Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 86
PROCEEDINGS OF THE ELEVENTH INTERNATIONAL CONFERENCE ON CIVIL, STRUCTURAL AND ENVIRONMENTAL ENGINEERING COMPUTING Edited by: B.H.V. Topping
Paper 214
A Discontinuous TaylorGalerkin Procedure for Solution of the BoltzmannBGK Equation in Two Dimensions B. Evans, K. Morgan and O. Hassan
Civil & Computational Engineering Centre, Swansea University, United Kingdom B. Evans, K. Morgan, O. Hassan, "A Discontinuous TaylorGalerkin Procedure for Solution of the BoltzmannBGK Equation in Two Dimensions", in B.H.V. Topping, (Editor), "Proceedings of the Eleventh International Conference on Civil, Structural and Environmental Engineering Computing", CivilComp Press, Stirlingshire, UK, Paper 214, 2007. doi:10.4203/ccp.86.214
Keywords: discontinuous, TaylorGalerkin, BoltzmannBGK.
Summary
In recent decades, the computational fluid dynamics community has largely focused on the application of the Euler and NavierStokes equation systems in finding solutions to fluid flow problems. This has been based on the continuum assumption that allows macroscopic variables such as pressure and temperature to be clearly defined. There are conditions, however, under which this continuum assumption breaks down. For example, nanoscale flows, rarefied gas dynamics and hypersonics are applications in which it is possible that insufficient molecular collisions are taking place to allow the gas to maintain equilibrium and for the continuum assumption to hold [1].
Continuum validity can be expressed in terms of the Knudsen number (K_{n}) which is the ratio of mean free molecular path to some typical lengthscale of the flow and is a measure of the degree of rarefaction of the gas. Above K_{n}=0.1, we enter the transition regime in which the standard conservation equations no longer form a closed set due to our inability to express shear stresses and heat flux in terms of the lower order macroscopic variables. Once we reach the limit of validity of the continuum equations, we must take a step of abstraction back and solve the Boltzmann equation from which the Euler and NavierStokes systems are derived. In this paper we present a discontinuous TaylorGalerkin [2] twostep finite element procedure for solution of the BoltzmannBGK equation [3]. This is a simplified version of the full Boltzmann equation. The term which poses the greatest mathematical difficulty in the full equation is the molecular collision term on the righthand side and it is this term which is modified in the BGK version. This is done by making the assumption that the effect of molecular collisions is to force a nonequilibrium gas back to equilibrium over some timescale, t. The BoltzmannBGK solver allows us to analyse significantly nonequilibrium problems without contravening any of the assumptions that make our equation systems valid. In this paper we present some 2D gas flow examples including a shock tube problem, flow over a ridge and transonic flow over an aerofoil. Results are compared with analytical solutions where they exist and results from NavierStokes solvers and Direct Simulation Monte Carlo (DSMC) methods [1]. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
