TY - JOUR

T1 - Stabilization of the Eulerian model for incompressible multiphase flow by artificial diffusion

AU - Vreman, A.W.

PY - 2011/2/20

Y1 - 2011/2/20

N2 - The commonly used Eulerian or continuum model for incompressible multiphase flow is known to be unstable to perturbations for all wavenumbers, even if viscosity terms are used in the momentum equations. In the present work the model is stabilized by adding explicit artificial diffusion to the mass equations. The artificial diffusion terms lead to improved stability properties: uniform flow becomes linearly stable for large wavenumbers, and above an analytically derived threshold for the artificial diffusivity, stability for all wavenumbers is achieved. The artificial diffusivity reappears in the momentum equations, in such a way that fundamental properties of the standard equations remain valid: Galilean invariance is maintained, total mass and momentum are conserved, decay of total kinetic energy is ensured in the absence of external forces, and a flow initially at rest at hydrostatic pressure remains unchanged, even if the spatial distribution of volume fractions is nonuniform. A staggered finite volume pressure correction method using central differencing (leading to energy conserving discretization of convective and pressure terms) is presented. Application of the method to one-dimensional two-phase flow of falling particles particles confirms that the equations are stable with and unstable without artificial diffusion in the volume fraction equation.

AB - The commonly used Eulerian or continuum model for incompressible multiphase flow is known to be unstable to perturbations for all wavenumbers, even if viscosity terms are used in the momentum equations. In the present work the model is stabilized by adding explicit artificial diffusion to the mass equations. The artificial diffusion terms lead to improved stability properties: uniform flow becomes linearly stable for large wavenumbers, and above an analytically derived threshold for the artificial diffusivity, stability for all wavenumbers is achieved. The artificial diffusivity reappears in the momentum equations, in such a way that fundamental properties of the standard equations remain valid: Galilean invariance is maintained, total mass and momentum are conserved, decay of total kinetic energy is ensured in the absence of external forces, and a flow initially at rest at hydrostatic pressure remains unchanged, even if the spatial distribution of volume fractions is nonuniform. A staggered finite volume pressure correction method using central differencing (leading to energy conserving discretization of convective and pressure terms) is presented. Application of the method to one-dimensional two-phase flow of falling particles particles confirms that the equations are stable with and unstable without artificial diffusion in the volume fraction equation.

KW - Artificial diffusion

KW - Multicomponent flow

KW - Multiphase flow

KW - Stability

KW - Staggered pressure correction methods

UR - http://www.scopus.com/inward/record.url?scp=78650559434&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2010.11.025

DO - 10.1016/j.jcp.2010.11.025

M3 - Article

AN - SCOPUS:78650559434

VL - 230

SP - 1639

EP - 1651

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 4

ER -