30 #include "twoPhaseSystem.H" 39 JohnsonJacksonParticleSlipFvPatchVectorField
53 partialSlipFvPatchVectorField(
p, iF),
54 specularityCoefficient_(
"specularityCoefficient",
dimless,
Zero)
67 partialSlipFvPatchVectorField(ptf,
p, iF, mapper),
68 specularityCoefficient_(ptf.specularityCoefficient_)
80 partialSlipFvPatchVectorField(
p, iF),
81 specularityCoefficient_(
"specularityCoefficient",
dimless,
dict)
85 (specularityCoefficient_.
value() < 0)
86 || (specularityCoefficient_.
value() > 1)
90 <<
"The specularity coefficient has to be between 0 and 1" 107 partialSlipFvPatchVectorField(ptf),
108 specularityCoefficient_(ptf.specularityCoefficient_)
119 partialSlipFvPatchVectorField(ptf, iF),
120 specularityCoefficient_(ptf.specularityCoefficient_)
131 partialSlipFvPatchVectorField::autoMap(m);
141 partialSlipFvPatchVectorField::rmap(ptf, addr);
153 const twoPhaseSystem&
fluid = db().lookupObject<twoPhaseSystem>
158 const phaseModel& phased
168 patch().lookupPatchField<volScalarField, scalar>
170 phased.volScalarField::name()
176 patch().lookupPatchField<volScalarField, scalar>
184 patch().lookupPatchField<volScalarField, scalar>
192 patch().lookupPatchField<volScalarField, scalar>
202 db().foundObject<volScalarField>(ThetaName)
203 ?
patch().lookupPatchField<volScalarField, scalar>(ThetaName)
213 .lookupObject<IOdictionary>
218 .subDict(
"kineticTheoryCoeffs")
227 *specularityCoefficient_.
value()
232 this->valueFraction() =
c/(
c +
patch().deltaCoeffs());
234 partialSlipFvPatchVectorField::updateCoeffs();
244 os.
writeEntry(
"specularityCoefficient", specularityCoefficient_);
245 writeEntry(
"value",
os);
const phaseModel & phase1() const
Return phase model 1.
const Type & value() const noexcept
Return const reference to value.
fvPatchField< vector > fvPatchVectorField
JohnsonJacksonParticleSlipFvPatchVectorField(const fvPatch &, const DimensionedField< vector, volMesh > &)
Construct from patch and internal field.
virtual void write(Ostream &) const
Write.
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.
virtual void updateCoeffs()
Update the coefficients.
friend Ostream & operator(Ostream &, const fvPatchField< Type > &)
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
virtual void rmap(const fvPatchVectorField &, const labelList &)
Reverse map the given fvPatchField onto this fvPatchField.
dimensionedScalar sqrt(const dimensionedScalar &ds)
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
Abstract base class with a fat-interface to all derived classes covering all possible ways in which t...
Ostream & writeEntry(const keyType &key, const T &value)
Write a keyword/value entry.
const dimensionSet dimless
Dimensionless.
virtual void write(Ostream &) const
Write.
Macros for easy insertion into run-time selection tables.
const phaseModel & phase2() const
Return phase model 2.
static word groupName(StringType base, const word &group)
Create dot-delimited name.group string.
dimensionedScalar alphaMax("alphaMax", dimless/dimTime, laminarTransport)
fvPatchField< scalar > fvPatchScalarField
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
virtual void autoMap(const fvPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.
A FieldMapper for finite-volume patch fields.
constexpr scalar pi(M_PI)
errorManip< error > abort(error &err)
An Ostream is an abstract base class for all output systems (streams, files, token lists...
OBJstream os(runTime.globalPath()/outputName)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
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.
const word & name() const
Field< vector > vectorField
Specialisation of Field<T> for vector.
Partial slip boundary condition for the particulate velocity.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)
static constexpr const zero Zero
Global zero (0)