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)
198 RASModel<EddyDiffusivity<phaseCompressibleTurbulenceModel>>
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();
247 return tmp<volSymmTensorField>
271 tmp<volScalarField> tpPrime
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_);
325 return tmp<volSymmTensorField>
367 const twoPhaseSystem&
fluid = refCast<const twoPhaseSystem>(phase_.fluid());
378 tmp<volScalarField> tda(phase_.d());
381 tmp<volTensorField> tgradU(
fvc::grad(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.0 + e_)*ThetaSqrt/sqrtPi;
401 rho*(2.0*nut_*
D + (lambda_ - (2.0/3.0)*nut_)*
tr(
D)*
I)
410 *
rho*gs0_*(1.0/da)*ThetaSqrt/sqrtPi
416 fluid.lookupSubModel<dragModel>
419 fluid.otherPhase(phase_)
431 *sqrtPi*(ThetaSqrt + ThetaSmallSqrt)
438 granularPressureModel_->granularPressureCoeff
448 kappa_ = conductivityModel_->kappa(
alpha, Theta_, gs0_,
rho, da, e_);
470 +
fvm::Sp(J2/(Theta_ + ThetaSmall), Theta_)
489 (sqrtPi/(3.0*(3.0 - e_)))
490 *(1.0 + 0.4*(1.0 + e_)*(3.0*e_ - 1.0)*
alpha*gs0_)
491 +1.6*
alpha*gs0_*(1.0 + e_)/sqrtPi
498 4.0*da*
rho*(1.0 + e_)*
alpha*gs0_/(3.0*sqrtPi) - 2.0*
K3/3.0
521 *(2.0*
K3*trD2 +
K2*tr2D)
530 kappa_ = conductivityModel_->kappa(
alpha, Theta_, gs0_,
rho, da, e_);
533 Theta_.clamp_range(0, 100);
537 nut_ = viscosityModel_->nu(
alpha, Theta_, gs0_,
rho, da, e_);
542 lambda_ = (4.0/3.0)*
sqr(
alpha)*da*gs0_*(1.0 + e_)*ThetaSqrt/sqrtPi;
547 frictionalStressModel_->frictionalPressure
555 nuFric_ = frictionalStressModel_->nu
567 nuFric_ =
min(nuFric_, maxNut_ - nut_);
573 Info<< typeName <<
':' <<
nl 574 <<
" max(Theta) = " <<
max(Theta_).value() <<
nl 575 <<
" max(nut) = " <<
max(nut_).value() <<
endl;
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
List of values generated by applying the access operation to each list item.
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
fvMatrix< scalar > fvScalarMatrix
ThermalDiffusivity< PhaseCompressibleTurbulenceModel< phaseModel > > phaseCompressibleTurbulenceModel
Typedef for phaseCompressibleTurbulenceModel.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
GeometricField< tensor, fvPatchField, volMesh > volTensorField
dimensionedSymmTensor sqr(const dimensionedVector &dv)
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
constexpr char nl
The newline '\n' character (0x0a)
virtual tmp< volSymmTensorField > R() const
Return the Reynolds stress tensor.
const dimensionSet dimViscosity
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
tmp< GeometricField< Type, fvPatchField, volMesh > > Sp(const volScalarField &sp, const GeometricField< Type, fvPatchField, volMesh > &vf)
GeometricBoundaryField< scalar, fvPatchField, volMesh > Boundary
Type of boundary fields.
Ignore writing from objectRegistry::writeObject()
const dimensionSet dimless
Dimensionless.
GeometricField< vector, fvPatchField, volMesh > volVectorField
CGAL::Exact_predicates_exact_constructions_kernel K
virtual bool read()
Re-read model coefficients if they have changed.
bool read(const char *buf, int32_t &val)
Same as readInt32.
#define forAll(list, i)
Loop across all elements in list.
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
virtual tmp< fvVectorMatrix > divDevRhoReff(volVectorField &U) const
Return the source term for the momentum equation.
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
static word groupName(StringType base, const word &group)
Create dot-delimited name.group string.
virtual void correct()
Solve the kinetic theory equations and correct the viscosity.
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for expressions::valueTypeCode::INVALID.
virtual ~kineticTheoryModel()
Destructor.
static const SymmTensor I
static const Identity< scalar > I
virtual tmp< surfaceScalarField > pPrimef() const
Return the face-phase-pressure'.
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
constexpr scalar pi(M_PI)
zeroField SuSp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
virtual void printCoeffs(const word &type)
Print model coefficients.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
zeroField Sp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
virtual tmp< volScalarField > omega() const
Return the specific dissipation rate.
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
virtual tmp< volScalarField > k() const
Return the turbulence kinetic energy.
int debug
Static debugging option.
virtual tmp< volSymmTensorField > devRhoReff() const
Return the effective stress tensor.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
RASModel< EddyDiffusivity< turbulenceModel > > RASModel
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.
void clamp_max(const Type &upper)
Impose upper (ceiling) clamp on the field values (in-place)
virtual tmp< volScalarField > pPrime() const
Return the phase-pressure'.
dimensionedSymmTensor dev2(const dimensionedSymmTensor &dt)
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
messageStream Info
Information stream (stdout output on master, null elsewhere)
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
const dimensionedScalar & D
A class for managing temporary objects.
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
static options & New(const fvMesh &mesh)
Construct fvOptions and register to database if not present.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
#define NotImplemented
Issue a FatalErrorIn for a function not currently implemented.
virtual tmp< volScalarField > epsilon() const
Return the turbulence kinetic energy dissipation rate.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
SymmTensor< Cmpt > devTwoSymm(const SymmTensor< Cmpt > &st)
Return the deviatoric part of twice the symmetric part of a SymmTensor.
static constexpr const zero Zero
Global zero (0)