Assessment of very high order of accuracy in implicit LES models

Andrew Mosedale, Dimitris Drikakis

Research output: Contribution to journalArticle

31 Citations (Scopus)

Abstract

This paper looks at the use of high-resolution and very high-order methods for implicit large-eddy simulation (ILES), with the specific example of simulating the multicomponent two-dimensional single-mode Richtmyer–Meshkov instability for which experimental data is available. The two gases are air and SF6, making stringent demands on the models used in the code. The interface between the two gases is initialized with a simple sinusoidal perturbation over a wavelength of 59mm, and a shock of strength Mach 1.3 is passed through this interface. The main comparison is between the second-order monotone upwind-centered scheme for conservation law methods of van Leer (1979, “Towards the Ultimate Conservative Difference Scheme ,” J. Comput. Phys.32, pp. 101–136) and the current state-of-the-art weighted essentially nonoscillatory interpolation, which is presented to ninth order, concentrating on the effect on resolution of the instability on coarse grids. The higher-order methods as expected provide better resolved and more physical features than the second-order methods on the same grid resolution. While it is not possible to make a definitive statement, the simulations have indicated that the extra time required for the higher-order reconstruction is less than the time saved by being able to obtain the same or better accuracy at lower computational cost (fewer grid points). It should also be noted that all simulations give a good representation of the growth rate of the instability, comparing very favorably to the experimental results, and as such far better than the currently existing theoretical models. This serves to further indicate that the ILES approach is capable of providing accurately physical information despite the lack of any formal subgrid model.
Original languageEnglish
Pages (from-to)1497-1503
Number of pages7
JournalJournal of Fluids Engineering
Volume129
Issue number12
DOIs
Publication statusPublished - 1 Dec 2007

Fingerprint

Large eddy simulation
Gases
Mach number
Conservation
Interpolation
Wavelength
Air
Costs

Keywords

  • computer simulation
  • gas dynamics
  • image reconstruction
  • information analysis
  • system stability
  • conservation law
  • experimental data
  • Richtmyer–Meshkov instability
  • second-order methods
  • large eddy simulation

Cite this

Mosedale, Andrew ; Drikakis, Dimitris. / Assessment of very high order of accuracy in implicit LES models. In: Journal of Fluids Engineering. 2007 ; Vol. 129, No. 12. pp. 1497-1503.
@article{2b0ee996569449158def12f8d869167a,
title = "Assessment of very high order of accuracy in implicit LES models",
abstract = "This paper looks at the use of high-resolution and very high-order methods for implicit large-eddy simulation (ILES), with the specific example of simulating the multicomponent two-dimensional single-mode Richtmyer–Meshkov instability for which experimental data is available. The two gases are air and SF6, making stringent demands on the models used in the code. The interface between the two gases is initialized with a simple sinusoidal perturbation over a wavelength of 59mm, and a shock of strength Mach 1.3 is passed through this interface. The main comparison is between the second-order monotone upwind-centered scheme for conservation law methods of van Leer (1979, “Towards the Ultimate Conservative Difference Scheme ,” J. Comput. Phys.32, pp. 101–136) and the current state-of-the-art weighted essentially nonoscillatory interpolation, which is presented to ninth order, concentrating on the effect on resolution of the instability on coarse grids. The higher-order methods as expected provide better resolved and more physical features than the second-order methods on the same grid resolution. While it is not possible to make a definitive statement, the simulations have indicated that the extra time required for the higher-order reconstruction is less than the time saved by being able to obtain the same or better accuracy at lower computational cost (fewer grid points). It should also be noted that all simulations give a good representation of the growth rate of the instability, comparing very favorably to the experimental results, and as such far better than the currently existing theoretical models. This serves to further indicate that the ILES approach is capable of providing accurately physical information despite the lack of any formal subgrid model.",
keywords = "computer simulation, gas dynamics, image reconstruction, information analysis, system stability, conservation law, experimental data, Richtmyer–Meshkov instability, second-order methods, large eddy simulation",
author = "Andrew Mosedale and Dimitris Drikakis",
year = "2007",
month = "12",
day = "1",
doi = "10.1115/1.2801374",
language = "English",
volume = "129",
pages = "1497--1503",
journal = "Journal of Fluids Engineering",
issn = "0098-2202",
number = "12",

}

