30 #include "twoPhaseSystem.H" 39 JohnsonJacksonParticleSlipFvPatchVectorField
50 const DimensionedField<vector, volMesh>& iF
53 partialSlipFvPatchVectorField(
p, iF),
54 specularityCoefficient_(
"specularityCoefficient",
dimless, 0)
61 const JohnsonJacksonParticleSlipFvPatchVectorField& ptf,
63 const DimensionedField<vector, volMesh>& iF,
64 const fvPatchFieldMapper& mapper
67 partialSlipFvPatchVectorField(ptf,
p, iF, mapper),
68 specularityCoefficient_(ptf.specularityCoefficient_)
76 const DimensionedField<vector, volMesh>& iF,
77 const dictionary&
dict 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" 101 const JohnsonJacksonParticleSlipFvPatchVectorField& ptf
104 partialSlipFvPatchVectorField(ptf),
105 specularityCoefficient_(ptf.specularityCoefficient_)
112 const JohnsonJacksonParticleSlipFvPatchVectorField& ptf,
113 const DimensionedField<vector, volMesh>& iF
116 partialSlipFvPatchVectorField(ptf, iF),
117 specularityCoefficient_(ptf.specularityCoefficient_)
125 const fvPatchFieldMapper& m
128 partialSlipFvPatchVectorField::autoMap(m);
138 partialSlipFvPatchVectorField::rmap(ptf, addr);
150 const twoPhaseSystem&
fluid = db().lookupObject<twoPhaseSystem>
155 const phaseModel& phased
166 phased.volScalarField::name()
177 patch().lookupPatchField<volScalarField>
185 patch().lookupPatchField<volScalarField>
195 db().foundObject<volScalarField>(ThetaName)
196 ?
patch().lookupPatchField<volScalarField>(ThetaName)
206 .lookupObject<IOdictionary>
211 .subDict(
"kineticTheoryCoeffs")
220 *specularityCoefficient_.
value()
225 this->valueFraction() =
c/(
c +
patch().deltaCoeffs());
227 partialSlipFvPatchVectorField::updateCoeffs();
237 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...
#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)
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.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
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.
constexpr scalar pi(M_PI)
errorManip< error > abort(error &err)
OBJstream os(runTime.globalPath()/outputName)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const dimensionedScalar c
Speed of light in a vacuum.
const std::string patch
OpenFOAM patch number as a std::string.
const word & name() const
List< label > labelList
A List of labels.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)