51 if (geometricD[iDir] == -1)
53 averageVol /= bounds.
span()[iDir];
58 scalar multMeanRadius =
dict.getOrDefault<scalar>(
"meanRadiusMult", 10);
59 Info<<
"Computed a mean radius of " << radius
60 <<
" and multiplying with " << multMeanRadius <<
endl;
61 return multMeanRadius*radius;
67 Foam::isotropic::isotropic
71 bool adjustWallThickness
81 dict_.getOrDefault<scalar>(
"radius", computeRadius(
dict))
86 if (adjustWallThickness)
88 const scalar wallThicknessMult =
89 dict.getOrDefault<scalar>(
"wallThicknessMult", 0.75);
91 "Adjusting wall thickness by " << wallThicknessMult <<
endl;
Base class for selecting the regulatisation radius.
A list of keyword definitions, which are a keyword followed by a number of values (eg...
An isotropic regularisationRadius (same in all spatial directions)
virtual void addRegularisationTerm(fvScalarMatrix &matrix, bool isTopoField) const
Add a Laplacian term with an isotropic diffusivity.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
A bounding box defined in terms of min/max extrema points.
dimensionedScalar radius_
Smoothing radius of the first regulatisation.
const DimensionedField< scalar, volMesh > & V() const
Return cell volumes.
label nGeometricD() const
Return the number of valid geometric dimensions in the mesh.
Macros for easy insertion into run-time selection tables.
#define forAll(list, i)
Loop across all elements in list.
const Vector< label > & geometricD() const
Return the vector of geometric directions in mesh.
const GeometricField< Type, fvPatchField, volMesh > & psi(const label i=0) const
Return psi.
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.
#define DebugInfo
Report an information message using Foam::Info.
defineTypeNameAndDebug(combustionModel, 0)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
Type gAverage(const FieldField< Field, Type > &f)
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Mesh data needed to do the Finite Volume discretisation.
vector span() const
The bounding box span (from minimum to maximum)
messageStream Info
Information stream (stdout output on master, null elsewhere)
const boundBox & bounds() const noexcept
Return mesh bounding box.
scalar computeRadius(const dictionary &dict)
Compute smoothing radius, if not directly given.
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)