2,302 research outputs found
Asymptotic-preserving exponential methods for the quantum Boltzmann equation with high-order accuracy
In this paper we develop high-order asymptotic-preserving methods for the
spatially inhomogeneous quantum Boltzmann equation. We follow the work in Li
and Pareschi, where asymptotic preserving exponential Runge-Kutta methods for
the classical inhomogeneous Boltzmann equation were constructed. A major
difficulty here is related to the non Gaussian steady states characterizing the
quantum kinetic behavior. We show that the proposed schemes work with
high-order accuracy uniformly in time for all Planck constants ranging from
classical regime to quantum regime, and all Knudsen numbers ranging from
kinetic regime to fluid regime. Computational results are presented for both
Bose gas and Fermi gas
Numerical study of Bose-Einstein condensation in the Kaniadakis-Quarati model for bosons
Kaniadakis and Quarati (1994) proposed a Fokker--Planck equation with
quadratic drift as a PDE model for the dynamics of bosons in the spatially
homogeneous setting. It is an open question whether this equation has solutions
exhibiting condensates in finite time. The main analytical challenge lies in
the continuation of exploding solutions beyond their first blow-up time while
having a linear diffusion term. We present a thoroughly validated time-implicit
numerical scheme capable of simulating solutions for arbitrarily large time,
and thus enabling a numerical study of the condensation process in the
Kaniadakis--Quarati model. We show strong numerical evidence that above the
critical mass rotationally symmetric solutions of the Kaniadakis--Quarati model
in 3D form a condensate in finite time and converge in entropy to the unique
minimiser of the natural entropy functional at an exponential rate. Our
simulations further indicate that the spatial blow-up profile near the origin
follows a universal power law and that transient condensates can occur for
sufficiently concentrated initial data.Comment: To appear in Kinet. Relat. Model
Matrix-valued Quantum Lattice Boltzmann Method
We devise a lattice Boltzmann method (LBM) for a matrix-valued quantum
Boltzmann equation, with the classical Maxwell distribution replaced by
Fermi-Dirac functions. To accommodate the spin density matrix, the distribution
functions become 2 x 2 matrix-valued. From an analytic perspective, the
efficient, commonly used BGK approximation of the collision operator is valid
in the present setting. The numerical scheme could leverage the principles of
LBM for simulating complex spin systems, with applications to spintronics.Comment: 18 page
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows IV: full Boltzmann and Model Equations
Fluid dynamic equations are valid in their respective modeling scales. With a
variation of the modeling scales, theoretically there should have a continuous
spectrum of fluid dynamic equations. In order to study multiscale flow
evolution efficiently, the dynamics in the computational fluid has to be
changed with the scales. A direct modeling of flow physics with a changeable
scale may become an appropriate approach. The unified gas-kinetic scheme (UGKS)
is a direct modeling method in the mesh size scale, and its underlying flow
physics depends on the resolution of the cell size relative to the particle
mean free path. The cell size of UGKS is not limited by the particle mean free
path. With the variation of the ratio between the numerical cell size and local
particle mean free path, the UGKS recovers the flow dynamics from the particle
transport and collision in the kinetic scale to the wave propagation in the
hydrodynamic scale.
The previous UGKS is mostly constructed from the evolution solution of
kinetic model equations. This work is about the further development of the UGKS
with the implementation of the full Boltzmann collision term in the region
where it is needed. The central ingredient of the UGKS is the coupled treatment
of particle transport and collision in the flux evaluation across a cell
interface, where a continuous flow dynamics from kinetic to hydrodynamic scales
is modeled. The newly developed UGKS has the asymptotic preserving (AP)
property of recovering the NS solutions in the continuum flow regime, and the
full Boltzmann solution in the rarefied regime. In the mostly unexplored
transition regime, the UGKS itself provides a valuable tool for the flow study
in this regime. The mathematical properties of the scheme, such as stability,
accuracy, and the asymptotic preserving, will be analyzed in this paper as
well
Quantum Lattice Boltzmann is a quantum walk
Numerical methods for the 1-D Dirac equation based on operator splitting and
on the quantum lattice Boltzmann (QLB) schemes are reviewed. It is shown that
these discretizations fall within the class of quantum walks, i.e. discrete
maps for complex fields, whose continuum limit delivers Dirac-like relativistic
quantum wave equations. The correspondence between the quantum walk dynamics
and these numerical schemes is given explicitly, allowing a connection between
quantum computations, numerical analysis and lattice Boltzmann methods. The QLB
method is then extended to the Dirac equation in curved spaces and it is
demonstrated that the quantum walk structure is preserved. Finally, it is
argued that the existence of this link between the discretized Dirac equation
and quantum walks may be employed to simulate relativistic quantum dynamics on
quantum computers.Comment: 18 pages, 3 figure
An Asymptotic Preserving Maxwell Solver Resulting in the Darwin Limit of Electrodynamics
In plasma simulations, where the speed of light divided by a characteristic
length is at a much higher frequency than other relevant parameters in the
underlying system, such as the plasma frequency, implicit methods begin to play
an important role in generating efficient solutions in these multi-scale
problems. Under conditions of scale separation, one can rescale Maxwell's
equations in such a way as to give a magneto static limit known as the Darwin
approximation of electromagnetics. In this work, we present a new approach to
solve Maxwell's equations based on a Method of Lines Transpose (MOL)
formulation, combined with a fast summation method with computational
complexity , where is the number of grid points (particles).
Under appropriate scaling, we show that the proposed schemes result in
asymptotic preserving methods that can recover the Darwin limit of
electrodynamics
Entropies from coarse-graining: convex polytopes vs. ellipsoids
We examine the Boltzmann/Gibbs/Shannon and the
non-additive Havrda-Charv\'{a}t / Dar\'{o}czy/Cressie-Read/Tsallis \
\ and the Kaniadakis -entropy \ \
from the viewpoint of coarse-graining, symplectic capacities and convexity. We
argue that the functional form of such entropies can be ascribed to a
discordance in phase-space coarse-graining between two generally different
approaches: the Euclidean/Riemannian metric one that reflects independence and
picks cubes as the fundamental cells and the symplectic/canonical one that
picks spheres/ellipsoids for this role. Our discussion is motivated by and
confined to the behaviour of Hamiltonian systems of many degrees of freedom. We
see that Dvoretzky's theorem provides asymptotic estimates for the minimal
dimension beyond which these two approaches are close to each other. We state
and speculate about the role that dualities may play in this viewpoint.Comment: 63 pages. No figures. Standard LaTe
- âŠ