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" 104 partialSlipFvPatchVectorField(ptf),
105 specularityCoefficient_(ptf.specularityCoefficient_)
116 partialSlipFvPatchVectorField(ptf, iF),
117 specularityCoefficient_(ptf.specularityCoefficient_)
128 partialSlipFvPatchVectorField::autoMap(m);
138 partialSlipFvPatchVectorField::rmap(ptf, addr);
150 const twoPhaseSystem&
fluid = db().lookupObject<twoPhaseSystem>
155 const phaseModel& phased
165 patch().lookupPatchField<volScalarField>
167 phased.volScalarField::name()
173 patch().lookupPatchField<volScalarField>
181 patch().lookupPatchField<volScalarField>
189 patch().lookupPatchField<volScalarField>
199 db().foundObject<volScalarField>(ThetaName)
200 ?
patch().lookupPatchField<volScalarField>(ThetaName)
210 .lookupObject<IOdictionary>
215 .subDict(
"kineticTheoryCoeffs")
224 *specularityCoefficient_.
value()
229 this->valueFraction() =
c/(
c +
patch().deltaCoeffs());
231 partialSlipFvPatchVectorField::updateCoeffs();
241 os.
writeEntry(
"specularityCoefficient", specularityCoefficient_);
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.
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.
Ostream & writeEntry(const keyType &key, const T &value)
Write a keyword/value entry.
void writeValueEntry(Ostream &os) const
Write *this field as a "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
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)