Stable computation of turbulent flows with a low-Reynolds-number k-∈ turbulence model and explicit solver

Research output: Contribution to journalArticle

15 Citations (Scopus)

Abstract

It is well known that the incorporation of two-equation turbulence models into explicit time-marching schemes, especially the low Reynolds versions, encounters various numerical difficulties. This paper sets out to review the numerical techniques presently used and to analyse the problems encountered. A numerical approach is proposed for the implementation of the low-Reynolds-number (LRN) k-∈ models into an explicit Runge-Kutta time-marching solver. In this method, a semi-implicit treatment of the source terms of the k and ∈ equations is introduced, and modified residuals of the two equations are used in order to adopt implicit residual smoothing. These techniques ensure the positivity of k and ∈ in the computation, even with uniform initial conditions and local time stepping. The asymptotic behaviour of the turbulence quantities during time marching is analysed, providing insight into the stability and convergence characteristics of the k and ∈ equations using explicit and implicit schemes with and without the damping terms. It is found that the time step size has a strong influence on the evolution of the turbulence quantities and the stability of the solution, whether implicit or explicit schemes are used. The time step size determined by the stability analysis is a function of k, ∈ and their damping terms. This analysis also provides a guide for the selection of time step size for the turbulence equations. It is found that the selection of time step size is an important factor in ensuring realistic values of k, ∈ and μt in the flow field during time marching from uniform initial conditions. Numerical solutions of several flow test cases are given and the computed results are compared with experimental measurements or theoretical solutions for validation of the method. It is shown that this method is reasonably efficient and robust.

Original languageEnglish
Pages (from-to)487-499
Number of pages13
JournalAdvances in Engineering Software
Volume28
Issue number8
Publication statusPublished - Nov 1997
Externally publishedYes

Fingerprint

Turbulence models
Turbulent flow
Reynolds number
Turbulence
Damping
Flow fields

ASJC Scopus subject areas

  • Software
  • Computational Mechanics

Cite this

@article{20bbf4118e534e1d8fd86cf19bfec907,
title = "Stable computation of turbulent flows with a low-Reynolds-number k-∈ turbulence model and explicit solver",
abstract = "It is well known that the incorporation of two-equation turbulence models into explicit time-marching schemes, especially the low Reynolds versions, encounters various numerical difficulties. This paper sets out to review the numerical techniques presently used and to analyse the problems encountered. A numerical approach is proposed for the implementation of the low-Reynolds-number (LRN) k-∈ models into an explicit Runge-Kutta time-marching solver. In this method, a semi-implicit treatment of the source terms of the k and ∈ equations is introduced, and modified residuals of the two equations are used in order to adopt implicit residual smoothing. These techniques ensure the positivity of k and ∈ in the computation, even with uniform initial conditions and local time stepping. The asymptotic behaviour of the turbulence quantities during time marching is analysed, providing insight into the stability and convergence characteristics of the k and ∈ equations using explicit and implicit schemes with and without the damping terms. It is found that the time step size has a strong influence on the evolution of the turbulence quantities and the stability of the solution, whether implicit or explicit schemes are used. The time step size determined by the stability analysis is a function of k, ∈ and their damping terms. This analysis also provides a guide for the selection of time step size for the turbulence equations. It is found that the selection of time step size is an important factor in ensuring realistic values of k, ∈ and μt in the flow field during time marching from uniform initial conditions. Numerical solutions of several flow test cases are given and the computed results are compared with experimental measurements or theoretical solutions for validation of the method. It is shown that this method is reasonably efficient and robust.",
author = "Yong Zhao",
year = "1997",
month = "11",
language = "English",
volume = "28",
pages = "487--499",
journal = "Advances in Engineering Software",
issn = "0965-9978",
publisher = "Elsevier",
number = "8",

}

TY - JOUR

T1 - Stable computation of turbulent flows with a low-Reynolds-number k-∈ turbulence model and explicit solver

AU - Zhao, Yong

PY - 1997/11

Y1 - 1997/11

N2 - It is well known that the incorporation of two-equation turbulence models into explicit time-marching schemes, especially the low Reynolds versions, encounters various numerical difficulties. This paper sets out to review the numerical techniques presently used and to analyse the problems encountered. A numerical approach is proposed for the implementation of the low-Reynolds-number (LRN) k-∈ models into an explicit Runge-Kutta time-marching solver. In this method, a semi-implicit treatment of the source terms of the k and ∈ equations is introduced, and modified residuals of the two equations are used in order to adopt implicit residual smoothing. These techniques ensure the positivity of k and ∈ in the computation, even with uniform initial conditions and local time stepping. The asymptotic behaviour of the turbulence quantities during time marching is analysed, providing insight into the stability and convergence characteristics of the k and ∈ equations using explicit and implicit schemes with and without the damping terms. It is found that the time step size has a strong influence on the evolution of the turbulence quantities and the stability of the solution, whether implicit or explicit schemes are used. The time step size determined by the stability analysis is a function of k, ∈ and their damping terms. This analysis also provides a guide for the selection of time step size for the turbulence equations. It is found that the selection of time step size is an important factor in ensuring realistic values of k, ∈ and μt in the flow field during time marching from uniform initial conditions. Numerical solutions of several flow test cases are given and the computed results are compared with experimental measurements or theoretical solutions for validation of the method. It is shown that this method is reasonably efficient and robust.

AB - It is well known that the incorporation of two-equation turbulence models into explicit time-marching schemes, especially the low Reynolds versions, encounters various numerical difficulties. This paper sets out to review the numerical techniques presently used and to analyse the problems encountered. A numerical approach is proposed for the implementation of the low-Reynolds-number (LRN) k-∈ models into an explicit Runge-Kutta time-marching solver. In this method, a semi-implicit treatment of the source terms of the k and ∈ equations is introduced, and modified residuals of the two equations are used in order to adopt implicit residual smoothing. These techniques ensure the positivity of k and ∈ in the computation, even with uniform initial conditions and local time stepping. The asymptotic behaviour of the turbulence quantities during time marching is analysed, providing insight into the stability and convergence characteristics of the k and ∈ equations using explicit and implicit schemes with and without the damping terms. It is found that the time step size has a strong influence on the evolution of the turbulence quantities and the stability of the solution, whether implicit or explicit schemes are used. The time step size determined by the stability analysis is a function of k, ∈ and their damping terms. This analysis also provides a guide for the selection of time step size for the turbulence equations. It is found that the selection of time step size is an important factor in ensuring realistic values of k, ∈ and μt in the flow field during time marching from uniform initial conditions. Numerical solutions of several flow test cases are given and the computed results are compared with experimental measurements or theoretical solutions for validation of the method. It is shown that this method is reasonably efficient and robust.

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

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

M3 - Article

AN - SCOPUS:0031279415

VL - 28

SP - 487

EP - 499

JO - Advances in Engineering Software

JF - Advances in Engineering Software

SN - 0965-9978

IS - 8

ER -