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()