Transient Heat Transfer Analysis Research Articles

Transient heat transfer analysis of airflow in a thermal water-bearing tunnel considering airflow turbulence and surrounding rock seepage effects

In deep tunnels traversing fractured water-rich strata, surface water descend, heats up at depth, and then rises into the tunnel, resulting in a ‘thermal water-bearing’ tunnel. This paper investigates the transient temperature evolution of such tunnels through laboratory-scale model experiment and numerical simulation, considering the convective heat transfer of thermal water in the surrounding rock. The experimental model was designed using similarity criteria derived from the equation analysis method. The 3D numerical model simulates the interaction between thermal water-saturated surrounding rock and tunnel airflow. The thermal field of the surrounding rock accounts for the coupling effects of the flow field, which follows Darcy's law. The tunnel airflow is governed by the RANS (Reynolds-Averaged Navier-Stokes) equations and the k−ω turbulence model, combined with the LRNM (Low-Reynolds Number Modeling) method for near-wall treatment. This approach is coupled with the energy balance equation to calculate the airflow temperature. The numerical model is validated using monitoring data from experimental measurements of grouted tunnel wall and airflow temperatures. Finally, parameter studies are conducted for both ungrouted and grouted tunnels, demonstrating the effectiveness of grouting in controlling airflow temperature by blocking thermal water.

Conservative Analysis of Transient Heat Transfer in Thermal Protection Systems with Interval Parameters

Conservative Analysis of Transient Heat Transfer in Thermal Protection Systems with Interval Parameters

Multi-scale modeling: Finite element analysis of the thermophysical properties of carbon/carbon composites considering manufacturing defects and porosity at high temperature

Carbon/carbon (C/C) composites have strict requirements for their thermal properties in extreme environments, especially at high temperature, where their properties are crucial for long-term life. The structure of C/C composites is exceptionally complex and exhibits multi-scale characteristics, and their thermophysical properties are closely related to temperature. Therefore, exploring their thermal response and heat transfer mechanisms through experimental method is relatively costly. This paper constructs a multi-scale finite element model to investigate the influence of structure and manufacturing defects on performance, and analyzes the thermophysical properties of the composites at High-temperature. By constructing a micro-representative volume element (Micro-RVE) that includes fibers, matrix, and pores, and using a steady-state heat transfer analysis method, the influence of material phase distribution on performance is studied. At the same time, a Meso-RVE reflecting the layering form, needle punching effect and manufacturing defects is established, and the transient heat transfer analysis method is used to investigate the influence of the structural form on performance. Finite element analysis shows that, the simulation values of the three C/C composites of unidirectional fiber bundles, short-chopped fiber felts, and needle punched C/C composites show good consistency with the experimental results in the published literature. This paper uses numerical methods to calculate the thermophysical properties of C/C composites from room temperature to 900°C and explores their variation with temperature. The constructed multi-scale model provides an accurate and effective method for predicting the thermophysical properties of C/C composites and also provides new perspectives and insights for thermal-mechanical coupling analysis and structural design.

Experimental analysis of transient and steady-state heat transfer from an impinging jet to a moving plate

Experimental analysis of transient and steady-state heat transfer from an impinging jet to a moving plate

Influence of Part Size on Microstructure during Austempering of Medium Carbon Low‐Alloy Steels

In this work, the influence of part dimensions on the microstructure of low‐alloy steel during its austempering has been studied. 0.41C‐0.24Si‐0.63Mn‐0.93Cr steel of different diameters ranging from 5 to 16 mm has been austempered at 350 °C for 10 min. Despite similar austempering conditions, an increase in the size and fraction of acicular/polygonal ferrite has been noted with an increase in the diameter of the samples. From the experimental observations and analytical calculations based on multi‐dimensional transient heat transfer analysis, the cooling rate during austempering has been observed to increase with a decrease in the sample diameter. Variation in ferrite grain size and fraction is due to different cooling rates encountered by different sized specimens during austempering. Close agreement between the experimental and theoretically estimated ferrite grain size due to variation in cooling rate during austempering has been noted.

Transient heat transfer analysis of a sandwich panel with a cracked honeycomb core

Sandwich structures with ceramic honeycomb cores are extensively employed in thermal protection systems owing to their exceptional ability to resist high temperatures. This work aims at exploring the effect of cracking on the transient thermal process of the sandwich panel subject to impulsive and cyclic thermal loadings. Both the conventional and re-entrant hexagonal alumina honeycombs are considered for the core material. By the integral transform method, combined with singular integral equations, the transient temperatures of the whole sandwich panel are determined from the semi-analytical solution. The straightforward temperature difference of the crack face’s midpoints is exploited to characterize the heat intensification near the crack. Parametric investigations are carried out for the internal cell angle, the relative density, crack length, crack position, and thickness of face sheets, which provides a better understanding of the honeycomb materials working in thermal protection systems.

