41 mixedFvPatchScalarField(
p, iF,
dict, false),
43 relax_(
dict.lookupOrDefault<scalar>(
"relax", 1))
46 refValue() = patchInternalField();
62 mixedFvPatchScalarField(psf,
p, iF, mapper),
75 mixedFvPatchScalarField(psf, iF),
92 const scalar q = q_->value(db().time().userTimeValue());
100 scalarField sumAlphaKappaEff(patch().size(), rootVSmall);
101 scalarField sumNotThisAlphaKappaEff(patch().size(), rootVSmall);
102 scalarField sumNotThisAlphaKappaEffT(patch().size(), rootVSmall);
109 if (&phase != &thisPhase)
117 sumAlphaKappaEff += alphaKappaEff;
118 sumNotThisAlphaKappaEff += alphaKappaEff;
119 sumNotThisAlphaKappaEffT += alphaKappaEff*
T.patchInternalField();
128 thisPhase.thermo().T().boundaryField()[
patchi];
130 sumAlphaKappaEff += alphaKappaEff;
131 sumNotThisAlphaKappaEff =
134 sumNotThisAlphaKappaEff,
137 sumNotThisAlphaKappaEffT =
140 sumNotThisAlphaKappaEffT,
141 rootSmall*kappaEff*
T.patchInternalField()
145 valueFraction() = sumNotThisAlphaKappaEff/sumAlphaKappaEff;
146 refValue() = sumNotThisAlphaKappaEffT/sumNotThisAlphaKappaEff;
147 refGrad() = q/
max(
alpha, rootSmall)/kappaEff;
153 valueFraction() = 1 - relax_*(1 -
f);
154 refValue() = (
f*relax_*refValue() + (1 - relax_)*
T)/valueFraction();
158 mixedFvPatchScalarField::updateCoeffs();
#define forAll(list, i)
Loop across all elements in list.
Macros for easy insertion into run-time selection tables.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
Run-time selectable general function of one variable.
const Boundary & boundaryField() const
Return const-reference to the boundary field.
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
virtual const volScalarField & T() const =0
Temperature [K].
A list of keyword definitions, which are a keyword followed by any number of values (e....
Foam::fvPatchFieldMapper.
Abstract base class with a fat-interface to all derived classes covering all possible ways in which t...
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
virtual tmp< scalarField > kappaEff(const label patchi) const =0
Effective thermal turbulent conductivity.
virtual const rhoThermo & thermo() const =0
Return the thermophysical model.
Class to represent a system of phases and model interfacial transfers between them.
static const word propertiesName
Default name of the phase properties dictionary.
const phaseModelList & phases() const
Return the phase models.
volScalarField alpha(IOobject("alpha", runTime.name(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
autoPtr< CompressibleMomentumTransportModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const viscosity &viscosity)
void write(std::ostream &os, const bool binary, List< floatScalar > &fField)
Write floats ascii or binary.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
tmp< fvMatrix< Type > > operator==(const fvMatrix< Type > &, const fvMatrix< Type > &)
void writeEntry(Ostream &os, const HashTable< T, Key, Hash > &ht)
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)