update flash (#1288)

* update flash

* remove output

* remove output in test

Author: EvenSol

src/main/java/neqsim/thermodynamicoperations/flashops/Flash.java

@@ -158,6 +158,8 @@ public void stabilityAnalysis() throws neqsim.util.exception.IsNaNException,
       }
       iterations = 0;
       fNorm = 1.0e10;
+      boolean acceleration = true;
+      double olderror = 1.0e10;
 
       do {
         iterations++;
@@ -191,7 +193,8 @@ public void stabilityAnalysis() throws neqsim.util.exception.IsNaNException,
               break;
             }
           }
-          if (iterations % 7 == 0 && fNorm < fNormOld && !secondOrderStabilityAnalysis) {
+          if (iterations % 7 == 0 && fNorm < fNormOld && !secondOrderStabilityAnalysis
+              && acceleration) {
             double vec1 = 0.0;
 
             double vec2 = 0.0;
@@ -211,6 +214,9 @@ public void stabilityAnalysis() throws neqsim.util.exception.IsNaNException,
               error[j] += Math.abs((logWi[i] - oldlogw[i]) / oldlogw[i]);
               Wi[j][i] = Math.exp(logWi[i]);
             }
+            if (error[j] > olderror) {
+              acceleration = false;
+            }
           } else {
             // succsessive substitution
             for (int i = 0; i < clonedSystem.getPhases()[0].getNumberOfComponents(); i++) {
@@ -268,14 +274,15 @@ public void stabilityAnalysis() throws neqsim.util.exception.IsNaNException,
         // logger.info("fnorm " + f.norm1() + " err " + error[j] + " iterations " +
         // iterations
         // + " phase " + j);
-      } while ((f.norm1() > 1e-6 && error[j] > 1e-9 && iterations < maxiterations)
+        // logger.info("error " + error[j]);
+        olderror = error[j];
+      } while ((f.norm1() > 1e-3 && error[j] > 1e-3 && iterations < maxiterations)
           || (iterations % 7) == 0 || iterations < 3);
       // (error[j]<oldErr && oldErr<oldOldErr) &&
       // logger.info("err " + error[j]);
       // logger.info("iterations " + iterations);
       // logger.info("f.norm1() " + f.norm1());
       if (iterations >= maxiterations) {
-        logger.error("err staability check " + error[j]);
         throw new neqsim.util.exception.TooManyIterationsException("too many iterations", null,
             maxiterations);
       }

src/test/java/neqsim/process/equipment/valve/ThrottlingValveTest.java

@@ -130,7 +130,6 @@ void testCalcCvLiquidPropnane() {
     valve1.setOutletPressure(19);
     valve1.run();
 
-    stream1.getFluid().prettyPrint();
     assertEquals(135.3602060, valve1.getCv(), 1e-2);
     assertEquals(135.3602060, valve1.getCv("SI"), 1e-2);
     assertEquals(2.46557752, valve1.getCv("US"), 1e-2);

src/test/java/neqsim/process/processmodel/LargeCombinedModelsTest.java

@@ -624,7 +624,7 @@ public void testExpanderProcess() {
     combinedProcess.run();
     combinedProcess.run();
 
-    ((Separator) sepprocessTrain1.getUnit("dew point scrubber")).getFluid().prettyPrint();
+    // ((Separator) sepprocessTrain1.getUnit("dew point scrubber")).getFluid().prettyPrint();
 
     Assertions.assertEquals(8.14523949005678,
         ((Separator) sepprocessTrain1.getUnit("dew point scrubber 2")).getGasOutStream()