Tune FF sonicLiquidFoam to converge

precice · Jul 29, 2020 · 91160d7 · 91160d7
1 parent c0e61f5
commit 91160d7
Show file tree

Hide file tree

Showing 5 changed files with 37 additions and 26 deletions.
diff --git a/tutorials/FF/pipe-pipe/sonicLiquidFoam-sonicLiquidFoam/Fluid1/system/fvSchemes b/tutorials/FF/pipe-pipe/sonicLiquidFoam-sonicLiquidFoam/Fluid1/system/fvSchemes
@@ -17,7 +17,8 @@ FoamFile
 
 ddtSchemes
 {
-    default         Euler;
+    //default         Euler;
+    default         CrankNicolson 0.45;
 }
 
 gradSchemes
@@ -48,4 +49,4 @@ snGradSchemes
 }
 
 
-// ************************************************************************* //
+// ************************************************************************* //
diff --git a/tutorials/FF/pipe-pipe/sonicLiquidFoam-sonicLiquidFoam/Fluid1/system/fvSolution b/tutorials/FF/pipe-pipe/sonicLiquidFoam-sonicLiquidFoam/Fluid1/system/fvSolution
@@ -19,34 +19,37 @@ solvers
 {
     "p.*"
     {
-        solver          smoothSolver;
-        smoother        symGaussSeidel;
+        solver          PBiCG;
+        preconditioner  DILU
+        //smoother        symGaussSeidel;
         tolerance       1e-06;
-        relTol          0;
+        relTol          0.01;
     }
 
     "U.*"
     {
         $p;
+        solver          PBiCG;
+        preconditioner  DILU;
         tolerance       1e-05;
-        relTol          0;
+        relTol          0.1;
     }
 
     "rho.*"
     {
         solver          PCG;
         preconditioner  DIC;
         tolerance       1e-05;
-        relTol          0;
+        relTol          0.1;
     }
 }
 
 PIMPLE
 {
-    nOuterCorrectors 2;
-    nCorrectors      1;
+    nOuterCorrectors 4;
+    nCorrectors      2;
     nNonOrthogonalCorrectors 0;
 }
 
 
-// ************************************************************************* //
+// ************************************************************************* //
diff --git a/tutorials/FF/pipe-pipe/sonicLiquidFoam-sonicLiquidFoam/Fluid2/system/fvSchemes b/tutorials/FF/pipe-pipe/sonicLiquidFoam-sonicLiquidFoam/Fluid2/system/fvSchemes
@@ -17,7 +17,8 @@ FoamFile
 
 ddtSchemes
 {
-    default         Euler;
+    //default         Euler;
+    default         CrankNicolson 0.45;
 }
 
 gradSchemes
@@ -48,4 +49,4 @@ snGradSchemes
 }
 
 
-// ************************************************************************* //
+// ************************************************************************* //
diff --git a/tutorials/FF/pipe-pipe/sonicLiquidFoam-sonicLiquidFoam/Fluid2/system/fvSolution b/tutorials/FF/pipe-pipe/sonicLiquidFoam-sonicLiquidFoam/Fluid2/system/fvSolution
@@ -19,34 +19,37 @@ solvers
 {
     "p.*"
     {
-        solver          smoothSolver;
-        smoother        symGaussSeidel;
+        solver          PBiCG;
+        preconditioner  DILU
+        //smoother        symGaussSeidel;
         tolerance       1e-06;
-        relTol          0;
+        relTol          0.01;
     }
 
     "U.*"
     {
         $p;
+        solver          PBiCG;
+        preconditioner  DILU;
         tolerance       1e-05;
-        relTol          0;
+        relTol          0.1;
     }
 
     "rho.*"
     {
         solver          PCG;
         preconditioner  DIC;
         tolerance       1e-05;
-        relTol          0;
+        relTol          0.1;
     }
 }
 
 PIMPLE
 {
-    nOuterCorrectors 2;
-    nCorrectors      1;
+    nOuterCorrectors 4;
+    nCorrectors      2;
     nNonOrthogonalCorrectors 0;
 }
 
 
-// ************************************************************************* //
+// ************************************************************************* //
diff --git a/tutorials/FF/pipe-pipe/sonicLiquidFoam-sonicLiquidFoam/precice-config.xml b/tutorials/FF/pipe-pipe/sonicLiquidFoam-sonicLiquidFoam/precice-config.xml
@@ -50,7 +50,7 @@
 
     <m2n:sockets from="Fluid1" to="Fluid2"/>
 
-    <coupling-scheme:serial-implicit>
+    <coupling-scheme:parallel-implicit>
         <time-window-size value="0.01"/>
         <max-time value="1.0"/>
         <participants first="Fluid1" second="Fluid2"/>
@@ -60,18 +60,21 @@
         <!--<exchange data="VelocityGradient" mesh="Fluid2-Mesh" from="Fluid2" to="Fluid1"/>-->
 
         <max-iterations value="100"/>
-        <relative-convergence-measure limit="1.0e-6" data="Pressure" mesh="Fluid2-Mesh"/>
+        <relative-convergence-measure limit="1.0e-3" data="Pressure" mesh="Fluid2-Mesh"/>
         <!--<relative-convergence-measure limit="1.0e-6" data="PressureGradient" mesh="Fluid2-Mesh"/>-->
-        <relative-convergence-measure limit="1.0e-6" data="Velocity" mesh="Fluid1-Mesh"/>
+        <relative-convergence-measure limit="1.0e-3" data="Velocity" mesh="Fluid1-Mesh"/>
 
         <acceleration:IQN-ILS>
             <data mesh="Fluid2-Mesh" name="Pressure" />
-            <initial-relaxation value="0.01" />
+            <data mesh="Fluid1-Mesh" name="Velocity" />
+            <preconditioner type="residual-sum" />
+            <initial-relaxation value="0.1" />
             <max-used-iterations value="80" />
             <time-windows-reused value="10" />
-            <filter type="QR1" limit="1e-8" />
+            <filter type="QR1" limit="1e-6" />
         </acceleration:IQN-ILS>
-    </coupling-scheme:serial-implicit>
+
+    </coupling-scheme:parallel-implicit>
 
     </solver-interface>
 </precice-configuration>