KTA (#1290)

* Added new viscosity model for PURE HELIUM

* Modified KTA-viscosity model with lower uncertainty

Author: akselrs

src/main/java/neqsim/physicalproperties/methods/commonphasephysicalproperties/viscosity/KTAViscosityMethod.java

@@ -0,0 +1,37 @@
+package neqsim.physicalproperties.methods.commonphasephysicalproperties.viscosity;
+
+import neqsim.physicalproperties.system.PhysicalProperties;
+
+/**
+ * <p>
+ * KTAViscosityMethod class.
+ * </p>
+ */
+
+public class KTAViscosityMethod extends Viscosity {
+  /**
+   * <p>
+   * Constructor for KTAViscosityMethod.
+   * </p>
+   *
+   * @param phase a {@link neqsim.physicalproperties.system.PhysicalProperties} object
+   */
+  public KTAViscosityMethod(PhysicalProperties phase) {
+    super(phase);
+  }
+
+  /** {@inheritDoc} */
+  @Override
+  public double calcViscosity() {
+    // Check if there are other components than helium
+    if (phase.getPhase().getNumberOfComponents() > 1 
+      || !phase.getPhase().getComponent(0).getName().equalsIgnoreCase("helium")) {
+      throw new Error("This method only supports PURE HELIUM.");
+    }
+
+    double T = phase.getPhase().getTemperature();
+    //Source of KTA-model https://www.sciencedirect.com/science/article/pii/S0149197024004670#bib28
+    double visc = 3.674 * Math.pow(10, -7) * Math.pow(T, 0.7);    //[Pa*s]
+    return visc;
+  }
+}

src/main/java/neqsim/physicalproperties/methods/commonphasephysicalproperties/viscosity/KTAViscosityMethodMod.java

@@ -0,0 +1,42 @@
+package neqsim.physicalproperties.methods.commonphasephysicalproperties.viscosity;
+
+import neqsim.physicalproperties.system.PhysicalProperties;
+
+/**
+ * <p>
+ * KTAViscosityMethodMod class.
+ * </p>
+ */
+public class KTAViscosityMethodMod extends Viscosity {
+  /**
+   * <p>
+   * Constructor for KTAViscosityMethodMod.
+   * </p>
+   *
+   * @param phase a {@link neqsim.physicalproperties.system.PhysicalProperties} object
+   */
+  public KTAViscosityMethodMod(PhysicalProperties phase) {
+    super(phase);
+  }
+
+  /** {@inheritDoc} */
+  @Override
+  public double calcViscosity() {
+    // Check if there are other components than helium
+    if (phase.getPhase().getNumberOfComponents() > 1 
+      || !phase.getPhase().getComponent(0).getName().equalsIgnoreCase("helium")) {
+      throw new Error("This method only supports PURE HELIUM.");
+    }
+
+    double T = phase.getPhase().getTemperature();
+    double P = phase.getPhase().getPressure();
+    double P_crit = 0.22832;    // [MPa] (Source: NIST)
+    double A = Math.pow((2 - T / 300), 5.05);
+    double B = Math.pow((2 - 300 / T), 2) - 1;
+    double viscosity = 1e-7 * (3.817 * Math.pow(T, 0.6938) 
+                        + Math.pow(P, A) / (T * P_crit) 
+                        + Math.exp(-Math.pow(T - 325, 2) / 1000) 
+                          * (Math.pow(P / 25, 2.7) - Math.pow(T, B)));
+    return viscosity;   // [Pa*s]
+  }
+}

src/main/java/neqsim/physicalproperties/system/PhysicalProperties.java

@@ -12,6 +12,8 @@
 import neqsim.physicalproperties.methods.commonphasephysicalproperties.conductivity.PFCTConductivityMethodMod86;
 import neqsim.physicalproperties.methods.commonphasephysicalproperties.diffusivity.CorrespondingStatesDiffusivity;
 import neqsim.physicalproperties.methods.commonphasephysicalproperties.viscosity.FrictionTheoryViscosityMethod;
+import neqsim.physicalproperties.methods.commonphasephysicalproperties.viscosity.KTAViscosityMethod;
+import neqsim.physicalproperties.methods.commonphasephysicalproperties.viscosity.KTAViscosityMethodMod;
 import neqsim.physicalproperties.methods.commonphasephysicalproperties.viscosity.LBCViscosityMethod;
 import neqsim.physicalproperties.methods.commonphasephysicalproperties.viscosity.PFCTViscosityMethodHeavyOil;
 import neqsim.physicalproperties.methods.commonphasephysicalproperties.viscosity.PFCTViscosityMethodMod86;
@@ -219,6 +221,10 @@ public void setViscosityModel(String model) {
       viscosityCalc = new PFCTViscosityMethodMod86(this);
     } else if ("PFCT-Heavy-Oil".equals(model)) {
       viscosityCalc = new PFCTViscosityMethodHeavyOil(this);
+    } else if ("KTA".equals(model)) {
+      viscosityCalc = new KTAViscosityMethod(this);
+    } else if ("KTA_mod".equals(model)) {
+      viscosityCalc = new KTAViscosityMethodMod(this);
     }
   }
 