Additive Manufacturing Process Simulation of Laser Powder Bed Fusion and Benchmarks

As the aerospace industry continues to adopt additively manufactured (AM) parts for flight hardware, process simulation becomes more attractive to improve the manufacturing process by understanding how process parameters affect part quality and performance. Process simulation can also be used to predict and prevent build failures before the printing process. The powder bed fusion process of Ti-6Al-4V material is modeled using commercial finite element software to simulate the selective laser melting process. The thermal history is obtained from transient heat transfer analysis. Both the inherent strain approach and a sequential thermal-mechanical approach are employed to predict residual stress and part distortion. A National Agency for Finite Element Methods and Standards (NAFEMS) benchmark problem is presented as a numerical example. It is a thin wall structure with geometry features that can lead to part defects due to thermal distortion. It is shown that both analysis approaches are able to capture the thin-member bridging behavior, stepping behavior, and general distortion contour plot as those published by NAFEMS.

Numerical Analysis of Transient Convective Heat Transfer in a Dental Implant

Hot substance consumption can have adverse effects on the neighbouring bone tissue in proximity to a dental implant. Elevated temperatures at the interface between the bone and implant could potentially disrupt the local cellular processes crucial for osteointegration. The primary goal of this study was to analyse the temperature and heat flux distributions within the implant body, surrounding bone, and bone-implant interface when the implant system subjected to a thermal load of transient nature. Transient thermal finite element analysis was utilised to analyse a three-dimensional model of dental implant with three different lengths – 6, 10, and 13 mm – placed in a mandible section. In order to obtain realistic results, thermal load was applied through convection on the outer surface of the prosthesis, simulating exposure to a hot liquid with the temperature and convection heat transfer coefficient of 67°C and 0.005 W/mm2°C, respectively. The temperature of the other components in the model was maintained at a constant 37°C. The results showed that increasing the implant length generally led to lower temperature and heat flux levels in the implant body, bone, and bone-implant interface. The highest temperature and heat flux values were concentrated in the superior region, gradually decreasing toward the inferior region. Importantly, all maximum temperature values remained below the limits associated with cellular bone necrosis and remodelling, thereby reducing the risk of osteoporosis. It is noteworthy that, when considering transient thermal load, shorter implants pose a significantly higher risk of implant failure compared to longer ones.

A new analytical model for transient heat transfer modeling of geothermal monotube heat exchangers

In this study, we present a new analytical model for the analysis of transient heat transfer (the model is solved after Laplace transformation) in geothermal monotube heat exchangers (GHE). The model accounts for various factors, including the heat exchanger's geometric parameters, the nature of the fluid (air and water), flow characteristics (laminar and turbulent), and heat transfer boundary conditions. Initially, the model's performance is validated by comparing its results with those obtained from computational fluid dynamics simulations using OpenFOAM. This validation process ensures the reliability and accuracy of the model. Additionally, an analysis of the heat transfer mechanisms and key modeling parameters identified from the comparison is provided. In the subsequent section, the analysis is extended by comparing the proposed model to four existing analytical models available in the literature. These reference models employ distinct approaches to simulate heat transfer to the ground. By examining a range of cases, the performance of the model is systematically evaluated against these conventional methods. This work closes by an analysis of the analytical models and the impact on the accuracy of the results highlighting the advantages of the proposed method in geothermal applications.

Transient Heat Transfer Analysis in Metal Plates with Variable Thickness

Nonlinear transient heat transfer via conduction–radiation is a dynamic topic of long-standing interest with applications ranging from aeronautical and mechanical engineering to industrial and civil security. To gain a better understanding of the performances of materials having thermal proprieties that change during nonlinear heat transfer, several studies using the finite element method (FEM) have been conducted. Such studies apply nonlinear thermal material characteristics to describe the complete system under different loading conditions in each region by adjusting the temperature values for the other three edges and the thickness parameter with Dirichlet boundary conditions. As a result, while modeling and simulating temperature distributions for such situations, nonlinearities generated by temperature-dependent thermal conductivity must be considered. In this work, we focus on the analysis of coupled transient heat transfer through two metal plates with temperature-dependent thermal characteristics in which the temperature is fixed along the bottom edge and heat is transferred from both the top and bottom faces of the two plates. FEM is employed to solve the nonlinear heat equation and compute the temperature as a function of time for variable thickness. The study examines the effect of modifying the thickness parameter values on the temperature distribution over time for various edge values over 5,000 s.

