Inverse analysis of forced convection in micro-channels with slip flow via integral transforms and Bayesian inference

The present work addresses the direct and inverse problems for convective heat transfer with incompressible laminar gas flow in micro-channels, within the range of validity of the slip-flow regime. The direct problem analysis combines the classical integral transform method and the generalized integral transform technique (GITT), by analytically solving the two-dimensional steady-state convection problem and finding a hybrid numerical-analytical solution for the required eigenvalue problem. The inverse problem analysis makes use of the accuracy and robustness of the direct problem solution and focus on the simultaneous identification of the momentum and thermal accommodation coefficients, related to gas flow and heat transfer within micro-channels, besides the usually unknown boundary condition parameters, here represented by the external Biot number. The inverse analysis is based on the availability solely of temperature measurements at the channel external wall, along its length, as obtained for instance via infrared camera thermography. A Bayesian inference approach is adopted in the solution of the identification problem based on the Monte Carlo Markov Chain method (MCMC) and the Metropolis-Hastings sampling algorithm. A typical example of slip flow in parallel-plates micro-channel is selected to illustrate both the direct and inverse problems solution approaches.     

Inverse Problem of a Piezothermoelastic Cylinder Subject to Transient Heating

An inverse problem of a radially polarized piezoelectric hollow circular cylinder of crystal class 6 mm is investigated. It is assumed that a voltage induced by the action of a time-varying temperature applied to the inner surface of the cylinder is measured on the outer surface. The inverse problem entails a determination of the heating temperature from knowledge of the measured voltage. First, an exact solution to the problem is found by solving the equations of equilibrium and electrostatics for the cylinder subject to the prescribed boundary conditions. Then a least-squares residual method that incorporates Lagrange multipliers for satisfaction of the boundary equations is employed to derive an approximate analytic solution. Both formulations are utilized in order to calculate the unknown heating temperature, and the corresponding temperature, displacement and stress fields in a cylinder of cadmium selenide. Numerical results based upon the least-squares residual formulation are found to compare favorably with those obtained by the exact analysis.     

A Markov Chain Monte Carlo-Metropolis Hastings Approach for the Simultaneous Estimation of Heat Generation and Heat Transfer Coefficient from a Teflon Cylinder

This paper reports the use of Markov Chain Monte Carlo (MCMC) and Metropolis Hastings (MH) approach, to solve an inverse heat transfer problem. Three-dimensional, steady state, conjugate heat transfer from a Teflon cylinder of dimensions 100 mm diameter and 100 mm length with uniform volumetric internal heat generation is considered. The goal is to estimate volumetric heat generation and heat transfer coefficient, given the temperature data at certain fixed location on the surface of the cylinder. The internal volumetric heat generation is specified as input and the temperature and heat transfer coefficient values are obtained by a numerical solution to the governing equation. The temperature values also depend on heat transfer coefficient which is obtained by solving Navier–Stokes equation to obtain flow information. In order to reduce the computational cost, a neural network is trained from the computational fluid dynamics simulations. This is posed as an inverse problem wherein volumetric heat generation and heat transfer coefficient are unknown but the temperature data is known by conducting experiments. The novelty of the paper is the simultaneous determination of volumetric heat generation and heat transfer coefficient for the experimentally measured steady-state temperatures from a Teflon cylinder using MCMC-MH as an inverse model in a Bayesian framework and finally, the estimates are reported in terms of mean, maximum a posteriori, and the standard deviation which is the uncertainty associated with the estimated parameters.     

GENERALIZED THERMOELASTIC PROBLEM FOR AN INFINITELY LONG HOLLOW CYLINDER FOR SHORT TIMES

In this work, we consider the one-dimensional problem for an infinitely long hollow cylinder in the context of the theory of generalized thermoelasticity with one relaxation time. The Laplace transform technique is used to solve the problem. The solution in the transformed domain is obtained by a direct approach. The inverse transforms are obtained in an approximate analytical manner using asymptotic expansions valid for small values of time. The temperature, displacement, and stress are computed and represented graphically.     

Application of optimization techniques and the enthalpy method to solve a 3D-inverse problem during a TIG welding process

