Non-Conformal Coupled (NCC): Conversion of tutorials from AMI to NCC

The following examples in the tutorials ($FOAM_TUTORIALS) directory have
been converted from using AMI to the new NCC system:

+ compressible/rhoPimpleFoam/RAS/annularThermalMixer
+ incompressible/pimpleFoam/RAS/propeller
+ lagrangian/particleFoam/mixerVessel2D (formerly mixerVesselAMI2D)
+ multiphase/interFoam/RAS/mixerVessel
+ multiphase/interFoam/RAS/propeller
+ multiphase/multiphaseEulerFoam/laminar/mixerVessel2D (formerly mixerVesselAMI2D)

The following tutorial has been converted from using ACMI:

+ incompressible/pimpleFoam/RAS/oscillatingInlet

The following tutorial has been converted from using Repeat AMI:

+ incompressible/pimpleFoam/RAS/impeller

The following tutorial has been added to demonstrate NCC's ability to
create a sufficiently conservative solution in a closed domain to
maintain phase fraction boundedness:

+ multiphase/interFoam/laminar/mixerVessel2D

The following tutorials have been added to demonstrate NCC's ability to
simulate partially overlapping couples on curved surfaces:

+ incompressible/pimpleFoam/RAS/ballValve
+ multiphase/compressibleInterFoam/RAS/ballValve

The following tutorial has been added to provide a simple comparison of
the conservation behaviour of AMI and NCC:

+ incompressible/pimpleFoam/laminar/nonConformalChannel

The following tutorial has been removed, as there were sufficiently many
examples involving this geometry:

+ incompressible/pimpleFoam/laminar/mixerVesselAMI2D

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/0/T

@@ -9,6 +9,7 @@ FoamFile
 {
     format      ascii;
     class       volScalarField;
+    location    "0";
     object      T;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -21,6 +22,8 @@ internalField   uniform $!outerInlet/T;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     innerInlet
     {
         type            fixedValue;
@@ -48,7 +51,10 @@ boundaryField
         $staticWalls;
     }
 
-    #includeEtc "caseDicts/setConstraintTypes"
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 // ************************************************************************* //

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/0/U

@@ -22,6 +22,8 @@ internalField   uniform (0 0 0);
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     innerInlet
     {
         type            fixedValue;
@@ -49,7 +51,10 @@ boundaryField
         $!movingWall/U;
     }
 
-    #includeEtc "caseDicts/setConstraintTypes"
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 // ************************************************************************* //

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/0/alphat

@@ -22,6 +22,8 @@ internalField   uniform 0;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     inlet
     {
         $!calculated/alphat;
@@ -42,8 +44,15 @@ boundaryField
         $movingWalls;
     }
 
-    #includeEtc "caseDicts/setConstraintTypes"
+    nonCouple
+    {
+        type            zeroGradient;
     }
 
+    nonCouple
+    {
+        type            zeroGradient;
+    }
+}
 
 // ************************************************************************* //

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/0/epsilon

@@ -22,6 +22,8 @@ internalField   uniform $!innerInlet/epsilon;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     innerInlet
     {
         type            fixedValue;
@@ -49,7 +51,15 @@ boundaryField
         $staticWalls;
     }
 
-    #includeEtc "caseDicts/setConstraintTypes"
+    nonCouple
+    {
+        type            zeroGradient;
+    }
+
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 // ************************************************************************* //

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/0/k

@@ -22,6 +22,8 @@ internalField   uniform $!innerInlet/k;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     innerInlet
     {
         type            fixedValue;
@@ -49,7 +51,10 @@ boundaryField
         $staticWalls;
     }
 
-    #includeEtc "caseDicts/setConstraintTypes"
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 // ************************************************************************* //

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/0/nut

@@ -22,6 +22,8 @@ internalField   uniform 0;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     inlet
     {
         $!calculated/nut;
@@ -42,7 +44,10 @@ boundaryField
         $movingWalls;
     }
 
-    #includeEtc "caseDicts/setConstraintTypes"
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 // ************************************************************************* //

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/0/p

@@ -22,6 +22,8 @@ internalField   uniform 1e5;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     inlet
     {
         $!inlet/p;
@@ -42,7 +44,10 @@ boundaryField
         $staticWalls;
     }
 
-    #includeEtc "caseDicts/setConstraintTypes"
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 // ************************************************************************* //

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/Allrun

@@ -11,6 +11,7 @@ runApplication snappyHexMesh -overwrite
 
 runApplication createBaffles -overwrite
 runApplication splitBaffles -overwrite
+runApplication createNonConformalCouples -overwrite nonCouple1 nonCouple2
 
 runApplication $(getApplication)
 

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/constant/dynamicMeshDict

@@ -24,7 +24,7 @@ mover
 
     motionSolver    solidBody;
 
-    cellZone        rotatingZone;
+    cellZone        rotating;
 
     solidBodyMotionFunction  rotatingMotion;
 

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/system/controlDict

@@ -32,7 +32,7 @@ writeInterval   0.02;
 
 purgeWrite      0;
 
-writeFormat     ascii;
+writeFormat     binary;
 
 writePrecision  6;
 

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/system/createBafflesDict

@@ -20,21 +20,21 @@ baffles
     baffleFaces
     {
         type        faceZone;
-        zoneName    rotatingZone;
+        zoneName    couple;
 
         patches
         {
             master
             {
-                name            AMI1;
+                name            nonCouple1;
-                type            cyclicAMI;
+                type            patch;
-                neighbourPatch  AMI2;
+                inGroups        (nonCouple);
             }
             slave
             {
-                $master;
+                name            nonCouple2;
-                name            AMI2;
+                type            patch;
-                neighbourPatch  AMI1;
+                inGroups        (nonCouple);
             }
         }
     }

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/system/fvSchemes

@@ -52,4 +52,5 @@ snGradSchemes
     default         corrected;
 }
 
+
 // ************************************************************************* //

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/system/fvSolution

@@ -19,9 +19,7 @@ solvers
     p
     {
         solver          GAMG;
-
         smoother        GaussSeidel;
-
         tolerance       1e-6;
         relTol          0.01;
     }
@@ -39,6 +37,14 @@ solvers
         relTol          0;
     }
 
+    MeshPhi
+    {
+        solver          smoothSolver;
+        smoother        symGaussSeidel;
+        tolerance       1e-2;
+        relTol          0;
+    }
+
     "(rho|U|h|k|epsilon|omega)"
     {
         solver          smoothSolver;
@@ -52,13 +58,13 @@ solvers
         $U;
         relTol          0;
     }
