CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

LES WallFunctions in O.F. 2.2.0 (nutUSpaldingWallFunction)

Register Blogs Community New Posts Updated Threads Search

Like Tree6Likes
  • 3 Post By Villo
  • 1 Post By RodriguezFatz
  • 2 Post By Villo

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 5, 2013, 07:25
Default LES WallFunctions in O.F. 2.2.0 (nutUSpaldingWallFunction)
  #1
New Member
 
Tommy V
Join Date: Nov 2013
Posts: 29
Rep Power: 12
Villo is on a distinguished road
Dear Foamers,
hi to all! new in the forum and, almost, in OpenFoam... so don`t complain too much
I`m simulating a transiet injection of a high density fluid into a geometry filled with a low density fluid, solving with twoLiquidMixingFoam with this turbulenceModel:
  1. LES model: Vreman
  2. LESdelta: cubeRootVol
I had to move from OF 2.0.0 to OF 2.2.0, otherwise the alpha field was over the upper limit and lower than the lower limit (also with "famous" numerical schemes that shows great div(alpha) limiting and stability as filteredLinear, vanLeer).
My mesh is fine but, maybe, not enough for the physics of the case. So i want to check differences in the solution using and not wall functions.
After this introduction let`s go to the questions:
  1. i found that nuSgsWallFunction is not still available in OF2.2.0, was it replaced from nutUSpaldingWallFunction?
  2. in OF 2.0, if wall functions are applied you have to set nuSgsWallFunction for 0/nuSgs and (if you want) nutSOMETHING for 0/nuTilda. How to treat it in OF 2.2.0, is nutUSpaldingWallFunction allowed for both viscosity definitions?
  3. trying to apply nutUSpaldingWallFunction i have follow this way: create the wallFunctionDict (Spalding law), run the pre-processing utilities wallFunctionTable, generate the uPlusWallFunctionData.gz archive. Then i have set in the nuSgs file (for wall patch) the following string

    wall
    {
    type nutUSpaldingWallFunction;
    value ???? ;
    }

    but which value ( ????) does i have to set? I image that i have to insert the uPlusWallFunctionData, how?
  4. for k field i think that i can continue to use kqRWallFunction, isn`t it?
Any suggestion?
Cheers
Tommy

p.s. i would ask about the new OF 2.2.0 settings for numerical schemes (http://www.openfoam.org/version2.2.0/numerics.php) to understand how does works a similar settings div(phi,h) bounded Gauss linearUpwind grad(h);... but maybe in a second time
Villo is offline   Reply With Quote

Old   December 9, 2013, 10:08
Default
  #2
New Member
 
Tommy V
Join Date: Nov 2013
Posts: 29
Rep Power: 12
Villo is on a distinguished road
No suggestions? also the comment "stupid questions cannot expect great answer" would be useful!
tariq, random_ran and Gang Wang like this.
Villo is offline   Reply With Quote

Old   December 10, 2013, 05:20
Default
  #3
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27
RodriguezFatz will become famous soon enough
to 3) -> As far as I know the "value" in such cases is just used for the very first iteration (of the first time-step) as initialization. After that it uses the calculated value of your boundary condition. At least I don't worry about that.
tariq likes this.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   December 10, 2013, 09:18
Default
  #4
New Member
 
Tommy V
Join Date: Nov 2013
Posts: 29
Rep Power: 12
Villo is on a distinguished road
Hi Philipp, thanks for your reply!
So after the first iteration the uPlusWallFunctionData.gz is automatically read?
I want to try to make you one more question: when this wall function is applied to 0/nuSgs how do you use to treat the boundary condition in 0/nuTilda?
Villo is offline   Reply With Quote

Old   December 10, 2013, 09:44
Default
  #5
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 27
RodriguezFatz will become famous soon enough
That's at least what I understood!

I am not 100% sure about the wall functions, because I don't catch on the whole code. But I think this is done by the "correctBoundaryConditions" command that comes after all the calculations in all the turbulence model files. Have a look at your LES model files! Maybe someone can confirm that?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   February 6, 2014, 04:31
Default
  #6
New Member
 
Tommy V
Join Date: Nov 2013
Posts: 29
Rep Power: 12
Villo is on a distinguished road
Hey there!
Playing around wallFunctions, in order to apply nutUSpaldingWallFunction is sufficient to set-up in 0/nuSgs:

bottom
{
type nutUSpaldingWallFunction;
value uniform 1e-4;
}


there is no need to use wallFunctionDict.
I have some doubts about the nuSgs because the solved field seems to be strictly related to the initial value (as showed in the nuSgs chart).
The effect of this wall function (about velocity) is showed in the attached file.
blue line: wall function on left, LES solving (wall refined) on right
red line: wall function on left and right

About wallFunctionDict: seems to be necessary when you set in 0/nuSgs (or 0/nut 0/nuTilda etc etc) for the patch wall
bottom
{
type nutTabulatedWallFunction;
uPlusTable uPlusWallFunctionData;
value uniform 1e-4;
}

(uPlusWallFunctionData have to be the same as setted-up in ./constant/wallFunctionDict)

Thene you have to choice on ./constant/wallFunctionDict the kind of tabulated wallFunction (as shown in the following output from compiling trial):
Valid tabulatedWallFunction types are:
2
(
SpaldingsLaw
general
)


Thene run on ./ wallFunctionTable and will be created a ./constant/uPlusWallFunctionData available for the computation. Just in this case the solver will read the generated table (i`m running a computation to confirm that)
Attached Images
File Type: png wall.PNG (50.8 KB, 217 views)
File Type: png nuSgs.PNG (22.8 KB, 181 views)
kovamaniac and dneves like this.

Last edited by Villo; February 6, 2014 at 06:24.
Villo is offline   Reply With Quote

Reply

Tags
les, spalding, twoliquidmixingfoam, wallfunctions


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
LES Filtering: how are the small and large scales equations solved? atmcfd Main CFD Forum 38 March 14, 2016 15:52
Beginning to work in LES - Suggestions needed atmcfd Main CFD Forum 5 July 20, 2012 23:16
Turbulence dampening due to magnetic field in LES and RAS eelcovv OpenFOAM 0 June 8, 2010 12:35
Differences between a laminar code and a les one ben Main CFD Forum 9 February 17, 2005 00:40
Some Questions about LES. Bin Li Main CFD Forum 2 February 20, 2004 10:58


All times are GMT -4. The time now is 23:55.