37namespace kineticTheoryModels
39namespace frictionalStressModels
56Foam::kineticTheoryModels::frictionalStressModels::
57JohnsonJacksonSchaeffer::JohnsonJacksonSchaeffer
63 coeffDict_(
dict.optionalSubDict(typeName +
"Coeffs")),
65 eta_(
"eta",
dimless, coeffDict_),
67 phi_(
"phi",
dimless, coeffDict_),
68 alphaDeltaMin_(
"alphaDeltaMin",
dimless, coeffDict_)
95 Fr_*
pow(
max(
alpha - alphaMinFriction, scalar(0)), eta_)
113 eta_*
pow(
max(
alpha - alphaMinFriction, scalar(0)), eta_ - 1)
115 + p_*
pow(
max(
alpha - alphaMinFriction, scalar(0)), eta_)
137 "JohnsonJacksonSchaeffer:nu",
150 0.5*pf[celli]*
sin(phi_.value())
171 mag(
U.boundaryField()[patchi].snGrad())
188 coeffDict_ <<= dict_.optionalSubDict(typeName +
"Coeffs");
190 Fr_.read(coeffDict_);
191 eta_.read(coeffDict_);
194 phi_.read(coeffDict_);
197 alphaDeltaMin_.read(coeffDict_);
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
const Mesh & mesh() const
Return mesh.
Boundary & boundaryFieldRef(const bool updateAccessTime=true)
Return a reference to the boundary field.
const Boundary & boundaryField() const
Return const-reference to the boundary field.
void correctBoundaryConditions()
Correct boundary field.
A list of pointers to objects of type <T>, with allocation/deallocation management of the pointers....
static autoPtr< Time > New()
Construct (dummy) Time - no functionObjects or libraries.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Dimension set for the base types, which can be used to implement rigorous dimension checking for alge...
const Type & value() const
Return const reference to value.
virtual tmp< volScalarField > frictionalPressure(const phaseModel &phase, const dimensionedScalar &alphaMinFriction, const dimensionedScalar &alphaMax) const
virtual tmp< volScalarField > frictionalPressurePrime(const phaseModel &phase, const dimensionedScalar &alphaMinFriction, const dimensionedScalar &alphaMax) const
virtual ~JohnsonJacksonSchaeffer()
Destructor.
Single incompressible phase derived from the phase-fraction. Used as part of the multiPhaseMixture fo...
Single incompressible phase derived from the phase-fraction. Used as part of the multiPhaseMixture fo...
A class for managing temporary objects.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
const polyBoundaryMesh & patches
bool coupled(solutionDict.getOrDefault("coupledEnergyField", false))
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
const dimensionSet dimless
Dimensionless.
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Cmpt invariantII(const SymmTensor< Cmpt > &st)
Return the 2nd invariant of a SymmTensor.
dimensionedScalar sin(const dimensionedScalar &ds)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
constexpr scalar degToRad() noexcept
Multiplication factor for degrees to radians conversion.
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const dimensionedScalar & D
dimensionedScalar alphaMax("alphaMax", dimless/dimTime, laminarTransport)
#define forAll(list, i)
Loop across all elements in list.
Unit conversion functions.