added more units for reading flow rate of phase

equinor · Jul 15, 2023 · a0383a9 · a0383a9
1 parent 537ad12
commit a0383a9
Showing 1 changed file with 13 additions and 0 deletions.
diff --git a/src/main/java/neqsim/thermo/phase/Phase.java b/src/main/java/neqsim/thermo/phase/Phase.java
@@ -10,6 +10,7 @@
 import org.apache.logging.log4j.LogManager;
 import org.apache.logging.log4j.Logger;
 import neqsim.physicalProperties.PhysicalPropertyHandler;
+import neqsim.thermo.ThermodynamicConstantsInterface;
 import neqsim.thermo.component.ComponentInterface;
 import neqsim.thermo.system.SystemInterface;
 import neqsim.util.exception.InvalidInputException;
@@ -2137,6 +2138,18 @@ public double getFlowRate(String flowunit) {
       return numberOfMolesInPhase * 60.0;
     } else if (flowunit.equals("mole/hr")) {
       return numberOfMolesInPhase * 3600.0;
+    } else if (flowunit.equals("Sm3/sec")) {
+      return numberOfMolesInPhase * ThermodynamicConstantsInterface.R
+          * ThermodynamicConstantsInterface.standardStateTemperature / 101325.0;
+    } else if (flowunit.equals("Sm3/hr")) {
+      return numberOfMolesInPhase * 3600.0 * ThermodynamicConstantsInterface.R
+          * ThermodynamicConstantsInterface.standardStateTemperature / 101325.0;
+    } else if (flowunit.equals("Sm3/day")) {
+      return numberOfMolesInPhase * 3600.0 * 24.0 * ThermodynamicConstantsInterface.R
+          * ThermodynamicConstantsInterface.standardStateTemperature / 101325.0;
+    } else if (flowunit.equals("MSm3/day")) {
+      return numberOfMolesInPhase * 3600.0 * 24.0 * ThermodynamicConstantsInterface.R
+          * ThermodynamicConstantsInterface.standardStateTemperature / 101325.0 / 1.0e6;
     } else {
       throw new RuntimeException("failed.. unit: " + flowunit + " not supported");
     }