37 namespace kineticTheoryModels
39 namespace frictionalStressModels
45 frictionalStressModel,
46 JohnsonJacksonSchaeffer,
59 const dictionary&
dict 62 frictionalStressModel(
dict),
63 coeffDict_(
dict.optionalSubDict(typeName +
"Coeffs")),
64 Fr_(
"Fr", dimensionSet(1, -1, -2, 0, 0), coeffDict_),
65 eta_(
"eta",
dimless, coeffDict_),
67 phi_(
"phi",
dimless, coeffDict_),
68 alphaDeltaMin_(
"alphaDeltaMin",
dimless, coeffDict_)
87 const phaseModel& phase,
95 Fr_*
pow(
max(
alpha - alphaMinFriction, scalar(0)), eta_)
104 const phaseModel& phase,
113 eta_*
pow(
max(
alpha - alphaMinFriction, scalar(0)), eta_ - 1.0)
115 + p_*
pow(
max(
alpha - alphaMinFriction, scalar(0)), eta_)
124 const phaseModel& phase,
140 auto& nuf = tnu.ref();
147 0.5*pf[celli]*
sin(phi_.value())
166 pf.boundaryField()[patchi]*
sin(phi_.value())
168 mag(
U.boundaryField()[patchi].snGrad())
176 nuf.correctBoundaryConditions();
185 coeffDict_ <<= dict_.optionalSubDict(typeName +
"Coeffs");
187 Fr_.read(coeffDict_);
188 eta_.read(coeffDict_);
191 phi_.read(coeffDict_);
194 alphaDeltaMin_.read(coeffDict_);
Cmpt invariantII(const SymmTensor< Cmpt > &st)
Return the 2nd invariant of a SymmTensor.
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
const Type & value() const noexcept
Return const reference to value.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Unit conversion functions.
dimensionedScalar sqrt(const dimensionedScalar &ds)
GeometricBoundaryField< scalar, fvPatchField, volMesh > Boundary
Type of boundary fields.
const dimensionSet dimless
Dimensionless.
GeometricField< vector, fvPatchField, volMesh > volVectorField
addToRunTimeSelectionTable(frictionalStressModel, JohnsonJackson, dictionary)
Macros for easy insertion into run-time selection tables.
#define forAll(list, i)
Loop across all elements in list.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
static word scopedName(const std::string &scope, const word &name)
Create scope:name or scope_name string.
dimensionedScalar alphaMax("alphaMax", dimless/dimTime, laminarTransport)
defineTypeNameAndDebug(JohnsonJackson, 0)
virtual tmp< volScalarField > nu(const phaseModel &phase, const dimensionedScalar &alphaMinFriction, const dimensionedScalar &alphaMax, const volScalarField &pf, const volSymmTensorField &D) const
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
static tmp< GeometricField< scalar, fvPatchField, volMesh > > New(const word &name, IOobjectOption::registerOption regOpt, const Mesh &mesh, const dimensionSet &dims, const word &patchFieldType=fvPatchField< scalar >::calculatedType())
Return tmp field (NO_READ, NO_WRITE) from name, mesh, dimensions and patch type. [Takes current timeN...
dimensionedScalar sin(const dimensionedScalar &ds)
virtual ~JohnsonJacksonSchaeffer()
Destructor.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const polyBoundaryMesh & patches
virtual tmp< volScalarField > frictionalPressurePrime(const phaseModel &phase, const dimensionedScalar &alphaMinFriction, const dimensionedScalar &alphaMax) const
virtual tmp< volScalarField > frictionalPressure(const phaseModel &phase, const dimensionedScalar &alphaMinFriction, const dimensionedScalar &alphaMax) const
const dimensionedScalar & D
A class for managing temporary objects.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
constexpr scalar degToRad(const scalar deg) noexcept
Conversion from degrees to radians.
JohnsonJacksonSchaeffer(const dictionary &dict)
Construct from components.
Do not request registration (bool: false)
PtrList< fvPatch > fvPatchList
Store lists of fvPatch as a PtrList.
static constexpr const zero Zero
Global zero (0)