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()),
116 IOobject::AUTO_WRITE,
126 IOobject::groupName(
"lambda", phase.
name()),
140 IOobject::groupName(
"gs0", phase.
name()),
154 IOobject::groupName(
"kappa", phase.
name()),
168 IOobject::groupName(
"nuFric", phase.
name()),
172 IOobject::AUTO_WRITE,
179 if (
type == typeName)
204 coeffDict().readEntry(
"equilibrium", equilibrium_);
205 e_.readIfPresent(coeffDict());
206 alphaMax_.readIfPresent(coeffDict());
207 alphaMinFriction_.readIfPresent(coeffDict());
209 viscosityModel_->read();
210 conductivityModel_->read();
211 radialModel_->read();
212 granularPressureModel_->read();
213 frictionalStressModel_->read();
270 *granularPressureModel_->granularPressureCoeffPrime
273 radialModel_->g0(alpha_, alphaMinFriction_, alphaMax_),
274 radialModel_->g0prime(alpha_, alphaMinFriction_, alphaMax_),
278 + frictionalStressModel_->frictionalPressurePrime
287 tpPrime.ref().boundaryFieldRef();
291 if (!bpPrime[patchi].
coupled())
293 bpPrime[patchi] == 0;
311 return devRhoReff(U_);
362 refCast<const twoPhaseSystem>(phase_.fluid()).otherPhase(phase_).U();
377 gs0_ = radialModel_->g0(
alpha, alphaMinFriction_, alphaMax_);
382 nut_ = viscosityModel_->nu(
alpha, Theta_, gs0_,
rho, da, e_);
387 lambda_ = (4.0/3.0)*
sqr(
alpha)*da*gs0_*(1 + e_)*ThetaSqrt/sqrtPi;
392 rho*(2*nut_*
D + (lambda_ - (2.0/3.0)*nut_)*
tr(
D)*
I)
401 *
rho*gs0_*(1.0/da)*ThetaSqrt/sqrtPi
407 refCast<const twoPhaseSystem>(phase_.fluid()).
Kd()
418 *sqrtPi*(ThetaSqrt + ThetaSmallSqrt)
425 granularPressureModel_->granularPressureCoeff
435 kappa_ = conductivityModel_->kappa(
alpha, Theta_, gs0_,
rho, da, e_);
457 +
fvm::Sp(J2/(Theta_ + ThetaSmall), Theta_)
476 (sqrtPi/(3*(3.0 - e_)))
477 *(1 + 0.4*(1 + e_)*(3*e_ - 1)*
alpha*gs0_)
478 +1.6*
alpha*gs0_*(1 + e_)/sqrtPi
485 4*da*
rho*(1 + e_)*
alpha*gs0_/(3*sqrtPi) - 2*
K3/3.0
508 *(2*
K3*trD2 +
K2*tr2D)
517 kappa_ = conductivityModel_->kappa(
alpha, Theta_, gs0_,
rho, da, e_);
520 Theta_.clamp_range(0, 100);
524 nut_ = viscosityModel_->nu(
alpha, Theta_, gs0_,
rho, da, e_);
529 lambda_ = (4.0/3.0)*
sqr(
alpha)*da*gs0_*(1 + e_)*ThetaSqrt/sqrtPi;
534 frictionalStressModel_->frictionalPressure
542 nuFric_ = frictionalStressModel_->nu
554 nuFric_ =
min(nuFric_, maxNut_ - nut_);
560 Info<< typeName <<
':' <<
nl 561 <<
" max(Theta) = " <<
max(Theta_).value() <<
nl 562 <<
" 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.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
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.
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)
const volScalarField Kd(fluid.Kd())
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< volScalarField > trho
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)
Generic dimensioned Type class.
const dimensionSet dimless
Dimensionless.
Eddy viscosity turbulence model base class.
GeometricField< vector, fvPatchField, volMesh > volVectorField
const dimensionSet dimDynamicViscosity
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.
Templated abstract base class for RAS turbulence models.
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.
const dimensionedScalar e
Elementary charge.
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)
A special matrix type and solver, designed for finite volume solutions of scalar equations. Face addressing is used to make all matrix assembly and solution loops vectorise.
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.
static tmp< GeometricField< symmTensor, fvPatchField, volMesh > > New(const word &name, IOobjectOption::registerOption regOpt, const Mesh &mesh, const dimensionSet &dims, const word &patchFieldType=fvPatchField< symmTensor >::calculatedType())
Return tmp field (NO_READ, NO_WRITE) from name, mesh, dimensions and patch type. [Takes current timeN...
SolverPerformance< Type > solve(const dictionary &)
Solve returning the solution statistics.
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 relax(const scalar alpha)
Relax matrix (for steady-state solution).
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.
Do not request registration (bool: false)
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)