46 psiThermo(mesh, word::null),
47 twoPhaseMixture(mesh, *this),
64 thermo1_->validate(phase1Name(),
"e");
65 thermo2_->validate(phase2Name(),
"e");
81 thermo1_->he() = thermo1_->he(
p_,
T_);
84 thermo2_->he() = thermo2_->he(
p_,
T_);
95 return thermo1_->incompressible() && thermo2_->incompressible();
101 return thermo1_->isochoric() && thermo2_->isochoric();
111 return alpha1()*thermo1_->he(p, T) +
alpha2()*thermo2_->he(p, T);
143 return alpha1()*thermo1_->hc() +
alpha2()*thermo2_->hc();
175 return alpha1()*thermo1_->Cp() +
alpha2()*thermo2_->Cp();
194 return alpha1()*thermo1_->Cv() +
alpha2()*thermo2_->Cv();
213 return alpha1()*thermo1_->gamma() +
alpha2()*thermo2_->gamma();
232 return alpha1()*thermo1_->Cpv() +
alpha2()*thermo2_->Cpv();
252 alpha1()*thermo1_->CpByCpv()
253 +
alpha2()*thermo2_->CpByCpv();
272 return mu()/(
alpha1()*thermo1_->rho() +
alpha2()*thermo2_->rho());
292 return alpha1()*thermo1_->kappa() +
alpha2()*thermo2_->kappa();
313 alpha1()*thermo1_->kappaEff(alphat)
314 +
alpha2()*thermo2_->kappaEff(alphat);
337 alpha1()*thermo1_->alphaEff(alphat)
338 +
alpha2()*thermo2_->alphaEff(alphat);
virtual bool isochoric() const
Return true if the equation of state is isochoric.
virtual bool incompressible() const
Return true if the equation of state is incompressible.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
virtual tmp< volScalarField > Cpv() const
Heat capacity at constant pressure/volume [J/kg/K].
virtual tmp< volScalarField > Cp() const
Heat capacity at constant pressure [J/kg/K].
virtual tmp< volScalarField > Cv() const
Heat capacity at constant volume [J/kg/K].
volScalarField T_
Temperature [K].
virtual tmp< volScalarField > CpByCpv() const
Heat capacity ratio [].
virtual tmp< volScalarField > kappaEff(const volScalarField &alphat) const
Effective thermal diffusivity of mixture [J/m/s/K].
const volScalarField & alpha1() const
Return the phase-fraction of phase 1.
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
virtual tmp< volScalarField > mu() const
Dynamic viscosity of mixture [kg/m/s].
GeometricField< scalar, fvPatchField, volMesh > volScalarField
virtual tmp< volScalarField > kappa() const
Thermal diffusivity for temperature of mixture [J/m/s/K].
virtual tmp< volScalarField > gamma() const
Gamma = Cp/Cv [].
const Boundary & boundaryField() const
Return const-reference to the boundary field.
volScalarField & p_
Pressure [Pa].
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
virtual ~twoPhaseMixtureThermo()
Destructor.
List< label > labelList
A List of labels.
virtual tmp< volScalarField > hc() const
Chemical enthalpy [J/kg].
virtual tmp< scalarField > THE(const scalarField &h, const scalarField &p, const scalarField &T0, const labelList &cells) const
Temperature from enthalpy/internal energy for cell-set.
Info<< "Predicted p max-min : "<< max(p).value()<< " "<< min(p).value()<< endl;rho==max(psi *p+alphal *rhol0+((alphav *psiv+alphal *psil)-psi)*pSat, rhoMin);#1"/home/ubuntu/OpenFOAM-4.1/applications/solvers/multiphase/cavitatingFoam/alphavPsi.H"1{alphav=max(min((rho-rholSat)/(rhovSat-rholSat), scalar(1)), scalar(0));alphal=1.0-alphav;Info<< "max-min alphav: "<< max(alphav).value()<< " "<< min(alphav).value()<< endl;psiModel-> correct()
defineTypeNameAndDebug(combustionModel, 0)
virtual volScalarField & he()
Enthalpy/Internal energy [J/kg].
const dimensionedScalar h
Planck constant.
volScalarField alpha_
Laminar thermal diffusuvity [kg/m/s].
const volScalarField & alpha2() const
Return the phase-fraction of phase 2.
virtual void correct()
Update properties.
volScalarField mu_
Dynamic viscosity [kg/m/s].
virtual tmp< volScalarField > alphaEff(const volScalarField &alphat) const
Effective thermal diffusivity of mixture [J/m/s/K].
volScalarField psi_
Compressibility [s^2/m^2].
A class for managing temporary objects.
#define NotImplemented
Issue a FatalErrorIn for a function not currently implemented.
twoPhaseMixtureThermo(const fvMesh &mesh)
Construct from mesh.
virtual tmp< volScalarField > nu() const
Kinematic viscosity of mixture [m^2/s].