Implicitly coupled phase fraction equations for the Eulerian multi-fluid model

Abstract In this work, the implementation, verification and validation of an implicitly coupled solution procedure for the phase fraction equations in the Eulerian multi-fluid model are presented. The model is implemented within the foam-extend toolbox, a community-driven fork of OpenFOAM. The implicitly coupled system for an arbitrary number of phases is based on the modified formulation of the phase fraction equation. This formulation takes advantage of the mixture divergence-free velocity and the cross-coupling with the remaining phase fraction equations via the non-linear relative velocity term. The linearised and implicitly coupled phase-fraction equations are solved simultaneously within a single block matrix. The model is tested for a bubbly air-water upward flow which forms a mixing layer inside a square duct. In the first test, the mesh verification analysis is performed on structured grids with different levels of refinement. The second test investigates the influence of the number of bubble phases on the flow solution for the same flow conditions. In the third test, the implemented model is validated against experimental data from the literature. The last test compares the performance of the implemented implicitly coupled solution procedure for the phase fraction equations against the standard segregated implementation. The proposed method shows good agreement with experimental data, and has proven to be consistent both in terms of the number of phases and grid refinement. Furthermore, the method improved the convergence of the solution for flows at higher bubble phase fraction.

[1]  S. Subramaniam.  Lagrangian-Eulerian methods for multiphase flows , 2013 .

[2]  A. D. Gosman,et al.  Multidimensional modeling of turbulent two‐phase flows in stirred vessels , 1992 .

[3]  V. Roig,et al.  Experimental investigation of a turbulent bubbly mixing layer , 1998 .

[4]  B. P. Leonard,et al.  Simple high-accuracy resolution program for convective modelling of discontinuities , 1988 .

[5]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[6]  B. Launder,et al.  The numerical computation of turbulent flows , 1990 .

[7]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[8]  Sergio Hoyas,et al.  Evaluation of the Eulerian-Lagrangian Spray Atomization (ELSA) model in spray simulations: 2D cases , 2013, Math. Comput. Model..

[9]  A. Gosman,et al.  High resolution NVD differencing scheme for arbitrarily unstructured meshes , 1999 .

[10]  Hrvoje Jasak,et al.  A coupled pressure based solution algorithm based on the volume-of-fluid approach for two or more immiscible fluids , 2010 .

[11]  Luiz Fernando L. R. Silva,et al.  Implementation of an implicit pressure–velocity coupling for the Eulerian multi-fluid model , 2019, Computers & Fluids.

[12]  Hrvoje Jasak,et al.  A tensorial approach to computational continuum mechanics using object-oriented techniques , 1998 .

[13]  H. Jasak,et al.  Block-selective algebraic multigrid for implicitly coupled pressure-velocity system , 2018 .

[14]  F. Durst,et al.  Eulerian and Lagrangian predictions of particulate two phase flows , 1984 .

[15]  A. Harten.  High Resolution Schemes for Hyperbolic Conservation Laws , 2017 .

[16]  Victor H. Ransom,et al.  Hyperbolic two-pressure models for two-phase flow , 1984 .

[17]  Luís Eça,et al.  A procedure for the estimation of the numerical uncertainty of CFD calculations based on grid refinement studies , 2014, J. Comput. Phys..

[18]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

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

[20]  Luiz Fernando L. R. Silva,et al.  Development and implementation of a polydispersed multiphase flow model in OpenFOAM , 2011, Comput. Chem. Eng..

[21]  J-M Buchlin,et al.  Numerical computation of multi-phase flows , 1981 .

[22]  H. Ström,et al.  On the dynamics of instabilities in two-fluid models for bubbly flows , 2017 .

[23]  M. Sommerfeld,et al.  Multiphase Flows with Droplets and Particles , 2011 .

[24]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[25]  M. B. Carver.  A Method of Limiting Intermediate Values of Volume Fraction in Iterative Two-Fluid Computations , 1982 .

[26]  R. I. Issa,et al.  Numerical aspects of an algorithm for the Eulerian simulation of two‐phase flows , 2003 .

[27]  J. M. Nóbrega,et al.  A coupled finite volume flow solver for the solution of incompressible viscoelastic flows , 2019, Journal of Non-Newtonian Fluid Mechanics.

[28]  C. Dopazo,et al.  On conditioned averages for intermittent turbulent flows , 1977, Journal of Fluid Mechanics.

[29]  Alberto Passalacqua,et al.  On the hyperbolicity of the two-fluid model for gas–liquid bubbly flows , 2018 .

[30]  M. Ishii,et al.  Thermo-Fluid Dynamics of Two-Phase Flow , 2007 .