Mathematical Modelling and Numerical Analysis an Entropy Preserving Finite-element/finite-volume Pressure Correction Scheme for the Drift-flux Model

We present in this paper a pressure correction scheme for the drift-flux model combining finite element and finite volume discretizations, which is shown to enjoy essential stability features of the continuous problem: the scheme is conservative, the unknowns are kept within their physical bounds and, in the homogeneous case (i.e. when the drift velocity vanishes), the discrete entropy of the system decreases; in addition, when using for the drift velocity a closure law which takes the form of a Darcy-like relation, the drift term becomes dissipative. Finally, the present algorithm preserves a constant pressure and a constant velocity through moving interfaces between phases. To ensure the stability as well as to obtain this latter property, a key ingredient is to couple the mass balance and the transport equation for the dispersed phase in an original pressure correction step. The existence of a solution to each step of the algorithm is proven; in particular, the existence of a solution to the pressure correction step is derived as a consequence of a more general existence result for discrete problems associated to the drift-flux model. Numerical tests show a near-first-order convergence rate for the scheme, both in time and space, and confirm its stability.

[1]  F. Browder.  Nonlinear functional analysis , 1970 .

[2]  P. Raviart,et al.  Conforming and nonconforming finite element methods for solving the stationary Stokes equations I , 1973 .

[3]  Gerhart Eigenberger,et al.  Gas—liquid flow in bubble columns and loop reactors: Part II. Comparison of detailed experiments and flow simulations , 1994 .

[4]  L. Quartapelle,et al.  A projection FEM for variable density incompressible flows , 2000 .

[5]  Miloslav Feistauer,et al.  Mathematical and Computational Methods for Compressible Flow , 2003 .

[6]  Gerhart Eigenberger,et al.  Applicability of the standard k–ε turbulence model to the dynamic simulation of bubble columns: Part I. Detailed numerical simulations , 1999 .

[7]  Jean-Marc Hérard,et al.  On the use of symmetrizing variables for vacuums , 2003 .

[8]  B. Larrouturou.  How to preserve the mass fractions positivity when computing compressible multi-component flows , 1991 .

[9]  G. Eigenberger,et al.  Simulation of Buoyancy Driven Bubbly Flow: Established Simplifications and Open Questions (Journal Review) , 2004 .

[10]  M. Darwish,et al.  A pressure-based algorithm for multi-phase flow at all speeds , 2003 .

[11]  Jean E. Roberts,et al.  Mixed and hybrid finite element methods , 1987 .

[12]  Jie Shen,et al.  An overview of projection methods for incompressible flows , 2006 .

[13]  A. A. Amsden,et al.  A numerical fluid dynamics calculation method for all flow speeds , 1971 .

[14]  R. Eymard,et al.  Error estimates for the approximate solutions of a nonlinear hyperbolic equation given by finite volume schemes , 1998 .

[15]  R. Rannacher,et al.  Simple nonconforming quadrilateral Stokes element , 1992 .

[16]  R. Temam,et al.  Navier-stokes equations: Theory and approximation , 1998 .

[17]  P. Wesseling.  Principles of Computational Fluid Dynamics , 2000 .

[18]  Hervé Guillard,et al.  A Darcy law for the drift velocity in a two-phase flow model , 2007, J. Comput. Phys..

[19]  Raphaele Herbin,et al.  An unconditionnally stable pressure correction scheme for compressible barotropic Navier-Stokes equations , 2007, 0710.2987.

[20]  Miloslav Feistauer,et al.  Error Estimates for Barycentric Finite Volumes Combined with Nonconforming Finite Elements Applied to Nonlinear Convection-Diffusion Problems , 2002 .

[21]  D. Kuzmin,et al.  NUMERICAL SIMULATION OF TURBULENT BUBBLY FLOWS , 2004 .

[22]  Thierry Gallouët,et al.  A convergent finite element-finite volume scheme for the compressible Stokes problem. Part I: The isothermal case , 2007, Math. Comput..

[23]  P. G. Ciarlet,et al.  Basic error estimates for elliptic problems , 1991 .

[24]  Philippe Angot,et al.  Analysis of a combined barycentric finite volume—nonconforming finite element method for nonlinear convection-diffusion problems , 1998 .