-
 }
 
 PIMPLE
 {
     momentumPredictor   yes;
     transonic           no;
+    correctMeshPhi      yes;
     nOuterCorrectors    3;
     nCorrectors         1;
     nNonOrthogonalCorrectors 0;

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/system/snappyHexMeshDict

@@ -19,10 +19,10 @@ addLayers       false;
 
 geometry
 {
-    AMI
+    couple
     {
         type        triSurfaceMesh;
-        file        "AMI.obj";
+        file        "couple.obj";
     }
 
     innerInlet
@@ -85,7 +85,7 @@ castellatedMeshControls
     features
     (
         {
-            file        "AMI.eMesh";
+            file        "couple.eMesh";
             level       2;
         }
         {
@@ -108,11 +108,11 @@ castellatedMeshControls
 
     refinementSurfaces
     {
-        AMI
+        couple
         {
             level       (2 2);
-            faceZone    rotatingZone;
+            faceZone    couple;
-            cellZone    rotatingZone;
+            cellZone    rotating;
             mode        inside;
         }
 

tutorials/compressible/rhoPimpleFoam/RAS/annularThermalMixer/system/surfaceFeaturesDict

@@ -15,7 +15,7 @@ FoamFile
 
 surfaces
 (
-    "AMI.obj"
+    "couple.obj"
     "shaft.obj"
     "wall.obj"
     "statorBlades.obj"

tutorials/incompressible/pimpleFoam/RAS/ballValve/0/U

@@ -0,0 +1,46 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volVectorField;
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            pressureInletOutletVelocity;
+        value           $internalField;
+    }
+    lowerOutlet
+    {
+        type            pressureInletOutletVelocity;
+        value           $internalField;
+    }
+    upperOutlet
+    {
+        type            pressureInletOutletVelocity;
+        value           $internalField;
+    }
+    wall
+    {
+        type            movingWallVelocity;
+        value           uniform (0 0 0);
+    }
+}
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/ballValve/0/epsilon

@@ -20,46 +20,31 @@ internalField   uniform 1.8e-3;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     inlet
     {
-        type            fixedValue;
+        type            turbulentMixingLengthDissipationRateInlet;
+        mixingLength    0.1;
         value           $internalField;
     }
-    outlet
+    lowerOutlet
     {
-        type            inletOutlet;
+        type            turbulentMixingLengthDissipationRateInlet;
-        inletValue      $internalField;
+        mixingLength    0.1;
         value           $internalField;
     }
-    walls
+    upperOutlet
+    {
+        type            turbulentMixingLengthDissipationRateInlet;
+        mixingLength    0.1;
+        value           $internalField;
+    }
+    wall
     {
         type            epsilonWallFunction;
         value           $internalField;
     }
-    defaultFaces
-    {
-        type            empty;
-    }
-    ACMI1_blockage
-    {
-        type            epsilonWallFunction;
-        value           $internalField;
-    }
-    ACMI1_couple
-    {
-        type            cyclicACMI;
-        value           $internalField;
-    }
-    ACMI2_blockage
-    {
-        type            epsilonWallFunction;
-        value           $internalField;
-    }
-    ACMI2_couple
-    {
-        type            cyclicACMI;
-        value           $internalField;
-    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/ballValve/0/k

@@ -0,0 +1,51 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      k;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 3.75e-3;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            turbulentIntensityKineticEnergyInlet;
+        intensity       0.05;
+        value           $internalField;
+    }
+    lowerOutlet
+    {
+        type            turbulentIntensityKineticEnergyInlet;
+        intensity       0.05;
+        value           $internalField;
+    }
+    upperOutlet
+    {
+        type            turbulentIntensityKineticEnergyInlet;
+        intensity       0.05;
+        value           $internalField;
+    }
+    wall
+    {
+        type            kqRWallFunction;
+        value           $internalField;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/ballValve/0/nut

@@ -9,49 +9,40 @@ FoamFile
 {
     format      ascii;
     class       volScalarField;
-    object      p;
+    location    "0";
+    object      nut;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-dimensions          [1 -1 -2 0 0 0 0];
+dimensions      [0 2 -1 0 0 0 0];
 
-internalField       uniform 1e5;
+internalField   uniform 0;
 
 boundaryField
 {
-    rotor
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
     {
         type            calculated;
         value           $internalField;
     }
-
+    lowerOutlet
-    stator
     {
         type            calculated;
         value           $internalField;
     }
-
+    upperOutlet
-    front
     {
-        type            empty;
+        type            calculated;
-    }
-
-    back
-    {
-        type            empty;
-    }
-
-    AMI1
-    {
-        type            cyclicAMI;
         value           $internalField;
     }
-
+    wall
-    AMI2
     {
-        type            cyclicAMI;
+        type            nutkWallFunction;
         value           $internalField;
     }
 }
 
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/ballValve/0/omega

@@ -0,0 +1,51 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      omega;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 -1 0 0 0 0];
+
+internalField   uniform 5.33333;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            turbulentMixingLengthFrequencyInlet;
+        mixingLength    0.1;
+        value           $internalField;
+    }
+    lowerOutlet
+    {
+        type            turbulentMixingLengthFrequencyInlet;
+        mixingLength    0.1;
+        value           $internalField;
+    }
+    upperOutlet
+    {
+        type            turbulentMixingLengthFrequencyInlet;
+        mixingLength    0.1;
+        value           $internalField;
+    }
+    wall
+    {
+        type            omegaWallFunction;
+        value           $internalField;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/ballValve/0/p

@@ -0,0 +1,46 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            totalPressure;
+        p0              uniform 10;
+    }
+    lowerOutlet
+    {
+        type            entrainmentPressure;
+        p0              uniform 0;
+    }
+    upperOutlet
+    {
+        type            entrainmentPressure;
+        p0              uniform 0;
+    }
+    wall
+    {
+        type            fixedFluxPressure;
+    }
+}
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/ballValve/0/s

@@ -0,0 +1,49 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      s;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+    lowerOutlet
+    {
+        type            inletOutlet;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+    upperOutlet
+    {
+        type            inletOutlet;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+    wall
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/ballValve/Allclean

@@ -0,0 +1,7 @@
+#!/bin/sh
+
+cd "${0%/*}" || exit 1
+
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase && rm -rf 0/reconstruct\(phi\) constant/geometry flowRate.eps error.eps

tutorials/incompressible/pimpleFoam/RAS/ballValve/Allrun

@@ -0,0 +1,31 @@
+#!/bin/sh
+
+cd ${0%/*} || exit 1
+
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+# Copy torus surface from resources directory
+mkdir constant/geometry && cp $FOAM_TUTORIALS/resources/geometry/ballValve-torus.obj.gz constant/geometry
+
+# Mesh
+runApplication blockMesh -dict $FOAM_TUTORIALS/resources/blockMesh/ballValve
+runApplication topoSet
+runApplication transformPoints -pointSet ball "Rz=-45"
+
+# Create couples
+runApplication createNonConformalCouples -overwrite pipeNonCouple ballNonCouple
+
+# Decompose
+runApplication decomposePar -cellDist
+
+# Set the value of the scalar transport field
+runParallel setFields
+
+# Run
+runParallel $(getApplication);
+
+# Reconstruct
+runApplication reconstructPar -newTimes
+
+# Plot the flow rates
+./createGraphs

tutorials/incompressible/pimpleFoam/RAS/ballValve/constant/dynamicMeshDict

@@ -22,13 +22,13 @@ mover
 
     motionSolver    solidBody;
 
-    cellZone        rotor;
+    cellZone        ball;
 
-    solidBodyMotionFunction  rotatingMotion;
+    solidBodyMotionFunction oscillatingRotatingMotion;
 
     origin        (0 0 0);
-    axis          (0 0 1);
+    amplitude     (0 0 45);
-    omega         6.2832; // rad/s
+    omega         6.283185307179586;
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/ballValve/constant/momentumTransport

@@ -0,0 +1,29 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      momentumTransport;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+simulationType  RAS;
+
+RAS
+{
+    model           kOmegaSST; //kEpsilon;
+
+    turbulence      on;
+
+    printCoeffs     on;
+}
+
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/ballValve/createGraphs

@@ -0,0 +1,30 @@
+#!/bin/sh
+
+cd ${0%/*} || exit 1
+
+inletFlowRate=postProcessing/inletFlowRate/0/surfaceFieldValue.dat
+lowerOutletFlowRate=postProcessing/lowerOutletFlowRate/0/surfaceFieldValue.dat
+upperOutletFlowRate=postProcessing/upperOutletFlowRate/0/surfaceFieldValue.dat
+
+gnuplot << EOF
+set terminal postscript eps color enhanced
+set output "flowRate.eps"
+set xlabel "Time (s)"
+set ylabel "Flow Rate (m^3/s)"
+plot \
+    "$inletFlowRate" us 1:(-\$2) every 10 w p t "Inlet", \
+    "$lowerOutletFlowRate" us 1:2 w l t "Lower Outlet", \
+    "$upperOutletFlowRate" us 1:2 w l t "Upper Outlet", \
+    "< paste $lowerOutletFlowRate $upperOutletFlowRate" \
+    us 1:(\$2+\$4) w l t "Total Outlet"
+EOF
+
+gnuplot << EOF
+set terminal postscript eps color enhanced
+set output "error.eps"
+set xlabel "Time (s)"
+set ylabel "Inlet-Outlet Conservation Error (%)"
+plot \
+    "< paste $inletFlowRate $lowerOutletFlowRate $upperOutletFlowRate" \
+    us 1:(100*(\$2+\$4+\$6)/\$2) w l t ''
+EOF

tutorials/incompressible/pimpleFoam/RAS/ballValve/system/controlDict

@@ -16,19 +16,19 @@ FoamFile
 
 application     pimpleFoam;
 
-startFrom       startTime;
+startFrom       latestTime;
 
 startTime       0;
 
 stopAt          endTime;
 
-endTime         5;
+endTime         1;
 
 deltaT          1e-3;
 
 writeControl    adjustableRunTime;
 
-writeInterval   0.05;
+writeInterval   0.01;
 
 purgeWrite      0;
 
@@ -50,37 +50,13 @@ maxCo           0.5;
 
 functions
 {
-    cartesianToCylindrical
+    #includeFunc patchFlowRate(funcName=inletFlowRate, patch=inlet)
-    {
+    #includeFunc patchFlowRate(funcName=lowerOutletFlowRate, patch=lowerOutlet)
-        type        cylindrical;
+    #includeFunc patchFlowRate(funcName=upperOutletFlowRate, patch=upperOutlet)
-        libs        ("libfieldFunctionObjects.so");
+    #includeFunc scalarTransport(s)
-
+    #includeFunc reconstruct(phi)
-        origin      (0 0 0);
+    #includeFunc streamlinesPatch(patch=inlet, nPoints=100, U=reconstruct(phi))
-        axis        (0 0 1);
+    #includeFunc wallsSurface
-
-        field       U;
-
-        writeControl    outputTime;
-        writeInterval   1;
-    }
-
-    #includeFunc fieldAverage(cylindrical(U))
-
-    cylindricalToCartesian
-    {
-        type        cylindrical;
-        libs        ("libfieldFunctionObjects.so");
-
-        origin      (0 0 0);
-        axis        (0 0 1);
-
-        field       cylindrical(U)Mean;
-        toCartesian true;
-        result      UMean;
-
-        writeControl    outputTime;
-        writeInterval   1;
-    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/ballValve/system/decomposeParDict

@@ -0,0 +1,28 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    object      decomposeParDict;
+}
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+numberOfSubdomains 8;
+
+method          scotch;
+
+hierarchicalCoeffs
+{
+    n               (2 2 2);
+    order           xyz;
+}
+
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/ballValve/system/fvSchemes

@@ -0,0 +1,64 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSchemes;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+ddtSchemes
+{
+    default         Euler;
+}
+
+gradSchemes
+{
+    default         Gauss linear;
+    grad(U)         cellLimited Gauss linear 1;
+    grad(k)         cellLimited Gauss linear 1;
+    grad(epsilon)   cellLimited Gauss linear 1;
+    grad(omega)     cellLimited Gauss linear 1;
+}
+
+divSchemes
+{
+    default         none;
+    div(phi,U)      Gauss linearUpwind grad(U);
+    div(phi,k)      Gauss limitedLinear 1;
+    div(phi,epsilon) Gauss limitedLinear 1;
+    div(phi,omega)  Gauss limitedLinear 1;
+    div(phi,s)      Gauss limitedLinear 1;
+    div((nuEff*dev2(T(grad(U))))) Gauss linear;
+}
+
+laplacianSchemes
+{
+    default         Gauss linear corrected;
+}
+
+interpolationSchemes
+{
+    default         linear;
+}
+
+snGradSchemes
+{
+    default         corrected;
+}
+
+wallDist
+{
+    method          meshWave;
+    correctWalls    false;
+}
+
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/ballValve/system/fvSolution

@@ -20,8 +20,6 @@ solvers
     {
         solver          GAMG;
         smoother        GaussSeidel;
-        cacheAgglomeration no;
-
         tolerance       0.02;
         relTol          0;
     }
@@ -29,37 +27,27 @@ solvers
     p
     {
         $pcorr;
-        tolerance       1e-06;
+        tolerance       1e-6;
         relTol          0.01;
     }
 
     pFinal
     {
         $p;
-        tolerance       1e-06;
         relTol          0;
     }
 
-    U
+    "(U|k|epsilon|omega|s)"
     {
-        solver          smoothSolver;
+        solver          PBiCGStab;
-        smoother        GaussSeidel;
+        preconditioner  DILU;
-        tolerance       1e-05;
+        tolerance       1e-5;
         relTol          0.01;
     }
 
-    UFinal
+    "(U|k|epsilon|omega|s)Final"
     {
         $U;
-        tolerance       1e-06;
-        relTol          0;
-    }
-
-    cellMotionUx
-    {
-        solver          PCG;
-        preconditioner  DIC;
-        tolerance       1e-08;
         relTol          0;
     }
 }
@@ -67,12 +55,10 @@ solvers
 PIMPLE
 {
     correctPhi          yes;
+    correctMeshPhi      no;
     nOuterCorrectors    2;
     nCorrectors         1;
     nNonOrthogonalCorrectors 0;
-
-    pRefCell            0;
-    pRefValue           0;
 }
 
 relaxationFactors

tutorials/incompressible/pimpleFoam/RAS/ballValve/system/setFieldsDict

@@ -0,0 +1,37 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    object      setFieldsDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+defaultFieldValues
+(
+    volScalarFieldValue s 0
+);
+
+regions
+(
+    cylinderToFace
+    {
+        point1 (-0.51 0 0);
+        point2 (-0.49 0 0);
+        radius 0.05;
+
+        fieldValues
+        (
+            volScalarFieldValue s 1
+        );
+    }
+);
+
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/ballValve/system/topoSetDict

@@ -15,16 +15,20 @@ FoamFile
 
 actions
 (
-    // Get both sides of ami
-    // ~~~~~~~~~~~~~~~~~~~~~
-
-    // Get all faces in cellSet
     {
-        name    AMI;
-        type    faceSet;
         action  new;
-        source  patchToFace;
+        type    cellSet;
-        patch   "AMI.*";
+        name    ball;
+        source  zoneToCell;
+        zone    ball;
+    }
+    {
+        action  new;
+        type    pointSet;
+        name    ball;
+        source  cellToPoint;
+        set     ball;
+        option  all;
     }
 );
 

tutorials/incompressible/pimpleFoam/RAS/ballValve/system/wallsSurface

@@ -0,0 +1,31 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+
+type            surfaces;
+libs            ("libsampling.so");
+
+writeControl    writeTime;
+
+surfaceFormat   vtk;
+
+interpolationScheme cell;
+interpolate false;
+
+fields      ();
+
+surfaces
+(
+    walls
+    {
+        type        patch;
+        interpolate $interpolate;
+        patches     (pipeWalls pipeNonCouple ballWalls ballNonCouple);
+    }
+);
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/impeller/0/U

@@ -28,16 +28,22 @@ boundaryField
         value           $internalField;
     }
 
-    stationaryWalls
+    "(statorDuct|nonCoupleIn1|nonCoupleOut2)"
     {
         type            noSlip;
     }
 
-    movingWalls
+    "(rotorDuct|rotorBlades|nonCoupleIn2|nonCoupleOut1)"
     {
         type            movingWallVelocity;
         value           $internalField;
     }
+
+    nonCouple
+    {
+        type            movingWallSlipVelocity;
+        value           $internalField;
+    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/impeller/0/epsilon

@@ -34,6 +34,11 @@ boundaryField
         type            epsilonWallFunction;
         value           $internalField;
     }
+
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/impeller/0/k

@@ -34,6 +34,11 @@ boundaryField
         type            kqRWallFunction;
         value           $internalField;
     }
+
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/impeller/0/nut

@@ -33,6 +33,12 @@ boundaryField
         type            nutkWallFunction;
         value           uniform 0;
     }
+
+    nonCouple
+    {
+        type            calculated;
+        value           uniform 0;
+    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/impeller/0/p

@@ -24,7 +24,7 @@ boundaryField
 
     "(inlet|outlet)"
     {
-        type            totalPressure;
+        type            entrainmentPressure;
         p0              $internalField;
         value           $internalField;
     }
@@ -33,6 +33,11 @@ boundaryField
     {
         type            fixedFluxPressure;
     }
+
+    nonCouple
+    {
+        type            fixedFluxPressure;
+    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/impeller/Allclean

@@ -0,0 +1,7 @@
+#!/bin/sh
+
+cd "${0%/*}" || exit 1
+
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase && rm -f 0/cylindrical\(U\)

tutorials/incompressible/pimpleFoam/RAS/impeller/Allrun

@@ -6,8 +6,6 @@ cd "${0%/*}" || exit 1
 
 runApplication blockMesh
 
-runApplication decomposePar
+runApplication createNonConformalCouples -overwrite
 
-runParallel $(getApplication)
+runApplication $(getApplication)
-
-runApplication reconstructPar

tutorials/incompressible/pimpleFoam/RAS/impeller/system/blockMeshDict

@@ -83,15 +83,9 @@ edges
     arc 15 31 (2 $r4 0)
 );
 
-defaultPatch
-{
-    name movingWalls;
-    type wall;
-}
-
 boundary
 (
-    stationaryWalls
+    statorDuct
     {
         type wall;
         faces
@@ -103,6 +97,34 @@ boundary
         );
     }
 
+    rotorDuct
+    {
+        type wall;
+        faces
+        (
+            (0 1 17 16)
+            (1 2 18 17)
+            (0 5 21 16)
+            (5 10 26 21)
+            (6 7 23 22)
+            (6 11 27 22)
+        );
+    }
+
+    rotorBlades
+    {
+        type wall;
+        faces
+        (
+            (0 1 6 5)
+            (1 2 7 6)
+            (5 6 11 10)
+            (16 17 22 21)
+            (17 18 23 22)
+            (21 22 27 26)
+        );
+    }
+
     inlet
     {
         type patch;
@@ -141,22 +163,20 @@ boundary
         );
     }
 
-    cyclicRepeatAMIIn1
+    nonCoupleIn1
     {
-        type cyclicRepeatAMI;
+        type patch;
-        neighbourPatch cyclicRepeatAMIIn2;
+        inGroups (nonCouple);
-        transformPatch cyclicIn1;
         faces
         (
             (3 8 24 19)
         );
     }
 
-    cyclicRepeatAMIIn2
+    nonCoupleIn2
     {
-        type cyclicRepeatAMI;
+        type patch;
-        neighbourPatch cyclicRepeatAMIIn1;
+        inGroups (nonCouple);
-        transformPatch cyclicIn1;
         faces
         (
             (2 7 23 18)
@@ -183,22 +203,20 @@ boundary
         );
     }
 
-    cyclicRepeatAMIOut1
+    nonCoupleOut1
     {
-        type cyclicRepeatAMI;
+        type patch;
-        neighbourPatch cyclicRepeatAMIOut2;
+        inGroups (nonCouple);
-        transformPatch cyclicIn1;
         faces
         (
             (10 11 27 26)
         );
     }
 
-    cyclicRepeatAMIOut2
+    nonCoupleOut2
     {
-        type cyclicRepeatAMI;
+        type patch;
-        neighbourPatch cyclicRepeatAMIOut1;
+        inGroups (nonCouple);
-        transformPatch cyclicIn1;
         faces
         (
             (12 13 29 28)

tutorials/incompressible/pimpleFoam/RAS/impeller/system/controlDict

@@ -50,10 +50,23 @@ maxCo           1.0;
 
 functions
 {
-    #includeFunc patchFlowRate(patch=inlet)
+    #includeFunc patchFlowRate(patch=inlet, funcName=inletFlowRate)
-    #includeFunc patchFlowRate(patch=outlet)
+    #includeFunc patchFlowRate(patch=outlet, funcName=outletFlowRate)
-    #includeFunc triSurfaceVolumetricFlowRate(name=surface0.stl)
+
-    #includeFunc triSurfaceVolumetricFlowRate(name=surface1.stl)
+    cartesianToCylindrical
+    {
+        type        cylindrical;
+        libs        ("libfieldFunctionObjects.so");
+
+        origin      (0 0 0);
+        axis        (1 0 0);
+
+        field       U;
+
+        writeControl    outputTime;
+        writeInterval   1;
+    }
 }
 
+
 // ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/impeller/system/createNonConformalCouplesDict

@@ -0,0 +1,69 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      createNonConformalCouplesDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+nonConformalCoupleIn0
+{
+    patches         (nonCoupleIn1 nonCoupleIn2);
+    transform       none;
+}
+
+nonConformalCoupleIn30
+{
+    $nonConformalCoupleIn0;
+    transform       rotational;
+    rotationAxis    (-1 0 0);
+    rotationCentre  (0 0 0);
+    rotationAngle   30;
+}
+
+nonConformalCoupleIn60 { $nonConformalCoupleIn30; rotationAngle 60; }
+nonConformalCoupleIn90 { $nonConformalCoupleIn30; rotationAngle 90; }
+nonConformalCoupleIn120 { $nonConformalCoupleIn30; rotationAngle 120; }
+nonConformalCoupleIn150 { $nonConformalCoupleIn30; rotationAngle 150; }
+nonConformalCoupleIn180 { $nonConformalCoupleIn30; rotationAngle 180; }
+nonConformalCoupleIn210 { $nonConformalCoupleIn30; rotationAngle 210; }
+nonConformalCoupleIn240 { $nonConformalCoupleIn30; rotationAngle 240; }
+nonConformalCoupleIn270 { $nonConformalCoupleIn30; rotationAngle 270; }
+nonConformalCoupleIn300 { $nonConformalCoupleIn30; rotationAngle 300; }
+nonConformalCoupleIn330 { $nonConformalCoupleIn30; rotationAngle 330; }
+
+nonConformalCoupleOut0
+{
+    patches         (nonCoupleOut1 nonCoupleOut2);
+    transform       none;
+}
+
+nonConformalCoupleOut30
+{
+    patches         (nonCoupleOut1 nonCoupleOut2);
+    transform       rotational;
+    rotationAxis    (-1 0 0);
+    rotationCentre  (0 0 0);
+    rotationAngle   -30;
+}
+
+nonConformalCoupleOut60 { $nonConformalCoupleOut30; rotationAngle -60; }
+nonConformalCoupleOut90 { $nonConformalCoupleOut30; rotationAngle -90; }
+nonConformalCoupleOut120 { $nonConformalCoupleOut30; rotationAngle -120; }
+nonConformalCoupleOut150 { $nonConformalCoupleOut30; rotationAngle -150; }
+nonConformalCoupleOut180 { $nonConformalCoupleOut30; rotationAngle -180; }
+nonConformalCoupleOut210 { $nonConformalCoupleOut30; rotationAngle -210; }
+nonConformalCoupleOut240 { $nonConformalCoupleOut30; rotationAngle -240; }
+nonConformalCoupleOut270 { $nonConformalCoupleOut30; rotationAngle -270; }
+nonConformalCoupleOut300 { $nonConformalCoupleOut30; rotationAngle -300; }
+nonConformalCoupleOut330 { $nonConformalCoupleOut30; rotationAngle -330; }
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/impeller/system/fvSolution

@@ -54,13 +54,11 @@ solvers
 
 PIMPLE
 {
-    correctPhi          yes;
+    correctPhi          no;
+    correctMeshPhi      no;
     nOuterCorrectors    2;
     nCorrectors         1;
     nNonOrthogonalCorrectors 0;
-
-    pRefCell            0;
-    pRefValue           0;
 }
 
 relaxationFactors

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/0/U

@@ -20,43 +20,30 @@ internalField   uniform (0 0 0);
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     inlet
     {
         type            fixedValue;
         value           uniform (1 0 0);
     }
+
     outlet
     {
         type            pressureInletOutletVelocity;
-        value           uniform (0 0 0);
+        value           $internalField;
     }
+
     walls
     {
-//        type            fixedValue;
         type            movingWallVelocity;
-        value           uniform (0 0 0);
+        value           $internalField;
     }
-    defaultFaces
+
+    "nonCouple.*"
     {
-        type            empty;
+        type            movingWallVelocity;
-    }
+        value           $internalField;
-    ACMI1_blockage
-    {
-        type            noSlip;
-    }
-    ACMI1_couple
-    {
-        type            cyclicACMI;
-        value           uniform (0 0 0);
-    }
-    ACMI2_blockage
-    {
-        type            noSlip;
-    }
-    ACMI2_couple
-    {
-        type            cyclicACMI;
-        value           uniform (0 0 0);
     }
 }
 

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/0/epsilon

@@ -16,43 +16,36 @@ FoamFile
 
 dimensions      [0 2 -3 0 0 0 0];
 
-internalField   uniform 20;
+internalField   uniform 1.8e-3;
 
 boundaryField
 {
-    rotor
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+
+    walls
     {
         type            epsilonWallFunction;
         value           $internalField;
     }
 
-    stator
+    "nonCouple.*"
     {
         type            epsilonWallFunction;
         value           $internalField;
     }
-
-    front
-    {
-        type            empty;
-    }
-
-    back
-    {
-        type            empty;
-    }
-
-    AMI1
-    {
-        type            cyclicAMI;
-        value           $internalField;
-    }
-
-    AMI2
-    {
-        type            cyclicAMI;
-        value           $internalField;
-    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/0/k

@@ -20,46 +20,32 @@ internalField   uniform 3.75e-3;
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     inlet
     {
         type            fixedValue;
         value           $internalField;
     }
+
     outlet
     {
         type            inletOutlet;
         inletValue      $internalField;
         value           $internalField;
     }
+
     walls
     {
         type            kqRWallFunction;
         value           $internalField;
     }
-    defaultFaces
+
-    {
+    "nonCouple.*"
-        type            empty;
-    }
-    ACMI1_blockage
     {
         type            kqRWallFunction;
         value           $internalField;
     }
-    ACMI1_couple
-    {
-        type            cyclicACMI;
-        value           $internalField;
-    }
-    ACMI2_blockage
-    {
-        type            kqRWallFunction;
-        value           $internalField;
-    }
-    ACMI2_couple
-    {
-        type            cyclicACMI;
-        value           $internalField;
-    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/0/nut

@@ -20,37 +20,30 @@ internalField   uniform 0;
 
 boundaryField
 {
-    rotor
+    #includeEtc "caseDicts/setConstraintTypes"
-    {
-        type            nutkWallFunction;
-        value           uniform 0;
-    }
 
-    stator
+    inlet
     {
-        type            nutkWallFunction;
+        type            fixedValue;
-        value           uniform 0;
-    }
-
-    front
-    {
-        type            empty;
-    }
-
-    back
-    {
-        type            empty;
-    }
-
-    AMI1
-    {
-        type            cyclicAMI;
         value           $internalField;
     }
 
-    AMI2
+    outlet
     {
-        type            cyclicAMI;
+        type            inletOutlet;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+
+    walls
+    {
+        type            nutkWallFunction;
+        value           $internalField;
+    }
+
+    "nonCouple.*"
+    {
+        type            nutkWallFunction;
         value           $internalField;
     }
 }

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/0/omega

@@ -0,0 +1,52 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      omega;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 -1 0 0 0 0];
+
+internalField   uniform 5.33333;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+
+    outlet
+    {
+        type            inletOutlet;
+        inletValue      $internalField;
+        value           $internalField;
+    }
+
+    walls
+    {
+        type            omegaWallFunction;
+        value           $internalField;
+    }
+
+    "nonCouple.*"
+    {
+        type            omegaWallFunction;
+        value           $internalField;
+    }
+}
+
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/0/p

@@ -9,6 +9,7 @@ FoamFile
 {
     format      ascii;
     class       volScalarField;
+    location    "0";
     object      p;
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -19,37 +20,29 @@ internalField   uniform 0;
 
 boundaryField
 {
-    rotor
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
     {
         type            zeroGradient;
     }
 
-    stator
+    outlet
+    {
+        type            fixedValue;
+        value           uniform 0;
+    }
+
+    walls
     {
         type            zeroGradient;
     }
 
-    front
+    "nonCouple.*"
     {
-        type            empty;
+        type            zeroGradient;
+    }
 }
 
-    back
-    {
-        type            empty;
-    }
-
-    AMI1
-    {
-        type            cyclicAMI;
-        value           $internalField;
-    }
-
-    AMI2
-    {
-        type            cyclicAMI;
-        value           $internalField;
-    }
-}
 
 // ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/Allrun

@@ -0,0 +1,11 @@
+#!/bin/sh
+
+cd ${0%/*} || exit 1
+
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+runApplication blockMesh
+
+runApplication createNonConformalCouples -overwrite nonCouple1 nonCouple2
+
+runApplication $(getApplication)

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/constant/momentumTransport

@@ -18,7 +18,7 @@ simulationType  RAS;
 
 RAS
 {
-    model           kEpsilon;
+    model           realizableKE;
 
     turbulence      on;
 

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/system/blockMeshDict

@@ -38,14 +38,13 @@ vertices
 
 blocks
 (
-    hex (0 1 2 3 4 5 6 7) (80 40 1) simpleGrading (1 1 1)
+    hex (0 1 2 3 4 5 6 7) inletChannel (80 40 1) simpleGrading (1 1 1)
     hex (8 9 10 11 12 13 14 15) (80 96 1) simpleGrading (1 1 1)
 );
 
-defaultPatch
+edges
-{
+(
-    type empty;
+);
-}
 
 boundary
 (
@@ -76,17 +75,17 @@ boundary
             (9 13 12 8)
         );
     }
-    couple1
+    nonCouple1
     {
-        type patch;
+        type wall;
         faces
         (
             (2 6 5 1)
         );
     }
-    couple2
+    nonCouple2
     {
-        type patch;
+        type wall;
         faces
         (
             (8 12 15 11)
@@ -94,5 +93,8 @@ boundary
     }
 );
 
+mergePatchPairs
+(
+);
 
 // ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/system/controlDict

@@ -48,5 +48,11 @@ adjustTimeStep  true;
 
 maxCo           0.5;
 
+functions
+{
+    #includeFunc patchFlowRate(patch=inlet, funcName=inletFlowRate)
+    #includeFunc patchFlowRate(patch=outlet, funcName=outletFlowRate)
+}
+
 
 // ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/system/fvSchemes

@@ -29,16 +29,16 @@ gradSchemes
 divSchemes
 {
     default         none;
-//    div(phi,U)      Gauss upwind;
     div(phi,U)      Gauss linearUpwind grad(U);
     div(phi,k)      Gauss upwind;
     div(phi,epsilon) Gauss upwind;
+    div(phi,omega)  Gauss upwind;
     div((nuEff*dev2(T(grad(U))))) Gauss linear;
 }
 
 laplacianSchemes
 {
-    default         Gauss linear limited corrected 0.33;
+    default         Gauss linear corrected;
 }
 
 interpolationSchemes
@@ -48,7 +48,7 @@ interpolationSchemes
 
 snGradSchemes
 {
-    default         limited corrected 0.33;
+    default         corrected;
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/oscillatingInlet/system/fvSolution

@@ -18,39 +18,35 @@ solvers
     "pcorr.*"
     {
         solver          GAMG;
+        smoother        GaussSeidel;
         tolerance       1e-2;
         relTol          0;
-        smoother        GaussSeidel;
-        cacheAgglomeration no;
-        maxIter         50;
     }
 
     p
     {
         $pcorr;
-        tolerance       1e-5;
+        tolerance       1e-6;
         relTol          0.01;
     }
 
     pFinal
     {
         $p;
-        tolerance       1e-6;
         relTol          0;
     }
 
-    "(U|k|epsilon)"
+    "(U|k|epsilon|omega)"
     {
-        solver          smoothSolver;
+        solver          PBiCGStab;
-        smoother        symGaussSeidel;
+        preconditioner  DILU;
         tolerance       1e-6;
         relTol          0.1;
     }
 
-    "(U|k|epsilon)Final"
+    "(U|k|epsilon|omega)Final"
     {
         $U;
-        tolerance       1e-6;
         relTol          0;
     }
 }
@@ -58,15 +54,14 @@ solvers
 PIMPLE
 {
     correctPhi          no;
+    correctMeshPhi      no;
     nOuterCorrectors    1;
     nCorrectors         2;
     nNonOrthogonalCorrectors 0;
 }
 
 relaxationFactors
-{
+{}
-//    "(U|k|epsilon).*"   1;
-}
 
 cache
 {

tutorials/incompressible/pimpleFoam/RAS/oscillatingInletACMI2D/0/p

@@ -1,69 +0,0 @@
-/*--------------------------------*- C++ -*----------------------------------*\
-  =========                 |
-  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
-   \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Version:  dev
-     \\/     M anipulation  |
-\*---------------------------------------------------------------------------*/
-FoamFile
-{
-    format      ascii;
-    class       volScalarField;
-    location    "0";
-    object      p;
-}
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-dimensions      [0 2 -2 0 0 0 0];
-
-internalField   uniform 0;
-
-boundaryField
-{
-    inlet
-    {
-        type            zeroGradient;
-    }
-    outlet
-    {
-        type            fixedValue;
-        value           uniform 0;
-    }
-    walls
-    {
-        type            zeroGradient;
-    }
-    couple1
-    {
-        type            zeroGradient;
-    }
-    couple2
-    {
-        type            zeroGradient;
-    }
-    defaultFaces
-    {
-        type            empty;
-    }
-    ACMI1_blockage
-    {
-        type            zeroGradient;
-    }
-    ACMI1_couple
-    {
-        type            cyclicACMI;
-        value           uniform 0;
-    }
-    ACMI2_blockage
-    {
-        type            zeroGradient;
-    }
-    ACMI2_couple
-    {
-        type            cyclicACMI;
-        value           uniform 0;
-    }
-}
-
-
-// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/oscillatingInletACMI2D/Allmesh

@@ -1,14 +0,0 @@
-#!/bin/sh
-cd ${0%/*} || exit 1    # Run from this directory
-
-# Source tutorial run functions
-. $WM_PROJECT_DIR/bin/tools/RunFunctions
-
-runApplication blockMesh
-
-runApplication topoSet -constant
-
-# Split the mesh to generate the ACMI coupled patches
-runApplication createBaffles -overwrite
-
-#------------------------------------------------------------------------------

tutorials/incompressible/pimpleFoam/RAS/oscillatingInletACMI2D/Allrun

@@ -1,9 +0,0 @@
-#!/bin/sh
-cd ${0%/*} || exit 1    # Run from this directory
-
-# Source tutorial run functions
-. $WM_PROJECT_DIR/bin/tools/RunFunctions
-
-./Allmesh
-
-runApplication $(getApplication)

tutorials/incompressible/pimpleFoam/RAS/oscillatingInletACMI2D/Allrun-parallel

@@ -1,13 +0,0 @@
-#!/bin/sh
-cd ${0%/*} || exit 1    # Run from this directory
-
-# Source tutorial run functions
-. $WM_PROJECT_DIR/bin/tools/RunFunctions
-
-./Allmesh
-
-runApplication decomposePar
-
-runParallel $(getApplication)
-
-runApplication reconstructPar

tutorials/incompressible/pimpleFoam/RAS/oscillatingInletACMI2D/README

@@ -1,131 +0,0 @@
-oscillatingInletACMI2D
-
-This tutorial case gives an example of the Arbitrarily Coupled Mesh Interface
-(ACMI) usage.  The mesh is composed of two mesh regions: an inlet channel which
-oscillates in the +/- Y-direction, and a fixed mesh region.
-
-Each ACMI patch requires the specification of a 'non-overlapping' patch.  In
-this example, the non-overlapping patches are described as walls, e.g. taken
-from the constant/polyMesh/boundary file:
-
-    1. First ACMI patch pair applied to the inlet channel outlet
-
-        ACMI1_couple
-        {
-            type            cyclicACMI;
-            nFaces          40;
-            startFace       43720;
-            matchTolerance  0.0001;
-            neighbourPatch  ACMI2_couple;
-            nonOverlapPatch ACMI1_blockage;
-        }
-        ACMI1_blockage
-        {
-            type            wall;
-            nFaces          40;
-            startFace       43680;
-        }
-
-
-    1. Second ACMI patch pair applied to the fixed mesh region inlet
-
-        ACMI2_couple
-        {
-            type            cyclicACMI;
-            nFaces          96;
-            startFace       43856;
-            matchTolerance  0.0001;
-            neighbourPatch  ACMI1_couple;
-            nonOverlapPatch ACMI2_blockage;
-        }
-        ACMI2_blockage
-        {
-            type            wall;
-            nFaces          96;
-            startFace       43760;
-        }
-
-
-In the above, the ACMI1_blockage and ACMI1_couple patches occupy the same space,
-with duplicate points, edges and faces.  The ACMI2_blockage and ACMI2_couple
-patches are created similarly.
-
-The duplicate patches are initially created using the createBaffles utility.
-Firstly, the original (non-duplicated) patch faces are collected into zones
-using the topoSet utility.
-
-Note: when specifying the ACMI patch-pairs, the coupled patch must be specified
-before its associated non-overlapping patch.
-
-Each ACMI/non-overlapping patch pair is specified using a master-slave approach.
-However, since we are generating boundary patches (which are always master
-patches) the slave patches are simply defined using 'dummy' entries, e.g.:
-
-    type        faceZone;
-    zoneName    couple1Faces;
-
-    patches
-    {
-        // create cyclic ACMI patch
-        master
-        {
-            //- Master side patch
-            name            ACMI1_couple;
-            type            cyclicACMI;
-            matchTolerance  0.0001;
-            neighbourPatch  ACMI2_couple;
-            nonOverlapPatch ACMI1_blockage;
-        }
-        slave // dummy entries only
-        {
-            //- Slave side patch
-            name            ACMI1_couple;
-            type            patch;
-        }
-        // create blockage patch
-        master2
-        {
-            //- Master side patch
-            name            ACMI1_blockage;
-            type            wall;
-        }
-        slave12 // dummy entries only
-        {
-            //- Slave side patch
-            name            ACMI1_blockage;
-            type            wall;
-        }
-
-    }
-
-Boundary conditions must then be applied to all geometric patches in the usual,
-manner, and the cases can be executed in parallel (as shown when running the
-Allrun-parallel script) without any special treatment, i.e. the case set-up is
-the same as when operating in serial mode.
-
-
-checkMesh
----------
-checkMesh will see the 'duplicate' boundary faces but does not know about
-the area factors so will complain:
-
-     ***Boundary openness (-0.0103092 2.3845e-17 3.80774e-17) possible hole in boundary description.
-     ***Open cells found, max cell openness: 0.333333, number of open cells 136
-      <<Writing 136 non closed cells to set nonClosedCells
-
-As long as these non-closed cells are on the ACMI they can be ignored.
-
-
-paraFoam
---------
-- display: it will not display the outside of the ACMI, there will be a hole.
-- cuttingPlanes: they leave out the duplicate faces (i.e. the cells on
-  the ACMI patches)
-- interpolation: the interpolation does not take into account the
-area-weights on the ACMI.
-
-
-pointFields
------------
-Same as paraFoam: the interpolation does not take into account the
-area-weights on the ACMI.

tutorials/incompressible/pimpleFoam/RAS/oscillatingInletACMI2D/system/createBafflesDict

@@ -1,106 +0,0 @@
-/*--------------------------------*- C++ -*----------------------------------*\
-  =========                 |
-  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
-   \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Version:  dev
-     \\/     M anipulation  |
-\*---------------------------------------------------------------------------*/
-FoamFile
-{
-    format      ascii;
-    class       dictionary;
-    object      createBafflesDict;
-}
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-// Whether to convert internal faces only (so leave boundary faces intact).
-// This is only relevant if your face selection type can pick up boundary
-// faces.
-internalFacesOnly false;
-
-
-// Baffles to create.
-baffles
-{
-    // NOTE: cyclicAMI patches MUST BE defined PRIOR to their associted
-    //       blockage patches
-
-    ACMI1
-    {
-        //- Use predefined faceZone to select faces and orientation.
-        type        faceZone;
-        zoneName    couple1Faces;
-
-        patches
-        {
-            master
-            {
-                //- Master side patch
-                name            ACMI1_couple;
-                type            cyclicACMI;
-                neighbourPatch  ACMI2_couple;
-                nonOverlapPatch ACMI1_blockage;
-            }
-            slave // not used since we're manipulating a boundary patch
-            {
-                //- Slave side patch
-                name            ACMI1_couple;
-                type            patch;
-            }
-
-            master2
-            {
-                //- Master side patch
-                name            ACMI1_blockage;
-                type            wall;
-            }
-            slave2 // not used since we're manipulating a boundary patch
-            {
-                //- Slave side patch
-                name            ACMI1_blockage;
-                type            wall;
-            }
-
-        }
-    }
-    ACMI2
-    {
-        //- Use predefined faceZone to select faces and orientation.
-        type        faceZone;
-        zoneName    couple2Faces;
-
-        patches
-        {
-            master
-            {
-                //- Master side patch
-                name            ACMI2_couple;
-                type            cyclicACMI;
-                neighbourPatch  ACMI1_couple;
-                nonOverlapPatch ACMI2_blockage;
-            }
-            slave // not used since we're manipulating a boundary patch
-            {
-                //- Slave side patch
-                name            ACMI2_couple;
-                type            patch;
-            }
-
-            master2
-            {
-                //- Master side patch
-                name            ACMI2_blockage;
-                type            wall;
-            }
-            slave2 // not used since we're manipulating a boundary patch
-            {
-                //- Slave side patch
-                name            ACMI2_blockage;
-                type            wall;
-            }
-        }
-    }
-}
-
-
-// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/oscillatingInletACMI2D/system/topoSetDict

@@ -1,70 +0,0 @@
-/*--------------------------------*- C++ -*----------------------------------*\
-  =========                 |
-  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
-   \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Version:  dev
-     \\/     M anipulation  |
-\*---------------------------------------------------------------------------*/
-FoamFile
-{
-    format      ascii;
-    class       dictionary;
-    object      topoSetDict;
-}
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-actions
-(
-    // Get both sides of ami
-    // ~~~~~~~~~~~~~~~~~~~~~
-
-    // Create faceZone for patch couple1
-    {
-        name    couple1FaceSet;
-        type    faceSet;
-        action  new;
-        source  patchToFace;
-        patch   couple1;
-    }
-    {
-        name    couple1Faces;
-        type    faceZoneSet;
-        action  new;
-        source  setToFaceZone;
-        faceSet couple1FaceSet;
-    }
-
-    // Create faceZone for patch couple2
-    {
-        name    couple2FaceSet;
-        type    faceSet;
-        action  new;
-        source  patchToFace;
-        patch   couple2;
-    }
-    {
-        name    couple2Faces;
-        type    faceZoneSet;
-        action  new;
-        source  setToFaceZone;
-        faceSet couple2FaceSet;
-    }
-
-    // Create cellZone for moving cells in inlet channel
-    {
-        name    inletChannelCellSet;
-        type    cellSet;
-        action  new;
-        source  boxToCell;
-        box     (-100 -100 -100) (1.0001 100 100);
-    }
-    {
-        name    inletChannel;
-        type    cellZoneSet;
-        action  new;
-        source  setToCellZone;
-        set     inletChannelCellSet;
-    }
-);
-
-// ************************************************************************* //

tutorials/incompressible/pimpleFoam/RAS/propeller/0/U

@@ -20,7 +20,6 @@ internalField   uniform (0 0 0);
 
 boundaryField
 {
-    //- Set patchGroups for constraint patches
     #includeEtc "caseDicts/setConstraintTypes"
 
     inlet
@@ -32,8 +31,8 @@ boundaryField
     outlet
     {
         type            inletOutlet;
-        inletValue      uniform (0 0 0);
+        inletValue      $internalField;
-        value           uniform (0 0 0);
+        value           $internalField;
     }
 
     outerCylinder
@@ -44,7 +43,13 @@ boundaryField
     "propeller.*"
     {
         type            movingWallVelocity;
-        value           uniform (0 0 0);
+        value           $internalField;
+    }
+
+    nonCouple
+    {
+        type            movingWallSlipVelocity;
+        value           $internalField;
     }
 }
 

tutorials/incompressible/pimpleFoam/RAS/propeller/0/epsilon

@@ -20,7 +20,6 @@ internalField   uniform 0.0495;
 
 boundaryField
 {
-    //- Set patchGroups for constraint patches
     #includeEtc "caseDicts/setConstraintTypes"
 
     inlet
@@ -41,6 +40,11 @@ boundaryField
         type            epsilonWallFunction;
         value           $internalField;
     }
+
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/propeller/0/k

@@ -20,7 +20,6 @@ internalField   uniform 0.06;
 
 boundaryField
 {
-    //- Set patchGroups for constraint patches
     #includeEtc "caseDicts/setConstraintTypes"
 
     inlet
@@ -41,6 +40,11 @@ boundaryField
         type            kqRWallFunction;
         value           $internalField;
     }
+
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/propeller/0/nut

@@ -20,25 +20,30 @@ internalField   uniform 0;
 
 boundaryField
 {
-    //- Set patchGroups for constraint patches
     #includeEtc "caseDicts/setConstraintTypes"
 
     inlet
     {
         type            calculated;
-        value           uniform 0;
+        value           $internalField;
     }
 
     outlet
     {
         type            calculated;
-        value           uniform 0;
+        value           $internalField;
     }
 
     wall
     {
         type            nutkWallFunction;
-        value           uniform 0;
+        value           $internalField;
+    }
+
+    nonCouple
+    {
+        type            calculated;
+        value           $internalField;
     }
 }
 

tutorials/incompressible/pimpleFoam/RAS/propeller/0/p

@@ -20,7 +20,6 @@ internalField   uniform 0;
 
 boundaryField
 {
-    //- Set patchGroups for constraint patches
     #includeEtc "caseDicts/setConstraintTypes"
 
     inlet
@@ -31,13 +30,18 @@ boundaryField
     outlet
     {
         type            fixedValue;
-        value           uniform 0;
+        value           $internalField;
     }
 
     wall
     {
         type            zeroGradient;
     }
+
+    nonCouple
+    {
+        type            zeroGradient;
+    }
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/propeller/Allclean

@@ -7,6 +7,9 @@ cd ${0%/*} || exit 1    # Run from this directory
 # Remove surfaces
 rm -f constant/geometry/propeller*.obj.gz
 
+# Remove post-processing fields
+rm -f 0/Q
+
 cleanCase
 cleanExplicitFeatures
 

tutorials/incompressible/pimpleFoam/RAS/propeller/Allmesh

@@ -7,16 +7,20 @@ cd ${0%/*} || exit 1    # Run from this directory
 # Copy propeller surfaces from resources directory
 cp $FOAM_TUTORIALS/resources/geometry/propeller*.obj.gz constant/geometry/
 
-# Meshing
+# Mesh in parallel
 runApplication blockMesh
 runApplication surfaceFeatures
-runApplication snappyHexMesh -overwrite
+runApplication decomposePar -noFields
+runParallel snappyHexMesh -overwrite
 
-# Create the inlet/outlet and AMI patches
+# Create the sliding patches
-runApplication createBaffles -overwrite
+runParallel createBaffles -overwrite
-runApplication splitBaffles -overwrite
+runParallel splitBaffles -overwrite
 
-# Renumbering
+# Renumber
-runApplication renumberMesh -noFields -overwrite
+runParallel renumberMesh -noFields -overwrite
+
+# Create non-conformal couples
+runParallel createNonConformalCouples -overwrite nonCouple1 nonCouple2
 
 #------------------------------------------------------------------------------

tutorials/incompressible/pimpleFoam/RAS/propeller/Allrun

@@ -8,10 +8,18 @@ cd ${0%/*} || exit 1    # Run from this directory
 if ! cloneMesh ../../../../multiphase/interFoam/RAS/propeller .
 then
     ./Allmesh
+else
+    runApplication decomposePar -noFields
 fi
 
-runApplication decomposePar
+# Copy the fields into the processor directories
+runApplication -a decomposePar -fields -copyZero
+
+# Run
 runParallel $(getApplication)
+
+# Reconstruct
+runApplication reconstructParMesh -constant
 runApplication reconstructPar
 
 #------------------------------------------------------------------------------

tutorials/incompressible/pimpleFoam/RAS/propeller/system/createBafflesDict

@@ -9,7 +9,7 @@ FoamFile
 {
     format      ascii;
     class       dictionary;
-    object      createPatchDict;
+    object      createBafflesDict;
 }
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -27,16 +27,16 @@ baffles
         {
             master
             {
-                name            AMI1;
+                name            nonCouple1;
-                type            cyclicAMI;
+                type            patch;
-                neighbourPatch  AMI2;
+                inGroups        (nonCouple);
             }
 
             slave
             {
-                name            AMI2;
+                name            nonCouple2;
-                type            cyclicAMI;
+                type            patch;
-                neighbourPatch  AMI1;
+                inGroups        (nonCouple);
             }
         }
     }

tutorials/incompressible/pimpleFoam/RAS/propeller/system/decomposeParDict

@@ -14,13 +14,13 @@ FoamFile
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-numberOfSubdomains 4;
+numberOfSubdomains 8;
 
 method          hierarchical;
 
 hierarchicalCoeffs
 {
-    n               (1 4 1);
+    n               (2 2 2);
     order           xyz;
 }
 

tutorials/incompressible/pimpleFoam/RAS/propeller/system/fvSchemes

@@ -22,14 +22,12 @@ ddtSchemes
 gradSchemes
 {
     default         Gauss linear;
-    grad(p)         Gauss linear;
     grad(U)         cellLimited Gauss linear 1;
 }
 
 divSchemes
 {
     default         none;
-//    div(phi,U)      Gauss upwind;
     div(phi,U)      Gauss linearUpwind grad(U);
     div(phi,k)      Gauss upwind;
     div(phi,epsilon) Gauss upwind;
@@ -38,7 +36,7 @@ divSchemes
 
 laplacianSchemes
 {
-    default         Gauss linear limited corrected 0.33;
+    default         Gauss linear corrected;
 }
 
 interpolationSchemes
@@ -48,7 +46,7 @@ interpolationSchemes
 
 snGradSchemes
 {
-    default         limited corrected 0.33;
+    default         corrected;
 }
 
 

tutorials/incompressible/pimpleFoam/RAS/propeller/system/fvSolution

@@ -19,45 +19,44 @@ solvers
     {
         solver          GAMG;
         smoother        DICGaussSeidel;
+        tolerance       1e-6;
+        relTol          0.01;
         cacheAgglomeration no;
         maxIter         50;
-
-        tolerance       1e-5;
-        relTol          0.01;
     }
 
     pFinal
     {
         $p;
-
-        tolerance       1e-6;
         relTol          0;
     }
 
     "pcorr.*"
     {
         $p;
-
         tolerance       1e-2;
         relTol          0;
+    }
 
+    MeshPhi
+    {
+        solver          smoothSolver;
+        smoother        symGaussSeidel;
+        tolerance       1e-2;
+        relTol          0;
     }
 
     "(U|k|epsilon)"
     {
         solver          smoothSolver;
         smoother        symGaussSeidel;
-
         tolerance       1e-6;
         relTol          0.1;
     }
 
     "(U|k|epsilon)Final"
     {
-        solver          smoothSolver;
+        $U;
-        smoother        symGaussSeidel;
-
-        tolerance       1e-6;
         relTol          0;
     }
 }
@@ -65,6 +64,7 @@ solvers
 PIMPLE
 {
     correctPhi          yes;
+    correctMeshPhi      yes;
     nOuterCorrectors    2;
     nCorrectors         1;
     nNonOrthogonalCorrectors 0;

tutorials/incompressible/pimpleFoam/laminar/mixerVesselAMI2D/Allrun

@@ -1,19 +0,0 @@
-#!/bin/sh
-cd ${0%/*} || exit 1    # Run from this directory
-
-# Source tutorial run functions
-. $WM_PROJECT_DIR/bin/tools/RunFunctions
-
-application=$(getApplication)
-
-runApplication blockMesh -dict $FOAM_TUTORIALS/resources/blockMesh/mixerVessel2D
-runApplication createBaffles -overwrite
-runApplication splitBaffles -overwrite
-runApplication topoSet
-
-#runApplication $application
-runApplication decomposePar
-runParallel $application
-runApplication reconstructPar
-
-#------------------------------------------------------------------------------

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/0/U

@@ -0,0 +1,47 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volVectorField;
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 1 -1 0 0 0 0];
+
+internalField   uniform (1 0 0);
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+
+    outlet
+    {
+        type            pressureInletOutletVelocity;
+        value           $internalField;
+    }
+
+    walls
+    {
+        type            noSlip;
+    }
+
+    nonCouple
+    {
+        type            slip;
+    }
+}
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/0/p

@@ -0,0 +1,47 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 2 -2 0 0 0 0];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    #includeEtc "caseDicts/setConstraintTypes"
+
+    inlet
+    {
+        type            fixedFluxPressure;
+    }
+
+    outlet
+    {
+        type            entrainmentPressure;
+        p0              $internalField;
+        value           $internalField;
+    }
+
+    walls
+    {
+        type            fixedFluxPressure;
+    }
+
+    nonCouple
+    {
+        type            fixedFluxPressure;
+    }
+}
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/Allclean

@@ -0,0 +1,7 @@
+#!/bin/sh
+
+cd "${0%/*}" || exit 1
+
+. $WM_PROJECT_DIR/bin/tools/CleanFunctions
+
+cleanCase && rm -rf constant/fvMesh && rm -f error.eps

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/Allrun

@@ -0,0 +1,46 @@
+#!/bin/sh
+
+cd "${0%/*}" || exit 1
+
+. $WM_PROJECT_DIR/bin/tools/RunFunctions
+
+mesh()
+{
+    runApplication -s $1 blockMesh -dict system/blockMeshDict.$1
+}
+
+run()
+{
+    runApplication -s $1 $(getApplication)
+    mv postProcessing/inletFlowRate postProcessing/inletFlowRate.$1
+    mv postProcessing/outletFlowRate postProcessing/outletFlowRate.$1
+}
+
+if isTest $@
+then
+    foamDictionary -entry startFrom -set startTime system/controlDict
+fi
+
+mesh AMI
+run AMI
+
+mesh nonConformalCyclic
+runApplication -s nonConformalCyclic createNonConformalCouples -overwrite nonCoupleUp nonCoupleDown
+run nonConformalCyclic
+
+inletFlowRateAMI=postProcessing/inletFlowRate.AMI/0/surfaceFieldValue.dat
+outletFlowRateAMI=postProcessing/outletFlowRate.AMI/0/surfaceFieldValue.dat
+inletFlowRateNonConformalCyclic=postProcessing/inletFlowRate.nonConformalCyclic/0/surfaceFieldValue.dat
+outletFlowRateNonConformalCyclic=postProcessing/outletFlowRate.nonConformalCyclic/0/surfaceFieldValue.dat
+
+gnuplot << EOF
+set terminal postscript eps color enhanced
+set output "error.eps"
+set xlabel "Time (s)"
+set ylabel "Inlet-Outlet Conservation Error (%)"
+plot \
+    "< paste $inletFlowRateAMI $outletFlowRateAMI" \
+    us 1:(100*(\$4+\$2)/\$2) every ::1 w l t "AMI" , \
+    "< paste $inletFlowRateNonConformalCyclic $outletFlowRateNonConformalCyclic" \
+    us 1:(100*(\$4+\$2)/\$2) every ::1 w l t "nonConformalCyclic"
+EOF

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/constant/momentumTransport

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

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/constant/physicalProperties

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

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/system/blockMeshDict.AMI

@@ -0,0 +1,44 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    object      blockMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#include "blockMeshDict.internal"
+
+boundary
+(
+    #include "blockMeshDict.boundary"
+
+    coupleUp
+    {
+        type cyclicAMI;
+        neighbourPatch coupleDown;
+        faces
+        (
+            (7 8 24 23)
+            (8 9 25 24)
+        );
+    }
+    coupleDown
+    {
+        type cyclicAMI;
+        neighbourPatch coupleUp;
+        faces
+        (
+            (10 11 27 26)
+            (11 12 28 27)
+        );
+    }
+);
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/system/blockMeshDict.boundary

@@ -0,0 +1,34 @@
+inlet
+{
+    type patch;
+    faces
+    (
+        (0 1 17 16)
+        (1 2 18 17)
+    );
+}
+outlet
+{
+    type patch;
+    faces
+    (
+        (13 14 30 29)
+        (14 15 31 30)
+    );
+}
+walls
+{
+    type wall;
+    faces
+    (
+        (0 3 19 16)
+        (3 5 21 19)
+        (5 4 20 21)
+        (4 7 23 20)
+        (10 13 29 26)
+
+        (2 6 18 22)
+        (6 9 25 22)
+        (12 15 31 28)
+    );
+}

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/system/blockMeshDict.internal

@@ -0,0 +1,47 @@
+convertToMeters 0.001;
+
+vertices
+(
+    (-150 0 -1)            (-150 25 -1) (-150 50 -1)
+    ( -50 0 -1) (-50 0 -1) ( -50 25 -1) ( -50 50 -1)
+    (   0 0 -1)            (   0 25 -1) (   0 50 -1)
+
+    (   0 0 -1)            (   0 25 -1) (   0 50 -1)
+    ( 150 0 -1)            ( 150 25 -1) ( 150 50 -1)
+
+    (-150 0  1)            (-150 25  1) (-150 50  1)
+    ( -50 0  1) (-50 0  1) ( -50 25  1) ( -50 50  1)
+    (   0 0  1)            (   0 25  1) (   0 50  1)
+
+    (   0 0  1)            (   0 25  1) (   0 50  1)
+    ( 150 0  1)            ( 150 25  1) ( 150 50  1)
+);
+
+blocks
+(
+    hex (0 3 5 1 16 19 21 17) (12 5 1) simpleGrading (1 1 1)
+    hex (1 5 6 2 17 21 22 18) (12 5 1) simpleGrading (1 1 1)
+    hex (4 7 8 5 20 23 24 21) (6 5 1) simpleGrading (1 1 1)
+    hex (5 8 9 6 21 24 25 22) (6 5 1) simpleGrading (1 1 1)
+    hex (10 13 14 11 26 29 30 27) (18 6 1) simpleGrading (1 1 1)
+    hex (11 14 15 12 27 30 31 28) (18 6 1) simpleGrading (1 1 1)
+);
+
+edges
+(
+    arc 7 8 30 (0 0 1)
+    arc 23 24 30 (0 0 1)
+    arc 8 9 30 (0 0 -1)
+    arc 24 25 30 (0 0 -1)
+
+    arc 10 11 30 (0 0 1)
+    arc 26 27 30 (0 0 1)
+    arc 11 12 30 (0 0 -1)
+    arc 27 28 30 (0 0 -1)
+);
+
+defaultPatch
+{
+    name frontAndBack;
+    type empty;
+}

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/system/blockMeshDict.nonConformalCyclic

@@ -0,0 +1,44 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    object      blockMeshDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#include "blockMeshDict.internal"
+
+boundary
+(
+    #include "blockMeshDict.boundary"
+
+    nonCoupleUp
+    {
+        type patch;
+        inGroups (nonCouple);
+        faces
+        (
+            (7 8 24 23)
+            (8 9 25 24)
+        );
+    }
+    nonCoupleDown
+    {
+        type patch;
+        inGroups (nonCouple);
+        faces
+        (
+            (10 11 27 26)
+            (11 12 28 27)
+        );
+    }
+);
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/system/controlDict

@@ -0,0 +1,58 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      controlDict;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+application     pimpleFoam;
+
+startFrom       startTime;
+
+startTime       0;
+
+stopAt          endTime;
+
+endTime         1;
+
+deltaT          0.005;
+
+writeControl    adjustableRunTime;
+
+writeInterval   0.1;
+
+purgeWrite      0;
+
+writeFormat     ascii;
+
+writePrecision  12;
+
+writeCompression off;
+
+timeFormat      general;
+
+timePrecision   6;
+
+runTimeModifiable yes;
+
+adjustTimeStep  no;
+
+maxCo           5;
+
+functions
+{
+    #includeFunc patchFlowRate(patch=inlet, funcName=inletFlowRate)
+    #includeFunc patchFlowRate(patch=outlet, funcName=outletFlowRate)
+}
+
+
+// ************************************************************************* //

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/system/fvSchemes

@@ -22,8 +22,6 @@ ddtSchemes
 gradSchemes
 {
     default         Gauss linear;
-    grad(p)         Gauss linear;
-    grad(U)         Gauss linear;
 }
 
 divSchemes
@@ -41,7 +39,6 @@ laplacianSchemes
 interpolationSchemes
 {
     default         linear;
-    interpolate(HbyA) linear;
 }
 
 snGradSchemes

tutorials/incompressible/pimpleFoam/laminar/nonConformalChannel/system/fvSolution

@@ -0,0 +1,56 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  dev
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvSolution;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+solvers
+{
+    p
+    {
+        solver          GAMG;
+        smoother        DICGaussSeidel;
+        preconditioner  DIC;
+        tolerance       1e-8;
+        relTol          0.01;
+    }
+
+    pFinal
+    {
+        $p;
+        relTol          0;
+    }
+
+    U
+    {
+        solver          smoothSolver;
+        smoother        symGaussSeidel;
+        tolerance       1e-5;
+        relTol          0.1;
+    }
+
+    UFinal
+    {
+        $U;
+        relTol          0;
+    }
+}
+
+PIMPLE
+{
+    nNonOrthogonalCorrectors 0;
+    nCorrectors         2;
+}
+
+
+// ************************************************************************* //

tutorials/lagrangian/particleFoam/mixerVessel2D/0/U

@@ -19,37 +19,23 @@ internalField   uniform (0 0 0);
 
 boundaryField
 {
+    #includeEtc "caseDicts/setConstraintTypes"
+
     rotor
     {
         type            movingWallVelocity;
-        value           uniform (0 0 0);
+        value           $internalField;
     }
 
     stator
     {
         type            movingWallVelocity;
-        value           uniform (0 0 0);
-    }
-
-    front
-    {
-        type            empty;
-    }
-
-    back
-    {
-        type            empty;
-    }
-
-    AMI1
-    {
-        type            cyclicAMI;
         value           $internalField;
     }
 
-    AMI2
+    nonCouple
     {
-        type            cyclicAMI;
+        type            movingWallSlipVelocity;
         value           $internalField;
     }
 }

tutorials/lagrangian/particleFoam/mixerVessel2D/Allrun

@@ -7,5 +7,10 @@ cd ${0%/*} || exit 1 # Run from this directory
 runApplication blockMesh -dict $FOAM_TUTORIALS/resources/blockMesh/mixerVessel2D
 runApplication createBaffles -overwrite
 runApplication splitBaffles -overwrite
+runApplication createNonConformalCouples -overwrite nonCouple1 nonCouple2
 
-runApplication $(getApplication)
+runApplication decomposePar
+
+runParallel $(getApplication)
+
+runApplication reconstructPar