26 #include "SchaefferFrictionalStress.H" 33 namespace kineticTheoryModels
35 namespace frictionalStressModels
41 frictionalStressModel,
54 const dictionary&
dict 57 frictionalStressModel(dict),
58 coeffDict_(dict.subDict(typeName +
"Coeffs")),
59 phi_(
"phi",
dimless, coeffDict_)
84 *
pow(
Foam::max(alpha1 - alphaMinFriction, scalar(0)), 10.0);
99 *
pow(
Foam::max(alpha1 - alphaMinFriction, scalar(0)), 9.0);
113 const scalar I2Dsmall = 1.0e-15;
117 tmp<volScalarField> tnu
124 alpha1.mesh().time().timeName(),
139 if (alpha1[celli] > alphaMinFriction.value())
144 sqrt(1.0/6.0*(
sqr(D[celli].xx() - D[celli].yy())
145 +
sqr(D[celli].yy() - D[celli].zz())
146 +
sqr(D[celli].zz() - D[celli].xx()))
147 +
sqr(D[celli].xy()) +
sqr(D[celli].xz())
148 +
sqr(D[celli].yz())) + I2Dsmall
154 nuf.correctBoundaryConditions();
164 phi_.
read(coeffDict_);
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
Schaeffer(const dictionary &dict)
Construct from components.
#define forAll(list, i)
Loop across all elements in list.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const dictionary & subDict(const word &) const
Find and return a sub-dictionary.
dimensionedScalar sqrt(const dimensionedScalar &ds)
const Type & value() const
Return const reference to value.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
virtual tmp< volScalarField > nu(const phaseModel &phase, const dimensionedScalar &alphaMinFriction, const dimensionedScalar &alphaMax, const volScalarField &pf, const volSymmTensorField &D) const
void read(const dictionary &)
Update the value of dimensioned<Type>
dimensionedScalar sin(const dimensionedScalar &ds)
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
defineTypeNameAndDebug(combustionModel, 0)
Internal & ref()
Return a reference to the dimensioned internal field.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
virtual ~Schaeffer()
Destructor.
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
const dictionary & dict_
Reference to higher-level dictionary for re-read.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
virtual tmp< volScalarField > frictionalPressurePrime(const phaseModel &phase, const dimensionedScalar &alphaMinFriction, const dimensionedScalar &alphaMax) const
A class for managing temporary objects.
virtual tmp< volScalarField > frictionalPressure(const phaseModel &phase, const dimensionedScalar &alphaMinFriction, const dimensionedScalar &alphaMax) const