Inverse Heat Conduction Problem Research Articles

The inverse problem of heat conduction in the case of non-uniqueness: A functional identification approach

The inverse problem of heat conduction in the case of non-uniqueness: A functional identification approach

Identification of an N‐valued heterogeneous conductivity profile in an inverse heat conduction problem

Identification of an <i>N</i>‐valued heterogeneous conductivity profile in an inverse heat conduction problem

An accelerated sequential function specification method based on LM gradient for transient inverse heat conduction problem

The sequential function specification method (SFSM) and the Levenberg-Marquardt method (LM) are two typical methods for inverse heat conduction problems (IHCP). However, these two methods have poor convergence and low convergence speed. Therefore, the combination of LM with SFSM is depended on to propose a modified Sequential Levenberg-Marquardt gradient method (SLM). The computation results of three shapes of heat fluxes with noise show the SLM is averagely improved over 16.3% more than SFSM in computation accuracy, but equally reduced over 39.9% less than in computation time. On other hand, when the emphasis is on the effect of the Aitken acceleration method, Tikhonov regularization and polynomial order on SLM, a neural network prediction model of the root mean square error is trained by regularization factor, polynomial order and the number of future time steps. Subsequently, the genetic algorithm optimization method is proposed for the minimum root mean square error so that it can improve the computation accuracy and efficiency of IHCP to the maximum extent. Finally, Second-order polynomial order and small order-of-magnitude regularization factor are highly recommended for parameter selection with multiple choices. The ASLM only needs to determine the optimal number of future time steps.

Noninvasive reconstruction of internal heat source in biological tissue using adaptive simulated annealing algorithm

The heat distribution information of human lesions is of great value for disease analysis, diagnosis, and treatment. It is a typical inverse problem of heat conduction that deriving the distribution of internal heat sources from the temperature distribution on the body surface. This paper transforms such an inverse problem of bio-heat transfer into a direct one, thereby avoiding complex boundary conditions and regularization processes. To noninvasively reconstruct the internal heat source and its corresponding 3D temperature field in biological tissue, the adaptive simulated annealing (ASA) algorithm is used in the simulation module, where the position P(x, y, z) of point heat source in biological tissue and its corresponding temperature T are set as the optimization variables. Under a certain optimized sample, one can obtain the simulated temperature distributing on the surface of the module, then subtract the simulated temperature from the measured temperature of the same surface which was measured using a thermal infrared imager. If the sum of absolute values of the difference is smaller, it indicates that the current sample is closer to the true location and temperature of the heat source. When the values of optimization variables are determined, the corresponding 3D temperature field is also confirmed. The simulation results show the simulated position and temperature of the heat source are very approximate with those of the real experimental module. The method presented in this paper has enormous potential and promising prospects in clinical research and application, such as tumor hyperthermia, disease thermal diagnosis technology, etc.

A Bayesian Spatiotemporal Modeling Approach to the Inverse Heat Conduction Problem

Abstract This study introduces a Bayesian spatiotemporal modeling approach to solve inverse heat conduction problems (IHCPs), employing penalized splines within a spatiotemporal forward model. The complexity and ill-posed nature of IHCPs, characterized by potential nonexistence, nonuniqueness, or instability of solutions, pose significant challenges for traditional methods. Addressing this, our study presents a spatiotemporal forward model that incorporates spatial, temporal, and interaction terms, accurately capturing the intricate dynamics inherent in IHCPs and using this information as a leverage to solve the inverse problem. We adopted a Bayesian inference framework for the subsequent parameter estimation problem and developed a Gibbs sampling algorithm to sample from the posterior distribution of the model's parameters, enhancing the estimation process. Through case studies on a one-dimensional (1D) heat simulation and a pool boiling experiment using multisensor thermocouple data for heat flux reconstruction, we demonstrate the model's superiority over traditional methods. The inclusion of the spatiotemporal interaction term significantly enhances model performance, indicating its potential for broader application in solving IHCPs. The application of this method in both simulated and real-world scenarios highlights its effectiveness in capturing the spatiotemporal complexities of IHCPs. This work contributes to the field by offering a robust methodology for addressing the spatial and temporal complexities inherent in IHCPs, supported by a comprehensive Bayesian inference framework and the use of a Gibbs sampling algorithm for parameter estimation.

Genetic Algorithm as the Solution of Non-Linear Inverse Heat Conduction Problems: A Novel Sequential Approach

