2,258 research outputs found
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
An Asymptotic Preserving Scheme for the ES-BGK model
In this paper, we study a time discrete scheme for the initial value problem
of the ES-BGK kinetic equation. Numerically solving these equations are
challenging due to the nonlinear stiff collision (source) terms induced by
small mean free or relaxation time. We study an implicit-explicit (IMEX) time
discretization in which the convection is explicit while the relaxation term is
implicit to overcome the stiffness. We first show how the implicit relaxation
can be solved explicitly, and then prove asymptotically that this time
discretization drives the density distribution toward the local Maxwellian when
the mean free time goes to zero while the numerical time step is held fixed.
This naturally imposes an asymptotic-preserving scheme in the Euler limit. The
scheme so designed does not need any nonlinear iterative solver for the
implicit relaxation term. Moreover, it can capture the macroscopic fluid
dynamic (Euler) limit even if the small scale determined by the Knudsen number
is not numerically resolved. We also show that it is consistent to the
compressible Navier-Stokes equations if the viscosity and heat conductivity are
numerically resolved. Several numerical examples, in both one and two space
dimensions, are used to demonstrate the desired behavior of this scheme
Flux Splitting for stiff equations: A notion on stability
For low Mach number flows, there is a strong recent interest in the
development and analysis of IMEX (implicit/explicit) schemes, which rely on a
splitting of the convective flux into stiff and nonstiff parts. A key
ingredient of the analysis is the so-called Asymptotic Preserving (AP)
property, which guarantees uniform consistency and stability as the Mach number
goes to zero. While many authors have focussed on asymptotic consistency, we
study asymptotic stability in this paper: does an IMEX scheme allow for a CFL
number which is independent of the Mach number? We derive a stability criterion
for a general linear hyperbolic system. In the decisive eigenvalue analysis,
the advective term, the upwind diffusion and a quadratic term stemming from the
truncation in time all interact in a subtle way. As an application, we show
that a new class of splittings based on characteristic decomposition, for which
the commutator vanishes, avoids the deterioration of the time step which has
sometimes been observed in the literature
On the Eulerian Large Eddy Simulation of disperse phase flows: an asymptotic preserving scheme for small Stokes number flows
In the present work, the Eulerian Large Eddy Simulation of dilute disperse
phase flows is investigated. By highlighting the main advantages and drawbacks
of the available approaches in the literature, a choice is made in terms of
modelling: a Fokker-Planck-like filtered kinetic equation proposed by Zaichik
et al. 2009 and a Kinetic-Based Moment Method (KBMM) based on a Gaussian
closure for the NDF proposed by Vie et al. 2014. The resulting Euler-like
system of equations is able to reproduce the dynamics of particles for small to
moderate Stokes number flows, given a LES model for the gaseous phase, and is
representative of the generic difficulties of such models. Indeed, it
encounters strong constraints in terms of numerics in the small Stokes number
limit, which can lead to a degeneracy of the accuracy of standard numerical
methods. These constraints are: 1/as the resulting sound speed is inversely
proportional to the Stokes number, it is highly CFL-constraining, and 2/the
system tends to an advection-diffusion limit equation on the number density
that has to be properly approximated by the designed scheme used for the whole
range of Stokes numbers. Then, the present work proposes a numerical scheme
that is able to handle both. Relying on the ideas introduced in a different
context by Chalons et al. 2013: a Lagrange-Projection, a relaxation formulation
and a HLLC scheme with source terms, we extend the approach to a singular flux
as well as properly handle the energy equation. The final scheme is proven to
be Asymptotic-Preserving on 1D cases comparing to either converged or
analytical solutions and can easily be extended to multidimensional
configurations, thus setting the path for realistic applications
A class of asymptotic preserving schemes for kinetic equations and related problems with stiff sources
In this paper, we propose a general framework to design asymptotic preserving
schemes for the Boltzmann kinetic kinetic and related equations. Numerically
solving these equations are challenging due to the nonlinear stiff collision
(source) terms induced by small mean free or relaxation time. We propose to
penalize the nonlinear collision term by a BGK-type relaxation term, which can
be solved explicitly even if discretized implicitly in time. Moreover, the
BGK-type relaxation operator helps to drive the density distribution toward the
local Maxwellian, thus natually imposes an asymptotic-preserving scheme in the
Euler limit. The scheme so designed does not need any nonlinear iterative
solver or the use of Wild Sum. It is uniformly stable in terms of the (possibly
small) Knudsen number, and can capture the macroscopic fluid dynamic (Euler)
limit even if the small scale determined by the Knudsen number is not
numerically resolved. It is also consistent to the compressible Navier-Stokes
equations if the viscosity and heat conductivity are numerically resolved. The
method is applicable to many other related problems, such as hyperbolic systems
with stiff relaxation, and high order parabilic equations
- …