Go to the documentation of this file.
46 const label patchi =
patch.index();
49 nutWallFunctionFvPatchScalarField::nutw(turbModel, patchi);
65 const scalar t = db().time().timeOutputValue();
69 for (
const scalar z : z0)
74 <<
"z0 field can only contain positive values. "
75 <<
"Please check input field z0."
88 const scalar w = cornerWeights[facei];
91 w*
sqrt(
k[celli])/(Cmu25*nutw.
kappa()*(
y[facei] + z0[facei]));
95 *(nutw[facei] + nuw[facei])
98 /(nutw.
kappa()*(
y[facei] + z0[facei]));
127 z0_(ptf.
z0_.clone(
p.patch()))
151 z0_(owfpsf.
z0_.clone(this->patch().patch()))
163 z0_(owfpsf.
z0_.clone(this->patch().patch()))
174 omegaWallFunctionFvPatchScalarField::autoMap(m);
185 omegaWallFunctionFvPatchScalarField::rmap(ptf, addr);
188 refCast<const atmOmegaWallFunctionFvPatchScalarField>(ptf);
190 z0_->rmap(atmpsf.
z0_(), addr);
virtual tmp< Field< Type > > snGrad() const
Return patch-normal gradient.
virtual void write(Ostream &) const
Write.
atmOmegaWallFunctionFvPatchScalarField(const fvPatch &, const DimensionedField< scalar, volMesh > &)
Construct from patch and internal field.
virtual tmp< volScalarField > nu() const =0
Return the laminar viscosity.
This boundary condition provides a wall constraint on the specific dissipation rate,...
scalar Cmu() const
Return Cmu.
A class for managing temporary objects.
This boundary condition provides a wall constraint on the specific dissipation rate (i....
const nearWallDist & y() const
Return the near wall distances.
virtual void rmap(const fvPatchScalarField &, const labelList &)
Reverse map the given fvPatchField onto this fvPatchField.
#define forAll(list, i)
Loop across all elements in list.
dimensionedScalar pow025(const dimensionedScalar &ds)
autoPtr< PatchFunction1< scalar > > z0_
Surface roughness length field [m].
virtual void calculate(const turbulenceModel &turbulence, const List< scalar > &cornerWeights, const fvPatch &patch, scalarField &G, scalarField &omega)
Calculate the omega and G.
const dimensionedScalar G0
Conductance quantum: default SI units: [S].
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
Abstract base class for turbulence models (RAS, LES and laminar).
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
OBJstream os(runTime.globalPath()/outputName)
Macros for easy insertion into run-time selection tables.
Top level data entry class for use in dictionaries. Provides a mechanism to specify a variable as a c...
virtual tmp< volScalarField > k() const =0
Return the turbulence kinetic energy.
virtual void autoMap(const fvPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.
errorManipArg< error, int > exit(error &err, const int errNo=1)
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
const std::string patch
OpenFOAM patch number as a std::string.
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
scalar kappa() const
Return kappa.
label k
Boltzmann constant.
void write(vtk::formatter &fmt, const Type &val, const label n=1)
Component-wise write of a value (N times)
Foam::fvPatchFieldMapper.
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)
The class nutWallFunction is a base class that parents the derived boundary conditions which provide ...
Smooth ATC in cells next to a set of patches supplied by type.
const Boundary & boundaryField() const
Return const-reference to the boundary field.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
const volVectorField & U() const
Access function to velocity field.