49 extrapolateProfile_(false)
64 rhoInlet_(
dict.getOrDefault<scalar>(
"rhoInlet", -VGREAT)),
68 dict.getOrDefault<
Switch>(
"extrapolateProfile", false)
80 dict.readIfPresent(
"rho", rhoName_);
88 <<
"Please supply either 'volumetricFlowRate' or" 89 <<
" 'massFlowRate' (optional: with 'rho')" <<
nl 94 if (!this->readValueEntry(
dict))
104 const flowRateInletVelocityFvPatchVectorField& ptf,
106 const DimensionedField<vector, volMesh>& iF,
107 const fvPatchFieldMapper& mapper
110 fixedValueFvPatchField<
vector>(ptf,
p, iF, mapper),
111 flowRate_(ptf.flowRate_.clone()),
112 rhoName_(ptf.rhoName_),
113 rhoInlet_(ptf.rhoInlet_),
114 volumetric_(ptf.volumetric_),
115 extrapolateProfile_(ptf.extrapolateProfile_)
126 flowRate_(ptf.flowRate_.clone()),
127 rhoName_(ptf.rhoName_),
128 rhoInlet_(ptf.rhoInlet_),
129 volumetric_(ptf.volumetric_),
130 extrapolateProfile_(ptf.extrapolateProfile_)
142 flowRate_(ptf.flowRate_.clone()),
143 rhoName_(ptf.rhoName_),
144 rhoInlet_(ptf.rhoInlet_),
145 volumetric_(ptf.volumetric_),
146 extrapolateProfile_(ptf.extrapolateProfile_)
152 template<
class RhoType>
153 void Foam::flowRateInletVelocityFvPatchVectorField::updateValues
158 const scalar t = db().time().timeOutputValue();
162 if (extrapolateProfile_)
173 nUp =
min(nUp, scalar(0));
175 const scalar flowRate = flowRate_->value(t);
176 const scalar estimatedFlowRate = -
gSum(
rho*(this->
patch().magSf()*nUp));
178 if (estimatedFlowRate > 0.5*flowRate)
180 nUp *= (
mag(flowRate)/
mag(estimatedFlowRate));
184 nUp -= ((flowRate - estimatedFlowRate)/
gSum(
rho*
patch().magSf()));
195 const scalar avgU = -flowRate_->value(t)/
gSum(
rho*
patch().magSf());
208 if (volumetric_ || rhoName_ ==
"none")
215 if (db().foundObject<volScalarField>(rhoName_))
228 <<
"Did not find registered density field " << rhoName_
229 <<
" and no constant density 'rhoInlet' specified" 233 updateValues(rhoInlet_);
237 fixedValueFvPatchVectorField::updateCoeffs();
244 flowRate_->writeData(
os);
250 if (extrapolateProfile_)
252 os.
writeEntry(
"extrapolateProfile", extrapolateProfile_);
265 flowRateInletVelocityFvPatchVectorField
This boundary condition supplies a fixed value constraint, and is the base class for a number of othe...
fvPatchField< vector > fvPatchVectorField
errorManipArg< error, int > exit(error &err, const int errNo=1)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
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)
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
A simple wrapper around bool so that it can be read as a word: true/false, on/off, yes/no, any/none. Also accepts 0/1 as a string and shortcuts t/f, y/n.
Ostream & writeEntry(const keyType &key, const T &value)
Write a keyword/value entry.
virtual void write(Ostream &) const
Write.
void writeValueEntry(Ostream &os) const
Write *this field as a "value" entry.
virtual void write(Ostream &) const
Write.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Type gSum(const FieldField< Field, Type > &f)
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
flowRateInletVelocityFvPatchVectorField(const fvPatch &, const DimensionedField< vector, volMesh > &)
Construct from patch and internal field.
virtual void updateCoeffs()
Update the coefficients associated with the patch field.
string evaluate(label fieldWidth, const std::string &s, size_t pos=0, size_t len=std::string::npos)
String evaluation with specified (positive, non-zero) field width.
Velocity inlet boundary condition either correcting the extrapolated velocity or creating a uniform v...
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Ostream & writeEntryIfDifferent(const word &key, const T &value1, const T &value2)
Write a keyword/value entry only when the two values differ.
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
OBJstream os(runTime.globalPath()/outputName)
static autoPtr< Function1< Type > > NewIfPresent(const word &entryName, const dictionary &dict, const word &redirectType, const objectRegistry *obrPtr=nullptr)
An optional selector, with fallback redirection.
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
A simple container of IOobject preferences. Can also be used for general handling of read/no-read/rea...
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
const std::string patch
OpenFOAM patch number as a std::string.
Field< vector > vectorField
Specialisation of Field<T> for vector.
tmp< faMatrix< Type > > operator==(const faMatrix< Type > &, const faMatrix< Type > &)
"buffered" : (MPI_Bsend, MPI_Recv)
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)
IOerror FatalIOError
Error stream (stdout output on all processes), with additional 'FOAM FATAL IO ERROR' header text and ...