72 for (
const auto& z : z0)
77 <<
"z0 field can only contain positive values. " 78 <<
"Please check input field z0." 91 const scalar
yPlus = Cmu25*
y[facei]*
sqrt(
k[celli])/nuw[facei];
93 const scalar w = cornerWeights[facei];
97 w*Cmu75*
pow(
k[celli], 1.5)/(
kappa*(
y[facei] + z0[facei]));
101 *(nutw[facei] + nuw[facei])
103 *Cmu25*
sqrt(
k[celli])
104 /(
kappa*(
y[facei] + z0[facei]));
108 epsilonc = w*2.0*
k[celli]*nuw[facei]/
sqr(
y[facei] + z0[facei]);
131 wallCoeffs_.writeEntries(
os);
159 z0_(ptf.z0_.clone(
p.
patch()))
224 refCast<const atmEpsilonWallFunctionFvPatchScalarField>(ptf);
227 z0_->rmap(atmpsf.z0_(), addr);
238 writeLocalEntries(
os);
250 atmEpsilonWallFunctionFvPatchScalarField
const objectRegistry & db() const
The associated objectRegistry.
scalar yPlusLam() const noexcept
Return the object: yPlusLam.
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.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
atmEpsilonWallFunctionFvPatchScalarField(const fvPatch &, const DimensionedField< scalar, volMesh > &)
Construct from patch and internal field.
virtual tmp< volScalarField > nu() const =0
Return the laminar viscosity.
virtual void write(Ostream &) const
Write.
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.
virtual void autoMap(const fvPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.
label k
Boltzmann constant.
Smooth ATC in cells next to a set of patches supplied by type.
virtual void write(Ostream &) const
Write.
const dimensionedScalar epsilon0
Electric constant: default SI units: [F/m].
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.
fvPatchField< scalar > fvPatchScalarField
virtual void calculate(const turbulenceModel &turbulence, const List< scalar > &cornerWeights, const fvPatch &patch, scalarField &G, scalarField &epsilon)
Calculate the epsilon and G.
const Time & time() const noexcept
Return time registry.
virtual const labelUList & faceCells() const
Return faceCells.
A FieldMapper for finite-volume patch fields.
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.
This boundary condition provides wall functions for the turbulent kinetic energy dissipation rate (i...
Ostream & writeEntryIfDifferent(const word &key, const T &value1, const T &value2)
Write a keyword/value entry only when the two values differ.
Top level data entry class for use in dictionaries. Provides a mechanism to specify a variable as a c...
void writeLocalEntries(Ostream &) const
Write local wall function variables.
OBJstream os(runTime.globalPath()/outputName)
scalar Cmu() const noexcept
Return the object: Cmu.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
This boundary condition provides a wall constraint on the turbulent kinetic energy dissipation rate (...
autoPtr< PatchFunction1< scalar > > z0_
Surface roughness length field [m].
scalar timeOutputValue() const
Return the current user-time value. (ie, after applying any timeToUserTime() conversion) ...
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.
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
The index of this patch in the boundary mesh.
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)
const bool lowReCorrection_
Apply low-Re correction term (default = no)
virtual void rmap(const fvPatchScalarField &, const labelList &)
Reverse map the given fvPatchField onto this fvPatchField.
const Boundary & boundaryField() const noexcept
Return const-reference to the boundary field.