An approximate method for solving rarefied and transitional flows using TDEFM with isotropic mesh adaptation

Matt-Hew Smith, H. M. Cave, J. S. Wu, M. N. Macrossan

Research output: Contribution to journalConference article

Abstract

DSMC [1] can become increasingly expensive when extended to the near-continuum regime. Because of the statistical nature of the results, long run times are required to build up samples of simulator particles large enough to reduce the statistical scatter to acceptable levels. Here we adapt a kinetic theory based flux method to produce a quick approximate solver for transition and near-continuum flows. The results have no statistical scatter. The CPU times are similar to those of traditional continuum (Navier-Stokes or Euler) solvers. The True Direction Equilibrium Flux Method (TDEFM) [2, 3] is a generalisation of Pullin's kinetic theory based EFM [4]. TDEFM can transfer fluxes of mass, momentum and energy in physically realistic directions from any source cell to any destination cell, even if the cells do not share an interface. TDEFM, as an Euler solver, has been shown to provide good results on a Cartesian grid for flows where standard continuum methods produce unphysical asymmetries apparently because the continuum fluxes are constrained (in one time step) to flow in the grid coordinate directions rather than the correct physical direction. [2, 3] The new method for rarefied flow does not try to produce the correct velocity distribution function, but does ensure that mass, momentum and energy are transported within the flow over the physically correct distances between "pseudo-collisions". To ensure this, (1) the time step is restricted so that mass, momentum and energy are exchanged between contiguous cells only in one time step, and (2) the cells sizes are adapted, as steady state is approached, to be approximately equal to the local mean free path. The results for Mach 5 flow over a flat plate for varying Knudsen numbers show an average difference (compared to DSMC) in the X-velocity profile near the surface of the plate of less than 6 percent. TDEFM, employing adaptive mesh refinement, required less than 9 percent of the computational time required by DSMC for the same flow. Thus the approximate method could be useful for quick "first-estimate" solutions of otherwise time consuming design problems.

Original languageEnglish
Pages (from-to)371-376
Number of pages6
JournalAIP Conference Proceedings
Volume1084
Publication statusPublished - 2009 Apr 13
Event26th International Symposium on Rarefied Gas Dynamics, RGD26 - Kyoto, Japan
Duration: 2008 Jul 202008 Jul 25

Fingerprint

mesh
continuums
cells
kinetic theory
momentum
velocity distribution
continuum flow
grids
Knudsen flow
flat plates
mean free path
simulators
energy
distribution functions
asymmetry
collisions
estimates

All Science Journal Classification (ASJC) codes

  • Physics and Astronomy(all)

Cite this

Smith, Matt-Hew ; Cave, H. M. ; Wu, J. S. ; Macrossan, M. N. / An approximate method for solving rarefied and transitional flows using TDEFM with isotropic mesh adaptation. In: AIP Conference Proceedings. 2009 ; Vol. 1084. pp. 371-376.
@article{888b972984424f1f9d4ade2cb63daa3f,
title = "An approximate method for solving rarefied and transitional flows using TDEFM with isotropic mesh adaptation",
abstract = "DSMC [1] can become increasingly expensive when extended to the near-continuum regime. Because of the statistical nature of the results, long run times are required to build up samples of simulator particles large enough to reduce the statistical scatter to acceptable levels. Here we adapt a kinetic theory based flux method to produce a quick approximate solver for transition and near-continuum flows. The results have no statistical scatter. The CPU times are similar to those of traditional continuum (Navier-Stokes or Euler) solvers. The True Direction Equilibrium Flux Method (TDEFM) [2, 3] is a generalisation of Pullin's kinetic theory based EFM [4]. TDEFM can transfer fluxes of mass, momentum and energy in physically realistic directions from any source cell to any destination cell, even if the cells do not share an interface. TDEFM, as an Euler solver, has been shown to provide good results on a Cartesian grid for flows where standard continuum methods produce unphysical asymmetries apparently because the continuum fluxes are constrained (in one time step) to flow in the grid coordinate directions rather than the correct physical direction. [2, 3] The new method for rarefied flow does not try to produce the correct velocity distribution function, but does ensure that mass, momentum and energy are transported within the flow over the physically correct distances between {"}pseudo-collisions{"}. To ensure this, (1) the time step is restricted so that mass, momentum and energy are exchanged between contiguous cells only in one time step, and (2) the cells sizes are adapted, as steady state is approached, to be approximately equal to the local mean free path. The results for Mach 5 flow over a flat plate for varying Knudsen numbers show an average difference (compared to DSMC) in the X-velocity profile near the surface of the plate of less than 6 percent. TDEFM, employing adaptive mesh refinement, required less than 9 percent of the computational time required by DSMC for the same flow. Thus the approximate method could be useful for quick {"}first-estimate{"} solutions of otherwise time consuming design problems.",
author = "Matt-Hew Smith and Cave, {H. M.} and Wu, {J. S.} and Macrossan, {M. N.}",
year = "2009",
month = "4",
day = "13",
language = "English",
volume = "1084",
pages = "371--376",
journal = "AIP Conference Proceedings",
issn = "0094-243X",
publisher = "American Institute of Physics Publising LLC",

}

