35 incompressibleMomentumTransportModel
42 namespace incompressible
52 RASincompressibleMomentumTransportModel,
65 return exp((-scalar(2.5) + Rt/20.0)/
pow3(scalar(1) + Rt/130.0));
70 exp(-6.0/
sqr(scalar(1) + Rt/50.0))
71 *(scalar(1) + 3.0*
exp(-Rt/10.0));
78 tmp<volScalarField> Rt =
q_*
k_/(2.0*nu()*
zeta_);
79 return scalar(1) - 0.3*
exp(-
sqr(Rt));
94 const geometricOneField&
alpha,
95 const geometricOneField&
rho,
103 eddyViscosity<incompressible::RASModel>
160 qMin_(
"qMin",
sqrt(kMin_)),
161 zetaMin_(
"zetaMin", epsilonMin_/(2*qMin_)),
167 this->groupName(
"k"),
180 this->groupName(
"epsilon"),
193 this->groupName(
"q"),
200 k_.boundaryField().types()
207 this->groupName(
"zeta"),
213 max(epsilon_, epsilonMin_)/(2.0*q_),
214 epsilon_.boundaryField().types()
217 bound(zeta_, zetaMin_);
219 if (
type == typeName)
263 tmp<fvScalarMatrix> zetaEqn
274 zetaEqn.ref().relax();
280 tmp<fvScalarMatrix> qEqn
Bound the given scalar field where it is below the specified minimum.
void correctBoundaryConditions()
Correct boundary field.
bool readIfPresent(const word &, const dictionary &)
Update the value of the Switch if it is found in the dictionary.
static Switch lookupOrAddToDict(const word &, dictionary &, const Switch &defaultValue=false)
Construct from dictionary, supplying default value so that if the.
static dimensioned< scalar > lookupOrAddToDict(const word &, dictionary &, const dimensionSet &dims=dimless, const scalar &defaultValue=pTraits< scalar >::zero)
Construct from dictionary, with default value.
bool readIfPresent(const dictionary &)
Update the value of dimensioned<Type> if found in the dictionary.
virtual bool read()=0
Re-read model coefficients if they have changed.
virtual void correct()=0
Solve the turbulence equations and correct the turbulence viscosity.
virtual void correctNut()
tmp< volScalarField > f2() const
qZeta(const geometricOneField &alpha, const geometricOneField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const viscosity &viscosity, const word &type=typeName)
Construct from components.
tmp< volScalarField > DqEff() const
Return the effective diffusivity for q.
dimensionedScalar sigmaZeta_
dimensionedScalar qMin_
Lower limit of q.
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
virtual bool read()
Read RASProperties dictionary.
dimensionedScalar zetaMin_
Lower limit of zeta.
tmp< volScalarField > fMu() const
tmp< volScalarField > DzetaEff() const
Return the effective diffusivity for epsilon.
Calculate the magnitude of the square of the gradient of the gradient of the given volField.
volScalarField alpha(IOobject("alpha", runTime.name(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
const dimensionedScalar G
Newtonian constant of gravitation.
tmp< volScalarField > magSqrGradGrad(const VolField< Type > &vf)
tmp< VolField< typename outerProduct< vector, Type >::type > > grad(const SurfaceField< Type > &ssf)
tmp< fvMatrix< Type > > laplacian(const VolField< Type > &vf, const word &name)
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const VolField< Type > &)
tmp< fvMatrix< Type > > SuSp(const volScalarField::Internal &, const VolField< Type > &)
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const VolField< Type > &vf, const word &name)
tmp< fvMatrix< Type > > ddt(const VolField< Type > &vf)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
dimensionedScalar exp(const dimensionedScalar &ds)
VolField< vector > volVectorField
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar pow3(const dimensionedScalar &ds)
SurfaceField< scalar > surfaceScalarField
dimensionedScalar sqrt(const dimensionedScalar &ds)
VolField< scalar > volScalarField
defineTypeNameAndDebug(combustionModel, 0)
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
bool bound(volScalarField &, const dimensionedScalar &min)
Bound the given scalar field where it is below the specified min value.
dimensioned< scalar > magSqr(const dimensioned< Type > &)
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
SolverPerformance< Type > solve(fvMatrix< Type > &, const word &)
Solve returning the solution statistics given convergence tolerance.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
makeMomentumTransportModelTypes(geometricOneField, geometricOneField, incompressibleMomentumTransportModel)
Info<< "Reading field U\n"<< endl;volVectorField U(IOobject("U", runTime.name(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Creating face flux\n"<< endl;surfaceScalarField phi(IOobject("phi", runTime.name(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedScalar(mesh.Sf().dimensions() *U.dimensions(), 0));autoPtr< viscosityModel > viscosity(viscosityModel::New(mesh))