From 7c19715f5b11ca3c32fad344782bdf00b10d2b88 Mon Sep 17 00:00:00 2001 From: Even Solbraa <41290109+EvenSol@users.noreply.github.com> Date: Sun, 26 Nov 2023 22:28:47 +0100 Subject: [PATCH] Update TPmultiflash.java --- .../flashOps/TPmultiflash.java | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/src/main/java/neqsim/thermodynamicOperations/flashOps/TPmultiflash.java b/src/main/java/neqsim/thermodynamicOperations/flashOps/TPmultiflash.java index ae13aa9c13..5393e0bdc0 100644 --- a/src/main/java/neqsim/thermodynamicOperations/flashOps/TPmultiflash.java +++ b/src/main/java/neqsim/thermodynamicOperations/flashOps/TPmultiflash.java @@ -718,8 +718,8 @@ public void stabilityAnalysis3() { int iter = 0; double errOld = 1.0e100; boolean useaccsubst = true; - int maxsucssubiter = 50; - int maxiter = 100; + int maxsucssubiter = 150; + int maxiter = 200; do { errOld = err; iter++; @@ -1190,9 +1190,9 @@ public void stabilityAnalysis2() { } } } while ((Math.abs(err) > 1e-9 || err > errOld) && iter < 200); - if (iter > 198) { - System.out.println("too many iterations....." + err + " temperature " - + system.getTemperature("C") + " C " + system.getPressure("bara") + " bara"); + if (iter > maxIter-2) { + // System.out.println("too many iterations....." + err + " temperature " + // + system.getTemperature("C") + " C " + system.getPressure("bara") + " bara"); new Exception("to many iterations "); } // logger.info("err: " + err + " ITER " + iter); @@ -1211,7 +1211,7 @@ public void stabilityAnalysis2() { xTrivialCheck0 += Math.abs(x[j][i] - system.getPhase(0).getComponent(i).getx()); xTrivialCheck1 += Math.abs(x[j][i] - system.getPhase(1).getComponent(i).getx()); } - if (iter >= 199) { + if (iter >= maxIter-1) { logger.info("iter > maxiter multiphase stability "); logger.info("error " + Math.abs(err)); logger.info("tm: " + tm[j]);