Inverse estimation of spatially and temporally varying heating boundary conditions of a two-dimensional object

In many dynamic heat transfer situations, the temperature at the heated boundary is not directly measurable and can be obtained by solving an inverse heat conduction problem (IHCP) based on measured temperature or/and heat flux at the accessible boundary. In this study, IHCP in a two-dimensional rectangular object is solved by using the conjugate gradient method (CGM) with temperature and heat flux measured at the boundary opposite to the heated boundary. The inverse problem is formulated in such a way that the heat flux at heated boundary is chosen as the unknown function to be recovered, and the temperature at the heated boundary is computed as a byproduct of the IHCP solution. The measurement data, i.e., the temperature and heat flux at the opposite boundary, are obtained by numerically solving a direct problem where the heated boundary of the object is subjected to spatially and temporally varying heat flux. The robustness of the formulated IHCP algorithm is tested for different profiles of heat fluxes along with different random errors of the measured heat flux at the opposite boundary. The effects of the uncertainties of the thermophysical properties and back-surface temperature measurement on inverse solutions are also examined.     

Two-Dimensional Inverse Heat Transfer Analysis of Functionally Graded Materials in Estimating Time-Dependent Surface Heat Flux

Two-dimensional transient inverse heat conduction problem (IHCP) of functionally graded materials (FGMs) is studied herein. A combination of the finite element (FE) and differential quadrature (DQ) methods as a simple, accurate, and efficient numerical method for FGMs transient heat transfer analysis is employed for solving the direct problem. In order to estimate the unknown boundary heat flux in solving the inverse problem, conjugate gradient method (CGM) in conjunction with adjoint problem is used. The results obtained show good accuracy for the estimation of boundary heat fluxes. The effects of measurement errors on the inverse solutions are also discussed.     

A SERIAL SOLUTION FOR THE 2-D INVERSE HEAT CONDUCTION PROBLEM FOR ESTIMATING MULTIPLE HEAT FLUX COMPONENTS

A serial algorithm for the inverse heat conduction problem (IHCP) has been developed to estimate the individual flux components, one by one, at the unknown boundary, based on the function specification method. The sensitivity coefficient defined in this algorithm brings out the influence of the heat flux components independent of each other. The objective function minimizes the difference in the measured temperature and the contribution of the individual flux component to the thermal field at the sensor location. The serial algorithm developed here could be used with data from both overspecified and underspecified sensors with respect to the number of flux components. The method was tested for delineating independent heat fluxes at the boundary of a two-dimensional solid for both space- and time-varying heat fluxes. Simulated thermal histories obtained from direct solution were used as inputs for the inverse problem for characterizing the new algorithm.

Three types of analyses were done on the results of the IHCP, focused on (1) the convergence of error in estimated temperatures at the different sensor locations, (2) overall error in estimated temperatures for the whole domain, and (3) the total heat energy transferred across the boundary. It is shown that the optimum configuration of independent unknown fluxes is given by the one with minimum energy estimates across the boundary, for both cases.     

A two-dimensional inverse heat conduction problem in estimating the fluid temperature in a pipeline

An inverse heat conduction problem (IHCP) was investigated in the two-dimensional section of a pipe elbow with thermal stratification to estimate the unknown transient fluid temperatures near the inner wall of the pipeline. An inverse algorithm based on the conjugate gradient method (CGM) was proposed to solve the IHCP using temperature measurements on the outer wall. In order to examine the accuracy of estimations, some comparisons have been made in this case. The temperatures obtained from the solution of the direct heat conduction problem (DHCP) using the finite element method (FEM) were pseudo-experimental input data on the outer wall for the IHCP. Comparisons of the estimated fluid temperatures with experimental fluid temperatures near the inner wall showed that the IHCP could accurately capture the actual temperature in form of the frequency of the temperature fluctuations. The analysis also showed that the IHCP needed at least 13 measurement points for the average absolute error to be dramatically reduced for the present IHCP with 37 nodes on each half of the pipe wall.     

Two-dimensional non-linear inverse heat conduction problem based on the singular value decomposition

