Immersed boundary tutorial: two interfering IB patches

tutorials/immersedBoundary/twoIbPatches/.gitignore

@@ -0,0 +1,52 @@
+# git-ls-files --others --exclude-from=.git/info/exclude
+# Lines that start with '#' are comments.
+
+# editor and misc backup files - anywhere
+*~
+.*~
+*.bak
+*.bak[0-9][0-9]
+*.orig
+*.orig[0-9][0-9]
+\#*\#
+
+# file-browser settings - anywhere
+.directory
+
+# CVS recovered versions - anywhere
+.#*
+
+# SVN directories - anywhere
+
+.svn/
+
+# OpenFOAM results
+
+[1-9]*/
+!/0/
+processor*
+*/polyMesh/*
+!*/polyMesh/blockMeshDict
+cellToRegion*
+log*
+
+# packages - anywhere
+
+*.tar.bz2
+*.tar.gz
+*.tar
+*.tgz
+*.gtgz
+
+# Pictures and movies
+
+*.png
+*.jpg
+*.jpeg
+*.bmp
+*.png
+*.avi
+*.mp4
+*.mpg
+
+#end-of-file

tutorials/immersedBoundary/twoIbPatches/0_org/U

@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volVectorField;
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (1 0 0);
+
+boundaryField
+{
+    ibCylinder
+    {
+        type movingImmersedBoundaryVelocity;
+        patchType immersedBoundary;
+        setDeadValue   yes;
+        deadValue      (0 0 0);
+        value uniform (0 0 0);
+    }
+    ibCylinder2
+    {
+        type movingImmersedBoundaryVelocity;
+        patchType immersedBoundary;
+        setDeadValue   yes;
+        deadValue      (0 0 0);
+        value uniform (0 0 0);
+    }
+    in
+    {
+        type fixedValue;
+        value uniform (1 0 0);
+    }
+    out
+    {
+        type inletOutlet;
+        inletValue uniform (0 0 0);
+        value uniform (1 0 0);
+    }
+    top
+    {
+        type fixedValue;
+        value uniform (0 0 0);
+    }
+    bottom
+    {
+        type fixedValue;
+        value uniform (0 0 0);
+    }
+    frontAndBack
+    {
+        type empty;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/0_org/epsilon

@@ -0,0 +1,74 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      omega;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -3 0 0 0 0];
+
+internalField   uniform 2e-4;
+
+boundaryField
+{
+    ibCylinder
+    {
+        type            immersedBoundaryEpsilonWallFunction;
+        patchType       immersedBoundary;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        setDeadValue    true;
+        deadValue       2e-4;
+
+        value           uniform 2e-4;
+    }
+    ibCylinder2
+    {
+        type            immersedBoundaryEpsilonWallFunction;
+        patchType       immersedBoundary;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+        setDeadValue    true;
+        deadValue       2e-4;
+
+        value           uniform 2e-4;
+    }
+    top
+    {
+        type            zeroGradient;
+    }
+    in
+    {
+        type            fixedValue;
+        value           uniform 2e-4;
+    }
+    out
+    {
+        type            inletOutlet;
+        inletValue      uniform 2e-4;
+        value           uniform 2e-4;
+    }
+    bottom
+    {
+        type            zeroGradient;
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/0_org/k

@@ -0,0 +1,66 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      k;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 1.5e-4;
+
+boundaryField
+{
+    ibCylinder
+    {
+        type            immersedBoundaryKqRWallFunction;
+        patchType       immersedBoundary;
+        setDeadValue    yes;
+        deadValue       1.5e-4;
+        value           uniform 1.5e-4;
+    }
+    ibCylinder2
+    {
+        type            immersedBoundaryKqRWallFunction;
+        patchType       immersedBoundary;
+        setDeadValue    yes;
+        deadValue       1.5e-4;
+        value           uniform 1.5e-4;
+    }
+    top
+    {
+        type            zeroGradient;
+    }
+    in
+    {
+        type            fixedValue;
+        value           uniform 1.5e-4;
+    }
+    out
+    {
+        type            inletOutlet;
+        inletValue      uniform 1.5e-4;
+        value           uniform 1.5e-4;
+    }
+    bottom
+    {
+        type            zeroGradient;
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/0_org/nut

@@ -0,0 +1,67 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      nut;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    ibCylinder
+    {
+        type            immersedBoundaryNutWallFunction;
+        patchType       immersedBoundary;
+        setDeadValue    true;
+        deadValue       1e-5;
+        value           uniform 0;
+    }
+    ibCylinder2
+    {
+        type            immersedBoundaryNutWallFunction;
+        patchType       immersedBoundary;
+        setDeadValue    true;
+        deadValue       1e-5;
+        value           uniform 0;
+    }
+    top
+    {
+        type            zeroGradient;
+        value           uniform 0;
+    }
+    in
+    {
+        type            zeroGradient;
+        value           uniform 0;
+    }
+    out
+    {
+        type            zeroGradient;
+        value           uniform 0;
+    }
+    bottom
+    {
+        type            zeroGradient;
+        value           uniform 0;
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/0_org/omega

@@ -0,0 +1,70 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      omega;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 -1 0 0 0 0];
+
+internalField   uniform 500;
+
+boundaryField
+{
+    ibCylinder
+    {
+        type            immersedBoundaryOmegaWallFunction;
+        patchType       immersedBoundary;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+
+        value           uniform 500;
+    }
+    ibCylinder2
+    {
+        type            immersedBoundaryOmegaWallFunction;
+        patchType       immersedBoundary;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+
+        value           uniform 500;
+    }
+    top
+    {
+        type            zeroGradient;
+    }
+    in
+    {
+        type            fixedValue;
+        value           uniform 500;
+    }
+    out
+    {
+        type            inletOutlet;
+        inletValue      uniform 500;
+        value           uniform 500;
+    }
+    bottom
+    {
+        type            zeroGradient;
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/0_org/p

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions [0 2 -2 0 0 0 0];
+
+internalField uniform 0;
+
+boundaryField
+{
+    ibCylinder
+    {
+        type mixedIb;
+        patchType immersedBoundary;
+        triValue uniform 0;
+        triGradient  uniform 0;
+        triValueFraction uniform 0;
+        setDeadValue   yes;
+        deadValue      0;
+        value uniform 0;
+    }
+    ibCylinder2
+    {
+        type mixedIb;
+        patchType immersedBoundary;
+        triValue uniform 0;
+        triGradient  uniform 0;
+        triValueFraction uniform 0;
+        setDeadValue   yes;
+        deadValue      0;
+        value uniform 0;
+    }
+    in
+    {
+        type zeroGradient;
+    }
+    out
+    {
+        type fixedValue;
+        value uniform 0;
+    }
+    top
+    {
+        type zeroGradient;
+    }
+    bottom
+    {
+        type zeroGradient;
+    }
+    frontAndBack
+    {
+        type empty;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/Allclean

@@ -0,0 +1,10 @@
+#!/bin/sh
+
+# Source tutorial clean functions
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase
+\rm -f *.dat
+\rm -rf 0
+\rm -rf constant/polyMesh
+\rm -rf postProcessing

tutorials/immersedBoundary/twoIbPatches/Allrun

@@ -0,0 +1,16 @@
+#!/bin/sh
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+# Get application name
+application="pimpleDyMIbFoam"
+
+runApplication blockMesh -dict save/blockMeshDict
+cp save/boundary constant/polyMesh/
+
+mkdir 0
+\cp -f 0_org/* 0/
+
+runApplication potentialFoam
+runApplication writeIbMasks
+runApplication $application

tutorials/immersedBoundary/twoIbPatches/constant/RASProperties

@@ -0,0 +1,25 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      RASProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+RASModel            kEpsilon;
+// RASModel            kOmegaSST;
+
+turbulence          on;
+
+printCoeffs         on;
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/constant/dynamicMeshDict

@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      dynamicMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+// dynamicFvMesh   staticFvMesh;
+dynamicFvMesh   immersedBoundarySolidBodyMotionFvMesh;
+
+immersedBoundarySolidBodyMotionFvMeshCoeffs
+{
+    motionFunctions
+    (
+
+        ibCylinder
+        {
+//             solidBodyMotionFunction     translation;
+//             translationCoeffs
+//             {
+//                 velocity (0.1 0 0);
+//             }
+
+            solidBodyMotionFunction     linearOscillation;
+            linearOscillationCoeffs
+            {
+                amplitude    (0.5 0 0);
+                period       2.5;
+            }
+        }
+
+        ibCylinder2
+        {
+            solidBodyMotionFunction     linearOscillation;
+            linearOscillationCoeffs
+            {
+                amplitude    (0 0.05 0);
+                period       2.5;
+            }
+        }
+
+
+    );
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/constant/transportProperties

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      transportProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+transportModel  Newtonian;
+nu              nu [0 2 -1 0 0 0 0] 1e-6;
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder.ftr

@@ -0,0 +1,205 @@
+1(
+cylinder
+patch)
+
+96
+(
+(0.176776 -0.176776 -0.05)
+(0.176776 -0.176776 0.05)
+(0.176776 0.176776 -0.05)
+(0.176776 0.176776 0.05)
+(-0.176776 0.176776 -0.05)
+(-0.176776 0.176776 0.05)
+(-0.176776 -0.176776 -0.05)
+(-0.176776 -0.176776 0.05)
+(0.241481 -0.064705 -0.05)
+(0.241481 -0.064705 0.05)
+(0.241481 0.064705 -0.05)
+(0.241481 0.064705 0.05)
+(0.064705 0.241481 -0.05)
+(0.064705 0.241481 0.05)
+(-0.064705 0.241481 -0.05)
+(-0.064705 0.241481 0.05)
+(-0.241481 0.064705 -0.05)
+(-0.241481 0.064705 0.05)
+(-0.241481 -0.064705 -0.05)
+(-0.241481 -0.064705 0.05)
+(0.064705 -0.241481 0.05)
+(0.064705 -0.241481 -0.05)
+(-0.064705 -0.241481 0.05)
+(-0.064705 -0.241481 -0.05)
+(0.216506 -0.125 -0.05)
+(0.216506 -0.125 0.05)
+(0.216506 0.125 -0.05)
+(0.216506 0.125 0.05)
+(0.125 0.216506 -0.05)
+(0.125 0.216506 0.05)
+(-0.125 0.216506 -0.05)
+(-0.125 0.216506 0.05)
+(-0.216506 0.125 -0.05)
+(-0.216506 0.125 0.05)
+(-0.216506 -0.125 -0.05)
+(-0.216506 -0.125 0.05)
+(0.125 -0.216506 0.05)
+(0.125 -0.216506 -0.05)
+(-0.125 -0.216506 0.05)
+(-0.125 -0.216506 -0.05)
+(0.247861 -0.0326315 -0.05)
+(0.247861 -0.0326315 0.05)
+(0.247861 0.0326315 -0.05)
+(0.247861 0.0326315 0.05)
+(0.0326315 0.247861 -0.05)
+(0.0326315 0.247861 0.05)
+(-0.0326315 0.247861 -0.05)
+(-0.0326315 0.247861 0.05)
+(-0.247861 0.0326315 -0.05)
+(-0.247861 0.0326315 0.05)
+(-0.247861 -0.0326315 -0.05)
+(-0.247861 -0.0326315 0.05)
+(0.0326315 -0.247861 0.05)
+(0.0326315 -0.247861 -0.05)
+(-0.0326315 -0.247861 0.05)
+(-0.0326315 -0.247861 -0.05)
+(0.25 -8.32665e-17 0.05)
+(0.25 8.32665e-17 -0.05)
+(8.32665e-17 0.25 0.05)
+(-8.32665e-17 0.25 -0.05)
+(-0.25 8.32665e-17 0.05)
+(-0.25 -8.32665e-17 -0.05)
+(8.32665e-17 -0.25 -0.05)
+(-8.32665e-17 -0.25 0.05)
+(0.198338 -0.152191 -0.05)
+(0.198338 -0.152191 0.05)
+(0.198338 0.152191 -0.05)
+(0.198338 0.152191 0.05)
+(0.152191 0.198338 -0.05)
+(0.152191 0.198338 0.05)
+(-0.152191 0.198338 -0.05)
+(-0.152191 0.198338 0.05)
+(-0.198338 0.152191 -0.05)
+(-0.198338 0.152191 0.05)
+(-0.198338 -0.152191 -0.05)
+(-0.198338 -0.152191 0.05)
+(0.152191 -0.198338 0.05)
+(0.152191 -0.198338 -0.05)
+(-0.152191 -0.198338 0.05)
+(-0.152191 -0.198338 -0.05)
+(0.23097 -0.095671 -0.05)
+(0.23097 -0.095671 0.05)
+(0.23097 0.095671 -0.05)
+(0.23097 0.095671 0.05)
+(0.095671 0.23097 -0.05)
+(0.095671 0.23097 0.05)
+(-0.095671 0.23097 -0.05)
+(-0.095671 0.23097 0.05)
+(-0.23097 0.095671 -0.05)
+(-0.23097 0.095671 0.05)
+(-0.23097 -0.095671 -0.05)
+(-0.23097 -0.095671 0.05)
+(0.095671 -0.23097 0.05)
+(0.095671 -0.23097 -0.05)
+(-0.095671 -0.23097 0.05)
+(-0.095671 -0.23097 -0.05)
+)
+
+
+96
+(
+((0 64 65) 0)
+((65 1 0) 0)
+((64 24 25) 0)
+((25 65 64) 0)
+((24 80 81) 0)
+((81 25 24) 0)
+((80 8 9) 0)
+((9 81 80) 0)
+((8 40 41) 0)
+((41 9 8) 0)
+((56 41 40) 0)
+((40 57 56) 0)
+((57 42 43) 0)
+((43 56 57) 0)
+((42 10 11) 0)
+((11 43 42) 0)
+((10 82 83) 0)
+((83 11 10) 0)
+((82 26 27) 0)
+((27 83 82) 0)
+((26 66 67) 0)
+((67 27 26) 0)
+((66 2 3) 0)
+((3 67 66) 0)
+((2 68 69) 0)
+((69 3 2) 0)
+((68 28 29) 0)
+((29 69 68) 0)
+((28 84 85) 0)
+((85 29 28) 0)
+((84 12 13) 0)
+((13 85 84) 0)
+((12 44 45) 0)
+((45 13 12) 0)
+((58 45 44) 0)
+((44 59 58) 0)
+((59 46 47) 0)
+((47 58 59) 0)
+((46 14 15) 0)
+((15 47 46) 0)
+((14 86 87) 0)
+((87 15 14) 0)
+((86 30 31) 0)
+((31 87 86) 0)
+((30 70 71) 0)
+((71 31 30) 0)
+((70 4 5) 0)
+((5 71 70) 0)
+((4 72 73) 0)
+((73 5 4) 0)
+((72 32 33) 0)
+((33 73 72) 0)
+((32 88 89) 0)
+((89 33 32) 0)
+((88 16 17) 0)
+((17 89 88) 0)
+((16 48 49) 0)
+((49 17 16) 0)
+((60 49 48) 0)
+((48 61 60) 0)
+((61 50 51) 0)
+((51 60 61) 0)
+((50 18 19) 0)
+((19 51 50) 0)
+((18 90 91) 0)
+((91 19 18) 0)
+((90 34 35) 0)
+((35 91 90) 0)
+((34 74 75) 0)
+((75 35 34) 0)
+((74 6 7) 0)
+((7 75 74) 0)
+((0 1 76) 0)
+((76 77 0) 0)
+((77 76 36) 0)
+((36 37 77) 0)
+((37 36 92) 0)
+((92 93 37) 0)
+((93 92 20) 0)
+((20 21 93) 0)
+((21 20 52) 0)
+((52 53 21) 0)
+((62 53 52) 0)
+((52 63 62) 0)
+((63 54 55) 0)
+((55 62 63) 0)
+((55 54 22) 0)
+((22 23 55) 0)
+((23 22 94) 0)
+((94 95 23) 0)
+((95 94 38) 0)
+((38 39 95) 0)
+((39 38 78) 0)
+((78 79 39) 0)
+((79 78 7) 0)
+((7 6 79) 0)
+)
+

tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder.stl

@@ -0,0 +1,674 @@
+solid cylinder
+  facet normal 0.751817 -0.659372 0
+    outer loop
+       vertex 0.176776 -0.176776 -0.05
+       vertex 0.198338 -0.152191 -0.05
+       vertex 0.198338 -0.152191 0.05
+    endloop
+  endfacet
+  facet normal 0.751817 -0.659372 0
+    outer loop
+       vertex 0.198338 -0.152191 0.05
+       vertex 0.176776 -0.176776 0.05
+       vertex 0.176776 -0.176776 -0.05
+    endloop
+  endfacet
+  facet normal 0.831476 -0.555561 0
+    outer loop
+       vertex 0.198338 -0.152191 -0.05
+       vertex 0.216506 -0.125 -0.05
+       vertex 0.216506 -0.125 0.05
+    endloop
+  endfacet
+  facet normal 0.831476 -0.555561 0
+    outer loop
+       vertex 0.216506 -0.125 0.05
+       vertex 0.198338 -0.152191 0.05
+       vertex 0.198338 -0.152191 -0.05
+    endloop
+  endfacet
+  facet normal 0.896866 -0.442302 0
+    outer loop
+       vertex 0.216506 -0.125 -0.05
+       vertex 0.23097 -0.095671 -0.05
+       vertex 0.23097 -0.095671 0.05
+    endloop
+  endfacet
+  facet normal 0.896866 -0.442302 0
+    outer loop
+       vertex 0.23097 -0.095671 0.05
+       vertex 0.216506 -0.125 0.05
+       vertex 0.216506 -0.125 -0.05
+    endloop
+  endfacet
+  facet normal 0.946935 -0.321425 0
+    outer loop
+       vertex 0.23097 -0.095671 -0.05
+       vertex 0.241481 -0.064705 -0.05
+       vertex 0.241481 -0.064705 0.05
+    endloop
+  endfacet
+  facet normal 0.946935 -0.321425 0
+    outer loop
+       vertex 0.241481 -0.064705 0.05
+       vertex 0.23097 -0.095671 0.05
+       vertex 0.23097 -0.095671 -0.05
+    endloop
+  endfacet
+  facet normal 0.980784 -0.195096 0
+    outer loop
+       vertex 0.241481 -0.064705 -0.05
+       vertex 0.247861 -0.0326315 -0.05
+       vertex 0.247861 -0.0326315 0.05
+    endloop
+  endfacet
+  facet normal 0.980784 -0.195096 0
+    outer loop
+       vertex 0.247861 -0.0326315 0.05
+       vertex 0.241481 -0.064705 0.05
+       vertex 0.241481 -0.064705 -0.05
+    endloop
+  endfacet
+  facet normal 0.997858 -0.0654098 0
+    outer loop
+       vertex 0.25 -8.32665e-17 0.05
+       vertex 0.247861 -0.0326315 0.05
+       vertex 0.247861 -0.0326315 -0.05
+    endloop
+  endfacet
+  facet normal 0.997858 -0.0654098 -1.11896e-16
+    outer loop
+       vertex 0.247861 -0.0326315 -0.05
+       vertex 0.25 8.32665e-17 -0.05
+       vertex 0.25 -8.32665e-17 0.05
+    endloop
+  endfacet
+  facet normal 0.997858 0.0654098 0
+    outer loop
+       vertex 0.25 8.32665e-17 -0.05
+       vertex 0.247861 0.0326315 -0.05
+       vertex 0.247861 0.0326315 0.05
+    endloop
+  endfacet
+  facet normal 0.997858 0.0654098 1.11896e-16
+    outer loop
+       vertex 0.247861 0.0326315 0.05
+       vertex 0.25 -8.32665e-17 0.05
+       vertex 0.25 8.32665e-17 -0.05
+    endloop
+  endfacet
+  facet normal 0.980784 0.195096 0
+    outer loop
+       vertex 0.247861 0.0326315 -0.05
+       vertex 0.241481 0.064705 -0.05
+       vertex 0.241481 0.064705 0.05
+    endloop
+  endfacet
+  facet normal 0.980784 0.195096 0
+    outer loop
+       vertex 0.241481 0.064705 0.05
+       vertex 0.247861 0.0326315 0.05
+       vertex 0.247861 0.0326315 -0.05
+    endloop
+  endfacet
+  facet normal 0.946935 0.321425 0
+    outer loop
+       vertex 0.241481 0.064705 -0.05
+       vertex 0.23097 0.095671 -0.05
+       vertex 0.23097 0.095671 0.05
+    endloop
+  endfacet
+  facet normal 0.946935 0.321425 0
+    outer loop
+       vertex 0.23097 0.095671 0.05
+       vertex 0.241481 0.064705 0.05
+       vertex 0.241481 0.064705 -0.05
+    endloop
+  endfacet
+  facet normal 0.896866 0.442302 0
+    outer loop
+       vertex 0.23097 0.095671 -0.05
+       vertex 0.216506 0.125 -0.05
+       vertex 0.216506 0.125 0.05
+    endloop
+  endfacet
+  facet normal 0.896866 0.442302 0
+    outer loop
+       vertex 0.216506 0.125 0.05
+       vertex 0.23097 0.095671 0.05
+       vertex 0.23097 0.095671 -0.05
+    endloop
+  endfacet
+  facet normal 0.831476 0.555561 0
+    outer loop
+       vertex 0.216506 0.125 -0.05
+       vertex 0.198338 0.152191 -0.05
+       vertex 0.198338 0.152191 0.05
+    endloop
+  endfacet
+  facet normal 0.831476 0.555561 0
+    outer loop
+       vertex 0.198338 0.152191 0.05
+       vertex 0.216506 0.125 0.05
+       vertex 0.216506 0.125 -0.05
+    endloop
+  endfacet
+  facet normal 0.751817 0.659372 0
+    outer loop
+       vertex 0.198338 0.152191 -0.05
+       vertex 0.176776 0.176776 -0.05
+       vertex 0.176776 0.176776 0.05
+    endloop
+  endfacet
+  facet normal 0.751817 0.659372 0
+    outer loop
+       vertex 0.176776 0.176776 0.05
+       vertex 0.198338 0.152191 0.05
+       vertex 0.198338 0.152191 -0.05
+    endloop
+  endfacet
+  facet normal 0.659372 0.751817 0
+    outer loop
+       vertex 0.176776 0.176776 -0.05
+       vertex 0.152191 0.198338 -0.05
+       vertex 0.152191 0.198338 0.05
+    endloop
+  endfacet
+  facet normal 0.659372 0.751817 0
+    outer loop
+       vertex 0.152191 0.198338 0.05
+       vertex 0.176776 0.176776 0.05
+       vertex 0.176776 0.176776 -0.05
+    endloop
+  endfacet
+  facet normal 0.555561 0.831476 0
+    outer loop
+       vertex 0.152191 0.198338 -0.05
+       vertex 0.125 0.216506 -0.05
+       vertex 0.125 0.216506 0.05
+    endloop
+  endfacet
+  facet normal 0.555561 0.831476 0
+    outer loop
+       vertex 0.125 0.216506 0.05
+       vertex 0.152191 0.198338 0.05
+       vertex 0.152191 0.198338 -0.05
+    endloop
+  endfacet
+  facet normal 0.442302 0.896866 0
+    outer loop
+       vertex 0.125 0.216506 -0.05
+       vertex 0.095671 0.23097 -0.05
+       vertex 0.095671 0.23097 0.05
+    endloop
+  endfacet
+  facet normal 0.442302 0.896866 0
+    outer loop
+       vertex 0.095671 0.23097 0.05
+       vertex 0.125 0.216506 0.05
+       vertex 0.125 0.216506 -0.05
+    endloop
+  endfacet
+  facet normal 0.321425 0.946935 0
+    outer loop
+       vertex 0.095671 0.23097 -0.05
+       vertex 0.064705 0.241481 -0.05
+       vertex 0.064705 0.241481 0.05
+    endloop
+  endfacet
+  facet normal 0.321425 0.946935 0
+    outer loop
+       vertex 0.064705 0.241481 0.05
+       vertex 0.095671 0.23097 0.05
+       vertex 0.095671 0.23097 -0.05
+    endloop
+  endfacet
+  facet normal 0.195096 0.980784 0
+    outer loop
+       vertex 0.064705 0.241481 -0.05
+       vertex 0.0326315 0.247861 -0.05
+       vertex 0.0326315 0.247861 0.05
+    endloop
+  endfacet
+  facet normal 0.195096 0.980784 0
+    outer loop
+       vertex 0.0326315 0.247861 0.05
+       vertex 0.064705 0.241481 0.05
+       vertex 0.064705 0.241481 -0.05
+    endloop
+  endfacet
+  facet normal 0.0654098 0.997858 0
+    outer loop
+       vertex 8.32665e-17 0.25 0.05
+       vertex 0.0326315 0.247861 0.05
+       vertex 0.0326315 0.247861 -0.05
+    endloop
+  endfacet
+  facet normal 0.0654098 0.997858 -1.11896e-16
+    outer loop
+       vertex 0.0326315 0.247861 -0.05
+       vertex -8.32665e-17 0.25 -0.05
+       vertex 8.32665e-17 0.25 0.05
+    endloop
+  endfacet
+  facet normal -0.0654098 0.997858 0
+    outer loop
+       vertex -8.32665e-17 0.25 -0.05
+       vertex -0.0326315 0.247861 -0.05
+       vertex -0.0326315 0.247861 0.05
+    endloop
+  endfacet
+  facet normal -0.0654098 0.997858 1.11896e-16
+    outer loop
+       vertex -0.0326315 0.247861 0.05
+       vertex 8.32665e-17 0.25 0.05
+       vertex -8.32665e-17 0.25 -0.05
+    endloop
+  endfacet
+  facet normal -0.195096 0.980784 0
+    outer loop
+       vertex -0.0326315 0.247861 -0.05
+       vertex -0.064705 0.241481 -0.05
+       vertex -0.064705 0.241481 0.05
+    endloop
+  endfacet
+  facet normal -0.195096 0.980784 0
+    outer loop
+       vertex -0.064705 0.241481 0.05
+       vertex -0.0326315 0.247861 0.05
+       vertex -0.0326315 0.247861 -0.05
+    endloop
+  endfacet
+  facet normal -0.321425 0.946935 0
+    outer loop
+       vertex -0.064705 0.241481 -0.05
+       vertex -0.095671 0.23097 -0.05
+       vertex -0.095671 0.23097 0.05
+    endloop
+  endfacet
+  facet normal -0.321425 0.946935 0
+    outer loop
+       vertex -0.095671 0.23097 0.05
+       vertex -0.064705 0.241481 0.05
+       vertex -0.064705 0.241481 -0.05
+    endloop
+  endfacet
+  facet normal -0.442302 0.896866 0
+    outer loop
+       vertex -0.095671 0.23097 -0.05
+       vertex -0.125 0.216506 -0.05
+       vertex -0.125 0.216506 0.05
+    endloop
+  endfacet
+  facet normal -0.442302 0.896866 0
+    outer loop
+       vertex -0.125 0.216506 0.05
+       vertex -0.095671 0.23097 0.05
+       vertex -0.095671 0.23097 -0.05
+    endloop
+  endfacet
+  facet normal -0.555561 0.831476 0
+    outer loop
+       vertex -0.125 0.216506 -0.05
+       vertex -0.152191 0.198338 -0.05
+       vertex -0.152191 0.198338 0.05
+    endloop
+  endfacet
+  facet normal -0.555561 0.831476 0
+    outer loop
+       vertex -0.152191 0.198338 0.05
+       vertex -0.125 0.216506 0.05
+       vertex -0.125 0.216506 -0.05
+    endloop
+  endfacet
+  facet normal -0.659372 0.751817 0
+    outer loop
+       vertex -0.152191 0.198338 -0.05
+       vertex -0.176776 0.176776 -0.05
+       vertex -0.176776 0.176776 0.05
+    endloop
+  endfacet
+  facet normal -0.659372 0.751817 0
+    outer loop
+       vertex -0.176776 0.176776 0.05
+       vertex -0.152191 0.198338 0.05
+       vertex -0.152191 0.198338 -0.05
+    endloop
+  endfacet
+  facet normal -0.751817 0.659372 0
+    outer loop
+       vertex -0.176776 0.176776 -0.05
+       vertex -0.198338 0.152191 -0.05
+       vertex -0.198338 0.152191 0.05
+    endloop
+  endfacet
+  facet normal -0.751817 0.659372 0
+    outer loop
+       vertex -0.198338 0.152191 0.05
+       vertex -0.176776 0.176776 0.05
+       vertex -0.176776 0.176776 -0.05
+    endloop
+  endfacet
+  facet normal -0.831476 0.555561 0
+    outer loop
+       vertex -0.198338 0.152191 -0.05
+       vertex -0.216506 0.125 -0.05
+       vertex -0.216506 0.125 0.05
+    endloop
+  endfacet
+  facet normal -0.831476 0.555561 0
+    outer loop
+       vertex -0.216506 0.125 0.05
+       vertex -0.198338 0.152191 0.05
+       vertex -0.198338 0.152191 -0.05
+    endloop
+  endfacet
+  facet normal -0.896866 0.442302 0
+    outer loop
+       vertex -0.216506 0.125 -0.05
+       vertex -0.23097 0.095671 -0.05
+       vertex -0.23097 0.095671 0.05
+    endloop
+  endfacet
+  facet normal -0.896866 0.442302 0
+    outer loop
+       vertex -0.23097 0.095671 0.05
+       vertex -0.216506 0.125 0.05
+       vertex -0.216506 0.125 -0.05
+    endloop
+  endfacet
+  facet normal -0.946935 0.321425 0
+    outer loop
+       vertex -0.23097 0.095671 -0.05
+       vertex -0.241481 0.064705 -0.05
+       vertex -0.241481 0.064705 0.05
+    endloop
+  endfacet
+  facet normal -0.946935 0.321425 0
+    outer loop
+       vertex -0.241481 0.064705 0.05
+       vertex -0.23097 0.095671 0.05
+       vertex -0.23097 0.095671 -0.05
+    endloop
+  endfacet
+  facet normal -0.980784 0.195096 0
+    outer loop
+       vertex -0.241481 0.064705 -0.05
+       vertex -0.247861 0.0326315 -0.05
+       vertex -0.247861 0.0326315 0.05
+    endloop
+  endfacet
+  facet normal -0.980784 0.195096 0
+    outer loop
+       vertex -0.247861 0.0326315 0.05
+       vertex -0.241481 0.064705 0.05
+       vertex -0.241481 0.064705 -0.05
+    endloop
+  endfacet
+  facet normal -0.997858 0.0654098 0
+    outer loop
+       vertex -0.25 8.32665e-17 0.05
+       vertex -0.247861 0.0326315 0.05
+       vertex -0.247861 0.0326315 -0.05
+    endloop
+  endfacet
+  facet normal -0.997858 0.0654098 -1.11896e-16
+    outer loop
+       vertex -0.247861 0.0326315 -0.05
+       vertex -0.25 -8.32665e-17 -0.05
+       vertex -0.25 8.32665e-17 0.05
+    endloop
+  endfacet
+  facet normal -0.997858 -0.0654098 0
+    outer loop
+       vertex -0.25 -8.32665e-17 -0.05
+       vertex -0.247861 -0.0326315 -0.05
+       vertex -0.247861 -0.0326315 0.05
+    endloop
+  endfacet
+  facet normal -0.997858 -0.0654098 1.11896e-16
+    outer loop
+       vertex -0.247861 -0.0326315 0.05
+       vertex -0.25 8.32665e-17 0.05
+       vertex -0.25 -8.32665e-17 -0.05
+    endloop
+  endfacet
+  facet normal -0.980784 -0.195096 0
+    outer loop
+       vertex -0.247861 -0.0326315 -0.05
+       vertex -0.241481 -0.064705 -0.05
+       vertex -0.241481 -0.064705 0.05
+    endloop
+  endfacet
+  facet normal -0.980784 -0.195096 0
+    outer loop
+       vertex -0.241481 -0.064705 0.05
+       vertex -0.247861 -0.0326315 0.05
+       vertex -0.247861 -0.0326315 -0.05
+    endloop
+  endfacet
+  facet normal -0.946935 -0.321425 0
+    outer loop
+       vertex -0.241481 -0.064705 -0.05
+       vertex -0.23097 -0.095671 -0.05
+       vertex -0.23097 -0.095671 0.05
+    endloop
+  endfacet
+  facet normal -0.946935 -0.321425 0
+    outer loop
+       vertex -0.23097 -0.095671 0.05
+       vertex -0.241481 -0.064705 0.05
+       vertex -0.241481 -0.064705 -0.05
+    endloop
+  endfacet
+  facet normal -0.896866 -0.442302 0
+    outer loop
+       vertex -0.23097 -0.095671 -0.05
+       vertex -0.216506 -0.125 -0.05
+       vertex -0.216506 -0.125 0.05
+    endloop
+  endfacet
+  facet normal -0.896866 -0.442302 0
+    outer loop
+       vertex -0.216506 -0.125 0.05
+       vertex -0.23097 -0.095671 0.05
+       vertex -0.23097 -0.095671 -0.05
+    endloop
+  endfacet
+  facet normal -0.831476 -0.555561 0
+    outer loop
+       vertex -0.216506 -0.125 -0.05
+       vertex -0.198338 -0.152191 -0.05
+       vertex -0.198338 -0.152191 0.05
+    endloop
+  endfacet
+  facet normal -0.831476 -0.555561 0
+    outer loop
+       vertex -0.198338 -0.152191 0.05
+       vertex -0.216506 -0.125 0.05
+       vertex -0.216506 -0.125 -0.05
+    endloop
+  endfacet
+  facet normal -0.751817 -0.659372 0
+    outer loop
+       vertex -0.198338 -0.152191 -0.05
+       vertex -0.176776 -0.176776 -0.05
+       vertex -0.176776 -0.176776 0.05
+    endloop
+  endfacet
+  facet normal -0.751817 -0.659372 0
+    outer loop
+       vertex -0.176776 -0.176776 0.05
+       vertex -0.198338 -0.152191 0.05
+       vertex -0.198338 -0.152191 -0.05
+    endloop
+  endfacet
+  facet normal 0.659372 -0.751817 0
+    outer loop
+       vertex 0.176776 -0.176776 -0.05
+       vertex 0.176776 -0.176776 0.05
+       vertex 0.152191 -0.198338 0.05
+    endloop
+  endfacet
+  facet normal 0.659372 -0.751817 0
+    outer loop
+       vertex 0.152191 -0.198338 0.05
+       vertex 0.152191 -0.198338 -0.05
+       vertex 0.176776 -0.176776 -0.05
+    endloop
+  endfacet
+  facet normal 0.555561 -0.831476 0
+    outer loop
+       vertex 0.152191 -0.198338 -0.05
+       vertex 0.152191 -0.198338 0.05
+       vertex 0.125 -0.216506 0.05
+    endloop
+  endfacet
+  facet normal 0.555561 -0.831476 0
+    outer loop
+       vertex 0.125 -0.216506 0.05
+       vertex 0.125 -0.216506 -0.05
+       vertex 0.152191 -0.198338 -0.05
+    endloop
+  endfacet
+  facet normal 0.442302 -0.896866 0
+    outer loop
+       vertex 0.125 -0.216506 -0.05
+       vertex 0.125 -0.216506 0.05
+       vertex 0.095671 -0.23097 0.05
+    endloop
+  endfacet
+  facet normal 0.442302 -0.896866 0
+    outer loop
+       vertex 0.095671 -0.23097 0.05
+       vertex 0.095671 -0.23097 -0.05
+       vertex 0.125 -0.216506 -0.05
+    endloop
+  endfacet
+  facet normal 0.321425 -0.946935 0
+    outer loop
+       vertex 0.095671 -0.23097 -0.05
+       vertex 0.095671 -0.23097 0.05
+       vertex 0.064705 -0.241481 0.05
+    endloop
+  endfacet
+  facet normal 0.321425 -0.946935 0
+    outer loop
+       vertex 0.064705 -0.241481 0.05
+       vertex 0.064705 -0.241481 -0.05
+       vertex 0.095671 -0.23097 -0.05
+    endloop
+  endfacet
+  facet normal 0.195096 -0.980784 0
+    outer loop
+       vertex 0.064705 -0.241481 -0.05
+       vertex 0.064705 -0.241481 0.05
+       vertex 0.0326315 -0.247861 0.05
+    endloop
+  endfacet
+  facet normal 0.195096 -0.980784 0
+    outer loop
+       vertex 0.0326315 -0.247861 0.05
+       vertex 0.0326315 -0.247861 -0.05
+       vertex 0.064705 -0.241481 -0.05
+    endloop
+  endfacet
+  facet normal 0.0654098 -0.997858 0
+    outer loop
+       vertex 8.32665e-17 -0.25 -0.05
+       vertex 0.0326315 -0.247861 -0.05
+       vertex 0.0326315 -0.247861 0.05
+    endloop
+  endfacet
+  facet normal 0.0654098 -0.997858 1.11896e-16
+    outer loop
+       vertex 0.0326315 -0.247861 0.05
+       vertex -8.32665e-17 -0.25 0.05
+       vertex 8.32665e-17 -0.25 -0.05
+    endloop
+  endfacet
+  facet normal -0.0654098 -0.997858 0
+    outer loop
+       vertex -8.32665e-17 -0.25 0.05
+       vertex -0.0326315 -0.247861 0.05
+       vertex -0.0326315 -0.247861 -0.05
+    endloop
+  endfacet
+  facet normal -0.0654098 -0.997858 -1.11896e-16
+    outer loop
+       vertex -0.0326315 -0.247861 -0.05
+       vertex 8.32665e-17 -0.25 -0.05
+       vertex -8.32665e-17 -0.25 0.05
+    endloop
+  endfacet
+  facet normal -0.195096 -0.980784 0
+    outer loop
+       vertex -0.0326315 -0.247861 -0.05
+       vertex -0.0326315 -0.247861 0.05
+       vertex -0.064705 -0.241481 0.05
+    endloop
+  endfacet
+  facet normal -0.195096 -0.980784 -0
+    outer loop
+       vertex -0.064705 -0.241481 0.05
+       vertex -0.064705 -0.241481 -0.05
+       vertex -0.0326315 -0.247861 -0.05
+    endloop
+  endfacet
+  facet normal -0.321425 -0.946935 0
+    outer loop
+       vertex -0.064705 -0.241481 -0.05
+       vertex -0.064705 -0.241481 0.05
+       vertex -0.095671 -0.23097 0.05
+    endloop
+  endfacet
+  facet normal -0.321425 -0.946935 -0
+    outer loop
+       vertex -0.095671 -0.23097 0.05
+       vertex -0.095671 -0.23097 -0.05
+       vertex -0.064705 -0.241481 -0.05
+    endloop
+  endfacet
+  facet normal -0.442302 -0.896866 0
+    outer loop
+       vertex -0.095671 -0.23097 -0.05
+       vertex -0.095671 -0.23097 0.05
+       vertex -0.125 -0.216506 0.05
+    endloop
+  endfacet
+  facet normal -0.442302 -0.896866 -0
+    outer loop
+       vertex -0.125 -0.216506 0.05
+       vertex -0.125 -0.216506 -0.05
+       vertex -0.095671 -0.23097 -0.05
+    endloop
+  endfacet
+  facet normal -0.555561 -0.831476 0
+    outer loop
+       vertex -0.125 -0.216506 -0.05
+       vertex -0.125 -0.216506 0.05
+       vertex -0.152191 -0.198338 0.05
+    endloop
+  endfacet
+  facet normal -0.555561 -0.831476 -0
+    outer loop
+       vertex -0.152191 -0.198338 0.05
+       vertex -0.152191 -0.198338 -0.05
+       vertex -0.125 -0.216506 -0.05
+    endloop
+  endfacet
+  facet normal -0.659372 -0.751817 0
+    outer loop
+       vertex -0.152191 -0.198338 -0.05
+       vertex -0.152191 -0.198338 0.05
+       vertex -0.176776 -0.176776 0.05
+    endloop
+  endfacet
+  facet normal -0.659372 -0.751817 -0
+    outer loop
+       vertex -0.176776 -0.176776 0.05
+       vertex -0.176776 -0.176776 -0.05
+       vertex -0.152191 -0.198338 -0.05
+    endloop
+  endfacet
+endsolid cylinder

tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder2.ftr

@@ -0,0 +1,205 @@
+1(
+cylinder
+patch)
+
+96
+(
+(0.45 8.32665e-17 0.05)
+(0.45 -8.32665e-17 -0.05)
+(0.452139 0.0326315 -0.05)
+(0.452139 0.0326315 0.05)
+(0.452139 -0.0326315 -0.05)
+(0.452139 -0.0326315 0.05)
+(0.458519 0.064705 -0.05)
+(0.458519 0.064705 0.05)
+(0.458519 -0.064705 -0.05)
+(0.458519 -0.064705 0.05)
+(0.46903 0.095671 -0.05)
+(0.46903 0.095671 0.05)
+(0.46903 -0.095671 -0.05)
+(0.46903 -0.095671 0.05)
+(0.483494 0.125 -0.05)
+(0.483494 0.125 0.05)
+(0.483494 -0.125 -0.05)
+(0.483494 -0.125 0.05)
+(0.501662 0.152191 -0.05)
+(0.501662 0.152191 0.05)
+(0.501662 -0.152191 -0.05)
+(0.501662 -0.152191 0.05)
+(0.523224 0.176776 -0.05)
+(0.523224 0.176776 0.05)
+(0.523224 -0.176776 -0.05)
+(0.523224 -0.176776 0.05)
+(0.547809 0.198338 -0.05)
+(0.547809 0.198338 0.05)
+(0.547809 -0.198338 0.05)
+(0.547809 -0.198338 -0.05)
+(0.575 0.216506 -0.05)
+(0.575 0.216506 0.05)
+(0.575 -0.216506 0.05)
+(0.575 -0.216506 -0.05)
+(0.604329 0.23097 -0.05)
+(0.604329 0.23097 0.05)
+(0.604329 -0.23097 0.05)
+(0.604329 -0.23097 -0.05)
+(0.635295 0.241481 -0.05)
+(0.635295 0.241481 0.05)
+(0.635295 -0.241481 0.05)
+(0.635295 -0.241481 -0.05)
+(0.667368 0.247861 -0.05)
+(0.667368 0.247861 0.05)
+(0.667368 -0.247861 0.05)
+(0.667368 -0.247861 -0.05)
+(0.7 0.25 0.05)
+(0.7 0.25 -0.05)
+(0.7 -0.25 -0.05)
+(0.7 -0.25 0.05)
+(0.732632 0.247861 -0.05)
+(0.732632 0.247861 0.05)
+(0.732632 -0.247861 0.05)
+(0.732632 -0.247861 -0.05)
+(0.764705 0.241481 -0.05)
+(0.764705 0.241481 0.05)
+(0.764705 -0.241481 0.05)
+(0.764705 -0.241481 -0.05)
+(0.795671 0.23097 -0.05)
+(0.795671 0.23097 0.05)
+(0.795671 -0.23097 0.05)
+(0.795671 -0.23097 -0.05)
+(0.825 0.216506 -0.05)
+(0.825 0.216506 0.05)
+(0.825 -0.216506 0.05)
+(0.825 -0.216506 -0.05)
+(0.852191 0.198338 -0.05)
+(0.852191 0.198338 0.05)
+(0.852191 -0.198338 0.05)
+(0.852191 -0.198338 -0.05)
+(0.876776 -0.176776 -0.05)
+(0.876776 -0.176776 0.05)
+(0.876776 0.176776 -0.05)
+(0.876776 0.176776 0.05)
+(0.898338 -0.152191 -0.05)
+(0.898338 -0.152191 0.05)
+(0.898338 0.152191 -0.05)
+(0.898338 0.152191 0.05)
+(0.916506 -0.125 -0.05)
+(0.916506 -0.125 0.05)
+(0.916506 0.125 -0.05)
+(0.916506 0.125 0.05)
+(0.93097 -0.095671 -0.05)
+(0.93097 -0.095671 0.05)
+(0.93097 0.095671 -0.05)
+(0.93097 0.095671 0.05)
+(0.941481 -0.064705 -0.05)
+(0.941481 -0.064705 0.05)
+(0.941481 0.064705 -0.05)
+(0.941481 0.064705 0.05)
+(0.947861 -0.0326315 -0.05)
+(0.947861 -0.0326315 0.05)
+(0.947861 0.0326315 -0.05)
+(0.947861 0.0326315 0.05)
+(0.95 -8.32665e-17 0.05)
+(0.95 8.32665e-17 -0.05)
+)
+
+
+96
+(
+((70 74 75) 0)
+((75 71 70) 0)
+((74 78 79) 0)
+((79 75 74) 0)
+((78 82 83) 0)
+((83 79 78) 0)
+((82 86 87) 0)
+((87 83 82) 0)
+((86 90 91) 0)
+((91 87 86) 0)
+((94 91 90) 0)
+((90 95 94) 0)
+((95 92 93) 0)
+((93 94 95) 0)
+((92 88 89) 0)
+((89 93 92) 0)
+((88 84 85) 0)
+((85 89 88) 0)
+((84 80 81) 0)
+((81 85 84) 0)
+((80 76 77) 0)
+((77 81 80) 0)
+((76 72 73) 0)
+((73 77 76) 0)
+((72 66 67) 0)
+((67 73 72) 0)
+((66 62 63) 0)
+((63 67 66) 0)
+((62 58 59) 0)
+((59 63 62) 0)
+((58 54 55) 0)
+((55 59 58) 0)
+((54 50 51) 0)
+((51 55 54) 0)
+((46 51 50) 0)
+((50 47 46) 0)
+((47 42 43) 0)
+((43 46 47) 0)
+((42 38 39) 0)
+((39 43 42) 0)
+((38 34 35) 0)
+((35 39 38) 0)
+((34 30 31) 0)
+((31 35 34) 0)
+((30 26 27) 0)
+((27 31 30) 0)
+((26 22 23) 0)
+((23 27 26) 0)
+((22 18 19) 0)
+((19 23 22) 0)
+((18 14 15) 0)
+((15 19 18) 0)
+((14 10 11) 0)
+((11 15 14) 0)
+((10 6 7) 0)
+((7 11 10) 0)
+((6 2 3) 0)
+((3 7 6) 0)
+((0 3 2) 0)
+((2 1 0) 0)
+((1 4 5) 0)
+((5 0 1) 0)
+((4 8 9) 0)
+((9 5 4) 0)
+((8 12 13) 0)
+((13 9 8) 0)
+((12 16 17) 0)
+((17 13 12) 0)
+((16 20 21) 0)
+((21 17 16) 0)
+((20 24 25) 0)
+((25 21 20) 0)
+((70 71 68) 0)
+((68 69 70) 0)
+((69 68 64) 0)
+((64 65 69) 0)
+((65 64 60) 0)
+((60 61 65) 0)
+((61 60 56) 0)
+((56 57 61) 0)
+((57 56 52) 0)
+((52 53 57) 0)
+((48 53 52) 0)
+((52 49 48) 0)
+((49 44 45) 0)
+((45 48 49) 0)
+((45 44 40) 0)
+((40 41 45) 0)
+((41 40 36) 0)
+((36 37 41) 0)
+((37 36 32) 0)
+((32 33 37) 0)
+((33 32 28) 0)
+((28 29 33) 0)
+((29 28 25) 0)
+((25 24 29) 0)
+)
+

tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder2.stl

@@ -0,0 +1,674 @@
+solid cylinder
+  facet normal 0.751817 -0.659372 0
+    outer loop
+       vertex 0.876776 -0.176776 -0.05
+       vertex 0.898338 -0.152191 -0.05
+       vertex 0.898338 -0.152191 0.05
+    endloop
+  endfacet
+  facet normal 0.751817 -0.659372 0
+    outer loop
+       vertex 0.898338 -0.152191 0.05
+       vertex 0.876776 -0.176776 0.05
+       vertex 0.876776 -0.176776 -0.05
+    endloop
+  endfacet
+  facet normal 0.831476 -0.555561 0
+    outer loop
+       vertex 0.898338 -0.152191 -0.05
+       vertex 0.916506 -0.125 -0.05
+       vertex 0.916506 -0.125 0.05
+    endloop
+  endfacet
+  facet normal 0.831476 -0.555561 0
+    outer loop
+       vertex 0.916506 -0.125 0.05
+       vertex 0.898338 -0.152191 0.05
+       vertex 0.898338 -0.152191 -0.05
+    endloop
+  endfacet
+  facet normal 0.896866 -0.442302 0
+    outer loop
+       vertex 0.916506 -0.125 -0.05
+       vertex 0.93097 -0.095671 -0.05
+       vertex 0.93097 -0.095671 0.05
+    endloop
+  endfacet
+  facet normal 0.896866 -0.442302 0
+    outer loop
+       vertex 0.93097 -0.095671 0.05
+       vertex 0.916506 -0.125 0.05
+       vertex 0.916506 -0.125 -0.05
+    endloop
+  endfacet
+  facet normal 0.946935 -0.321425 0
+    outer loop
+       vertex 0.93097 -0.095671 -0.05
+       vertex 0.941481 -0.064705 -0.05
+       vertex 0.941481 -0.064705 0.05
+    endloop
+  endfacet
+  facet normal 0.946935 -0.321425 0
+    outer loop
+       vertex 0.941481 -0.064705 0.05
+       vertex 0.93097 -0.095671 0.05
+       vertex 0.93097 -0.095671 -0.05
+    endloop
+  endfacet
+  facet normal 0.980784 -0.195096 0
+    outer loop
+       vertex 0.941481 -0.064705 -0.05
+       vertex 0.947861 -0.0326315 -0.05
+       vertex 0.947861 -0.0326315 0.05
+    endloop
+  endfacet
+  facet normal 0.980784 -0.195096 0
+    outer loop
+       vertex 0.947861 -0.0326315 0.05
+       vertex 0.941481 -0.064705 0.05
+       vertex 0.941481 -0.064705 -0.05
+    endloop
+  endfacet
+  facet normal 0.997858 -0.0654098 0
+    outer loop
+       vertex 0.95 -8.32665e-17 0.05
+       vertex 0.947861 -0.0326315 0.05
+       vertex 0.947861 -0.0326315 -0.05
+    endloop
+  endfacet
+  facet normal 0.997858 -0.0654098 -1.07752e-16
+    outer loop
+       vertex 0.947861 -0.0326315 -0.05
+       vertex 0.95 8.32665e-17 -0.05
+       vertex 0.95 -8.32665e-17 0.05
+    endloop
+  endfacet
+  facet normal 0.997858 0.0654098 0
+    outer loop
+       vertex 0.95 8.32665e-17 -0.05
+       vertex 0.947861 0.0326315 -0.05
+       vertex 0.947861 0.0326315 0.05
+    endloop
+  endfacet
+  facet normal 0.997858 0.0654098 1.07752e-16
+    outer loop
+       vertex 0.947861 0.0326315 0.05
+       vertex 0.95 -8.32665e-17 0.05
+       vertex 0.95 8.32665e-17 -0.05
+    endloop
+  endfacet
+  facet normal 0.980784 0.195096 0
+    outer loop
+       vertex 0.947861 0.0326315 -0.05
+       vertex 0.941481 0.064705 -0.05
+       vertex 0.941481 0.064705 0.05
+    endloop
+  endfacet
+  facet normal 0.980784 0.195096 0
+    outer loop
+       vertex 0.941481 0.064705 0.05
+       vertex 0.947861 0.0326315 0.05
+       vertex 0.947861 0.0326315 -0.05
+    endloop
+  endfacet
+  facet normal 0.946935 0.321425 0
+    outer loop
+       vertex 0.941481 0.064705 -0.05
+       vertex 0.93097 0.095671 -0.05
+       vertex 0.93097 0.095671 0.05
+    endloop
+  endfacet
+  facet normal 0.946935 0.321425 0
+    outer loop
+       vertex 0.93097 0.095671 0.05
+       vertex 0.941481 0.064705 0.05
+       vertex 0.941481 0.064705 -0.05
+    endloop
+  endfacet
+  facet normal 0.896866 0.442302 0
+    outer loop
+       vertex 0.93097 0.095671 -0.05
+       vertex 0.916506 0.125 -0.05
+       vertex 0.916506 0.125 0.05
+    endloop
+  endfacet
+  facet normal 0.896866 0.442302 0
+    outer loop
+       vertex 0.916506 0.125 0.05
+       vertex 0.93097 0.095671 0.05
+       vertex 0.93097 0.095671 -0.05
+    endloop
+  endfacet
+  facet normal 0.831476 0.555561 0
+    outer loop
+       vertex 0.916506 0.125 -0.05
+       vertex 0.898338 0.152191 -0.05
+       vertex 0.898338 0.152191 0.05
+    endloop
+  endfacet
+  facet normal 0.831476 0.555561 0
+    outer loop
+       vertex 0.898338 0.152191 0.05
+       vertex 0.916506 0.125 0.05
+       vertex 0.916506 0.125 -0.05
+    endloop
+  endfacet
+  facet normal 0.751817 0.659372 0
+    outer loop
+       vertex 0.898338 0.152191 -0.05
+       vertex 0.876776 0.176776 -0.05
+       vertex 0.876776 0.176776 0.05
+    endloop
+  endfacet
+  facet normal 0.751817 0.659372 0
+    outer loop
+       vertex 0.876776 0.176776 0.05
+       vertex 0.898338 0.152191 0.05
+       vertex 0.898338 0.152191 -0.05
+    endloop
+  endfacet
+  facet normal 0.659372 0.751817 0
+    outer loop
+       vertex 0.876776 0.176776 -0.05
+       vertex 0.852191 0.198338 -0.05
+       vertex 0.852191 0.198338 0.05
+    endloop
+  endfacet
+  facet normal 0.659372 0.751817 0
+    outer loop
+       vertex 0.852191 0.198338 0.05
+       vertex 0.876776 0.176776 0.05
+       vertex 0.876776 0.176776 -0.05
+    endloop
+  endfacet
+  facet normal 0.555561 0.831476 0
+    outer loop
+       vertex 0.852191 0.198338 -0.05
+       vertex 0.825 0.216506 -0.05
+       vertex 0.825 0.216506 0.05
+    endloop
+  endfacet
+  facet normal 0.555561 0.831476 0
+    outer loop
+       vertex 0.825 0.216506 0.05
+       vertex 0.852191 0.198338 0.05
+       vertex 0.852191 0.198338 -0.05
+    endloop
+  endfacet
+  facet normal 0.442302 0.896866 0
+    outer loop
+       vertex 0.825 0.216506 -0.05
+       vertex 0.795671 0.23097 -0.05
+       vertex 0.795671 0.23097 0.05
+    endloop
+  endfacet
+  facet normal 0.442302 0.896866 0
+    outer loop
+       vertex 0.795671 0.23097 0.05
+       vertex 0.825 0.216506 0.05
+       vertex 0.825 0.216506 -0.05
+    endloop
+  endfacet
+  facet normal 0.321425 0.946935 0
+    outer loop
+       vertex 0.795671 0.23097 -0.05
+       vertex 0.764705 0.241481 -0.05
+       vertex 0.764705 0.241481 0.05
+    endloop
+  endfacet
+  facet normal 0.321425 0.946935 0
+    outer loop
+       vertex 0.764705 0.241481 0.05
+       vertex 0.795671 0.23097 0.05
+       vertex 0.795671 0.23097 -0.05
+    endloop
+  endfacet
+  facet normal 0.195096 0.980784 0
+    outer loop
+       vertex 0.764705 0.241481 -0.05
+       vertex 0.732632 0.247861 -0.05
+       vertex 0.732632 0.247861 0.05
+    endloop
+  endfacet
+  facet normal 0.195096 0.980784 0
+    outer loop
+       vertex 0.732632 0.247861 0.05
+       vertex 0.764705 0.241481 0.05
+       vertex 0.764705 0.241481 -0.05
+    endloop
+  endfacet
+  facet normal 0.0654098 0.997858 0
+    outer loop
+       vertex 0.7 0.25 0.05
+       vertex 0.732632 0.247861 0.05
+       vertex 0.732632 0.247861 -0.05
+    endloop
+  endfacet
+  facet normal 0.0654098 0.997858 -1.49195e-16
+    outer loop
+       vertex 0.732632 0.247861 -0.05
+       vertex 0.7 0.25 -0.05
+       vertex 0.7 0.25 0.05
+    endloop
+  endfacet
+  facet normal -0.0654098 0.997858 0
+    outer loop
+       vertex 0.7 0.25 -0.05
+       vertex 0.667368 0.247861 -0.05
+       vertex 0.667368 0.247861 0.05
+    endloop
+  endfacet
+  facet normal -0.0654098 0.997858 1.49195e-16
+    outer loop
+       vertex 0.667368 0.247861 0.05
+       vertex 0.7 0.25 0.05
+       vertex 0.7 0.25 -0.05
+    endloop
+  endfacet
+  facet normal -0.195096 0.980784 0
+    outer loop
+       vertex 0.667368 0.247861 -0.05
+       vertex 0.635295 0.241481 -0.05
+       vertex 0.635295 0.241481 0.05
+    endloop
+  endfacet
+  facet normal -0.195096 0.980784 0
+    outer loop
+       vertex 0.635295 0.241481 0.05
+       vertex 0.667368 0.247861 0.05
+       vertex 0.667368 0.247861 -0.05
+    endloop
+  endfacet
+  facet normal -0.321425 0.946935 0
+    outer loop
+       vertex 0.635295 0.241481 -0.05
+       vertex 0.604329 0.23097 -0.05
+       vertex 0.604329 0.23097 0.05
+    endloop
+  endfacet
+  facet normal -0.321425 0.946935 0
+    outer loop
+       vertex 0.604329 0.23097 0.05
+       vertex 0.635295 0.241481 0.05
+       vertex 0.635295 0.241481 -0.05
+    endloop
+  endfacet
+  facet normal -0.442302 0.896866 0
+    outer loop
+       vertex 0.604329 0.23097 -0.05
+       vertex 0.575 0.216506 -0.05
+       vertex 0.575 0.216506 0.05
+    endloop
+  endfacet
+  facet normal -0.442302 0.896866 0
+    outer loop
+       vertex 0.575 0.216506 0.05
+       vertex 0.604329 0.23097 0.05
+       vertex 0.604329 0.23097 -0.05
+    endloop
+  endfacet
+  facet normal -0.555561 0.831476 0
+    outer loop
+       vertex 0.575 0.216506 -0.05
+       vertex 0.547809 0.198338 -0.05
+       vertex 0.547809 0.198338 0.05
+    endloop
+  endfacet
+  facet normal -0.555561 0.831476 0
+    outer loop
+       vertex 0.547809 0.198338 0.05
+       vertex 0.575 0.216506 0.05
+       vertex 0.575 0.216506 -0.05
+    endloop
+  endfacet
+  facet normal -0.659372 0.751817 0
+    outer loop
+       vertex 0.547809 0.198338 -0.05
+       vertex 0.523224 0.176776 -0.05
+       vertex 0.523224 0.176776 0.05
+    endloop
+  endfacet
+  facet normal -0.659372 0.751817 0
+    outer loop
+       vertex 0.523224 0.176776 0.05
+       vertex 0.547809 0.198338 0.05
+       vertex 0.547809 0.198338 -0.05
+    endloop
+  endfacet
+  facet normal -0.751817 0.659372 0
+    outer loop
+       vertex 0.523224 0.176776 -0.05
+       vertex 0.501662 0.152191 -0.05
+       vertex 0.501662 0.152191 0.05
+    endloop
+  endfacet
+  facet normal -0.751817 0.659372 0
+    outer loop
+       vertex 0.501662 0.152191 0.05
+       vertex 0.523224 0.176776 0.05
+       vertex 0.523224 0.176776 -0.05
+    endloop
+  endfacet
+  facet normal -0.831476 0.555561 0
+    outer loop
+       vertex 0.501662 0.152191 -0.05
+       vertex 0.483494 0.125 -0.05
+       vertex 0.483494 0.125 0.05
+    endloop
+  endfacet
+  facet normal -0.831476 0.555561 0
+    outer loop
+       vertex 0.483494 0.125 0.05
+       vertex 0.501662 0.152191 0.05
+       vertex 0.501662 0.152191 -0.05
+    endloop
+  endfacet
+  facet normal -0.896866 0.442302 0
+    outer loop
+       vertex 0.483494 0.125 -0.05
+       vertex 0.46903 0.095671 -0.05
+       vertex 0.46903 0.095671 0.05
+    endloop
+  endfacet
+  facet normal -0.896866 0.442302 0
+    outer loop
+       vertex 0.46903 0.095671 0.05
+       vertex 0.483494 0.125 0.05
+       vertex 0.483494 0.125 -0.05
+    endloop
+  endfacet
+  facet normal -0.946935 0.321425 0
+    outer loop
+       vertex 0.46903 0.095671 -0.05
+       vertex 0.458519 0.064705 -0.05
+       vertex 0.458519 0.064705 0.05
+    endloop
+  endfacet
+  facet normal -0.946935 0.321425 0
+    outer loop
+       vertex 0.458519 0.064705 0.05
+       vertex 0.46903 0.095671 0.05
+       vertex 0.46903 0.095671 -0.05
+    endloop
+  endfacet
+  facet normal -0.980784 0.195096 0
+    outer loop
+       vertex 0.458519 0.064705 -0.05
+       vertex 0.452139 0.0326315 -0.05
+       vertex 0.452139 0.0326315 0.05
+    endloop
+  endfacet
+  facet normal -0.980784 0.195096 0
+    outer loop
+       vertex 0.452139 0.0326315 0.05
+       vertex 0.458519 0.064705 0.05
+       vertex 0.458519 0.064705 -0.05
+    endloop
+  endfacet
+  facet normal -0.997858 0.0654098 0
+    outer loop
+       vertex 0.45 8.32665e-17 0.05
+       vertex 0.452139 0.0326315 0.05
+       vertex 0.452139 0.0326315 -0.05
+    endloop
+  endfacet
+  facet normal -0.997858 0.0654098 -1.07752e-16
+    outer loop
+       vertex 0.452139 0.0326315 -0.05
+       vertex 0.45 -8.32665e-17 -0.05
+       vertex 0.45 8.32665e-17 0.05
+    endloop
+  endfacet
+  facet normal -0.997858 -0.0654098 0
+    outer loop
+       vertex 0.45 -8.32665e-17 -0.05
+       vertex 0.452139 -0.0326315 -0.05
+       vertex 0.452139 -0.0326315 0.05
+    endloop
+  endfacet
+  facet normal -0.997858 -0.0654098 1.07752e-16
+    outer loop
+       vertex 0.452139 -0.0326315 0.05
+       vertex 0.45 8.32665e-17 0.05
+       vertex 0.45 -8.32665e-17 -0.05
+    endloop
+  endfacet
+  facet normal -0.980784 -0.195096 0
+    outer loop
+       vertex 0.452139 -0.0326315 -0.05
+       vertex 0.458519 -0.064705 -0.05
+       vertex 0.458519 -0.064705 0.05
+    endloop
+  endfacet
+  facet normal -0.980784 -0.195096 0
+    outer loop
+       vertex 0.458519 -0.064705 0.05
+       vertex 0.452139 -0.0326315 0.05
+       vertex 0.452139 -0.0326315 -0.05
+    endloop
+  endfacet
+  facet normal -0.946935 -0.321425 0
+    outer loop
+       vertex 0.458519 -0.064705 -0.05
+       vertex 0.46903 -0.095671 -0.05
+       vertex 0.46903 -0.095671 0.05
+    endloop
+  endfacet
+  facet normal -0.946935 -0.321425 0
+    outer loop
+       vertex 0.46903 -0.095671 0.05
+       vertex 0.458519 -0.064705 0.05
+       vertex 0.458519 -0.064705 -0.05
+    endloop
+  endfacet
+  facet normal -0.896866 -0.442302 0
+    outer loop
+       vertex 0.46903 -0.095671 -0.05
+       vertex 0.483494 -0.125 -0.05
+       vertex 0.483494 -0.125 0.05
+    endloop
+  endfacet
+  facet normal -0.896866 -0.442302 0
+    outer loop
+       vertex 0.483494 -0.125 0.05
+       vertex 0.46903 -0.095671 0.05
+       vertex 0.46903 -0.095671 -0.05
+    endloop
+  endfacet
+  facet normal -0.831476 -0.555561 0
+    outer loop
+       vertex 0.483494 -0.125 -0.05
+       vertex 0.501662 -0.152191 -0.05
+       vertex 0.501662 -0.152191 0.05
+    endloop
+  endfacet
+  facet normal -0.831476 -0.555561 0
+    outer loop
+       vertex 0.501662 -0.152191 0.05
+       vertex 0.483494 -0.125 0.05
+       vertex 0.483494 -0.125 -0.05
+    endloop
+  endfacet
+  facet normal -0.751817 -0.659372 0
+    outer loop
+       vertex 0.501662 -0.152191 -0.05
+       vertex 0.523224 -0.176776 -0.05
+       vertex 0.523224 -0.176776 0.05
+    endloop
+  endfacet
+  facet normal -0.751817 -0.659372 0
+    outer loop
+       vertex 0.523224 -0.176776 0.05
+       vertex 0.501662 -0.152191 0.05
+       vertex 0.501662 -0.152191 -0.05
+    endloop
+  endfacet
+  facet normal 0.659372 -0.751817 0
+    outer loop
+       vertex 0.876776 -0.176776 -0.05
+       vertex 0.876776 -0.176776 0.05
+       vertex 0.852191 -0.198338 0.05
+    endloop
+  endfacet
+  facet normal 0.659372 -0.751817 0
+    outer loop
+       vertex 0.852191 -0.198338 0.05
+       vertex 0.852191 -0.198338 -0.05
+       vertex 0.876776 -0.176776 -0.05
+    endloop
+  endfacet
+  facet normal 0.555561 -0.831476 0
+    outer loop
+       vertex 0.852191 -0.198338 -0.05
+       vertex 0.852191 -0.198338 0.05
+       vertex 0.825 -0.216506 0.05
+    endloop
+  endfacet
+  facet normal 0.555561 -0.831476 0
+    outer loop
+       vertex 0.825 -0.216506 0.05
+       vertex 0.825 -0.216506 -0.05
+       vertex 0.852191 -0.198338 -0.05
+    endloop
+  endfacet
+  facet normal 0.442302 -0.896866 0
+    outer loop
+       vertex 0.825 -0.216506 -0.05
+       vertex 0.825 -0.216506 0.05
+       vertex 0.795671 -0.23097 0.05
+    endloop
+  endfacet
+  facet normal 0.442302 -0.896866 0
+    outer loop
+       vertex 0.795671 -0.23097 0.05
+       vertex 0.795671 -0.23097 -0.05
+       vertex 0.825 -0.216506 -0.05
+    endloop
+  endfacet
+  facet normal 0.321425 -0.946935 0
+    outer loop
+       vertex 0.795671 -0.23097 -0.05
+       vertex 0.795671 -0.23097 0.05
+       vertex 0.764705 -0.241481 0.05
+    endloop
+  endfacet
+  facet normal 0.321425 -0.946935 0
+    outer loop
+       vertex 0.764705 -0.241481 0.05
+       vertex 0.764705 -0.241481 -0.05
+       vertex 0.795671 -0.23097 -0.05
+    endloop
+  endfacet
+  facet normal 0.195096 -0.980784 0
+    outer loop
+       vertex 0.764705 -0.241481 -0.05
+       vertex 0.764705 -0.241481 0.05
+       vertex 0.732632 -0.247861 0.05
+    endloop
+  endfacet
+  facet normal 0.195096 -0.980784 0
+    outer loop
+       vertex 0.732632 -0.247861 0.05
+       vertex 0.732632 -0.247861 -0.05
+       vertex 0.764705 -0.241481 -0.05
+    endloop
+  endfacet
+  facet normal 0.0654098 -0.997858 0
+    outer loop
+       vertex 0.7 -0.25 -0.05
+       vertex 0.732632 -0.247861 -0.05
+       vertex 0.732632 -0.247861 0.05
+    endloop
+  endfacet
+  facet normal 0.0654098 -0.997858 1.49195e-16
+    outer loop
+       vertex 0.732632 -0.247861 0.05
+       vertex 0.7 -0.25 0.05
+       vertex 0.7 -0.25 -0.05
+    endloop
+  endfacet
+  facet normal -0.0654098 -0.997858 0
+    outer loop
+       vertex 0.7 -0.25 0.05
+       vertex 0.667368 -0.247861 0.05
+       vertex 0.667368 -0.247861 -0.05
+    endloop
+  endfacet
+  facet normal -0.0654098 -0.997858 -1.49195e-16
+    outer loop
+       vertex 0.667368 -0.247861 -0.05
+       vertex 0.7 -0.25 -0.05
+       vertex 0.7 -0.25 0.05
+    endloop
+  endfacet
+  facet normal -0.195096 -0.980784 0
+    outer loop
+       vertex 0.667368 -0.247861 -0.05
+       vertex 0.667368 -0.247861 0.05
+       vertex 0.635295 -0.241481 0.05
+    endloop
+  endfacet
+  facet normal -0.195096 -0.980784 -0
+    outer loop
+       vertex 0.635295 -0.241481 0.05
+       vertex 0.635295 -0.241481 -0.05
+       vertex 0.667368 -0.247861 -0.05
+    endloop
+  endfacet
+  facet normal -0.321425 -0.946935 0
+    outer loop
+       vertex 0.635295 -0.241481 -0.05
+       vertex 0.635295 -0.241481 0.05
+       vertex 0.604329 -0.23097 0.05
+    endloop
+  endfacet
+  facet normal -0.321425 -0.946935 -0
+    outer loop
+       vertex 0.604329 -0.23097 0.05
+       vertex 0.604329 -0.23097 -0.05
+       vertex 0.635295 -0.241481 -0.05
+    endloop
+  endfacet
+  facet normal -0.442302 -0.896866 0
+    outer loop
+       vertex 0.604329 -0.23097 -0.05
+       vertex 0.604329 -0.23097 0.05
+       vertex 0.575 -0.216506 0.05
+    endloop
+  endfacet
+  facet normal -0.442302 -0.896866 -0
+    outer loop
+       vertex 0.575 -0.216506 0.05
+       vertex 0.575 -0.216506 -0.05
+       vertex 0.604329 -0.23097 -0.05
+    endloop
+  endfacet
+  facet normal -0.555561 -0.831476 0
+    outer loop
+       vertex 0.575 -0.216506 -0.05
+       vertex 0.575 -0.216506 0.05
+       vertex 0.547809 -0.198338 0.05
+    endloop
+  endfacet
+  facet normal -0.555561 -0.831476 -0
+    outer loop
+       vertex 0.547809 -0.198338 0.05
+       vertex 0.547809 -0.198338 -0.05
+       vertex 0.575 -0.216506 -0.05
+    endloop
+  endfacet
+  facet normal -0.659372 -0.751817 0
+    outer loop
+       vertex 0.547809 -0.198338 -0.05
+       vertex 0.547809 -0.198338 0.05
+       vertex 0.523224 -0.176776 0.05
+    endloop
+  endfacet
+  facet normal -0.659372 -0.751817 -0
+    outer loop
+       vertex 0.523224 -0.176776 0.05
+       vertex 0.523224 -0.176776 -0.05
+       vertex 0.547809 -0.198338 -0.05
+    endloop
+  endfacet
+endsolid cylinder

tutorials/immersedBoundary/twoIbPatches/constant/turbulenceProperties

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      turbulenceProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType  RASModel;
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/save/blockMeshDict

@@ -0,0 +1,70 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      blockMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+convertToMeters 1;
+
+vertices
+(
+    (-1   -0.5 -0.04)
+    ( 2   -0.5 -0.04)
+    ( 2    0.5 -0.04)
+    (-1    0.5 -0.04)
+    (-1   -0.5  0.04)
+    ( 2   -0.5  0.04)
+    ( 2    0.5  0.04)
+    (-1    0.5  0.04)
+);
+
+blocks
+(
+    hex (0 1 2 3 4 5 6 7) (75 25 1) simpleGrading (1 1 1)
+);
+
+edges
+(
+);
+
+patches
+(
+    patch in
+    (
+        (0 4 7 3)
+    )
+    patch out
+    (
+        (2 6 5 1)
+    )
+    patch top
+    (
+        (3 7 6 2)
+    )
+    patch bottom
+    (
+        (1 5 4 0)
+    )
+    empty frontAndBack
+    (
+        (0 3 2 1)
+        (4 5 6 7)
+    )
+);
+
+mergePatchPairs
+(
+);
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/save/boundary

@@ -0,0 +1,71 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.1                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       polyBoundaryMesh;
+    location    "constant/polyMesh";
+    object      boundary;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+7
+(
+    ibCylinder
+    {
+        type            immersedBoundary;
+        nFaces          0;
+        startFace       3650;
+
+        internalFlow    no;
+        isWall          yes;
+    }
+    ibCylinder2
+    {
+        type            immersedBoundary;
+        nFaces          0;
+        startFace       3650;
+
+        internalFlow    no;
+        isWall          yes;
+    }
+    in
+    {
+        type            patch;
+        nFaces          25;
+        startFace       3650;
+    }
+    out
+    {
+        type            patch;
+        nFaces          25;
+        startFace       3675;
+    }
+    top
+    {
+        type            patch;
+        nFaces          75;
+        startFace       3700;
+    }
+    bottom
+    {
+        type            patch;
+        nFaces          75;
+        startFace       3775;
+    }
+    frontAndBack
+    {
+        type            empty;
+        nFaces          3750;
+        startFace       3850;
+    }
+)
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/save/k

@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      k;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 0.375;
+
+boundaryField
+{
+    ibCylinder
+    {
+        type            immersedBoundaryKqRWallFunction;
+        patchType       immersedBoundary;
+        value           uniform 0.375;
+    }
+    top
+    {
+        type            zeroGradient;
+    }
+    in
+    {
+        type            fixedValue;
+        value           uniform 0.375;
+    }
+    out
+    {
+        type            inletOutlet;
+        inletValue      uniform 0.375;
+        value           uniform 0.375;
+    }
+    bottom
+    {
+        type            zeroGradient;
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/save/nut

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      nut;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -1 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    ibCylinder
+    {
+        type            immersedBoundaryNutWallFunction;
+        patchType       immersedBoundary;
+        value           uniform 0;
+    }
+    top
+    {
+        type            zeroGradient;
+        value           uniform 0;
+    }
+    in
+    {
+        type            zeroGradient;
+        value           uniform 0;
+    }
+    out
+    {
+        type            zeroGradient;
+        value           uniform 0;
+    }
+    bottom
+    {
+        type            zeroGradient;
+        value           uniform 0;
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/save/omega

@@ -0,0 +1,60 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      omega;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 -1 0 0 0 0];
+
+internalField   uniform 500;
+
+boundaryField
+{
+    ibCylinder
+    {
+        type            immersedBoundaryOmegaWallFunction;
+        patchType       immersedBoundary;
+        Cmu             0.09;
+        kappa           0.41;
+        E               9.8;
+
+        value           uniform 500;
+    }
+    top
+    {
+        type            zeroGradient;
+    }
+    in
+    {
+        type            fixedValue;
+        value           uniform 500;
+    }
+    out
+    {
+        type            inletOutlet;
+        inletValue      uniform 500;
+        value           uniform 500;
+    }
+    bottom
+    {
+        type            zeroGradient;
+    }
+    frontAndBack
+    {
+        type            empty;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/system/controlDict

@@ -0,0 +1,152 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+DebugSwitches
+{
+
+  //     immersedBoundary 1;
+  //    immersedBoundaryFvPatch 10;
+    // polyMesh 2;
+    // fvMesh   2;
+  /*
+    immersedBoundary 2;
+    immersedBoundaryCorrectedMeshFields 2;
+    immersedBoundaryEpsilonWallFunction 2;
+    immersedBoundaryKqRWallFunction 2;
+    immersedBoundaryNutWallFunction 2;
+    immersedBoundaryOmegaWallFunction 2;
+    foo 2;
+  */
+}
+
+application     pimpleDyMIbFoam;
+
+startFrom       startTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         5;
+
+deltaT          0.01;
+
+// writeControl    runTime;
+// writeInterval   0.2;
+
+writeControl    timeStep;
+writeInterval   5;
+
+purgeWrite      0;
+
+writeFormat     binary;
+
+writePrecision  6;
+
+writeCompression uncompressed;
+
+timeFormat      general;
+
+timePrecision   6;
+
+runTimeModifiable yes;
+
+adjustTimeStep  no;
+
+maxCo           0.2;
+
+libs
+(
+    "libimmersedBoundary.so"
+    "libimmersedBoundaryDynamicMesh.so"
+    //   "libimmersedBoundaryTurbulence.so"
+"libincompressibleImmersedBoundaryTurbulence.so"
+);
+
+functions
+(
+    forces
+    {
+        type               forces;
+        functionObjectLibs ("libforces.so");
+
+        outputControl timeStep;
+        outputInterval 1;
+        patches (  ibCylinder2);
+
+        pName       p;
+        UName       U;
+        rhoName     rhoInf;
+        rhoInf      1;
+
+        log         true;
+        CofR        ( 0 0 0 );
+
+        Aref 0.05;
+        Uref 1;
+    }
+    /*
+    probes
+    {
+        type            probes;
+        functionObjectLibs ("libsampling.so");
+        outputControl   timeStep;
+        outputInterval  1;
+        probeLocations
+        (
+            ( 1.5 0 0)
+            ( 1.75 0 0)
+        );
+        fields
+        (
+            p
+        );
+    }
+    minMaxP
+    {
+        type minMaxField;
+        functionObjectLibs ("libfieldFunctionObjects.so");
+        name p;
+    }
+    minMaxU
+    {
+        type minMaxField;
+        functionObjectLibs ("libfieldFunctionObjects.so");
+        name U;
+    }
+    minMaxNut
+    {
+        type minMaxField;
+        functionObjectLibs ("libfieldFunctionObjects.so");
+        name nut;
+    }
+    minMaxK
+    {
+        type minMaxField;
+        functionObjectLibs ("libfieldFunctionObjects.so");
+        name k;
+    }
+    minMaxNut
+    {
+        type minMaxField;
+        functionObjectLibs ("libfieldFunctionObjects.so");
+        name epsilon;
+    }
+    */
+);
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/system/decomposeParDict

@@ -0,0 +1,34 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      decomposeParDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+numberOfSubdomains 3;
+
+method scotch;
+
+simpleCoeffs
+{
+    n (5 2 1);
+    delta           0.001;
+}
+
+distributed     no;
+
+roots
+(
+);
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/system/fvSchemes

@@ -0,0 +1,54 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      fvSchemes;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+ddtSchemes
+{
+    default Euler;
+}
+
+gradSchemes
+{
+    default         cellLimited Gauss linear 1;
+}
+
+divSchemes
+{
+    default         none;
+    div(phi,U)      Gauss linearUpwindV grad(U);//Gauss upwind;
+    div(phi,k)      Gauss upwind;
+    div(phi,epsilon) Gauss upwind;
+    div(phi,omega)  Gauss upwind;
+
+    div((nuEff*dev(T(grad(U))))) Gauss linear;
+}
+
+laplacianSchemes
+{
+    default         Gauss linear limited 0.5;
+}
+
+interpolationSchemes
+{
+    default         linear;
+}
+
+snGradSchemes
+{
+    default         limited 0.5;
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/system/fvSolution

@@ -0,0 +1,78 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      fvSolution;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solvers
+{
+    "p|pFinal|pcorr"
+    {
+        solver          CG;
+        preconditioner  Cholesky;
+        minIter         1;
+        maxIter         1000;
+        tolerance       1e-7;
+        relTol          0.0;
+    }
+    "U|UFinal|k|epsilon|omega"
+    {
+        solver           BiCGStab;
+        preconditioner   ILU0;
+
+        minIter          1;
+        maxIter          1000;
+        tolerance        1e-08;
+        relTol           0;
+    }
+}
+
+SIMPLE
+{
+    nNonOrthogonalCorrectors 1;
+
+    pRefPoint (0 -0.45 0.05);
+    pRefValue 0;
+}
+
+PIMPLE
+{
+    nOuterCorrectors 4;
+    nCorrectors     2;
+    nNonOrthogonalCorrectors 0;
+
+    pRefPoint (0 -0.45 0.05);
+    pRefValue 0;
+    checkMeshCourantNo no;
+    correctPhi yes;
+}
+
+PISO
+{
+    nNonOrthogonalCorrectors 1;
+}
+
+relaxationFactors
+{
+    equations
+    {
+        U         1;
+        UFinal    1;
+    }
+    fields
+    {
+    }
+}
+
+
+// ************************************************************************* //

tutorials/immersedBoundary/twoIbPatches/system/mapFieldsDict

@@ -0,0 +1,31 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | foam-extend: Open Source CFD                    |
+|  \\    /   O peration     | Version:     4.0                                |
+|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    object      mapFieldsDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+// List of pairs of source/target patches for mapping
+patchMap
+(
+//     lid movingWall
+);
+
+// List of target patches cutting the source domain (these need to be
+// handled specially e.g. interpolated from internal values)
+cuttingPatches
+(
+//     fixedWalls
+);
+
+
+// ************************************************************************* //