We present a minimal theoretical model for self-sustained oscillations of a thin elastic sheet on a hot plate,induced by thermomechanical coupling.As the plate temperature increases,the sheet’s static deflection beco...We present a minimal theoretical model for self-sustained oscillations of a thin elastic sheet on a hot plate,induced by thermomechanical coupling.As the plate temperature increases,the sheet’s static deflection becomes unstable via a Hopf bifurcation at a critical temperature TC,giving rise to spontaneous periodic motion.Linear stability analysis yields analytical expressions for the critical oscillation temperature TC and the oscillation period at onset.Numerical simulations of the nonlinear equations confirm the bifurcation and reveal how key parameters(stiffness,thermal softening,thermal coupling,etc.)govern the oscillation amplitude and waveform.Finally,we demonstrate that the self-oscillating sheet can perform mechanical work as a heat engine,and we compare its performance to the Carnot efficiency limit.This work provides design principles for thermally driven selfoscillators with potential applications in soft robotics,adaptive structures,and thermal energy harvesting.展开更多
With the miniaturization of devices and the development of modern heating technologies,the generalization of heat conduction and thermoelastic coupling has become crucial,effectively emulating the thermodynamic behavi...With the miniaturization of devices and the development of modern heating technologies,the generalization of heat conduction and thermoelastic coupling has become crucial,effectively emulating the thermodynamic behavior of materials in ultrashort time scales.Theoretically,generalized heat conductive models are considered in this work.By analogy with mechanical viscoelastic models,this paper further enriches the heat conduction models and gives their one-dimensional physical expression.Numerically,the transient thermoelastic response of the slim strip material under thermal shock is investigated by applying the proposed models.First,the analytical solution in the Laplace domain is obtained by the Laplace transform.Then,the numerical results of the transient responses are obtained by the numerical inverse Laplace transform.Finally,the transient responses of different models are analyzed and compared,and the effects of material parameters are discussed.This work not only opens up new research perspectives on generalized heat conductive and thermoelastic coupling theories,but also is expected to be beneficial for the deeper understanding of the heat wave theory.展开更多
The paper proposes an approximate solution to the classical (parabolic) multidimensional 2D and 3D heat conduction equation for a 5 × 5 cm aluminium plate and a 5 × 5 × 5 cm aluminum cube. An approximat...The paper proposes an approximate solution to the classical (parabolic) multidimensional 2D and 3D heat conduction equation for a 5 × 5 cm aluminium plate and a 5 × 5 × 5 cm aluminum cube. An approximate solution of the generalized (hyperbolic) 2D and 3D equation for the considered plate and cube is also proposed. Approximate solutions were obtained by applying calculus of variations and Euler-Lagrange equations. In order to verify the correctness of the proposed approximate solutions, they were compared with the exact solutions of parabolic and hyperbolic equations. The paper also presents the research on the influence of time parameters τ as well as the relaxation times τ ∗ to the variation of the profile of the temperature field for the considered aluminum plate and cube.展开更多
The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element ...The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.展开更多
Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using t...Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.展开更多
It is known that Fourier’s heat equation, which is parabolic, implies an infinite velocity propagation, or, in other words, that the mechanism of heat conduction is established instantaneously under all conditions. T...It is known that Fourier’s heat equation, which is parabolic, implies an infinite velocity propagation, or, in other words, that the mechanism of heat conduction is established instantaneously under all conditions. This is unacceptable on physical grounds in spite of the fact that Fourier’s law agrees well with experiment. However, discrepancies are likely to occur when extremely short distances or extremely short time intervals are considered, as they must in some modern problems of aero-thermodynamics. Cattaneo and independently Vernotte proved that such process can be described by Heaviside’s telegraph equation. This paper shows that this fact can be derived using calculus of variations, by application of the Euler-Lagrange equation. So, we proved that the equation of heat conduction with finite velocity propagation of the thermal disturbance can be obtained as a solution to one variational problem.展开更多
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.展开更多
This paper consists of two parts. (1) For a hollow sphere with sudden temperature changes on its inner and outer surfaces, the hyperbolic heat conduction equation is employed to describe this extreme thermal case and...This paper consists of two parts. (1) For a hollow sphere with sudden temperature changes on its inner and outer surfaces, the hyperbolic heat conduction equation is employed to describe this extreme thermal case and an analytical expression of its temperature distribution is obtained. According to the expression, the non-Fourier heat conduction behavior that will appear in the hollow sphere is studied and some qualitative conditions that will result in distinct non-Fourier behavior in the medium is ultimately attained. (2) A novel experiment to observe non-Fourier heat conduction behavior in porous material (mainly ordinary duplicating paper) heated by a microsecond laser pulse is presented. The conditions for observing distinct non-Fourier heat conduction behavior in the experimental sample agree well with the theoretical results qualitatively.展开更多
The field measurements of decay rates and time lags of heat conduction in a building construction taken in Nanjing during the summer of 2001 are presented.The decay rates and time lags are calculated according to the ...The field measurements of decay rates and time lags of heat conduction in a building construction taken in Nanjing during the summer of 2001 are presented.The decay rates and time lags are calculated according to the frequency responses of the heat absorbed by the room＇s internal surfaces,inside surface temperature,indoor air temperature and outdoor synthetic temperature.The measured results match very well with the theoretical results of the zeroth and the first order values of the decay rates and time lags of heat conduction in the building construction,but the difference between the measured values and the theoretical values for the second order is too great to be accepted.It is therefore difficult to accurately test the second order value.However,it is still advisable to complete the analysis using the zeroth-and the first-orders values of the decay rates and time lags of heat conduction in building construction under field conditions,because in these cases the decay rates of heat conduction reach twenty which meets the requirements of engineering plans.展开更多
According to the axisymmetric heat conduction of monolayer cylinder, a general method was deduced to calculate the axisymmetric temperature of linear heat conduction multilayer cylinder. Four types of boundary conditi...According to the axisymmetric heat conduction of monolayer cylinder, a general method was deduced to calculate the axisymmetric temperature of linear heat conduction multilayer cylinder. Four types of boundary conditions were summarized and formulas for each type were derived. Then, a general calculating program was developed. Four temperature formulas could be expressed by a uniform equation, and the intermediate interface temperatures of axisymmetrical linear conduction multilayer cylinder satisfied tridiagonal linear and nonlinear systems of equations, which could be solved with the pursuit method and the Newton＇s method, respectively. With the calculating program, the temperature at any point of linear heat conduction multilayer cylinder could be obtained.展开更多
A new type of hybrid finite element formulation with fundamental solutions as internal interpolation functions, named as HFS-FEM, is presented in this paper and used for solving two dimensional heat conduction problem...A new type of hybrid finite element formulation with fundamental solutions as internal interpolation functions, named as HFS-FEM, is presented in this paper and used for solving two dimensional heat conduction problems in single and multi-layer materials. In the proposed approach, a new variational functional is firstly constructed for the proposed HFS-FE model and the related existence of extremum is presented. Then, the assumed internal potential field constructed by the linear combination of fundamental solutions at points outside the elemental domain under consideration is used as the internal interpolation function, which analytically satisfies the governing equation within each element. As a result, the domain integrals in the variational functional formulation can be converted into the boundary integrals which can significantly simplify the calculation of the element stiffness matrix. The independent frame field is also introduced to guarantee the inter-element continuity and the stationary condition of the new variational functional is used to obtain the final stiffness equations. The proposed method inherits the advantages of the hybrid Trefftz finite element method （HT-FEM） over the conventional finite element method （FEM） and boundary element method （BEM）, and avoids the difficulty in selecting appropriate terms of T-complete functions used in HT-FEM, as the fundamental solutions contain usually one term only, rather than a series containing infinitely many terms. Further, the fundamental solutions of a problem are, in general, easier to derive than the T-complete functions of that problem. Finally, several examples are presented to assess the performance of the proposed method, and the obtained numerical results show good numerical accuracy and remarkable insensitivity to mesh distortion.展开更多
In order to simulate thermal strains, thermal stresses, residual stresses and microstructure of the steel during gas quenching by means of the numerical method, it is necessary to obtain an accurate boundary condition...In order to simulate thermal strains, thermal stresses, residual stresses and microstructure of the steel during gas quenching by means of the numerical method, it is necessary to obtain an accurate boundary condition of temperature field. The surface heat transfer coefficient is a key parameter. The explicit finite difference method, nonlinear estimation method and the experimental relation between temperature and time during gas quenching have been used to solve the inverse problem of heat conduction. The relationship between surface temperature and surface heat transfer coefficient of a cylinder has been given. The nonlinear surface heat transfer coefficients include the coupled effects between martensitic phase transformation and temperature.展开更多
The study of transient heat conduction in multilayered slabs is widely used in various engineering fields. In this paper, the transient heat conduction in multilayered slabs with general boundary conditions and arbitr...The study of transient heat conduction in multilayered slabs is widely used in various engineering fields. In this paper, the transient heat conduction in multilayered slabs with general boundary conditions and arbitrary heat generations is analysed. The boundary conditions are general and include various combinations of Dirichlet, Neumann or Robin boundary conditions at either surface. Moreover, arbitrary heat generations in the slabs are taken into account. The solutions are derived by basic methods, including the superposition method, separation variable method and orthogonal expansion method. The simplified double-layered analytical solution is validated by a numerical method and applied to predicting the temporal and spatial distribution of the temperature inside a landfill. It indicates the ability of the proposed analytical solutions for solving the wide range of applied transient heat conduction problems.展开更多
A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to t...A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.展开更多
Inverse heat conduction method (IHCM) is one of the most effective approaches to obtaining the boiling heat transfer coefficient from measured results. This paper focuses on its application in cryogenic boiling heat t...Inverse heat conduction method (IHCM) is one of the most effective approaches to obtaining the boiling heat transfer coefficient from measured results. This paper focuses on its application in cryogenic boiling heat transfer. Experiments were conducted on the heat transfer of a stainless steel block in a liquid nitrogen bath, with the assumption of a 1D conduction condition to realize fast acquisition of the temperature of the test points inside the block. With the inverse-heat conduction theory and the explicit finite difference model, a solving program was developed to calculate the heat flux and the boiling heat transfer coefficient of a stainless steel block in liquid nitrogen bath based on the temperature acquisition data. Considering the oscillating data and some unsmooth transition points in the inverse-heat-conduction calculation result of the heat-transfer coefficient, a two-step data-fitting procedure was proposed to obtain the expression for the boiling heat transfer coefficients. The coefficient was then verified for accuracy by a comparison between the simulation results using this expression and the verifying experimental results of a stainless steel block. The maximum error with a revised segment fitting is around 6%, which verifies the feasibility of using IHCM to measure the boiling heat transfer coefficient in liquid nitrogen bath.展开更多
With rapid developments in the field of very large-scale integrated circuits,heat dissipation has emerged as a significant factor that restricts the high-density integration of chips.Due to their high thermal conducti...With rapid developments in the field of very large-scale integrated circuits,heat dissipation has emerged as a significant factor that restricts the high-density integration of chips.Due to their high thermal conductivity and low thermal expansion coefficient,diamond/Cu composites have attracted considerable attention as a promising thermal management material.In this study,a surface tungsten carbide gradient layer coating of diamond particles has been realized using comprehensive magnetron sputtering technology and a heat treatment process.Diamond/Cu composites were prepared using high-temperature and high-pressure technology.The results show that,by adjusting the heat treatment process,tungsten carbide and di-tungsten carbide are generated by an in situ reaction at the tungsten–diamond interface,and W–WC–W_(2)C gradient layer-coated diamond particles were obtained.The diamond/Cu composites were sintered by high-temperature and high-pressure technology,and the density of surface-modified diamond/Cu composites was less than 4 g cm^(-3).The W–WC–W_(2)C@diamond/Cu composites have a thermal diffusivity as high as 331 mm^(2)s^(-1),and their thermal expansion coefficient is as low as 1.76×10^(-6)K^(-1).The interface coherent structure of the gradient layer-coated diamond/copper composite can effectively improve the interface heat transport efficiency.展开更多
In this paper, based on the improved complex variable moving least-square （ICVMLS） approximation, a new complex variable meshless method （CVMM） for two-dimensional （2D） transient heat conduction problems is pres...In this paper, based on the improved complex variable moving least-square （ICVMLS） approximation, a new complex variable meshless method （CVMM） for two-dimensional （2D） transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.展开更多
A new numerical technique named as fuzzy finite difference method is proposed to solve the heat conduction problems with fuzzy uncertainties in both the phys- ical parameters and initial/boundary conditions. In virtue...A new numerical technique named as fuzzy finite difference method is proposed to solve the heat conduction problems with fuzzy uncertainties in both the phys- ical parameters and initial/boundary conditions. In virtue of the level-cut method, the difference discrete equations with fuzzy parameters are equivalently transformed into groups of interval equations. New stability analysis theory suited to fuzzy difference schemes is developed. Based on the parameter perturbation method, the interval ranges of the uncertain temperature field can be approximately predicted. Subsequently, fuzzy solutions to the original difference equations are obtained by the fuzzy resolution theorem. Two numerical examples are given to demonstrate the feasibility and efficiency of the presented method for solving both steady-state and transient heat conduction problems.展开更多
Honeycomb sandwich materials have been widely used in aerospace industry as secondary structural materials or interior panels. Recently the research of full thermoplastic honeycomb sandwich materials has been of inter...Honeycomb sandwich materials have been widely used in aerospace industry as secondary structural materials or interior panels. Recently the research of full thermoplastic honeycomb sandwich materials has been of interest due to their recyclability. In their production, they usually demand an adhesive-free process, namely fusion bonding, to connect thermoplastic honeycomb core and face sheets. It is a heat induced process where a parameter of temperature should be well controlled to guarantee the product quality. This article presents a mathematical model of heat conduction to analyze the transient temperature distribution from heating tools towards inner part of the core under fusion bonding conditions. In order to simplify the complexity of 3D honeycomb geometry, a homogenization method is used to obtain average thermal properties of the honeycomb along the major heat flux direction. The model is validated by comparing with the results of in-situ temperature measurement during fusion bonding. The presented model can also be applied to analyzing general out-of-plane heat conduction through honeycomb sandwich structures made from other materials.展开更多
In this work, laser heat conduction lap welding(LHCLW) of AZ31B magnesium alloy sheet and DP780galvanized steel sheet was carried out by the defocused laser beam. The effects of laser power on the microstructure and m...In this work, laser heat conduction lap welding(LHCLW) of AZ31B magnesium alloy sheet and DP780galvanized steel sheet was carried out by the defocused laser beam. The effects of laser power on the microstructure and mechanical properties of the joint were studied. The pros and cons of the joint were identified and evaluated by measuring the tensile shear strength, microhardness and microstructure observation. The formation mechanism of various phases at the Mg/steel interface was analyzed. The results indicated that the galvanized layer could promote the metallurgical bonding between magnesium alloy and steel by improving the diffusion ability of molten magnesium alloy at the steel interface and reacting with Mg, so as to enhance the strength of the joint. A continuous dense layered eutectic structure(α-Mg+MgZn) was formed at the interface of the joint, while MgZn_(2)and MgZn phase was formed at the weld edge zone and heat affective zone(HAZ), whereas no reaction layer was generated between the uncoated steel and magnesium alloy. A sound joint could be obtained at 2.5 kW, and the corresponding tensile shear strength reached the maximum value of 42.9 N/mm. The strength was slightly reduced at 2.6 kW due to the existence of microcracks in the eutectic reaction layer.展开更多
基金supported by the Guangdong Basic and Applied Basic Research Foundation(Grant Nos.2025B1515020077 and 2024A15150301-39)the National Natural Science Foundation of China(Grant No.12205138)the Shenzhen Science and Technology Innovation Committee(Grant No.JCYJ2022-0530113206015).
文摘We present a minimal theoretical model for self-sustained oscillations of a thin elastic sheet on a hot plate,induced by thermomechanical coupling.As the plate temperature increases,the sheet’s static deflection becomes unstable via a Hopf bifurcation at a critical temperature TC,giving rise to spontaneous periodic motion.Linear stability analysis yields analytical expressions for the critical oscillation temperature TC and the oscillation period at onset.Numerical simulations of the nonlinear equations confirm the bifurcation and reveal how key parameters(stiffness,thermal softening,thermal coupling,etc.)govern the oscillation amplitude and waveform.Finally,we demonstrate that the self-oscillating sheet can perform mechanical work as a heat engine,and we compare its performance to the Carnot efficiency limit.This work provides design principles for thermally driven selfoscillators with potential applications in soft robotics,adaptive structures,and thermal energy harvesting.
基金Project supported by the Guangdong Basic and Applied Basic Research Foundation of China(No.2023A1515012809)the Natural Science Foundation of Shaanxi Province of China(No.2023-JC-YB-073)the Fundamental Research Funds for the Central Universities of China(No.D5000230066)。
文摘With the miniaturization of devices and the development of modern heating technologies,the generalization of heat conduction and thermoelastic coupling has become crucial,effectively emulating the thermodynamic behavior of materials in ultrashort time scales.Theoretically,generalized heat conductive models are considered in this work.By analogy with mechanical viscoelastic models,this paper further enriches the heat conduction models and gives their one-dimensional physical expression.Numerically,the transient thermoelastic response of the slim strip material under thermal shock is investigated by applying the proposed models.First,the analytical solution in the Laplace domain is obtained by the Laplace transform.Then,the numerical results of the transient responses are obtained by the numerical inverse Laplace transform.Finally,the transient responses of different models are analyzed and compared,and the effects of material parameters are discussed.This work not only opens up new research perspectives on generalized heat conductive and thermoelastic coupling theories,but also is expected to be beneficial for the deeper understanding of the heat wave theory.
文摘The paper proposes an approximate solution to the classical (parabolic) multidimensional 2D and 3D heat conduction equation for a 5 × 5 cm aluminium plate and a 5 × 5 × 5 cm aluminum cube. An approximate solution of the generalized (hyperbolic) 2D and 3D equation for the considered plate and cube is also proposed. Approximate solutions were obtained by applying calculus of variations and Euler-Lagrange equations. In order to verify the correctness of the proposed approximate solutions, they were compared with the exact solutions of parabolic and hyperbolic equations. The paper also presents the research on the influence of time parameters τ as well as the relaxation times τ ∗ to the variation of the profile of the temperature field for the considered aluminum plate and cube.
基金Project supported by the National Natural Science Foundation of China (Nos. 12102043, 12072375U2241240)the Natural Science Foundation of Hunan Province of China (Nos. 2023JJ40698 and 2021JJ40710)。
文摘The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.
基金This work was financially supported by the Key Science and Technology Project of Longmen Laboratory(No.LMYLKT-001)Innovation and Entrepreneurship Training Program for College Students of Henan Province(No.202310464050)。
文摘Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.
文摘It is known that Fourier’s heat equation, which is parabolic, implies an infinite velocity propagation, or, in other words, that the mechanism of heat conduction is established instantaneously under all conditions. This is unacceptable on physical grounds in spite of the fact that Fourier’s law agrees well with experiment. However, discrepancies are likely to occur when extremely short distances or extremely short time intervals are considered, as they must in some modern problems of aero-thermodynamics. Cattaneo and independently Vernotte proved that such process can be described by Heaviside’s telegraph equation. This paper shows that this fact can be derived using calculus of variations, by application of the Euler-Lagrange equation. So, we proved that the equation of heat conduction with finite velocity propagation of the thermal disturbance can be obtained as a solution to one variational problem.
文摘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.
基金Supported by the Chinese Academy of Sciences (No. KJ 951-B1-704), the National Natural Science Foundation of China (No. 59736130) and the State Key Fundamental Research Plan of China (No. G2000026305).
文摘This paper consists of two parts. (1) For a hollow sphere with sudden temperature changes on its inner and outer surfaces, the hyperbolic heat conduction equation is employed to describe this extreme thermal case and an analytical expression of its temperature distribution is obtained. According to the expression, the non-Fourier heat conduction behavior that will appear in the hollow sphere is studied and some qualitative conditions that will result in distinct non-Fourier behavior in the medium is ultimately attained. (2) A novel experiment to observe non-Fourier heat conduction behavior in porous material (mainly ordinary duplicating paper) heated by a microsecond laser pulse is presented. The conditions for observing distinct non-Fourier heat conduction behavior in the experimental sample agree well with the theoretical results qualitatively.
基金The Advance Research Projects of Southeast Universityfor the National Natural Science Foundation of China(No.XJ0701262)the National Key Technologies R&D Program of China during the 11th Five-Year Plan Period(No.2008BAJ12B04,2008BAJ12B05,2006BAJ03A04)
文摘The field measurements of decay rates and time lags of heat conduction in a building construction taken in Nanjing during the summer of 2001 are presented.The decay rates and time lags are calculated according to the frequency responses of the heat absorbed by the room＇s internal surfaces,inside surface temperature,indoor air temperature and outdoor synthetic temperature.The measured results match very well with the theoretical results of the zeroth and the first order values of the decay rates and time lags of heat conduction in the building construction,but the difference between the measured values and the theoretical values for the second order is too great to be accepted.It is therefore difficult to accurately test the second order value.However,it is still advisable to complete the analysis using the zeroth-and the first-orders values of the decay rates and time lags of heat conduction in building construction under field conditions,because in these cases the decay rates of heat conduction reach twenty which meets the requirements of engineering plans.
基金Item Sponsored by National Natural Science Foundation of China (50474014)Provincial Key Technologies Research and Development Program of Liaoning of China(2008216005)
文摘According to the axisymmetric heat conduction of monolayer cylinder, a general method was deduced to calculate the axisymmetric temperature of linear heat conduction multilayer cylinder. Four types of boundary conditions were summarized and formulas for each type were derived. Then, a general calculating program was developed. Four temperature formulas could be expressed by a uniform equation, and the intermediate interface temperatures of axisymmetrical linear conduction multilayer cylinder satisfied tridiagonal linear and nonlinear systems of equations, which could be solved with the pursuit method and the Newton＇s method, respectively. With the calculating program, the temperature at any point of linear heat conduction multilayer cylinder could be obtained.
文摘A new type of hybrid finite element formulation with fundamental solutions as internal interpolation functions, named as HFS-FEM, is presented in this paper and used for solving two dimensional heat conduction problems in single and multi-layer materials. In the proposed approach, a new variational functional is firstly constructed for the proposed HFS-FE model and the related existence of extremum is presented. Then, the assumed internal potential field constructed by the linear combination of fundamental solutions at points outside the elemental domain under consideration is used as the internal interpolation function, which analytically satisfies the governing equation within each element. As a result, the domain integrals in the variational functional formulation can be converted into the boundary integrals which can significantly simplify the calculation of the element stiffness matrix. The independent frame field is also introduced to guarantee the inter-element continuity and the stationary condition of the new variational functional is used to obtain the final stiffness equations. The proposed method inherits the advantages of the hybrid Trefftz finite element method （HT-FEM） over the conventional finite element method （FEM） and boundary element method （BEM）, and avoids the difficulty in selecting appropriate terms of T-complete functions used in HT-FEM, as the fundamental solutions contain usually one term only, rather than a series containing infinitely many terms. Further, the fundamental solutions of a problem are, in general, easier to derive than the T-complete functions of that problem. Finally, several examples are presented to assess the performance of the proposed method, and the obtained numerical results show good numerical accuracy and remarkable insensitivity to mesh distortion.
基金This work has been supported by the National Natural Science Foundation of China (10162002) and Foundation for University Key Teacher by the Ministry of Education and The Yunnan Foundation of Natural Science (1999A0023M).
文摘In order to simulate thermal strains, thermal stresses, residual stresses and microstructure of the steel during gas quenching by means of the numerical method, it is necessary to obtain an accurate boundary condition of temperature field. The surface heat transfer coefficient is a key parameter. The explicit finite difference method, nonlinear estimation method and the experimental relation between temperature and time during gas quenching have been used to solve the inverse problem of heat conduction. The relationship between surface temperature and surface heat transfer coefficient of a cylinder has been given. The nonlinear surface heat transfer coefficients include the coupled effects between martensitic phase transformation and temperature.
基金Projects(41530637,41877222,41702290)supported by the National Natural Science Foundation of China
文摘The study of transient heat conduction in multilayered slabs is widely used in various engineering fields. In this paper, the transient heat conduction in multilayered slabs with general boundary conditions and arbitrary heat generations is analysed. The boundary conditions are general and include various combinations of Dirichlet, Neumann or Robin boundary conditions at either surface. Moreover, arbitrary heat generations in the slabs are taken into account. The solutions are derived by basic methods, including the superposition method, separation variable method and orthogonal expansion method. The simplified double-layered analytical solution is validated by a numerical method and applied to predicting the temporal and spatial distribution of the temperature inside a landfill. It indicates the ability of the proposed analytical solutions for solving the wide range of applied transient heat conduction problems.
基金supported by the National Natural Science Foundation of China (Nos. 10902020 and 10721062)
文摘A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.
基金supported by the National Natural Sciences Foundation of China (No. 50776075)
文摘Inverse heat conduction method (IHCM) is one of the most effective approaches to obtaining the boiling heat transfer coefficient from measured results. This paper focuses on its application in cryogenic boiling heat transfer. Experiments were conducted on the heat transfer of a stainless steel block in a liquid nitrogen bath, with the assumption of a 1D conduction condition to realize fast acquisition of the temperature of the test points inside the block. With the inverse-heat conduction theory and the explicit finite difference model, a solving program was developed to calculate the heat flux and the boiling heat transfer coefficient of a stainless steel block in liquid nitrogen bath based on the temperature acquisition data. Considering the oscillating data and some unsmooth transition points in the inverse-heat-conduction calculation result of the heat-transfer coefficient, a two-step data-fitting procedure was proposed to obtain the expression for the boiling heat transfer coefficients. The coefficient was then verified for accuracy by a comparison between the simulation results using this expression and the verifying experimental results of a stainless steel block. The maximum error with a revised segment fitting is around 6%, which verifies the feasibility of using IHCM to measure the boiling heat transfer coefficient in liquid nitrogen bath.
基金National Natural Science Foundation of China(Grant No.52072327)the China National Key R&D Program(2021YFB3701802)+6 种基金Scientific and Technological Projects of Henan Province(No.232102231050)the Higher Education and Teaching Reformation Project(2014SJGLX064)the Project for Work-station of Zhongyuan scholars of Henan Province(Nos.214400510002,224400510023)the Science and Technology Major Project of Henan Province(No.221100230300)the Postgraduate Education Reform and QualityAcademic Degrees&Graduate Education Reform Project of Henan Province(No.2021SJGLX060Y)the Postgraduate Education Reform and Quality Improvement Project of Henan Province(No.YJS2022JD34)the Science and Technology on Plasma Physics Laboratory(Grant No.JCKYS2021212010).
文摘With rapid developments in the field of very large-scale integrated circuits,heat dissipation has emerged as a significant factor that restricts the high-density integration of chips.Due to their high thermal conductivity and low thermal expansion coefficient,diamond/Cu composites have attracted considerable attention as a promising thermal management material.In this study,a surface tungsten carbide gradient layer coating of diamond particles has been realized using comprehensive magnetron sputtering technology and a heat treatment process.Diamond/Cu composites were prepared using high-temperature and high-pressure technology.The results show that,by adjusting the heat treatment process,tungsten carbide and di-tungsten carbide are generated by an in situ reaction at the tungsten–diamond interface,and W–WC–W_(2)C gradient layer-coated diamond particles were obtained.The diamond/Cu composites were sintered by high-temperature and high-pressure technology,and the density of surface-modified diamond/Cu composites was less than 4 g cm^(-3).The W–WC–W_(2)C@diamond/Cu composites have a thermal diffusivity as high as 331 mm^(2)s^(-1),and their thermal expansion coefficient is as low as 1.76×10^(-6)K^(-1).The interface coherent structure of the gradient layer-coated diamond/copper composite can effectively improve the interface heat transport efficiency.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Shanghai Leading Academic Discipline Project,China(Grant No.S30106)the Innovation Fund for Graduate Student of Shanghai University of China (Grant No.SHUCX120125)
文摘In this paper, based on the improved complex variable moving least-square （ICVMLS） approximation, a new complex variable meshless method （CVMM） for two-dimensional （2D） transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.
基金supported by the National Special Fund for Major Research Instrument Development(2011YQ140145)111 Project(B07009)+1 种基金the National Natural Science Foundation of China(11002013)Defense Industrial Technology Development Program(A2120110001 and B2120110011)
文摘A new numerical technique named as fuzzy finite difference method is proposed to solve the heat conduction problems with fuzzy uncertainties in both the phys- ical parameters and initial/boundary conditions. In virtue of the level-cut method, the difference discrete equations with fuzzy parameters are equivalently transformed into groups of interval equations. New stability analysis theory suited to fuzzy difference schemes is developed. Based on the parameter perturbation method, the interval ranges of the uncertain temperature field can be approximately predicted. Subsequently, fuzzy solutions to the original difference equations are obtained by the fuzzy resolution theorem. Two numerical examples are given to demonstrate the feasibility and efficiency of the presented method for solving both steady-state and transient heat conduction problems.
文摘Honeycomb sandwich materials have been widely used in aerospace industry as secondary structural materials or interior panels. Recently the research of full thermoplastic honeycomb sandwich materials has been of interest due to their recyclability. In their production, they usually demand an adhesive-free process, namely fusion bonding, to connect thermoplastic honeycomb core and face sheets. It is a heat induced process where a parameter of temperature should be well controlled to guarantee the product quality. This article presents a mathematical model of heat conduction to analyze the transient temperature distribution from heating tools towards inner part of the core under fusion bonding conditions. In order to simplify the complexity of 3D honeycomb geometry, a homogenization method is used to obtain average thermal properties of the honeycomb along the major heat flux direction. The model is validated by comparing with the results of in-situ temperature measurement during fusion bonding. The presented model can also be applied to analyzing general out-of-plane heat conduction through honeycomb sandwich structures made from other materials.
基金Projects(51705219, 51905227) supported by the National Natural Science Foundation of ChinaProject(BK20200915) supported by the Natural Science Foundation of Jiangsu Province,ChinaProject(19KJB460013) supported by the General University Science Research Project of Jiangsu Province,China。
文摘In this work, laser heat conduction lap welding(LHCLW) of AZ31B magnesium alloy sheet and DP780galvanized steel sheet was carried out by the defocused laser beam. The effects of laser power on the microstructure and mechanical properties of the joint were studied. The pros and cons of the joint were identified and evaluated by measuring the tensile shear strength, microhardness and microstructure observation. The formation mechanism of various phases at the Mg/steel interface was analyzed. The results indicated that the galvanized layer could promote the metallurgical bonding between magnesium alloy and steel by improving the diffusion ability of molten magnesium alloy at the steel interface and reacting with Mg, so as to enhance the strength of the joint. A continuous dense layered eutectic structure(α-Mg+MgZn) was formed at the interface of the joint, while MgZn_(2)and MgZn phase was formed at the weld edge zone and heat affective zone(HAZ), whereas no reaction layer was generated between the uncoated steel and magnesium alloy. A sound joint could be obtained at 2.5 kW, and the corresponding tensile shear strength reached the maximum value of 42.9 N/mm. The strength was slightly reduced at 2.6 kW due to the existence of microcracks in the eutectic reaction layer.
