41 template<
class BasicMomentumTransportModel>
49 const viscosity& viscosity,
64 gasTurbulencePtr_(
nullptr),
106 if (type == typeName)
108 this->printCoeffs(type);
115 template<
class BasicMomentumTransportModel>
120 alphaInversion_.readIfPresent(this->coeffDict());
121 Cp_.readIfPresent(this->coeffDict());
122 C4_.readIfPresent(this->coeffDict());
123 Cmub_.readIfPresent(this->coeffDict());
134 template<
class BasicMomentumTransportModel>
138 if (!gasTurbulencePtr_)
143 refCast<const phaseModel>(this->properties());
145 const phaseModel& gas = fluid.otherPhase(liquid);
155 momentumTransportModel::typeName,
161 return *gasTurbulencePtr_;
165 template<
class BasicMomentumTransportModel>
169 this->gasTurbulence();
173 const phaseModel& gas = fluid.otherPhase(liquid);
176 this->Cmu_*
sqr(this->k_)/this->epsilon_
177 + Cmub_*gas.
d()*gasTurbulence.
alpha()
178 *(
mag(this->U_ - gasTurbulence.
U()));
180 this->nut_.correctBoundaryConditions();
185 template<
class BasicMomentumTransportModel>
189 this->gasTurbulence();
193 const phaseModel& gas = fluid.otherPhase(liquid);
217 template<
class BasicMomentumTransportModel>
229 max(alphaInversion_ - alpha, scalar(0))
236 template<
class BasicMomentumTransportModel>
243 this->gasTurbulence();
245 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
249 + phaseTransferCoeff*gasTurbulence.
k()
250 -
fvm::Sp(phaseTransferCoeff, this->k_);
254 template<
class BasicMomentumTransportModel>
262 this->gasTurbulence();
264 const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());
267 alpha*rho*this->C4_*this->epsilon_*bubbleG()/this->k_
268 + phaseTransferCoeff*gasTurbulence.
epsilon()
269 -
fvm::Sp(phaseTransferCoeff, this->epsilon_);
273 template<
class BasicMomentumTransportModel>
Generic thermophysical properties class for a liquid in which the functions and coefficients for each...
Templated abstract base class for multiphase compressible turbulence models.
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
virtual bool read()
Read model coefficients if they have changed.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
const alphaField & alpha() const
Access function to phase fraction.
Info<< "Reading strained laminar flame speed field Su\"<< endl;volScalarField Su(IOobject("Su", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field betav\"<< endl;volScalarField betav(IOobject("betav", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field Lobs\"<< endl;volScalarField Lobs(IOobject("Lobs", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field CT\"<< endl;volSymmTensorField CT(IOobject("CT", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field Nv\"<< endl;volScalarField Nv(IOobject("Nv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field nsv\"<< endl;volSymmTensorField nsv(IOobject("nsv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);IOdictionary PDRProperties(IOobject("PDRProperties", runTime.constant(), mesh, IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE));autoPtr< PDRDragModel > drag
const word & name() const
Generic dimensioned Type class.
tmp< volScalarField > d() const
Return the Sauter-mean diameter.
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
Continuous-phase k-epsilon model including bubble-generated turbulence.
LaheyKEpsilon(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const viscosity &viscosity, const word &type=typeName)
Construct from components.
Standard k-epsilon turbulence model for incompressible and compressible flows including rapid distort...
Class to represent a system of phases and model interfacial transfers between them.
A class for handling words, derived from string.
const rhoThermo & thermo() const
Return const-access to phase rhoThermo.
static word groupName(Name name, const word &group)
Info<< "Reading field p\"<< endl;volScalarField p(IOobject("p", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);pressureReference pressureReference(p, simple.dict());mesh.schemes().setFluxRequired(p.name());Info<< "Reading field pa\"<< endl;volScalarField pa(IOobject("pa", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field Ua\"<< endl;volVectorField Ua(IOobject("Ua", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);# 65 "/home/ubuntu/OpenFOAM-10/applications/solvers/incompressible/adjointShapeOptimisationFoam/createFields.H" 2label paRefCell=0;scalar paRefValue=0.0;setRefCell(pa, simple.dict(), paRefCell, paRefValue);mesh.schemes().setFluxRequired(pa.name());autoPtr< viscosityModel > viscosity(viscosityModel::New(mesh))
BasicMomentumTransportModel::alphaField alphaField
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
Class to represent a interface between phases where one phase is considered dispersed within the othe...
virtual tmp< volScalarField > CdRe() const =0
Drag coefficient.
static autoPtr< dictionary > New(Istream &)
Construct top-level dictionary on freestore from Istream.
tmp< volScalarField > phaseTransferCoeff() const
virtual tmp< volScalarField > k() const =0
Return the turbulence kinetic energy.
Abstract base class for turbulence models (RAS, LES and laminar).
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
virtual tmp< fvScalarMatrix > kSource() const
dimensionedScalar pow3(const dimensionedScalar &ds)
const volVectorField & U() const
Access function to velocity field.
const Time & time() const
Return time.
const phaseSystem & fluid() const
Return the system to which this phase belongs.
BasicMomentumTransportModel::rhoField rhoField
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
virtual tmp< volScalarField > nu() const
Kinematic viscosity of mixture [m^2/s].
tmp< volScalarField > bubbleG() const
virtual void correctNut()
dimensioned< scalar > mag(const dimensioned< Type > &)
virtual tmp< fvScalarMatrix > epsilonSource() const
A class for managing temporary objects.
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
const objectRegistry & db() const
Return the local objectRegistry.
virtual tmp< volScalarField > epsilon() const =0
Return the turbulence kinetic energy dissipation rate.
Single incompressible phase derived from the phase-fraction. Used as part of the multiPhaseMixture fo...
dimensionedScalar deltaT() const
Return time step.