[Geovany]

Hi everyone,

I'm currently working with a 3D cavity-based geometry using rhoSimpleFoam solver. I'd like to get the u_rms value of my simulations. I've already done my simulation in a 5M mesh and I've validated my results.

Reading another thread I know that I have to re-run my simulation and use fieldAverage function object to get UPrime2Mean and then u_rms,x=sqrt(UPrime2Meanxx) using calculator in paraview.

My question is: In a steady simulations, 'cause I used a steady solver, the answer of u_rms,x that I will get is correct? Or it is just averaging through every iteration and this is not useful, since the solution has not yet converged until the last iteration?
In case that the approach of using fieldAverage isn't correct, someone can tell me how I could get the u_rms using a steady solver?

Thanks in advance!

[HPE]

No, don't think you can use `fieldAverage` on steady simulations. There can't be any `UPrime2Mean` because there is no `UPrime` (i.e. fluctuating part) for steady simulations.

The object you are looking for is `R` for steady cases. You can use turbulenceFields function object to calculate it, and I don't think you need to rerun your completed RANS simulation. Postprocessing the simulation should be enough.

[Geovany]

Thank you very much for your help.

I used the function that you mentioned and then I was able to plot the turbulence properties in Matlab. But I have another question, I hope that you could help me. How can I compute u'_rms from the Reynolds Stress Tensor?
=
or
=

I think that in the latter option it is assuming that the velocity fluctuations are equal to and that only applies to isotropic turbulence, right? Or am I mixing basic concepts?

My goal is to plot cross-stream profiles of turbulent fluctations

[HPE]

Yes, "turbulenceFields" function object assumes isotropic turbulence for the computation of the turbulence intensity.

This assumption is fine for the linear eddy-viscosity turbulence models, because one of the derivation assumptions of these models is the turbulence isotropy in freestream.

Having said that, if you still want to use the Reynolds stress tensor (you can use "k" safely for the isotropic configs), you need to have access to "R" field - and pick the orthogonal components.

For a given "R" field, the order of components in OpenFOAM is like the following:

(Rxx Rxy Rxz Ryy Ryz Rzz)

Hope this helps.