Combining Integral Transforms and Bayesian Inference in the Simultaneous Identification of Variable Thermal Conductivity and Thermal Capacity in Heterogeneous Media

2011 ◽  
Vol 133 (11) ◽  
Author(s):  
Carolina P. Naveira-Cotta ◽  
Helcio R. B. Orlande ◽  
Renato M. Cotta

This work presents the combined use of the integral transform method, for the direct problem solution, and of Bayesian inference, for the inverse problem analysis, in the simultaneous estimation of spatially variable thermal conductivity and thermal capacity for 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), offering a fast, precise, and robust solution for the transient temperature field. The inverse problem analysis employs a Markov chain Monte Carlo (MCMC) method, through the implementation of the Metropolis-Hastings sampling algorithm. Instead of seeking the functions estimation in the form of local values for the thermal conductivity and capacity, an alternative approach is employed based on the eigenfunction expansion of the thermophysical properties themselves. Then, the unknown parameters become the corresponding series coefficients for the properties eigenfunction expansions. Simulated temperatures obtained via integral transforms are used in the inverse analysis, for a prescribed concentration distribution of the dispersed phase in a heterogeneous media such as particle filled composites. Available correlations for the thermal conductivity and theory of mixtures relations for the thermal capacity are employed to produce the simulated results with high precision in the direct problem solution, while eigenfunction expansions with reduced number of terms are employed in the inverse analysis itself, in order to avoid the inverse crime. Gaussian distributions were used as priors for the parameter estimation procedure. In addition, simulated results with different randomly generated errors were employed in order to test the inverse analysis robustness.

Author(s):  
Carolina P. Naveira-Cotta ◽  
Helcio R. B. Orlande ◽  
Renato M. Cotta ◽  
Jeziel S. Nunes

This work deals with the simultaneous estimation of spatially variable thermal conductivity and diffusivity for one-dimensional heat conduction in 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 by Bayesian inference through a Markov Chain Monte Carlo (MCMC) method. Instead of seeking the function estimation in the form of a sequence of local values for the thermal properties, an alternative approach is utilized here, which is based on the eigenfunction expansion of the thermal conductivity and diffusivity themselves. Then, the unknown parameters become the corresponding expansion coefficients. In addition, the inverse analysis is performed on the transformed temperature field, instead of employing the actual local temperature measurements, thus promoting a significant data reduction through the integral transformation of the experimental measurements. A demonstration experiment is built involving a partially heated thin bakelite plate. Temperature measurements obtained via infrared thermography are used in the inverse analysis.


Author(s):  
D.K. Durdiev ◽  
J.Z. Nuriddinov

The inverse problem of determining a multidimensional kernel of an integral term depending on a time variable $t$ and $ (n-1)$-dimensional spatial variable $x'=\left(x_1,\ldots, x_ {n-1}\right)$ in the $n$-dimensional heat equation with a variable coefficient of thermal conductivity is investigated. The direct problem is the Cauchy problem for this equation. The integral term has the time convolution form of kernel and direct problem solution. As additional information for solving the inverse problem, the solution of the direct problem on the hyperplane $x_n = 0$ is given. At the beginning, the properties of the solution to the direct problem are studied. For this, the problem is reduced to solving an integral equation of the second kind of Volterra-type and the method of successive approximations is applied to it. Further the stated inverse problem is reduced to two auxiliary problems, in the second one of them an unknown kernel is included in an additional condition outside integral. Then the auxiliary problems are replaced by an equivalent closed system of Volterra-type integral equations with respect to unknown functions. Applying the method of contraction mappings to this system in the Hölder class of functions, we prove the main result of the article, which is a local existence and uniqueness theorem of the inverse problem solution.


Author(s):  
Kleber Marques Lisboa ◽  
Renato Machado Cotta

A hybrid numerical–analytical solution is developed for laminar flow development in a parallel plate duct partially filled with porous media. The integral transform method is employed in combination with a single domain reformulation strategy for representing the heterogeneous media within the channel. A novel eigenfunction expansion basis is proposed, including abrupt spatial variations of physical properties due to the domain transitions. The introduction of the new basis allows for a solution with similar convergence rates as in previous applications with simpler formulations, as demonstrated through a careful convergence analysis of the expansions. The inherent automatic error control characteristic of the integral transforms approach then provides benchmark results for the developing velocity profile. Moreover, a physical analysis further verifies the consistency of both the proposed expansion and the mixed symbolic–numerical code developed. A detailed verification with a finite-element commercial code is also performed.


Sign in / Sign up

Export Citation Format

Share Document