[HelioVillanueva]

Hi all.

I implemented the Partially-Averaged Navier-Stokes (PANS) methodology in OpenFOAM-3.0.x for the kEpsilon and kOmegaSST turbulence models. I didn't rigorously validate the code, just tested on the pitzDaily tutorial for LES and the results seems to agree well for the mean velocity but with the big advantage in processing time. While the LES tutorial took 1hr and 24 mins, the kEpsilonPANS and kOmegaSSTPANS took aprox 15 mins. Although it's a simple tutorial case, the results are promising. Here is a video comparing the models.

Since the implementation was done using a lot of information available online (credits below), I hope that the validation and improvement of the code could also be done with the community help. For now, there is the need to copy the k and epsilon/omega files in 0 folder and rename to kU and epsilonU/omegaU to run a case. I think it would be nice if don't, but I couldn't manage to initialize the kU and epsilonU/omegaU fields with the NO_READ option in a way that was possible to use the turbulentMixingLenghtFrequencyInlet and turbulentMixingLenghtDissipationRateInlet bcs.

Another thing is the fK value. In the papers some use a constant value or a fK field. The second approach is implemented but it is possible to use a constant fK by specifying a narrow limit in constant/turbulenceProperties. See the case example provided with the code.

Code: here
it should compile just coping the directories to your OF user folder and run the Allmake script.

Hakan Nillson's phd cfd course at Chalmers. Report by Guglielmo Minelli
http://www.tfd.chalmers.se/~hani/kur.../OF_report.pdf

Abdol-Hamid, Girimaji 2005
http://ntrs.nasa.gov/archive/nasa/ca...0050182919.pdf

Lakshmipathy, Girimaji 2006
http://arc.aiaa.org/doi/abs/10.2514/6.2006-119

SUNIL LAKSHMIPATHY PHD Thesis:
http://oaktrust.library.tamu.edu/bit...pdf?sequence=2

About implementing new turbulence models in OF-3.0.x:
http://www.cfd-online.com/Forums/ope...tml#post578372