Large-eddy simulation of turbulent channel flows with conservative IDO scheme

Naoyuki Onodera, Takayuki Aoki, Hiromichi Kobayashi

Research output: Contribution to journalArticle

7 Citations (Scopus)

Abstract

The resolution of a numerical scheme in both physical and Fourier spaces is one of the most important requirements to calculate turbulent flows. A conservative form of the interpolated differential operator (IDO-CF) scheme is a multi-moment Eulerian scheme in which point values and integrated average values are separately defined in one cell. Since the IDO-CF scheme using high-order interpolation functions is constructed with compact stencils, the boundary conditions are able to be treated as easy as the 2nd-order finite difference method (FDM). It is unique that the first-order spatial derivative of the point value is derived from the interpolation function with 4th-order accuracy and the volume averaged value is based on the exact finite volume formulation, so that the IDO-CF scheme has higher spectral resolution than conventional FDMs with 4th-order accuracy. The computational cost to calculate the first-order spatial derivative with non-uniform grid spacing is one-third of the 4th-order FDM. For a large-eddy simulation (LES), we use the coherent structure model (CSM) in which the model coefficient is locally obtained from a turbulent structure extracted from a second invariant of the velocity gradient tensor, and the model coefficient correctly satisfies asymptotic behaviors to walls.The results of the IDO-CF scheme with the CSM for turbulent channel flows are compared to the FDM with the CSM and dynamic Smagorinsky model as well as the direct numerical simulation (DNS) by Moser et al. Adding the sub-grid scale stress tensor of LES to the IDO-CF scheme improves the profile of the mean velocity in comparison with an implicit eddy viscosity of the IDO-CF upwind scheme. The IDO-CF scheme with the CSM gives better turbulent intensities than conventional FDMs with the same number of grid points. The turbulent statistics calculated by IDO-CF scheme are in good agreement with the DNS at the various values of Reynolds number Reτ=180,395, and 590. It is found that the IDO-CF scheme is suitable for the turbulent flow computation and improves the turbulent statistics with compact stencils.

Original languageEnglish
Pages (from-to)5787-5805
Number of pages19
JournalJournal of Computational Physics
Volume230
Issue number14
DOIs
Publication statusPublished - 2011 Jun 20

Fingerprint

Large eddy simulation
large eddy simulation
channel flow
Channel flow
Model structures
Finite difference method
Frequency division multiplexing
Direct numerical simulation
Turbulent flow
Tensors
Interpolation
grids
Statistics
Derivatives
direct numerical simulation
turbulent flow
interpolation
Spectral resolution
statistics
Dynamic models

Keywords

  • Coherent structure model
  • IDO scheme
  • Large-eddy simulation
  • Multi-moment scheme

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy (miscellaneous)

Cite this

Large-eddy simulation of turbulent channel flows with conservative IDO scheme. / Onodera, Naoyuki; Aoki, Takayuki; Kobayashi, Hiromichi.

In: Journal of Computational Physics, Vol. 230, No. 14, 20.06.2011, p. 5787-5805.

Research output: Contribution to journalArticle

@article{676593d5a688453a8ef7e9104c56b7a1,
title = "Large-eddy simulation of turbulent channel flows with conservative IDO scheme",
abstract = "The resolution of a numerical scheme in both physical and Fourier spaces is one of the most important requirements to calculate turbulent flows. A conservative form of the interpolated differential operator (IDO-CF) scheme is a multi-moment Eulerian scheme in which point values and integrated average values are separately defined in one cell. Since the IDO-CF scheme using high-order interpolation functions is constructed with compact stencils, the boundary conditions are able to be treated as easy as the 2nd-order finite difference method (FDM). It is unique that the first-order spatial derivative of the point value is derived from the interpolation function with 4th-order accuracy and the volume averaged value is based on the exact finite volume formulation, so that the IDO-CF scheme has higher spectral resolution than conventional FDMs with 4th-order accuracy. The computational cost to calculate the first-order spatial derivative with non-uniform grid spacing is one-third of the 4th-order FDM. For a large-eddy simulation (LES), we use the coherent structure model (CSM) in which the model coefficient is locally obtained from a turbulent structure extracted from a second invariant of the velocity gradient tensor, and the model coefficient correctly satisfies asymptotic behaviors to walls.The results of the IDO-CF scheme with the CSM for turbulent channel flows are compared to the FDM with the CSM and dynamic Smagorinsky model as well as the direct numerical simulation (DNS) by Moser et al. Adding the sub-grid scale stress tensor of LES to the IDO-CF scheme improves the profile of the mean velocity in comparison with an implicit eddy viscosity of the IDO-CF upwind scheme. The IDO-CF scheme with the CSM gives better turbulent intensities than conventional FDMs with the same number of grid points. The turbulent statistics calculated by IDO-CF scheme are in good agreement with the DNS at the various values of Reynolds number Reτ=180,395, and 590. It is found that the IDO-CF scheme is suitable for the turbulent flow computation and improves the turbulent statistics with compact stencils.",
keywords = "Coherent structure model, IDO scheme, Large-eddy simulation, Multi-moment scheme",
author = "Naoyuki Onodera and Takayuki Aoki and Hiromichi Kobayashi",
year = "2011",
month = "6",
day = "20",
doi = "10.1016/j.jcp.2011.04.004",
language = "English",
volume = "230",
pages = "5787--5805",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "14",

}