Tungsten inert Gas (TIG) welding takes place in an atmosphere of inert gas and uses a tungsten electrode. In this process heat input identification is a complex task and represents an important role in the optimization of the welding process. The technique used to estimate the heat flux is based on solution of an inverse three-dimensional transient heat conduction model with moving heat sources. The thermal fields at any region of the plate or at any instant are determined from the estimation of the heat rate delivered to the workpiece. The direct problem is solved by an implicit finite difference method. The system of linear algebraic equations is solved by Successive Over Relaxation method (SOR) and the inverse problem is solved using the Golden Section technique. The golden section technique minimizes an error square function based on the difference of theoretical and experimental temperature. The temperature measurements are obtained using thermocouples at accessible regions of the workpiece surface while the theoretical temperatures are calculated from the 3D transient thermal model.     

Optimization of two-dimensional nonstationary thermal stresses and displacements in a half-space through the use of internal heat sources

A solution to the problem on optimal control for either a displacement-vector or a stress-tensor component in a section of an elastic half-space is presented for the case of plane strain. The power of internal heat sources located in a plane, which is parallel to the boundary surface, is selected to be the control function. Assuming the existence of a control function, which provides an exact lower limit for the uniform deviation of the controlled stress or displacement from their desired or given values, the optimization problem is reduced to an inverse thermoelasticity problem governed by an integral equation of the first kind. The latter equation has been reduced to Volterra integral equations of the first and second kind, whose solutions were found approximately by making use of the linear spline technique.     

Application of the method of fundamental solutions and radial basis functions for inverse heat source problem in case of steady-state

The paper deals with the inverse determination of heat sources in steady 2-D heat conduction problem. The problem is described by Poisson equation in which the function of the right hand side is unknown. The identification of the strength of a heat source is given by using the boundary condition and a known value of temperature in chosen points placed inside the domain. For the solution of the inverse problem of identification of the heat source the method of fundamental solution with radial basis functions is proposed. The accurate results have been obtained for five test problems where the analytical solutions were available.     

Practical application of inverse heat conduction for wall condition estimation on a rotating cylinder

The solution of the linear, inverse, transient heat conduction problem (IHCP) in a cylindrical geometry is analysed. The rotating cylinder under investigation is experiencing boiling convection induced by the impingement of a water jet. The initial temperature is known, additional temperature measurements in time are taken with sensors positioned at a constant radius within the solid material, and the estimation of the wall heat flux at the external radius is sought. First, simulated temperature measurements inside the cylinder are processed in order to be used to estimate the wall heat flux. When noise is present in the data, some of the simulated results obtained using the least squares method exhibit oscillatory behavior, but these large oscillations are substantially reduced by the implementation of a regularization technique. Real experimental data are also used for the wall condition estimation and for the subsequent building of local boiling curves are plotted and discussed. The question of the possible effect of a temperature dependent conductivity on the reconstructed wall condition is also considered.     

Inverse problem of coupled thermoelasticity for prediction of heat flux and thermal stresses in an annular cylinder

This study presents a means of solving the inverse problem of coupled thermoelasticity in an annular cylinder. While knowing the temperature history at any point of the body, the boundary time-varying heat flux can be computed, and subsequently the distributions of temperature and thermal stresses in the cylinder can be determined as well. The accuracy of the inverse analysis is examined by using simulated exact and inexact measurements obtained on the surface and at an interior location of the cylinder. Numerical results demonstrate that excellent estimation for the heat flux and thermal stresses distributions can be obtained for all the test cases considered in this study.     

Thermoelastic interaction in functionally graded thick hollow cylinder with temperature-dependent properties

The prediction of thermoelastic behavior induced by transient thermal shock is important to evaluate the durability of functionally graded materials. The purpose of this article is to study the axisymmetric thermoelastic interaction in a functionally graded thick hollow cylinder by an asymptotic approach. The governing equations with variable material properties, which are spatially graded and temperature dependent, are proposed based on the generalized theory of thermoelasticity with one relaxation time (L–S theory). The Laplace transform technique is used to derive the general solutions with the cylinder divided into thin cylinders and material properties assumed constant in each thin cylinder. The inverse Laplace transform is then conducted analytically by some approximations in the time domain, and the short-time solution of the problem with its interior boundary subjected to a sudden temperature rise and the outer surface maintained at constant temperature are obtained. Utilizing these asymptotic solutions, the propagation of thermal and thermoelastic waves are studied, which display dependence of each wave’s propagation upon the relaxation time, volume fraction parameter and temperature. The distributions of the radial displacement, temperature and stresses are also plotted and discussed. These results reveal effects of these variable material properties with spatial position and temperature on thermoelastic behavior.     

