[ravicfdd]

Hello everyone, I am using OpenFOAM-9 and I have a question regarding lewis number models. Starting from OpenFOAM-8 there was a redesign to the multi-component diffusion modelling. New Transport Models that does not use unity Lewis number approximation have been introduced.

I am performing combustion simulations and I am running a case file and when the case runs the solver outputs:

Creating turbulence model
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity
Creating reaction model

I notice that it selects a unityLewisEddyDiffusivity model.


My turbulenceProperties/momentumTransport file look like this:

simulationType RAS;

RAS
{
RASModel kOmegaSST;
turbulence on;
printCoeffs no;

kOmegaSSTCoeffs
{
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.55555;
gamma2 0.5;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1.0;
c1 10.0;
F3 false;
}
}

I have found another thread similar to my question : reactingFoam non unity Lewis number

But according to this it looks like I can not use the kOmegaSST model while wanting a non unity lewis number. Any solution to this will be of great help! Thanks!

[ravicfdd]

Update: I have solved it myself, an additional file "thermophysicalTransport" is required in the constant folder to specify nonUnityLewis.


It looks like this:

FoamFile
{
format ascii;
class dictionary;
location "constant";
object thermophysicalTransport;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

RAS
{
model nonUnityLewisEddyDiffusivity;

Prt 0.85;
Sct 0.7;
}


// ************************************************** *********************** //