Thermal diffusivity identification based on an iterative regularization method

This article deals with the identification of a space and time dependent material thermal diffusivity. Such parameter is involved in heat transfers described by partial differential equations. An iterative regularization method based on a conjugate gradient algorithm is implemented. Such approach is attractive in order to efficiently deal with measurement noises and model errors. Numerical results are illustrated according to several simulations

Calibration and Rescaling Principles for Nonlinear Inverse Heat Conduction and Parameter Estimation Problems

This dissertation provides a systematic method for resolving nonlinear inverse heat conduction problems based on a calibration formulation and its accompanying principles. It is well-known that inverse heat conduction problems are ill-posed and hence subject to stability and uniqueness issues. Regularization methods are required to extract the best prediction based on a family of solutions. To date, most studies require sophisticated and combined numerical methods and regularization schemes for producing predictions. All thermophysical and geometrical properties must be provided in the simulations. The successful application of the numerical methods relies on the accuracy of the related system parameters as previously described. Due to the existence of uncertainties in the system parameters, these numerical methods possess bias of varying magnitudes. The calibration based approaches are proposed to minimize the systematic errors since system parameters are implicitly included in the mathematical formulation based on several calibration tests. To date, most calibration inverse studies have been based on the assumption of constant thermophysical properties. In contrast, this dissertation focuses on accounting for temperature-dependent thermophysical properties that produces a nonlinear heat equation. A novel rescaling principle is introduced for linearzing the system. This concept generates a mathematical framework similar to that of the linear formulation. Unlike the linear formulation, the present approach does require knowledge of thermophysical properties. However, all geometrical properties and sensor characterization are completely removed from the system. In this dissertation, a linear one-probe calibration method is first introduced as background. After that, the calibration method is generalized to the one-probe and two-probe, one-dimensional thermal system based on the assumption of temperature-dependent thermophysical properties. All previously proposed calibration equations are expressed in terms of a Volterra integral equation of the first kind for the unknown surface (net) heat flux and hence requires regularization owning to the ill-posed nature of first kind equations. A new strategy is proposed for determining the optimal regularization parameter that is independent of the applied regularization approach. As a final application, the described calibration principle is used for estimating unknown thermophysical properties above room temperature

Inverse problems and optimal experiment design in unsteady heat transfer processes identification

Experimental-computational methods for estimating characteristics of unsteady heat transfer processes are analyzed. The methods are based on the principles of distributed parameter system identification. The theoretical basis of such methods is the numerical solution of nonlinear ill-posed inverse heat transfer problems and optimal experiment design problems. Numerical techniques for solving problems are briefly reviewed. The results of the practical application of identification methods are demonstrated when estimating effective thermophysical characteristics of composite materials and thermal contact resistance in two-layer systems

An inverse problem of finding the time-dependent diffusion coefficient from an integral condition

We consider the inverse problem of determining the time-dependent diffusivity in one-dimensional heat equation with periodic boundary conditions and nonlocal over-specified data. The problem is highly nonlinear and it serves as a mathematical model for the technological process of external guttering applied in cleaning admixtures from silicon chips. First, the well-posedness conditions for the existence, uniqueness, and continuous dependence upon the data of the classical solution of the problem are established. Then, the problem is discretized using the finite-difference method and recasts as a nonlinear least-squares minimization problem with a simple positivity lower bound on the unknown diffusivity. Numerically, this is effectively solved using the lsqnonlin routine from the MATLAB toolbox. In order to investigate the accuracy, stability, and robustness of the numerical method, results for a few test examples are presented and discussed

Novel Numerical Approaches for the Resolution of Direct and Inverse Heat Transfer Problems

This dissertation describes an innovative and robust global time approach which has been developed for the resolution of direct and inverse problems, specifically in the disciplines of radiation and conduction heat transfer. Direct problems are generally well-posed and readily lend themselves to standard and well-defined mathematical solution techniques. Inverse problems differ in the fact that they tend to be ill-posed in the sense of Hadamard, i.e., small perturbations in the input data can produce large variations and instabilities in the output. The stability problem is exacerbated by the use of discrete experimental data which may be subject to substantial measurement error. This tendency towards ill-posedness is the main difficulty in developing a suitable prediction algorithm for most inverse problems. Previous attempts to overcome the inherent instability have involved the utilization of smoothing techniques such as Tikhonov regularization and sequential function estimation (Beck’s future information method). As alternatives to the existing methodologies, two novel mathematical schemes are proposed. They are the Global Time Method (GTM) and the Function Decomposition Method (FDM). Both schemes are capable of rendering time and space in a global fashion thus resolving the temporal and spatial domains simultaneously. This process effectively treats time elliptically or as a fourth spatial dimension. AWeighted Residuals Method (WRM) is utilized in the mathematical formulation wherein the unknown function is approximated in terms of a finite series expansion. Regularization of the solution is achieved by retention of expansion terms as opposed to smoothing in the classical Tikhonov sense. In order to demonstrate the merit and flexibility of these approaches, the GTM and FDM have been applied to representative problems of direct and inverse heat transfer. Those chosen are a direct problem of radiative transport, a parameter estimation problem found in Differential Scanning Calorimetry (DSC) and an inverse heat conduction problem (IHCP). The IHCP is resolved for the cases of diagnostic deduction (discrete temperature data at the boundary) and thermal design (prescribed functional data at the boundary). Both methods are shown to provide excellent results for the conditions under which they were tested. Finally, a number of suggestions for future work are offered

An automatic methodology for estimating eddy diffusivities from experimental data

A technique for estimating eddy diffusivities in a turbulent atmospheric layer is presented; the scheme adopted is based on an inverse-problem methodology. The inverse problem is formulated as a nonlinear constrained optimization problem, where the objective function is defined through the square differences between experimental and model data. The direct mathematical model is given by the advection-diffusion equation, which is solved by second-order finite-difference method. In the presence of noise it is necessary to use some regularization term; the Tikhonov function and an entropic regularization of zeroth, first and second orders are used in this paper. In addition, two inversion strategies are used: alternate and simultaneous eddy diffusivities estimation. Numerical experiments show a good performance of the proposed methodology