Numerical study of Bose-Einstein condensation in the Kaniadakis-Quarati model for bosons
arXiv:1902.06266 · doi:10.3934/krm.2020017
Abstract
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.
arXiv commentTo appear in Kinet. Relat. Models
References in corpus (11)
- Mathematical theory and numerical methods for Bose-Einstein condensation
- A Finite-Volume Method for Nonlinear Nonlocal Equations with a Gradient Flow Structure
- Kinetics of Bose-Condensation
- Condensation of bosons in kinetic regime
- Convergence of a variational Lagrangian scheme for a nonlinear drift diffusion equation
- Dynamics of the Bose-Einstein condensation: analogy with the collapse dynamics of a classical self-gravitating Brownian gas
- Numerical simulation of nonlinear continuity equations by evolving diffeomorphisms
- A Lagrangian scheme for the solution of nonlinear diffusion equations using moving simplex meshes
- The Fokker-Planck equation for bosons in 2D: well-possedness and asymptotic analysis
- On the singularity formation and relaxation to equilibrium in 1D Fokker-Planck model with superlinear drift
- Blow-up rate estimates for the solutions of the bosonic Boltzmann-Nordheim equation