An improved Quiet Direct Simulation method for Eulerian fluids using a second-order scheme

Matt-Hew Smith, H. M. Cave, J. S. Wu, M. C. Jermy, Y. S. Chen

Research output: Contribution to journalArticle

27 Citations (Scopus)

Abstract

In this paper, a second-order scheme for the Quiet Direct Simulation (QDS) of Eulerian fluids is proposed. The QDS method replaces the random sampling method used in Direct Simulation Monte Carlo (DSMC) methods with a technique whereby particles are moved, have their properties distributed onto a mesh, are destroyed and then are recreated deterministically from the properties stored on the mesh using Gauss-Hermite quadrature weights and abscissas. Particles are permitted to move in physically realistic directions so flux exchange is not limited to cells sharing an adjacent interface as in conventional, direction decoupled finite volume solvers. In this paper the method is extended by calculating the fluxes of mass, momentum and energy between cells assuming a linear variation of density, temperature and velocity in each cell and using these fluxes to update the mass, velocity and internal energy carried by each particle. This Euler solver has several advantages including large dynamic range, no statistical scatter in the results, true direction fluxes to all nearby neighbors and is computationally inexpensive. The second-order method is found to reduce the numerical diffusion of QDS as demonstrated in several verification studies. These include unsteady shock tube flow, a two-dimensional blast wave and of the development of Mach 3 flow over a forward facing step in a wind tunnel, which are compared with previous results from the literature wherever is possible. Finally the implementation of QUIETWAVE, a rapid method of simulating blast events in urban environments, is introduced and the results of a test case are presented.

Original languageEnglish
Pages (from-to)2213-2224
Number of pages12
JournalJournal of Computational Physics
Volume228
Issue number6
DOIs
Publication statusPublished - 2009 Apr 1

Fingerprint

Fluxes
Fluids
fluids
blasts
mesh
forward facing steps
simulation
cells
random sampling
Shock tubes
Pipe flow
shock tubes
wind tunnels
internal energy
quadratures
Mach number
Wind tunnels
dynamic range
Monte Carlo method
energy

All Science Journal Classification (ASJC) codes

  • Numerical Analysis
  • Modelling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Cite this

Smith, Matt-Hew ; Cave, H. M. ; Wu, J. S. ; Jermy, M. C. ; Chen, Y. S. / An improved Quiet Direct Simulation method for Eulerian fluids using a second-order scheme. In: Journal of Computational Physics. 2009 ; Vol. 228, No. 6. pp. 2213-2224.
@article{0a542629eb064075b26b76c3427ecd22,
title = "An improved Quiet Direct Simulation method for Eulerian fluids using a second-order scheme",
abstract = "In this paper, a second-order scheme for the Quiet Direct Simulation (QDS) of Eulerian fluids is proposed. The QDS method replaces the random sampling method used in Direct Simulation Monte Carlo (DSMC) methods with a technique whereby particles are moved, have their properties distributed onto a mesh, are destroyed and then are recreated deterministically from the properties stored on the mesh using Gauss-Hermite quadrature weights and abscissas. Particles are permitted to move in physically realistic directions so flux exchange is not limited to cells sharing an adjacent interface as in conventional, direction decoupled finite volume solvers. In this paper the method is extended by calculating the fluxes of mass, momentum and energy between cells assuming a linear variation of density, temperature and velocity in each cell and using these fluxes to update the mass, velocity and internal energy carried by each particle. This Euler solver has several advantages including large dynamic range, no statistical scatter in the results, true direction fluxes to all nearby neighbors and is computationally inexpensive. The second-order method is found to reduce the numerical diffusion of QDS as demonstrated in several verification studies. These include unsteady shock tube flow, a two-dimensional blast wave and of the development of Mach 3 flow over a forward facing step in a wind tunnel, which are compared with previous results from the literature wherever is possible. Finally the implementation of QUIETWAVE, a rapid method of simulating blast events in urban environments, is introduced and the results of a test case are presented.",
author = "Matt-Hew Smith and Cave, {H. M.} and Wu, {J. S.} and Jermy, {M. C.} and Chen, {Y. S.}",
year = "2009",
month = "4",
day = "1",
doi = "10.1016/j.jcp.2008.12.013",
language = "English",
volume = "228",
pages = "2213--2224",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "6",

}

