70 for (
const scalar z : z0)
75 <<
"z0 field can only contain positive values. " 76 <<
"Please check input field z0." 89 const scalar w = cornerWeights[facei];
92 w*
sqrt(
k[celli])/(Cmu25*
kappa*(
y[facei] + z0[facei]));
96 *(nutw[facei] + nuw[facei])
99 /(
kappa*(
y[facei] + z0[facei]));
114 wallCoeffs_.writeEntries(
os);
142 z0_(ptf.z0_.clone(
p.
patch()))
207 refCast<const atmOmegaWallFunctionFvPatchScalarField>(ptf);
211 z0_->rmap(atmpsf.z0_(), addr);
222 writeLocalEntries(
os);
223 writeEntry(
"value",
os);
234 atmOmegaWallFunctionFvPatchScalarField
const objectRegistry & db() const
The associated objectRegistry.
virtual void autoMap(const fvPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.
virtual tmp< Field< Type > > snGrad() const
Return patch-normal gradient.
const dimensionedScalar G0
Conductance quantum: default SI units: [S].
errorManipArg< error, int > exit(error &err, const int errNo=1)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const volVectorField & U() const
Access function to velocity field.
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...
const fvPatch & patch() const noexcept
Return the patch.
virtual void rmap(const fvPatchScalarField &, const labelList &)
Reverse map the given fvPatchField onto this fvPatchField.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
atmOmegaWallFunctionFvPatchScalarField(const fvPatch &, const DimensionedField< scalar, volMesh > &)
Construct from patch and internal field.
This boundary condition provides a wall constraint on the specific dissipation rate (i...
virtual tmp< volScalarField > nu() const =0
Return the laminar viscosity.
dimensionedScalar sqrt(const dimensionedScalar &ds)
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
dimensionedScalar pow025(const dimensionedScalar &ds)
Abstract base class with a fat-interface to all derived classes covering all possible ways in which t...
label k
Boltzmann constant.
Smooth ATC in cells next to a set of patches supplied by type.
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
virtual void write(Ostream &) const
Write.
void writeLocalEntries(Ostream &) const
Write local wall function variables.
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.
virtual tmp< volScalarField > k() const =0
Return the turbulence kinetic energy.
#define forAll(list, i)
Loop across all elements in list.
virtual void autoMap(const fvPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.
fvPatchField< scalar > fvPatchScalarField
virtual void write(Ostream &) const
Write.
const Time & time() const noexcept
Return time registry.
virtual const labelUList & faceCells() const
Return faceCells.
A FieldMapper for finite-volume patch fields.
virtual void calculate(const turbulenceModel &turbulence, const List< scalar > &cornerWeights, const fvPatch &patch, scalarField &G, scalarField &omega)
Calculate the omega and G.
virtual void rmap(const fvPatchField< scalar > &, const labelList &)
Reverse map the given fvPatchField onto this fvPatchField.
const nearWallDist & y() const
Return the near wall distances.
Top level data entry class for use in dictionaries. Provides a mechanism to specify a variable as a c...
OBJstream os(runTime.globalPath()/outputName)
wallFunctionCoefficients wallCoeffs_
Wall-function coefficients.
scalar Cmu() const noexcept
Return the object: Cmu.
scalar timeOutputValue() const
Return current time value.
This boundary condition provides a wall function for the specific dissipation rate (i...
autoPtr< PatchFunction1< scalar > > z0_
Surface roughness length field [m].
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.
virtual tmp< volScalarField > nut() const =0
Return the turbulence viscosity.
List< label > labelList
A List of labels.
scalar kappa() const noexcept
Return the object: kappa.
A class for managing temporary objects.
label index() const noexcept
Return the index of this patch in the fvBoundaryMesh.
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)
const Boundary & boundaryField() const noexcept
Return const-reference to the boundary field.