Hostname: page-component-7b7567c8c-mnklg Total loading time: 0 Render date: 2026-11-12T19:33:52.906Z Has data issue: false hasContentIssue false

A numerical scheme for the quantum Boltzmann equation withstiff collision terms⋆

Published online by Cambridge University Press:  24 October 2011

Francis Filbet
Affiliation:
Universitéde Lyon, Université Lyon I, CNRS UMR 5208, Institut Camille Jordan, 43 boulevard du 11 Novembre 1918, 69622 Villeurbanne cedex, France. filbet@math.univ-lyon1.fr
Jingwei Hu
Affiliation:
Department of Mathematics, University of Wisconsin-Madison, 480 Lincoln Drive, Madison, 53706 WI, USA; hu@ices.utexas.edu ; jin@math.wisc.edu ;Institute for Computational Engineering and Sciences (ICES), The University of Texas at Austin, 1 University Station C0200, Austin, 78712 TX, USA
Shi Jin
Affiliation:
Department of Mathematics, University of Wisconsin-Madison, 480 Lincoln Drive, Madison, 53706 WI, USA; hu@ices.utexas.edu ; jin@math.wisc.edu ;

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the 'Save PDF' action button.

Numerically solving the Boltzmann kinetic equations with the small Knudsen number ischallenging due to the stiff nonlinear collision terms. A class of asymptotic-preservingschemes was introduced in [F. Filbet and S. Jin,J. Comput. Phys. 229 (2010)7625–7648] to handle this kind of problems. The idea is to penalize the stiff collisionterm by a BGK type operator. This method, however, encounters its own difficulty whenapplied to the quantum Boltzmann equation. To define the quantum Maxwellian (Bose-Einsteinor Fermi-Dirac distribution) at each time step and every mesh point, one has to invert anonlinear equation that connects the macroscopic quantity fugacity with density andinternal energy. Setting a good initial guess for the iterative method is troublesome inmost cases because of the complexity of the quantum functions (Bose-Einstein orFermi-Dirac function). In this paper, we propose to penalize the quantum collision term bya ‘classical’ BGK operator instead of the quantum one. This is based on the observationthat the classical Maxwellian, with the temperature replaced by the internal energy, hasthe same first five moments as the quantum Maxwellian. The scheme so designed avoids theaforementioned difficulty, and one can show that the density distribution is still driventoward the quantum equilibrium. Numerical results are presented to illustrate theefficiency of the new scheme in both the hydrodynamic and kinetic regimes. We also developa spectral method for the quantum collision operator.

Information

Type
Research Article
Copyright
© EDP Sciences, SMAI, 2011