34 namespace compressible
48 const geometricOneField&
alpha,
53 const fluidThermo& thermophysicalModel,
70 dimensioned<scalar>::lookupOrAddToDict
82 PDRkEpsilon::~PDRkEpsilon()
92 C4_.readIfPresent(coeffDict_);
107 nut_ = Cmu_*
sqr(k_)/epsilon_;
108 nut_.correctBoundaryConditions();
131 epsilon_.boundaryFieldRef().updateCoeffs();
150 tmp<fvScalarMatrix> epsEqn
156 C1_*betav*
G*epsilon_/k_
157 + 1.5*
pow(Cmu_, 3.0/4.0)*GR*
sqrt(k_)/LI
162 epsEqn.ref().relax();
164 epsEqn.ref().boundaryManipulate(epsilon_.boundaryFieldRef());
167 bound(epsilon_, epsilonMin_);
172 tmp<fvScalarMatrix> kEqn
188 nut_ = Cmu_*
sqr(k_)/epsilon_;
189 nut_.correctBoundaryConditions();
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
tmp< fvMatrix< Type > > SuSp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
dimensionedScalar sqrt(const dimensionedScalar &ds)
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
GeometricField< vector, fvPatchField, volMesh > volVectorField
virtual bool read()
Read model coefficients if they have changed.
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
PDRkEpsilon(const geometricOneField &alpha, const volScalarField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const fluidThermo &thermophysicalModel, const word &modelName=typeName)
Construct from components.
const compressible::RASModel & turbulence_
Macros for easy insertion into run-time selection tables.
const volScalarField & rho_
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
const volVectorField & U_
GeometricField< scalar, fvPatchField, volMesh > volScalarField
CompressibleMomentumTransportModel< fluidThermo > momentumTransportModel
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
const surfaceScalarField & phi_
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Info<< "Predicted p max-min : "<< max(p).value()<< " "<< min(p).value()<< endl;rho==max(rho0+psi *p, rhoMin);# 1 "/home/ubuntu/OpenFOAM-8/applications/solvers/multiphase/cavitatingFoam/alphavPsi.H" 1{ alphav=max(min((rho - rholSat)/(rhovSat - rholSat), scalar(1)), scalar(0));alphal=1.0 - alphav;Info<< "max-min alphav: "<< max(alphav).value()<< " "<< min(alphav).value()<< endl;psiModel-> correct()
defineTypeNameAndDebug(combustionModel, 0)
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
PDRDragModel(const dictionary &PDRProperties, const compressible::RASModel &turbulence, const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi)
Construct from components.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
volScalarField & bound(volScalarField &, const dimensionedScalar &lowerBound)
Bound the given scalar field if it has gone unbounded.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
SolverPerformance< Type > solve(fvMatrix< Type > &, const word &)
Solve returning the solution statistics given convergence tolerance.
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const objectRegistry & db() const
Return the local objectRegistry.
const dimensionedScalar & G
Newtonian constant of gravitation.
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
RASModel< momentumTransportModel > RASModel
Typedefs for turbulence, RAS and LES models for compressible flow based on the standard laminar trans...