added tutorial for foam-extend

quickstart/fluid-openfoam/0/U

@@ -22,11 +22,13 @@ boundaryField
     top
     {
         type            noSlip;
+        value           uniform (0 0 0);
     }
 
     bottom
     {
         type            noSlip;
+        value           uniform (0 0 0);
     }
 
     inlet

quickstart/fluid-openfoam/constant/dynamicMeshDict

@@ -14,6 +14,8 @@ motionSolverLibs ("libfvMotionSolvers.so");
 solver      displacementLaplacian;
 // OpenFOAM9 or newer: rename "solver" to "motionSolver"
 
+diffusivity uniform;
+
 displacementLaplacianCoeffs {
     diffusivity quadratic inverseDistance (flap);
 }

quickstart/fluid-openfoam/system/controlDict

@@ -3,13 +3,13 @@ FoamFile
     version     2.0;
     format      ascii;
     class       dictionary;
+    location    "system";
     object      controlDict;
 }
 
 application     pimpleFoam;       // latest OpenFOAM
 // application     pimpleDyMFoam; // OpenFOAM v1712, OpenFOAM 5.x, or older
 
-
 startFrom           startTime;
 
 startTime           0;
@@ -36,14 +36,23 @@ timeFormat          general;
 
 timePrecision       8;
 
-libs ("libpreciceAdapterFunctionObject.so");
+libs
+(
+    "libpreciceAdapterFunctionObject.so"
+);
 functions
 {
     forces
     {
         type                forces;
         libs                ( "libforces.so" );
+        functionObjectLibs  ( "libforces.so" );
+        outputControl       timeStep;
+        outputInterval      1;
         patches             (flap);
+        pName               p;
+        UName               U;
+        rhoName             rhoInf;
         rho                 rhoInf;
         log                 true;
         rhoInf              10;

quickstart/fluid-openfoam/system/fvSchemes

@@ -21,6 +21,7 @@ divSchemes
 {
     default         none;
     div(phi,U)      Gauss linearUpwind grad(U);
+    div((nuEff*dev(T(grad(U))))) Gauss linear;
     div((nuEff*dev2(T(grad(U))))) Gauss linear;
 }
 

quickstart/fluid-openfoam/system/fvSolution

@@ -14,6 +14,9 @@ solvers
         tolerance        1e-6;
         relTol           1e-4;
         smoother         DICGaussSeidel;
+        agglomerator    faceAreaPair;
+        nCellsInCoarsestLevel 10;
+        mergeLevels 1;
     }
 
     pFinal
@@ -78,3 +81,16 @@ potentialFlow
 {
     nNonOrthogonalCorrectors 1;
 }
+
+relaxationFactors
+{
+    U          0.7;
+    UFinal     1;
+}
+
+fieldBounds
+{
+    p -1e5 1e5;
+    U 100;
+}
+

quickstart/fluid-openfoam/system/preciceDict

@@ -36,4 +36,5 @@ interfaces
 FSI
 {
   rho rho [1 -3 0 0 0 0 0] 1000;
+  nu nu [0 2 -1 0 0 0 0] 0.001;
 }

tools/run-foam-extend.sh

@@ -0,0 +1,41 @@
+#!/usr/bin/env sh
+set -e # Not setting -u as it gets triggered by the OpenFOAM RunFunctions
+
+# Prepare an (intentionally empty) .foam file for the ParaView OpenFOAM reader
+CASENAME="$(pwd | xargs basename)"
+touch "$CASENAME.foam"
+
+# Modify code for foam-extend
+echo "modifying everything now"
+sed -i "s/noSlip;/noSlipWall;/g" 0/U
+sed -i "s,application     pimpleFoam;,//application     pimpleFoam;,g" system/controlDict
+sed -i "s,// application     pimpleDyMFoam;,application     pimpleDyMFoam;,g" system/controlDict
+sed -i '41i\ \ \ \ "liblduSolvers.so"' system/controlDict
+sed -i '41i\ \ \ \ "libforces.so"' system/controlDict
+sed -i "s,writeCompression    off,writeCompression    uncompressed,g" system/controlDict
+
+sed -i "s/libfvMotionSolvers\./libfvMotionSolver\./g" constant/dynamicMeshDict
+
+# OpenFOAM run functions: getApplication, getNumberOfProcessors
+# shellcheck disable=SC1090 # This is an OpenFOAM file which we don't need to check
+. "${WM_PROJECT_DIR}/bin/tools/RunFunctions"
+solver=$(getApplication | cut -f 1 -d " " | sed '\~//~d')
+if [ "${1:-}" = "-parallel" ]; then
+    procs=$(getNumberOfProcessors)
+    decomposePar -force
+    mpirun -np "${procs}" "${solver}" -parallel
+    reconstructPar
+else
+    ${solver}
+fi
+
+# Reverse code for OpenFOAM
+#rm -rf constant/polyMesh
+#sed -i "s/noSlipWall;/noSlip;/g" 0/U
+#sed -i "s,application     pimpleDyMFoam;,// application     pimpleDyMFoam;,g" system/controlDict
+#sed -i "s,//application     pimpleFoam;,application     pimpleFoam;,g" system/controlDict
+#sed -i '/   "liblduSolvers.so"/d' system/controlDict
+#sed -i '/   "libforces.so/d' system/controlDict
+#sed -i "s,writeCompression    uncompressed,writeCompression    off,g" system/controlDict
+#
+#sed -i "s/libfvMotionSolver\./libfvMotionSolvers\./g" constant/dynamicMeshDict