LES investigations of turbulent heat transfer structures with exponential power escalation

The BORAX-type accident is a key scenario in the safety design of pool-type experimental nuclear reactors. It considers a reactivity insertion large enough to initiate an exponential power escalation with a period as small as tens of milliseconds, necessitating effective heat removal by the coolant system. This study presents a computational analysis of single-phase transient heat transfer under thermal-hydraulic conditions relevant to such reactors, focusing on turbulent channel flow between planar fuel plates with highly subcooled water. Our investigation uses Large Eddy Simulation (LES) performed with TrioCFD, the open-source CFD code developed by the French Atomic Energy Commission (CEA). Several modeling options are tested: incompressible flow with temperature as a passive scalar, incompressible flow with varying viscosity, and quasicompressible flow with temperature-dependent viscosity and mass density. To evaluate the accuracy of our physical models, all simulations were performed under uniform conditions, using the same experimental parameters, mesh, and numerical strategies. The specific examined scenario involves a highly subcooled turbulent channel flow at moderate pressure, experiencing an exponential power increase. We ensure detailed turbulence resolution at the Batchelor scale in the wall normal direction close to the heating wall. This communication presents a comparative study of LES simulations of the single-phase flow under exponential power excursion at the wall. Subsequent analyses focused on the impact of different velocity-temperature couplings on the prediction of wall heat transfer, incorporating comparisons with experimental infrared thermometry measurements. Early in the transient phase, temperature peaks are confined within turbulent streaks, as reported in the existing literature. However, as the exponential power transient unfolds, the viscous layer and streak structures destabilize, leading to a more dispersed distribution of hot spots. Preliminary findings suggest a transition in heat transfer modes from turbulence-driven to conduction during the final stages of the single-phase transient.

A methodology for the optimal placement of conformal cooling channels in injection molds: 2D transient heat transfer analysis

A methodology for the optimal placement of conformal cooling channels in injection molds: 2D transient heat transfer analysis

Numerical investigation and optimisation of solid–solid phase change material composite-based plate-fin heat sink for thermal management of electronic package

Phase Change Material based Thermal Management Systems show promise in electronic packaging systems for safe operation, acting as a thermal storage element to prevent sudden temperature surges, enhancing operational time and mitigating electronic system failures. Therefore, the meticulous design of a solid–solid phase change material (SSPCM) based heat sink with fins is imperative for effective electronic thermal management by improving its safe operational time. This study performs a numerical analysis of two-dimensional transient heat transfer for the SSPCM-based plate-fin heat sink using the modified enthalpy porosity method, where the plate fins serve as Thermal Conductivity Enhancers (TCEs). Neopentyl glycol with 6% (w/w) Copper oxide nanoparticle is used as the SSPCM composite. Enhancement ratio based on safe operational time at different power levels (2–6 W) with different percentages of TCE and fin count recorded variation from 1.32 to 13.7. The TCE integration is more effective at higher power levels. The optimisation of SSPCM-based heat sink design for maximum safe operational time is accomplished through Response Surface Methodology. The significance of the TCE fraction, number of fins, heat generation rate, and interaction effects over the safe operational time are observed. The correlation for safe operational time as a function of TCE volume fraction, fin count and heat generation rate is developed. TCE volume fraction possesses the highest influencing factor over safe operational time at lower and higher power levels compared to fin count. At 2 W and 6 W, a maximum change in safe operational time is observed from 72.7% at 2 W to 83.45% at 6 W. The heat sink with 0.262 volume fraction of TCE and 11 fins is the optimised design that can provide the desirable safe operational time for the tested range of power levels.

Thermal boundary conditions for heat transfer analysis of bridges considering non-uniform distribution of internal air temperature by computational fluid dynamics

Heat transfer analysis has been used to calculate the temperature distribution in bridges. Thermal boundary conditions play a critical role in this analysis. However, existing studies on thermal boundary conditions simplify the air temperature inside the bridge deck as uniform, which is not realistic and thus causes inaccurate simulation results. This study proposes a new approach to thermal boundary conditions in the heat transfer analysis of bridges. For the first time, computational fluid dynamics is used to calculate non-uniform air temperatures inside the bridge deck. In addition, non-approximate heat exchange equations for long-wave radiation are also incorporated into the approach. The techniques are applied to the 1377-m main span Tsing Ma Suspension Bridge to calculate the internal air temperatures of a deck segment. Transient heat transfer analysis is then conducted to calculate the time-dependent temperature distribution of the segment. As compared with the field monitoring results, the proposed approach can simulate the temperature distribution of the bridge with an average discrepancy of 0.88 °C and is more accurately than other existing approaches.

