|
[Sponsors] |
How to set initial conditions for LES 3D case? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 20, 2018, 21:29 |
How to set initial conditions for LES 3D case?
|
#1 |
New Member
Derek
Join Date: Mar 2018
Posts: 10
Rep Power: 8 |
Still doing unsteady flow past 3D cylinder. Can someone advise how to set initial conditions for LES 3D case?? Or if you could provide any similar tutorial. Much appreciate!
Initial conditions includes: k, nut, nuTilda, p and u (btw, is "s" necessary for LES? |
|
March 21, 2018, 01:12 |
|
#2 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
I use:
k: Fixed value 0.1 nearly everywhere, kqRWallFunction at a surface of interest nut: calculated nearly eveyrwhere, nutkWallFunction at a surface of interest p: zeroGradient nearly everywhere, a fixed value at one of the patches I don't need much more for the dynamicKEqn model.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
March 21, 2018, 23:07 |
|
#3 |
New Member
Derek
Join Date: Mar 2018
Posts: 10
Rep Power: 8 |
Thank you Sir!
|
|
March 30, 2018, 11:38 |
|
#4 |
New Member
Derek
Join Date: Mar 2018
Posts: 10
Rep Power: 8 |
sir, will the initial k and nut significant affect the calculation results?
|
|
March 30, 2018, 12:25 |
|
#5 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
In any case it influences how fast and stable the simulation starts. But you may try this out for your case.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
April 22, 2018, 02:41 |
|
#6 | |
New Member
JD
Join Date: May 2017
Posts: 24
Rep Power: 9 |
Quote:
So, with these wall functions in place, do you use 30<y+<300 ? |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Problem with chtMultiregionFoam radiation boundary condition | baran_foam | OpenFOAM Running, Solving & CFD | 10 | December 17, 2019 17:36 |
Foam::error::printStack after large bounding epsilon events | StoneCFD727 | OpenFOAM Running, Solving & CFD | 1 | January 1, 2018 17:04 |
simpleFoam error - "Floating point exception" | mbcx4jc2 | OpenFOAM Running, Solving & CFD | 12 | August 4, 2015 02:20 |
should Courant number always be kept below 1? | wc34071209 | OpenFOAM Running, Solving & CFD | 16 | March 9, 2014 19:31 |
pimpleFoam: turbulence->correct(); is not executed when using residualControl | hfs | OpenFOAM Running, Solving & CFD | 3 | October 29, 2013 08:35 |