16,542 research outputs found
Unified Gas-kinetic Wave-Particle Methods III: Multiscale Photon Transport
In this paper, we extend the unified gas-kinetic wave-particle (UGKWP) method
to the multiscale photon transport. In this method, the photon free streaming
and scattering processes are treated in an un-splitting way. The duality
descriptions, namely the simulation particle and distribution function, are
utilized to describe the photon. By accurately recovering the governing
equations of the unified gas-kinetic scheme (UGKS), the UGKWP preserves the
multiscale dynamics of photon transport from optically thin to optically thick
regime. In the optically thin regime, the UGKWP becomes a Monte Carlo type
particle tracking method, while in the optically thick regime, the UGKWP
becomes a diffusion equation solver. The local photon dynamics of the UGKWP, as
well as the proportion of wave-described and particle-described photons are
automatically adapted according to the numerical resolution and transport
regime. Compared to the -type UGKS, the UGKWP requires less memory cost
and does not suffer ray effect. Compared to the implicit Monte Carlo (IMC)
method, the statistical noise of UGKWP is greatly reduced and computational
efficiency is significantly improved in the optically thick regime. Several
numerical examples covering all transport regimes from the optically thin to
optically thick are computed to validate the accuracy and efficiency of the
UGKWP method. In comparison to the -type UGKS and IMC method, the UGKWP
method may have several-order-of-magnitude reduction in computational cost and
memory requirement in solving some multsicale transport problems.Comment: 27 pages, 15 figures. arXiv admin note: text overlap with
arXiv:1810.0598
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
Exponential Runge-Kutta methods for stiff kinetic equations
We introduce a class of exponential Runge-Kutta integration methods for
kinetic equations. The methods are based on a decomposition of the collision
operator into an equilibrium and a non equilibrium part and are exact for
relaxation operators of BGK type. For Boltzmann type kinetic equations they
work uniformly for a wide range of relaxation times and avoid the solution of
nonlinear systems of equations even in stiff regimes. We give sufficient
conditions in order that such methods are unconditionally asymptotically stable
and asymptotic preserving. Such stability properties are essential to guarantee
the correct asymptotic behavior for small relaxation times. The methods also
offer favorable properties such as nonnegativity of the solution and entropy
inequality. For this reason, as we will show, the methods are suitable both for
deterministic as well as probabilistic numerical techniques
ORB5: a global electromagnetic gyrokinetic code using the PIC approach in toroidal geometry
This paper presents the current state of the global gyrokinetic code ORB5 as
an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177
409 (2007)]. The ORB5 code solves the electromagnetic Vlasov-Maxwell system of
equations using a PIC scheme and also includes collisions and strong flows. The
code assumes multiple gyrokinetic ion species at all wavelengths for the
polarization density and drift-kinetic electrons. Variants of the physical
model can be selected for electrons such as assuming an adiabatic response or a
``hybrid'' model in which passing electrons are assumed adiabatic and trapped
electrons are drift-kinetic. A Fourier filter as well as various control
variates and noise reduction techniques enable simulations with good
signal-to-noise ratios at a limited numerical cost. They are completed with
different momentum and zonal flow-conserving heat sources allowing for
temperature-gradient and flux-driven simulations. The code, which runs on both
CPUs and GPUs, is well benchmarked against other similar codes and analytical
predictions, and shows good scalability up to thousands of nodes
Kinetic Solvers with Adaptive Mesh in Phase Space
An Adaptive Mesh in Phase Space (AMPS) methodology has been developed for
solving multi-dimensional kinetic equations by the discrete velocity method. A
Cartesian mesh for both configuration (r) and velocity (v) spaces is produced
using a tree of trees data structure. The mesh in r-space is automatically
generated around embedded boundaries and dynamically adapted to local solution
properties. The mesh in v-space is created on-the-fly for each cell in r-space.
Mappings between neighboring v-space trees implemented for the advection
operator in configuration space. We have developed new algorithms for solving
the full Boltzmann and linear Boltzmann equations with AMPS. Several recent
innovations were used to calculate the discrete Boltzmann collision integral
with dynamically adaptive mesh in velocity space: importance sampling,
multi-point projection method, and the variance reduction method. We have
developed an efficient algorithm for calculating the linear Boltzmann collision
integral for elastic and inelastic collisions in a Lorentz gas. New AMPS
technique has been demonstrated for simulations of hypersonic rarefied gas
flows, ion and electron kinetics in weakly ionized plasma, radiation and light
particle transport through thin films, and electron streaming in
semiconductors. We have shown that AMPS allows minimizing the number of cells
in phase space to reduce computational cost and memory usage for solving
challenging kinetic problems
- …