A gas-kinetic BGK solver for two-dimensional turbulent compressible flow
In this paper, a gas kinetic solver is developed for the Reynolds Average Navier-Stokes (RANS) equations in two-space dimensions. To our best knowledge, this is the first attempt to extend the application of the BGK (Bhatnagaar-Gross-Krook) scheme to solve RANS equations with a turbulence model usin...
Main Authors: | , , , |
---|---|
Format: | Conference or Workshop Item |
Language: | English |
Published: |
2008
|
Subjects: | |
Online Access: | http://irep.iium.edu.my/6125/ http://irep.iium.edu.my/6125/1/BGKsolver.pdf |
id |
iium-6125 |
---|---|
recordtype |
eprints |
spelling |
iium-61252011-11-21T20:52:28Z http://irep.iium.edu.my/6125/ A gas-kinetic BGK solver for two-dimensional turbulent compressible flow Ong, Jiunn Chit Omar, Ashraf Ali Asrar, Waqar Ismail, Ahmad Faris TL500 Aeronautics In this paper, a gas kinetic solver is developed for the Reynolds Average Navier-Stokes (RANS) equations in two-space dimensions. To our best knowledge, this is the first attempt to extend the application of the BGK (Bhatnagaar-Gross-Krook) scheme to solve RANS equations with a turbulence model using finite difference method. The convection flux terms which appear on the left hand side of the RANS equations are discretized by a semi-discrete finite difference method. Then, the resulting inviscid flux functions are approximated by gas-kinetic BGK scheme which is based on the BGK model of the approximate collisional Boltzmann equation. The cell interface values required by the inviscid flux functions are reconstructed to higher-order spatial accuracy via the MUSCL (Monotone Upstream-Centered Schemes for Conservation Laws) variable interpolation method coupled with a minmod limiter. As for the diffusion flux terms, they are discretized by a second-order central difference scheme. To account for the turbulence effect, a combined k-ε / k-ω SST (Shear-Stress Transport) two-equation turbulence model is used in the solver. An explicit-type time integration method known as the modified fourth-order Runge-Kutta method is used to compute steady-state solutions. The computed results for a supersonic flow past a flat plate where the transition is artificially triggered at 50% of plate length are presented in this paper. Validating the computed results against existing analytical solutions and also comparing them with results from other well-known numerical schemes show that a very good agreement is obtained. 2008-07 Conference or Workshop Item PeerReviewed application/pdf en http://irep.iium.edu.my/6125/1/BGKsolver.pdf Ong, Jiunn Chit and Omar, Ashraf Ali and Asrar, Waqar and Ismail, Ahmad Faris (2008) A gas-kinetic BGK solver for two-dimensional turbulent compressible flow. In: International Conference on Scientific Computing , July 14-17, 2008, Las Vegas, Nevada, USA. |
repository_type |
Digital Repository |
institution_category |
Local University |
institution |
International Islamic University Malaysia |
building |
IIUM Repository |
collection |
Online Access |
language |
English |
topic |
TL500 Aeronautics |
spellingShingle |
TL500 Aeronautics Ong, Jiunn Chit Omar, Ashraf Ali Asrar, Waqar Ismail, Ahmad Faris A gas-kinetic BGK solver for two-dimensional turbulent compressible flow |
description |
In this paper, a gas kinetic solver is developed for the Reynolds Average Navier-Stokes (RANS) equations in two-space dimensions. To our best knowledge, this is the first attempt to extend the application of the BGK (Bhatnagaar-Gross-Krook) scheme to solve RANS equations with a turbulence model using finite difference method. The convection flux terms which appear on the left hand side of the RANS equations are discretized by a semi-discrete finite difference method. Then, the resulting inviscid flux functions are approximated by gas-kinetic BGK scheme which is based on the BGK model of the approximate collisional Boltzmann equation. The cell interface values required by the inviscid flux functions are reconstructed to higher-order spatial accuracy via the MUSCL (Monotone Upstream-Centered Schemes for Conservation Laws) variable interpolation method coupled with a minmod limiter. As for the diffusion flux terms, they are discretized by a second-order central difference scheme. To account for the turbulence effect, a combined k-ε / k-ω SST (Shear-Stress Transport) two-equation turbulence model is used in the solver. An explicit-type time integration method known as the modified fourth-order Runge-Kutta method is used to compute steady-state solutions. The computed results for a supersonic flow past a flat plate where the transition is artificially triggered at 50% of plate length are presented in this paper. Validating the computed results against existing analytical solutions and also comparing them with results from other well-known numerical schemes show that a very good agreement is obtained. |
format |
Conference or Workshop Item |
author |
Ong, Jiunn Chit Omar, Ashraf Ali Asrar, Waqar Ismail, Ahmad Faris |
author_facet |
Ong, Jiunn Chit Omar, Ashraf Ali Asrar, Waqar Ismail, Ahmad Faris |
author_sort |
Ong, Jiunn Chit |
title |
A gas-kinetic BGK solver for two-dimensional turbulent compressible flow |
title_short |
A gas-kinetic BGK solver for two-dimensional turbulent compressible flow |
title_full |
A gas-kinetic BGK solver for two-dimensional turbulent compressible flow |
title_fullStr |
A gas-kinetic BGK solver for two-dimensional turbulent compressible flow |
title_full_unstemmed |
A gas-kinetic BGK solver for two-dimensional turbulent compressible flow |
title_sort |
gas-kinetic bgk solver for two-dimensional turbulent compressible flow |
publishDate |
2008 |
url |
http://irep.iium.edu.my/6125/ http://irep.iium.edu.my/6125/1/BGKsolver.pdf |
first_indexed |
2023-09-18T20:14:57Z |
last_indexed |
2023-09-18T20:14:57Z |
_version_ |
1777407725213319168 |