Abstract Direct measurement of surface heat flux could be extremely challenging, if not impossible, in numerous applications. In such cases, the use of temperature measurement data from subsurface locations can facilitate the determination of surface heat flux and temperature through the solution of the inverse heat conduction problem (IHCP). Several different techniques have been developed over the years for solving IHCPs, with different levels of complexity and accuracy. The filter coefficient technique has proved to be a promising approach for solving IHCPs. Inspired by the filter coefficient approach, a novel method is presented in this paper for solving one-dimensional IHCPs in a domain with temperature-dependent material properties. A test case is developed in COMSOL Multiphysics where the front side of a slab is subject to known transient heat flux and the temperature distributions within the domain are calculated through numerical simulation. The IHCP solution in the form of filter coefficients is defined and a genetic algorithm (GA) is used for the calculation of filter matrix. The number of significant filter coefficients required to evaluate surface heat flux at each time-step is determined through trial and error and the optimal number is selected for evaluating the solution. The structure of the filter matrix is assessed and discussed. The resulting filter coefficients are then used to evaluate the surface heat flux for several different test cases and the performance of the proposed approach is assessed in detail. The results showed that the presented approach is robust and can result in finding optimal filter coefficients that can accurately estimate various types of surface heat flux profiles using temperature data from a limited number of time steps and in a near real-time fashion.

Restoration of unsteady heat flow from a thermal energy accumulator by solving the inverse heat conduction problem

Restoration of unsteady heat flow from a thermal energy accumulator by solving the inverse heat conduction problem

Oil and Gas Structures: Forecasting the Fire Resistance of Steel Structures with Fire Protection under Hydrocarbon Fire Conditions

The hydrocarbon temperature–time curve is widely used instead of the standard curve to describe the temperature in the environment of structural surfaces exposed to fire in oil and gas chemical facilities and tunnels. This paper presents calculations of the ratio of time to reach critical temperatures at different nominal fire curves for steel structures such as bulkheads and columns with different types of fireproofing. The thermophysical properties of the fireproofing materials were obtained by solving the inverse heat conduction problem using computer simulation. It was found that the time interval for reaching critical temperatures in structures with different types of fireproofing in a hydrocarbon fire decreased, on average, by a factor of 1.2–1.7 compared to the results of standard fire tests. For example, for decks and bulkheads with mineral wool fireproofing, the K-factor of the ratio of the time for reaching the critical temperature of steel under the standard curve to the hydrocarbon curve was 1.30–1.62; for plaster, it was 1.56; for cement boards, it was 1.34; for non-combustible coatings, it was 1.38–2.0; and, for epoxy paints, it was 1.71. The recommended values of the K-factor for fire resistance up to 180 min (incl.) were 1.7 and, after 180 min, 1.2. The obtained dependencies would allow fireproofing manufacturers to predict the insulation thickness for expensive hydrocarbon fire experiments if the results of fire tests under standard (cellulosic) conditions are known.

An Inverse Problem for Estimating Spatially and Temporally Dependent Surface Heat Flux with Thermography Techniques

In this study, an inverse conjugate heat transfer problem is examined to estimate temporally and spatially the dependent unknown surface heat flux using thermography techniques with a thermal camera in a three-dimensional domain. Thermography techniques encompass a broad set of methods and procedures used for capturing and analyzing thermal data, while thermal cameras are specific tools used within those techniques to capture thermal images. In the present study, the interface conditions of the plate and air domains are obtained using perfect thermal contact conditions, and therefore we define the problem studied as an inverse conjugate heat transfer problem. Achieving the simultaneous solution of the continuity, Navier–Stokes, and energy equations within the air domain, alongside the heat conduction equation in the plate domain, presents a more intricate challenge compared to conventional inverse heat conduction problems. The validity of our inverse solutions was verified through numerical simulations, considering various inlet air velocities and plate thicknesses. Notably, it was found that due to the singularity of the gradient of the cost function at the final time point, the estimated results near the final time must be discarded, and exact measurements consistently produce accurate boundary heat fluxes under thin-plate conditions, with air velocity exhibiting no significant impact on the estimates. Additionally, an analysis of measurement errors and their influence on the inverse solutions was conducted. The numerical results conclusively demonstrated that the maximum error when estimating heat flux consistently remained below 3% and higher measurement noise resulted in the accuracy of the heat flux estimation decreasing. This underscores the inherent challenges associated with inverse problems and highlights the importance of obtaining accurate measurement data in the problem domain.

Parametric investigation of flow and heat transfer on the reflooding quenching of a cylinder rod in forced convection

