initial commit

.gitignore

@@ -0,0 +1,7 @@
+constant/polyMesh
+processor*
+log.*
+dynamicCode
+postProcessing
+temp
+anim

cell-simple/0/U

@@ -0,0 +1,34 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    inlet  { type cyclicAMI; };
+    outlet  { type cyclicAMI; };
+    walls { type noSlip; }
+    cperiodic-west { type cyclicAMI; }
+    cperiodic-east { type cyclicAMI; }
+
+    cperiodic-down { type cyclicAMI; }
+    cperiodic-up { type cyclicAMI; }
+}
+
+
+// ************************************************************************* //

cell-simple/0/alpha.oil

@@ -0,0 +1,57 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      alpha.water;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    inlet
+    {
+        type            cyclicAMI;
+    }
+    outlet
+    {
+        type            cyclicAMI;
+    }
+    cperiodic-west
+    {
+        type            cyclicAMI;
+    }
+    cperiodic-east
+    {
+        type            cyclicAMI;
+    }
+    cperiodic-down
+    {
+        type            cyclicAMI;
+    }
+    cperiodic-up
+    {
+        type            cyclicAMI;
+    }
+    walls
+    {
+        type            constantAlphaContactAngle;
+        theta0          180;
+        limit           gradient;
+        value           uniform 0;
+    }
+}
+
+
+// ************************************************************************* //

cell-simple/0/p_rgh

@@ -0,0 +1,48 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    object      p_rgh;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+UTranslationDirectionNbTimes    1;
+unitCellPressureJump            #calc "22765.0";
+pressureJump                    #calc "-$UTranslationDirectionNbTimes * $unitCellPressureJump";
+rho0                            1000.0;
+unitCellSize                    2.828427e-05;
+
+unitCellPassageTime             #calc "$unitCellSize / sqrt(2.0 * abs($unitCellPressureJump) / $rho0)";
+
+dimensions      [1 -1 -2 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    inlet  {
+        type uniformJumpAMI;
+        patchType cyclicAMI;
+        jumpTable table
+        (
+            (0 0.0)
+            ($unitCellPassageTime $pressureJump) //0.6500712769900825)
+        );
+        value $internalField;
+    }
+    outlet  { type uniformJumpAMI; patchType cyclicAMI; value $internalField; }
+    walls { type fixedFluxPressure; }
+    cperiodic-west { type cyclicAMI; }
+    cperiodic-east { type cyclicAMI; }
+    cperiodic-down { type cyclicAMI; }
+    cperiodic-up { type cyclicAMI; }
+}
+
+// ************************************************************************* //

cell-simple/Allclean

@@ -0,0 +1,17 @@
+#!/bin/sh
+
+# Run from this directory
+cd "${0%/*}" || exit 1
+
+# Source tutorial clean functions
+. "$WM_PROJECT_DIR/bin/tools/CleanFunctions"
+
+# Delete the following:
+# - time directories
+# - constant/polyMesh directory
+# - postProcessing and VTK directories
+# - log files
+# - field files with a ".orig" backup
+cleanCase
+
+#------------------------------------------------------------------------------

cell-simple/Allrun

@@ -0,0 +1,25 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+./Allclean
+
+# Source tutorial run functions
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+# Get application directory
+application=$(getApplication)
+
+runApplication ideasUnvToFoam simple_orig.unv
+runApplication transformPoints "scale=(1e-5 1e-5 1e-5)"
+#runApplication createPatch -overwrite
+cp -f logs/boundary constant/polyMesh/
+
+
+#- Run parallel
+#runApplication decomposePar #-cellDist
+runApplication setFields
+runApplication renumberMesh
+runApplication $application
+#runApplication reconstructPar
+
+#------------------------------------------------------------------------------

cell-simple/constant/g

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       uniformDimensionedVectorField;
+    location    "constant";
+    object      g;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -2 0 0 0 0];
+value           (0 0 -9.81);
+
+
+// ************************************************************************* //

cell-simple/constant/momentumTransport

@@ -0,0 +1,20 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      momentumTransport;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType  laminar;
+
+
+// ************************************************************************* //

cell-simple/constant/phaseProperties

@@ -0,0 +1,21 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      phaseProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+phases (oil water);
+
+sigma            0.046; //0.07;
+
+// ************************************************************************* //

