doublereal Phase::massFraction(const std::string& nameSpec) const
{
size_t iloc = speciesIndex(nameSpec);
if (iloc != npos) {
return massFractions()[iloc];
} else {
return 0.0;
}
}
doublereal Phase::massFraction(std::string name) const {
int iloc = speciesIndex(name);
if (iloc >= 0) return massFractions()[iloc];
else return 0.0;
}
void
PorousFlowWaterNCG::thermophysicalProperties(Real pressure,
Real temperature,
Real /* Xnacl */,
Real Z,
unsigned int qp,
std::vector<FluidStateProperties> & fsp) const
{
FluidStateProperties & liquid = fsp[_aqueous_phase_number];
FluidStateProperties & gas = fsp[_gas_phase_number];
// Check whether the input temperature is within the region of validity
checkVariables(temperature);
// Clear all of the FluidStateProperties data
clearFluidStateProperties(fsp);
FluidStatePhaseEnum phase_state;
massFractions(pressure, temperature, Z, phase_state, fsp);
switch (phase_state)
{
case FluidStatePhaseEnum::GAS:
{
// Set the gas saturations
gas.saturation = 1.0;
// Calculate gas properties
gasProperties(pressure, temperature, fsp);
break;
}
case FluidStatePhaseEnum::LIQUID:
{
// Calculate the liquid properties
Real liquid_pressure = pressure - _pc.capillaryPressure(1.0, qp);
liquidProperties(liquid_pressure, temperature, fsp);
break;
}
case FluidStatePhaseEnum::TWOPHASE:
{
// Calculate the gas properties
gasProperties(pressure, temperature, fsp);
// Calculate the saturation
saturationTwoPhase(pressure, temperature, Z, fsp);
// Calculate the liquid properties
Real liquid_pressure = pressure - _pc.capillaryPressure(1.0 - gas.saturation, qp);
liquidProperties(liquid_pressure, temperature, fsp);
break;
}
}
// Liquid saturations can now be set
liquid.saturation = 1.0 - gas.saturation;
liquid.dsaturation_dp = -gas.dsaturation_dp;
liquid.dsaturation_dT = -gas.dsaturation_dT;
liquid.dsaturation_dZ = -gas.dsaturation_dZ;
// Save pressures to FluidStateProperties object
gas.pressure = pressure;
liquid.pressure = pressure - _pc.capillaryPressure(liquid.saturation, qp);
}
