46 const label patchi =
patch().index();
53 internalField().
group()
66 auto&
nutw = tnutw.ref();
71 const scalar t = db().time().timeOutputValue();
75 for (
const scalar z : z0)
80 <<
"z0 field can only contain positive values. " 81 <<
"Please check input field z0." 94 const scalar uStar = Cmu25*
sqrt(
k[celli]);
95 const scalar
yPlus = uStar*
y[facei]/nuw[facei];
96 const scalar Edash = (
y[facei] + z0[facei])/z0[facei];
129 const DimensionedField<scalar, volMesh>& iF
147 boundNut_(ptf.boundNut_),
148 z0_(ptf.z0_.clone(
p.
patch()))
160 boundNut_(
dict.getOrDefault<bool>(
"boundNut", false)),
171 boundNut_(rwfpsf.boundNut_),
183 boundNut_(rwfpsf.boundNut_),
195 nutkWallFunctionFvPatchScalarField::autoMap(m);
210 nutkWallFunctionFvPatchScalarField::rmap(ptf, addr);
212 const auto& nrwfpsf =
213 refCast<const atmNutkWallFunctionFvPatchScalarField>(ptf);
217 z0_->rmap(nrwfpsf.z0_(), addr);
235 atmNutkWallFunctionFvPatchScalarField
virtual void write(Ostream &) const
Write.
dimensionedScalar log(const dimensionedScalar &ds)
errorManipArg< error, int > exit(error &err, const int errNo=1)
This boundary condition provides a wall constraint on the turbulent viscosity (i.e. nut) based on the turbulent kinetic energy (i.e. k) for atmospheric boundary layer modelling. It is designed to be used in conjunction with the atmBoundaryLayerInletVelocity boundary condition.
virtual tmp< scalarField > calcNut() const
Calculate the turbulent viscosity.
error FatalError
Error stream (stdout output on all processes), with additional 'FOAM FATAL ERROR' header text and sta...
A list of keyword definitions, which are a keyword followed by a number of values (eg...
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedScalar sqrt(const dimensionedScalar &ds)
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
dimensionedScalar pow025(const dimensionedScalar &ds)
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
void writeValueEntry(Ostream &os) const
Write *this field as a "value" entry.
label k
Boltzmann constant.
Smooth ATC in cells next to a set of patches supplied by type.
const dimensionedScalar kappa
Coulomb constant: default SI units: [N.m2/C2].
Abstract base class for turbulence models (RAS, LES and laminar).
Macros for easy insertion into run-time selection tables.
#define forAll(list, i)
Loop across all elements in list.
constexpr const char *const group
Group name for atomic constants.
static word groupName(StringType base, const word &group)
Create dot-delimited name.group string.
const dimensionedScalar e
Elementary charge.
This boundary condition provides a wall function for the turbulent viscosity (i.e. nut) based on the turbulent kinetic energy, (i.e. k) for for low- and high-Reynolds number applications.
fvPatchField< scalar > fvPatchScalarField
static const word propertiesName
Default name of the turbulence properties dictionary.
A FieldMapper for finite-volume patch fields.
static tmp< T > New(Args &&... args)
Construct tmp with forwarding arguments.
virtual tmp< scalarField > yPlus() const
Calculate and return the yPlus at the boundary.
Ostream & writeEntryIfDifferent(const word &key, const T &value1, const T &value2)
Write a keyword/value entry only when the two values differ.
virtual void rmap(const fvPatchScalarField &, const labelList &)
Reverse map the given fvPatchField onto this fvPatchField.
Top level data entry class for use in dictionaries. Provides a mechanism to specify a variable as a c...
virtual void write(Ostream &) const
Write.
Bound the given scalar field if it has gone unbounded.
OBJstream os(runTime.globalPath()/outputName)
scalar Cmu() const noexcept
Return the object: Cmu.
void writeLocalEntries(Ostream &) const
Write local wall function variables.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
const std::string patch
OpenFOAM patch number as a std::string.
wallFunctionCoefficients wallCoeffs_
Wall-function coefficients.
List< label > labelList
A List of labels.
scalar kappa() const noexcept
Return the object: kappa.
A class for managing temporary objects.
static const nutWallFunctionFvPatchScalarField & nutw(const turbulenceModel &turbModel, const label patchi)
Return the nut patchField for the given wall patch.
atmNutkWallFunctionFvPatchScalarField(const fvPatch &, const DimensionedField< scalar, volMesh > &)
Construct from patch and internal field.
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)
virtual void autoMap(const fvPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.