Main CfdTm Forum Saturne Aster ATAAC TWikiEditRaw editAttachPrint version

TWiki> Main Web>TWikiUsers>NeilAshton>SSTSAS (2011-06-19, NeilAshton)

TWiki> Main Web>TWikiUsers>NeilAshton>SSTSAS (2011-06-19, NeilAshton)

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

See Menter F. R.,Egorov Y. (2005) A Scale-adaptive simulation model using two-equation models. AIAA paper 2005-1095.

The coefficients are:

Where is the modelled turbulence length scale

The von Karman Length Scale is as follows:

The second velocity derivative is:

Limiter for the :

- include and define the index, ISAS, which turns off the SAS source term if it is equals to zero and enables it when it is one.
- include and define and .

- Extra memory is required for the calculations of the magnitude of the gradients of and . The call headers need to be changed to include extra working arrays, W10 and W11.

- W10 and W11 defined as double precision.

The SAS model will run using an index argument, where if ISAS.EQ.1 then the SAS source term will be used and if ISAS.EQ.0 then the standard SST will be run.

- Working array W10 is

- Working array W11 is

- Working array W11 is max (W10/),W11/)

The above operation frees up W10 working array for the following calculations

- Strain invariant is

- Calculate the laplacian of U in part 7 of turkw.F

- include IF statement for ISAS

- Take VISCCFA and ITRGRP to calculate the laplacian, W2

- Repeat for U, V and W

- The work arrays for ITRGRP are W3,W4,W5,W6,W7,W8

- U, W9 = -divgrad U
- V, W9 = -divgrad V
- W, W9 = divgrad W

- W10 is used for the square root of the sum of the gradients.

**Source terms**

- Calculate $L_{vk}^{min}$.

- Calculate $L_{vk}$, and then setup max statement

- Calculate $L$

- Calculate $Q_{SAS}$

- Divide working array W10 by the volume

- Add an if statement to the source terms, where if ISAS.Eq.1 the source term has QSAS added to it

- Change memkom call to include extra working arrays W10 and W11

- Change call to turbkw, to include extra working arrays W10 and W11

- Add extra working arrays, IW10, IW11

- Change definitions, integers for IW10, IW11

- Change memory pointers to include new working arrays

- set default value, ISAS = 0

- include incsas.h

- include default values for and

- ISAS = 0
- = 0.11
- = 3.51
- include incsas.h

- include incsas.h
- Include printout for ISAS
- include printout for

- include incsas.h
- include an IF statement to check that DDES and SAS are not both running.

I | Attachment | Action | Size | Date | Who | Comment |
---|---|---|---|---|---|---|

csas_sst_32.pdf | manage | 151.2 K | 2011-06-19 - 22:17 | NeilAshton | csas calibration for the SST-SAS model |

Edit | Attach | Print version | History: r1 | Backlinks | Raw View | Raw edit | More topic actions

Topic revision: r1 - 2011-06-19 - 22:17:31 - NeilAshton

Code_Saturne

KNOO

FLOMANIA

Photo Wall

Upcoming Events

Add Event

Copyright © by the contributing authors. Unless noted otherwise, all material on this web site is the property of the contributing authors.