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_trapping():
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)
air = OpenPNM.Phases.Air(network=pn)
phys = OpenPNM.Physics.GenericPhysics(network=pn,
geometry=geom,
phase=water)
f = OpenPNM.Physics.models.capillary_pressure.washburn
phys.models.add(propname='throat.capillary_pressure', model=f)
OP = OrdinaryPercolation(network=pn,
invading_phase=water,
defending_phase=air)
OP.run(inlets=pn.pores('left'), outlets=pn.pores('right'), trapping=True)
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()
def test_OP_trapping():
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)
air = OpenPNM.Phases.Air(network=pn)
phys = OpenPNM.Physics.GenericPhysics(network=pn,
geometry=geom,
phase=water)
f = OpenPNM.Physics.models.capillary_pressure.washburn
phys.models.add(propname='throat.capillary_pressure',
model=f)
OP = OrdinaryPercolation(network=pn,
invading_phase=water,
defending_phase=air)
OP.run(inlets=pn.pores('left'), outlets=pn.pores('right'), trapping=True)
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()
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_old_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)
OP_1 = OrdinaryPercolation(network=pn, invading_phase=water)
Ps = pn.pores(labels=['left'])
OP_1.run(inlets=Ps)
OP_1.return_results(Pc=7000)
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_1.keys())) == a
ctrl.clear()
def test_OP_old_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)
OP_1 = OrdinaryPercolation(network=pn, invading_phase=water)
Ps = pn.pores(labels=['left'])
OP_1.run(inlets=Ps)
OP_1.return_results(Pc=7000)
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_1.keys())) == a
mgr.clear()