Transient magnetohydrodynamic heat and mass transfer analysis of chemically reacting fluid flow over oscillatory permeable media

This study focuses on the numerical observation of the convective motion of chemically active magnetohydrodynamic (MHD) fluid through a vertically oriented permeable medium, incorporating variations in mass and heat transfer. The fluid type is assumed to be incompressible, chemically strongly ionized and viscous with some mass infusibility. The model associated with this problem is solved by a highly stable Implicit Finite Difference Method (IFDM). The method is used for small and large deflection of the physical parameters, which results in a noticeable fluid flow behavior. Numerical configuration is graphically depicted to scrutinize the fluid behavior. The momentum, energy, concentration diffusion, skin friction, Nusselt number and Sherwood number are investigated for numerous factors such as magnetic field, permeability and chemical reaction rate. The current study unveils significant findings, demonstrating that a heightened rate of chemical reaction in the presence of magnetic effects, coupled with specific porosity, diminishes ionization energy, resulting in a concurrent decrease in the concentration and momentum profiles of the fluid flow. The rise in the viscous diffusion rate is attributed to escalating values of the Schmidt number, causing an augmentation in dynamic viscosity and consequently resulting in an overall reduction in the momentum of the fluid flow.

Numerical Modeling and Analysis of Transient and Three-Dimensional Heat Transfer in 3D Printing via Fused-Deposition Modeling (FDM)

Fused-Deposition Modeling (FDM) is a commonly used 3D printing method for rapid prototyping and the fabrication of plastic components. The history of temperature variation during the FDM process plays a crucial role in the degree of bonding between layers. This study presents research on the thermal analysis of the 3D printing process using a developed simulation code. The code employs numerical discretization methods with an implicit scheme and an effective heat transfer coefficient for cooling. The computational model is validated by comparing the results with analytical solutions, demonstrating an agreement of more than 99%. The code is then utilized to perform thermal analyses for the 3D printing process. Interlayer and intralayer reheating effects, sensitivity to printing parameters, and realistic printing patterns are investigated. It is shown that concentric and zigzag paths yield similar peaks at different time intervals. Nodal temperatures can fall below the glass transition temperature (Tg) during the printing process, especially at the outer nodes of the domain and under conditions where the cooling period is longer and the printed volume per unit time is smaller. The article suggests future work to calculate welding time at different conditions and locations for the estimation of the degree of bonding.

Prediction of latent heat storage transient thermal performance for integrated solar combined cycle using machine learning techniques

The objective of this work is to develop a prediction methodology for latent heat storage thermal performance to reduce computational efforts. Thermal storage can effectively compensate for the limitations of solar energy by extending operation time and leveling fluctuations. Although latent heat storage is advantageous compared to sensible heat storage for its large energy density and stable output temperature, its thermal performance has strong nonlinearity. With the conventional numerical simulations, it takes large computational time to obtain the one year or longer thermal performance of the latent heat storage. In this study, the exit steam enthalpy of latent heat storage for an integrated solar combined cycle (ISCC) is predicted using machine learning techniques. As latent heat storage is used to store solar thermal energy, the inlet steam properties such as enthalpy, pressure, and flow rate are continuously altered. Heat charged-discharged history in the regression model inputs significantly improves the prediction accuracy. The regression models are trained with transient heat transfer analysis results based on meteorological data of Bawean Island, Indonesia. Seven regression algorithms (Gaussian process regression (GPR), kernel approximation, ensemble of trees, support vector machines (SVM), artificial neural networks (ANN), linear regression, and regression trees) are compared and GPR is selected for its root mean square error (RMSE) of 16.9 kJ/kg or less. The RMSE values are equivalent to Bawean results even when the GPR is applied to two other locations: Kupang and Cikarang. Meanwhile, an ensemble classifier, which is selected from nine classification algorithms (Ensemble classifiers, ANN, SVM, decision trees, nearest neighbor classifiers, naive bayes classifiers, kernel approximation classifiers, discriminant analysis, and logistic regression classifiers), is used to determine when to terminate the latent heat storage operation. Owing to its short calculation time, the developed methodology is advantageous at the early design stages of latent heat storage, especially when multiple locations need to be compared.

