# Implementation of the - DDES model in Code_Saturne

This page contains the equations necessary for the implementation of the DDES model into Code Saturne v2.0.1.

See A hybrid numerical scheme for a new formulation of delayed detached-eddy simulation (DDES) based on elliptic relaxation, N. Ashton, A. Revell, R. Prosser, Journal of Physics: Conference Series, 13th European Turbulence Conference, Vol. 318, (2011)

See 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).

Laurence, D L, Uribe, J C & Utyuzhinkov, S V 2004 A robust formulation of the v2 − f model. Flow Turbulence and Combustion 73, 169–185.

## Source term for equation

Destruction term in the equation

 (1)

Where the function is:

 (2)

The length scales are:

 (3)

Where, , and is the distance to the wall

## Modified Code_Saturne v2.0.1 subroutines

### incddes.h

• include and define the index values , IDDES, which turns off the change in K transport equation if it is equals to zero and enables it when it is 1 or 2 depending on the formulation.
• include and define .

### turbke.F

• Extra memory is required for the storage of the wall distance and . The call headers need to be changed to include extra working arrays, W10 and W11.

• W10 and W11 defined as double precision.

The DDES model will run using an index argument, where if IDDES.EQ.2 then the DDES modification will be used.

• 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.2 then the becomes

### tridim.F

• Change memkep call to include extra working arrays W10 and W11

• Change call to turbke, to include extra working arrays W10 and W11

### memkep.F

• Add extra working arrays, IW10, IW11

• Change definitions, integers for IW10, IW11

• Change memory pointers to include new working arrays

### modini.F

• Change INEEDY IF statement to include the v2f model. i.e ITURB = 50 to calculate the distance to the wall

### iniini.F

• set default value, IDDES = 0

• include incddes.h

• include default values for

### usini1.F

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

### impini.F

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

### verini.F

• include incddes.h
• include an IF statement to check that SAS and DDES cannot run at the same time.

Current Tags:
create new tag
, view all tags
Topic attachments
I Attachment Action Size Date Who Comment
pdf cddes_phif_etc.pdf manage 151.2 K 2011-06-21 - 20:47 NeilAshton Cddes calibration for the phif DDES model
Topic revision: r4 - 2012-04-28 - 21:27:03 - NeilAshton
Main Web
18 Mar 2019

Manchester CfdTm
Code_Saturne

ATAAC
KNOO

DESider
FLOMANIA