INVERSE FORCED CONVECTION PROBLEM OF SIMULTANEOUS ESTIMATION OF TWO BOUNDARY HEAT FLUXES IN IRREGULARLY SHAPED CHANNELS

This article deals with the use of the conjugate gradient method of function estimation for the simultaneous identification of two unknown boundary heat fluxes in channels with laminar flows. The irregularly shaped channel in the physical domain is transformed into a parallel plate channel in the computational domain by using an elliptic scheme of numerical grid generation. The direct problem, as well as the auxiliary problems and the gradient equations, required for the solution of the inverse problem with the conjugate gradient method are formulated in terms of generalized boundary-fitted coordinates. Therefore, the solution approach presented here can be readily applied to forced convection boundary inverse problems in channels of any shape. Direct and auxiliary problems are solved with finite volumes. The numerical solution for the direct problem is validated by comparing the results obtained here with benchmark solutions for smoothly expanding channels. Simulated temperature measurements containing random errors are used in the inverse analysis for strict cases involving functional forms with discontinuities and sharp corners for the unknown functions. The estimation of three different types of inverse problems are addressed in the paper: (i) time-dependent heat fluxes; (ii) spatially dependent heat fluxes; and (iii) time and spatially dependent heat fluxes.     

THREE-DIMENSIONAL TRANSIENT THERMAL STRESS ANALYSIS OF A NONHOMOGENEOUS HOLLOW CIRCULAR CYLINDER DUE TO A MOVING HEAT SOURCE IN THE AXIAL DIRECTION

In this paper, a theoretical analysis of a three-dimensional transient thermal stress problem is developed for a nonhomogeneous hollow circular cylinder due to a moving heat source in the axial direction from the inner and /or outer surfaces. Assuming that the hollow circular cylinder has nonhomogeneous thermal and mechanical material properties in the radial direction, the heat conduction problem and the associated thermoelastic behaviors for such nonhomogeneous medium are developed by introducing the theory of laminated composites as one of theoretical approximation. The transient heat conduction problem is treated with the help of the methods of Fourier cosine transformation and Laplace transformation, and the associated thermoelastic field is analyzed making use of the thermoelastic displacement potential, Michell's function, and the Boussinesq's function. Some numerical results for the temperature change and the stress distributions are shown in figures, and the effect of relaxing the thermal stress in the nonhomogeneous hollow circular cylinder and the influence of the velocity of a moving heat source are briefly discussed     

Fractional Order Thermoelastic Problem for an Infinitely Long Solid Circular Cylinder

In this study, we consider an infinitely long solid circular cylinder, whose lateral surface is traction free and subject to a known surrounding temperature. The problem is in the context of the fractional order thermoelasticity theory. The medium is assumed to be initially quiescent. The solution is obtained by a direct approach without the customary use of potential functions. Laplace transform technique is used to obtain the general solution for any set of boundary conditions. The inversion process is carried out using a numerical method based on Fourier series expansions. Numerical results are given and represented graphically.     

Inverse analysis with integral transformed temperature fields: Identification of thermophysical properties in heterogeneous media

