53 if (isA<epsilonWallFunctionFvPatchScalarField>(bf[
patchi]))
100 if (isA<epsilonWallFunctionFvPatchScalarField>(bf[
patchi]))
102 epsilonPatches.append(patchi);
107 weights[faceCells[i]]++;
133 const volScalarField::Boundary& bf = epsilon.
boundaryField();
136 refCast<const epsilonWallFunctionFvPatchScalarField>(bf[
patchi]);
202 const scalar Cmu75 =
pow(nutw.
Cmu(), 0.75);
209 const scalar
yPlus = Cmu25*y[facei]*
sqrt(k[celli])/nuw[facei];
211 const scalar w = cornerWeights[facei];
216 w*Cmu75*
pow(k[celli], 1.5)/(nutw.
kappa()*y[facei]);
220 *(nutw[facei] + nuw[facei])
222 *Cmu25*
sqrt(k[celli])
223 /(nutw.
kappa()*y[facei]);
227 epsilon0[celli] += w*2.0*k[celli]*nuw[facei]/
sqr(y[facei]);
369 momentumTransportModel::typeName,
388 const_cast<FieldType&
> 399 G[celli] = G0[celli];
400 epsilon[celli] = epsilon0[celli];
422 momentumTransportModel::typeName,
441 const_cast<FieldType&
> 453 scalar w = weights[facei];
459 G[celli] = (1.0 - w)*G[celli] + w*G0[celli];
460 epsilon[celli] = (1.0 - w)*epsilon[celli] + w*epsilon0[celli];
461 epsilonf[facei] = epsilon[celli];
503 label nConstrainedCells = 0;
513 label celli = faceCells[facei];
515 constraintCells.
append(celli);
516 constraintEpsilon.append(epsilon[celli]);
523 <<
": number of constrained cells = " << nConstrainedCells
virtual tmp< volScalarField > nu() const =0
Return the laminar viscosity.
const char *const group
Group name for atomic constants.
bool changing() const
Is mesh changing (topology changing and/or moving)
#define forAll(list, i)
Loop across all elements in list.
bool empty() const
Return true if the UList is empty (ie, size() is zero)
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
virtual tmp< Field< Type > > snGrad() const
Return patch-normal gradient.
label index() const
Return the index of this patch in the fvBoundaryMesh.
virtual void calculate(const momentumTransportModel &turbulence, const List< scalar > &cornerWeights, const fvPatch &patch, scalarField &G, scalarField &epsilon)
Calculate the epsilon and G.
const dimensionedScalar & epsilon0
Electric constant: default SI units: [F/m].
A list of keyword definitions, which are a keyword followed by any number of values (e...
bool updated() const
Return true if the boundary condition has already been updated.
const Boundary & boundaryField() const
Return const-reference to the boundary field.
virtual void calculateTurbulenceFields(const momentumTransportModel &turbulence, scalarField &G0, scalarField &epsilon0)
Main driver to calculate the turbulence fields.
scalarField G_
Local copy of turbulence G field.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
void size(const label)
Override size to be inconsistent with allocated storage.
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
dimensionedScalar pow025(const dimensionedScalar &ds)
virtual void manipulateMatrix(fvMatrix< scalar > &matrix)
Manipulate matrix.
static scalar tolerance_
Tolerance used in weighted calculations.
Abstract base class with a fat-interface to all derived classes covering all possible ways in which t...
const nearWallDist & y() const
Return the near wall distances.
static word timeName(const scalar, const int precision=precision_)
Return time name of given scalar time.
label k
Boltzmann constant.
const Time & time() const
Return the top-level database.
bool manipulatedMatrix() const
Return true if the matrix has already been manipulated.
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
scalar Cmu() const
Return Cmu.
Macros for easy insertion into run-time selection tables.
virtual void operator==(const fvPatchField< Type > &)
static const nutWallFunctionFvPatchScalarField & nutw(const momentumTransportModel &turbModel, const label patchi)
Return the nut patchField for the given wall patch.
label master_
Master patch ID.
This boundary condition provides a turbulent kinematic viscosity condition when using wall functions...
virtual label & master()
Return non-const access to the master patch ID.
const dimensionedScalar & e
Elementary charge.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
virtual const labelUList & faceCells() const
Return faceCells.
static word groupName(Name name, const word &group)
void setValues(const labelUList &cells, const UList< Type > &values)
Set solution in given cells to the specified values.
Foam::fvPatchFieldMapper.
DynamicList< T, SizeInc, SizeMult, SizeDiv > & append(const T &)
Append an element at the end of the list.
scalarField epsilon_
Local copy of turbulence epsilon field.
A special matrix type and solver, designed for finite volume solutions of scalar equations. Face addressing is used to make all matrix assembly and solution loops vectorise.
virtual void manipulateMatrix(fvMatrix< Type > &matrix)
Manipulate matrix.
This boundary condition provides a turbulence dissipation wall constraint for low- and high-Reynolds ...
A 1D vector of objects of type <T>, where the size of the vector is known and can be used for subscri...
const fvPatch & patch() const
Return patch.
scalarField & G(bool init=false)
Return non-const access to the master's G field.
virtual void updateCoeffs()
Update the coefficients associated with the patch field.
virtual label size() const
Return size.
const dimensionedScalar & G0
Conductance quantum: default SI units: [S].
virtual epsilonWallFunctionFvPatchScalarField & epsilonPatch(const label patchi)
Helper function to return non-const access to an epsilon patch.
virtual tmp< Field< Type > > patchInternalField() const
Return internal field next to patch as patch field.
const Mesh & mesh() const
Return mesh.
virtual tmp< volScalarField > k() const =0
Return the turbulence kinetic energy.
Abstract base class for turbulence models (RAS, LES and laminar).
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
void setSize(const label)
Reset size of List.
const volVectorField & U() const
Access function to velocity field.
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
List< List< scalar > > cornerWeights_
List of averaging corner weights.
virtual void updateWeightedCoeffs(const scalarField &weights)
Update the coefficients associated with the patch field.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
prefixOSstream Pout(cout, "Pout")
virtual void updateCoeffs()
Update the coefficients associated with the patch field.
Mesh data needed to do the Finite Volume discretisation.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
const objectRegistry & db() const
Return local objectRegistry.
dimensioned< scalar > mag(const dimensioned< Type > &)
word GName() const
Helper function to return the name of the turbulence G field.
virtual const word & name() const
Return name.
virtual void setMaster()
Set the master patch - master is responsible for updating all.
bool initialised_
Initialised flag.
static scalar yPlusLam(const scalar kappa, const scalar E)
Calculate the Y+ at the edge of the laminar sublayer.
A class for managing temporary objects.
virtual void createAveragingWeights()
Create the averaging weights for cells which are bounded by.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
epsilonWallFunctionFvPatchScalarField(const fvPatch &, const DimensionedField< scalar, volMesh > &)
Construct from patch and internal field.
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)
label size() const
Return the number of elements in the UList.
scalar kappa() const
Return kappa.
const DimensionedField< Type, volMesh > & internalField() const
Return dimensioned internal field reference.
scalarField & epsilon(bool init=false)
Return non-const access to the master's epsilon field.