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 long 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 \begin{document}$ 3 $\end{document} D form a condensate in finite time and converge in entropy to the unique minimiser of the natural entropy functional. 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.