Search results

This paper aims to review the existing bond-based peridynamic (PD) and state-based PD heat conduction models, and further propose a refined bond-based PD thermal conduction model by using the PD differential operator.

The general refined bond-based PD is established by replacing the local spatial derivatives in the classical heat conduction equations with their corresponding nonlocal integral forms obtained by the PD differential operator. This modeling approach is representative of the state-based PD models, whereas the resulting governing equations appear as the bond-based PD models.

The refined model can be reduced to the existing bond-based PD heat conduction models by specifying particular influence functions. Also, the refined model does not require any calibration procedure unlike the bond-based PD. A systematic explicit dynamic solver is introduced to validate 1 D, 2 D and 3 D heat conduction in domains with and without a crack subjected to a combination of Dirichlet, Neumann and convection boundary conditions. All of the PD predictions are in excellent agreement with the classical solutions and demonstrate the nonlocal feature and advantage of PD in dealing with heat conduction in discontinuous domains.

The existing PD heat conduction models are reviewed. A refined bond-based PD thermal conduction model by using PD differential operator is proposed and 3 D thermal conduction in intact or cracked structures is simulated.

The effect of the tube wall heat conduction on the natural convection in a tilted long cylindrical envelope with constant, but different temperature of the two ends and an adiabatic outer surface was numerically investigated. The envelope is supposed to be a simplified model for the pulse tube in a pulse tube refrigerator when the pulse tube is positioned at different orientations. It is found that the cylindrical envelope lateral wall heat conduction can enhance the heat transfer from the hot end to the cold end, not only because of the increase in pure heat conduction in the wall, but more importantly, also the intensification of the natural convection within the enclosure. This enhancement is resulted from the big temperature difference between the tube wall and the adjacent fluid near the hot and cold ends. Adoption of low thermal conductivity tube can effectively reduce such additional heat transfers from hot to cold end, thus reducing the loss of cooling capacity for the pulse tube refrigerator.

This paper aims to develop an efficient numerical method for nonlinear transient heat conduction problems with local radiation boundary conditions and nonlinear heat sources.

Based on the physical characteristic of the transient heat conduction and the distribution characteristic of the Green’s function, a quasi-superposition principle is presented for the transient heat conduction problems with local nonlinearities. Then, an efficient method is developed, which indicates that the solution of the original nonlinear problem can be derived by solving some nonlinear problems with small structures and a linear problem with the original structure. These problems are independent of each other and can be solved simultaneously by the parallel computing technique.

Within a small time step, the nonlinear thermal loads can only induce significant temperature responses of the regions near the positions of the nonlinear thermal loads, whereas the temperature responses of the remaining regions are very close to zero. According to the above physical characteristic, the original nonlinear problem can be transformed into some nonlinear problems with small structures and a linear problem with the original structure.

An efficient and accurate numerical method is presented for transient heat conduction problems with local nonlinearities, and some numerical examples demonstrate the high efficiency and accuracy of the proposed method.

The purpose of this paper is to present a new approach of meshless local B-spline based finite difference (FD) method for solving two dimensional transient heat conduction problems.

In the present method, any governing equations are discretized by B-spline approximation which is implemented in the spirit of FD technique using a local B-spline collocation scheme. The key aspect of the method is that any derivative is stated as neighbouring nodal values based on B-spline interpolants. The set of neighbouring nodes are allowed to be randomly distributed thus enhanced flexibility in the numerical simulation can be obtained. The method requires no mesh connectivity at all for either field variable approximation or integration. Time integration is performed by using the Crank-Nicolson implicit time stepping technique.

Several heat conduction problems in complex domains which represent for extended surfaces in industrial applications are examined to demonstrate the effectiveness of the present approach. Comparison of the obtained results with solutions from other numerical method available in literature is given. Excellent agreement with reference numerical method has been found.

The method is presented for 2D problems. Nevertheless, it would be also applicable for 3D problems.

A transient two dimensional heat conduction in complex domains which represent for extended surfaces in industrial applications is presented.

The presented new meshless local method is simple and accurate, while it is also suitable for analysis in domains of arbitrary geometries.

3D steady heat conduction analysis considering heat source is conducted on the fundamental of the fast multipole method (FMM)-accelerated line integration boundary element method (LIBEM).

Due to considering the heat source, domain integral is generated in the traditional heat conduction boundary integral equation (BIE), which will counteract the well-known merit of the BEM, namely, boundary-only discretization. To avoid volume discretization, the enhanced BEM, the LIBEM with dimension reduction property is introduced to transfer the domain integral into line integrals. Besides, owing to the unsatisfactory performance of the LIBEM when it comes to large-scale structures requiring massive computation, the FMM-accelerated LIBEM (FM-LIBEM) is proposed to improve the computation efficiency further.

Assuming N and M are the numbers of nodes and integral lines, respectively, the FM-LIBEM can reduce the time complexity from O(NM) to about O(N+ M), and a full discussion and verification of the advantage are done based on numerical examples under heat conduction.

(1) The LIBEM is applied to 3D heat conduction analysis with heat source. (2) The domain integrals can be transformed into boundary integrals with straight line integrals by the LIM. (3) A FM-LIBEM is proposed and can reduce the time complexity from O(NM) to O(N+ M). (4) The FM-LIBEM with high computational efficiency is exerted to solve 3D heat conduction analysis with heat source in massive computation successfully.