TY - JOUR

T1 - Large-eddy simulation of turbulent channel flows with conservative IDO scheme

AU - Onodera, Naoyuki

AU - Aoki, Takayuki

AU - Kobayashi, Hiromichi

PY - 2011/6/20

Y1 - 2011/6/20

N2 - The resolution of a numerical scheme in both physical and Fourier spaces is one of the most important requirements to calculate turbulent flows. A conservative form of the interpolated differential operator (IDO-CF) scheme is a multi-moment Eulerian scheme in which point values and integrated average values are separately defined in one cell. Since the IDO-CF scheme using high-order interpolation functions is constructed with compact stencils, the boundary conditions are able to be treated as easy as the 2nd-order finite difference method (FDM). It is unique that the first-order spatial derivative of the point value is derived from the interpolation function with 4th-order accuracy and the volume averaged value is based on the exact finite volume formulation, so that the IDO-CF scheme has higher spectral resolution than conventional FDMs with 4th-order accuracy. The computational cost to calculate the first-order spatial derivative with non-uniform grid spacing is one-third of the 4th-order FDM. For a large-eddy simulation (LES), we use the coherent structure model (CSM) in which the model coefficient is locally obtained from a turbulent structure extracted from a second invariant of the velocity gradient tensor, and the model coefficient correctly satisfies asymptotic behaviors to walls.The results of the IDO-CF scheme with the CSM for turbulent channel flows are compared to the FDM with the CSM and dynamic Smagorinsky model as well as the direct numerical simulation (DNS) by Moser et al. Adding the sub-grid scale stress tensor of LES to the IDO-CF scheme improves the profile of the mean velocity in comparison with an implicit eddy viscosity of the IDO-CF upwind scheme. The IDO-CF scheme with the CSM gives better turbulent intensities than conventional FDMs with the same number of grid points. The turbulent statistics calculated by IDO-CF scheme are in good agreement with the DNS at the various values of Reynolds number Reτ=180,395, and 590. It is found that the IDO-CF scheme is suitable for the turbulent flow computation and improves the turbulent statistics with compact stencils.

AB - The resolution of a numerical scheme in both physical and Fourier spaces is one of the most important requirements to calculate turbulent flows. A conservative form of the interpolated differential operator (IDO-CF) scheme is a multi-moment Eulerian scheme in which point values and integrated average values are separately defined in one cell. Since the IDO-CF scheme using high-order interpolation functions is constructed with compact stencils, the boundary conditions are able to be treated as easy as the 2nd-order finite difference method (FDM). It is unique that the first-order spatial derivative of the point value is derived from the interpolation function with 4th-order accuracy and the volume averaged value is based on the exact finite volume formulation, so that the IDO-CF scheme has higher spectral resolution than conventional FDMs with 4th-order accuracy. The computational cost to calculate the first-order spatial derivative with non-uniform grid spacing is one-third of the 4th-order FDM. For a large-eddy simulation (LES), we use the coherent structure model (CSM) in which the model coefficient is locally obtained from a turbulent structure extracted from a second invariant of the velocity gradient tensor, and the model coefficient correctly satisfies asymptotic behaviors to walls.The results of the IDO-CF scheme with the CSM for turbulent channel flows are compared to the FDM with the CSM and dynamic Smagorinsky model as well as the direct numerical simulation (DNS) by Moser et al. Adding the sub-grid scale stress tensor of LES to the IDO-CF scheme improves the profile of the mean velocity in comparison with an implicit eddy viscosity of the IDO-CF upwind scheme. The IDO-CF scheme with the CSM gives better turbulent intensities than conventional FDMs with the same number of grid points. The turbulent statistics calculated by IDO-CF scheme are in good agreement with the DNS at the various values of Reynolds number Reτ=180,395, and 590. It is found that the IDO-CF scheme is suitable for the turbulent flow computation and improves the turbulent statistics with compact stencils.

KW - Coherent structure model

KW - IDO scheme

KW - Large-eddy simulation

KW - Multi-moment scheme

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

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

U2 - 10.1016/j.jcp.2011.04.004

DO - 10.1016/j.jcp.2011.04.004

M3 - Article

VL - 230

SP - 5787

EP - 5805

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 14

ER -