Guilin rice noodles, a unique cuisine from Guilin, Guangxi, is renowned both domestically and internationally as one of the top ten “Guilin Classics”. Utilizing a heat conduction model, this study explores the effec...Guilin rice noodles, a unique cuisine from Guilin, Guangxi, is renowned both domestically and internationally as one of the top ten “Guilin Classics”. Utilizing a heat conduction model, this study explores the effectiveness of the cooking process in sterilizing Guilin rice noodles before consumption. The model assumes that a large pot is filled with boiling water which is maintained at a constant high temperature heat resource through continuous gentle heating. And the room temperature is set as the initial temperature for the preheating process and the final temperature for the cooling process. The objective is to assess whether the cooking process achieves satisfactory sterilization results. The temperature distribution function of rice noodle with time is analytically obtained using the separation of variables method in the three-dimensional cylindrical coordinate system. Meanwhile, the thermal diffusion coefficient of Guilin rice noodles is obtained in terms of Riedel’ theory. By analyzing the elimination characteristics of Pseudomonas cocovenenans subsp. farinofermentans, this study obtains the optimal time required for effective sterilization at the core of Guilin rice noodles. The results show that the potential Pseudomonas cocovenenans subsp. farinofermentans will be completely eliminated through continuously preheating more than 31 seconds during the cooking process before consumption. This study provides a valuable reference of food safety standards in the cooking process of Guilin rice noodles, particularly in ensuring the complete inactivation of potentially harmful strains such as Pseudomonas cocovenenans subsp. farinofermentans.展开更多
A 1D radially self-consistent model in helicon plasmas has been established to investigate the influence of radial heat conduction on plasma transport and wave propagation.Two kinds of 1D radial fluid models,with and ...A 1D radially self-consistent model in helicon plasmas has been established to investigate the influence of radial heat conduction on plasma transport and wave propagation.Two kinds of 1D radial fluid models,with and without considering heat conduction,have been developed to couple the 1D plasma-wave interaction model,and self-consistent solutions have been obtained.It is concluded that in the low magnetic field range the radial heat conduction plays a moderate role in the transport of helicon plasmas and the importance depends on the application of the helicon source.It influences the local energy balance leading to enhancement of the electron temperature in the bulk region and a decrease in plasma density.The power deposition in the plasma is mainly balanced by collisional processes and axial diffusion,whereas it is compensated by heat conduction in the bulk region and consumed near the boundary.The role of radial heat conduction in the large magnetic field regime becomes negligible and the two fluid models show consistency.The local power balance,especially near the wall,is improved when conductive heat is taken into account.展开更多
Aerogel nanoporous materials possess high porosity, high specific surface area, and extremely low density due to their unique nanoscale network structure. Moreover, their effective thermal conductivity is very low, ma...Aerogel nanoporous materials possess high porosity, high specific surface area, and extremely low density due to their unique nanoscale network structure. Moreover, their effective thermal conductivity is very low, making them a new type of lightweight and highly efficient nanoscale super-insulating material. However, prediction of their effective thermal conductivity is challenging due to their uneven pore size distribution. To investigate the internal heat transfer mechanism of aerogel nanoporous materials, this study constructed a cross-aligned and cubic pore model(CACPM) based on the actual pore arrangement of SiO_(2) aerogel. Based on the established CACPM, the effective thermal conductivity expression for the aerogel was derived by simultaneously considering gas-phase heat conduction, solid-phase heat conduction, and radiative heat transfer. The derived expression was then compared with available experimental data and the Wei structure model. The results indicate that, according to the model established in this study for the derived thermal conductivity formula of silica aerogel, for powdery silica aerogel under the conditions of T = 298 K, a_(2)= 0.85, D_(1)= 90 μm, ρ = 128 kg/m^(3), within the pressure range of 0–10^(5)Pa, the average deviation between the calculated values and experimental values is 10.51%. In the pressure range of 10^(3)–10^(4)Pa, the deviation between calculated values and experimental values is within 4%. Under these conditions, the model has certain reference value in engineering verification. This study also makes a certain contribution to the research of aerogel thermal conductivity heat transfer models and calculation formulae.展开更多
Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an applic...Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an application pertaining to the safety of light water nuclear reactors.Postulating a core meltdown accident,the behaviour of the core melt(aka corium)into a steel vessel is of tremendous importance when evaluating the vessel integrity.Evaluating correctly the heat fluxes requires the numerical simulation of the interaction between the liquid material and its solid counterpart which forms during the solidification process,but also may melt back.To simulate this configuration,encoun-tered in various industrial applications,one considers a bi-phase model constituted by a liquid phase in contact and interaction with its solid phase.The liquid phase may solidify in presence of low energetic source,while the solid phase may melt due to an intense heat flux from the high-temperature liquid.In the frame of the in-house legacy code,several simplifying assumptions(0D multi-layer discretization,instantaneous heat transfer via a quadratic temperature profile in solids)are made for the modelling of such phase changes.In the present work,these shortcomings are illustrated and further overcome by solving a 2D heat conduction model in the solid by a mixed Raviart-Thomas finite element method coupled to the liquid phase due to heat and mass exchanges through Stefan condition.The liquid phase is modeled with a 0D multi-layer approach.The 0D-liquid and 2D-solid mod-els are coupled by a Stefan like phase change interface model.Several sanity checks are performed to assess the validity of the approach on 1D and 2D academical configurations for which exact or reference solutions are available.Then more advanced situations(genu-ine multi-dimensional phase changes and an"industrial-like scenario")are simulated to verify the appropriate behavior of the obtained coupled simulation scheme.展开更多
A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh rat...A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.展开更多
We propose a new concept, the centre of energy, to study energy diffusion and heat conduction in a one-dimensional hard-point model. For the diatom model, we find an anomalous energy diffusion as (x2) - tβ with β ...We propose a new concept, the centre of energy, to study energy diffusion and heat conduction in a one-dimensional hard-point model. For the diatom model, we find an anomalous energy diffusion as (x2) - tβ with β = 1.33, which is independent of initial condition and mass rate. The present model can be viewed as the model composed by independent quasi-particles, the centre of energy. In this way, heat current can be calculated. Based on the theory of dynamic billiard, the divergent exponent of heat conductivity is estimated to be α = 0.33, which is confirmed by a simple numerical calculation.展开更多
The phonon relaxation and heat conduction in one-dimensional Fermi Pasta-Ulam (FPU) β lattices are studied by using molecular dynamics simulations. The phonon relaxation rate, which dominates the length dependence ...The phonon relaxation and heat conduction in one-dimensional Fermi Pasta-Ulam (FPU) β lattices are studied by using molecular dynamics simulations. The phonon relaxation rate, which dominates the length dependence of the FPU β lattice, is first calculated from the energy autoeorrelation function for different modes at various temperatures through equilibrium molecular dynamics simulations. We find that the relaxation rate as a function of wave number k is proportional to k^1.688, which leads to a N^0.41 divergence of the thermal conductivity in the framework of Green-Kubo relation. This is also in good agreement with the data obtained by non-equilibrium molecular dynamics simulations which estimate the length dependence exponent of the thermal conductivity as 0.415. Our results confirm the N^2/5 divergence in one-dimensional FPU β lattices. The effects of the heat flux on the thermal conductivity are also studied by imposing different temperature differences on the two ends of the lattices. We find that the thermal conductivity is insensitive to the heat flux under our simulation conditions. It implies that the linear response theory is applicable towards the heat conduction in one-dimensional FPU β lattices.展开更多
The welding mechanism of laser-TIG hybrid welding process is analyzed. Withthe variation of arc current, the welding process is divided into two patterns: deep-penetrationwelding and heat conductive welding. The heat ...The welding mechanism of laser-TIG hybrid welding process is analyzed. Withthe variation of arc current, the welding process is divided into two patterns: deep-penetrationwelding and heat conductive welding. The heat flow model of hybrid welding is presented. As todeep-penetration welding, the heat source includes a surface heat flux and a volume heat flux. Theheat source of heat conductive welding is composed of two Gaussian distribute surface heat sources.With this heat source model, a temperature field is calculated. The finite element code MARC isemployed for this purpose. The calculation results show a good agreement with the experimental data.展开更多
This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials (FGMs). First, a Laplace transform approach is used to handle the time variable. Then, a f...This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials (FGMs). First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest's algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.展开更多
Heat conduction through conventional and interlocking building bricks with cavities was studied in this work. Heat transfer analysis was carried out using MATLAB? partial differential equation toolbox. Regular and sta...Heat conduction through conventional and interlocking building bricks with cavities was studied in this work. Heat transfer analysis was carried out using MATLAB? partial differential equation toolbox. Regular and staggered hole arrangements were studied. Results showed that four staggered holed interlocking bricks were effective in thermal resistance into the bricks and increasing the holes beyond four did not give any thermal resistance advantage. For the conventional bricks staggered holes did not give any thermal resistance advantage but the four-holed bricks were also adjudged to be effective in thermal resistance into the brick surface. Increasing the number of holes beyond four in conventional bricks did give some thermal resistivity advantage but very minimal. Structural strengths of holed bricks were not considered in this study.展开更多
Studies on heat conduction are so far mainly focused on regular systems such as the one-dimensional(1D) and twodimensional(2D) lattices where atoms are regularly connected and temperatures of atoms are homogeneous...Studies on heat conduction are so far mainly focused on regular systems such as the one-dimensional(1D) and twodimensional(2D) lattices where atoms are regularly connected and temperatures of atoms are homogeneously distributed.However, realistic systems such as the nanotube/nanowire networks are not regular but heterogeneously structured, and their heat conduction remains largely unknown. We present a model of quasi-physical networks to study heat conduction in such physical networks and focus on how the network structure influences the heat conduction coefficient κ. In this model,we for the first time consider each link as a 1D chain of atoms instead of a spring in the previous studies. We find that κ is different from link to link in the network, in contrast to the same constant in a regular 1D or 2D lattice. Moreover, for each specific link, we present a formula to show how κ depends on both its link length and the temperatures on its two ends.These findings show that the heat conduction in physical networks is not a straightforward extension of 1D and 2D lattices but seriously influenced by the network structure.展开更多
Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditi...Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.展开更多
We build a fractional dual-phase-lag model and the corresponding bioheat transfer equation, which we use to interpret the experiment results for processed meat that have been explained by applying the hyperbolic condu...We build a fractional dual-phase-lag model and the corresponding bioheat transfer equation, which we use to interpret the experiment results for processed meat that have been explained by applying the hyperbolic conduction. Analytical solutions expressed by H-functions are obtained by using the Laplace and Fourier transforms method. The inverse fractional dual-phase-lag heat conduction problem for the simultaneous estimation of two relaxation times and orders of fractionality is solved by applying the nonlinear least-square method. The estimated model parameters are given. Finally, the measured and the calculated temperatures versus time are compared and discussed. Some numerical examples are also given and discussed.展开更多
In this study, we consider the heat-induced withdrawal reflex caused by exposure to an electromagnetic beam. We propose a concise dose-response relation for predicting the occurrence of withdrawal reflex from a given ...In this study, we consider the heat-induced withdrawal reflex caused by exposure to an electromagnetic beam. We propose a concise dose-response relation for predicting the occurrence of withdrawal reflex from a given spatial temperature profile. Our model is distilled from sub-step components in the ADT CHEETEH-E model developed at the Institute for Defense Analyses. Our model has only two parameters: the activation temperature of nociceptors and the critical threshold on the activated volume. When the spatial temperature profile is measurable, the two parameters can be determined from test data. We connect this dose-response relation to a temperature evolution model for electromagnetic heating. The resulting composite model governs the process from the electromagnetic beam deposited on the skin to the binary outcome of subject’s reflex response. We carry out non-dimensionalization in the time evolution model. The temperature solution of the non-dimensional system is the product of the applied power density and a parameter-free function. The effects of physical parameters are contained in non-dimensional time and depth. Scaling the physical temperature distribution into a parameter-free function greatly simplifies the analytical solution, and helps to pinpoint the effects of beam spot area and applied power density. With this formulation, we study the theoretical behaviors of the system, including the time of reflex, effect of heat conduction, biological latency in observed reflex, energy consumption by the time of reflex, and the strategy of selecting test conditions in experiments for the purpose of inferring model parameters from test data.展开更多
Taking the Lindemann model as a sample system in which there exist chemical reactions, diffusion and heat conduction, we found the theoretical framework of linear stability analysis for a unidimensional nonhomogeneous...Taking the Lindemann model as a sample system in which there exist chemical reactions, diffusion and heat conduction, we found the theoretical framework of linear stability analysis for a unidimensional nonhomogeneous two-variable system with one end subject to Dirichlet conditions, while the other end no-flux conditions. Furthermore, the conditions for the emergence of temperature waves are found out by the linear stability analysis and verified by a diagram for successive steps of evolution of spatial profile of temperature during a period that is plotted by numerical simulations on a computer. Without doubt, these results are in favor of the heat balance in chemical reactor designs.展开更多
This paper presents an improved understanding of coupled hydro-thermo-mechanical(HTM) hydraulic fracturing of quasi-brittle rock using the bonded particle model(BPM) within the discrete element method(DEM). BPM has be...This paper presents an improved understanding of coupled hydro-thermo-mechanical(HTM) hydraulic fracturing of quasi-brittle rock using the bonded particle model(BPM) within the discrete element method(DEM). BPM has been recently extended by the authors to account for coupled convective econductive heat flow and transport, and to enable full hydro-thermal fluidesolid coupled modeling.The application of the work is on enhanced geothermal systems(EGSs), and hydraulic fracturing of hot dry rock(HDR) is studied in terms of the impact of temperature difference between rock and a flowing fracturing fluid. Micro-mechanical investigation of temperature and fracturing fluid effects on hydraulic fracturing damage in rocks is presented. It was found that fracture is shorter with pronounced secondary microcracking along the main fracture for the case when the convectiveeconductive thermal heat exchange is considered. First, the convection heat exchange during low-viscosity fluid infiltration in permeable rock around the wellbore causes significant rock cooling, where a finger-like fluid infiltration was observed. Second, fluid infiltration inhibits pressure rise during pumping and delays fracture initiation and propagation. Additionally, thermal damage occurs in the whole area around the wellbore due to rock cooling and cold fluid infiltration. The size of a damaged area around the wellbore increases with decreasing fluid dynamic viscosity. Fluid and rock compressibility ratio was found to have significant effect on the fracture propagation velocity.展开更多
Thermal performance of a loop heat pipe with two evaporators and two condensers was examined using a lumped network model analysis. Thermosyphon-type vertical loop heat pipe and capillary-pump-type horizontal loop hea...Thermal performance of a loop heat pipe with two evaporators and two condensers was examined using a lumped network model analysis. Thermosyphon-type vertical loop heat pipe and capillary-pump-type horizontal loop heat pipe were calculated by examining the change of heating rate of two evaporators. Calculation results showed that the vapor and liquid flow rates in the loop heat pipe and the thermal conductance of the heat pipe changed significantly depending on the distribution ratio of the heating rate of the multiple evaporators. The thermal performance of the vertical loop heat pipe with two evaporators was also examined and experimental results of flow direction and thermal conductance of the heat pipe agreed with the analytical results. The lumped network model analysis is therefore considered accurate and preferable for the practical design of a loop heat pipe with multiple evaporators.展开更多
We investigate full counting statistics of quantum heat transfer in a collective-qubit system constructed by multiqubits interacting with two thermal baths. The nonequilibrium polaron-transformed Redfield approach emb...We investigate full counting statistics of quantum heat transfer in a collective-qubit system constructed by multiqubits interacting with two thermal baths. The nonequilibrium polaron-transformed Redfield approach embedded with an auxiliary counting field is applied to obtain the steady state heat current and fluctuations, which enables us to study the impact of the qubit–bath interaction in a wide regime. The heat current, current noise, and skewness are all found to clearly unify the limiting results in the weak and strong couplings. Moreover, the superradiant heat transfer is clarified as a system-size-dependent effect, and large number of qubits dramatically suppress the nonequilibrium superradiant signature.展开更多
Analytically solving a three-dimensional (3-D) bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green func...Analytically solving a three-dimensional (3-D) bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.展开更多
文摘Guilin rice noodles, a unique cuisine from Guilin, Guangxi, is renowned both domestically and internationally as one of the top ten “Guilin Classics”. Utilizing a heat conduction model, this study explores the effectiveness of the cooking process in sterilizing Guilin rice noodles before consumption. The model assumes that a large pot is filled with boiling water which is maintained at a constant high temperature heat resource through continuous gentle heating. And the room temperature is set as the initial temperature for the preheating process and the final temperature for the cooling process. The objective is to assess whether the cooking process achieves satisfactory sterilization results. The temperature distribution function of rice noodle with time is analytically obtained using the separation of variables method in the three-dimensional cylindrical coordinate system. Meanwhile, the thermal diffusion coefficient of Guilin rice noodles is obtained in terms of Riedel’ theory. By analyzing the elimination characteristics of Pseudomonas cocovenenans subsp. farinofermentans, this study obtains the optimal time required for effective sterilization at the core of Guilin rice noodles. The results show that the potential Pseudomonas cocovenenans subsp. farinofermentans will be completely eliminated through continuously preheating more than 31 seconds during the cooking process before consumption. This study provides a valuable reference of food safety standards in the cooking process of Guilin rice noodles, particularly in ensuring the complete inactivation of potentially harmful strains such as Pseudomonas cocovenenans subsp. farinofermentans.
基金National Natural Science Foundation of China(No.51907039)Shenzhen Technology Project(Nos.JCYJ20190806142603534 and ZDSYS201707280904031)+1 种基金ESPEOS project(No.PID2019108034RB-I00/AEI/10.13039/501100011033)funded by the Agencia Estatal de Investigacion(Spanish National Research Agency)。
文摘A 1D radially self-consistent model in helicon plasmas has been established to investigate the influence of radial heat conduction on plasma transport and wave propagation.Two kinds of 1D radial fluid models,with and without considering heat conduction,have been developed to couple the 1D plasma-wave interaction model,and self-consistent solutions have been obtained.It is concluded that in the low magnetic field range the radial heat conduction plays a moderate role in the transport of helicon plasmas and the importance depends on the application of the helicon source.It influences the local energy balance leading to enhancement of the electron temperature in the bulk region and a decrease in plasma density.The power deposition in the plasma is mainly balanced by collisional processes and axial diffusion,whereas it is compensated by heat conduction in the bulk region and consumed near the boundary.The role of radial heat conduction in the large magnetic field regime becomes negligible and the two fluid models show consistency.The local power balance,especially near the wall,is improved when conductive heat is taken into account.
基金supported by the National Natural Science Foundation of China (Grant Nos. 51764046 and 52160013)the Inner Mongolia Autonomous Region Postgraduate Research Innovation Project of China (Grant No. S20231165Z)the Research Program of Science and Technology at Universities of Inner Mongolia Autonomous Region of China (Grant Nos. 2023RCTD016 and 2024RCTD008)。
文摘Aerogel nanoporous materials possess high porosity, high specific surface area, and extremely low density due to their unique nanoscale network structure. Moreover, their effective thermal conductivity is very low, making them a new type of lightweight and highly efficient nanoscale super-insulating material. However, prediction of their effective thermal conductivity is challenging due to their uneven pore size distribution. To investigate the internal heat transfer mechanism of aerogel nanoporous materials, this study constructed a cross-aligned and cubic pore model(CACPM) based on the actual pore arrangement of SiO_(2) aerogel. Based on the established CACPM, the effective thermal conductivity expression for the aerogel was derived by simultaneously considering gas-phase heat conduction, solid-phase heat conduction, and radiative heat transfer. The derived expression was then compared with available experimental data and the Wei structure model. The results indicate that, according to the model established in this study for the derived thermal conductivity formula of silica aerogel, for powdery silica aerogel under the conditions of T = 298 K, a_(2)= 0.85, D_(1)= 90 μm, ρ = 128 kg/m^(3), within the pressure range of 0–10^(5)Pa, the average deviation between the calculated values and experimental values is 10.51%. In the pressure range of 10^(3)–10^(4)Pa, the deviation between calculated values and experimental values is within 4%. Under these conditions, the model has certain reference value in engineering verification. This study also makes a certain contribution to the research of aerogel thermal conductivity heat transfer models and calculation formulae.
基金funded by CEA,EDF and Framatomefinancial and scientific support of CEA Cadarache.
文摘Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an application pertaining to the safety of light water nuclear reactors.Postulating a core meltdown accident,the behaviour of the core melt(aka corium)into a steel vessel is of tremendous importance when evaluating the vessel integrity.Evaluating correctly the heat fluxes requires the numerical simulation of the interaction between the liquid material and its solid counterpart which forms during the solidification process,but also may melt back.To simulate this configuration,encoun-tered in various industrial applications,one considers a bi-phase model constituted by a liquid phase in contact and interaction with its solid phase.The liquid phase may solidify in presence of low energetic source,while the solid phase may melt due to an intense heat flux from the high-temperature liquid.In the frame of the in-house legacy code,several simplifying assumptions(0D multi-layer discretization,instantaneous heat transfer via a quadratic temperature profile in solids)are made for the modelling of such phase changes.In the present work,these shortcomings are illustrated and further overcome by solving a 2D heat conduction model in the solid by a mixed Raviart-Thomas finite element method coupled to the liquid phase due to heat and mass exchanges through Stefan condition.The liquid phase is modeled with a 0D multi-layer approach.The 0D-liquid and 2D-solid mod-els are coupled by a Stefan like phase change interface model.Several sanity checks are performed to assess the validity of the approach on 1D and 2D academical configurations for which exact or reference solutions are available.Then more advanced situations(genu-ine multi-dimensional phase changes and an"industrial-like scenario")are simulated to verify the appropriate behavior of the obtained coupled simulation scheme.
文摘A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10605020)the Natural Science Foundation of Zhejiang Province of China (Grant No. Y605376.)
文摘We propose a new concept, the centre of energy, to study energy diffusion and heat conduction in a one-dimensional hard-point model. For the diatom model, we find an anomalous energy diffusion as (x2) - tβ with β = 1.33, which is independent of initial condition and mass rate. The present model can be viewed as the model composed by independent quasi-particles, the centre of energy. In this way, heat current can be calculated. Based on the theory of dynamic billiard, the divergent exponent of heat conductivity is estimated to be α = 0.33, which is confirmed by a simple numerical calculation.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.50976052,51136001,and 50730006)the Program for New Century Excellent Talents in University,China+1 种基金the Tsinghua University Initiative Scientific Research Program,Chinathe Tsinghua National Laboratory for Information Science and Technology TNList Cross-discipline Foundation,China
文摘The phonon relaxation and heat conduction in one-dimensional Fermi Pasta-Ulam (FPU) β lattices are studied by using molecular dynamics simulations. The phonon relaxation rate, which dominates the length dependence of the FPU β lattice, is first calculated from the energy autoeorrelation function for different modes at various temperatures through equilibrium molecular dynamics simulations. We find that the relaxation rate as a function of wave number k is proportional to k^1.688, which leads to a N^0.41 divergence of the thermal conductivity in the framework of Green-Kubo relation. This is also in good agreement with the data obtained by non-equilibrium molecular dynamics simulations which estimate the length dependence exponent of the thermal conductivity as 0.415. Our results confirm the N^2/5 divergence in one-dimensional FPU β lattices. The effects of the heat flux on the thermal conductivity are also studied by imposing different temperature differences on the two ends of the lattices. We find that the thermal conductivity is insensitive to the heat flux under our simulation conditions. It implies that the linear response theory is applicable towards the heat conduction in one-dimensional FPU β lattices.
文摘The welding mechanism of laser-TIG hybrid welding process is analyzed. Withthe variation of arc current, the welding process is divided into two patterns: deep-penetrationwelding and heat conductive welding. The heat flow model of hybrid welding is presented. As todeep-penetration welding, the heat source includes a surface heat flux and a volume heat flux. Theheat source of heat conductive welding is composed of two Gaussian distribute surface heat sources.With this heat source model, a temperature field is calculated. The finite element code MARC isemployed for this purpose. The calculation results show a good agreement with the experimental data.
文摘This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials (FGMs). First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest's algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.
文摘Heat conduction through conventional and interlocking building bricks with cavities was studied in this work. Heat transfer analysis was carried out using MATLAB? partial differential equation toolbox. Regular and staggered hole arrangements were studied. Results showed that four staggered holed interlocking bricks were effective in thermal resistance into the bricks and increasing the holes beyond four did not give any thermal resistance advantage. For the conventional bricks staggered holes did not give any thermal resistance advantage but the four-holed bricks were also adjudged to be effective in thermal resistance into the brick surface. Increasing the number of holes beyond four in conventional bricks did give some thermal resistivity advantage but very minimal. Structural strengths of holed bricks were not considered in this study.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11135001 and 11375066)the National Basic Research Program of China(Grant No.2013CB834100)
文摘Studies on heat conduction are so far mainly focused on regular systems such as the one-dimensional(1D) and twodimensional(2D) lattices where atoms are regularly connected and temperatures of atoms are homogeneously distributed.However, realistic systems such as the nanotube/nanowire networks are not regular but heterogeneously structured, and their heat conduction remains largely unknown. We present a model of quasi-physical networks to study heat conduction in such physical networks and focus on how the network structure influences the heat conduction coefficient κ. In this model,we for the first time consider each link as a 1D chain of atoms instead of a spring in the previous studies. We find that κ is different from link to link in the network, in contrast to the same constant in a regular 1D or 2D lattice. Moreover, for each specific link, we present a formula to show how κ depends on both its link length and the temperatures on its two ends.These findings show that the heat conduction in physical networks is not a straightforward extension of 1D and 2D lattices but seriously influenced by the network structure.
基金supported by National Natural Science Foundation of China(11101244,11271231)National Tackling Key Problems Program(20050200069)Doctorate Foundation of the Ministry of Education of China(20030422047)
文摘Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.
基金supported by the National Natural Science Foundation of China(Grant Nos.11102102,11472161,and 91130017)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2014AQ015)the Independent Innovation Foundation of Shandong University,China(Grant No.2013ZRYQ002)
文摘We build a fractional dual-phase-lag model and the corresponding bioheat transfer equation, which we use to interpret the experiment results for processed meat that have been explained by applying the hyperbolic conduction. Analytical solutions expressed by H-functions are obtained by using the Laplace and Fourier transforms method. The inverse fractional dual-phase-lag heat conduction problem for the simultaneous estimation of two relaxation times and orders of fractionality is solved by applying the nonlinear least-square method. The estimated model parameters are given. Finally, the measured and the calculated temperatures versus time are compared and discussed. Some numerical examples are also given and discussed.
文摘In this study, we consider the heat-induced withdrawal reflex caused by exposure to an electromagnetic beam. We propose a concise dose-response relation for predicting the occurrence of withdrawal reflex from a given spatial temperature profile. Our model is distilled from sub-step components in the ADT CHEETEH-E model developed at the Institute for Defense Analyses. Our model has only two parameters: the activation temperature of nociceptors and the critical threshold on the activated volume. When the spatial temperature profile is measurable, the two parameters can be determined from test data. We connect this dose-response relation to a temperature evolution model for electromagnetic heating. The resulting composite model governs the process from the electromagnetic beam deposited on the skin to the binary outcome of subject’s reflex response. We carry out non-dimensionalization in the time evolution model. The temperature solution of the non-dimensional system is the product of the applied power density and a parameter-free function. The effects of physical parameters are contained in non-dimensional time and depth. Scaling the physical temperature distribution into a parameter-free function greatly simplifies the analytical solution, and helps to pinpoint the effects of beam spot area and applied power density. With this formulation, we study the theoretical behaviors of the system, including the time of reflex, effect of heat conduction, biological latency in observed reflex, energy consumption by the time of reflex, and the strategy of selecting test conditions in experiments for the purpose of inferring model parameters from test data.
文摘Taking the Lindemann model as a sample system in which there exist chemical reactions, diffusion and heat conduction, we found the theoretical framework of linear stability analysis for a unidimensional nonhomogeneous two-variable system with one end subject to Dirichlet conditions, while the other end no-flux conditions. Furthermore, the conditions for the emergence of temperature waves are found out by the linear stability analysis and verified by a diagram for successive steps of evolution of spatial profile of temperature during a period that is plotted by numerical simulations on a computer. Without doubt, these results are in favor of the heat balance in chemical reactor designs.
基金Financial support provided by the U.S. Department of Energy under DOE Grant No. DE-FE0002760
文摘This paper presents an improved understanding of coupled hydro-thermo-mechanical(HTM) hydraulic fracturing of quasi-brittle rock using the bonded particle model(BPM) within the discrete element method(DEM). BPM has been recently extended by the authors to account for coupled convective econductive heat flow and transport, and to enable full hydro-thermal fluidesolid coupled modeling.The application of the work is on enhanced geothermal systems(EGSs), and hydraulic fracturing of hot dry rock(HDR) is studied in terms of the impact of temperature difference between rock and a flowing fracturing fluid. Micro-mechanical investigation of temperature and fracturing fluid effects on hydraulic fracturing damage in rocks is presented. It was found that fracture is shorter with pronounced secondary microcracking along the main fracture for the case when the convectiveeconductive thermal heat exchange is considered. First, the convection heat exchange during low-viscosity fluid infiltration in permeable rock around the wellbore causes significant rock cooling, where a finger-like fluid infiltration was observed. Second, fluid infiltration inhibits pressure rise during pumping and delays fracture initiation and propagation. Additionally, thermal damage occurs in the whole area around the wellbore due to rock cooling and cold fluid infiltration. The size of a damaged area around the wellbore increases with decreasing fluid dynamic viscosity. Fluid and rock compressibility ratio was found to have significant effect on the fracture propagation velocity.
文摘Thermal performance of a loop heat pipe with two evaporators and two condensers was examined using a lumped network model analysis. Thermosyphon-type vertical loop heat pipe and capillary-pump-type horizontal loop heat pipe were calculated by examining the change of heating rate of two evaporators. Calculation results showed that the vapor and liquid flow rates in the loop heat pipe and the thermal conductance of the heat pipe changed significantly depending on the distribution ratio of the heating rate of the multiple evaporators. The thermal performance of the vertical loop heat pipe with two evaporators was also examined and experimental results of flow direction and thermal conductance of the heat pipe agreed with the analytical results. The lumped network model analysis is therefore considered accurate and preferable for the practical design of a loop heat pipe with multiple evaporators.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11874011 and 11704093)
文摘We investigate full counting statistics of quantum heat transfer in a collective-qubit system constructed by multiqubits interacting with two thermal baths. The nonequilibrium polaron-transformed Redfield approach embedded with an auxiliary counting field is applied to obtain the steady state heat current and fluctuations, which enables us to study the impact of the qubit–bath interaction in a wide regime. The heat current, current noise, and skewness are all found to clearly unify the limiting results in the weak and strong couplings. Moreover, the superradiant heat transfer is clarified as a system-size-dependent effect, and large number of qubits dramatically suppress the nonequilibrium superradiant signature.
基金Project supported by the National Natural Science Foundation of China (No. 50776097)
文摘Analytically solving a three-dimensional (3-D) bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.