diff --git a/applications/solvers/incompressible/pimpleFoam/UEqn.H b/applications/solvers/incompressible/pimpleFoam/UEqn.H
index 3574c6553a57a632d982afa53c3be42de27a62d4..3655677979ca0e81a0dea53b8e019a667c304d02 100644
--- a/applications/solvers/incompressible/pimpleFoam/UEqn.H
+++ b/applications/solvers/incompressible/pimpleFoam/UEqn.H
@@ -9,7 +9,7 @@ tmp<fvVectorMatrix> UEqn
 
 UEqn().relax();
 
-volScalarField rUA = 1.0/UEqn().A();
+volScalarField rAU = 1.0/UEqn().A();
 
 if (momentumPredictor)
 {
@@ -17,6 +17,6 @@ if (momentumPredictor)
 }
 else
 {
-    U = rUA*(UEqn().H() - fvc::grad(p));
+    U = rAU*(UEqn().H() - fvc::grad(p));
     U.correctBoundaryConditions();
 }
diff --git a/applications/solvers/incompressible/pimpleFoam/pEqn.H b/applications/solvers/incompressible/pimpleFoam/pEqn.H
index b5dc23a58cf17caa05252a4470fe174f07cacc27..acd86ed4b088d80c437fd17d2cc8ed9c75cbf8df 100644
--- a/applications/solvers/incompressible/pimpleFoam/pEqn.H
+++ b/applications/solvers/incompressible/pimpleFoam/pEqn.H
@@ -1,4 +1,4 @@
-U = rUA*UEqn().H();
+U = rAU*UEqn().H();
 
 if (nCorr <= 1)
 {
@@ -6,7 +6,7 @@ if (nCorr <= 1)
 }
 
 phi = (fvc::interpolate(U) & mesh.Sf())
-    + fvc::ddtPhiCorr(rUA, U, phi);
+    + fvc::ddtPhiCorr(rAU, U, phi);
 
 adjustPhi(phi, U, p);
 
@@ -16,7 +16,7 @@ for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
     // Pressure corrector
     fvScalarMatrix pEqn
     (
-        fvm::laplacian(rUA, p) == fvc::div(phi)
+        fvm::laplacian(rAU, p) == fvc::div(phi)
     );
 
     pEqn.setReference(pRefCell, pRefValue);
@@ -47,5 +47,5 @@ for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
 // Explicitly relax pressure for momentum corrector
 p.relax();
 
-U -= rUA*fvc::grad(p);
+U -= rAU*fvc::grad(p);
 U.correctBoundaryConditions();
diff --git a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/createFields.H b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/createFields.H
index bc798ba988df54a0fe3518443c6382ec5375d818..16b3bd977d07b6979d2cb84e5f12100cc8910921 100644
--- a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/createFields.H
+++ b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/createFields.H
@@ -27,7 +27,7 @@
         mesh
     );
 
-#   include "createPhi.H"
+    #include "createPhi.H"
 
 
     label pRefCell = 0;
diff --git a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H
index 155f43a6328d8acaf6aac801960f909d31784d4b..74045cabb7dfc05bfa3cac28d510f9ffa6bf4260 100644
--- a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H
+++ b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H
@@ -46,6 +46,7 @@ for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
 
 #include "continuityErrs.H"
 
+// Explicitly relax pressure for momentum corrector
 p.relax();
 
 // Make the fluxes relative to the mesh motion
diff --git a/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C b/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C
index 8f29b9e9343d81bcb1a6815d63196fa22b461e7e..4187da2d937cb7f2e64273bbeb75ce1b0d5b32a4 100644
--- a/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C
+++ b/src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C
@@ -27,7 +27,10 @@ License
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
-const Foam::scalar Foam::polyMeshTetDecomposition::minTetQuality = 1e-12;
+// Note: the use of this tolerance is ad-hoc, there may be extreme
+// cases where the resulting tetrahedra still have particle tracking
+// problems.
+const Foam::scalar Foam::polyMeshTetDecomposition::minTetQuality = SMALL;
 
 
 // * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //
diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/controlDict b/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/controlDict
index 590a47c55a6124c258398c7e6ce82f659aed0aea..86119ee8dae1ae16ff705f1252d960d210bf6791 100644
--- a/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/controlDict
+++ b/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/controlDict
@@ -15,8 +15,6 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-libs            ("libcompressibleTurbulenceModel.so" "libcompressibleRASModels.so");
-
 application     chtMultiRegionFoam;
 
 startFrom       latestTime;
@@ -54,5 +52,11 @@ maxDi           10.0;
 
 adjustTimeStep  yes;
 
+libs
+(
+    "libOpenFOAM.so"
+    "libcompressibleTurbulenceModel.so"
+    "libcompressibleRASModels.so"
+);
 
 // ************************************************************************* //
diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/controlDict b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/controlDict
index ec9f3455bd051faa1da52f6088d977381ec873b1..1d46884f63f1620675110ebc74a814515ed3937d 100644
--- a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/controlDict
+++ b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/controlDict
@@ -15,8 +15,6 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-libs            ("libcompressibleTurbulenceModel.so" "libcompressibleRASModels.so");
-
 application     chtMultiRegionFoam;
 
 startFrom       latestTime;
@@ -53,5 +51,11 @@ maxDi          10.0;
 
 adjustTimeStep  yes;
 
+libs
+(
+    "libOpenFOAM.so"
+    "libcompressibleTurbulenceModel.so"
+    "libcompressibleRASModels.so"
+);
 
 // ************************************************************************* //
diff --git a/tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/controlDict b/tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/controlDict
index 32e60a1dd1db695e64f63501af2c5e5de42616cc..f6ba2fe5fb575bd406412d259556087b2309bccd 100644
--- a/tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/controlDict
+++ b/tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/controlDict
@@ -15,8 +15,6 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-libs            ("libcompressibleTurbulenceModel.so" "libcompressibleRASModels.so");
-
 application     chtMultiRegionSimpleFoam;
 
 startFrom       startTime;
@@ -47,4 +45,11 @@ timePrecision   6;
 
 runTimeModifiable true;
 
+libs
+(
+    "libOpenFOAM.so"
+    "libcompressibleTurbulenceModel.so"
+    "libcompressibleRASModels.so"
+);
+
 // ************************************************************************* //
diff --git a/tutorials/incompressible/boundaryFoam/boundaryWallFunctionsProfile/system/controlDict b/tutorials/incompressible/boundaryFoam/boundaryWallFunctionsProfile/system/controlDict
index 8df3f4af88c6eba5c58b19cb7ddb7ab0a56821d3..7aafcc5b170d6e68634d9357670c0244ed09f55f 100644
--- a/tutorials/incompressible/boundaryFoam/boundaryWallFunctionsProfile/system/controlDict
+++ b/tutorials/incompressible/boundaryFoam/boundaryWallFunctionsProfile/system/controlDict
@@ -17,13 +17,6 @@ FoamFile
 
 application     boundaryFoam;
 
-libs
-(
-    "libOpenFOAM.so"
-    "libinterpolationTables_16x.so"
-    "libtabulatedWallFunctionFvPatchFields_16x.so"
-);
-
 startFrom       startTime;
 
 startTime       0;
@@ -54,5 +47,11 @@ runTimeModifiable true;
 
 graphFormat     raw;
 
+libs
+(
+    "libOpenFOAM.so"
+    "libinterpolationTables_16x.so"
+    "libtabulatedWallFunctionFvPatchFields_16x.so"
+);
 
 // ************************************************************************* //
