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

Abstract 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.

[1]  Robert D. Skeel,et al.  An explicit Runge-Kutta-Nystro¨m method is canonical if and only if its adjoint is explicit , 1992 .

[2]  Inverse estimation of the boundary condition in three-dimensional heat conduction , 1997 .

[3]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[4]  Gad Hetsroni,et al.  Heat transfer to a liquid—solid mixture in a flume , 1994 .

[5]  M. Krízek,et al.  On the maximum angle condition for linear tetrahedral elements , 1992 .

[6]  O. Alifanov.  Inverse heat transfer problems , 1994 .

[7]  H. Engl,et al.  Regularization of Inverse Problems , 1996 .

[8]  B. Blackwell,et al.  Comparison of some inverse heat conduction methods using experimental data , 1996 .

[9]  P. Hansen.  Rank-Deficient and Discrete Ill-Posed Problems: Numerical Aspects of Linear Inversion , 1987 .

[10]  I. Mudawar,et al.  Statistical investigation of the relationship between interfacial waviness and sensible heat transfer to a falling liquid film , 1991 .

[11]  W. Marquardt.  MODEL-BASED EXPERIMENTAL ANALYSIS : A SYSTEMS APPROACH TO MECHANISTIC MODELING OF KINETIC PHENOMENA , 2004 .

[12]  A. Reusken,et al.  The DROPS package for numerical simulations of incompressible flows using parallel adaptive multigrid techniques , 2022 .

[13]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[14]  I. Babuska,et al.  ON THE ANGLE CONDITION IN THE FINITE ELEMENT METHOD , 1976 .

[15]  Wolfgang Marquardt,et al.  Model-Based Experimental Analysis of Kinetic Phenomena in Multi-Phase Reactive Systems , 2005 .

[16]  Cheng-Hung Huang,et al.  A three-dimensional inverse heat conduction problem in estimating surface heat flux by conjugate gradient method , 1999 .

[17]  Wolfgang Marquardt,et al.  DESIGN, FORMULATION, AND SOLUTION OF MULTIDIMENSIONAL INVERSE HEAT CONDUCTION PROBLEMS , 2005 .

[18]  U. Renz,et al.  Local and instantaneous distribution of heat transfer rates through wavy films , 2002 .

[19]  Joël Bertrand,et al.  EditorialThe Changing Face of Chemical Engineering Research , 2005 .

[20]  H. Auracher,et al.  Heat and mass transfer in wavy falling films of binary mixtures , 1999 .

[21]  E. Haber,et al.  A multigrid method for distributed parameter estimation problems. , 2003 .

[22]  Somchart Chantasiriwan,et al.  An algorithm for solving multidimensional inverse heat conduction problem , 2001 .

[23]  Albert Mosyak,et al.  Heat transfer to two-phase flow in inclined tubes , 2003 .

[24]  Ulrich Renz,et al.  Measurement of instantaneous heat transfer using a hot-foil infrared technique , 2001, SPIE Defense + Commercial Sensing.