In this paper an efficient sequential method is developed in order to estimate the unknown boundary condition on the surface of a body from transient temperature measurements inside the solid. This numerical approach for solving an inverse heat conduction problem (IHCP) takes into account two-dimensional problems, planar or axisymmetric cylindrical, composite materials with irregular boundaries and temperature-dependent thermal properties. The unknown surface condition is assumed to have abrupt changes at unknown times. The regularization procedure used for the solution of the IHCP is based on the singular value decomposition technique. An overall estimate of error is defined in order to find the optimal estimation in the 2D IHCP (linear and non-linear). The stability and accuracy of the scheme presented is evaluated by comparison with the Function Specification Method. This comparative study has been carried out using numerically simulated data, and the parameters considered include shape of input, noise level of measurement, size of time step and temperature-dependent thermal properties. A good agreement was found between both methods. Beside this, the slight differences on estimations and number of future temperatures are discussed in this paper.     

The Estimation of Fluid Temperatures Through an Inverse Heat Conduction Technique

Transient, internal temperatures within an instrumented probe are considered as part of an inverse heat conduction problem (IHCP) to compute the temperature of the surrounding fluid. A linear scheme is used where the exchange coefficients are treated as known parameters.Input data to the IHCP have been generated numerically. When these are uncorrupted, the inverse algorithm works well without stabilization. However, in practice the algorithm must be stabilized, as it is shown that noise is amplified substantially. It becomes necessary both to parameterize spatial variations in the fluid temperature and to utilize a functional specification method to address the noncausal solution.     

Boundary reconstruction in two-dimensional steady state anisotropic heat conduction using a regularized meshless method

We study the stable numerical identification of an unknown portion of the boundary on which either a Dirichlet or a Robin boundary condition is provided, while additional Cauchy data are given on the remaining known part of the boundary of a two-dimensional domain, in the case of steady state anisotropic heat conduction problems. This inverse geometric problem is solved using the method of fundamental solutions (MFS) in conjunction with the Tikhonov regularization method [53]. The optimal value for the regularization parameter is chosen according to Hansen’s L-curve criterion [17]. The stability, convergence, accuracy and efficiency of the proposed method are investigated by considering several examples in both smooth and piecewise smooth geometries.     

Solution of inverse heat conduction problem using the lattice Boltzmann method

In this paper the D2Q9 lattice Boltzmann method (LBM) was utilized for the solution of a two-dimensional inverse heat conduction (IHCP) problem. The accuracy of the LBM results was validated against those obtained from prevalent numerical methods using a common benchmark problem. The conjugate gradient method was used in order to estimate the heat flux test case. A complete error analysis was performed. As the LBM is attuned to parallel computations, its use is recommended in conjugation with IHCP solution methods.     

Two dimensional steady-state inverse heat conduction problems

In this work we estimate the surface temperature in two dimensional steady-state in a rectangular region by two different methods, the singular value decomposition (SVD) with boundary element method (BEM) and the least-squares approach with integral transform method (ITM). The BEM method is efficient for solving inverse heat conduction problems (IHCP) because only the boundary of the region needs to be discretized. Furthermore, both temperature and heat flux at the unknown boundary are estimated at the same time. The least-squares technique involves solving the equations constructed from the measured temperature and the exact solution. The measured data are simulated by adding random errors to the exact solution of the direct problem. The effects of random errors on the accuracy of the predictions are examined. The sensitivity coefficients are also presented to illustrate the effect of sensor location on the estimated surface conditions. Numerical experiments are given to demonstrate the accuracy of the present approaches.     

Heat flux estimation of a plasma rocket helicon source by solution of the inverse heat conduction problem

A solution of the inverse heat conduction problem (IHCP) by the steepest descent method is carried out in order to determine the waste heat flux from a helicon plasma discharge using transient surface temperature measurements obtained from infrared thermography. The infrared camera data is calibrated against thermocouple data and mapped to real locations on the observed surface. The magnitude and distribution of the heat flux to the gas containment tube in the helicon is investigated as the applied power, gas flow rate, magnetic field distribution and neutral gas are varied.     

A NONLINEAR INDIRECT MEASUREMENT PROBLEM FOR A MULTITHERMOCOUPLE PROBE IMMERSED IN A LIQUID

The problem described herein concerns the processing of the time-dependent, internal temperatures within a multithermocouple probe. These are used to compute the temperature of the surrounding fluid, as part of an inverse heat conduction problem (IHCP). The novel achievement in this work is that the exchange coefficients do not have to be supplied a priori, but instead are an additional solution output. Consequently the IHCP is nonlinear and requires significant stabilization. Four methods are applied successively, until a satisfactory solution is found: the parameterization of spatial variations in fluid temperatures and exchange coefficients; a functional specification method (using future time data) to address the noncausal nature of the solution; a lower bound on the exchange coefficient; and a maximum number of iterations at each time step (in accordance with the discrepancy principle).     