cell-simple/constant/physicalProperties.oil

@@ -0,0 +1,26 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      physicalProperties.oil;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+
+viscosityModel constant;
+
+nu 1e-5;
+rho 900;
+
+//nu [0 2 -1 0 0 0 0] 1e-5;
+//rho [1 -3 0 0 0 0 0] 900;
+
+// ************************************************************************* //

cell-simple/constant/physicalProperties.water

@@ -0,0 +1,27 @@
+
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      physicalProperties.water;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+
+viscosityModel constant;
+
+nu 1e-06;
+rho 1000;
+
+//nu [0 2 -1 0 0 0 0] 1e-06;
+//rho [1 -3 0 0 0 0 0] 1000;
+
+// ************************************************************************* //

cell-simple/constant/polyMesh/boundary

@@ -0,0 +1,115 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       polyBoundaryMesh;
+    location    "constant/polyMesh";
+    object      boundary;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+7
+(
+    inlet
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.001;
+neighbourPatch  outlet;
+transformType   translational;
+separation      (-2.828427e-05 0 0);
+method          faceAreaWeightAMI;
+
+        nFaces          2444;
+        startFace       473972;
+    }
+    outlet
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.001;
+neighbourPatch  inlet;
+transformType   translational;
+separation      (2.828427e-05 0 0);
+method          faceAreaWeightAMI;
+
+        nFaces          2476;
+        startFace       476416;
+    }
+
+    cperiodic-west
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.0001;
+neighbourPatch  cperiodic-east;
+transformType   translational;
+separation      (0 -2.828427e-5 0);
+method          faceAreaWeightAMI;
+
+        nFaces          2042;
+        startFace       478892;
+    }
+    cperiodic-east
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.0001;
+neighbourPatch  cperiodic-west;
+transformType   translational;
+separation      (0 2.828427e-5 0);
+method          faceAreaWeightAMI;
+
+        nFaces          2092;
+        startFace       480934;
+    }
+
+    cperiodic-down
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.001;
+neighbourPatch  cperiodic-up;
+transformType   translational;
+separation      (0 0 -2e-5);
+method          faceAreaWeightAMI;
+
+        nFaces          860;
+        startFace       483026;
+    }
+    cperiodic-up
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.001;
+neighbourPatch  cperiodic-down;
+transformType   translational;
+separation      (0 0 2e-5);
+method          faceAreaWeightAMI;
+
+        nFaces          858;
+        startFace       483886;
+    }
+    walls
+    {
+        type            wall;
+
+inGroups        List<word> 1(wall);
+        nFaces          15888;
+        startFace       484744;
+    }
+)
+
+// ************************************************************************* //

cell-simple/constant/transportProperties

@@ -0,0 +1,45 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  9
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      transportProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+phases (oil water);
+
+water
+{
+    transportModel  Newtonian;
+    nu              1e-06;
+    rho             1000;
+}
+
+oil
+{
+    transportModel  Newtonian;
+    nu              1e-5; //1.48e-05;
+    rho             900;
+}
+
+sigma            0.046;
+
+
+thetaAngle                      thetaAngle                      [0 0 0 0 0 0 0]     0;
+
+UTranslationDirectionNbTimes    UTranslationDirectionNbTimes    [0 0 0 0 0 0 0]     0;
+VTranslationDirectionNbTimes    VTranslationDirectionNbTimes    [0 0 0 0 0 0 0]     0;
+WTranslationDirectionNbTimes    WTranslationDirectionNbTimes    [0 0 0 0 0 0 0]     0;
+
+unitCellSize                    unitCellSize                    [0 1 0 0 0 0 0]     0;
+L                               L                               [0 1 0 0 0 0 0]     0;
+
+// ************************************************************************* //

cell-simple/logs/boundary

