35namespace PDRDragModels
47 const dictionary& PDRProperties,
49 const volScalarField&
rho,
50 const volVectorField&
U,
51 const surfaceScalarField&
phi
55 Csu(
"Csu",
dimless, PDRDragModelCoeffs_),
56 Csk(
"Csk",
dimless, PDRDragModelCoeffs_),
63 U_.
mesh().facesInstance(),
76 U_.
mesh().facesInstance(),
96 tmp<volSymmTensorField> tDragDcu
103 U_.mesh().time().constant(),
121 (0.5*rho_)*CR_*
mag(U_) + (Csu*
I)*
betav*turbulence_.muEff()*
sqr(Aw_);
130 tmp<volScalarField> tGk
137 U_.mesh().time().constant(),
158 (0.5*rho_)*
mag(U_)*(U_ & CT & U_)
170 PDRDragModelCoeffs_.readEntry(
"Csu", Csu.value());
171 PDRDragModelCoeffs_.readEntry(
"Csk", Csk.value());
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
Internal & ref(const bool updateAccessTime=true)
Return a reference to the dimensioned internal field.
const objectRegistry & db() const noexcept
Return the local objectRegistry.
virtual bool read()
Inherit read from regIOobject.
Basic sub-grid obstacle drag model. Details supplied by J Puttock 2/7/06.
virtual ~basic()
Destructor.
virtual tmp< volSymmTensorField > Dcu() const
Return the momentum drag coefficient.
virtual tmp< volScalarField > Gk() const
Return the momentum drag turbulence generation rate.
void writeFields() const
Write fields.
virtual bool read()
Re-read model coefficients if they have changed.
const Type & lookupObject(const word &name, const bool recursive=false) const
A class for managing temporary objects.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
compressible::turbulenceModel & turbulence
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const dimensionSet dimless
Dimensionless.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensionedScalar pow3(const dimensionedScalar &ds)
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
static const Identity< scalar > I
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
dimensioned< symmTensor > dimensionedSymmTensor
Dimensioned tensor obtained from generic dimensioned type.
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
static constexpr const zero Zero
Global zero (0)
const dimensionSet dimVolume(pow3(dimLength))
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
const volScalarField & betav