a little different form to keep the derivatives

This commit is contained in:
Kai Bao 2022-11-17 22:57:42 +01:00
parent 3ddb830ae9
commit 5e35477ce1

View File

@ -888,7 +888,7 @@ pressureDropSpiralICD(const int seg) const
const EvalWell liquid_fraction = water_fraction + oil_fraction;
// viscosity contribution from the liquid
const EvalWell liquid_viscosity_fraction = liquid_fraction < 1.e-30 ? 0 :
const EvalWell liquid_viscosity_fraction = liquid_fraction < 1.e-30 ? oil_fraction * oil_viscosity + water_fraction * water_viscosity :
liquid_fraction * mswellhelpers::emulsionViscosity(water_fraction, water_viscosity, oil_fraction, oil_viscosity, sicd);
const EvalWell mixture_viscosity = liquid_viscosity_fraction + gas_fraction * gas_viscosity;