Finite Difference-Peridynamic Differential Operator for Solving Transient Heat Conduction Problems

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.

Galerkin-based quasi-smooth manifold element(QSME)method for anisotropic heat conduction problems in composites with complex geometry

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.

Variational Approach to 2D and 3D Heat Conduction Modeling

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.

Variational Approach to Heat Conduction Modeling

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.

Sterilization Effect of Cooking Process for Guilin Rice Noodles Based on Heat Conduction Model

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.

Layout Design of Conductive Heat Channel by Emulating Natural Branch Systems

To design effective and easy-to-manufacture conductive heat channels, a heuristic method by emulating the natural branch systems is suggested. The design process of the method is divided into two steps, which are the principal channel design and the lateral channel design. During the process, the width of each channel is controlled by the bifurcation law, and the end point of the channel is located at the point with the maximum temperature while the start points of the principal channel and the lateral channel are respectively determined by the location of the heat sink and the law of the minimum thermal resistance. Four design examples with different boundary conditions are studied by the suggested method, and the design results are compared with that of the traditional structural topology optimization method. Not only lower maximum temperature and relatively uniform distribution of temperature are obtained by the suggested method, but also straight channels are achieved without gray element, which is easy to manufacture. The suggested method inspired by the natural branch systems can provide an effective solution for heat channel design in the heat dissipation structures.

Transient Radiation-Conductive Heat transfer Problems:"The Quadrupole Method"

This paper presents a statement of the works performed in L.E.M.T.A by the members of the thermal and mechanical heterogeneous media research group during the last six years concerning the solving of coupled conductive and radiative heat transfers within a multilayer and semi-transparent“wall”.Out of the authors,this paper allows to take inspiration from the works of D.Maillet,M.Lazard and V.Manias[19,20,21].The aim of these works is to represent in a macroscopic way,with the minimum number of thermophysical parameters,the heat transfers in a plane system composed of semi-transparent media.The approach we propose is semi-analytic(Kernel substitution technique,Laplace transformation)and allow to obtain in the Laplace domain an analytical solution that can be easily used.This method can be applied in two main scopes of applications:the estimation of thermophysical properties(phononic conductivity,optical thickness,Planck number for instance)of semi-transparent materials(glasses,crystals,glass wool,semi-conductors,synthetic diamonds,vitroceramics and so on)and the modelling of processes with semitransparent walls(for instance bottles forming,flat glass production,drying of paper).The method will be first presented and validated and two examples of applications will be then given.This method can be applied to semitransparent walls that emit,absorb and scatter the radiant energy(participating medium).It appears from the principle of a Kernel substitution technique applied to the radiative flux expression and initially introduced by Lick[1]that allows to change the character of the governing heat equation from the integro-differential form to a purely differential one.In the case of limiting cases of purely scattering and purely absorbing media,the solution of the radiative transfer equation is exact.In the general case,we make a two-flux approximation.In all cases,we assume a linear transfer and use the Laplace transform.The method can be applied to grey or grey by bands media,with isotropic or anisotropic scattering.The advantage of the method is fast computational times for good precision.

Lightweight diamond/Cu interface tuning for outstanding heat conduction

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.

NONLINEAR STABILITY OF RAREFACTION WAVES TO THE COMPRESSIBLE NAVIER-STOKES EQUATIONS FOR A REACTING MIXTURE WITH ZERO HEAT CONDUCTIVITY

In this paper,we study the time-asymptotically nonlinear stability of rarefaction waves for the Cauchy problem of the compressible Navier-Stokes equations for a reacting mixture with zero heat conductivity in one dimension.If the corresponding Riemann problem for the compressible Euler system admits the solutions consisting of rarefaction waves only,it is shown that its Cauchy problem has a unique global solution which tends time-asymptotically towards the rarefaction waves,while the initial perturbation and the strength of rarefaction waves are suitably small.

Investigation of radial heat conduction with 1D self-consistent model in helicon plasmas

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.

A Novel Localized Meshless Method for Solving Transient Heat Conduction Problems in Complicated Domains

This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by a finite difference scheme,and thus the governing equation of transient heat transfer is transformed into a non-homogeneous modified Helmholtz equation.Secondly,the solution of the non-homogeneous modified Helmholtz equation is decomposed into a particular solution and a homogeneous solution.And then,the DRM and LKM are used to solve the particular solution of the non-homogeneous equation and the homogeneous solution of the modified Helmholtz equation,respectively.The LKM is a recently proposed local radial basis function collocationmethod with themerits of being simple,accurate,and free ofmesh and integration.Compared with the traditional domain-type and boundary-type schemes,the present coupling algorithm could be treated as a really good alternative for the analysis of transient heat conduction on high-dimensional and complicated domains.Numerical experiments,including two-and three-dimensional heat transfer models,demonstrated the effectiveness and accuracy of the new methodology.

Research on the heat conduction mechanism of induction brazecoating

Induction brazecoating technology is an important means to improve the surface properties of materials.In this paper,copper plate and corundum are selected as substrates for induction brazecoating respectively.The temperature variation of powder and paste coating is systematically studied,and the heat conduction mode and path in the brazecoating process are analyzed.The results show that rise of the coating temperature mainly depends on the heat absorption from the substrate.The liquid-solid interface conducts heat violently and advances step by step,which promotes the melting spread of metal filler metal.The powdery brazecoating material is in a free state,and there is a gas insulation film between the powder particles and the diamond,making it difficult for substrate to conduct heat to the coating.The binder not only assists forming,shell making and oxygen isolation,but also plays an important role during melting.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Dynamic response of a thermal transistor to time-varying signals

Thermal transistor,the thermal analog of an electronic transistor,is one of the most important thermal devices for microscopic-scale heat manipulating.It is a three-terminal device,and the heat current flowing through two terminals can be largely controlled by the temperature of the third one.Dynamic response plays an important role in the application of electric devices and also thermal devices,which represents the devices’ability to treat fast varying inputs.In this paper,we systematically study two typical dynamic responses of a thermal transistor,i.e.,the response to a step-function input(a switching process)and the response to a square-wave input.The role of the length L of the control segment is carefully studied.It is revealed that when L is increased,the performance of the thermal transistor worsens badly.Both the relaxation time for the former process and the cutoff frequency for the latter one follow the power-law dependence on L quite well,which agrees with our analytical expectation.However,the detailed power exponents deviate from the expected values noticeably.This implies the violation of the conventional assumptions that we adopt.

Decay rates and time lags of heat conduction in building construction under field conditions

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.

Solution of an Inverse Problem of Heat Conduction of 45 Steel with Martensite Phase Transformation in High Pressure during Gas Quenching

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.

HYBRID FEM WITH FUNDAMENTAL SOLUTIONS AS TRIAL FUNCTIONS FOR HEAT CONDUCTION SIMULATION

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.

Analytical solutions of transient heat conduction in multilayered slabs and application to thermal analysis of landfills

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.

Fast precise integration method for hyperbolic 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 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.

A Numerical Method on Inverse Determination of Heat Transfer Coefficient Based on Thermographic Temperature Measurement

The heat transfer coefficient in a multidimensional heat conduction problem is obtained from the solution of the inverse heat conduction problem based on the thermographic temperature measurement. The modified one-dimensional correction method （MODCM）, along with the finite volume method, is employed for both two- and three-dimensional inverse problems. A series of numerical experiments are conducted in order to verify the effectiveness of the method. In addition, the effect of the temperature measurement error, the ending criterion of the iteration, etc. on the result of the inverse problem is investigated. It is proved that the method is a simple, stable and accurate one that can solve successfully the inverse heat conduction problem.