Assessment of very high order of accuracy in implicit LES models. / Mosedale, Andrew; Drikakis, Dimitris.

In: Journal of Fluids Engineering, Vol. 129, No. 12, 01.12.2007, p. 1497-1503.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Assessment of very high order of accuracy in implicit LES models

AU - Mosedale, Andrew

AU - Drikakis, Dimitris

PY - 2007/12/1

Y1 - 2007/12/1

N2 - This paper looks at the use of high-resolution and very high-order methods for implicit large-eddy simulation (ILES), with the specific example of simulating the multicomponent two-dimensional single-mode Richtmyer–Meshkov instability for which experimental data is available. The two gases are air and SF6, making stringent demands on the models used in the code. The interface between the two gases is initialized with a simple sinusoidal perturbation over a wavelength of 59mm, and a shock of strength Mach 1.3 is passed through this interface. The main comparison is between the second-order monotone upwind-centered scheme for conservation law methods of van Leer (1979, “Towards the Ultimate Conservative Difference Scheme ,” J. Comput. Phys.32, pp. 101–136) and the current state-of-the-art weighted essentially nonoscillatory interpolation, which is presented to ninth order, concentrating on the effect on resolution of the instability on coarse grids. The higher-order methods as expected provide better resolved and more physical features than the second-order methods on the same grid resolution. While it is not possible to make a definitive statement, the simulations have indicated that the extra time required for the higher-order reconstruction is less than the time saved by being able to obtain the same or better accuracy at lower computational cost (fewer grid points). It should also be noted that all simulations give a good representation of the growth rate of the instability, comparing very favorably to the experimental results, and as such far better than the currently existing theoretical models. This serves to further indicate that the ILES approach is capable of providing accurately physical information despite the lack of any formal subgrid model.

AB - This paper looks at the use of high-resolution and very high-order methods for implicit large-eddy simulation (ILES), with the specific example of simulating the multicomponent two-dimensional single-mode Richtmyer–Meshkov instability for which experimental data is available. The two gases are air and SF6, making stringent demands on the models used in the code. The interface between the two gases is initialized with a simple sinusoidal perturbation over a wavelength of 59mm, and a shock of strength Mach 1.3 is passed through this interface. The main comparison is between the second-order monotone upwind-centered scheme for conservation law methods of van Leer (1979, “Towards the Ultimate Conservative Difference Scheme ,” J. Comput. Phys.32, pp. 101–136) and the current state-of-the-art weighted essentially nonoscillatory interpolation, which is presented to ninth order, concentrating on the effect on resolution of the instability on coarse grids. The higher-order methods as expected provide better resolved and more physical features than the second-order methods on the same grid resolution. While it is not possible to make a definitive statement, the simulations have indicated that the extra time required for the higher-order reconstruction is less than the time saved by being able to obtain the same or better accuracy at lower computational cost (fewer grid points). It should also be noted that all simulations give a good representation of the growth rate of the instability, comparing very favorably to the experimental results, and as such far better than the currently existing theoretical models. This serves to further indicate that the ILES approach is capable of providing accurately physical information despite the lack of any formal subgrid model.

KW - computer simulation

KW - gas dynamics

KW - image reconstruction

KW - information analysis

KW - system stability

KW - conservation law

KW - experimental data

KW - Richtmyer–Meshkov instability

KW - second-order methods

KW - large eddy simulation

UR - http://www.scopus.com/inward/record.url?eid=2-s2.0-38149074129&partnerID=40&md5=b5ef12ab62f4d5c775acf510c72acf19

U2 - 10.1115/1.2801374

DO - 10.1115/1.2801374

M3 - Article

VL - 129

SP - 1497

EP - 1503

JO - Journal of Fluids Engineering

JF - Journal of Fluids Engineering

SN - 0098-2202

IS - 12

ER -