def test_OP_plotting():
ctrl.clear()
pn = OpenPNM.Network.Cubic(shape=[30, 30, 1], spacing=0.01)
geom = OpenPNM.Geometry.Toray090(network=pn, pores=pn.Ps, throats=pn.Ts)
water = OpenPNM.Phases.Water(network=pn)
OpenPNM.Physics.Standard(network=pn, phase=water, geometry=geom)
inlets = pn.pores('left')
OP = OrdinaryPercolation(network=pn)
OP.setup(invading_phase=water)
OP.set_inlets(pores=inlets)
OP.run(npts=25)
a = OP.plot_drainage_curve()
assert isinstance(a, matplotlib.figure.Figure)
a = OP.plot_primary_drainage_curve()
assert isinstance(a, matplotlib.figure.Figure)
def test_OP_plotting():
mgr.clear()
pn = OpenPNM.Network.Cubic(shape=[30, 30, 1], spacing=0.01)
geom = OpenPNM.Geometry.Toray090(network=pn, pores=pn.Ps, throats=pn.Ts)
water = OpenPNM.Phases.Water(network=pn)
OpenPNM.Physics.Standard(network=pn, phase=water, geometry=geom)
inlets = pn.pores('left')
OP = OrdinaryPercolation(network=pn)
OP.setup(invading_phase=water)
OP.set_inlets(pores=inlets)
OP.run(npts=25)
a = OP.plot_drainage_curve()
assert isinstance(a, matplotlib.figure.Figure)
a = OP.plot_primary_drainage_curve()
assert isinstance(a, matplotlib.figure.Figure)
def test_OP_new_approach():
mgr.clear()
pn = OpenPNM.Network.Cubic(shape=[30, 30, 1], spacing=0.01)
geom = OpenPNM.Geometry.Toray090(network=pn, pores=pn.Ps, throats=pn.Ts)
water = OpenPNM.Phases.Water(network=pn)
phys = OpenPNM.Physics.Standard(network=pn, phase=water, geometry=geom)
inlets = pn.pores('left')
OP = OrdinaryPercolation(network=pn)
OP.setup(invading_phase=water)
OP.set_inlets(pores=inlets)
OP.run(npts=25)
a = ['pore.all', 'pore.inlets', 'pore.inv_Pc', 'pore.inv_sat',
'pore.inv_seq', 'throat.all', 'throat.entry_pressure',
'throat.inv_Pc', 'throat.inv_sat', 'throat.inv_seq']
assert sorted(list(OP.keys())) == a
V_inv = sp.sum(pn['pore.volume'][OP['pore.inv_Pc'] < sp.inf])
V_tot = sp.sum(pn['pore.volume'])
assert V_inv/V_tot == 1.0
mgr.clear()
def test_OP_new_approach():
ctrl.clear()
pn = OpenPNM.Network.Cubic(shape=[30, 30, 1], spacing=0.01)
geom = OpenPNM.Geometry.Toray090(network=pn, pores=pn.Ps, throats=pn.Ts)
water = OpenPNM.Phases.Water(network=pn)
phys = OpenPNM.Physics.Standard(network=pn, phase=water, geometry=geom)
inlets = pn.pores('left')
OP = OrdinaryPercolation(network=pn)
OP.setup(invading_phase=water)
OP.set_inlets(pores=inlets)
OP.run(npts=25)
a = [
'pore.all', 'pore.inlets', 'pore.inv_Pc', 'pore.inv_sat',
'pore.inv_seq', 'throat.all', 'throat.entry_pressure', 'throat.inv_Pc',
'throat.inv_sat', 'throat.inv_seq'
]
assert sorted(list(OP.keys())) == a
V_inv = sp.sum(pn['pore.volume'][OP['pore.inv_Pc'] < sp.inf])
V_tot = sp.sum(pn['pore.volume'])
assert V_inv / V_tot == 1.0
ctrl.clear()
