1,168 research outputs found
A Multiscale Kinetic-Fluid Solver with Dynamic Localization of Kinetic Effects
This paper collects the efforts done in our previous works [P. Degond, S.
Jin, L. Mieussens, A Smooth Transition Between Kinetic and Hydrodynamic
Equations, J. Comp. Phys., 209 (2005) 665--694.],[P.Degond, G. Dimarco, L.
Mieussens, A Moving Interface Method for Dynamic Kinetic-fluid Coupling, J.
Comp. Phys., Vol. 227, pp. 1176-1208, (2007).],[P. Degond, J.G. Liu, L.
Mieussens, Macroscopic Fluid Model with Localized Kinetic Upscaling Effects,
SIAM Multi. Model. Sim. 5(3), 940--979 (2006)] to build a robust multiscale
kinetic-fluid solver. Our scope is to efficiently solve fluid dynamic problems
which present non equilibrium localized regions that can move, merge, appear or
disappear in time. The main ingredients of the present work are the followings
ones: a fluid model is solved in the whole domain together with a localized
kinetic upscaling term that corrects the fluid model wherever it is necessary;
this multiscale description of the flow is obtained by using a micro-macro
decomposition of the distribution function [P. Degond, J.G. Liu, L. Mieussens,
Macroscopic Fluid Model with Localized Kinetic Upscaling Effects, SIAM Multi.
Model. Sim. 5(3), 940--979 (2006)]; the dynamic transition between fluid and
kinetic descriptions is obtained by using a time and space dependent transition
function; to efficiently define the breakdown conditions of fluid models we
propose a new criterion based on the distribution function itself. Several
numerical examples are presented to validate the method and measure its
computational efficiency.Comment: 34 page
Multiphysics simulations of collisionless plasmas
Collisionless plasmas, mostly present in astrophysical and space
environments, often require a kinetic treatment as given by the Vlasov
equation. Unfortunately, the six-dimensional Vlasov equation can only be solved
on very small parts of the considered spatial domain. However, in some cases,
e.g. magnetic reconnection, it is sufficient to solve the Vlasov equation in a
localized domain and solve the remaining domain by appropriate fluid models. In
this paper, we describe a hierarchical treatment of collisionless plasmas in
the following way. On the finest level of description, the Vlasov equation is
solved both for ions and electrons. The next courser description treats
electrons with a 10-moment fluid model incorporating a simplified treatment of
Landau damping. At the boundary between the electron kinetic and fluid region,
the central question is how the fluid moments influence the electron
distribution function. On the next coarser level of description the ions are
treated by an 10-moment fluid model as well. It may turn out that in some
spatial regions far away from the reconnection zone the temperature tensor in
the 10-moment description is nearly isotopic. In this case it is even possible
to switch to a 5-moment description. This change can be done separately for
ions and electrons. To test this multiphysics approach, we apply this full
physics-adaptive simulations to the Geospace Environmental Modeling (GEM)
challenge of magnetic reconnection.Comment: 13 pages, 5 figure
The Moment Guided Monte Carlo method for the Boltzmann equation
In this work we propose a generalization of the Moment Guided Monte Carlo
method developed in [11]. This approach permits to reduce the variance of the
particle methods through a matching with a set of suitable macroscopic moment
equations. In order to guarantee that the moment equations provide the correct
solutions, they are coupled to the kinetic equation through a non equilibrium
term. Here, at the contrary to the previous work in which we considered the
simplified BGK operator, we deal with the full Boltzmann operator. Moreover, we
introduce an hybrid setting which permits to entirely remove the resolution of
the kinetic equation in the limit of infinite number of collisions and to
consider only the solution of the compressible Euler equation. This
modification additionally reduce the statistical error with respect to our
previous work and permits to perform simulations of non equilibrium gases using
only a few number of particles. We show at the end of the paper several
numerical tests which prove the efficiency and the low level of numerical noise
of the method.Comment: arXiv admin note: text overlap with arXiv:0908.026
Asymptotic-Preserving Monte Carlo methods for transport equations in the diffusive limit
We develop a new Monte Carlo method that solves hyperbolic transport
equations with stiff terms, characterized by a (small) scaling parameter. In
particular, we focus on systems which lead to a reduced problem of parabolic
type in the limit when the scaling parameter tends to zero. Classical Monte
Carlo methods suffer of severe time step limitations in these situations, due
to the fact that the characteristic speeds go to infinity in the diffusion
limit. This makes the problem a real challenge, since the scaling parameter may
differ by several orders of magnitude in the domain. To circumvent these time
step limitations, we construct a new, asymptotic-preserving Monte Carlo method
that is stable independently of the scaling parameter and degenerates to a
standard probabilistic approach for solving the limiting equation in the
diffusion limit. The method uses an implicit time discretization to formulate a
modified equation in which the characteristic speeds do not grow indefinitely
when the scaling factor tends to zero. The resulting modified equation can
readily be discretized by a Monte Carlo scheme, in which the particles combine
a finite propagation speed with a time-step dependent diffusion term. We show
the performance of the method by comparing it with standard (deterministic)
approaches in the literature
Fluid Simulations with Localized Boltzmann Upscaling by Direct Simulation Monte-Carlo
In the present work, we present a novel numerical algorithm to couple the
Direct Simulation Monte Carlo method (DSMC) for the solution of the Boltzmann
equation with a finite volume like method for the solution of the Euler
equations. Recently we presented in [14],[16],[17] different methodologies
which permit to solve fluid dynamics problems with localized regions of
departure from thermodynamical equilibrium. The methods rely on the
introduction of buffer zones which realize a smooth transition between the
kinetic and the fluid regions. In this paper we extend the idea of buffer zones
and dynamic coupling to the case of the Monte Carlo methods. To facilitate the
coupling and avoid the onset of spurious oscillations in the fluid regions
which are consequences of the coupling with a stochastic numerical scheme, we
use a new technique which permits to reduce the variance of the particle
methods [11]. In addition, the use of this method permits to obtain estimations
of the breakdowns of the fluid models less affected by fluctuations and
consequently to reduce the kinetic regions and optimize the coupling. In the
last part of the paper several numerical examples are presented to validate the
method and measure its computational performances
A Hierarchy of Hybrid Numerical Methods for Multi-Scale Kinetic Equations
In this paper, we construct a hierarchy of hybrid numerical methods for
multi-scale kinetic equations based on moment realizability matrices, a concept
introduced by Levermore, Morokoff and Nadiga. Following such a criterion, one
can consider hybrid scheme where the hydrodynamic part is given either by the
compressible Euler or Navier-Stokes equations, or even with more general
models, such as the Burnett or super-Burnett systems.Comment: 27 pages, edit: typo and metadata chang
- …