
[Sponsors] 
pisoFOAM (LES)  internal pipe flow  convergence 

LinkBack  Thread Tools  Search this Thread  Display Modes 
July 31, 2018, 09:03 
pisoFOAM (LES)  internal pipe flow  convergence

#1  
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 
Hello,
First, I'm using OpenFOAM 5.0 (LES  Smagorinsky  Re=5300 with R = 0,1 and L=10D)... I start the field in RANS (steadyState  simpleFoam) and later I mapped (mapFields) it to a new folder with intention to start the simulations in LES (pisoFoam). The problem lies in my results and in the way I am evaluating the statistically stationary state of turbulent flow. A piece of my log.pisoFoam: Quote:
Quote:
Unfortunately this is my first case in turbulence, I do not know how to evaluate if the flow is in a statistically stationary state (... kindly ask those who have experience, share with me). I tried to monitor the 'wallShearStress' on the wall, but without success (similar to how it is done in Fluent), could anyone tell me if this is possible? probes?! The maximum that I did was to consider the region of entrance of the flow, through the relation: L = 0.06 * Re * D. With possession of this value (L), and because it is a periodic domain with flow velocity equal to 1,472 m/s (flow in the Z direction)... it was determined that the time to reach the region of 'developed flow' was 43.20 seconds ... for that reason the simulation was set in 60 seconds. The fluid was perturbed by the code (perturbUCylinder) and evaluating the velocities at points in the domain, there was no laminarization (attached image  pt4: point near the wall). Thanks to all who can help. 

July 31, 2018, 13:00 

#2 
Senior Member
Join Date: Aug 2015
Posts: 433
Rep Power: 11 
For the convergence, statistical stationarity refers to how they're (the statistics) changing. Simply put, your statistics will not change if you have obtained convergence (provided that you are averaging over enough data). Given that you're modeling turbulent flow, the least you'll want to do is make sure your average of interest is not changing. The period/time you average over depends on the problem, but generally longer is better. This will of course increase as you become interested in higher order statistics. It is hard to tell how your results could be off without seeing the DNS results on the plots in addition to your own.
Caelan 

July 31, 2018, 14:14 

#3  
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 
Quote:
Thanks for your answer. The results that I am looking for to compare can be found in: "A direct numerical simulation study on the mean velocity characteristics in turbulent pipe flow  Wu (2008)", I will attach image of the RMS profiles that I am seeking to obtain, these are the ones that are most I'm intrigued by my simulation because they are totally different and are not varying over time (with longer simulation time). I wanted to understand what could be causing this ... could you help me? On the statistically stationary question, I understand the concept ... I want to understand what parameters I should evaluate in my simulation using the OF. As I mentioned, using Fluent just monitor the 'wallShearStress' on the wall ... however, I have not found out if this is possible in OF. 

July 31, 2018, 15:03 

#4 
Senior Member
Join Date: Aug 2015
Posts: 433
Rep Power: 11 
Can you plot your results on top of the DNS values?
For the statistics, here's a nice web page with a description : https://people.duke.edu/~rnau/411diff.htm. So you'd at least look at the timeaverage of your quantity of interest, but higher order statistics such as variance are also important. I am not sure exactly what fluent does, but wall shear stress is a function of velocity gradients near the wall as well as the effective viscosity  ie it will depend on turbulence quantities and treatment of the boundary layer. If you want to match that in openfoam, there's a function object called wallShearStress. So you could include it in your controlDict to use at runTime or as a post processing step. In terms of matching DNS with LES, there are many possible ways that the two could differ  for example, your choice of LES model, resolution of the boundary layer, resolution of the domain. So, how are you treating the boundary layer? Caelan 

July 31, 2018, 15:36 

#5  
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 
Quote:
The wallShearStress function is being used in my controlDict, what I can not do is monitor this function throughout the simulation to evaluate 'convergence'. The resolution on the wall is y + <1. 

July 31, 2018, 15:50 

#6 
Senior Member
Join Date: Aug 2015
Posts: 433
Rep Power: 11 
They are plotting against (1r+) on the x axis  where are they defining r from? Make sure you're using the same independent variable. If your first cell is below y+ of 1 you should have adequate resolution in the boundary layer  make sure you aren't applying any wall modeling.
Depending on what you have available, you don't need to use native linux tools to plot. I always read data into matlab for post processing. Try running postProcess func wallShearStress. This will create a series of text files for you to read into your software of choice (in a folder called postProcessing  this might actually be created already when you include the function object in the controlDict). Otherwise, there are various posts on the forum about pulling residual data from a log file to plot with gnuplot. Caelan 

