Volume 46, Number 2, November-December 2012
|Page(s)||443 - 463|
|Published online||24 October 2011|
A numerical scheme for the quantum Boltzmann equation with stiff collision terms⋆
Universitéde Lyon, Université Lyon I, CNRS UMR 5208, Institut
Camille Jordan, 43 boulevard du 11
Novembre 1918, 69622
2 Department of Mathematics, University of Wisconsin-Madison, 480 Lincoln Drive, Madison, 53706 WI, USA ;
3 Institute for Computational Engineering and Sciences (ICES), The University of Texas at Austin, 1 University Station C0200, Austin, 78712 TX, USA
Revised: 20 June 2011
Numerically solving the Boltzmann kinetic equations with the small Knudsen number is challenging due to the stiff nonlinear collision terms. A class of asymptotic-preserving schemes 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 collision term by a BGK type operator. This method, however, encounters its own difficulty when applied to the quantum Boltzmann equation. To define the quantum Maxwellian (Bose-Einstein or Fermi-Dirac distribution) at each time step and every mesh point, one has to invert a nonlinear equation that connects the macroscopic quantity fugacity with density and internal energy. Setting a good initial guess for the iterative method is troublesome in most cases because of the complexity of the quantum functions (Bose-Einstein or Fermi-Dirac function). In this paper, we propose to penalize the quantum collision term by a ‘classical’ BGK operator instead of the quantum one. This is based on the observation that the classical Maxwellian, with the temperature replaced by the internal energy, has the same first five moments as the quantum Maxwellian. The scheme so designed avoids the aforementioned difficulty, and one can show that the density distribution is still driven toward the quantum equilibrium. Numerical results are presented to illustrate the efficiency of the new scheme in both the hydrodynamic and kinetic regimes. We also develop a spectral method for the quantum collision operator.
Mathematics Subject Classification: 35Q20 / 65L04 / 76Y05
Key words: Quantum Boltzmann equation / Bose/Fermi gas / asymptotic-preserving schemes / fluid dynamic limit
This work was partially supported by NSF grant DMS-0608720 and NSF FRG grant DMS-0757285. FF was also supported by the ERC Starting Grant Project NuSiKiMo, project 239983-NuSiKiMo. SJ was also supported by a Van Vleck Distinguished Research Prize and a Vilas Associate Award from the University of Wisconsin-Madison.
© EDP Sciences, SMAI, 2011
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.