35 const volScalarField&
rho,
36 const volVectorField&
U,
37 const surfaceScalarField&
phi,
38 const surfaceScalarField&
rhoPhi,
40 const twoPhaseMixtureThermo&
mixture
43 twoPhaseTransport_(false),
49 IOdictionary turbulenceProperties
53 turbulenceModel::propertiesName,
61 const word simulationType
63 turbulenceProperties.get<word>(
"simulationType")
66 if (simulationType ==
"twoPhaseTransport")
68 twoPhaseTransport_ =
true;
72 if (twoPhaseTransport_)
84 IOobject::groupName(
"alphaRhoPhi",
alpha1.group()),
85 fvc::interpolate(
rho1)*alphaPhi10_
93 IOobject::groupName(
"alphaRhoPhi",
alpha2.group()),
94 fvc::interpolate(
rho2)*(phi_ - alphaPhi10_)
100 ThermalDiffusivity<PhaseCompressibleTurbulenceModel<fluidThermo>>
114 ThermalDiffusivity<PhaseCompressibleTurbulenceModel<fluidThermo>>
128 turbulence_ = compressible::turbulenceModel::New
136 turbulence_->validate();
146 if (twoPhaseTransport_)
151 turbulence1_->alphat()
155 turbulence2_->alphat()
160 return mixture_.
alphaEff(turbulence_->alphat());
171 if (twoPhaseTransport_)
174 turbulence1_->divDevRhoReff(
U)
175 + turbulence2_->divDevRhoReff(
U);
179 return turbulence_->divDevRhoReff(
U);
186 if (twoPhaseTransport_)
199 if (twoPhaseTransport_)
201 turbulence1_->correct();
202 turbulence2_->correct();
206 turbulence_->correct();
const volScalarField & alpha1
const volScalarField & alpha2
Internal & ref(const bool updateAccessTime=true)
Return a reference to the dimensioned internal field.
virtual tmp< volScalarField > alphaEff(const volScalarField &alphat) const =0
Effective thermal diffusivity of mixture [kg/m/s].
Transport model selection class for the compressibleInterFoam family of solvers.
void correct()
Correct the phase or mixture transport models.
tmp< volScalarField > alphaEff() const
Return the effective temperature transport coefficient.
void correctPhasePhi()
Correct the phase mass-fluxes.
tmp< volVectorField > divDevRhoReff()
Return the effective viscous stress (laminar + turbulent).
A class for managing temporary objects.
const rhoThermo & thermo1() const
const rhoThermo & thermo2() const
virtual tmp< volScalarField > alphaEff(const volScalarField &alphat) const
Effective thermal diffusivity of mixture [J/m/s/K].
const volScalarField & alpha2() const
Return the phase-fraction of phase 2.
const volScalarField & alpha1() const
Return the phase-fraction of phase 1.
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
GeometricField< vector, fvPatchField, volMesh > volVectorField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture