implement reactive distillation; improve inside-out method

Author: yoelcortes

biosteam/_tea.py

@@ -403,7 +403,7 @@ def total_costs(self):
                 return sum([i[1] for i in self.data])
             else:
                 N = len(self.data[0])
-                return [sum([i[index] for i in self.data]) for index in range(1, N)]
+                return np.array([sum([i[index] for i in self.data]) for index in range(1, N)])
 
         def table(self):
             names = self.names

biosteam/process_tools/utils.py

@@ -426,8 +426,8 @@ def get_OSBL(units):
     
     """
     return [i for i in units
-            if isinstance(i, (bst.Facility, bst.StorageTank))
-            or 'storage' in i.line.lower() 
+            if isinstance(i, (bst.Facility, bst.StorageTank)) and not 'surge' in i.line.lower()
+            or 'storage' in i.line.lower()
             or 'wastewater' in i.line.lower()]
 
 def heat_exchanger_operation(units):

biosteam/units/distillation.py

@@ -2765,7 +2765,7 @@ class MESHDistillation(MultiStageEquilibrium, new_graphics=False):
 
     def _init(self, 
             LHK, N_stages, feed_stages, 
-            reflux=None, boilup=None, 
+            reflux=None, boilup=None, bottoms_product_to_feed=None,
             P=101325, 
             vapor_side_draws=None, liquid_side_draws=None,
             stage_reactions=None,
@@ -2806,6 +2806,8 @@ def _init(self,
             stage_specifications[0] = ('Reflux', reflux)
         if boilup is not None:
             stage_specifications[-1] = ('Boilup', boilup)
+        if bottoms_product_to_feed is not None:
+            stage_specifications[-1] = ('Flow', bottoms_product_to_feed)
         super()._init(N_stages=N_stages, feed_stages=feed_stages,
                       top_side_draws=vapor_side_draws, 
                       bottom_side_draws=liquid_side_draws,              
@@ -2988,7 +2990,7 @@ def update_liquid_holdup(self):
                 Ks = Ua * sqrt(rho_V / (rho_L - rho_V)) # Capacity parameter    
             Phi_e = exp(-4.257 * Ks**0.91) # Effective relative froth density 
             Lw = 0.73 * diameter * 12 # Weir length [in] assuming Ad/A = 0.1
-            # TODO: Compute weir length or other Ad/A
+            # TODO: Compute weir length for other Ad/A
             qL = liquid.get_total_flow('gal/min')
             CL = 0.362 + 0.317 * exp(-3.5 * weir_height)
             hL = Phi_e * (hw + CL * (qL / (Lw * Phi_e)) ** b) # equivalent height of clear liquid holdup [in]