A Hybrid Large Eddy Simulation Algorithm Based on the Implicit Domain Decomposition

The resolution of small near-wall eddies encountered in high-Reynolds number flows using large eddy simulation (LES) requires very fine meshes that may be computationally prohibitive. As a result, the use of wall-modeled LES as an alternative is becoming more popular. In this paper, the near-wall do...

Full description

Bibliographic Details
Main Authors: Amir E. Fard, Sergey Utyuzhnikov
Format: Article
Language:English
Published: MDPI AG 2023-10-01
Series:Mathematics
Subjects:
Online Access:https://www.mdpi.com/2227-7390/11/20/4340
_version_ 1797573120709623808
author Amir E. Fard
Sergey Utyuzhnikov
author_facet Amir E. Fard
Sergey Utyuzhnikov
author_sort Amir E. Fard
collection DOAJ
description The resolution of small near-wall eddies encountered in high-Reynolds number flows using large eddy simulation (LES) requires very fine meshes that may be computationally prohibitive. As a result, the use of wall-modeled LES as an alternative is becoming more popular. In this paper, the near-wall domain decomposition (NDD) approach that was originally developed for Reynolds-averaged Navier–Stokes simulations (RANSs) is extended to the hybrid RANS/LES zonal decomposition. The algorithm is implemented in two stages. First, the solution is computed everywhere with LES on a coarse grid using a new non-local slip boundary condition for the instantaneous velocity at the wall. The solution is then recomputed in the near-wall region with RANS. The slip boundary conditions used in the first stage guarantee that the composite solution is smooth at the inner/outer region interface. Another advantage of the model is that the turbulent viscosity in the inner region is computed based on the corresponding RANS velocity. This shows improvement over those hybrid models that have only one velocity field in the whole domain obtained from LES. The model is realized in the open source code OpenFOAM with different approximations of turbulent viscosity and is applied to the planar channel flow at frictional Reynolds numbers of <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mi>Re</mi><mi>τ</mi></msub><mo>=</mo><mn>950</mn></mrow></semantics></math></inline-formula>, 2000, and 4200. Mean streamwise velocity and Reynolds stress intensities are predicted reasonably well in comparison to the solutions obtained with unresolved LES and available DNS benchmarks. No additional forcing at the interface is required, while the log–layer mismatch is essentially reduced in all cases.
first_indexed 2024-03-10T21:05:11Z
format Article
id doaj.art-dc30c487d25144ea9fb8e73208a29df8
institution Directory Open Access Journal
issn 2227-7390
language English
last_indexed 2024-03-10T21:05:11Z
publishDate 2023-10-01
publisher MDPI AG
record_format Article
series Mathematics
spelling doaj.art-dc30c487d25144ea9fb8e73208a29df82023-11-19T17:14:38ZengMDPI AGMathematics2227-73902023-10-011120434010.3390/math11204340A Hybrid Large Eddy Simulation Algorithm Based on the Implicit Domain DecompositionAmir E. Fard0Sergey Utyuzhnikov1School of Engineering, University of Manchester, Manchester M13 9PL, UKSchool of Engineering, University of Manchester, Manchester M13 9PL, UKThe resolution of small near-wall eddies encountered in high-Reynolds number flows using large eddy simulation (LES) requires very fine meshes that may be computationally prohibitive. As a result, the use of wall-modeled LES as an alternative is becoming more popular. In this paper, the near-wall domain decomposition (NDD) approach that was originally developed for Reynolds-averaged Navier–Stokes simulations (RANSs) is extended to the hybrid RANS/LES zonal decomposition. The algorithm is implemented in two stages. First, the solution is computed everywhere with LES on a coarse grid using a new non-local slip boundary condition for the instantaneous velocity at the wall. The solution is then recomputed in the near-wall region with RANS. The slip boundary conditions used in the first stage guarantee that the composite solution is smooth at the inner/outer region interface. Another advantage of the model is that the turbulent viscosity in the inner region is computed based on the corresponding RANS velocity. This shows improvement over those hybrid models that have only one velocity field in the whole domain obtained from LES. The model is realized in the open source code OpenFOAM with different approximations of turbulent viscosity and is applied to the planar channel flow at frictional Reynolds numbers of <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mi>Re</mi><mi>τ</mi></msub><mo>=</mo><mn>950</mn></mrow></semantics></math></inline-formula>, 2000, and 4200. Mean streamwise velocity and Reynolds stress intensities are predicted reasonably well in comparison to the solutions obtained with unresolved LES and available DNS benchmarks. No additional forcing at the interface is required, while the log–layer mismatch is essentially reduced in all cases.https://www.mdpi.com/2227-7390/11/20/4340hybrid RANS/LESturbulencechannel flowlog–layer mismatchimplicit near-wall domain decomposition
spellingShingle Amir E. Fard
Sergey Utyuzhnikov
A Hybrid Large Eddy Simulation Algorithm Based on the Implicit Domain Decomposition
Mathematics
hybrid RANS/LES
turbulence
channel flow
log–layer mismatch
implicit near-wall domain decomposition
title A Hybrid Large Eddy Simulation Algorithm Based on the Implicit Domain Decomposition
title_full A Hybrid Large Eddy Simulation Algorithm Based on the Implicit Domain Decomposition
title_fullStr A Hybrid Large Eddy Simulation Algorithm Based on the Implicit Domain Decomposition
title_full_unstemmed A Hybrid Large Eddy Simulation Algorithm Based on the Implicit Domain Decomposition
title_short A Hybrid Large Eddy Simulation Algorithm Based on the Implicit Domain Decomposition
title_sort hybrid large eddy simulation algorithm based on the implicit domain decomposition
topic hybrid RANS/LES
turbulence
channel flow
log–layer mismatch
implicit near-wall domain decomposition
url https://www.mdpi.com/2227-7390/11/20/4340
work_keys_str_mv AT amirefard ahybridlargeeddysimulationalgorithmbasedontheimplicitdomaindecomposition
AT sergeyutyuzhnikov ahybridlargeeddysimulationalgorithmbasedontheimplicitdomaindecomposition
AT amirefard hybridlargeeddysimulationalgorithmbasedontheimplicitdomaindecomposition
AT sergeyutyuzhnikov hybridlargeeddysimulationalgorithmbasedontheimplicitdomaindecomposition