Add perpendicular-flap case solid-solids4foam and solid-openfoam

perpendicular-flap/README.md

@@ -37,6 +37,8 @@ Solid participant:
 
 * DUNE. For more information, have a look at the [experimental DUNE adapter](https://github.com/precice/dune-adapter) and send us your feedback.
 
+* solids4foam. Like for CalculuX, the geometrically linear solver is used by default. For more information, see the [solids4foam documentation](https://bitbucket.org/philip_cardiff/solids4foam-release/src/master/documentation/overview.md)
+
 ## Running the Simulation
 
 All listed solvers can be used in order to run the simulation. OpenFOAM can be executed in parallel using `run.sh -parallel`. The default setting uses 4 MPI ranks. Open two separate terminals and start the desired fluid and solid participant by calling the respective run script `run.sh` located in the participant directory. For example:

perpendicular-flap/solid-solids4foam/0/D

@@ -0,0 +1,33 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      D;
+}
+
+dimensions      [0 1 0 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    interface
+    {
+        type            solidForce;
+        forceField      solidForce;
+        value           uniform (0 0 0);
+    }
+    bottom
+    {
+        type            fixedDisplacement;
+        value           uniform (0 0 0);
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+

perpendicular-flap/solid-solids4foam/0/solidForce

@@ -0,0 +1,32 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      solidForce;
+}
+
+dimensions      [1 2 2 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    interface
+    {
+        type            calculated;
+        value           uniform (0 0 0);
+    }
+    bottom
+    {
+        type            calculated;
+        value           uniform (0 0 0);
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+

perpendicular-flap/solid-solids4foam/clean.sh

@@ -0,0 +1,6 @@
+#!/bin/sh
+set -e -u
+
+. ../../tools/cleaning-tools.sh
+
+clean_openfoam .

perpendicular-flap/solid-solids4foam/constant/dynamicMeshDict

@@ -0,0 +1,10 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      dynamicMeshDict;
+}
+
+dynamicFvMesh   staticFvMesh;
+

perpendicular-flap/solid-solids4foam/constant/g

@@ -0,0 +1,12 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       uniformDimensionedVectorField;
+    location    "constant";
+    object      g;
+}
+
+dimensions      [0 1 -2 0 0 0 0];
+value           ( 0 0 0 );
+

perpendicular-flap/solid-solids4foam/constant/mechanicalProperties

@@ -0,0 +1,21 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      mechanicalProperties;
+}
+
+planeStress no;
+
+mechanical
+(
+    rubber
+    {
+        type        linearElastic;
+        rho         rho [1 -3 0 0 0 0 0] 3000;
+        E           E [1 -1 -2 0 0 0 0] 4e6;
+        nu          nu [0 0 0 0 0 0 0] 0.3;
+    }
+);
+

perpendicular-flap/solid-solids4foam/constant/physicsProperties

@@ -0,0 +1,12 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      physicsProperties;
+}
+
+type    solid;
+
+

perpendicular-flap/solid-solids4foam/constant/solidProperties

@@ -0,0 +1,25 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      solidProperties;
+}
+
+solidModel    linearGeometryTotalDisplacement;
+
+linearGeometryTotalDisplacementCoeffs
+{
+    // Maximum number of momentum correctors
+    nCorrectors     1000;
+
+    // Solution tolerance for displacement
+    solutionTolerance 1e-08;
+
+    // Alternative solution tolerance for displacement
+    alternativeTolerance 1e-08;
+
+    // Write frequency for the residuals
+    infoFrequency   100;
+}
+

perpendicular-flap/solid-solids4foam/run.sh

@@ -0,0 +1,14 @@
+#!/bin/bash
+
+# Currently, the case has only been checked with OpenFOAM-v2012
+if [[ "${WM_PROJECT}" != "OpenFOAM" || "${WM_PROJECT_VERSION}" != "v2012" ]]
+then
+    echo; echo "This case currently only runs in OpenFOAM-v2012"; echo
+    exit 0
+fi
+
+blockMesh
+touch fluid-openfoam.foam
+
+../../tools/run-openfoam.sh "$@"
+. ../../tools/openfoam-remove-empty-dirs.sh && openfoam_remove_empty_dirs

