MEA solvent viscosity calculation correction (#785)

Co-authored-by: Andrew Lee <[email protected]>
IDAES · Apr 5, 2022 · 8310071 · 8310071
1 parent 1d433e9
commit 8310071
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 4 deletions.
diff --git a/idaes/models_extra/column_models/properties/MEA_solvent.py b/idaes/models_extra/column_models/properties/MEA_solvent.py
@@ -517,7 +517,9 @@ def build_parameters(pobj):
     def return_expression(blk, phase):
         pobj = blk.params.get_phase(phase)
 
-        r = blk.mass_frac_phase_comp_apparent["Liq", "MEA"] * 100
+        r = (blk.mass_frac_phase_comp_apparent["Liq", "MEA"] /
+             (blk.mass_frac_phase_comp_apparent["Liq", "MEA"] +
+              blk.mass_frac_phase_comp_apparent["Liq", "H2O"])) * 100
         T = blk.temperature
         alpha = (
             blk.mole_frac_phase_comp_apparent["Liq", "CO2"]

diff --git a/idaes/models_extra/column_models/properties/tests/test_mea_solvent.py b/idaes/models_extra/column_models/properties/tests/test_mea_solvent.py
@@ -135,7 +135,7 @@ def test_solution(self, model):
             model.props[1].vol_mol_phase["Liq"]
         )
 
-        assert pytest.approx(5.26403457e-4, rel=1e-8) == value(
+        assert pytest.approx(5.58038484e-4, rel=1e-8) == value(
             model.props[1].visc_d_phase["Liq"]
         )
 
@@ -163,9 +163,9 @@ def test_solution(self, model):
         assert pytest.approx(4.47017415e-09, rel=1e-8) == value(
             model.props[1].diffus_phase_comp_true["Liq", "MEA"]
         )
-        assert pytest.approx(2.27073817e-09, rel=1e-8) == value(
+        assert pytest.approx(2.17984326e-09, rel=1e-8) == value(
             model.props[1].diffus_phase_comp_true["Liq", "MEACOO_-"]
         )
-        assert pytest.approx(2.27073817e-09, rel=1e-8) == value(
+        assert pytest.approx(2.17984326e-09, rel=1e-8) == value(
             model.props[1].diffus_phase_comp_true["Liq", "MEA_+"]
         )