In this work we present a non stationary domain decomposition algorithm for multiscale hydrodynamic-kinetic problems, in which the Knudsen number may span from equilib- rium to highly rarefied regimes. Our approach is characterized by using the full Boltzmann equation for the kinetic regime, the Compressible Euler equations for equilibrium, with a buffer zone in which the BGK-ES equation is used to represent the transition between fully kinetic to equilibrium flows. In this fashion, the Boltzmann solver is used only when the collision integral is non-stiff, and the mean free path is of the same order as the mesh size needed to capture variations in macroscopic quantities. Thus, in principle, the same mesh size and time steps can be used in the whole computation. Moreover, the time step is limited only by convective terms. Since the Boltzmann solver is applied only in wholly kinetic regimes, we use the reduced noise DSMC scheme we have proposed in Part I of the present work. This ensures a smooth exchange of information across the different domains, with a natural way to construct interface numerical fluxes. Several tests comparing our hybrid scheme with full Boltzmann DSMC computations show the good agreement between the two solutions, on a wide range of Knudsen numbers.

A hybrid method for hydrodynamic-kinetic flow – Part II – Coupling of hydrodynamic and kinetic models / Alaia, A.; PUPPO, GABRIELLA ANNA. - In: JOURNAL OF COMPUTATIONAL PHYSICS. - ISSN 0021-9991. - 231:(2012), pp. 5217-5242. [10.1016/j.jcp.2012.02.022]

A hybrid method for hydrodynamic-kinetic flow – Part II – Coupling of hydrodynamic and kinetic models

PUPPO, GABRIELLA ANNA
2012

Abstract

In this work we present a non stationary domain decomposition algorithm for multiscale hydrodynamic-kinetic problems, in which the Knudsen number may span from equilib- rium to highly rarefied regimes. Our approach is characterized by using the full Boltzmann equation for the kinetic regime, the Compressible Euler equations for equilibrium, with a buffer zone in which the BGK-ES equation is used to represent the transition between fully kinetic to equilibrium flows. In this fashion, the Boltzmann solver is used only when the collision integral is non-stiff, and the mean free path is of the same order as the mesh size needed to capture variations in macroscopic quantities. Thus, in principle, the same mesh size and time steps can be used in the whole computation. Moreover, the time step is limited only by convective terms. Since the Boltzmann solver is applied only in wholly kinetic regimes, we use the reduced noise DSMC scheme we have proposed in Part I of the present work. This ensures a smooth exchange of information across the different domains, with a natural way to construct interface numerical fluxes. Several tests comparing our hybrid scheme with full Boltzmann DSMC computations show the good agreement between the two solutions, on a wide range of Knudsen numbers.
2012
Numerical methods for kinetic flows, domain decomposition strategies, hybrid schemes
01 Pubblicazione su rivista::01a Articolo in rivista
A hybrid method for hydrodynamic-kinetic flow – Part II – Coupling of hydrodynamic and kinetic models / Alaia, A.; PUPPO, GABRIELLA ANNA. - In: JOURNAL OF COMPUTATIONAL PHYSICS. - ISSN 0021-9991. - 231:(2012), pp. 5217-5242. [10.1016/j.jcp.2012.02.022]
File allegati a questo prodotto
File Dimensione Formato  
Alaia_A-hybrid-method_2012.pdf

solo gestori archivio

Tipologia: Versione editoriale (versione pubblicata con il layout dell'editore)
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 1.61 MB
Formato Adobe PDF
1.61 MB Adobe PDF   Contatta l'autore

I documenti in IRIS 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/11573/1280452
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 16
  • ???jsp.display-item.citation.isi??? 14
social impact