31#include "twoPhaseSystem.H"
38 const volScalarField&
alpha,
39 const volScalarField&
rho,
40 const volVectorField&
U,
41 const surfaceScalarField& alphaRhoPhi,
42 const surfaceScalarField&
phi,
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
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;
Templated abstract base class for single-phase compressible turbulence models.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
static word groupName(StringType base, const word &group)
Create dot-delimited name.group string.
Templated abstract base class for RAS turbulence models.
Kinetic theory particle phase RAS model.
virtual tmp< volScalarField > epsilon() const
Return the turbulence kinetic energy dissipation rate.
virtual tmp< volSymmTensorField > R() const
Return the Reynolds stress tensor.
virtual void correct()
Solve the kinetic theory equations and correct the viscosity.
virtual tmp< volScalarField > k() const
Return the turbulence kinetic energy.
virtual tmp< volScalarField > pPrime() const
Return the phase-pressure'.
virtual tmp< surfaceScalarField > pPrimef() const
Return the face-phase-pressure'.
virtual tmp< volScalarField > omega() const
Return the specific dissipation rate.
virtual tmp< volSymmTensorField > devRhoReff() const
Return the effective stress tensor.
virtual ~kineticTheoryModel()
Destructor.
virtual bool read()
Re-read model coefficients if they have changed.
static const SymmTensor I
static autoPtr< Time > New()
Construct (dummy) Time - no functionObjects or libraries.
Generic dimensioned Type class.
Eddy viscosity turbulence model base class.
tmp< volVectorField > divDevRhoReff()
Return the effective viscous stress (laminar + turbulent).
A special matrix type and solver, designed for finite volume solutions of scalar equations....
void relax(const scalar alpha)
Relax matrix (for steady-state solution).
SolverPerformance< Type > solve(const dictionary &)
Solve returning the solution statistics.
A class for managing temporary objects.
bool coupled(solutionDict.getOrDefault("coupledEnergyField", false))
#define NotImplemented
Issue a FatalErrorIn for a function not currently implemented.
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
RASModel< EddyDiffusivity< turbulenceModel > > RASModel
constexpr scalar pi(M_PI)
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.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh > > grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
tmp< GeometricField< Type, fvPatchField, volMesh > > Sp(const volScalarField &sp, const GeometricField< Type, fvPatchField, volMesh > &vf)
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
zeroField SuSp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
zeroField Sp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
dimensionedSymmTensor dev2(const dimensionedSymmTensor &dt)
const dimensionSet dimViscosity
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
const dimensionSet dimless
Dimensionless.
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
bool read(const char *buf, int32_t &val)
Same as readInt32.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
messageStream Info
Information stream (stdout output on master, null elsewhere)
static const Identity< scalar > I
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
Ostream & endl(Ostream &os)
Add newline and flush stream.
dimensionedScalar sqrt(const dimensionedScalar &ds)
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
ThermalDiffusivity< PhaseCompressibleTurbulenceModel< phaseModel > > phaseCompressibleTurbulenceModel
Typedef for phaseCompressibleTurbulenceModel.
const dimensionSet dimDynamicViscosity
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh > > &tdf1, const word &name, const dimensionSet &dimensions)
Global function forwards to reuseTmpDimensionedField::New.
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
constexpr char nl
The newline '\n' character (0x0a)
const dimensionedScalar & D
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
#define forAll(list, i)
Loop across all elements in list.