August 1, 2018, 13:39 

#7  
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 
Quote:
Some things I'm taking, I do not know if they will work. On the question of the wall, as I said, the resolution is of y + <1. As soon as I get the new results, I'll let you know. 

August 9, 2018, 13:16 

#8 
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 
I was thinking that the problem I am facing may be related to the disturbance being imposed on the velocity field .... OR being related to my fvSchenes, because the field only keeps fluctuations when I use the 'Gauss cubic' model for velocity.... what to me is very strange, 'Gauss limitedLinear 1.0' and 'Gauss LinearUpwind grad (U)' are laminarizing my fluid.


August 9, 2018, 13:36 

#9 
Senior Member
Join Date: Aug 2015
Posts: 433
Rep Power: 11 
That is a good observation  limitedLinear and linearUpwind are nominally first/second order and second order upwind, respectively. From a highlevel perspective, this means that (possibly undue) numerical dissipation will be introduced by the schemes when compared to central differencing (linear) for the same grid. Generally for LES you'll want to use Gauss linear for your divergence terms as long as stability is not a concern.
Caelan 

August 9, 2018, 13:48 

#10  
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 
Quote:


August 9, 2018, 13:57 

#11 
Senior Member
Join Date: Aug 2015
Posts: 433
Rep Power: 11 
What temporal integration scheme are you using? For LES 2nd order in time is also desirable.
Caelan 

August 9, 2018, 15:06 

#12  
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 
Quote:
I'm alternating to see if I find '' differences '' between one setup and another. I have used 'backward', but I also did tests using 'CrankNicolson 0.9'. Is my problem, because it is a pipe, directly related to the results with the cartesian coodenates? Is it necessary to convert them to cylinder coodenates in order to "get" the actual values of my problem (if they exist)? Guilherme 

August 9, 2018, 15:11 

#13 
Senior Member
Join Date: Aug 2015
Posts: 433
Rep Power: 11 
Good, it sounds like your time integration should be fine. When you say there's dissipation, are you seeing this when you visualize, or when you plot? Related to your latest question, I'd recommend solving/plotting like the DNS data that you are comparing to. I am unfamiliar with using cylindrical coordinates in OF, though.
Caelan 

August 9, 2018, 15:27 

#14  
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 
Quote:


August 9, 2018, 15:28 

#15 
Senior Member
Join Date: Aug 2015
Posts: 433
Rep Power: 11 
Have you visualized the flow to confirm this?
Caelan 

August 9, 2018, 16:22 

#16 
Senior Member
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 310
Rep Power: 12 
If you have non zero velocity variance (autoccorelations, R. Stress) you have turbulence. Your plots are not appropriately scaled, you need to normalize by the shear velocity if you want to compare. Besides, you have to average not only in time but in directions of homogeneity, that is, in polar coordinates along radial lines.


August 9, 2018, 17:33 

#17 
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 

August 9, 2018, 17:35 

#18  
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 
Quote:
Thanks for the answer, Guilherme *EDIT1: Do you find it necessary to convert from cartesian coordinates to cylindrical coordinates? 

August 9, 2018, 18:09 

#19 
Senior Member
Join Date: Aug 2015
Posts: 433
Rep Power: 11 
I think we've established that your case looks to be fine numerically. That leaves the setup and analysis/post processing. As mentioned earlier, you should be analyzing your data so you can plot the same way the DNS data is plotted.
As for setup, here is a similar thread that will hopefully provide useful insight : question regarding LES of pipe flow  pimpleFoam Caelan 

August 10, 2018, 08:00 

#20 
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 180
Rep Power: 5 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
question regarding LES of pipe flow  pimpleFoam  Dan1788  OpenFOAM Running, Solving & CFD  37  December 26, 2017 15:42 
Velocity curve in channel flow by pisoFoam and LES  mxylondon  OpenFOAM Verification & Validation  9  June 21, 2015 08:28 
LES transitional pipe entrance flow  hasanduz  Main CFD Forum  2  October 11, 2013 18:59 
setup problems  LES pipe flow with cyclic BC (1) and direct mapped inlet (2)  florian_krause  OpenFOAM  22  June 13, 2013 22:25 
convergence problem when use pisoFoam, LES for wind tunnel case  Forrest_Lei  OpenFOAM  3  July 19, 2011 07:00 