def test_skjaeveland_pc():
"""Simple test of Skjæveland capillary pressure correlation"""
wateroil = WaterOil(h=0.3, swl=0.2)
wateroil.add_skjaeveland_pc(swr=0.1, cw=0.1, co=-0.1, aw=0.1, ao=0.1)
check_table(wateroil.table)
# Add with wrong numbers
wateroil = WaterOil(h=0.3, swl=0.2, sorw=0.3)
wateroil.add_skjaeveland_pc(swr=0.8, cw=-0.1, co=0.1, aw=-0.1, ao=-0.1)
# (the code returns None when errors occur, no Exception)
assert wateroil.pccomment == ""
assert "pc" not in wateroil.table
def test_skjaeveland_pc():
"""Simple test of Skjæveland capillary pressure correlation"""
wateroil = WaterOil(h=0.3, swl=0.2, swirr=0.1)
wateroil.add_corey_water()
wateroil.add_skjaeveland_pc(swr=0.1, cw=0.1, co=-0.1, aw=0.1, ao=0.1)
check_table(wateroil.table)
swfn = wateroil.SWFN()
# Defaulting swr:
wateroil.add_skjaeveland_pc(swr=0.1, cw=0.1, co=-0.1, aw=0.1, ao=0.1)
wateroil.add_skjaeveland_pc(swr=0.1, cw=0.1, co=-0.1, aw=0.1, ao=0.1)
assert wateroil.SWFN() == swfn
# Add with wrong numbers
wateroil = WaterOil(h=0.3, swl=0.2, sorw=0.3)
with pytest.raises(ValueError, match="cw must be larger"):
wateroil.add_skjaeveland_pc(swr=0.1, cw=-0.1, co=-0.1, aw=0.1, ao=0.1)
with pytest.raises(ValueError, match="co must be less"):
wateroil.add_skjaeveland_pc(swr=0.1, cw=0.1, co=0.1, aw=0.1, ao=0.1)
with pytest.raises(ValueError, match="aw must be larger"):
wateroil.add_skjaeveland_pc(swr=0.1, cw=0.1, co=-0.1, aw=-0.1, ao=0.1)
with pytest.raises(ValueError, match="ao must be larger"):
wateroil.add_skjaeveland_pc(swr=0.1,
sor=0.3,
cw=0.1,
co=-0.1,
aw=0.1,
ao=-0.1)
with pytest.raises(ValueError,
match=r"swr \(swirr\) must be less than 1 - sor"):
wateroil.add_skjaeveland_pc(swr=0.8,
sor=0.3,
cw=0.1,
co=-0.1,
aw=0.1,
ao=0.1)
with pytest.raises(ValueError, match="swr must be contained in"):
wateroil.add_skjaeveland_pc(swr=-1.2,
sor=0.3,
cw=0.1,
co=-0.1,
aw=0.1,
ao=0.1)
