use rigorous VLE

Author: yoelcortes

How2STEAM

@@ -1481,7 +1481,7 @@ def _mass_and_energy_balance_specifications_table(self):
         return df
 
     def results(self, with_units=True, include_utilities=True,
-                include_total_cost=True, include_installed_cost=False,
+                include_total_cost=True, include_installed_cost=True,
                 include_zeros=True, external_utilities=None, key_hook=None,
                 basis=None):
         """

biosteam/_unit.py

@@ -325,7 +325,8 @@ def load_auxiliaries(self):
         )
 
     def _run_vent(self, vent, effluent):
-        vent.receive_vent(effluent, energy_balance=False, ideal=True)
+        stream = bst.MultiStream.from_streams([effluent, vent])
+        stream.vle(T=stream.T, P=stream.P)
 
     def _run(self):
         air = self.air
@@ -767,7 +768,8 @@ def load_auxiliaries(self):
         )
 
     def _run_vent(self, vent, effluent):
-        vent.receive_vent(effluent, energy_balance=False, ideal=True)
+        stream = bst.MultiStream.from_streams([effluent, vent])
+        stream.vle(T=stream.T, P=stream.P)
 
     def get_SURs(self, effluent):
         F_vol = effluent.ivol['Water'] # m3 / hr

biosteam/units/aerated_bioreactor.py

@@ -16,6 +16,7 @@
     & Ng, M. K. (2017). Product and Process Design Principles. Wiley.
 
 """
+import biosteam as bst
 from .abstract_stirred_tank_reactor import AbstractStirredTankReactor
 
 __all__ = (
@@ -190,7 +191,8 @@ def vent(self):
 
     def _run_vent(self, vent, effluent):
         vent.phase = 'g'
-        vent.receive_vent(effluent, energy_balance=False, ideal=True)
+        stream = bst.MultiStream.from_streams([effluent, vent])
+        stream.vle(T=stream.T, P=stream.P)
 
     def _run(self):
         vent, effluent = self.outs

biosteam/units/anaerobic_bioreactor.py

@@ -47,6 +47,7 @@
 from math import ceil
 from scipy.integrate import odeint
 from thermosteam.reaction import Reaction, ParallelReaction
+import biosteam as bst
 
 __all__ = (
     'NRELAnaerobicBatchBioreactor', 
@@ -130,7 +131,8 @@ def _get_design_info(self):
                 ('Working volume fraction', self.V_wf, ''))
 
     def _init(self, tau=None, N=None, V_max=None, T=305.15, P=101325,
-              Nmin=2, Nmax=36, loading_time=None):
+              Nmin=2, Nmax=36, loading_time=None, reactions=None, 
+              tau_0=None, V_wf=None):
         self._N = N; self._V_max = V_max
 
         #: [float] Reaction time [hr].
@@ -158,6 +160,12 @@ def _init(self, tau=None, N=None, V_max=None, T=305.15, P=101325,
         #: If not given, it will assume each vessel is constantly being filled.
         self.loading_time = loading_time
 
+        #: [Reaction|ReactionSystem] Stoichiomentric reactions.
+        self.reactions = reactions
+        
+        if tau_0 is not None: self.tau_0 = tau_0
+        if V_wf is not None: self.V_wf = V_wf
+        
     def _setup(self):
         super()._setup()
         vent, effluent = self.outs
@@ -210,6 +218,13 @@ def tau(self):
     def tau(self, tau):
         self._tau = tau
 
+    def _run(self):
+        vent, effluent = self.outs
+        effluent.mix_from(self.ins)
+        self.reactions(effluent)
+        stream = bst.MultiStream.from_streams([effluent, vent])
+        stream.vle(T=stream.T, P=stream.P)
+        
     def _design(self):
         effluent = self.effluent
         v_0 = effluent.F_vol