The purpose of this paper is to address unsteady heat conduction in two subsets of ordinary bodies. One subset consists of a large plane wall, a long cylinder and a sphere in one dimension. The other subset consists of a short cylinder and a large rectangular bar in two dimensions. The prevalent assumptions in the two subsets are: constant initial temperature, uniform surface heat flux and thermo-physical properties invariant with temperature. The engineering applications of the unsteady heat conduction deal with the determination of temperature–time histories in the two subsets using electric resistance heating, radiative heating and fire pool heating.

To this end, a novel numerical procedure named the enhanced method of discretization in time (EMDT) transforms the linear one-dimensional unsteady, heat conduction equations with non-homogeneous boundary conditions into equivalent nonlinear “quasi–steady” heat conduction equations having the time variable embedded as a time parameter. The equivalent nonlinear “quasi–steady” heat conduction equations are solved with a finite difference method.

Based on the numerical computations, it is demonstrated that the approximate temperature–time histories in the simple subset of ordinary bodies (large plane wall, long cylinder and sphere) exhibit a perfect matching over the entire time domain 0 < t < ∞ when compared against the rigorous exact temperature–time histories expressed by classical infinite series. Furthermore, using the method of superposition of solutions in the convoluted subset (short cylinder and large rectangular crossbar), the same level of agreement in the approximate temperature–time histories in the simple subset of ordinary bodies is evident.

The performance of the proposed EMDT coupled with a finite difference method is exhaustively assessed in the solution of the unsteady, one-dimensional heat conduction equations with prescribed surface heat flux for: a subset of one-dimensional bodies (plane wall, long cylinder and spheres) and a subset of two-dimensional bodies (short cylinder and large rectangular bar).

In this study, the purpose was to introduce two‐dimensional hyperbolic heat conduction equations in order to simulate the fast precooling process of a cylindrically shaped food product with internal heat generation. A modified model for internal heat generation due to respiration in the food product was proposed to take the effect of relaxation time into account. The obtained governing equations were solved numerically using an efficient finite difference technique. The influence of Biot number and heat generation parameters on thermal characteristics was examined and discussed. The results based on hyperbolic model were compared with the classical parabolic heat diffusion model. The present numerical code was validated via comparison with analytical solution and a good agreement was found.

The obtained governing equations were solved numerically using an efficient finite difference technique.

The influence of Biot number and heat generation parameters on thermal characteristics was examined and discussed. The results based on hyperbolic model were compared with the classical parabolic heat diffusion model. The present numerical code was validated via comparison with analytical solution and a good agreement was found.

Two‐dimensional analysis of fast precooling of cylindrical food product based on hyperbolic heat conduction model has not been investigated yet.

This investigation deals with the equilibration of heat conduction simulation in a very thin film using molecular dynamics. Two parameters, the positional order parameter and the kinetic H‐function, are employed simultaneously to monitor the evolution to the equilibrium. With the different boundary conditions, material parameters, and molecular lattice configurations, the results of the simulation show that the combination of the two parameters can give a consistent prediction to the approach of the equilibrium. At the equilibrium state, the process of heat conduction in a thin film is studied to understand the macroscopic behaviour from the standpoint of molecular dynamic motions. The method used can be applied to solve other microscopic flow and heat transfer problems using molecular dynamic simulation.

The present paper describes the applicability of hybrid transfinite element modelling/analysis formulations for non‐linear heat conduction problems involving phase change. The methodology is based on application of transform approaches and classical Galerkin schemes with finite element formulations to maintain the modelling versatility and numerical features for computational analysis. In addition, in conjunction with the above, the effects due to latent heat are modelled using enthalpy formulations to enable a physically realistic approximation to be effectively dealt computationally for materials exhibiting phase change within a narrow band of temperatures. Pertinent details of the approach and computational scheme adapted are described in technical detail. Numerical test cases of comparative nature are presented to demonstrate the applicability of the proposed formulations for numerical modelling/analysis of non‐linear heat conduction problems involving phase change.

The purpose of this study is to address one-dimensional, unsteady heat conduction in a large plane wall exchanging heat convection with a nearby fluid under “small time” conditions.

The Transversal Method of Lines (TMOL) was used to reformulate the unsteady, one-dimensional heat conduction equation in the space coordinate and time into a transformed “quasi-steady”, one-dimensional heat conduction equation in the space coordinate housing the time as an embedded parameter. The resulting ordinary differential equation of second order with heat convection boundary conditions is solved analytically with the method of undetermined coefficients.

Semi-analytical TMOL dimensionless temperature profiles of compact form with/without regressed terms are obtained for the whole spectrum of Biot number (0 < Bi < ∞) in the “small time” sub-domain. In addition, a new “large time” sub-domain is redefined, that is, setting a smaller critical dimensionless time or critical Fourier number τ_cr = 0.18.

The computed dimensionless center, surface and mean temperature profiles in the large plane wall accounting for all Biot number (0 < Bi < ∞) in the “small time” sub-domain τ < τ_cr = 0.18 exhibit excellent quality while carrying reasonable relative errors for engineering applications. The exemplary level of accuracy indicates that the traditional evaluation of the center, surface and mean temperatures with the standard infinite series retaining a large number of terms is no longer necessary.