29 #include "phasePair.H" 49 const dictionary&
dict,
71 sqr(pair_.dispersed().d())
72 /pair_.continuous().nu()
78 min(
Re).value() < 1200.0
79 ||
max(
Re).value() > 18800.0
80 ||
min(sqrSr).value() < 0.0016
81 ||
max(sqrSr).value() > 0.04
85 <<
"Re and/or Sr are out of the range of applicability of the " 86 <<
"Moraga model. Clamping to range bounds" 91 sqrSr.clamp_range(0.0016, 0.04);
93 return 0.2*
exp(-
Re*sqrSr/3.6e5 - 0.12)*
exp(
Re*sqrSr/3.0e7);
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
addToRunTimeSelectionTable(liftModel, constantLiftCoefficient, dictionary)
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)
defineTypeNameAndDebug(constantLiftCoefficient, 0)
Ostream & endl(Ostream &os)
Add newline and flush stream.
void clamp_range(const Type &lower, const Type &upper)
Clamp field values (in-place) to the specified range.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
virtual tmp< volScalarField > Cl() const
Lift coefficient.
dimensionedScalar exp(const dimensionedScalar &ds)
Calculate the gradient of the given field.
Moraga(const dictionary &dict, const phasePair &pair)
Construct from a dictionary and a phase pair.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
virtual ~Moraga()
Destructor.
scalarField Re(const UList< complex > &cmplx)
Extract real component.
#define WarningInFunction
Report a warning using Foam::Warning.
A class for managing temporary objects.