59 lines
1.8 KiB
Python
59 lines
1.8 KiB
Python
|
import gaspype as gp
|
||
|
|
# import gaspype.membrane as mb
|
||
|
|
import pytest
|
||
|
|
|
||
|
|
|
||
|
|
def test_oxygen_partial_pressure():
|
||
|
|
# Compare equalibrium calculation with oxygen_partial_pressure function
|
||
|
|
fs = gp.fluid_system('CO, CO2, H2, O2, H2O, N2')
|
||
|
|
fl = gp.fluid({'H2O': 0.99, 'H2': 0.01}, fs)
|
||
|
|
|
||
|
|
t = 2500
|
||
|
|
|
||
|
|
ox1 = gp.equilibrium(fl, t, 1e5).get_x('O2') * 1e5 # Pa
|
||
|
|
print(ox1)
|
||
|
|
|
||
|
|
ox2 = gp.oxygen_partial_pressure(fl, t, 1e5)
|
||
|
|
print(ox2)
|
||
|
|
|
||
|
|
assert ox1 == pytest.approx(ox2, abs=1e-1) # type: ignore
|
||
|
|
|
||
|
|
|
||
|
|
"""
|
||
|
|
def test_voltage_potential():
|
||
|
|
fs = gp.fluid_system('CO, CO2, H2, O2, H2O, N2')
|
||
|
|
fl1 = gp.fluid({'H2O': 0.5, 'H2': 0.5}, fs)
|
||
|
|
fl2 = gp.fluid({'O2': 0.2, 'N2': 0.8}, fs)
|
||
|
|
|
||
|
|
potential = mb.voltage_potential(fl1, fl2, 300, 1e5) # V
|
||
|
|
#assert potential == pytest.approx(1.1736986, abs=1e-4)
|
||
|
|
|
||
|
|
potential = mb.voltage_potential(fl1, fl2, 273.15 + 800, 1e5) # V
|
||
|
|
assert potential == pytest.approx(0.9397330, abs=1e-4)
|
||
|
|
|
||
|
|
fl1 = gp.fluid({'H2O': 0.01, 'H2': 0.99}, fs)
|
||
|
|
fl2 = gp.fluid({'O2': 0.99, 'N2': 0.01}, fs)
|
||
|
|
|
||
|
|
potential = mb.voltage_potential(fl1, fl2, 273.15 + 800, 1e5) # V
|
||
|
|
#assert potential == pytest.approx(1.18918092, abs=1e-4)
|
||
|
|
|
||
|
|
fl1 = gp.fluid({'H2O': 0.90, 'H2': 0.10}, fs)
|
||
|
|
fl2 = gp.fluid({'O2': 0.20, 'N2': 0.80}, fs)
|
||
|
|
|
||
|
|
potential = mb.voltage_potential(fl1, fl2, 273.15 + 800, 1e5) # V
|
||
|
|
assert potential == pytest.approx(0.83813680, abs=1e-4)
|
||
|
|
|
||
|
|
|
||
|
|
|
||
|
|
def test_voltage_potential_without_o2():
|
||
|
|
fl1 = gp.fluid({'H2O': 0.5, 'CO': 0, 'CO2': 0, 'H2': 0.5})
|
||
|
|
fs = gp.fluid_system('CO, CO2, H2, O2, H2O, N2')
|
||
|
|
fl2 = gp.fluid({'O2': 0.2, 'N2': 0.8}, fs)
|
||
|
|
|
||
|
|
potential = mb.voltage_potential(fl1, fl2, 300, 1e5) # V
|
||
|
|
#assert potential == pytest.approx(1.1736986, abs=1e-4)
|
||
|
|
|
||
|
|
potential = mb.voltage_potential(fl1, fl2, 273.15 + 800, 1e5) # V
|
||
|
|
#assert potential == pytest.approx(0.9397330, abs=1e-4)
|
||
|
|
"""
|