Heat Conduction Research Articles

Penetration based lunar regolith thermal conductivity inversion: Method and verification

Based on the temperature change data of the penetrator, to achieve in-situ detection of the thermal conductivity of the lunar regolith profile, it is necessary to establish a heat conduction model between the penetrator and the lunar regolith. This study simplifies the heat conduction model of the complex shaped penetrator through simulation analysis results.Then,we proposed a thermal conductivity inversion method based on the transient thermal cylinder source model. The thermal diffusion test was carried out under normal temperature and pressure on a standard reference object with known thermal conductivity and the thermal conductivity inversion work was completed, which verified the feasibility of the inversion method. Then, we completed the thermal diffusion test of the penetrator and the lunar regolith simulant under the simulated low-temperature vacuum environment of the lunar surface, and carried out the thermal conductivity inversion of the lunar regolith simulant based on the thermal diffusion test data, which proved that the proposed thermal inversion method is applicable for the lunar regolith under the low-temperature vacuum. Finally, the thermal conductivity test was conducted on icy lunar regolith simulant with different water contents and a thermal diffusion test in a vacuum low-temperature environment was carried out. The thermal conductivity inversion was completed using the same method, which proved that the inversion method is suitable for icy lunar regolith simulant.

Functionalization of boron nitride with poly(catechol-polyamine) and bis(γ-triethoxysilylpropyl)tetrasulfide for preparing epoxy nanocomposites with enhanced thermal conductivity