A reflooding quenching boiling experiment in forced convection was conducted to investigate the flow and heat transfer performance of a FeCrAl cylindrical rod under various subcooling degrees and flow rates. Employing image processing and inverse heat conduction problem (IHCP) methods, the research provides insights into vapor film thickness, quenching front propagation, wall temperature cooling rate and minimum film boiling temperature. Key observations include the presence of numerous extremely small bubbles dispersing around the rod and a third quenching front appears in the middle of the rod under higher subcooling condition at 25 °C, attributed to a thinner vapor film, smaller bubbles resulting from film rupture and the inertia of fluid flow. Within the experimental range, the impact of the subcooling degree on the fluid flow and heat transfer performance during the quenching boiling process is more pronounced than that of the coolant flow rate. The evolution of the vapor film indicates a substantial reduction of 78.76 % in film thickness and a significant increase of 185.71 % in propagation velocity with an increase in subcooling degree, while an elevation in flow rate lead to a 29.84 % decrease in film thickness and a 45 % increase in propagation velocity. Heat transfer performance benefits from increased subcooling degrees throughout the entire quenching boiling process, whereas flow rate impacts the heat transfer performance of film boiling alone. Furthermore, by incorporating the impact of flow velocity, a predictive formula for Tmin is established with an error margin within ±5 %.

Determination of boundary conditions from the solution of the inverse heat conduction problem in the gas nitriding process

Two non-linear unsteady methods for solving the inverse heat conduction problem are discussed and compared in this paper. The first method is based on the finite element method, uses the variational calculus and the BFGS (Broyden-Fletcher–Goldfarb–Shanno) optimization algorithm. The second is an analytical and numerical method based on the approximation of the solution in the form of a linear combination of Chebyshev polynomials. In numerical tests, the stability of the temperature, heat flux density, and heat transfer coefficient obtained from both methods was analyzed. On the basis of experimental data, oscillations of the measured gas temperature and temperature in the component were analyzed. Oscillations of the gas temperature obtained from the experiment were taken into consideration during tests. The first method to solve the inverse problem was used to determine the boundary conditions for the entire gas nitriding process. As a result, stable values of temperature, heat flux density, and the heat transfer coefficient on the surface of the component under treatment in a real gas nitriding process, so far unpublished, were obtained.

Determination of transient heat transfer by cooling channel in high-pressure die casting using inverse method

Complex shapes of aluminum castings are typically manufactured during the short cycle process known as the high-pressure die casting (HPDC). High productivity is ensured by introducing die cooling through a system of channels, die inserts or jet coolers. Die cooling can also effectively help in reducing internal porosity in cast components. Accurate simulations based on sophisticated numerical models require accurate input data such as material properties, initial and boundary conditions. Although the heat is dominantly dissipated through die cooling, indicating the importance of knowing precise thermal boundary conditions, open literature lacks a detailed information about the spatial distribution of heat transfer coefficient. This study presents an inverse method to determine accurate heat transfer coefficients of a die insert based on temperature measurements in multiple points by 0.5 mm K-type thermocouples and a subsequent solution of the two-dimensional inverse heat conduction problem. The solver was built in the open-source CFD code OpenFOAM and the free library for nonlinear optimization NLopt. The results are presented for the commonly used 10 mm die insert with a hemispherical tip and coolant flow rates ranging from 100 l/h to 200 l/h. Heat transfer coefficients reach values well above 50 kW/m2K in the hemispherical tip, which is followed by a secondary peak and then a gradual drop to values around 1 kW/m2K further downstream.

Reconstruction of internal heat source in biological tissue using parallel particle swarm optimization

With the rapid development of engineering thermophysics, researches on human biological heat transfer phenomena has gradually shifted from qualitative to quantitative. The heat distribution information of human lesions is of great value for disease analysis, diagnosis, and treatment. It is a typical inverse problem of heat conduction that deriving the distribution of internal heat sources from the temperature distribution on the body surface. Differing from traditional numerical methods for solving heat conduction, this paper transforms such an inverse problem of bio-heat transfer into a direct one, thereby avoiding complex boundary conditions and regularization processes. To noninvasively reconstruct the internal heat source and its corresponding 3D temperature field in biological tissue, the parallel particle swarm optimization (PPSO) is used in the simulation module, where the position P(x, y, z) of point heat source in biological tissue and its corresponding temperature T are set as the optimization variables. Under a certain optimized sample, one can obtain the simulated temperature distributing on the surface of the module, then subtract the simulated temperature from the measured temperature of the same surface which was measured using a thermal infrared imager. If the absolute value of the difference is smaller, it indicates that the current sample is closer to the true location and temperature of the heat source. When the values of optimization variables are determined, the corresponding 3D temperature field is also confirmed. The simulation results show the simulated position and temperature of the heat source are very approximate with those of the real experimental module. The method presented in this paper has enormous potential and promising prospects in clinical research and application, such as tumor hyperthermia, disease thermal diagnosis technology, etc.

