26 #include "JohnsonJacksonParticleSlipFvPatchVectorField.H" 28 #include "twoPhaseSystem.H" 37 JohnsonJacksonParticleSlipFvPatchVectorField
48 const DimensionedField<vector, volMesh>& iF
51 partialSlipFvPatchVectorField(p, iF),
52 specularityCoefficient_(
"specularityCoefficient",
dimless, 0)
59 const JohnsonJacksonParticleSlipFvPatchVectorField& ptf,
61 const DimensionedField<vector, volMesh>& iF,
62 const fvPatchFieldMapper& mapper
65 partialSlipFvPatchVectorField(ptf, p, iF, mapper),
66 specularityCoefficient_(ptf.specularityCoefficient_)
74 const DimensionedField<vector, volMesh>& iF,
75 const dictionary& dict
78 partialSlipFvPatchVectorField(p, iF),
79 specularityCoefficient_
81 "specularityCoefficient",
83 dict.
lookup(
"specularityCoefficient")
88 (specularityCoefficient_.value() < 0)
89 || (specularityCoefficient_.value() > 1)
93 <<
"The specularity coefficient has to be between 0 and 1" 107 const JohnsonJacksonParticleSlipFvPatchVectorField& ptf
110 partialSlipFvPatchVectorField(ptf),
111 specularityCoefficient_(ptf.specularityCoefficient_)
118 const JohnsonJacksonParticleSlipFvPatchVectorField& ptf,
119 const DimensionedField<vector, volMesh>& iF
122 partialSlipFvPatchVectorField(ptf, iF),
123 specularityCoefficient_(ptf.specularityCoefficient_)
131 const fvPatchFieldMapper& m
134 partialSlipFvPatchVectorField::autoMap(m);
144 partialSlipFvPatchVectorField::rmap(ptf, addr);
156 const twoPhaseSystem& fluid = db().lookupObject<twoPhaseSystem>
161 const phaseModel& phased
163 fluid.phase1().name() == internalField().group()
171 patch().lookupPatchField<volScalarField, scalar>
173 phased.volScalarField::name()
179 patch().lookupPatchField<volScalarField, scalar>
187 patch().lookupPatchField<volScalarField, scalar>
195 patch().lookupPatchField<volScalarField, scalar>
205 db().foundObject<volScalarField>(ThetaName)
206 ? patch().lookupPatchField<volScalarField, scalar>(ThetaName)
216 .lookupObject<IOdictionary>
221 .subDict(
"kineticTheoryCoeffs")
231 *specularityCoefficient_.
value()
236 this->valueFraction() =
c/(
c + patch().deltaCoeffs());
238 partialSlipFvPatchVectorField::updateCoeffs();
248 os.writeKeyword(
"specularityCoefficient")
250 writeEntry(
"value", os);
virtual void autoMap(const fvPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.
fvPatchField< vector > fvPatchVectorField
JohnsonJacksonParticleSlipFvPatchVectorField(const fvPatch &, const DimensionedField< vector, volMesh > &)
Construct from patch and internal field.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
friend Ostream & operator(Ostream &, const fvPatchField< Type > &)
virtual void rmap(const fvPatchVectorField &, const labelList &)
Reverse map the given fvPatchField onto this fvPatchField.
dimensionedScalar sqrt(const dimensionedScalar &ds)
const Type & value() const
Return const reference to value.
Macros for easy insertion into run-time selection tables.
makePatchTypeField(fvPatchVectorField, SRFFreestreamVelocityFvPatchVectorField)
stressControl lookup("compactNormalStress") >> compactNormalStress
fvPatchField< scalar > fvPatchScalarField
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
static word groupName(Name name, const word &group)
List< label > labelList
A List of labels.
errorManip< error > abort(error &err)
dimensionedScalar alphaMax(laminarTransport.lookup("alphaMax"))
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const dimensionedScalar c
Speed of light in a vacuum.
virtual void updateCoeffs()
Update the coefficients.
Field< vector > vectorField
Specialisation of Field<T> for vector.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
virtual void write(Ostream &) const
Write.
virtual void write(Ostream &) const
Write.