44 namespace functionObjects
62 if (!foundObject<volScalarField>(fieldName_))
76 store(fieldName_, tfldPtr);
78 if (phaseName_ !=
"none")
84 return lookupObjectRef<volScalarField>(fieldName_);
94 const word Dname(
"D" +
s.name());
103 mesh_.time().timeName(),
113 if (nutName_ !=
"none")
124 findObject<incompressible::turbulenceModel>
134 alphaD_ *
turb->nu() + alphaDt_ *
turb->nut()
142 findObject<compressible::turbulenceModel>
152 alphaD_ *
turb->mu() + alphaDt_ *
turb->mut()
163 mesh_.time().timeName(),
176 Foam::functionObjects::scalarTransport::scalarTransport
184 fieldName_(
dict.getOrDefault<
word>(
"field",
"s")),
185 phiName_(
dict.getOrDefault<
word>(
"phi",
"phi")),
186 rhoName_(
dict.getOrDefault<
word>(
"rho",
"rho")),
187 nutName_(
dict.getOrDefault<
word>(
"nut",
"none")),
188 phaseName_(
dict.getOrDefault<
word>(
"phase",
"none")),
189 phasePhiCompressedName_
191 dict.getOrDefault<
word>(
"phasePhiCompressed",
"alphaPhiUn")
196 resetOnStartUp_(false),
197 schemesField_(
"unknown-schemesField"),
199 bounded01_(
dict.getOrDefault(
"bounded01", true))
226 dict.readIfPresent(
"phi", phiName_);
227 dict.readIfPresent(
"rho", rhoName_);
228 dict.readIfPresent(
"nut", nutName_);
229 dict.readIfPresent(
"phase", phaseName_);
230 dict.readIfPresent(
"bounded01", bounded01_);
232 schemesField_ =
dict.getOrDefault(
"schemesField", fieldName_);
233 constantD_ =
dict.readIfPresent(
"D", D_);
234 alphaD_ =
dict.getOrDefault<scalar>(
"alphaD", 1);
235 alphaDt_ =
dict.getOrDefault<scalar>(
"alphaDt", 1);
237 dict.readIfPresent(
"nCorr", nCorr_);
238 dict.readIfPresent(
"resetOnStartUp", resetOnStartUp_);
240 if (
dict.found(
"fvOptions"))
242 fvOptions_.reset(
dict.subDict(
"fvOptions"));
261 word divScheme(
"div(phi," + schemesField_ +
")");
262 word laplacianScheme(
"laplacian(" +
D.name() +
"," + schemesField_ +
")");
265 scalar relaxCoeff = 0.0;
266 if (mesh_.relaxEquation(schemesField_))
268 relaxCoeff = mesh_.equationRelaxationFactor(schemesField_);
272 if (phaseName_ !=
"none")
283 D.dimensions().reset(limitedPhiAlpha.dimensions()/
dimLength);
286 tmp<surfaceScalarField> tTPhiUD;
287 for (label i = 0; i <= nCorr_; i++)
292 +
fvm::div(limitedPhiAlpha,
s, divScheme)
298 sEqn.relax(relaxCoeff);
299 fvOptions_.constrain(sEqn);
300 sEqn.solve(mesh_.solverDict(schemesField_));
302 tTPhiUD = sEqn.flux();
322 for (label i = 0; i <= nCorr_; i++)
334 sEqn.relax(relaxCoeff);
336 fvOptions_.constrain(sEqn);
338 sEqn.solve(mesh_.solverDict(schemesField_));
343 for (label i = 0; i <= nCorr_; i++)
354 sEqn.relax(relaxCoeff);
356 fvOptions_.constrain(sEqn);
358 sEqn.solve(mesh_.solverDict(schemesField_));
364 <<
"Incompatible dimensions for phi: " <<
phi.dimensions() <<
nl 378 <<
tab << fieldName_ <<
nl
defineTypeNameAndDebug(ObukhovLength, 0)
fvMatrix< scalar > fvScalarMatrix
virtual ~scalarTransport()
Destructor.
errorManipArg< error, int > exit(error &err, const int errNo=1)
error FatalError
Error stream (stdout output on all processes), with additional 'FOAM FATAL ERROR' header text and sta...
A list of keyword definitions, which are a keyword followed by a number of values (eg...
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
constexpr char nl
The newline '\n' character (0x0a)
compressible::turbulenceModel & turb
Ostream & endl(Ostream &os)
Add newline and flush stream.
Calculate the matrix for the laplacian of the field.
constexpr char tab
The tab '\t' character(0x09)
Ignore writing from objectRegistry::writeObject()
const dimensionSet dimless
Dimensionless.
const Time & time() const
Return the top-level database.
Class to control time during OpenFOAM simulations that is also the top-level objectRegistry.
Macros for easy insertion into run-time selection tables.
void setFluxRequired(const word &name) const
Get flux-required for given name, or default.
virtual bool execute()
Calculate the scalarTransport.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
dimensionedScalar pos(const dimensionedScalar &ds)
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
const dimensionSet dimVolume(pow3(dimLength))
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
static const word propertiesName
Default name of the turbulence properties dictionary.
A class for handling words, derived from Foam::string.
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
bool store(word &fieldName, const tmp< ObjectType > &tfield, bool cacheable=false)
Store the field in the (sub) objectRegistry under the given name.
Calculate the matrix for the first temporal derivative.
virtual bool read(const dictionary &)
Read the scalarTransport data.
static tmp< T > New(Args &&... args)
Construct tmp with forwarding arguments.
static word timeName(const scalar t, const int precision=precision_)
Return time name of given scalar time formatted with the given precision.
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
addToRunTimeSelectionTable(functionObject, ObukhovLength, dictionary)
Calculate the matrix for the divergence of the given field and flux.
virtual bool write()
Do nothing.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
virtual bool read(const dictionary &dict)
Read optional controls.
Specialization of Foam::functionObject for an Foam::fvMesh, providing a reference to the Foam::fvMesh...
const dimensionedScalar & D
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
A class for managing temporary objects.
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
void explicitSolve(const RdeltaTType &rDeltaT, const RhoType &rho, volScalarField &psi, const surfaceScalarField &phiPsi, const SpType &Sp, const SuType &Su)
const fvMesh & mesh_
Reference to the fvMesh.
CMULES: Multidimensional universal limiter for explicit corrected implicit solution.
Calculate the finiteVolume matrix for implicit and explicit sources.
static constexpr const zero Zero
Global zero (0)