fix solid calcs (#976)

* fix solid flash * update * update * update * add chack for component in temp database * remove output
equinor · Apr 1, 2024 · 7204d95 · 7204d95
1 parent ae2bc5f
commit 7204d95
Show file tree

Hide file tree

Showing 6 changed files with 217 additions and 124 deletions.
diff --git a/src/main/java/neqsim/thermo/component/ComponentSolid.java b/src/main/java/neqsim/thermo/component/ComponentSolid.java
@@ -8,6 +8,7 @@
 
 import neqsim.thermo.phase.PhaseInterface;
 import neqsim.thermo.phase.PhaseType;
+import neqsim.util.database.NeqSimDataBase;
 
 /**
  * <p>
@@ -240,7 +241,13 @@ public void setSolidRefFluidPhase(PhaseInterface phase) {
       refPhase.setTemperature(273.0);
       refPhase.setPressure(1.0);
       try {
-        refPhase.addComponent(componentName, 10.0, 10.0, 0);
+        if (NeqSimDataBase.hasComponent(componentName)
+            || NeqSimDataBase.hasTempComponent(componentName)) {
+          refPhase.addComponent(componentName, 10.0, 10.0, 0);
+        } else {
+          refPhase.addComponent("methane", 10.0, 10.0, 0);
+          refPhase.getComponent("methane").setComponentName(componentName);
+        }
       } catch (Exception ex) {
         logger.error("error occured in setSolidRefFluidPhase ", ex);
         refPhase.addComponent("methane", 10.0, 10.0, 0);

diff --git a/src/main/java/neqsim/thermo/phase/PhasePureComponentSolid.java b/src/main/java/neqsim/thermo/phase/PhasePureComponentSolid.java
@@ -24,7 +24,6 @@ public class PhasePureComponentSolid extends PhaseSolid {
    */
   public PhasePureComponentSolid() {
     super();
-    setType(PhaseType.SOLID);
   }
 
   /** {@inheritDoc} */

diff --git a/src/main/java/neqsim/thermo/phase/PhaseSolid.java b/src/main/java/neqsim/thermo/phase/PhaseSolid.java
@@ -27,6 +27,7 @@ public abstract class PhaseSolid extends PhaseSrkEos {
   public PhaseSolid() {
     super();
     setType(PhaseType.SOLID);
+    calcMolarVolume = false;
   }
 
   /** {@inheritDoc} */
@@ -48,6 +49,7 @@ public void init(double totalNumberOfMoles, int numberOfComponents, int initType
       double beta) {
     try {
       super.init(totalNumberOfMoles, numberOfComponents, initType, pt, beta);
+      getDensityTemp();
     } catch (Exception ex) {
       logger.error(ex.getMessage());
     }