perpendicular-flap/solid-solids4foam/system/blockMeshDict

@@ -0,0 +1,55 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      blockMeshDict;
+}
+
+convertToMeters 0.01;
+
+vertices
+(
+    (-5 0 0)
+    (5 0 0)
+    (5 100 0)
+    (-5 100 0)
+
+    (-5 0 100)
+    (5 0 100)
+    (5 100 100)
+    (-5 100 100)
+);
+
+blocks
+(
+    hex (0 1 2 3 4 5 6 7) (6 15 1) simpleGrading (1 1 1)
+);
+
+edges
+(
+);
+
+patches
+(
+    patch interface
+    (
+        (4 7 3 0)
+        (7 6 2 3)
+        (1 2 6 5)
+    )
+    patch bottom
+    (
+        (0 1 5 4)
+    )
+    empty frontAndBack
+    (
+        (3 2 1 0)
+        (4 5 6 7)
+    )
+);
+
+mergePatchPairs
+(
+);
+

perpendicular-flap/solid-solids4foam/system/controlDict

@@ -0,0 +1,48 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      controlDict;
+}
+
+application     solids4Foam;
+
+startFrom       startTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         10;
+
+deltaT          0.01;
+
+writeControl    adjustableRunTime;
+
+writeInterval   0.1;
+
+purgeWrite      0;
+
+writeFormat     ascii;
+
+writePrecision  6;
+
+writeCompression off;
+
+timeFormat      general;
+
+timePrecision   6;
+
+functions
+{
+    preCICE_Adapter
+    {
+        type preciceAdapterFunctionObject;
+        libs ("libpreciceAdapterFunctionObject.so");
+    }
+}
+
+
+

perpendicular-flap/solid-solids4foam/system/decomposeParDict

@@ -0,0 +1,10 @@
+numberOfSubdomains 4;
+
+method          simple;
+
+simpleCoeffs
+{
+    n               (1 4 1);
+    delta           0.001;
+}
+

perpendicular-flap/solid-solids4foam/system/fvSchemes

@@ -0,0 +1,51 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      fvSchemes;
+}
+
+d2dt2Schemes
+{
+    default            backward;
+}
+
+ddtSchemes
+{
+    default            backward;
+}
+
+gradSchemes
+{
+    default            extendedLeastSquares 0;
+}
+
+divSchemes
+{
+    default            Gauss linear;
+}
+
+laplacianSchemes
+{
+    default            none;
+    laplacian(DD,D)    Gauss linear corrected;
+    laplacian(DDD,DD)  Gauss linear corrected;
+}
+
+snGradSchemes
+{
+    default            none;
+    snGrad(D)          corrected;
+    snGrad(DD)         corrected;
+}
+
+interpolationSchemes
+{
+    default            none;
+    interpolate(impK)  linear;
+    interpolate(grad(D)) linear;
+    interpolate(grad(DD)) linear;
+    interpolate(sigma0) linear;
+}
+

perpendicular-flap/solid-solids4foam/system/fvSolution

@@ -0,0 +1,38 @@
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSolution;
+}
+
+solvers
+{
+    "D|DD"
+    {
+        solver          PCG;
+        preconditioner  FDIC;
+        tolerance       1e-09;
+        relTol          0.1;
+    }
+}
+
+relaxationFactors
+{
+    // Under-relaxating the D equation by 0.99-0.9999 can improve convergence in
+    // some cases, in particular when there are solidContact boundaries
+    equations
+    {
+        //D   0.999;
+    }
+
+    // Under-relaxating the D field by 0.1-0.9 can improve convergence in some
+    // cases
+    fields
+    {
+        //D      0.9;
+    }
+}
+
+