Improve logic for liquid vapourization enthalpy

davidschlachter · Nov 15, 2017 · 1aed95d · 1aed95d
1 parent 7df175b
commit 1aed95d
Showing 1 changed file with 9 additions and 2 deletions.
diff --git a/Enthalpy.java b/Enthalpy.java
@@ -93,8 +93,15 @@ public double testFunction(double testTemp){
           } else {
             lBound = outletBubbleTemperature;
           }
-          result += inletLiquidFlowRate*HeatCapacity.integrate(this.inlet.getFlowSpecies().get(i),  initialTemperature, lBound, "liquid");
-          result += outletVapourFlowRate*HeatCapacity.integrate(this.inlet.getFlowSpecies().get(i), lBound, finalTemperature, "vapour");
+          // Liquid that was vapourized
+          result += (inletLiquidFlowRate-outletLiquidFlowRate)*HeatCapacity.integrate(this.inlet.getFlowSpecies().get(i),  initialTemperature, lBound, "liquid");
+          // Liquid that was heated but not vapourized
+          result += outletLiquidFlowRate*HeatCapacity.integrate(this.inlet.getFlowSpecies().get(i),  initialTemperature, finalTemperature, "liquid");
+          // Vapour that was formed
+          result += (outletVapourFlowRate-inletVapourFlowRate)*HeatCapacity.integrate(this.inlet.getFlowSpecies().get(i), lBound, finalTemperature, "vapour");
+          // Vapour that did not change state (just heated)
+          result += inletVapourFlowRate*HeatCapacity.integrate(this.inlet.getFlowSpecies().get(i), lBound, finalTemperature, "vapour");
+          // Heat of vapourization for vapourized liquid
           result += (inletLiquidFlowRate-outletLiquidFlowRate)*heatofVapourization;
         }
         // If some of the vapour has condensed