An approximate method for solving rarefied and transitional flows using TDEFM with isotropic mesh adaptation. / Smith, Matt-Hew; Cave, H. M.; Wu, J. S.; Macrossan, M. N.

In: AIP Conference Proceedings, Vol. 1084, 13.04.2009, p. 371-376.

Research output: Contribution to journalConference article

TY - JOUR

T1 - An approximate method for solving rarefied and transitional flows using TDEFM with isotropic mesh adaptation

AU - Smith, Matt-Hew

AU - Cave, H. M.

AU - Wu, J. S.

AU - Macrossan, M. N.

PY - 2009/4/13

Y1 - 2009/4/13

N2 - DSMC [1] can become increasingly expensive when extended to the near-continuum regime. Because of the statistical nature of the results, long run times are required to build up samples of simulator particles large enough to reduce the statistical scatter to acceptable levels. Here we adapt a kinetic theory based flux method to produce a quick approximate solver for transition and near-continuum flows. The results have no statistical scatter. The CPU times are similar to those of traditional continuum (Navier-Stokes or Euler) solvers. The True Direction Equilibrium Flux Method (TDEFM) [2, 3] is a generalisation of Pullin's kinetic theory based EFM [4]. TDEFM can transfer fluxes of mass, momentum and energy in physically realistic directions from any source cell to any destination cell, even if the cells do not share an interface. TDEFM, as an Euler solver, has been shown to provide good results on a Cartesian grid for flows where standard continuum methods produce unphysical asymmetries apparently because the continuum fluxes are constrained (in one time step) to flow in the grid coordinate directions rather than the correct physical direction. [2, 3] The new method for rarefied flow does not try to produce the correct velocity distribution function, but does ensure that mass, momentum and energy are transported within the flow over the physically correct distances between "pseudo-collisions". To ensure this, (1) the time step is restricted so that mass, momentum and energy are exchanged between contiguous cells only in one time step, and (2) the cells sizes are adapted, as steady state is approached, to be approximately equal to the local mean free path. The results for Mach 5 flow over a flat plate for varying Knudsen numbers show an average difference (compared to DSMC) in the X-velocity profile near the surface of the plate of less than 6 percent. TDEFM, employing adaptive mesh refinement, required less than 9 percent of the computational time required by DSMC for the same flow. Thus the approximate method could be useful for quick "first-estimate" solutions of otherwise time consuming design problems.

AB - DSMC [1] can become increasingly expensive when extended to the near-continuum regime. Because of the statistical nature of the results, long run times are required to build up samples of simulator particles large enough to reduce the statistical scatter to acceptable levels. Here we adapt a kinetic theory based flux method to produce a quick approximate solver for transition and near-continuum flows. The results have no statistical scatter. The CPU times are similar to those of traditional continuum (Navier-Stokes or Euler) solvers. The True Direction Equilibrium Flux Method (TDEFM) [2, 3] is a generalisation of Pullin's kinetic theory based EFM [4]. TDEFM can transfer fluxes of mass, momentum and energy in physically realistic directions from any source cell to any destination cell, even if the cells do not share an interface. TDEFM, as an Euler solver, has been shown to provide good results on a Cartesian grid for flows where standard continuum methods produce unphysical asymmetries apparently because the continuum fluxes are constrained (in one time step) to flow in the grid coordinate directions rather than the correct physical direction. [2, 3] The new method for rarefied flow does not try to produce the correct velocity distribution function, but does ensure that mass, momentum and energy are transported within the flow over the physically correct distances between "pseudo-collisions". To ensure this, (1) the time step is restricted so that mass, momentum and energy are exchanged between contiguous cells only in one time step, and (2) the cells sizes are adapted, as steady state is approached, to be approximately equal to the local mean free path. The results for Mach 5 flow over a flat plate for varying Knudsen numbers show an average difference (compared to DSMC) in the X-velocity profile near the surface of the plate of less than 6 percent. TDEFM, employing adaptive mesh refinement, required less than 9 percent of the computational time required by DSMC for the same flow. Thus the approximate method could be useful for quick "first-estimate" solutions of otherwise time consuming design problems.

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

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

M3 - Conference article

AN - SCOPUS:63849136263

VL - 1084

SP - 371

EP - 376

JO - AIP Conference Proceedings

JF - AIP Conference Proceedings

SN - 0094-243X

ER -