The objective of this work is to introduce the use of integral transformed temperature measured data for the solution of inverse heat transfer problems, instead of the common local transient temperature measurements. The proposed approach is capable of significantly compressing the measured data through the integral transformation, without losing the information contained in the measurements and required for the solution of the inverse problem. The data compression is of special interest for modern measurement techniques, such as the infrared thermography, that allows for fine spatial resolutions and large frequencies, possibly resulting on a very large amount of measured data. In order to critically address the use of integral transformed measurements, we examine in this paper the simultaneous estimation of spatially variable thermal conductivity and thermal diffusivity in one-dimensional heat conduction within heterogeneous media. The direct problem solution is analytically obtained via integral transforms and the related eigenvalue problem is solved by the Generalized Integral Transform Technique (GITT). The inverse problem is handled with Bayesian inference by employing a Markov Chain Monte Carlo (MCMC) method. The unknown functions appearing in the formulation are expanded in terms of eigenfunctions as well, so that the unknown parameters become the corresponding series coefficients. Such projection of the functions in an infinite dimensional space onto a parametric space of finite dimension also permits that several quantities appearing in the solution of the direct problem be analytically computed. Simulated measurements are used in the inverse analysis; they are assumed to be additive, uncorrelated, normally distributed, with zero means and known covariances. Both Gaussian and non-informative uniform distributions are used as priors for demonstrating the robustness of the estimation procedure.     

THERMOELASTIC TRANSIENT RESPONSE OF MULTILAYERED HOLLOW CYLINDER WITH INITIAL INTERFACE PRESSURE

This article deals with the transient response of one-dimensional axisymmetric quasi-static coupled thermoelastic problems with initial interface pressure. The initial interface pressure in a multilayered cylinder caused by the heat-assembling method is considered as an initial condition for the thermoelastic equilibrium problem. The Laplace transform and finite difference methods are used to analyze problems. Using the Laplace transform with respect to time, the general solutions of the governing equations are obtained in the transform domain. The solution is obtained using the matrix similarity transformation and inverse Laplace transform. We obtained solutions for the temperature and thermal stress distributions in a transient state. Moreover, the computational procedures established in this article can solve the generalized thermoelasticity problem for a multilayered hollow cylinder.     

Estimation of the time-dependent heat flux using the temperature distribution at a point by conjugate gradient method

In this paper, the conjugate gradient method coupled with adjoint problem is used in order to solve the inverse heat conduction problem and estimation of the time-dependent heat flux using the temperature distribution at a point. Also, the effects of noisy data and position of measured temperature on final solution are studied. The numerical solution of the governing equations is obtained by employing a finite-difference technique. For solving this problem the general coordinate method is used. We solve the inverse heat conduction problem of estimating the transient heat flux, applied on part of the boundary of an irregular region. The irregular region in the physical domain (r,z) is transformed into a rectangle in the computational domain (ξ,η). The present formulation is general and can be applied to the solution of boundary inverse heat conduction problems over any region that can be mapped into a rectangle. The obtained results for few selected examples show the good accuracy of the presented method. Also the solutions have good stability even if the input data includes noise and that the results are nearly independent of sensor position.     

Estimation of heat flux on the surface of an initially hot cylinder cooled by a laminar confined impinging jet

In this study, an inverse algorithm based on the conjugate gradient method and the discrepancy principle is applied to estimate the unknown space-and time-dependent heat flux at the surface of an initially hot cylinder cooled by a laminar confined slot impinging jet from the knowledge of temperature measurements taken on the cylinder’s surface. It is assumed that no prior information is available on the functional form of the unknown heat flux; hence the procedure is classified as the function estimation in inverse calculation. The temperature data obtained from the direct problem are used to simulate the temperature measurements, and the effect of the errors in these measurements upon the precision of the estimated results is also considered. The results show that an excellent estimation on the space-and time-dependent heat flux can be obtained even the distributions of thermal properties inside the cylinder is unknown.     

Transient coupled magnetothermoelastic problem of multilayered hollow cylinder subjected to magnetic and vapor fields

This paper deals three-dimensional axisymmetric quasistatic-coupled magnetothermoelastic problems for time-dependent boundary condition. The water vapor temperature and pressure relation is assumed for the inner boundary. The water vapor temperature and pressure data were obtained from a thermodynamic steam table. Laplace transform and finite difference methods are used to analyze problems. The solution is obtained by using the matrix similarity transformation and inverse Laplace transform. We obtain solutions for the temperature and thermal deformation distributions in a transient and steady state. Moreover, the computational procedures established in this thesis, can solve the generalized magnetothermoelasticity problem for multilayered hollow cylinder with nonhomogeneous materials.     

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.