40 Foam::pressurePIDControlInletVelocityFvPatchVectorField::facePressure()
const 42 const word pfName(pName_ +
"f");
75 upstreamName_(
word::null),
76 downstreamName_(
word::null),
105 upstreamName_(ptf.upstreamName_),
106 downstreamName_(ptf.downstreamName_),
107 deltaP_(ptf.deltaP_),
108 shapeFactor_(ptf.shapeFactor_),
110 phiName_(ptf.phiName_),
111 rhoName_(ptf.rhoName_),
117 errorIntegral_(ptf.errorIntegral_),
119 oldError_(ptf.oldError_),
120 oldErrorIntegral_(ptf.oldErrorIntegral_),
121 timeIndex_(ptf.timeIndex_)
136 deltaP_(
dict.
get<scalar>(
"deltaP")),
137 shapeFactor_(
dict.getOrDefault<scalar>(
"shapeFactor", 0)),
138 pName_(
dict.getOrDefault<
word>(
"p",
"p")),
139 phiName_(
dict.getOrDefault<
word>(
"phi",
"phi")),
140 rhoName_(
dict.getOrDefault<
word>(
"rho",
"none")),
145 error_(
dict.getOrDefault<scalar>(
"error", 0)),
146 errorIntegral_(
dict.getOrDefault<scalar>(
"errorIntegral", 0)),
149 oldErrorIntegral_(0),
161 upstreamName_(ptf.upstreamName_),
162 downstreamName_(ptf.downstreamName_),
163 deltaP_(ptf.deltaP_),
164 shapeFactor_(ptf.shapeFactor_),
166 phiName_(ptf.phiName_),
167 rhoName_(ptf.rhoName_),
173 errorIntegral_(ptf.errorIntegral_),
175 oldError_(ptf.oldError_),
176 oldErrorIntegral_(ptf.oldErrorIntegral_),
177 timeIndex_(ptf.timeIndex_)
189 upstreamName_(ptf.upstreamName_),
190 downstreamName_(ptf.downstreamName_),
191 deltaP_(ptf.deltaP_),
192 shapeFactor_(ptf.shapeFactor_),
194 phiName_(ptf.phiName_),
195 rhoName_(ptf.rhoName_),
201 errorIntegral_(ptf.errorIntegral_),
203 oldError_(ptf.oldError_),
204 oldErrorIntegral_(ptf.oldErrorIntegral_),
205 timeIndex_(ptf.timeIndex_)
222 const scalar deltaT(db().time().deltaTValue());
227 db().lookupObject<surfaceScalarField>(phiName_)
231 if (timeIndex_ != db().time().
timeIndex())
233 timeIndex_ = db().time().timeIndex();
236 oldErrorIntegral_ = errorIntegral_;
247 const fvPatchField<scalar>& rhoField =
255 <<
"The dimensions of the field " << phiName_
256 <<
"are not recognised. The dimensions are " <<
phi.dimensions()
258 <<
" for an incompressible case, or " 264 const scalar patchA =
gSum(
patch().magSf());
270 faceZoneAverage(upstreamName_,
mesh.
Cf(), Aa, xa);
271 faceZoneAverage(downstreamName_,
mesh.
Cf(), Ab, xb);
272 const scalar
L =
mag(xa - xb);
273 const scalar LbyALinear =
L/(Aa - Ab)*
log(Aa/Ab);
274 const scalar LbyAStep =
L/2*(1/Aa + 1/Ab);
275 const scalar LbyA = (1 - shapeFactor_)*LbyALinear + shapeFactor_*LbyAStep;
281 scalar deltaP = deltaP_;
285 faceZoneAverage(upstreamName_, facePressure(), Aa, pa);
286 faceZoneAverage(downstreamName_, facePressure(), Ab, pb);
292 <<
"The pressure field name, 'p' is \"" << pName_ <<
"\", " 293 <<
"but a field of that name was not found. The inlet velocity " 294 <<
"will be set to an analytical value calculated from the " 295 <<
"specified pressure drop. No PID control will be done and " 296 <<
"transient effects will be ignored. This behaviour is designed " 297 <<
"to be appropriate for potentialFoam solutions. If you are " 298 <<
"getting this warning from another solver, you have probably " 299 <<
"specified an incorrect pressure name." 304 scalar QTarget, QMeasured;
305 const scalar a = (1/
sqr(Ab) - 1/
sqr(Aa))/(2*
rho);
308 const scalar
b = LbyA/deltaT;
309 const scalar
c = - LbyA/deltaT*oldQ_ ;
310 QTarget = (-
b +
sqrt(
sqr(
b) - 4*a*(
c - deltaP_)))/(2*a);
311 QMeasured = (-
b +
sqrt(
sqr(
b) - 4*a*(
c - deltaP)))/(2*a);
315 QTarget =
sqrt(deltaP_/a);
316 QMeasured =
sqrt(deltaP/a);
320 error_ = QTarget - QMeasured;
321 errorIntegral_ = oldErrorIntegral_ + 0.5*(error_ + oldError_);
322 const scalar errorDifferential = oldError_ - error_;
332 + D_*errorDifferential
340 const scalar error = deltaP/deltaP_ - 1;
341 const scalar newQ = -
gSum(*
this &
patch().Sf());
342 Info<<
"pressurePIDControlInletVelocityFvPatchField " <<
patch().name()
347 <<
" (" <<
mag(error)*100 <<
"% " 348 << (error < 0 ?
"below" :
"above") <<
" the target)" <<
endl;
374 writeEntry(
"value",
os);
385 pressurePIDControlInletVelocityFvPatchVectorField
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
List of values generated by applying the access operation to each list item.
This boundary condition supplies a fixed value constraint, and is the base class for a number of othe...
virtual void write(Ostream &) const
Write.
fvPatchField< vector > fvPatchVectorField
dimensionedScalar log(const dimensionedScalar &ds)
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...
const vector L(dict.get< vector >("L"))
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
const surfaceVectorField & Cf() const
Return face centres as surfaceVectorField.
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > linearInterpolate(const GeometricField< Type, fvPatchField, volMesh > &vf)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
Ostream & writeEntry(const keyType &key, const T &value)
Write a keyword/value entry.
virtual void write(Ostream &) const
Write.
Lookup type of boundary radiation properties.
pressurePIDControlInletVelocityFvPatchVectorField(const fvPatch &, const DimensionedField< vector, volMesh > &)
Construct from patch and internal field.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Type gSum(const FieldField< Field, Type > &f)
const dimensionedScalar b
Wien displacement law constant: default SI units: [m.K].
A class for handling words, derived from Foam::string.
A FieldMapper for finite-volume patch fields.
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...
int debug
Static debugging option.
OBJstream os(runTime.globalPath()/outputName)
virtual void updateCoeffs()
Update the coefficients associated with the patch field.
const dimensionSet dimDensity
#define WarningInFunction
Report a warning using Foam::Warning.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
virtual void updateCoeffs()
Update the coefficients associated with the patch field.
const dimensionedScalar c
Speed of light in a vacuum.
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.
messageStream Info
Information stream (stdout output on master, null elsewhere)
bool steady() const noexcept
True if default ddtScheme is steady-state.
This boundary condition tries to generate an inlet velocity that maintains a specified pressure drop ...
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)
const dimensionSet dimArea(sqr(dimLength))
const dimensionSet dimVelocity