diff --git a/tutorials/incompressible/pimpleDyMFoam/movingCone/0/U b/tutorials/incompressible/pimpleDyMFoam/movingCone/0/U
index b3bc15c145b210bb2952e7ccdaadb14b7a9a89a0..ddeb5b36b82d69bdaacdaee86bae3560829c94d7 100644
--- a/tutorials/incompressible/pimpleDyMFoam/movingCone/0/U
+++ b/tutorials/incompressible/pimpleDyMFoam/movingCone/0/U
@@ -20,47 +20,47 @@ internalField   uniform (0 0 0);
 
 boundaryField
 {
-    movingWall      
+    movingWall
     {
         type            fixedValue;
         value           uniform (1 0 0);
     }
 
-    farFieldMoving  
+    farFieldMoving
     {
         type            fixedValue;
         value           uniform (0 0 0);
     }
 
-    fixedWall       
+    fixedWall
     {
         type            fixedValue;
         value           uniform (0 0 0);
     }
 
-    axis            
+    axis
     {
         type            symmetryPlane;
     }
 
-    left            
+    left
     {
         type            pressureInletOutletVelocity;
         value           uniform (0 0 0);
     }
 
-    farField        
+    farField
     {
         type            fixedValue;
         value           uniform (0 0 0);
     }
 
-    back            
+    back
     {
         type            wedge;
     }
 
-    front           
+    front
     {
         type            wedge;
     }
diff --git a/tutorials/incompressible/pimpleDyMFoam/movingCone/0/p b/tutorials/incompressible/pimpleDyMFoam/movingCone/0/p
index 112e30ef2f56ac5c29a2435e810cc31e839cf970..9f13ffc4a40a8341d1c6b632673707473c48cfe5 100644
--- a/tutorials/incompressible/pimpleDyMFoam/movingCone/0/p
+++ b/tutorials/incompressible/pimpleDyMFoam/movingCone/0/p
@@ -20,27 +20,27 @@ internalField   uniform 0;
 
 boundaryField
 {
-    movingWall      
+    movingWall
     {
         type            zeroGradient;
     }
 
-    farFieldMoving  
+    farFieldMoving
     {
         type            zeroGradient;
     }
 
-    fixedWall       
+    fixedWall
     {
         type            zeroGradient;
     }
 
-    axis            
+    axis
     {
         type            symmetryPlane;
     }
 
-    left            
+    left
     {
         type            totalPressure;
         p0              uniform 0;
@@ -52,17 +52,17 @@ boundaryField
         value           uniform 0;
     }
 
-    farField        
+    farField
     {
         type            zeroGradient;
     }
 
-    back            
+    back
     {
         type            wedge;
     }
 
-    front           
+    front
     {
         type            wedge;
     }
diff --git a/tutorials/incompressible/pimpleDyMFoam/movingCone/system/fvSolution b/tutorials/incompressible/pimpleDyMFoam/movingCone/system/fvSolution
index 6bc425e2f66d4d7349ba458f1e775c110537f9bb..f93592537c32e22a50581e80cfacfce27fa43228 100644
--- a/tutorials/incompressible/pimpleDyMFoam/movingCone/system/fvSolution
+++ b/tutorials/incompressible/pimpleDyMFoam/movingCone/system/fvSolution
@@ -21,7 +21,7 @@ solvers
     {
         solver          PCG;
         preconditioner  DIC;
-        tolerance       0.01;
+        tolerance       0.02;
         relTol          0;
     }
 
@@ -30,7 +30,7 @@ solvers
         solver          PCG;
         preconditioner  DIC;
         tolerance       1e-06;
-        relTol          0.05;
+        relTol          0.01;
     }
 
     pFinal
@@ -66,10 +66,16 @@ solvers
 
 PIMPLE
 {
+    correctPhi          yes;
     nOuterCorrectors    2;
     nCorrectors         1;
     nNonOrthogonalCorrectors 0;
 }
 
+relaxationFactors
+{
+    U        1;
+    UFinal   1;
+}
 
 // ************************************************************************* //
diff --git a/tutorials/incompressible/pimpleDyMFoam/wingMotion/wingMotion2D_pimpleDyMFoam/system/controlDict b/tutorials/incompressible/pimpleDyMFoam/wingMotion/wingMotion2D_pimpleDyMFoam/system/controlDict
index df1315c0430788819ec733bcd87a87d4970cd0b1..db20cc97149d6a051aa00fd3cbc7b601caf57ce2 100644
--- a/tutorials/incompressible/pimpleDyMFoam/wingMotion/wingMotion2D_pimpleDyMFoam/system/controlDict
+++ b/tutorials/incompressible/pimpleDyMFoam/wingMotion/wingMotion2D_pimpleDyMFoam/system/controlDict
@@ -49,14 +49,18 @@ adjustTimeStep  yes;
 
 maxCo           0.9;
 
-libs ("libforces.so");
+libs
+(
+    "libOpenFOAM.so"
+    "libforces.so"
+);
 
 functions
 {
     forces
     {
         type                forces;
-        functionObjectLibs  ( "libforces.so" );
+        functionObjectLibs  ("libforces.so");
         outputControl       timeStep;
         outputInterval      10;
         patches             (wing);
diff --git a/tutorials/incompressible/pimpleFoam/t-junction-with-fan/system/controlDict b/tutorials/incompressible/pimpleFoam/t-junction-with-fan/system/controlDict
index 700875183870f4b0fbe48b4571dd3a7146456071..9dc31d5561fb3c86c93e949d18a3dd5b39d41977 100644
--- a/tutorials/incompressible/pimpleFoam/t-junction-with-fan/system/controlDict
+++ b/tutorials/incompressible/pimpleFoam/t-junction-with-fan/system/controlDict
@@ -15,8 +15,6 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-libs            ("libOpenFOAM.so" "libincompressibleTurbulenceModel.so" "libincompressibleRASModels.so");
-
 application     pimpleFoam;
 
 startFrom       startTime;
@@ -51,6 +49,13 @@ adjustTimeStep  yes;
 
 maxCo           5;
 
+libs
+(
+    "libOpenFOAM.so"
+    "libincompressibleTurbulenceModel.so"
+    "libincompressibleRASModels.so"
+);
+
 functions
 {
     probes
diff --git a/tutorials/incompressible/pimpleFoam/t-junction/system/controlDict b/tutorials/incompressible/pimpleFoam/t-junction/system/controlDict
index 700875183870f4b0fbe48b4571dd3a7146456071..e5b21cbaea53cd956a2b12d85d813a2bbd2ff0b9 100644
--- a/tutorials/incompressible/pimpleFoam/t-junction/system/controlDict
+++ b/tutorials/incompressible/pimpleFoam/t-junction/system/controlDict
@@ -15,8 +15,6 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-libs            ("libOpenFOAM.so" "libincompressibleTurbulenceModel.so" "libincompressibleRASModels.so");
-
 application     pimpleFoam;
 
 startFrom       startTime;
@@ -51,6 +49,13 @@ adjustTimeStep  yes;
 
 maxCo           5;
 
+libs
+(
+    "libOpenFOAM.so"
+    "libinterpolationTables_16x.so"
+    "libtabulatedWallFunctionFvPatchFields_16x.so"
+);
+
 functions
 {
     probes
diff --git a/tutorials/incompressible/simpleFoam/motorBike/system/controlDict b/tutorials/incompressible/simpleFoam/motorBike/system/controlDict
index 7be247e0791833a6ab0b8c108cd0bac4fec1b4c9..c81f8fd5885ea5463ed2cb12e3032fe154880fb1 100644
--- a/tutorials/incompressible/simpleFoam/motorBike/system/controlDict
+++ b/tutorials/incompressible/simpleFoam/motorBike/system/controlDict
@@ -14,8 +14,6 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-libs            ("libincompressibleRASModels.so");
-
 application     simpleFoam;
 
 startFrom       startTime;
@@ -46,6 +44,11 @@ timePrecision   6;
 
 runTimeModifiable true;
 
+libs
+(
+    "libOpenFOAM.so"
+    "libincompressibleRASModels.so"
+);
 
 functions
 {
diff --git a/tutorials/lagrangian/porousExplicitSourceReactingParcelFoam/filter/system/controlDict b/tutorials/lagrangian/porousExplicitSourceReactingParcelFoam/filter/system/controlDict
index 8aa86473df34c01b00b51c8e12e0fe1952455525..6190c66a0722fedd86b6b711d2df393480c11013 100644
--- a/tutorials/lagrangian/porousExplicitSourceReactingParcelFoam/filter/system/controlDict
+++ b/tutorials/lagrangian/porousExplicitSourceReactingParcelFoam/filter/system/controlDict
@@ -15,8 +15,6 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-libs            ("libOpenFOAM.so" "libcompressibleTurbulenceModel.so" "libcompressibleRASModels.so" "libradiation.so");
-
 application     porousExplicitSourceReactingParcelFoam;
 
 startFrom       startTime;
@@ -53,5 +51,12 @@ maxCo           1.0;
 
 maxDeltaT       1;
 
+libs
+(
+    "libOpenFOAM.so"
+    "libcompressibleTurbulenceModel.so"
+    "libcompressibleRASModels.so"
+    "libradiation.so"
+);
 
 // ************************************************************************* //
diff --git a/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/system/controlDict b/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/system/controlDict
index f2ddb38b90087b21b555d31a543d51ac8e197158..caf18048c1960b2f8846ca1f7f592475fb49f319 100644
--- a/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/system/controlDict
+++ b/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/system/controlDict
@@ -15,15 +15,6 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-libs
-(
-    "libOpenFOAM.so"
-    "libcompressibleTurbulenceModel.so"
-    "libcompressibleRASModels.so"
-    "libSLGThermo.so"
-    "libsurfaceFilmModels.so"
-);
-
 application     reactingParcelFoam;
 
 startFrom       latestTime; // startTime;
@@ -60,5 +51,13 @@ maxCo           0.5;
 
 maxDeltaT       0.005;
 
+libs
+(
+    "libOpenFOAM.so"
+    "libcompressibleTurbulenceModel.so"
+    "libcompressibleRASModels.so"
+    "libSLGThermo.so"
+    "libsurfaceFilmModels.so"
+);
 
 // ************************************************************************* //
diff --git a/tutorials/mesh/moveDynamicMesh/simpleHarmonicMotion/system/controlDict b/tutorials/mesh/moveDynamicMesh/simpleHarmonicMotion/system/controlDict
index 2d71827bb43c766f572f270e302d476e6e2b55c2..9e9f696faa2ee580de662d150aaaa6ba6d101f50 100644
--- a/tutorials/mesh/moveDynamicMesh/simpleHarmonicMotion/system/controlDict
+++ b/tutorials/mesh/moveDynamicMesh/simpleHarmonicMotion/system/controlDict
@@ -53,16 +53,9 @@ maxDeltaT       0.025;
 
 libs
 (
+    "libOpenFOAM.so"
     "libincompressibleRASModels.so"
     "libforces.so"
 );
 
-
-// libs
-// (
-//     "libgenericPatchFields.so"
-//     "libincompressibleRASModels.so"
-//     "libforces.so"
-// );
-
 // ************************************************************************* //