LES calculation of EPSILON (turbulence dissipation)
I am using pisoFOAM with LES Smagorinsky model to simulate a fully developed turbulent flow in a rotating pipe.
Now i want to evaluate the turbulence dissipation epsilon.
Is is it possible to calculate or postprocess it with openfoam? Is there anyone who did some kind of postprocessing before?
Any help would be appreciated.
Hi Simon and everyone;
Did you find the solution to your problem?
I am also trying to calculate de dissipation of an LES simulation, but I am not getting the correct value.
I am doing a decaying isotropic turbulence, then the derivative of the energy with respect to time should be equal to my -dissipation: dE/dt=-epsilon.
If I calculate dE/dt in my simulations, my results are in agreement with a DNS simulation (A. Wray in the AGARD database). However if I calculate epsilon directly from my velocity field, it seems that I have a constant "shift" as you can see in my plot.
I am calculating my disipation as:
volScalarField epsilon=2*(nu+nuSgs)*(S && S); #epsilon total
Any suggestions or ideas of what I might be doing wrong ? Am I missing something?
Thanks for your help and have a nice day,
|All times are GMT -4. The time now is 12:37.|