3 research outputs found
Phase appearance or disappearance in two-phase flows
This paper is devoted to the treatment of specific numerical problems which
appear when phase appearance or disappearance occurs in models of two-phase
flows. Such models have crucial importance in many industrial areas such as
nuclear power plant safety studies. In this paper, two outstanding problems are
identified: first, the loss of hyperbolicity of the system when a phase appears
or disappears and second, the lack of positivity of standard shock capturing
schemes such as the Roe scheme. After an asymptotic study of the model, this
paper proposes accurate and robust numerical methods adapted to the simulation
of phase appearance or disappearance. Polynomial solvers are developed to avoid
the use of eigenvectors which are needed in usual shock capturing schemes, and
a method based on an adaptive numerical diffusion is designed to treat the
positivity problems. An alternate method, based on the use of the hyperbolic
tangent function instead of a polynomial, is also considered. Numerical results
are presented which demonstrate the efficiency of the proposed solutions
An Asymptotic-Preserving all-speed scheme for the Euler and Navier-Stokes equations
We present an Asymptotic-Preserving 'all-speed' scheme for the simulation of
compressible flows valid at all Mach-numbers ranging from very small to order
unity. The scheme is based on a semi-implicit discretization which treats the
acoustic part implicitly and the convective and diffusive parts explicitly.
This discretization, which is the key to the Asymptotic-Preserving property,
provides a consistent approximation of both the hyperbolic compressible regime
and the elliptic incompressible regime. The divergence-free condition on the
velocity in the incompressible regime is respected, and an the pressure is
computed via an elliptic equation resulting from a suitable combination of the
momentum and energy equations. The implicit treatment of the acoustic part
allows the time-step to be independent of the Mach number. The scheme is
conservative and applies to steady or unsteady flows and to general equations
of state. One and Two-dimensional numerical results provide a validation of the
Asymptotic-Preserving 'all-speed' properties