---

Ich würde gerne ein neues Abbruchkriterium einbauen:

Die Rechnung soll abbrechen, wenn eine Strömungsgröße auf einen Wert hin konvergiert ist. Beispielsweise soll die Rechnung aufhören, wenn die Werte der Strömungsgröße 100 Iterationen sich nicht mehr als 0.001 unterscheiden.

Ich muss dazu die Dateien convergenceCheck.H und initConvergenceCheck.H verwendet.
Mein Ziel ist es den Massenstrom durch den Auslass zu bestimmen, und diesen von Iterationsschritt zu Iterationsschritt vergleichen zu können.
Wie kann ich denn die einzelnen Iterationschritte aufrufen und diese miteinander vergleichen?

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

So, ich poste mal was ich bisher geschafft habe:

Der Löser berechnet mir den Massenstrom durch den Auslass.
Ich kann in der fvSolution Datei Grenzwerte definieren. Wenn der Massenstrom innerhalb dieser Ober- und Untergrenze liegt, stoppt die Rechnung, Konvergenz ist erreicht. Anbei mein Code.

createFields.H

Code:

---

scalar flux = 0;
mesh.readUpdate();
label patchi = mesh.boundaryMesh().findPatchID("outlet");

---

my_rhoSimpleFoam.C

Code:

---

flux = sum(phi.boundaryField()[patchi]);

---

initConvergenceCheck.H

Code:

---

scalar convergenceCriterion1 = 1000;
scalar convergenceCriterion2 = 1000;

simple.readIfPresent("convergence1", convergenceCriterion1);
simple.readIfPresent("convergence2", convergenceCriterion2);

---

convergenceCheck.H

Code:

---

  if (flux < convergenceCriterion1 && flux > convergenceCriterion2)
{
    Info<< "reached convergence criterion: " << convergenceCriterion2 << " < flux < " << convergenceCriterion1 << endl;
    runTime.writeAndEnd();
    Info<< "latestTime = " << runTime.timeName() << endl;
}

---

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP