def test_LET_pc_imb():
wo = WaterOil(swirr=0.1)
wo.add_LET_pc_imb(Ls=1, Es=1, Ts=1, Lf=1, Ef=1, Tf=1, Pcmax=10, Pcmin=-10, Pct=3)
assert np.isclose(wo.table["pc"].max(), 10)
assert np.isclose(wo.table["pc"].min(), -10)
wo = WaterOil(swirr=0.1)
wo.add_LET_pc_imb(Ls=5, Es=5, Ts=5, Lf=5, Ef=5, Tf=5, Pcmax=5, Pcmin=1, Pct=4)
assert np.isclose(wo.table["pc"].max(), 5)
assert np.isclose(wo.table["pc"].min(), 1)
wo = WaterOil(swirr=0.1, sorw=0.3)
wo.add_LET_pc_imb(Ls=5, Es=5, Ts=5, Lf=5, Ef=5, Tf=5, Pcmax=5, Pcmin=1, Pct=4)
assert np.isclose(wo.table["pc"].max(), 5)
assert np.isclose(wo.table["pc"].min(), 1)
def test_let_pc_imb():
"""Test the LET formulation for imbibition capillary pressures"""
wateroil = WaterOil(swirr=0.1)
wateroil.add_LET_pc_imb(Ls=1,
Es=1,
Ts=1,
Lf=1,
Ef=1,
Tf=1,
Pcmax=10,
Pcmin=-10,
Pct=3)
assert np.isclose(wateroil.table["pc"].max(), 10)
assert np.isclose(wateroil.table["pc"].min(), -10)
wateroil = WaterOil(swirr=0.1)
wateroil.add_LET_pc_imb(Ls=5,
Es=5,
Ts=5,
Lf=5,
Ef=5,
Tf=5,
Pcmax=5,
Pcmin=1,
Pct=4)
assert np.isclose(wateroil.table["pc"].max(), 5)
assert np.isclose(wateroil.table["pc"].min(), 1)
wateroil = WaterOil(swirr=0.1, sorw=0.3)
wateroil.add_LET_pc_imb(Ls=5,
Es=5,
Ts=5,
Lf=5,
Ef=5,
Tf=5,
Pcmax=5,
Pcmin=1,
Pct=4)
assert np.isclose(wateroil.table["pc"].max(), 5)
assert np.isclose(wateroil.table["pc"].min(), 1)
wateroil.add_corey_water()
wateroil.add_corey_oil()
sat_table_str_ok(wateroil.SWOF())
