2,302 research outputs found

    Asymptotic-preserving exponential methods for the quantum Boltzmann equation with high-order accuracy

    Get PDF
    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

    Full text link
    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

    Full text link
    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

    Full text link
    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

    Get PDF
    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

    Full text link
    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 (MOLT^T) formulation, combined with a fast summation method with computational complexity O(Nlog⁥N)O(N\log{N}), where NN 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

    Full text link
    We examine the Boltzmann/Gibbs/Shannon SBGS\mathcal{S}_{BGS} and the non-additive Havrda-Charv\'{a}t / Dar\'{o}czy/Cressie-Read/Tsallis \ Sq\mathcal{S}_q \ and the Kaniadakis Îș\kappa-entropy \ SÎș\mathcal{S}_\kappa \ 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
    • 

    corecore