Go to the documentation of this file.
31 #include "twoPhaseSystem.H"
36 Foam::RASModels::kineticTheoryModel::kineticTheoryModel
43 const transportModel& phase,
44 const word& propertiesName,
67 kineticTheoryModels::viscosityModel::
New
74 kineticTheoryModels::conductivityModel::
New
81 kineticTheoryModels::radialModel::
New
86 granularPressureModel_
88 kineticTheoryModels::granularPressureModel::
New
93 frictionalStressModel_
95 kineticTheoryModels::frictionalStressModel::
New
101 equilibrium_(coeffDict_.
get<
bool>(
"equilibrium")),
103 alphaMax_(
"alphaMax",
dimless, coeffDict_),
104 alphaMinFriction_(
"alphaMinFriction",
dimless, coeffDict_),
105 residualAlpha_(
"residualAlpha",
dimless, coeffDict_),
112 IOobject::groupName(
"Theta", phase.
name()),
125 IOobject::groupName(
"lambda", phase.
name()),
139 IOobject::groupName(
"gs0", phase.
name()),
153 IOobject::groupName(
"kappa", phase.
name()),
167 IOobject::groupName(
"nuFric", phase.
name()),
177 if (
type == typeName)
202 coeffDict().readEntry(
"equilibrium", equilibrium_);
203 e_.readIfPresent(coeffDict());
204 alphaMax_.readIfPresent(coeffDict());
205 alphaMinFriction_.readIfPresent(coeffDict());
207 viscosityModel_->read();
208 conductivityModel_->read();
209 radialModel_->read();
210 granularPressureModel_->read();
211 frictionalStressModel_->read();
274 *granularPressureModel_->granularPressureCoeffPrime
277 radialModel_->g0(alpha_, alphaMinFriction_, alphaMax_),
278 radialModel_->g0prime(alpha_, alphaMinFriction_, alphaMax_),
282 + frictionalStressModel_->frictionalPressurePrime
290 volScalarField::Boundary& bpPrime =
291 tpPrime.
ref().boundaryFieldRef();
295 if (!bpPrime[patchi].
coupled())
297 bpPrime[patchi] == 0;
315 return devRhoReff(U_);
372 refCast<const twoPhaseSystem>(phase_.fluid()).otherPhase(phase_).U();
386 gs0_ = radialModel_->g0(
alpha, alphaMinFriction_, alphaMax_);
391 nut_ = viscosityModel_->nu(
alpha, Theta_, gs0_,
rho, da, e_);
396 lambda_ = (4.0/3.0)*
sqr(
alpha)*da*gs0_*(1 + e_)*ThetaSqrt/sqrtPi;
401 rho*(2*nut_*
D + (lambda_ - (2.0/3.0)*nut_)*
tr(
D)*
I)
410 *
rho*gs0_*(1.0/da)*ThetaSqrt/sqrtPi
416 refCast<const twoPhaseSystem>(phase_.fluid()).
Kd()
427 *sqrtPi*(ThetaSqrt + ThetaSmallSqrt)
434 granularPressureModel_->granularPressureCoeff
444 kappa_ = conductivityModel_->kappa(
alpha, Theta_, gs0_,
rho, da, e_);
466 +
fvm::Sp(J2/(Theta_ + ThetaSmall), Theta_)
485 (sqrtPi/(3*(3.0 - e_)))
486 *(1 + 0.4*(1 + e_)*(3*e_ - 1)*
alpha*gs0_)
487 +1.6*
alpha*gs0_*(1 + e_)/sqrtPi
494 4*da*
rho*(1 + e_)*
alpha*gs0_/(3*sqrtPi) - 2*
K3/3.0
517 *(2*
K3*trD2 +
K2*tr2D)
526 kappa_ = conductivityModel_->kappa(
alpha, Theta_, gs0_,
rho, da, e_);
534 nut_ = viscosityModel_->nu(
alpha, Theta_, gs0_,
rho, da, e_);
539 lambda_ = (4.0/3.0)*
sqr(
alpha)*da*gs0_*(1 + e_)*ThetaSqrt/sqrtPi;
544 frictionalStressModel_->frictionalPressure
552 nuFric_ = frictionalStressModel_->nu
563 nuFric_ =
min(nuFric_, maxNut_ - nut_);
569 Info<< typeName <<
':' <<
nl
570 <<
" max(Theta) = " <<
max(Theta_).value() <<
nl
571 <<
" max(nut) = " <<
max(nut_).value() <<
endl;
int debug
Static debugging option.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
virtual bool read()
Re-read model coefficients if they have changed.
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
A class for managing temporary objects.
virtual tmp< volScalarField > k() const
Return the turbulence kinetic energy.
static constexpr const zero Zero
Global zero (0)
Templated abstract base class for RAS turbulence models.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to database if not present.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
bool read(const char *buf, int32_t &val)
Same as readInt32.
ThermalDiffusivity< PhaseCompressibleTurbulenceModel< phaseModel > > phaseCompressibleTurbulenceModel
Typedef for phaseCompressibleTurbulenceModel.
virtual tmp< surfaceScalarField > pPrimef() const
Return the face-phase-pressure'.
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Ostream & endl(Ostream &os)
Add newline and flush stream.
const volScalarField Kd(fluid.Kd())
dimensionedSymmTensor dev2(const dimensionedSymmTensor &dt)
void min(const dimensioned< Type > &dt)
Use the minimum of the field and specified value.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
#define forAll(list, i)
Loop across all elements in list.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
#define NotImplemented
Issue a FatalErrorIn for a function not currently implemented.
tmp< fvMatrix< Type > > SuSp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
messageStream Info
Information stream (stdout output on master, null elsewhere)
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
static const SymmTensor I
void relax(const scalar alpha)
Relax matrix (for steady-state solution).
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
virtual void correct()
Solve the kinetic theory equations and correct the viscosity.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
const dimensionSet dimViscosity
const dimensionSet dimDynamicViscosity
Generic dimensioned Type class.
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.
virtual tmp< volScalarField > epsilon() const
Return the turbulence kinetic energy dissipation rate.
virtual tmp< volScalarField > pPrime() const
Return the phase-pressure'.
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
Global function forwards to reuseTmpDimensionedField::New.
GeometricField< vector, fvPatchField, volMesh > volVectorField
virtual tmp< volSymmTensorField > devRhoReff() const
Return the effective stress tensor.
virtual ~kineticTheoryModel()
Destructor.
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
RASModel< EddyDiffusivity< turbulenceModel > > RASModel
constexpr scalar pi(M_PI)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
dimensionedScalar sqrt(const dimensionedScalar &ds)
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
const dimensionedScalar e
Elementary charge.
const dimensionedScalar & D
A special matrix type and solver, designed for finite volume solutions of scalar equations....
virtual void printCoeffs(const word &type)
Print model coefficients.
SolverPerformance< Type > solve(const dictionary &)
Solve returning the solution statistics.
Eddy viscosity turbulence model base class.
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
virtual tmp< volScalarField > omega() const
Return the specific dissipation rate.
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
static word groupName(StringType base, const word &group)
Create dot-delimited name.group string.
bool coupled(solutionDict.getOrDefault("coupledEnergyField", false))
tmp< GeometricField< Type, fvPatchField, volMesh > > Sp(const volScalarField &sp, const GeometricField< Type, fvPatchField, volMesh > &vf)
virtual tmp< volSymmTensorField > R() const
Return the Reynolds stress tensor.
virtual tmp< fvVectorMatrix > divDevRhoReff(volVectorField &U) const
Return the source term for the momentum equation.
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
static const Identity< scalar > I
const dimensionSet dimless
Dimensionless.
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)