@@ -0,0 +1,115 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       polyBoundaryMesh;
+    location    "constant/polyMesh";
+    object      boundary;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+7
+(
+    inlet
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.001;
+neighbourPatch  outlet;
+transformType   translational;
+separation      (-2.828427e-05 0 0);
+method          faceAreaWeightAMI;
+
+        nFaces          2444;
+        startFace       473972;
+    }
+    outlet
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.001;
+neighbourPatch  inlet;
+transformType   translational;
+separation      (2.828427e-05 0 0);
+method          faceAreaWeightAMI;
+
+        nFaces          2476;
+        startFace       476416;
+    }
+
+    cperiodic-west
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.0001;
+neighbourPatch  cperiodic-east;
+transformType   translational;
+separation      (0 -2.828427e-5 0);
+method          faceAreaWeightAMI;
+
+        nFaces          2042;
+        startFace       478892;
+    }
+    cperiodic-east
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.0001;
+neighbourPatch  cperiodic-west;
+transformType   translational;
+separation      (0 2.828427e-5 0);
+method          faceAreaWeightAMI;
+
+        nFaces          2092;
+        startFace       480934;
+    }
+
+    cperiodic-down
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.001;
+neighbourPatch  cperiodic-up;
+transformType   translational;
+separation      (0 0 -2e-5);
+method          faceAreaWeightAMI;
+
+        nFaces          860;
+        startFace       483026;
+    }
+    cperiodic-up
+    {
+
+type            cyclicAMI;
+inGroups        List<word> 1(cyclicAMI);
+matchTolerance  0.001;
+neighbourPatch  cperiodic-down;
+transformType   translational;
+separation      (0 0 2e-5);
+method          faceAreaWeightAMI;
+
+        nFaces          858;
+        startFace       483886;
+    }
+    walls
+    {
+        type            wall;
+
+inGroups        List<word> 1(wall);
+        nFaces          15888;
+        startFace       484744;
+    }
+)
+
+// ************************************************************************* //

cell-simple/system/controlDict

@@ -0,0 +1,92 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     interFoam;
+
+startFrom       latestTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         1.0;
+
+deltaT          1e-8;
+
+writeControl    adjustableRunTime;
+
+writeInterval   1e-7;
+
+purgeWrite      0;
+
+writeFormat     ascii;
+
+writePrecision  6;
+
+writeCompression off;
+
+timeFormat      general;
+
+timePrecision   6;
+
+runTimeModifiable yes;
+
+adjustTimeStep  yes;
+
+maxCo           0.2;
+maxAlphaCo      0.1;
+
+maxDeltaT       1;
+
+functions
+{
+    residuals
+    {
+        type            residuals;
+        libs            ("libutilityFunctionObjects.so");
+        fields  (p_rgh U);
+        writeResidualFields yes;
+        writeControl    timeStep;
+        writeInterval   1;
+    }
+    saturation1
+    {
+        type            saturation;
+        libs            ("multiphaseFunctionObjects.so");
+        writeControl    timeStep;
+        writeInterval   10;
+    }
+    volumetricFlowRate1
+    {
+        type volumetricFlowRate;
+        libs ("multiphaseFunctionObjects.so");
+        writeControl    timeStep;
+        writeInterval   10;
+        phi phi;
+        patches (inlet outlet cperiodic-west cperiodic-east cperiodic-down cperiodic-up);
+    }
+    massFlowRate1
+    {
+        type massFlowRate;
+        libs ("multiphaseFunctionObjects.so");
+        writeControl    timeStep;
+        writeInterval   10;
+        phi phi;
+        patches (inlet outlet cperiodic-west cperiodic-east cperiodic-down cperiodic-up);
+    }
+}
+
+// ************************************************************************* //

cell-simple/system/createPatchDict

