[pm11dt]

I'm relatively new to openFoam and require some help with a specific issue I'm having when dealing with the post processing of LES cases.

I need to calculate k and epsilon fields for my LES case in order to compare with a RANS K-epsilon case I ran previously. In order to achieve this I have been using the following 3rd party utility:

postAverageTurbulenceFields found here: Utilities: post average turbulence fields and create turbulence fields for LES

However, I believe this utility only calls back and writes values for k from the transport equation solved during simulation. This would be fine when used with RANS but not for LES. I used the oneEqEddy viscosity sgs model which has a transport equation for ksgs and writes ksgs as k in the time directory.

When i use this tool with oneEqEddy LES model it will write ksgs as k and is probably writing epsilon from the sgs based on the same assumptions.

I need a way of calculating the total k and total epsilon for my LES case by either modifying postAverageTurbulenceFields utility or using some other utility to calculate these values from the mean velocity information I have.

k=0.5<U'_i><U'_i>
epsilon= nu(d<U'_i>/d<x_j>)^2

nu = nuSgs + nu of material ..?

Any help would be very much appreciated, but I should mention I have very little programming experience i.e. actually modifying source code. I am using and older version of openFoam, v2.3.1.

Thanks!

Daniel.