/**
* Compute average hydrocarbon pressure and maximum
* dissolution and evaporation at average hydrocarbon
* pressure in all regions in field.
*
* Fluid properties are evaluated at average hydrocarbon
* pressure for purpose of conversion from surface rate to
* reservoir voidage rate.
*
* \param[in] state Dynamic reservoir state.
*/
void
defineState(const BlackoilState& state)
{
averagePressure(state);
averageTemperature(state);
calcRmax();
}
scalar pressureSwirlInjector::averageVelocity
(
const label i
) const
{
const injectorType& it = sm_.injectors()[i].properties();
scalar dt = it.teoi() - it.tsoi();
scalar injectedMassFlow = it.mass()/(it.teoi()-it.tsoi());
scalar injectionPressure = averagePressure(i);
scalar Tav = it.integrateTable(it.T())/dt;
scalar rhoFuel = sm_.fuels().rho(sm_.ambientPressure(), Tav, it.X());
scalar kV = kv(i, injectedMassFlow, injectionPressure);
return kV*sqrt(2.0*(injectionPressure-sm_.ambientPressure())/rhoFuel);
}
