For the space homogeneous Boltzmann equation, we formulate a hybrid Monte Carlo method that is robust in the fluid dynamic limit. This method is based on an analytic representation of the solution over a single time step and involves implicit time differencing derived from a suitable power series expansion of the solution (a generalized Wild expansion). A class of implicit, yet explicitly implementable, numerical schemes is obtained by substituting a Maxwellian distribution in place of the high order terms in the expansion. The numerical solution is represented as a convex combination of a non-equilibrium particle distribution and a Maxwellian. The hybrid distribution is then evolved by Monte Carlo using the implicit formulation for the time evolution. Computational simulations of spatially homogeneous problems by our method are presented here for the Kac model and for the variable hard sphere model (including Maxwell molecules). Comparison to exact solutions and to direct simulation Monte Carlo (DSMC) computations shows the robustness and the efficiency of the new method.

An implicit Monte arlo method for rarefied gas dynamics - 1. The space homogeneous case

PARESCHI, Lorenzo;
1999

Abstract

For the space homogeneous Boltzmann equation, we formulate a hybrid Monte Carlo method that is robust in the fluid dynamic limit. This method is based on an analytic representation of the solution over a single time step and involves implicit time differencing derived from a suitable power series expansion of the solution (a generalized Wild expansion). A class of implicit, yet explicitly implementable, numerical schemes is obtained by substituting a Maxwellian distribution in place of the high order terms in the expansion. The numerical solution is represented as a convex combination of a non-equilibrium particle distribution and a Maxwellian. The hybrid distribution is then evolved by Monte Carlo using the implicit formulation for the time evolution. Computational simulations of spatially homogeneous problems by our method are presented here for the Kac model and for the variable hard sphere model (including Maxwell molecules). Comparison to exact solutions and to direct simulation Monte Carlo (DSMC) computations shows the robustness and the efficiency of the new method.
1999
Pareschi, Lorenzo; Caflisch, Re
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in SFERA sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11392/1207320
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 23
  • ???jsp.display-item.citation.isi??? 53
social impact