An improved Quiet Direct Simulation method for Eulerian fluids using a second-order scheme. / Smith, Matt-Hew; Cave, H. M.; Wu, J. S.; Jermy, M. C.; Chen, Y. S.

In: Journal of Computational Physics, Vol. 228, No. 6, 01.04.2009, p. 2213-2224.

Research output: Contribution to journalArticle

TY - JOUR

T1 - An improved Quiet Direct Simulation method for Eulerian fluids using a second-order scheme

AU - Smith, Matt-Hew

AU - Cave, H. M.

AU - Wu, J. S.

AU - Jermy, M. C.

AU - Chen, Y. S.

PY - 2009/4/1

Y1 - 2009/4/1

N2 - In this paper, a second-order scheme for the Quiet Direct Simulation (QDS) of Eulerian fluids is proposed. The QDS method replaces the random sampling method used in Direct Simulation Monte Carlo (DSMC) methods with a technique whereby particles are moved, have their properties distributed onto a mesh, are destroyed and then are recreated deterministically from the properties stored on the mesh using Gauss-Hermite quadrature weights and abscissas. Particles are permitted to move in physically realistic directions so flux exchange is not limited to cells sharing an adjacent interface as in conventional, direction decoupled finite volume solvers. In this paper the method is extended by calculating the fluxes of mass, momentum and energy between cells assuming a linear variation of density, temperature and velocity in each cell and using these fluxes to update the mass, velocity and internal energy carried by each particle. This Euler solver has several advantages including large dynamic range, no statistical scatter in the results, true direction fluxes to all nearby neighbors and is computationally inexpensive. The second-order method is found to reduce the numerical diffusion of QDS as demonstrated in several verification studies. These include unsteady shock tube flow, a two-dimensional blast wave and of the development of Mach 3 flow over a forward facing step in a wind tunnel, which are compared with previous results from the literature wherever is possible. Finally the implementation of QUIETWAVE, a rapid method of simulating blast events in urban environments, is introduced and the results of a test case are presented.

AB - In this paper, a second-order scheme for the Quiet Direct Simulation (QDS) of Eulerian fluids is proposed. The QDS method replaces the random sampling method used in Direct Simulation Monte Carlo (DSMC) methods with a technique whereby particles are moved, have their properties distributed onto a mesh, are destroyed and then are recreated deterministically from the properties stored on the mesh using Gauss-Hermite quadrature weights and abscissas. Particles are permitted to move in physically realistic directions so flux exchange is not limited to cells sharing an adjacent interface as in conventional, direction decoupled finite volume solvers. In this paper the method is extended by calculating the fluxes of mass, momentum and energy between cells assuming a linear variation of density, temperature and velocity in each cell and using these fluxes to update the mass, velocity and internal energy carried by each particle. This Euler solver has several advantages including large dynamic range, no statistical scatter in the results, true direction fluxes to all nearby neighbors and is computationally inexpensive. The second-order method is found to reduce the numerical diffusion of QDS as demonstrated in several verification studies. These include unsteady shock tube flow, a two-dimensional blast wave and of the development of Mach 3 flow over a forward facing step in a wind tunnel, which are compared with previous results from the literature wherever is possible. Finally the implementation of QUIETWAVE, a rapid method of simulating blast events in urban environments, is introduced and the results of a test case are presented.

UR - http://www.scopus.com/inward/record.url?scp=59149095608&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=59149095608&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2008.12.013

DO - 10.1016/j.jcp.2008.12.013

M3 - Article

AN - SCOPUS:59149095608

VL - 228

SP - 2213

EP - 2224

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 6

ER -