An EM algorithm for dynamic estimation of interfacial boundary in stratified flow of immiscible liquids using EIT

Abstract Estimation of interfacial boundary between two immiscible liquids in two-phase flows through pipe line provides information about the flow characteristics and thus can aid in design and monitoring of the flow process. The interfacial boundary can be represented in several ways, one such method is the front point approach. Front points describe the location and the shape of the interfacial boundary separating the immiscible liquids. During the flow process, due to fluctuations the interfacial boundary and so the front points which describe the boundary changes with time. The time-varying interfacial boundary can be estimated using dynamic inverse algorithms based on Kalman filter. However, algorithms based on Kalman filter require complete knowledge of model parameters (initial states, state transition matrix, and noise covariance matrices) for implementation. In processes involving complex flow pattern such as two-phase flows, it is difficult to represent the model parameters in a prior form. This uncertainty in model parameters causes suboptimal performance of the Kalman type filters. In this paper, we employ expectation maximization algorithm (EM) to estimate model parameters along with the interfacial boundary using electrical impedance tomography (EIT). The estimation of model parameters reduces the modeling uncertainty and thus results in improving the tracking of interfacial boundary. Numerical and experimental studies are performed to validate the performance of the proposed method.

[1]  Ville Kolehmainen,et al.  Free‐surface and admittivity estimation in electrical impedance tomography , 2006 .

[2]  David K. Han,et al.  Regular Article: A Shape Decomposition Technique in Electrical Impedance Tomography , 1999 .

[3]  Dan Nettleton,et al.  Convergence properties of the EM algorithm in constrained parameter spaces , 1999 .

[4]  Roger T. Bonnecaze,et al.  Inverse method for imaging a free surface using electrical impedance tomography , 2000 .

[5]  Alfred O. Hero,et al.  Time-delay estimation for filtered Poisson processes using an EM-type algorithm , 1994, IEEE Trans. Signal Process..

[6]  William J. Byrne,et al.  Alternating minimization and Boltzmann machine learning , 1992, IEEE Trans. Neural Networks.

[7]  Geoffrey E. Hinton,et al.  Parameter estimation for linear dynamical systems , 1996 .

[8]  Sin Kim,et al.  Estimation of void boundaries in flow field using expectation–maximization algorithm , 2011 .

[9]  D. Isaacson,et al.  Electrode models for electric current computed tomography , 1989, IEEE Transactions on Biomedical Engineering.

[10]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[11]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[12]  C. Jiang,et al.  The use of mixture models to detect effects of major genes on quantitative characters in a plant breeding experiment. , 1994, Genetics.

[13]  R. Mehra,et al.  Computational aspects of maximum likelihood estimation and reduction in sensitivity function calculations , 1974 .

[14]  Sin Kim,et al.  Nonstationary phase boundary estimation in electrical impedance tomography using unscented Kalman filter , 2008, J. Comput. Phys..

[15]  T. Savolainen,et al.  Estimating shapes and free surfaces with electrical impedance tomography , 2004 .

[16]  New York Dover,et al.  ON THE CONVERGENCE PROPERTIES OF THE EM ALGORITHM , 1983 .

[17]  Yuri V. Fairuzov,et al.  Numerical simulation of transient flow of two immiscible liquids in pipeline , 2000 .

[18]  Yaakov Bar-Shalom,et al.  Estimation and Tracking: Principles, Techniques, and Software , 1993 .

[19]  Umer Zeeshan Ijaz,et al.  Electrical resistance imaging of a time-varying interface in stratified flows using an unscented Kalman filter , 2008 .

[20]  A. Unnikrishnan,et al.  EM Based Extended Kalman Filter for Estimation of Rotor Time-Constant of Induction Motor , 2006, 2006 IEEE International Symposium on Industrial Electronics.

[21]  E. Somersalo,et al.  Existence and uniqueness for electrode models for electric current computed tomography , 1992 .

[22]  Sin Kim,et al.  Unscented Kalman filter approach to tracking a moving interfacial boundary in sedimentation processes using three-dimensional electrical impedance tomography , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[23]  Umer Zeeshan Ijaz,et al.  Phase boundary estimation in electrical impedance tomography using the Hooke and Jeeves pattern search method , 2010 .

[24]  Eric T. Chung,et al.  Electrical impedance tomography using level set representation and total variational regularization , 2005 .

[25]  Bong Seok Kim,et al.  An EKF based estimation scheme for sedimentation processes in vessels using EIT-type measurement data , 2010 .

[26]  Dong Liu,et al.  EM algorithm applied for estimating non-stationary region boundaries using electrical impedance tomography , 2010 .

[27]  Simon R. Arridge,et al.  RECOVERY OF REGION BOUNDARIES OF PIECEWISE CONSTANT COEFFICIENTS OF AN ELLIPTIC PDE FROM BOUNDARY DATA , 1999 .

[28]  T. Moon The expectation-maximization algorithm , 1996, IEEE Signal Process. Mag..

[29]  D. Rubin,et al.  On Jointly Estimating Parameters and Missing Data by Maximizing the Complete-Data Likelihood , 1983 .

[30]  C. Striebel,et al.  On the maximum likelihood estimates for linear dynamic systems , 1965 .

[31]  Umer Zeeshan Ijaz,et al.  Moving interfacial boundary estimation in stratified flow of two immiscible liquids using electrical resistance tomography , 2007 .

[32]  Bernhard Brandstätter,et al.  A FEM‐BEM approach using level‐sets in electrical capacitance tomography , 2005 .

[33]  R. Shumway,et al.  AN APPROACH TO TIME SERIES SMOOTHING AND FORECASTING USING THE EM ALGORITHM , 1982 .