Numerical Analysis of Convective Heat Transfer in Quenching Treatments of Boron Steel under Different Configurations of Immersed Water Jets and Its Effects on Microstructure

In this work, the effects of jet impact angle and water flow on the heat-transfer coefficient in boron steel probes were analyzed. Angles of 90°, 75° and 60° were used with stirring flows of 33 l·min−1, 25 l·min−1, 13 l·min−1 and 6 l·min−1. The aim consisted of determining the heat-extraction rates by analyzing the correlation programmed in the Ansys Fluent 2020R2 software when different cooling conditions are used, avoiding many experiments, and establishing quenching conditions free of surface defects on the workpiece. This process is currently used in heavy machinery, requiring high hardness and wear resistance. The fluid-dynamic field was validated using a scale physical model using the particle image velocimetry technique, PIV. In contrast, the thermal field was validated with transient state experiments solving the inverse heat conduction problem, IHCP. The results show that for high flows (33 l·min−1), the jets with an angle of 90° impact the entire surface of the piece, but their cooling rate is slower compared to the other angles, being 243.61 K·s−1, and 271.70 K·s−1, 329.56 K·s−1 for 75° and 60°, respectively. However, for low flows (6 l·min−1), the impact velocities are very similar for the three cases, promoting more homogeneous cooling rates of 58.47 K·s−1, 73.58 K·s−1 and 63.98 K s−1 for angles of 90°, 75° and 60°, respectively. Likewise, through the use of CCT diagrams, it was determined that regardless of the cooling rate, the final structure will always be a mixture of martensite–bainite due to the effect of boron as determined experimentally, which implies a more significant proportion of martensite at higher cooling rates.

Analysis of the effect of surface coating and texturization on the estimation of effective thermophysical properties

This article aims to investigate the influence of surface modifications on the effective thermophysical properties of textured carbon steel and coated hard metal. The novelty of this work lies in the simultaneous estimation of thermal diffusivity and thermal conductivity via measurements using one active heating surface and Bayesian inference to solve the inverse heat conduction problem. To achieve this, an experimental setup was developed where samples were partially heated on an active surface, and the temperature was measured at two different points on the surface to estimate the thermal properties. The method was applied to two different thermal models using the same set of experimental data. Thus, the thermophysical properties were determined with and without the Markov Chain Monte Carlo technique using the Metropolis-Hastings sampling algorithm. Both techniques proved suitable for estimating the thermophysical properties. It is noteworthy that texturization and coating did not significantly modify the effective thermal diffusivity of the analyzed samples. However, there were considerable changes in the effective thermal conductivity, with an increase of approximately 13 % in the textured carbon steel sample and a reduction of approximately 11 % in the coated hard metal sample. The expanded uncertainty of the estimated effective thermophysical properties was less than 14 %. These results reflect low dispersion and good accuracy of the experimental technique used.

Band Heater Heat Flux Characterization Using Inverse Heat Conduction Problem Models

Abstract An experimental apparatus was constructed to correlate water flowrate and temperature rise under an external band heater. Due to the physical characteristics of the band heater, its transient heating behavior is unknown. This paper investigates the application of inverse heat conduction problem (IHCP) methods to characterize the heat flux from the band heater. Three experiments with different heating times (5, 10, and 20 s) and no flowrate were conducted to measure the transient temperature under the 400 W band heater. Type-T thermocouples measure surface temperature at the centerline of the band heater. The experimental results are computed with five different heat conduction models. The models are chosen to identify how the heat flux response varies from a simplified to a realistic model. Additionally, the results of the experimental heat flux are compared to the manufacturer band heater data (58.9 kW/m2) for each model. The minimum time needed for the heater to fully energize the system is from 10 to 12 s. The residuals for each model are analyzed and used to evaluate the appropriateness of the five different models. The results show that the use of simpler models can achieve results similar to those of complex models, with less time and computational cost.

Investigation of operational limit of a pulsating heat pipe by estimating local heat transfer