Reconstruction of local heat fluxes in pool boiling experiments along the entire boiling curve from high resolution transient temperature measurements

In this paper, we consider a transient inverse heat conduction problem (IHCP) defined on an irregular three-dimensional (3D) domain in pool boiling experiments. Heat input to a circular copper heater of 35 mm diameter and 7 mm thickness is provided by a resistance heating foil pressed to the bottom of the heater. The heat flux at the inaccessible boiling side is estimated from a number of temperature readings in the heater volume. These temperatures are measured by some high-resolution microthermocouples, which are mounted 3.6 μm below the surface in the test heater. The IHCP is formulated as a mathematical optimization problem and solved by the conjugate gradient (CG) method. The arising partial differential equations (PDEs) are solved using the software package DROPS. A simulation case study is used to validate the performance of the solution approach. Finally, we apply the solution approach to the IHCP in pool boiling experiments. The procedure enables the reconstruction of local instantaneous heat flux distribution on the heater surface at different locations along the boiling curve.     

Thermocouple Data in the Inverse Heat Conduction Problem

The presence of thermocouples inside a heat-conducting body will distort the temperature field in the body and may lead to significant bias in the temperature measurement. If temperature histories obtained from thermocouples are used in the inverse heat conduction problem (IHCP), errors are propagated into the IHCP results. The bias in the thermocouple measurements can be removed through use of appropriate detailed thermocouple models to account for the dynamics of the sensor measurement. The results of these models can be used to generate correction kernels to eliminate bias in the thermocouple reading, or can be applied as sensitivity coefficients in the IHCP directly. Three-dimensional and axisymmetric models are compared and contrasted and a simple sensitivity study is conducted to evaluate the significance of thermal property selection on the temperature correction and subsequent heat flux estimation. In this paper, a high-fidelity thermocouple model is used to account for thermocouple bias in an experiment to measure heat fluxes from solidifying aluminum to a sand mold. Correction kernels are obtained that are used to demonstrate the magnitude of the temperature measurement bias created by the thermocouples. The corrected temperatures are used in the IHCP to compute the surface heat flux. A comparison to IHCP results using uncorrected temperatures shows the impact of the bias correction on the computed heat fluxes.     

Identification of boundary heat fluxes in a falling film experiment using high resolution temperature measurements

In this paper, we consider a three-dimensional inverse heat conduction problem (IHCP) in a falling film experiment. The wavy film is heated electrically by a thin constantan foil and the temperature on the back side of this foil is measured by high resolution infrared (IR) thermography. The transient heat flux at the inaccessible film side of the foil is determined from the IR data and the electrical heating power. The IHCP is formulated as a mathematical optimization problem, which is solved with the conjugate gradient method. In each step of the iterative process two direct transient heat conduction problems must be solved. We apply a one step θ-method and piecewise linear finite elements on a tetrahedral grid for the time and space discretization, respectively. The resulting large sparse system of equations is solved with a preconditioned Krylov subspace method. We give results of simulated experiments, which illustrate the performance and tuning of the solution method, and finally present the estimation results from temperature measurements obtained during falling film experiments.     

Inverse problem on the identification of temperature and thermal stresses in an FGM hollow cylinder by the surface displacements

A problem on the identification of time-dependent temperature on one of the limiting surfaces of a radially inhomogeneous hollow cylinder is formulated and solved under the temperature and radial displacement given on the other limiting surface. The analysis of temperature and thermal stress distribution in the cylinder is performed. The solution has been constructed by the reduction to an inverse thermoelasticity problem. By making use of the finite difference method, a stable solution algorithm is suggested for the analysis of inverse problem. The solution technique is verified numerically by making use of the solution to a relevant direct problem. It is shown that the proposed technique can be e?ciently used for the identification of a heat flux or unknown parameters (the surrounding temperature or the heat-exchange coe?cient) in the third-kind boundary conditions.     

Inverse heat conduction estimation of inner wall temperature fluctuations under turbulent penetration

Turbulent penetration can occur when hot and cold fluids mix in a horizontal T-junction pipe at nuclear plants. Caused by the unstable turbulent penetration, temperature fluctuations with large amplitude and high frequency can lead to time-varying wall thermal stress and even thermal fatigue on the inner wall. Numerous cases, however, exist where inner wall temperatures cannot be measured and only outer wall temperature measurements are feasible. Therefore, it is one of the popular research areas in nuclear science and engineering to estimate temperature fluctuations on the inner wall from measurements of outer wall temperatures without damaging the structure of the pipe. In this study, both the one-dimensional(1D) and the two-dimensional(2D) inverse heat conduction problem(IHCP) were solved to estimate the temperature fluctuations on the inner wall. First, numerical models of both the 1D and the 2D direct heat conduction problem(DHCP) were structured in MATLAB, based on the finite difference method with an implicit scheme. Second, both the 1D IHCP and the 2D IHCP were solved by the steepest descent method(SDM), and the DHCP results of temperatures on the outer wall were used to estimate the temperature fluctuations on the inner wall. Third, we compared the temperature fluctuations on the inner wall estimated by the 1D IHCP with those estimated by the 2D IHCP in four cases:(1) when the maximum disturbance of temperature of fluid inside the pipe was 3℃,(2) when the maximum disturbance of temperature of fluid inside the pipe was 30℃,(3) when the maximum disturbance of temperature of fluid inside the pipe was 160℃, and(4) when the fluid temperatures inside the pipe were random from 50℃ to 210℃.     

An inversion approach for the inverse heat conduction problems

The inverse heat conduction problems (IHCP) analysis method provides an efficient approach for estimating the thermophysical properties of materials, the boundary conditions, or the initial conditions. Successful applications of the IHCP method depend mainly on the efficiency of the inversion algorithms. In this paper, a generalized objective functional, which has been developed using a generalized stabilizing functional and a combinational estimation that integrates the advantages of the least trimmed squares (LTS) estimation and the M-estimation, is proposed. The objective functional unifies the regularized M-estimation, the regularized least squares (LS) estimation, the regularized LTS estimation, the regularized combinational estimation of the LTS estimation and the M-estimation, and the regularized combinational estimation of the LS estimation and the M-estimation into a concise formula. The filled function method, which is coupled with the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm, is developed for searching a possible global optimal solution. Numerical simulations are implemented to evaluate the feasibility and effectiveness of the proposed algorithm. Favorable numerical performances and satisfactory results are observed, which indicates that the proposed algorithm is successful in solving the IHCP.     

A data method using the inner temperature difference to improve the stability of the inverse heat conduction problems

Abstract

This article proposes a method to construct two series systems for improving the stability of the inverse heat conduction problems (IHCP) in a finite slab. The transfer function between the surface heat flux or temperature and the inner temperature difference is respectively obtained by Laplace transform technique firstly. Then the series systems which can solve IHCP based on the inner temperature difference are constructed by replacing the unsuitable zero and pole points of the transfer function approximated by è approximation. Finally the effects of the series systems are evaluated by a typical example. The results of the evaluation show that this method can obtain the surface heat flux and temperature by the inner temperature difference, and enhance the response speed of the measurement system at the same time. In addition this method can also improve the signal to noise ratio (SNR) of the inverse solutions by selectively amplifying the high SNR parts of the inner temperature difference. The present work provides an effective method to improve the stability of IHCP.     

Effect of pressure on heat transfer coefficient at the metal/mold interface of A356 aluminum alloy

The aim of this paper is to correlate interfacial heat transfer coefficient (IHTC) to applied external pressure, in which IHTC at the interface between A356 aluminum alloy and metallic mold during the solidification of casting under different pressures were obtained using the inverse heat conduction problem (IHCP) method. The method covers the expedient of comparing theoretical and experimental thermal histories. Temperature profiles obtained from thermocouples were used in a finite difference heat flow program to estimate the transient heat transfer coefficients. The new simple formula was presented for correlation between external pressure and heat transfer coefficient. Acceptable agreement with data in literature shows the accuracy of the proposed formula.     

Numerically solving twofold ill-posed inverse problems of heat equation by the adjoint Trefftz method

The inverse problem endowing with multiple unknown functions gradually becomes an important topic in the field of numerical heat transfer, and one fundamental problem is how to use limited minimal data to solve the inverse problem. With this in mind, in the present article we search the solution of a general inverse heat conduction problem when two boundary data on the space-time boundary are missing and recover two unknown temperature functions with the help of a few extra measurements of temperature data polluted by random noise. This twofold ill-posed inverse heat conduction problem is more difficult than the backward heat conduction problem and the sideways heat conduction problem, both with one unknown function to be recovered. Based on a stable adjoint Trefftz method, we develop a global boundary integral equation method, which together with the compatibility conditions and some measured data can be used to retrieve two unknown temperature functions. Several numerical examples demonstrate that the present method is effective and stable, even for those of strongly ill-posed ones under quite large noises.