Preliminary analysis of maximum hypothetical accident for a solid core space nuclear reactor power system

In order to study the safety characteristics of the solid core space nuclear reactor power (SNRP) system under the maximum hypothetical accident, a two-dimensional entire core transient heat transfer analysis model was established, and the key parameters response characteristics of the ultra-small lithium-cooled SNRP under two maximum hypothetical accidents, namely, loss of heat sink (LOHS) and loss of coolant accident (LOCA), were calculated and analyzed. In the maximum hypothetical accidents, the coolant cooling capacity is lost, thus the core decay power is discharged into space only through the radiation of the residual heat removal system and the side surface of the reactor vessel. During the accidents, the heat transfer of the core deteriorated, and the fuel temperature may rise to the melting point, resulting in radioactive leakage. The results show that: (1) In the LOHS accident, the maximum fuel temperature reaches 2150 K at 550 s, and the pressure of the primary system volume accumulator continues to increase to the set system pressure safety limit of 2 MPa, resulting the primary loop overpressure failure. And the fuel matrix temperature is close to the set cladding limit temperature of 2200 K; (2) In the LOCA, the deterioration of heat transfer in the core makes the temperature increase rapidly, reaching a maximum of 3016 K at about 630 s, which is very close to the melting temperature of UN fuel 3123 K. As the decay power decreases, the maximum core temperature decreases to less than 1600 K after 24 h of the accident. The auxiliary cooling system of the solid core SNRP system under the maximum hypothetical accident is optimized, and the design parameters of the auxiliary cooling system meeting the safety requirements are obtained.

Nonlinear transient heat transfer analysis of multilayered thermal protection systems by incremental differential quadrature element method

As a first attempt, the incremental differential quadrature element method (IDQEM) is used for the one-dimensional nonlinear transient heat transfer analysis of multilayered thermal protection systems (TPSs) under time-dependent aerodynamic heating. All thermophysical properties are considered to be temperature-dependent (TD). Based on the IDQEM, the multilayered TPS is divided into several spatial sub-domains along the layer interfaces, and the entire heating process is also divided into several temporal sub-domains. For each temporal sub-domain, the governing equations are discretized by the differential quadrature method and then solved by the Newton-Raphson method. The initial condition of each temporal sub-domain is defined by the temperature results at the end of the previous sub-domain. The temperature profile of the TPS during the entire heating process can be obtained by repeating the calculation procedure from the first temporal sub-domain to the last one. Convergence and accuracy of the method are confirmed. Numerical examples are carried out to show the effects of TD thermophysical properties, thermal boundary, and layer thickness on the temperature profile of a three-layer TPS. The results show that the temperature dependency of thermophysical properties is worthy of consideration in predicting the temperature profile of the TPSs.

Transient heat transfer study of small-caliber barrel in continuous fire

In this paper, the barrel of small-caliber automatic rifle is taken as the research object. An experiment was taken on the gun barrel and an infrared thermal imager was used to measure the change of barrel temperature field during the continuous shooting, and the finite difference model of barrel heat transfer was established. The formula of the correction coefficient of heat transfer of gunpowder gas on the inner surface of the body tube was obtained by fitting with the experimental results. The barrel transient heat transfer analysis calculation software with independent intellectual property rights was developed by C++, and the characteristics of temperature field were analyzed. The experimental results show that the temperature of the outer surface of the barrel rises and then falls along the axial direction from the bottom of the barrel to the muzzle, and the rising slope will decrease in the next 30 rounds of continuous firing. The calculation results of the small-caliber barrel with chromium layer show that the temperature of the outer surface of the chromium layer of each shot rises instantaneously and then drops rapidly, in the form of pulse, and the temperature of the steel matrix rises continuously, but it is significantly lower than the peak temperature of the chromium layer, and the presence of the chromium layer provides buffer protection for the barrel matrix material. There is still an influence region in the matrix behind the interface between the chromium layer and the barrel matrix, and the temperature of the steel matrix in this region fluctuates. The thickness of the influence region is determined by analysis. The temperature of the inner surface of the barrel reaches a maximum at the tail of the barrel along the axial axis and then slowly decreases in the direction of the muzzle, the radial temperature distribution of the barrel is steeper at the tail end and the muzzle area is relatively flat. The calculation results of the model in this paper are basically consistent with the experimental results, which provides a rapid calculation model and tool for the thermal design and analysis of the gun barrel.