
[Sponsors] 
grid length scale in DES SpalartAllmaras  definition of dTilda 

LinkBack  Thread Tools  Display Modes 
February 2, 2011, 17:14 
grid length scale in DES SpalartAllmaras  definition of dTilda

#1 
New Member
Join Date: Mar 2009
Posts: 18
Rep Power: 9 
Does someone know how is calculate the grid length scale used in dTilda? In the original model, Spalart proposed to use the maximum length of the cell, but in the lesdelta options I was unable to find something similar. I am using the cube root delta.


May 16, 2012, 18:50 

#2  
New Member
Join Date: Nov 2011
Posts: 13
Rep Power: 7 
Quote:
One of the other options is maxDeltaxyz: 'Maximum of x, y and z; for structured hex cells only'. Is this equivalent to the DES Delta parameter, and if so does it work with snappyHexMesh grids? 

May 17, 2012, 04:53 

#3 
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 364
Rep Power: 12 
Until OF2.1.0, the SpalartAllmaras DES delta definition (max DeltaX, DeltaY, DeltaZ) was not available, so the only reasonable choice appeared to be the classical cubic root Delta definition (which, however, does not return the expected behavior in the RANS zone of the original SA DES model). Now it seems that someone has detected (at last...) this deficiency, so the maxDeltaxyz is the consistent choice for the DES model in its original form (although "consistent" does not mean "best", especially in the LES region of the flow with highly anisotropic or unstructured grids)
Regards V. 

May 17, 2012, 06:35 

#4 
New Member
Join Date: Nov 2011
Posts: 13
Rep Power: 7 
Thanks for your advice, I'll try both with DDES and see what the differences are.


July 17, 2012, 05:07 

#5 
Senior Member
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 168
Rep Power: 8 
Dear all,
I have a question concerning the computation of (max DeltaX, DeltaY, DeltaZ) in src/turbulenceModels/LES/LESdeltas/maxDeltaxyz. Here are some lines of the code HTML Code:
forAll(cells,cellI) { scalar deltaMaxTmp = 0.0; const labelList& cFaces = mesh().cells()[cellI]; const point& centrevector = mesh().cellCentres()[cellI]; forAll(cFaces, cFaceI) { label faceI = cFaces[cFaceI]; const point& facevector = mesh().faceCentres()[faceI]; scalar tmp = mag(facevector  centrevector); if (tmp > deltaMaxTmp) { deltaMaxTmp = tmp; } } hmax()[cellI] = deltaCoeff_*deltaMaxTmp; } But in order to obtain the desired (max DeltaX, DeltaY, DeltaZ), isn't it necessary to multiply the obtained value with a factor 2? Is it necessary to set deltaCoeff to 2 in constant/LESProperties? As far as I can see the default value is 1. I would appreciate any answers. Kind Regards Anne 

July 24, 2012, 05:07 

#6 
Member
David GISEN
Join Date: Jul 2009
Location: Germany
Posts: 51
Rep Power: 9 
Hi Anne,
if you are using Hexamesh, you have to set deltaCoeff to 2 to obtain the correct value. As far as I understand, it is implemented in the present way to account for Polymeshes. There is a bug report (ID 0000474) at http://www.openfoam.org/bugs/ including this topic. Regards, David 

July 25, 2012, 10:23 

#7 
Senior Member
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 168
Rep Power: 8 
Hey David,
thank you for this hint. I found this thread a while ago, but could not find it anymore. Kind Regards Anne 

August 1, 2012, 05:30 

#8  
Senior Member
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 168
Rep Power: 8 
Hello,
I have one small further question. Citation from the Bug report site Quote:
Quote:
Kind Regards Anne 

May 27, 2013, 16:28 

#9 
Member
Join Date: Jul 2010
Posts: 39
Rep Power: 8 
Hi everyone;
Does anyone know the answer to Ann's previous question about how to define deltaCoeff for hex meshes? I believe that when using DES or DDES: maxDeltaxyzCoeffs { deltaCoeff 2; } but for IDDES the deltaCoeff should be defined as 1. So: IDDESDeltaCoeffs { deltaCoeff 1; cw 0.15; } Is that correct? Thanks for your help. Last edited by MaryBau; May 29, 2013 at 11:48. 

December 11, 2014, 05:48 

#10 
New Member
Join Date: Nov 2011
Posts: 28
Rep Power: 7 
This question is actual for me. I use LES turbulence model and snappyhexmesh. What is the delta should I choose 1 or 2? Because there is the problem of the flow past cylinder. If Deltacoeffs = 1 then Cd = 0.33, if deltaCoeffs = 2 then Cd = 1.2 as in the experiment. But I have coarse mesh y+ ~ 40, V=72m\s
Thanks. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
How to use DES well?  Daniel  Main CFD Forum  0  October 26, 2008 06:59 
grid generation for Detached eddy simulation (DES)  ROOZBEH  Main CFD Forum  1  June 18, 2007 12:02 
Kolomogorov length scale  Bobby  Main CFD Forum  8  November 23, 2002 04:54 
Structured Grid Definition  craig shores  Main CFD Forum  3  March 21, 2001 15:48 
Grid Independent Solution  Chuck Leakeas  Main CFD Forum  2  May 26, 2000 11:18 