Real PorousFlowCapillaryPressureBC::d2CapillaryPressureCurve(Real saturation, unsigned /*qp*/) const { Real seff = effectiveSaturationFromSaturation(saturation); return PorousFlowBrooksCorey::d2CapillaryPressure(seff, _pe, _lambda, _pc_max) * _dseff_ds * _dseff_ds; }
Real PorousFlowCapillaryPressureBC::capillaryPressureCurve(Real saturation) const { Real seff = effectiveSaturationFromSaturation(saturation); return PorousFlowBrooksCorey::capillaryPressure(seff, _pe, _lambda, _pc_max); }