Go to the documentation of this file.
46 void Foam::fv::atmBuoyancyTurbSource::calcB()
61 Foam::fv::atmBuoyancyTurbSource::calcC3
84 neg0(Rig)*(1.0 - (1.0 + (C2_ - 1.0)/(C2_ - C1_))*LbyLmax)
85 +
pos(Rig)*(1.0 - LbyLmax);
88 return (C1_ - C2_)*alphaB + 1.0;
93 Foam::fv::atmBuoyancyTurbSource::calcC3
119 +
pos(Rig)*(1.0 - LbyLmax);
130 const word& sourceName,
131 const word& modelType,
138 rhoName_(coeffs_.getOrDefault<
word>(
"rho",
"rho")),
144 coeffs_.getCheckOrDefault<scalar>
148 [&](
const scalar Lmax){ return Lmax > SMALL; }
157 coeffs_.getCheckOrDefault<scalar>
161 [&](
const scalar
x){ return x > SMALL; }
179 mesh.time().timeName(),
188 const auto* turbPtr =
197 <<
"Unable to find a turbulence model."
201 fieldNames_.resize(2);
206 if (!tepsilon.
isTmp())
208 fieldNames_[0] = tepsilon().name();
210 const dictionary& turbDict = turbPtr->coeffDict();
212 Cmu_.
read(
"Cmu", turbDict);
213 C1_.read(
"C1", turbDict);
214 C2_.read(
"C2", turbDict);
216 else if (!tomega.
isTmp())
219 fieldNames_[0] = tomega().name();
221 const dictionary& turbDict = turbPtr->coeffDict();
223 Cmu_.
read(
"betaStar", turbDict);
228 <<
"Unable to find neither epsilon nor omega field." <<
nl
229 <<
"atmBuoyancyTurbSource needs either epsilon or omega field."
233 fieldNames_[1] = turbPtr->k()().
name();
237 Log <<
" Applying atmBuoyancyTurbSource to: "
238 << fieldNames_[0] <<
" and " << fieldNames_[1]
253 atmBuoyancyTurbSourceK
267 atmBuoyancyTurbSourceEpsilon
277 atmBuoyancyTurbSourceOmega
324 atmBuoyancyTurbSourceK(
alpha,
rho, eqn, fieldi);
332 atmBuoyancyTurbSourceEpsilon(
alpha,
rho, eqn, fieldi);
336 atmBuoyancyTurbSourceOmega(
alpha,
rho, eqn, fieldi);
Defines the attributes of an object for which implicit objectRegistry management is supported,...
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
A class for handling words, derived from Foam::string.
Intermediate abstract class for handling cell-set options for the derived fvOptions.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
const dimensionedScalar G
Newtonian constant of gravitation.
A class for managing temporary objects.
static constexpr const zero Zero
Global zero (0)
void resetApplied()
Resize/reset applied flag list for all fieldNames_ entries.
virtual void addSup(fvMatrix< scalar > &eqn, const label fieldi)
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
dimensionedScalar neg0(const dimensionedScalar &ds)
static const word propertiesName
Default name of the turbulence properties dictionary.
A class representing the concept of a GeometricField of 1 used to avoid unnecessary manipulations for...
bool isTmp() const noexcept
Identical to is_pointer()
Ostream & endl(Ostream &os)
Add newline and flush stream.
const fvMesh & mesh_
Reference to the mesh database.
dimensionedScalar pow025(const dimensionedScalar &ds)
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
dimensionedScalar pow3(const dimensionedScalar &ds)
DimensionedField< scalar, volMesh > Internal
Type of the internal field from which this GeometricField is derived.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
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,...
Macros for easy insertion into run-time selection tables.
Mesh data needed to do the Finite Volume discretisation.
errorManip< error > abort(error &err)
bool read(Istream &is)
Read dictionary from Istream. Discards the header.
Type & lookupObjectRef(const word &name, const bool recursive=false) const
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
label k
Boltzmann constant.
A special matrix type and solver, designed for finite volume solutions of scalar equations....
const dimensionSet dimTemperature(0, 0, 0, 1, 0, 0, 0)
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
defineTypeNameAndDebug(atmAmbientTurbSource, 0)
addToRunTimeSelectionTable(option, atmAmbientTurbSource, dictionary)
const dimensionSet & dimensions() const
Return const reference to dimensions.
atmBuoyancyTurbSource(const word &sourceName, const word &modelType, const dictionary &dict, const fvMesh &mesh)
Construct from explicit source name and mesh.
static const gravity & New(const Time &runTime)
Construct on Time.
const dimensionSet dimless
Dimensionless.
dimensionedScalar pos(const dimensionedScalar &ds)