Hexagonal boron nitride (hBN) has garnered significant interest as inorganic thermally conductive material in electronic packaging due to its high-efficiency thermal management capability and electrical insulation properties. However, the relatively inert chemical environment of the BN surface hinders phonon transfer at the interface. Additionally, the high content of inorganic components in the matrix leads to the deterioration of the composite’s mechanical properties. In this research, h-BN surface is modified using poly(catechol-polyamine) (PCPA) and bis (γ-triethoxysi-lylpropyl)tetrasulphide (Si69) (denoted as BN-PCPA-Si69). The noncovalent PCPA modification protects the endothermic properties of BN, Si69 provides covalent bonding between thermally conductive fillers and bisphenol A-type epoxy resin (Ep), and the interfacial compatibility between BN and epoxy resin is improved. The heat conductivity of final BN/Ep composite achieved 0.817 W/(m·K) containing 30 wt% BN. Correspondingly, the heat conductivity of the composite (mBN/Ep) achieved 0.881 W/(m·K) containing 30 wt% BN-PCPA-Si69. Compared to pure epoxy resin (0.200 W/(m·K), the thermal conductivity (TC) of the mBN/Ep composite increased by 364% while maintaining good electrical insulation properties. The proposed BN surface functionalization approach has potential applications in fields such as microelectronic packaging, where materials need to balance high thermal conductivity and electrical insulation.

Effects of OpenCL-Based Parallelization Methods on Explicit Numerical Methods to Solve the Heat Equation

In recent years, the need for high-performance computing solutions has increased due to the growing complexity of computational tasks. The use of parallel processing techniques has become essential to address this demand. In this study, an Open Computing Language (OpenCL)-based parallelization algorithm is implemented for the Constant Neighbors (CNe) and CNe with Predictor–Corrector (CpC) numerical methods, which are recently developed explicit and stable numerical algorithms to solve the heat conduction equation. The CPU time and error rate performance of these two methods are compared with the sequential implementation and Euler’s explicit method. The results demonstrate that the parallel version’s CPU time remains nearly constant under the examined circumstances, regardless of the number of spatial mesh points. This leads to a remarkable speed advantage over the sequential version for larger data point counts. Furthermore, the impact of the number of timesteps on the crossover point where the parallel version becomes faster than the sequential one is investigated.

Thermodynamically Consistent Evolution Equations in Continuum Mechanics

This paper addresses the modelling of material behaviour in terms of differential (or rate) equations. To comply with the objectivity principle, recourse is made to invariant fields in the Lagrangian description or to objective time derivatives in the Eulerian description. The thermodynamic consistency is investigated in terms of the Clausius–Duhem inequality with two unusual features. Firstly, the (non-negative) entropy production is viewed as a constitutive function per se. Secondly, the inequality is viewed as a constraint on the pertinent fields and it is solved by using a representation formula, which allows for the the admissibility of a class of models. For definiteness, models of heat conduction are established, within Lagrangian descriptions, while models of the Navier–Stokes–Voigt fluid are investigated within Eulerian descriptions. In connection with thermo-viscous fluids, evolution equations are investigated within the Eulerian description. It is shown that the thermodynamic consistency is compatible with both objective and non-objective evolution equations.

Enhancing thermoelectric generation: Integrating passive radiative cooling and concentrated solar heating with consideration of parasitic heat conduction

Thermoelectric generators (TEGs) integrated with solar energy and radiative cooling offer a promising approach for generating power. Concentrated solar energy enhances generation by increasing the solar flux density. However, the relationship between thermoelectric generation and concentration ratio remains not well understood. In this study, the finite element method is employed to investigate the temperature difference across TEGs and the change in generation performance with varying concentration ratios. The results show that due to the mismatch among cooling, heating, and parasitic heat conduction powers, the optimal power generation does not occur at the maximum concentration ratio. As the concentration ratio increases, the temperature difference across the single TEG and its output power initially increase and then decrease. To further improve the power generation performance at high concentration ratios, this study introduces TEG stacking strategy to enhance waste heat recovery by increasing thermal resistance. Outdoor experiments with self-made solar absorptive coating and radiative cooling coating validate this approach, achieving a temperature difference of 64.17 °C and a power density of 135.57 W m−2 by stacking three TEGs, outperforming a single TEG by over three times. This study offers a potential method to continuously power remote off-grid communities and low-power devices.

Investigation into the mechanism of ignition of magnesium alloy plates by high-temperature molten droplet

The ignition and spread of magnesium alloy fires are mainly attributed to the ignition of high-temperature molten metal droplets, but the underlying ignition mechanism remains unclear. This study addresses this knowledge gap by employing metallic aluminum and copper particles heated to temperatures over 1200 °C. Systematically, these heated particles were released at a fixed rate onto a magnesium alloy plate located inside a controlled heating furnace that maintained a constant temperature of 500 °C. Experimental results show that whether there are irregular pores and cracks in the resulting molten product can be used as a criterion for judging whether the magnesium alloy sample ignites. Furthermore, it is worth noting that molten copper droplets of the same size exhibit a higher heat-carrying capacity than their aluminum counterparts when both are heated to the same temperature. In addition, this study expanded the research scope through numerical simulation and analysis, using the PATO (Porous-Material Analysis Toolbox) solver to clarify the transient heat conduction process between the high-temperature droplet and the magnesium alloy plate. Simulations show that heat transfer upon droplet impact is primarily longitudinal, directed toward the base of the material. This results in a significantly higher temperature in the center of the affected magnesium alloy plate than at its periphery. In summary, this study helps to improve the understanding of the dangers of high-temperature molten droplets, and it is of great significance to investigate fire accidents caused by high-temperature molten droplets.

Three-phase lag model for thermal conductivity of a thermo-viscoelastic porous medium

This paper investigates the propagation of waves and the thermal behavior in a thermo-viscoelastic porous medium using the three-phase-lag (TPL) model which accounts for phase lags in the heat flux vector, temperature gradient and thermal displacement gradient. The study aims to capture the interactions between thermal, mechanical, and structural properties of the thermo-viscoelastic porous, isotropic, homogeneous medium. The governing equations, incorporating the TPL heat conduction law and constitutive relations for the viscoelastic porous material, are derived and analyzed. The impact of relaxation times, porosity, and viscoelastic parameters on heat propagation is investigated through analytical solutions by normal modal analysis. Analytical representations of various physical quantities (stresses components, displacements, temperature,…) in the material's domain are derived. These expressions are then assessed numerically for a specific material, with the outcomes depicted graphically. A comparison is made between the predictions of the three-phase-lag (TPL), the dual-phase-lag (DPL) and Lord–Shulman (L–S) theories in the absence and presence of voids.

Thermal modeling and its role in management of bifurcation aneurysms and associated ischemic complications in the middle cerebral artery

Bifurcated aneurysm of cerebral artery is a common cerebrovascular disease, which can easily lead to serious ischemic complications. Traditional treatments face challenges of risk and effectiveness. Therefore, it is particularly important to explore new therapeutic strategies. In recent years, the application of thermal modeling techniques in medical biological systems has provided a new perspective for the analysis of hemodynamics and tissue thermal response. The purpose of this study was to investigate the role of heat model in the treatment of bifurcated cerebral artery aneurysms and their associated ischemic complications, and to evaluate its potential value in understanding disease mechanisms and optimizing treatment protocols. In this study, a three-dimensional thermal model was used to simulate the heat conduction and blood flow characteristics of the bifurcated parts of cerebral arteries. Combined with the clinical data of patients, the influence of the thermal model on the formation mechanism of aneurysms and the prognosis evaluation of ischemic complications were analyzed through numerical calculation and experimental verification. Compare the performance of different treatment methods in thermal models. The thermal model shows that the temperature distribution at the aneurysm site is closely related to the blood flow velocity, and the local temperature increase may promote thrombosis and lead to ischemic complications. The treatment regimen optimized by the thermal model showed significant improvement in the simulation, reducing the risk of complications and improving hemodynamic parameters. Heat model can reveal the dynamic relationship between blood flow and temperature in the study of cerebral artery bifurcation aneurysm, and provide theoretical support for clinical treatment.

Fluid flow and heat transfer of a novel passive cooling system for gearless wind turbines with a power range of 3–12 MW

Today, the gearless horizontal axis wind turbines are mainstream in wind energy industry. High demands of electric power led to bigger systems and active cooling reduces the overall efficiency of the turbines. Passive cooling systems have been examined for the first time for a gearless wind energy generator with power range of 3–12 MW. With further developed heat conductors, it is possible to operate a wind generator in a larger power class with passive cooling components. This is accompanied by enormous cost savings due to elimination of costs for active cooling elements such as the use of fans, pumps, and heat exchangers. An additional factor is the significant reduction in necessary maintenance due to the minimized incidence of corrosion and wear. Design of the cooling fins and an ideal position of the generator within the housing of the wind turbine has been objectives. Mandatory is that maximum temperatures of the generator fins should stay under 155 °C, which could be reached with several designs for different heat exchanger geometries.

Improving a shell-tube latent heat thermal energy storage unit for building hot water demand using metal foam inserts at a constant pumping power

The phase transition heat transfer during the melting and solidification processes of phase change materials (PCMs) was modeled in a shell-tube thermal energy storage unit. For the first time, special attention was dedicated to the heat transfer enhancement on the heat transfer fluid (HTF) side by placing metal foam (MF) inserts in the HTF tube. The MF inserts on the HTF side were placed next to the fin bases to produce the maximum heat transfer enhancement. A two-temperature local thermal non-equilibrium model was applied to accurately model the transient heat transfer in the MF domains. To perform a fair heat transfer analysis and consider the hydrodynamic aspects of heat transfer, the pumping power for driving the HTF fluid was kept fixed. The finite element method with automatic time-step control was used to integrate the model equations. Results show that the insertion of the MF in the HTF tube enhances conduction heat transfer but reduces flow for a constant pump power. Using a small amount of MF provides a good flow rate and a convective heat transfer while using a large amount of MF provides good conduction heat transfer but a poor flow rate and convection. A moderate amount of MF inserts results in poor conduction and a poor flow rate, notably increasing the phase transition time. A 5 % MF insert could provide about a 13 % shorter solidification time compared to a 25 % MF insert. Thus, a well-designed and engineered MF insert configuration on the HTF side can reduce the phase transition time and improve the energy storage power.

Microwave-powered integrated carbon capture and utilisation (ICCU) for methane production with rapid temperature response from minimal thermal inertia

Integrated CO2 capture and utilisation (ICCU), emerging as an effective approach for achieving global carbon neutrality goals, demonstrates its potential through the transformation of captured CO2 into valuable methane and the concurrent hydrogen storage via methanation reaction. However, traditional ICCU strategies are heavily reliant on conventional thermal conduction heating, and the in-situ methanation process is constrained by reaction thermodynamics, posing significant challenges. In this study, we introduce for the first time a novel microwave-powered system, employing a clean and efficient energy source that offers rapid heating and on-demand operation and acts as an external field enhancement method in the ICCU-methanation procedure. The designed blend of Ni-loaded CaO and carbon nanotubes (CNTs) integrates the functions of CO2 adsorption, hydrogenation catalysis, and microwave absorbing ability. This novel approach achieved a CO2 capture capacity of 3.31 mmol gDFM-1 with CO2 conversion and CH4 selectivity reaching 71.01 % and 98.17 %, respectively, at a nickel loading of 20 wt%. We also propose a possible rate enhancement mechanism for the microwave-powered ICCU-methanation cyclic process compared to conventional methods. The microwave-powered system can substantially enhance overall efficiency by swiftly responding to temperature swings, thus reducing processing time in ICCU-methanation cycles.

Low Thermal Conductivity of Epidote and Its Cooling Effect on the Oceanic Crust

AbstractEpidote is a potential water carrier in subducting oceanic crust, and it is involved in temperature records from blueschist and eclogite exhumed from subduction zones. Thermal conductivity (κ) and thermal diffusivity (D) of epidote were measured at 303–1,473 K and 0.5–3.0 GPa using the transient plane‐source method. κ and D values decreased with temperature before dehydration, but increase anomalously with temperature after dehydration, and finally decreased with crystallization of the decomposition products. Thermal conductivity and thermal diffusivity of epidote exhibited a positive pressure dependence. The obtained κ and D values of epidote were significantly smaller than those of nominally anhydrous minerals and other hydrous minerals. X‐ray computed tomography analysis revealed that the pores inside the recovered samples were partially connected. κ and D values were dominated by conductive heat transport before dehydration, while the anomalous increase in the thermal parameters after dehydration may have been caused by the fluid inside the connected pores. Using variable thermal conductivity obtained in this study, the thermal structures of warm subducting slabs associated with epidote‐blueschist and epidote‐eclogite facies were modeled. Our estimations revealed that the temperatures of warm slabs below the forearc with the epidote‐blueschist and epidote‐eclogite facies may be lower 7–10 and ∼40–50 K compared to those determined using a constant thermal conductivity value. This implies that the presence of epidote‐blueschist and epidote‐eclogite could reduce the temperature of the slab, which may alter the dehydration depth of the slab.

Numerical simulations of conductive heating vacuum membrane distillation: Quantitative analysis of shunted heat flows and preventive strategy of salt crystallization

Numerical simulations of conductive heating vacuum membrane distillation: Quantitative analysis of shunted heat flows and preventive strategy of salt crystallization

Analytical modelling of transient conduction heat transfer in tubes for industrial applications

Analytical modelling of transient conduction heat transfer in tubes for industrial applications

Modeling and Simulation of the Aging Behavior of a Zinc Die Casting Alloy

While zinc die-casting alloy Zamak is widely used in vehicles and machines, its solidified state has yet to be thoroughly investigated experimentally or mathematically modeled. The material behavior is characterized by temperature and rate sensitivity, aging, and long-term influences under external loads. Thus, we model the thermo-mechanical behavior of Zamak in the solid state for a temperature range from −40 °C to 85 °C, and the aging state up to one year. The finite strain thermo-viscoplasticity model is derived from an extensive experimental campaign. This campaign involved tension, compression, and torsion tests at various temperatures and aging states. Furthermore, the thermo-physical properties of temperature- and aging-dependent heat capacity and heat conductivity are considered. One significant challenge is related to the multiplicative decompositions of the deformation gradient, which affects strain and stress measures relative to different intermediate configurations. The entire model is implemented into an implicit finite element program and validation examples at more complex parts are provided so that the predicability for complex parts is available, which has not been possible so far. Validation experiments using digital image correlation confirm the accuracy of the thermo-mechanically consistent constitutive equations for complex geometrical shapes. Moroever, validation measures are introduced and applied for a complex geometrical shape of a zinc die casting specimen. This provides a measure of the deformation state for complex components under real operating conditions.

Bivariate Jacobi polynomials depending on four parameters and their effect on solutions of time-fractional Burgers’ equations

The utilization of time-fractional Burgers’ equations is widespread, employed in modeling various phenomena such as heat conduction, acoustic wave propagation, gas turbulence, and the propagation of chaos in non-linear Markov processes. This study introduces a novel pseudo-operational collocation method, leveraging two-variable Jacobi polynomials. These polynomials are obtained through the Kronecker product of their one-variable counterparts, concerning both spatial (x) and temporal (t) domains. The study explores the impact of four parameters (θ,ϑ,σ,ς>−1) on the accuracy of resulting approximate solutions, marking the first examination of such influence. Collocation nodes in a tensor approach are constructed employing the roots of one-variable Jacobi polynomials of varying degrees in x and t. The study delves into analyzing how the distribution of these roots affects the outcomes. Consequently, pseudo-operational matrices are devised to integrate both integer and fractional orders, presenting a novel methodological advancement. By employing these matrices and appropriate approximations, the governing equations transform into an algebraic system, facilitating computational analysis. Furthermore, the existence and uniqueness of the equations under study are investigated and the study estimates error bounds within a Jacobi-weighted space for the obtained approximate solutions. Numerical simulations underscore the simplicity, applicability, and efficiency of the proposed matrix spectral scheme.

Застосування методу двобічних наближень до розв’язання першої крайової задачі для одновимірного рівняння теплопровідності з експоненціально нелінійним коефіцієнтом теплопровідності

The problem of heat conduction in nonlinear media reduces to solving boundary value problems for nonlinear heat conduction equations, where the coefficients of the equation or the function of the power of heat sources depend on temperature according to some law. Among the numerical methods for solving problems for nonlinear equations of mathematical physics, one can distinguish finite difference methods, finite element methods, variational and projection methods, as well as iterative methods. Among the latter group, the two-sided approximation method is particularly attractive due to its ability to provide a convenient estimate for the error of the approximate solution and to prove the existence of a solution to the original problem. The theory of linear partially ordered spaces was developed by L. V. Kantorovich in the second half of the 1930s. Further development of this theory is associated with the works of M. A. Krasnosel’skii, H. Amann, V. I. Opojtsev, N. S. Kurpel, B. A. Shuvar, A. I. Kolosov etc. The aim of this article is to develop a two-sided approximation method based on the use of Green's functions for solving the first boundary value problem for a one-dimensional nonlinear heat conduction equation and to investigate its performance when solving test problems. To achieve this goal, the unknown function was replaced, and the boundary value problem was reduced to a Hammerstein integral equation, which was considered as a nonlinear operator equation in a partially ordered Banach space. Conditions for the existence of a unique positive solution to the problem and conditions for two-sided convergence of successive approximations to it were obtained. The developed method was implemented in software and tested on solving test problems. The results of the computational experiment are illustrated with graphical and tabular information

An experimental setup for studying conduction and radiation heat transfer in a rotary drum

This study presents the development of a novel heating system that aims to bridge the experimental gap for analyzing both radiation and conduction heat transfer in a rotary drum. The experiments were conducted by loading the drum with 4 mm silica glass beads at varying fill levels (10 %, 17.5 %, and 25 %) and rotation rates (2, 5, and 10 RPM), and monitoring the temperature evolution with an infrared camera. Higher fill levels resulted in lower average particle bed temperatures, while rotation rate influenced the contact time between particles and the drum wall and the exposure to radiative heat. A rotation rate of 5 RPM resulted in the highest average bed temperature due to optimal contact and exposure time for conduction and absorption of radiative heat, respectively. Thus, a balance between adequate particle mixing and optimal contact/exposure time is crucial for maximizing heat transfer efficiency in rotary drums exhibiting conduction and radiation.

Recent Advances in High-Entropy Ceramics: Synthesis Methods, Properties, and Emerging Applications

High-entropy ceramics (HECs) represent an emerging class of materials composed of at least five different cations or anions in near-equiatomic proportions, garnering significant attention due to their extraordinary functional and structural properties. While multi-component ceramics have played a crucial role for many years, the concept of high-entropy materials was first introduced eighteen years ago with the synthesis of high-entropy alloys, and the first high-entropy nitride films were reported in 2014. These newly developed materials exhibit superior properties over traditional ceramics, such as enhanced thermal stability, hardness, and chemical resistance, making them suitable for a wide range of applications. High-entropy carbides, borides, oxides, oxi-carbides, oxi-borides, and other systems fall within the HEC category, typically occupying unique positions within phase diagrams that lead to novel properties. HECs are particularly well suited for high-temperature coatings, for tribological applications where low thermal conductivity and similar heat coefficients are critical, as well as for energy storage and dielectric uses. Computational tools like CALPHAD streamline the element selection process for designing HECs, while innovative, energy-efficient synthesis methods are being explored for producing dense specimens. This paper provides an in-depth analysis of the current state of the compositional design, the fabrication techniques, and the diverse applications of HECs, emphasizing their transformative potential in various industrial domains.

Bio-Based Phase Change Materials for Sustainable Development.

The increasing global population has intensified the demand for energy and food, leading to significant greenhouse gas (GHG) emissions from both sectors. To mitigate these impacts and achieve Sustainable Development Goals (SDGs), passive thermal storage methods, particularly using phase change materials (PCMs), have become crucial for enhancing energy efficiency and reducing GHG emissions across various industries. This paper discusses the state of the art of bio-based phase change materials (bio-PCMs), derived from animal fats and plant oils as sustainable alternatives to traditional paraffin-based PCMs, while addressing the challenges of developing bio-PCMs with suitable phase change properties for practical applications. A comprehensive process is proposed to convert bacon fats to bio-PCMs, which offer advantages such as non-toxicity, availability, cost-effectiveness, and stability, aligning with multiple SDGs. The synthesis process involves hydrolysis to break down fat molecules obtained from the extracted lipid, followed by three additional independent processes to further tune the phase change properties of PCMs. The esterification significantly decreases the phase transition temperatures while slightly improving latent heat; the UV-crosslinking moderately raises both the phase transition temperature and latent heat; the crystallization remarkably increases the both. The future research and guidelines are discussed to develop the large scale manufacturing with cost effectiveness, to optimize synthesis process by multiscale modeling, and to improve thermal conductivity and latent heat capacities at the same time.