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 purpose of this paper is to predict the plate heat exchanger performance when axial conduction in plates and in flow channels are present and fluids' viscosities are temperature dependent.

The approach to achieve the objective of the paper is deriving the governing equations and developing a computer program based on finite differences to solve them. The governing equations become dimensionless defining reference values and then discretized using FTBCS and FTCS methods. To solve the governing equations, the flow channel is divided into small elements in axial direction. Physical properties are constant for each element, while viscosity changes from one element to another one.

The effect of axial conduction in plates as well as in flow channels on temperature distributions, are studied individually and simultaneously. The program is run under four different conditions, namely: no axial conduction, axial conduction in the plates and in the flow channels, axial conduction in the plates only, and axial conduction in the flow channels only.In all the above cases, temperature distributions are achieved and characteristic curves are plotted. The numerical results are validated by comparing them with those published in an established reference carried out ignoring the effect of axial conduction, using the same plate geometry and flow details.

This paper gives valuable information and offers practical help to plate heat exchanger design engineer in order to choose the proper material for the plates as well as the right service and product fluids.

High temperature gradient induces high residual stress, producing an important effect on the part manufacturing during laser powder bed fusion (LPBF). The purpose of this study is to investigate the effect of the molten pool mode on the thermal stress of Ti-6Al-4V alloy during different deposition processes.

A coupled thermal-mechanical finite element model was built. The developed model was validated by comparing the numerical results with the experimental data in the maximum molten pool temperature, the molten pool dimension and the residual stress described in the previous work.

For the single-track process, the keyhole mode caused an increase in both the maximum stress and the high-stress area compared with the conduction mode. For the multitrack process, a lower tensile stress around the scanning track and a higher compressive stress below the scanning track were found in the keyhole mode. For the multilayer process, the stress along the scanning direction at the middle of the part changed from tensile stress to compressive stress with the increase in the deposition layer number. As the powder layer number increased, the stress along the scanning direction near the top surface of the part decreased while the stress along the deposition direction obviously increased, indicating that the stress along the deposition direction became the dominant stress. The keyhole mode can reduce the residual stress near the top of the part, and the conduction mode was more likely to produce a low residual stress near the bottom of the part.

The results provide a systematic understanding of thermal stress during the LPBF process.

The solution of radiative heat transfer problems in participating media is often obtained using the standard discrete ordinates method (SDOM). This method produces anomalies caused by ray effects if radiative boundary conditions have discontinuities or abrupt changes. Ray effects may be mitigated using the modified discrete ordinates method (MDOM), which is based on superposition of the solutions obtained by considering separately radiation from the walls and radiation from the medium. The purpose of this paper is to study the role of ray effects in combined conduction‐radiation problems and investigate the superiority of the MDOM over SDOM.

The MDOM has been used to calculate radiative heat transfer in irregular geometries using body‐fitted coordinates. Here, the blocked‐off region concept, originally developed in computational fluid dynamics, is used along with the finite volume method and SDOM or MDOM to solve combined conduction‐radiation heat transport problems in irregular geometries. Enclosures with an absorbing, emitting and isotropically or anisotropically scattering medium are analyzed.

The results confirm the capability of the MDOM to minimize the anomalies due to ray effects in combined heat transfer problems, and demonstrate that MDOM is more computationally efficient than SDOM.

The paper demonstrates the application of MDOM to combined conduction‐radiation heat transfer problems in irregular geometries using blocked‐off method.

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.

This study aims to perform a comprehensive investigation to model the thermal characteristics of a coupled conduction-radiation heat transfer in a two-dimensional irregular enclosure including a triangular-shaped heat source.

For this purpose, a promising hybrid technique based on the concepts of blocked-off method, FVM and DOM is developed. The enclosure consists of several horizontal, vertical and oblique walls, and thermal conductivity within the enclosure varies directly with temperature and indirectly with position. To simplify the complex geometry, a promising mathematical model is introduced using blocked-off method. Emitting, absorbing and non-isotropic scattering gray are assumed as the main radiative characteristics of the steady medium.

DOM and FVM are, respectively, applied for solving radiative transfer equation (RTE) and the energy equation, which includes conduction, radiation and heat source terms. The temperature and heat flux distributions are calculated inside the enclosure. For validation, results are compared with previous data reported in the literature under the same conditions. Results and comparisons show that this approach is highly efficient and reliable for complex geometries with coupled conduction-radiation heat transfer. Finally, the effects of thermo-radiative parameters including surface emissivity, extinction coefficient, scattering albedo, asymmetry factor and conduction-radiation parameter on temperature and heat flux distributions are studied.

In this paper, a hybrid numerical method is used to analyze coupled conduction-radiation heat transfer in an irregular geometry. Varying thermal conductivity is included in this analysis. By applying the method, results obtained for temperature and heat flux distributions are presented and also validated by the data provided by several previous papers.

Present study aims to extend the lattice Boltzmann method (LBM) to simulate radiation in geometries with curved boundaries, as the first step to simulate radiation in complex porous media. In recent years, researchers have increasingly explored the use of porous media to improve the heat transfer processes. The lattice Boltzmann method (LBM) is one of the most effective techniques for simulating heat transfer in such media. However, the application of the LBM to study radiation in complex geometries that contain curved boundaries, as found in many porous media, has been limited.

The numerical evaluation of the effect of the radiation-conduction parameter and extinction coefficient on temperature and incident radiation distributions demonstrates that the proposed LBM algorithm provides highly accurate results across all cases, compared to those found in the literature or those obtained using the finite volume method (FVM) with the discrete ordinates method (DOM) for radiative information.

For the case with a conduction-radiation parameter equal to 0.01, the maximum relative error is 1.9% in predicting temperature along vertical central line. The accuracy improves with an increase in the conduction-radiation parameter. Furthermore, the comparison between computational performances of two approaches reveals that the LBM-LBM approach performs significantly faster than the FVM-DOM solver.

The difficulty of radiative modeling in combined problems involving irregular boundaries has led to alternative approaches that generally increase the computational expense to obtain necessary radiative details. To address the limitations of existing methods, this study presents a new approach involving a coupled lattice Boltzmann and first-order blocked-off technique to efficiently model conductive-radiative heat transfer in complex geometries with participating media. This algorithm has been developed using the parallel lattice Boltzmann solver.

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.