In a recent paper we presented a new ultra efficient numerical method for solving kinetic equations of the Boltzmann type (Dimarco and Loubere, 2013) . The key idea, on which the method relies, is to solve the collision part on a grid and then to solve exactly the transport part by following the characteristics backward in time. On the contrary to classical semi-Lagrangian methods one does not need to reconstruct the distribution function at each time step. This allows to tremendously reduce the computational cost and to perform efficient numerical simulations of kinetic equations up to the six-dimensional case without parallelization. However, the main drawback of the method developed was the loss of spatial accuracy close to the fluid limit. In the present work, we modify the scheme in such a way that it is able to preserve the high order spatial accuracy for extremely rarefied and fluid regimes. In particular, in the fluid limit, the method automatically degenerates into a high order method for the compressible Euler equations. Numerical examples are presented which validate the method, show the higher accuracy with respect to the previous approach and measure its efficiency with respect to well-known schemes (Direct Simulation Monte Carlo, Finite Volume, MUSCL, WENO).

Towards an ultra efficient kinetic scheme. Part II: The high order case

DIMARCO, Giacomo
Primo
;
2013

Abstract

In a recent paper we presented a new ultra efficient numerical method for solving kinetic equations of the Boltzmann type (Dimarco and Loubere, 2013) . The key idea, on which the method relies, is to solve the collision part on a grid and then to solve exactly the transport part by following the characteristics backward in time. On the contrary to classical semi-Lagrangian methods one does not need to reconstruct the distribution function at each time step. This allows to tremendously reduce the computational cost and to perform efficient numerical simulations of kinetic equations up to the six-dimensional case without parallelization. However, the main drawback of the method developed was the loss of spatial accuracy close to the fluid limit. In the present work, we modify the scheme in such a way that it is able to preserve the high order spatial accuracy for extremely rarefied and fluid regimes. In particular, in the fluid limit, the method automatically degenerates into a high order method for the compressible Euler equations. Numerical examples are presented which validate the method, show the higher accuracy with respect to the previous approach and measure its efficiency with respect to well-known schemes (Direct Simulation Monte Carlo, Finite Volume, MUSCL, WENO).
2013
Dimarco, Giacomo; R., Loubere
File in questo prodotto:
File Dimensione Formato  
JCP1.pdf

solo gestori archivio

Descrizione: versione editoriale
Tipologia: Full text (versione editoriale)
Licenza: NON PUBBLICO - Accesso privato/ristretto
Dimensione 1.56 MB
Formato Adobe PDF
1.56 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
1212.0204.pdf

accesso aperto

Descrizione: versione pre print
Tipologia: Full text (versione editoriale)
Licenza: PUBBLICO - Pubblico con Copyright
Dimensione 1.79 MB
Formato Adobe PDF
1.79 MB Adobe PDF Visualizza/Apri

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/1937412
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 19
  • ???jsp.display-item.citation.isi??? 19
social impact