# Implementation of EBRSM - DDES model in Code_Saturne

This page contains the equations necessary for the implementation of the EBRSM - DDES model into Code Saturne.

Spalart, P.R. et al. A New Version of Detached-eddy Simulation, Resistant to Ambiguous Grid Densities. Theoretical and Computational Fluid Dynamics 20, 181-195(2006).

Manceau R et al. Elliptic Blending model: A new near-wall Reynolds-stress turbulence closure. Physics of Fluids, 14, 744-754(2001)

## Source terms for equation

The dissipation term of the equations becomes:

The length scale used in a DDES formulation is a switch between the RANS length scale and a LES length scale formulation:

The LES Length scale is:

The RANS length scale is defined as:

Where is:

The function is:

Where, and is the distance to the wall

## Modified Code_Saturne v1.4 subroutines

### incddes.h

• include and define the index values , IDDES which turns off the change in dissipation transport equation if it is equals to zero and enables it when it is one.
• include and define

### usini1.F

• IDDES = 0
• = 0.65
• include incddes.h

### iniini.F

• set default value, IDDES= 0,

• include incddes.h

• include default values for

### resebm.F

• allocate new memory for gradient calculations, W10 and W11

• Define W10,W11 as double precision

• Working array W10 is

• Working array W11 is

Source terms

• Calculate .

• Calculate

• Calculate

• Calculate

• Add an if statement to the equation, where if IDDES.Eq.1 the change in will be switched on

### turrij.F

• Include new working arrays, W10, W11 in the call to resebm
• Define W10 and W11 as double precision

### modini.F

• Change INEEDY IF statement to include ITURB = 32 to calculate the distance to the wall

### impini.F

• include incddes.h
• Include printout for IDDES
• include printout for

Current Tags:
create new tag
, view all tags
Topic revision: r6 - 2010-10-27 - 14:00:26 - NeilAshton
Main Web
21 Mar 2018

Manchester CfdTm
Code_Saturne

ATAAC
KNOO

DESider
FLOMANIA