Commit baa02e69 authored by Henry Weller's avatar Henry Weller

src/combustionModels: Updated LES model lookup

parent 142ff927
......@@ -184,7 +184,7 @@ public:
// Member Functions
//- Read LESProperties dictionary
//- Read model coefficients if they have changed
virtual bool read();
//- Return the effective diffusivity for nuTilda
......
......@@ -146,7 +146,7 @@ public:
// Member Functions
//- Read LESProperties dictionary
//- Read model coefficients if they have changed
virtual bool read();
};
......
......@@ -189,7 +189,10 @@ void FSD<CombThermoType, ThermoType>::calculateSourceNorm()
// Calculation of the mixture fraction variance (ftVar)
const compressible::LESModel& lesModel =
YO2.db().lookupObject<compressible::LESModel>("LESProperties");
YO2.db().lookupObject<compressible::LESModel>
(
turbulenceModel::propertiesName
);
const volScalarField& delta = lesModel.delta();
const volScalarField ftVar(Cv_*sqr(delta)*sqr(mgft));
......
......@@ -52,7 +52,7 @@ Description
is large (>1e-04) then a beta pdf is used for filtering.
At the moment the flame area combustion model is only fit to work in a LES
frame work. In RAS the subgrid fluctiuation has to be solved by an extra
frame work. In RAS the subgrid fluctuation has to be solved by an extra
transport equation.
SourceFiles
......
......@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
......@@ -75,7 +75,6 @@ void Foam::reactionRateFlameAreaModels::relaxation::correct
const volScalarField& sigma
)
{
dimensionedScalar omega0
(
"omega0",
......@@ -97,13 +96,19 @@ void Foam::reactionRateFlameAreaModels::relaxation::correct
1e-4
);
const compressible::LESModel& lesModel =
omega_.db().lookupObject<compressible::LESModel>("LESProperties");
dimensionedScalar kMin
(
"kMin",
sqr(dimVelocity),
SMALL
);
const compressibleTurbulenceModel& turbulence = combModel_.turbulence();
// Total strain : resolved and sub-grid (just LES for now)
// Total strain
const volScalarField sigmaTotal
(
sigma + alpha_*lesModel.epsilon()/(lesModel.k() + lesModel.kMin())
sigma + alpha_*turbulence.epsilon()/(turbulence.k() + kMin)
);
const volScalarField omegaInf(correlation_.omega0Sigma(sigmaTotal));
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment