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

van Driest delta kills performance?

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 18, 2016, 05:11
Default van Driest delta kills performance?
  #1
Senior Member
 
Timofey Mukha
Join Date: Mar 2012
Location: Stockholm, Sweden
Posts: 118
Rep Power: 14
tiam is on a distinguished road
Hello!

I have been running LES of channel flow for the meshes of following sizes: 1, 2, 4 and 8 million.
For each mesh I have 2 case which differ only in one aspect: one uses cubeRootVol delta, and the other one van Driest.

I decompose the cases in such a way that each processor has around 30k cells, which as far as I know a reasonably large amount.

For the cubeRootVol cases I have very nice scaling in the sense that I can compute (almost) the same simulation time for all the mesh sizes for a given wall-time.

However, for the van Driest cases I have problems. For the 1 and 2 millions meshes the performance is pretty much on par with cubeRootVol. But for 4 and 8 million, the difference becomes very significant.

For the 4 million case I have 0.52 hours per simulation-second for the cubeRootVol and for the van Driest 0.75 hours, so about 1.45 times slower.

For the 8 million case I have 0.47 h/sim-sec for the cubeRootVol (better than the 4 mill case !) and for the van Driest delta I have 1.76 h/sim-sec!!

As I said the cases are identical besides for the choice delta. Also by observing the log as the simulation runs I see that I have a pause in the end of each time-step, after the U and p solvers are done. As far as I know that is where turbulent quantities get updated.

I am using version 2.3.1.

Has anyone encountered similar issues? Is this to be expected? As I understand the van Driest delta has to calculate the wall distance. This should be done only once for a static mesh. Maybe the issue is that it is recalculated every time?

Best,
Timofey
tiam is offline   Reply With Quote

Old   January 18, 2016, 12:02
Default
  #2
Senior Member
 
Timofey Mukha
Join Date: Mar 2012
Location: Stockholm, Sweden
Posts: 118
Rep Power: 14
tiam is on a distinguished road
As a follow up.
Indeed, the distance is evaluated at every time-step.
This is necessary, since we are interested in considering cells where y+<500 or some similar value.
Since y+ gets updated every time-step we have to also recalculate the distance.

I am wondering if this is needed, since the form of the damping function makes it dissappear away from the wall anyway.
So it should be possible to just caluclate the whole distance field once and then use it, whithout worrying about the damping acting only near the wall.
Or am I missing something?
tiam is offline   Reply With Quote

Old   January 25, 2016, 03:27
Default
  #3
Senior Member
 
Timofey Mukha
Join Date: Mar 2012
Location: Stockholm, Sweden
Posts: 118
Rep Power: 14
tiam is on a distinguished road
An update.

The behaviour has now also been confirmed using a larger mesh (69 mil.) and a standard scaling test, i.e. doubling the amount of cores and comparing to the increase in the simulated time.

Fixing this is not very easy it seems. The damping function depends on y+, calculated for the interiour cells. On a general geometry it is not trivial which u_tau to pick up to do the scaling for each cell. And one is forced to recalculate this field every time-step.
tiam is offline   Reply With Quote

Old   February 10, 2018, 06:21
Default
  #4
Member
 
Join Date: Nov 2013
Posts: 35
Rep Power: 12
beatlejuice is on a distinguished road
Hello tiam,

I have exactly the same problem!

I have to identical LES Smagorinsky cases which differ only in the choice of delta: cubeRootVol & van Driest.

I also observe a very long pause after one time step is solved (after solving for U & p and preparing for next time step) for the van Driest method. This pause does not exist in the cubeRootVol simulation.

Did you find any work around or does anybody has any solution to this.

I use OF 2.3.1

cheers!
beatlejuice is offline   Reply With Quote

Old   February 10, 2018, 06:28
Default
  #5
Senior Member
 
Timofey Mukha
Join Date: Mar 2012
Location: Stockholm, Sweden
Posts: 118
Rep Power: 14
tiam is on a distinguished road
Quote:
Originally Posted by beatlejuice View Post
Hello tiam,

I have exactly the same problem!

I have to identical LES Smagorinsky cases which differ only in the choice of delta: cubeRootVol & van Driest.

I also observe a very long pause after one time step is solved (after solving for U & p and preparing for next time step) for the van Driest method. This pause does not exist in the cubeRootVol simulation.

Did you find any work around or does anybody has any solution to this.

I use OF 2.3.1

cheers!
Hi! No, there is now workaround I think. I would consider using a different SGS model, WALE for instance. This should be looked i to by the devs.
tiam is offline   Reply With Quote

Old   February 10, 2018, 13:50
Default
  #6
Member
 
Join Date: Nov 2013
Posts: 35
Rep Power: 12
beatlejuice is on a distinguished road
Hey tiam,

thx for your fast reply!

Did you use the WALE model? Were the results satisfying?

Do you know what the field B is?
What boundary condition did you use for the walls, inlet, outlet for the fields B, k, nut, nuTilda?

I would appreciate yout help!
Thanks! Beatle
beatlejuice is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Smoothing of IDDES delta Fuchs OpenFOAM Running, Solving & CFD 2 November 17, 2020 12:17
van Driest 1956 , classic paper Ferramis Main CFD Forum 3 January 18, 2017 10:16
Incomp. LES in pisoFOAM, how to set up Smagorinsky model with van Driest dampin Jan.Östh OpenFOAM Running, Solving & CFD 21 August 10, 2016 13:03
looking for the paper "van driest transformation and compressible wall bounded flows" hnemati Main CFD Forum 0 February 4, 2015 08:07
Implementation of Van Driest Damping Function MaLa OpenFOAM Running, Solving & CFD 1 May 25, 2014 22:15


All times are GMT -4. The time now is 16:43.