Relative roughness in openfoam
How to incorporate relative roughness in openfoam?
I am using pisoFoam to solve a simple turbulent flow through a pipe. I want to validate the result using Moody's chart, where i specify reynolds number and relative roughness to openfoam and check if the frcition factor given by openfoam matches with the value from moody's chart. Do we have to use nutkRoughWallFunction? If yes then what do Ks and Cs mean, and where do i find their values? Also, what is the default relative roughness used by openfoam? |
Hi,
I strongly recommend you to study this article: Blocken, B., Stathopoulos, T., Carmeliet, J., 2007. CFD simulation of the atmospheric boundary layer: wall function problems. Atmospheric Environment 41 (2), 238–252. Also check this thread out: http://www.cfd-online.com/Forums/ope...h-surface.html You can see nutkRoughWallFunction source file to see how are they defined. Cs values is something between 0 to 1, but typically it is considered 0.5. Ks (It is called equivalent sand-gran roughness height) differs from case to case, it is often calculated by: Ks = (20 to 30)*y0 in which, y0 is aerodynamic roughness length. a rule of thumb for values of y0 is 0.1 of the real roughness value. for example if height of the roughness is 0.1 m then y0 would be 0.01 m. I hope it helps a bit. Best |
Okay i just realised that my way of validating turbulent flow is wrong.
Can you tell me how to validate a turbulent flow using some analytical results, like we have poiseuille's law for laminar flow. |
Quote:
hope it helps. |
Quote:
How do i calculate relative roughness from Ks and Cs values? And how do i get the friction factor from openfoam? I read the nutkRoughWallFunction files and understood what Ks and Cs mean but there is no mention of how to obtain relative roughness from them. |
Quote:
OK, first of all, what do you mean by relative roughness? do you mean y0? I haven't worked with nutkRoughWallFunction for a while and forgot some, I think you define Ks values as an input data. Well friction factor can be obtained using wall shear stress, as you know: Cf=Tau/(0.5*rho*U^2) In which Tau is wall shear stress. You can calculate wall shear stress with a utility named as "wallShearStress" in openfoam. later you can calcluate Cf in paraview using its calculator according to rho and U values. |
Quote:
Also, Cf is skin friction factor. I believe it is not the same as the friction factor f which is on the left hand scale of moody's chart. According to texts - friction factor is calculated experimentally. So the approach i am using is wrong. |
1 Attachment(s)
There is one more problem though.
This is the test case i am trying out. Just to learn about turbulence. But even at 10^6 reynolds number i am getting laminar streamlines. Can you look up and tell me what is wrong with the problem. Assume smooth walls. Attachment 24567 |
Quote:
I couldn't get why your approach is wrong. |
Quote:
Now my conclusion is - Obtain Fanning friction coefficient Cf using wallShearStress and then multiply it by 4 to get Darcy friction factor and use Moodys chart with Reynolds number and relative pipe roughness to verify the friction factor. Thanks a lot. I think this solves my problem. Can you also take a look at another problem which i posted right before your current reply? |
Quote:
|
Quote:
Thanks in advance :) |
1 Attachment(s)
Quote:
The results look similar. The velocity is not blowing. Input bein 1 m/s i am getting a maximum velocity of around 1.98 m/s. The results still look laminar - Attachment 24777 |
Quote:
|
1 Attachment(s)
Quote:
These are my yPLus values Attachment 24826 Also, I am using k-epsilon turbulence model. |
Quote:
You have got to use kOmegaSST model alongside with a finer grid to be able to resolve the boundary layer. Try take a look at these: http://www.computationalfluiddynamic...t-cell-height/ http://www.computationalfluiddynamic...oundary-layer/ http://www.computationalfluiddynamic...oundary-layer/ http://www.computationalfluiddynamic...ds-number-cfd/ http://www.computationalfluiddynamic...-requirements/ http://www.computationalfluiddynamic...nce-modelling/ |
Quote:
Also, for kOmegaSST what boundary conditions should i use? Is it possible for you to simulate the problem and send me the files? I have been trying this problem for 3 weeks now and not able to understand anything. |
Quote:
you have to make a new initial boundary file named omega in your 0 directory and use omegaWallFunction for your wall boundaries. For example: wall { type omegaWallFunction; value uniform 0; } outlet { type zeroGradient; } inlet { type fixedValue; value uniform 0.0001; } try having a look at this thread: http://www.cfd-online.com/Forums/ope...omega-sst.html I hope it helps a bit :), Best |
Quote:
One more question, do the values of k and omega (or epsilon) affect the final solution? I was trying out various values but found that for k>epsilon i get floating point exception on the courant number and for k<<epsilon (whatever be the value of k and epsilon) I obtain solutions which are very very close to each other. from what i have read y+ should be between 30-300 right? Also I was looking at pitzDaily example from pisoFoam>LES, there you can see vortices but on running the same simulation with say kE you do not get the vortices. I read somewhere on the forum that LES is better at capturing vortices as with kE or kOmega the viscous forces dampen the vortices. |
Quote:
Quote:
Quote:
Therefore by using kOmegaSST model, you have got to lower your y+ value to something lower than 1. Quote:
Right now I have found this article which is quite useful: www.engmech.cz/2012/proceedings/pdf/195_Furst_J-FT.pdf It uses openFOAM and a new turbulence model which is best for laminar turbulent transition modeling. And plus take a look at this: http://www.cfd-online.com/Forums/ope...el-low-re.html I am studying about this. You have got to give me some time. Dig them out, maybe you can get what I can't. ;) I hope it helps a bit, Best, Mojtaba |
1 Attachment(s)
Thanks.
I will have a look at the links. Quote:
I do not want to go into the complications of LES since I will be using kOmegaSST only. The problem i am facing now is that in the pitzDaily case does not show continuous vortices as it showed with LES. I am using kOmegaSST turbulence model. I am using wall functions and my flow has a very high reynolds number. Attachment 24891 |
Quote:
http://www.engmech.cz/2012/proceedin...Furst_J-FT.pdf Quote:
Provided links has a comprehensive procedure to calculate first cell height. Right now I am working on kkLOmega model which is best for capturing laminar turbulent transition. I will report if I get anything useful. |
les nut
Quote:
|
All times are GMT -4. The time now is 06:45. |