42 for (
int i=0; i<10; i++)
64 const scalar Cmu25 =
pow025(Cmu_);
65 const scalar Cmu75 =
pow(Cmu_, 0.75);
91 scalar
yPlus = Cmu25*
sqrt(k[celli])*y[facei]/nuw[facei];
93 scalar w = cornerWeights[facei];
95 if (yPlus > yPlusLam_)
97 epsilon0[celli] += w*Cmu75*
pow(k[celli], 1.5)/(kappa_*y[facei]);
101 *(nutw[facei] + nuw[facei])
103 *Cmu25*
sqrt(k[celli])
108 epsilon0[celli] += w*2.0*k[celli]*nuw[facei]/
sqr(y[facei]);
109 G0[celli] += G[celli];
125 yPlusLam_(yPlusLam(kappa_, E_))
152 yPlusLam_(yPlusLam(kappa_, E_))
word GName() const
Helper function to return the name of the turbulence G field.
scalar yPlusLam(const scalar kappa, const scalar E)
Calculate the Y+ at the edge of the laminar sublayer.
#define forAll(list, i)
Loop across all elements in list.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
dimensionedScalar log(const dimensionedScalar &ds)
A list of keyword definitions, which are a keyword followed by any number of values (e...
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
const dimensionedScalar G
Newtonian constant of gravitation.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
virtual tmp< volScalarField > nu() const =0
Return the laminar viscosity.
const volVectorField & U() const
Access function to velocity field.
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.
virtual void calculate(const turbulenceModel &turbulence, const List< scalar > &cornerWeights, const fvPatch &patch, scalarField &G, scalarField &epsilon)
Calculate the epsilon and G.
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.
makePatchTypeField(fvPatchVectorField, SRFFreestreamVelocityFvPatchVectorField)
const Boundary & boundaryField() const
Return const-reference to the boundary field.
Foam::fvPatchFieldMapper.
This boundary condition provides a turbulence dissipation wall function condition for high Reynolds n...
virtual tmp< Field< Type > > snGrad() const
Return patch-normal gradient.
virtual const labelUList & faceCells() const
Return faceCells.
epsilonLowReWallFunctionFvPatchScalarField(const fvPatch &, const DimensionedField< scalar, volMesh > &)
Construct from patch and internal field.
label index() const
Return the index of this patch in the fvBoundaryMesh.
scalar yPlusLam_
Y+ at the edge of the laminar sublayer.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
This boundary condition provides a turbulence dissipation wall function condition for low- and high-R...
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
dimensioned< scalar > mag(const dimensioned< Type > &)
virtual tmp< volScalarField > nut() const =0
Return the turbulence viscosity.
A class for managing temporary objects.
const nearWallDist & y() const
Return the near wall distances.