@@ -0,0 +1,147 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version: 10                                      |
+|   \\  /    A nd           |                                                 |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    class dictionary;
+    format ascii;
+    location "system";
+    object createPatchDict;
+    version 2.0;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+// Do a synchronisation of coupled points after creation of any patches.
+// Note: this does not work with points that are on multiple coupled patches
+//       with transformations (i.e. cyclics).
+pointSync false;
+patches
+(
+    /*{
+        constructFrom patches;
+        name outlet;
+        patchInfo
+        {
+            inGroups (outlet);
+            type patch;
+        }
+        patches (periodic-west);
+    }
+
+    {
+        constructFrom patches;
+        name inlet;
+        patchInfo
+        {
+            inGroups (inlet);
+            type patch;
+        }
+        patches (periodic-east);
+    }*/
+
+    {
+        constructFrom patches;
+        name inlet;
+        patchInfo
+        {
+            inGroups (cyclicAMI);
+            type cyclicAMI;
+            transformType translational;
+            neighbourPatch outlet;
+            matchTolerance  1e-3;
+            separation      (2.828427e-5 0 0);
+        }
+        patches (periodic-north);
+    }
+
+    {
+        constructFrom patches;
+        name outlet;
+        patchInfo
+        {
+            inGroups (cyclicAMI);
+            type cyclicAMI;
+            transformType translational;
+            neighbourPatch inlet;
+            matchTolerance  1e-3;
+            separation      (-2.828427e-5 0 0);
+        }
+        patches (periodic-south);
+    }
+
+    {
+        constructFrom patches;
+        name walls;
+        patchInfo
+        {
+            inGroups (wall);
+            type wall;
+        }
+        patches (wall);
+    }
+
+    {
+        constructFrom patches;
+        name cperiodic-west;
+        patchInfo
+        {
+            inGroups (cyclicAMI);
+            type cyclicAMI;
+            transformType translational;
+            neighbourPatch cperiodic-east;
+            matchTolerance  1e-4;
+            separation      (0 -2.828427e-5 0);
+        }
+        patches (periodic-west);
+    }
+
+    {
+        constructFrom patches;
+        name cperiodic-east;
+        patchInfo
+        {
+            inGroups (cyclicAMI);
+            type cyclicAMI;
+            transformType translational;
+            neighbourPatch cperiodic-west;
+            matchTolerance  1e-4;
+            separation      (0 2.828427e-5 0);
+        }
+        patches (periodic-east);
+    }
+
+    {
+        constructFrom patches;
+        name cperiodic-down;
+        patchInfo
+        {
+            inGroups (cyclicAMI);
+            type cyclicAMI;
+            transformType translational;
+            neighbourPatch cperiodic-up;
+            matchTolerance  1e-3;
+            separation      (0 0 -2e-5);
+        }
+        patches (periodic-down);
+    }
+
+    {
+        constructFrom patches;
+        name cperiodic-up;
+        patchInfo
+        {
+            inGroups (cyclicAMI);
+            type cyclicAMI;
+            transformType translational;
+            neighbourPatch cperiodic-down;
+            matchTolerance  1e-3;
+            separation      (0 0 2e-5);
+        }
+        patches (periodic-up);
+    }
+);
+
+// ************************************************************************* //

cell-simple/system/decomposeParDict

@@ -0,0 +1,49 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      decomposeParDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+numberOfSubdomains 4;
+
+method          scotch;
+
+simpleCoeffs
+{
+    n               (2 2 1);
+    delta           0.001;
+}
+
+hierarchicalCoeffs
+{
+    n               ( 2 2 1 );
+    delta           0.001;
+    order           xyz;
+}
+
+manualCoeffs
+{
+    dataFile        "";
+}
+
+scotchCoeffs
+{
+    processorWeights ( );
+}
+
+distributed     no;
+
+roots           ( );
+
+
+// ************************************************************************* //

cell-simple/system/fvSchemes

@@ -0,0 +1,52 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSchemes;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+ddtSchemes
+{
+    default         Euler;
+}
+
+gradSchemes
+{
+    default         Gauss linear;
+}
+
+divSchemes
+{
+    default none;
+    div(rhoPhi,U)  bounded Gauss linearUpwind grad(U);
+    div(phi,alpha)  Gauss interfaceCompression vanLeer 1;
+    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
+    //div(nonlinearStress) Gauss linear ;
+}
+
+laplacianSchemes
+{
+    default         Gauss linear corrected;
+}
+
+interpolationSchemes
+{
+    default         linear;
+}
+
+snGradSchemes
+{
+    default         corrected;
+}
+
+
+// ************************************************************************* //

cell-simple/system/fvSolution

