A CLASS OF TWO-LEVEL EXPLICIT DIFFERENCE SCHEMES FOR SOLVING THREE DIMENSIONAL HEAT CONDUCTION EQUATION

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 new complex variable meshless method for transient heat conduction problems

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.

Fuzzy finite difference method for heat conduction analysis with uncertain parameters

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.

Isogeometric Boundary Element Analysis for 2D Transient Heat Conduction Problem with Radial Integration Method

This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.

Analysis of steady heat conduction for 3D axisymmetric functionally graded circular plate

The thermal conduction behavior of the three-dimensional axisymmetric functionally graded circular plate was studied under thermal loads on its top and bottom surfaces. Material properties were taken to be arbitrary distribution functions of the thickness. A temperature function that satisfies thermal boundary conditions at the edges and the variable separation method were used to reduce equation governing the steady state heat conduction to an ordinary differential equation （ODE） in the thickness coordinate which was solved analytically. Next, resulting variable coefficients ODE due to arbitrary distribution of material properties along thickness coordinate was also solved by the Peano-Baker series. Some numerical examples were given to demonstrate the accuracy, efficiency of the present model, mad to investigate the influence of different distributions of material properties on the temperature field. The numerical results confirm that the influence of different material distributions, gradient indices and thickness of plate to temperature field in plate can not be ignored.

Meshless analysis of three-dimensional steady-state heat conduction problems

Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson＇s equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method （RKPM）. A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.