How pulsating heat pipes (PHPs) reach their operational limit has not yet been fully understood. This study aims to provide a complete picture of the termination mechanism of the self-oscillation of vapor and liquid. Experimental studies on a 10-turn PHP with HFC-134a were conducted and the filling ratio (FR) was from 20% to 80%. The thermo-fluid behavior in the PHP was investigated by temperature measurements with a high-resolution and high-speed infrared camera and estimation of fluid-to-wall heat flux distributions by solving inverse heat conduction problems. The results suggested that the PHP, increasing heat load, reached the operational limit due to different mechanisms depending on the filling ratio: at a high FR (80%), the liquid volume ratio increased with the increase of the operating temperature, resulting in the compressed liquid phase. At a low FR (20%), when a large amount of heat was applied, the fluid in the evaporator dried out and became a superheated vapor. The PHP with an optimum FR (50%) transferred the maximum heat under the same evaporator temperature, as the fluid in the PHP was able to keep the saturated two-phase state until the evaporator temperature exceeded the critical temperature.

Fire Protection of Steel Structures of Oil and Gas Facilities: Multilayer, Removable, Non-Combustible Covers

Fire protection is required to protect metal structures of oil and gas facilities from fires. Such fire protection should provide high fire resistance limits: 60, 90, 120 and more minutes. Specialists of LLC “RPC PROMIZOL ” developed a multilayer, removable type of fire protection made of superfine basalt fibre and ceramic materials for operation in Arctic conditions. Five experimental studies were carried out in standard and hydrocarbon fire regimes. The fire protection effectiveness of the products for I20 beams without load was obtained: a 50 mm thick coating provided 130 min of a standard fire regime; a 15 mm thick coating provided 60 min. The 15 mm thick coating provided 30 min of a hydrocarbon fire regime and the 50 mm thick coating provided 93 min of a hydrocarbon fire regime. The I40 beam under a load of 19.9 tf showed an R243 for the standard fire regime. The coefficients of effective thermal conductivity and specific heat capacity of fire-retardant compositions were determined by solving the inverse heat conduction problem. The problem was solved by modelling using the QuickField 7.0 software package, which implements FEM. Modelling showed that for obtaining the fire resistance limit R120 under the standard fire regime for the sample steel structure from an I40 beam, it is enough to apply fire protection with a thickness of 25 mm instead of 50 mm, which agrees with the experimental data. For the hydrocarbon regime, it is predicted that R120 can be obtained at a thickness of 45 mm instead of 50 mm.

Conjugate gradient method with regularization in estimating mold surface heat flux during continuous casting

The reconstruction of the boundary heat flux between solidifying steel and water-cooled mold wall using measured temperature data is recognized as an inverse heat conduction problem (IHCP). A Tikhonov spatial regularization approach to enhance the smoothness of the Conjugate Gradient Method (CGMr) was proposed which incorporates three different orders of spatial regularization: zeroth, first, and second order. The optimal regularization parameter was selected using a modified L-curve method. The accuracy of the CGMr is investigated with Case 1: a triangular time-spatial variation heat flux, and Case 2: a step change in the form of rectangular variation heat flux. The effects of measurement noise, grid points and time-step size are investigated. The results show that the minimum relative error (eRMS) of the predicted Case 1 heat flux is 1.98%, 7.64%, and 8.02% for zeroth-, first-, and second-order spatial regularization, respectively. The corresponding values for the predicted Case 2 heat flux are 3.82%, 13.42%, and 14.91%. Subsequently, CGMr algorithm is applied to calculate the heat flux in a mold simulator experiment. By comparing the relationship between heat fluxes reconstructed by CGMr after different iteration numbers, it is observed that the recovered heat flux of 2.14 MW/m2 with zeroth regularization remains highly stable.

Determination of interfacial heat transfer coefficient at the frozen sand mold casting process of ZL101 alloy

Compared to the resin sand mold casting process, frozen casting is more environmentally friendly, providing a better working environment and enhanced supercooling degree. The interfacial heat transfer coefficient (IHTC) between frozen sand mold and metal is an important parameter that significantly influences the final mechanical properties and microstructure of the castings. This paper solved the inverse heat conduction problem using the finite difference method (FDM). In addition, the conjugate gradient method (CGM) was adopted to calculate the temperature distribution and heat flux in the molten metal. At the same time, the particle swarm optimization algorithm (PSO) was used in temperature distribution determination in frozen sand mold. The interfacial heat transfer coefficient (IHTC) was estimated during the solidification of ZL101. The results showed a good agreement between calculated and experimental data, obtaining accurate casting interface temperature Tm, frozen sand mold interface temperature Ts, heat flux q, and IHTC. The analysis of the IHTC variation revealed a water content value within the range of 4 wt.% to 5 wt.% resulted in IHTC in two types of interpretation, called ‘fluctuation type’ and ‘turning type’.