79 const scalar a =
aZone_[i];
88 for (label celli : cz)
95 kCoeff.correctBoundaryConditions();
111 typeName +
":epsilonCoeff",
112 mesh_.time().timeName(),
125 const scalar a = aZone_[i];
126 const scalar
N = NZone_[i];
127 const scalar Cep = CepZone_[i];
128 const scalar Cd = CdZone_[i];
130 for (label zonei : zoneIDs_[i])
132 const cellZone& cz = mesh_.cellZones()[zonei];
134 for (label celli : cz)
136 epsilonCoeff[celli] = Cep*Cd*a*
N*
mag(
U[celli]);
143 return tepsilonCoeff;
152 const word& modelType,
165 epsilonName_(
"epsilon")
181 if (eqn.
psi().name() == epsilonName_)
189 else if (eqn.
psi().name() == kName_)
209 if (eqn.
psi().name() == epsilonName_)
217 else if (eqn.
psi().name() == kName_)
232 if (!coeffs_.readIfPresent(
"epsilonNames", fieldNames_))
234 if (coeffs_.found(
"epsilon"))
236 fieldNames_.resize(1);
237 coeffs_.readEntry(
"epsilon", fieldNames_.first());
241 fieldNames_.resize(2);
242 fieldNames_[0] =
"epsilon";
243 fieldNames_[1] =
"k";
249 const dictionary& regionsDict(coeffs_.subDict(
"regions"));
263 const word zoneName(modelDict.
get<
word>(
"cellZone"));
265 zoneIDs_[i] = mesh_.cellZones().indices(zoneName);
266 if (zoneIDs_[i].empty())
269 <<
"Unable to find cellZone " << zoneName <<
nl
270 <<
"Valid cellZones are:" << mesh_.cellZones().names()
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.
void correctBoundaryConditions()
Correct boundary field.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
virtual bool read()
Re-read model coefficients if they have changed.
static autoPtr< Time > New()
Construct (dummy) Time - no functionObjects or libraries.
static word timeName(const scalar t, const int precision=precision_)
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
T get(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
const dictionary & subDict(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
Find and return a sub-dictionary.
wordList toc() const
Return the table of contents.
bool readEntry(const word &keyword, T &val, enum keyType::option matchOpt=keyType::REGEX, bool mandatory=true) const
A special matrix type and solver, designed for finite volume solutions of scalar equations....
const GeometricField< Type, fvPatchField, volMesh > & psi(const label i=0) const
Return psi.
Mesh data needed to do the Finite Volume discretisation.
const Time & time() const
Return the top-level database.
labelListList zoneIDs_
Zone indices.
tmp< volScalarField > kCoeff(const volVectorField &U) const
Return the k coefficient.
virtual void addSup(fvMatrix< scalar > &eqn, const label fieldi)
Add implicit contribution to momentum equation.
virtual bool read(const dictionary &dict)
Read dictionary.
scalarList aZone_
Width of the vegetation element.
scalarList NZone_
Number of plants per unit of area.
tmp< volScalarField > epsilonCoeff(const volVectorField &U) const
Return the epsilon coefficient.
scalarList CdZone_
Drag coefficient.
Base abstract class for handling finite volume options (i.e. fvOption).
const fvMesh & mesh_
Reference to the mesh database.
void resetApplied()
Resize/reset applied flag list for all fieldNames_ entries.
const cellZoneMesh & cellZones() const noexcept
Return cell zone mesh.
A class for managing temporary objects.
A class for handling words, derived from Foam::string.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
Foam::word regionName(Foam::polyMesh::defaultRegion)
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
A special matrix type and solver, designed for finite volume solutions of scalar equations.
Calculate the finiteVolume matrix for implicit and explicit sources.
zeroField Sp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
const dimensionSet dimless
Dimensionless.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
static constexpr const zero Zero
Global zero (0)
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
errorManipArg< error, int > exit(error &err, const int errNo=1)
constexpr char nl
The newline '\n' character (0x0a)
#define forAll(list, i)
Loop across all elements in list.
const Vector< label > N(dict.get< Vector< label > >("N"))