@@ -0,0 +1,117 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSolution;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solvers
+{
+    "alpha.oil.*"
+    {
+        nAlphaCorr      2;
+        nAlphaSubCycles 2;
+        cAlpha          1;
+        MULESCorr       yes;
+        nLimiterIter    5;
+
+        solver          smoothSolver;
+        smoother        symGaussSeidel;
+        tolerance       1e-8;
+        relTol          0;
+    }
+
+    "pcorr.*"
+    {
+        solver          GAMG;//PCG;
+        //preconditioner  DIC;
+        tolerance       1e-8;
+        relTol          1e-3;
+        smoother        DICGaussSeidel;
+        cacheAgglomeration true;
+        nCellInCoarsestLevel 10;
+        agglomerator faceAreaPair;
+        mergeLevels 1;
+    }
+
+    p_rgh
+    {
+        solver          GAMG;//PCG;
+        //preconditioner  DIC;
+        tolerance       1e-8;
+        relTol          1e-3;
+        smoother        DICGaussSeidel;
+        cacheAgglomeration true;
+        nCellInCoarsestLevel 10;
+        agglomerator faceAreaPair;
+        mergeLevels 1;
+    }
+
+    p_rghFinal
+    {
+        $p_rgh;
+        tolerance 1e-8;
+        relTol          0;
+    }
+
+    U
+    {
+        //nSweeps 2;
+        solver          smoothSolver;
+        smoother        symGaussSeidel;
+        tolerance       1e-08;
+        relTol          0;
+    }
+
+    UFinal
+    {
+        $U;
+        relTol 0;
+    }
+}
+
+PIMPLE
+{
+    momentumPredictor   no;
+    nOuterCorrectors    5;
+    nCorrectors         2;
+    nNonOrthogonalCorrectors 2;
+    pRefCell 0;
+    pRefValue 0;
+
+    outerCorrectorResidualControl
+    {
+        p_rgh
+        {
+            tolerance 1e-5;
+            relTol 0;
+        }
+    }
+
+}
+
+relaxationFactors
+{
+    fields
+    {
+        p_rgh 1;
+        p_rghFinal 1;
+    }
+
+    equations
+    {
+        ".*" 1;
+    }
+}
+
+
+// ************************************************************************* //

cell-simple/system/residuals

@@ -0,0 +1,19 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+Description
+    For specified fields, writes out the initial residuals for the first
+    solution of each time step; for non-scalar fields (e.g. vectors), writes
+    the largest of the residuals for each component (e.g. x, y, z).
+
+\*---------------------------------------------------------------------------*/
+
+#includeEtc "residuals.cfg"
+
+fields  (<fieldNames>);
+writeResidualFields yes;
+// ************************************************************************* //

cell-simple/system/residuals.cfg

@@ -0,0 +1,15 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+
+type            solverInfo;
+libs            ("libutilityFunctionObjects.so");
+
+writeControl    timeStep;
+writeInterval   1;
+
+// ************************************************************************* //

cell-simple/system/saturation

@@ -0,0 +1,16 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+
+type            saturation;
+libs            ("multiphaseFunctionObjects.so");
+
+writeControl    timeStep;
+writeInterval   1;
+
+
+// ************************************************************************* //

cell-simple/system/saturation.cfg

@@ -0,0 +1,15 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+
+type            saturation;
+libs            ("multiphaseFunctionObjects.so");
+
+writeControl    timeStep;
+writeInterval   1;
+
+// ************************************************************************* //

cell-simple/system/setFieldsDict

@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      setFieldsDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+defaultFieldValues
+(
+    volScalarFieldValue alpha.water 1
+    volScalarFieldValue alpha.oil 0
+);
+
+regions
+(
+    /*boxToCell
+    {
+        box (-0.001414 -0.001414 -0.001) (0.001414 0.001414 0.001);
+        fieldValues
+        (
+            volScalarFieldValue alpha.water 1
+        );
+    }*/
+    sphereToCell
+    {
+        centre (-1.414e-05 0 0); //(-1.0e-5 0 0);//(-1.414e-5 0 0);
+        radius 0.7e-5; // 0.722e-5
+        fieldValues
+        (
+            volScalarFieldValue alpha.water 0
+            volScalarFieldValue alpha.oil 1
+        );
+    }
+    sphereToCell
+    {
+        centre (1.414e-05 0 0); //(-1.0e-5 0 0);//(-1.414e-5 0 0);
+        radius 0.7e-5;
+        fieldValues
+        (
+            volScalarFieldValue alpha.water 0
+            volScalarFieldValue alpha.oil 1
+        );
+    }
+);
+
+
+// ************************************************************************* //