src/test/java/neqsim/physicalproperties/methods/commonphasephysicalproperties/viscosity/KTAViscosityMethodModTest.java

@@ -0,0 +1,25 @@
+package neqsim.physicalproperties.methods.commonphasephysicalproperties.viscosity;
+
+import static org.junit.jupiter.api.Assertions.assertEquals;
+import org.junit.jupiter.api.Test;
+import neqsim.thermodynamicoperations.ThermodynamicOperations;
+
+public class KTAViscosityMethodModTest {
+  static neqsim.thermo.system.SystemInterface testSystem = null;
+    @Test
+    void testCalcViscosity() {
+      double T = 273.15;
+      double P = 20.265; //Pressure in MPa
+      testSystem = new neqsim.thermo.system.SystemSrkEos(T, P);
+      testSystem.addComponent("helium", 1.0);
+      //testSystem.addComponent("hydrogen", 0.6);
+      testSystem.setMixingRule(2);
+      ThermodynamicOperations testOps = new ThermodynamicOperations(testSystem);
+      testOps.TPflash();
+      testSystem.getPhase("gas").getPhysicalProperties().setViscosityModel("KTA_mod");
+      testSystem.initProperties();
+      assertEquals(18.878e-6, testSystem.getPhase(0).getPhysicalProperties().getViscosity(), 1e-9);
+      //double viscosity = testSystem.getPhase(0).getPhysicalProperties().getViscosity() * Math.pow(10, 6);
+      //System.out.println("Viscosity: " + viscosity);
+    }
+}

src/test/java/neqsim/physicalproperties/methods/commonphasephysicalproperties/viscosity/KTAViscosityMethodTest.java

@@ -0,0 +1,26 @@
+package neqsim.physicalproperties.methods.commonphasephysicalproperties.viscosity;
+
+import static org.junit.jupiter.api.Assertions.assertEquals;
+import org.junit.jupiter.api.Test;
+import neqsim.thermodynamicoperations.ThermodynamicOperations;
+
+public class KTAViscosityMethodTest {
+  static neqsim.thermo.system.SystemInterface testSystem = null;
+
+  @Test
+  void testCalcViscosity() {
+    double T = 273.15;
+    double P = 20.265; //Pressure in MPa
+    testSystem = new neqsim.thermo.system.SystemSrkEos(T, P);
+    testSystem.addComponent("helium", 1.0);
+    //testSystem.addComponent("hydrogen", 0.6);
+    testSystem.setMixingRule(2);
+    ThermodynamicOperations testOps = new ThermodynamicOperations(testSystem);
+    testOps.TPflash();
+    testSystem.getPhase("gas").getPhysicalProperties().setViscosityModel("KTA");
+    testSystem.initProperties();
+    assertEquals(18.647e-6, testSystem.getPhase(0).getPhysicalProperties().getViscosity(), 1e-9);
+    //double viscosity = testSystem.getPhase(0).getPhysicalProperties().getViscosity() * Math.pow(10, 6);
+    //System.out.println("Viscosity: " + viscosity);
+  }
+}