Loading…
Hybrid LES/URANS simulation of turbulent natural convection by BEM
In this paper we have developed a hybrid LES/URANS turbulent model for a BEM based turbulent fluid flow solver. We employed the unified LES/URANS approach, where the interface between the LES and URANS regions is defined using a physical quantity, which dynamically changes during numerical simulatio...
Saved in:
Published in: | Engineering analysis with boundary elements 2015-12, Vol.61, p.16-26 |
---|---|
Main Authors: | , , |
Format: | Article |
Language: | English |
Subjects: | |
Citations: | Items that this one cites Items that cite this one |
Online Access: | Get full text |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
Summary: | In this paper we have developed a hybrid LES/URANS turbulent model for a BEM based turbulent fluid flow solver. We employed the unified LES/URANS approach, where the interface between the LES and URANS regions is defined using a physical quantity, which dynamically changes during numerical simulation. The main characteristic of the unified hybrid model is that only one set of governing equations is used for fluid flow simulation in both the LES and URANS regions. Regions where turbulent kinetic energy is calculated by LES and URANS models are determined using a switching criterion. We used the Reynolds number based on turbulent kinetic energy and the Reynolds number based on total turbulent kinetic energy to establish the LES/URANS interface switching criterion. Depending on flow characteristics and with the use of switching criterion, we chose between sub-grid scale viscosity (SGS) and URANS effective viscosity. The SGS or URANS effective viscosity is used in the transport equation for turbulent kinetic energy and in governing equations for fluid flow. The developed numerical algorithm was tested by simulating turbulent natural convection within a square cavity. The hybrid turbulent model was implemented within a numerical algorithm based on the boundary element method, where single domain and sub-domain approaches are used. The governing equations are written in velocity–vorticity formulation. We used the false transient time scheme for the kinematics equation. |
---|---|
ISSN: | 0955-7997 1873-197X |
DOI: | 10.1016/j.enganabound.2015.06.005 |