diff --git a/applications/utilities/finiteArea/checkFaMesh/faMeshWriteVTK.H b/applications/utilities/finiteArea/checkFaMesh/faMeshWriteVTK.H
index c6edf85a3e..5a138055d1 100644
--- a/applications/utilities/finiteArea/checkFaMesh/faMeshWriteVTK.H
+++ b/applications/utilities/finiteArea/checkFaMesh/faMeshWriteVTK.H
@@ -23,7 +23,7 @@ Description
         // vtk::formatType::INLINE_ASCII,
         fileName
         (
-            aMesh.mesh().time().globalPath() / "finiteArea"
+            aMesh.time().globalPath() / "finiteArea"
         )
     );
 
@@ -66,7 +66,7 @@ Description
         // vtk::formatType::INLINE_ASCII,
         fileName
         (
-            aMesh.mesh().time().globalPath() / "finiteArea-edges"
+            aMesh.time().globalPath() / "finiteArea-edges"
         )
     );
 
@@ -102,7 +102,7 @@ Description
         // vtk::formatType::INLINE_ASCII,
         fileName
         (
-            aMesh.mesh().time().globalPath() / "finiteArea-edgesCentres"
+            aMesh.time().globalPath() / "finiteArea-edgesCentres"
         )
     );
 
diff --git a/applications/utilities/finiteArea/makeFaMesh/faMeshWriteVTK.H b/applications/utilities/finiteArea/makeFaMesh/faMeshWriteVTK.H
index 1df8b4d312..ed988e7413 100644
--- a/applications/utilities/finiteArea/makeFaMesh/faMeshWriteVTK.H
+++ b/applications/utilities/finiteArea/makeFaMesh/faMeshWriteVTK.H
@@ -23,7 +23,7 @@ Description
         // vtk::formatType::INLINE_ASCII,
         fileName
         (
-            aMesh.mesh().time().globalPath() / "finiteArea"
+            aMesh.time().globalPath() / "finiteArea"
         )
     );
 
diff --git a/src/dynamicFaMesh/interfaceTrackingFvMesh/interfaceTrackingFvMesh.C b/src/dynamicFaMesh/interfaceTrackingFvMesh/interfaceTrackingFvMesh.C
index 8642f515db..ce798d638b 100644
--- a/src/dynamicFaMesh/interfaceTrackingFvMesh/interfaceTrackingFvMesh.C
+++ b/src/dynamicFaMesh/interfaceTrackingFvMesh/interfaceTrackingFvMesh.C
@@ -213,8 +213,8 @@ void Foam::interfaceTrackingFvMesh::makeUs() const
         IOobject
         (
             "Us",
-            mesh().time().timeName(),
-            mesh(),
+            aMesh().time().timeName(),
+            aMesh().thisDb(),
             IOobject::NO_READ,
             IOobject::AUTO_WRITE
         ),
@@ -271,8 +271,8 @@ void Foam::interfaceTrackingFvMesh::makeFsNetPhi() const
         IOobject
         (
             "fsNetPhi",
-            mesh().time().timeName(),
-            mesh(),
+            aMesh().time().timeName(),
+            aMesh().thisDb(),
             IOobject::NO_READ,
             IOobject::AUTO_WRITE
         ),
@@ -462,8 +462,8 @@ void Foam::interfaceTrackingFvMesh::makePhis()
         IOobject
         (
             "phis",
-            mesh().time().timeName(),
-            mesh(),
+            aMesh().time().timeName(),
+            aMesh().thisDb(),
             IOobject::NO_READ,
             IOobject::NO_WRITE
         ),
@@ -494,7 +494,7 @@ void Foam::interfaceTrackingFvMesh::makeSurfactConc() const
                 mesh().time().startTime().value()
             ),
             // mesh().time().timeName(),
-            mesh(),
+            aMesh().thisDb(),
             IOobject::MUST_READ,
             IOobject::AUTO_WRITE
         ),
@@ -525,7 +525,7 @@ void Foam::interfaceTrackingFvMesh::makeBulkSurfactConc() const
                 mesh().time().startTime().value()
             ),
             // mesh().time().timeName(),
-            mesh(),
+            aMesh().thisDb(),
             IOobject::MUST_READ,
             IOobject::AUTO_WRITE
         ),
@@ -559,8 +559,8 @@ void Foam::interfaceTrackingFvMesh::makeSurfaceTension() const
         IOobject
         (
             "surfaceTension",
-            mesh().time().timeName(),
-            mesh(),
+            aMesh().time().timeName(),
+            aMesh().thisDb(),
             IOobject::NO_READ,
             IOobject::NO_WRITE
         ),
@@ -604,28 +604,17 @@ void Foam::interfaceTrackingFvMesh::makeContactAngle()
     IOobject contactAngleHeader
     (
         "contactAngle",
-        mesh().time().timeName(),
-        mesh(),
-        IOobject::MUST_READ
+        aMesh().time().timeName(),
+        aMesh().thisDb(),
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
     );
 
     if (contactAngleHeader.typeHeaderOk<areaScalarField>())
     {
         Info<< "Reading contact angle field" << endl;
 
-        contactAnglePtr_ =
-            new areaScalarField
-            (
-                IOobject
-                (
-                    "contactAngle",
-                    mesh().time().timeName(),
-                    mesh(),
-                    IOobject::MUST_READ,
-                    IOobject::AUTO_WRITE
-                ),
-                aMesh()
-            );
+        contactAnglePtr_ = new areaScalarField(contactAngleHeader, aMesh());
     }
 }
 
@@ -886,8 +875,8 @@ void Foam::interfaceTrackingFvMesh::updateSurfactantConcentration()
                 IOobject
                 (
                     "Cb",
-                    mesh().time().timeName(),
-                    mesh(),
+                    aMesh().time().timeName(),
+                    aMesh().thisDb(),
                     IOobject::NO_READ,
                     IOobject::NO_WRITE
                 ),
@@ -1303,8 +1292,8 @@ void Foam::interfaceTrackingFvMesh::correctContactLinePointNormals()
             IOobject
             (
                 "tangent",
-                mesh().time().timeName(),
-                mesh(),
+                aMesh().time().timeName(),
+                aMesh().thisDb(),
                 IOobject::NO_READ,
                 IOobject::NO_WRITE
             ),
@@ -2051,8 +2040,8 @@ Foam::interfaceTrackingFvMesh::surfaceTensionGrad()
         IOobject
         (
             "surfaceTensionGrad",
-            mesh().time().timeName(),
-            mesh(),
+            aMesh().time().timeName(),
+            aMesh().thisDb(),
             IOobject::NO_READ,
             IOobject::NO_WRITE
         ),
diff --git a/src/finiteArea/faMesh/faPatches/constraint/processor/processorFaPatch.C b/src/finiteArea/faMesh/faPatches/constraint/processor/processorFaPatch.C
index 3312cb8173..d19a16e143 100644
--- a/src/finiteArea/faMesh/faPatches/constraint/processor/processorFaPatch.C
+++ b/src/finiteArea/faMesh/faPatches/constraint/processor/processorFaPatch.C
@@ -135,12 +135,11 @@ void Foam::processorFaPatch::makeNonGlobalPatchPoints() const
     // create a 1->1 map
 
     // Can not use faGlobalMeshData at this point yet
+    // - use polyMesh globalData instead
 
-    if
-    (
-        !Pstream::parRun()
-     || !boundaryMesh().mesh()().globalData().nGlobalPoints()
-    )
+    const auto& pMeshGlobalData = boundaryMesh().mesh().mesh().globalData();
+
+    if (!Pstream::parRun() || !pMeshGlobalData.nGlobalPoints())
     {
         // 1 -> 1 mapping
         nonGlobalPatchPointsPtr_.reset(new labelList(identity(nPoints())));
@@ -151,8 +150,7 @@ void Foam::processorFaPatch::makeNonGlobalPatchPoints() const
         labelList& ngpp = *nonGlobalPatchPointsPtr_;
 
         // Get reference to shared points
-        const labelList& sharedPoints =
-            boundaryMesh().mesh()().globalData().sharedPointLabels();
+        const labelList& sharedPoints = pMeshGlobalData.sharedPointLabels();
 
         const labelList& faMeshPatchPoints = pointLabels();
 
diff --git a/src/finiteArea/finiteArea/d2dt2Schemes/EulerFaD2dt2Scheme/EulerFaD2dt2Scheme.C b/src/finiteArea/finiteArea/d2dt2Schemes/EulerFaD2dt2Scheme/EulerFaD2dt2Scheme.C
index 41a085789f..9ab1276224 100644
--- a/src/finiteArea/finiteArea/d2dt2Schemes/EulerFaD2dt2Scheme/EulerFaD2dt2Scheme.C
+++ b/src/finiteArea/finiteArea/d2dt2Schemes/EulerFaD2dt2Scheme/EulerFaD2dt2Scheme.C
@@ -75,8 +75,8 @@ EulerFaD2dt2Scheme<Type>::facD2dt2
     IOobject d2dt2IOobject
     (
         "d2dt2("+dt.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -134,8 +134,8 @@ EulerFaD2dt2Scheme<Type>::facD2dt2
     IOobject d2dt2IOobject
     (
         "d2dt2("+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -212,8 +212,8 @@ EulerFaD2dt2Scheme<Type>::facD2dt2
     IOobject d2dt2IOobject
     (
         "d2dt2("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -297,8 +297,8 @@ EulerFaD2dt2Scheme<Type>::facD2dt2
     IOobject d2dt2IOobject
     (
         "d2dt2("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
diff --git a/src/finiteArea/finiteArea/ddtSchemes/EulerFaDdtScheme/EulerFaDdtScheme.C b/src/finiteArea/finiteArea/ddtSchemes/EulerFaDdtScheme/EulerFaDdtScheme.C
index 254cbc7570..cd75c7d092 100644
--- a/src/finiteArea/finiteArea/ddtSchemes/EulerFaDdtScheme/EulerFaDdtScheme.C
+++ b/src/finiteArea/finiteArea/ddtSchemes/EulerFaDdtScheme/EulerFaDdtScheme.C
@@ -53,8 +53,8 @@ EulerFaDdtScheme<Type>::facDdt
     IOobject ddtIOobject
     (
         "ddt("+dt.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -100,8 +100,8 @@ EulerFaDdtScheme<Type>::facDdt0
     IOobject ddtIOobject
     (
         "ddt("+dt.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -138,8 +138,8 @@ EulerFaDdtScheme<Type>::facDdt
     IOobject ddtIOobject
     (
         "ddt("+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -191,8 +191,8 @@ EulerFaDdtScheme<Type>::facDdt0
     IOobject ddtIOobject
     (
         "ddt0("+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -238,8 +238,8 @@ EulerFaDdtScheme<Type>::facDdt
     IOobject ddtIOobject
     (
         "ddt("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -291,8 +291,8 @@ EulerFaDdtScheme<Type>::facDdt0
     IOobject ddtIOobject
     (
         "ddt0("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -340,8 +340,8 @@ EulerFaDdtScheme<Type>::facDdt
     IOobject ddtIOobject
     (
         "ddt("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -397,8 +397,8 @@ EulerFaDdtScheme<Type>::facDdt0
     IOobject ddtIOobject
     (
         "ddt0("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
diff --git a/src/finiteArea/finiteArea/ddtSchemes/backwardFaDdtScheme/backwardFaDdtScheme.C b/src/finiteArea/finiteArea/ddtSchemes/backwardFaDdtScheme/backwardFaDdtScheme.C
index 6256f61088..da042c24e5 100644
--- a/src/finiteArea/finiteArea/ddtSchemes/backwardFaDdtScheme/backwardFaDdtScheme.C
+++ b/src/finiteArea/finiteArea/ddtSchemes/backwardFaDdtScheme/backwardFaDdtScheme.C
@@ -84,8 +84,8 @@ backwardFaDdtScheme<Type>::facDdt
     IOobject ddtIOobject
     (
         "ddt("+dt.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -139,8 +139,8 @@ backwardFaDdtScheme<Type>::facDdt0
     IOobject ddtIOobject
     (
         "ddt("+dt.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -186,8 +186,8 @@ backwardFaDdtScheme<Type>::facDdt
     IOobject ddtIOobject
     (
         "ddt("+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -259,8 +259,8 @@ backwardFaDdtScheme<Type>::facDdt0
     IOobject ddtIOobject
     (
         "ddt0("+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -330,8 +330,8 @@ backwardFaDdtScheme<Type>::facDdt
     IOobject ddtIOobject
     (
         "ddt("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -403,8 +403,8 @@ backwardFaDdtScheme<Type>::facDdt0
     IOobject ddtIOobject
     (
         "ddt0("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -474,8 +474,8 @@ backwardFaDdtScheme<Type>::facDdt
     IOobject ddtIOobject
     (
         "ddt("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -551,8 +551,8 @@ backwardFaDdtScheme<Type>::facDdt0
     IOobject ddtIOobject
     (
         "ddt0("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
diff --git a/src/finiteArea/finiteArea/ddtSchemes/boundedBackwardFaDdtScheme/boundedBackwardFaDdtScheme.C b/src/finiteArea/finiteArea/ddtSchemes/boundedBackwardFaDdtScheme/boundedBackwardFaDdtScheme.C
index 22c50c9c3c..140e452e3e 100644
--- a/src/finiteArea/finiteArea/ddtSchemes/boundedBackwardFaDdtScheme/boundedBackwardFaDdtScheme.C
+++ b/src/finiteArea/finiteArea/ddtSchemes/boundedBackwardFaDdtScheme/boundedBackwardFaDdtScheme.C
@@ -67,8 +67,8 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt
     IOobject ddtIOobject
     (
         "ddt("+dt.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -122,8 +122,8 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt0
     IOobject ddtIOobject
     (
         "ddt("+dt.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -167,8 +167,8 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt
     IOobject ddtIOobject
     (
         "ddt("+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -264,8 +264,8 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt0
     IOobject ddtIOobject
     (
         "ddt0("+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -359,8 +359,8 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt
     IOobject ddtIOobject
     (
         "ddt("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -456,8 +456,8 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt0
     IOobject ddtIOobject
     (
         "ddt0("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -551,8 +551,8 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt
     IOobject ddtIOobject
     (
         "ddt("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -652,8 +652,8 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt0
     IOobject ddtIOobject
     (
         "ddt0("+rho.name()+','+vf.name()+')',
-        mesh()().time().timeName(),
-        mesh()(),
+        mesh().time().timeName(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
diff --git a/src/finiteArea/finiteArea/ddtSchemes/steadyStateFaDdtScheme/steadyStateFaDdtScheme.C b/src/finiteArea/finiteArea/ddtSchemes/steadyStateFaDdtScheme/steadyStateFaDdtScheme.C
index 1a44f226c1..d430680d88 100644
--- a/src/finiteArea/finiteArea/ddtSchemes/steadyStateFaDdtScheme/steadyStateFaDdtScheme.C
+++ b/src/finiteArea/finiteArea/ddtSchemes/steadyStateFaDdtScheme/steadyStateFaDdtScheme.C
@@ -53,8 +53,8 @@ steadyStateFaDdtScheme<Type>::facDdt
         IOobject
         (
             "ddt("+dt.name()+')',
-            mesh()().time().timeName(),
-            mesh()()
+            mesh().time().timeName(),
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(dt.dimensions()/dimTime, Zero)
@@ -74,8 +74,8 @@ steadyStateFaDdtScheme<Type>::facDdt0
         IOobject
         (
             "ddt("+dt.name()+')',
-            mesh()().time().timeName(),
-            mesh()()
+            mesh().time().timeName(),
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(dt.dimensions()/dimTime, Zero)
@@ -95,8 +95,8 @@ steadyStateFaDdtScheme<Type>::facDdt
         IOobject
         (
             "ddt("+vf.name()+')',
-            mesh()().time().timeName(),
-            mesh()()
+            mesh().time().timeName(),
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(vf.dimensions()/dimTime, Zero)
@@ -116,8 +116,8 @@ steadyStateFaDdtScheme<Type>::facDdt0
         IOobject
         (
             "ddt0("+vf.name()+')',
-            mesh()().time().timeName(),
-            mesh()()
+            mesh().time().timeName(),
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(vf.dimensions()/dimTime, Zero)
@@ -138,8 +138,8 @@ steadyStateFaDdtScheme<Type>::facDdt
         IOobject
         (
             "ddt("+rho.name()+','+vf.name()+')',
-            mesh()().time().timeName(),
-            mesh()()
+            mesh().time().timeName(),
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
@@ -159,8 +159,8 @@ steadyStateFaDdtScheme<Type>::facDdt0
         IOobject
         (
             "ddt0("+rho.name()+','+vf.name()+')',
-            mesh()().time().timeName(),
-            mesh()()
+            mesh().time().timeName(),
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
@@ -181,8 +181,8 @@ steadyStateFaDdtScheme<Type>::facDdt
         IOobject
         (
             "ddt("+rho.name()+','+vf.name()+')',
-            mesh()().time().timeName(),
-            mesh()()
+            mesh().time().timeName(),
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
@@ -203,8 +203,8 @@ steadyStateFaDdtScheme<Type>::facDdt0
         IOobject
         (
             "ddt0("+rho.name()+','+vf.name()+')',
-            mesh()().time().timeName(),
-            mesh()()
+            mesh().time().timeName(),
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
diff --git a/src/finiteArea/interpolation/edgeInterpolation/edgeInterpolation.C b/src/finiteArea/interpolation/edgeInterpolation/edgeInterpolation.C
index b57f0e75f2..a75f322260 100644
--- a/src/finiteArea/interpolation/edgeInterpolation/edgeInterpolation.C
+++ b/src/finiteArea/interpolation/edgeInterpolation/edgeInterpolation.C
@@ -205,8 +205,8 @@ void Foam::edgeInterpolation::makeLPN() const
         IOobject
         (
             "lPN",
-            faMesh_.time().constant(),
-            faMesh_(),
+            mesh().time().constant(),
+            mesh().thisDb(),
             IOobject::NO_READ,
             IOobject::NO_WRITE,
             false
@@ -286,8 +286,8 @@ void Foam::edgeInterpolation::makeWeights() const
         IOobject
         (
             "weightingFactors",
-            mesh()().pointsInstance(),
-            mesh()(),
+            mesh().pointsInstance(),
+            mesh().thisDb(),
             IOobject::NO_READ,
             IOobject::NO_WRITE,
             false
@@ -366,8 +366,8 @@ void Foam::edgeInterpolation::makeDeltaCoeffs() const
         IOobject
         (
             "differenceFactors_",
-            mesh()().pointsInstance(),
-            mesh()(),
+            mesh().pointsInstance(),
+            mesh().thisDb(),
             IOobject::NO_READ,
             IOobject::NO_WRITE,
             false
@@ -460,8 +460,8 @@ void Foam::edgeInterpolation::makeCorrectionVectors() const
         IOobject
         (
             "correctionVectors",
-            mesh()().pointsInstance(),
-            mesh()(),
+            mesh().pointsInstance(),
+            mesh().thisDb(),
             IOobject::NO_READ,
             IOobject::NO_WRITE,
             false
@@ -571,8 +571,8 @@ void Foam::edgeInterpolation::makeSkewCorrectionVectors() const
         IOobject
         (
             "skewCorrectionVectors",
-            mesh()().pointsInstance(),
-            mesh()(),
+            mesh().pointsInstance(),
+            mesh().thisDb(),
             IOobject::NO_READ,
             IOobject::NO_WRITE,
             false
diff --git a/src/finiteVolume/cfdTools/general/MRF/MRFZoneList.C b/src/finiteVolume/cfdTools/general/MRF/MRFZoneList.C
index 54ec416cde..c3b4084131 100644
--- a/src/finiteVolume/cfdTools/general/MRF/MRFZoneList.C
+++ b/src/finiteVolume/cfdTools/general/MRF/MRFZoneList.C
@@ -195,7 +195,7 @@ Foam::tmp<Foam::volVectorField> Foam::MRFZoneList::DDt
             (
                 "MRFZoneList:acceleration",
                 U.mesh().time().timeName(),
-                U.mesh()
+                U.mesh().thisDb()
             ),
             U.mesh(),
             dimensionedVector(U.dimensions()/dimTime, Zero)
diff --git a/src/finiteVolume/finiteVolume/d2dt2Schemes/EulerD2dt2Scheme/EulerD2dt2Scheme.C b/src/finiteVolume/finiteVolume/d2dt2Schemes/EulerD2dt2Scheme/EulerD2dt2Scheme.C
index d27011b387..dfd383475c 100644
--- a/src/finiteVolume/finiteVolume/d2dt2Schemes/EulerD2dt2Scheme/EulerD2dt2Scheme.C
+++ b/src/finiteVolume/finiteVolume/d2dt2Schemes/EulerD2dt2Scheme/EulerD2dt2Scheme.C
@@ -55,7 +55,7 @@ EulerD2dt2Scheme<Type>::fvcD2dt2
     (
         "d2dt2("+vf.name()+')',
         mesh().time().timeName(),
-        mesh(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
@@ -133,7 +133,7 @@ EulerD2dt2Scheme<Type>::fvcD2dt2
     (
         "d2dt2("+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh(),
+        mesh().thisDb(),
         IOobject::NO_READ,
         IOobject::NO_WRITE
     );
diff --git a/src/finiteVolume/finiteVolume/d2dt2Schemes/steadyStateD2dt2Scheme/steadyStateD2dt2Scheme.C b/src/finiteVolume/finiteVolume/d2dt2Schemes/steadyStateD2dt2Scheme/steadyStateD2dt2Scheme.C
index 3313ba1519..0ff77f5722 100644
--- a/src/finiteVolume/finiteVolume/d2dt2Schemes/steadyStateD2dt2Scheme/steadyStateD2dt2Scheme.C
+++ b/src/finiteVolume/finiteVolume/d2dt2Schemes/steadyStateD2dt2Scheme/steadyStateD2dt2Scheme.C
@@ -56,7 +56,7 @@ steadyStateD2dt2Scheme<Type>::fvcD2dt2
             (
                 "d2dt2("+vf.name()+')',
                 mesh().time().timeName(),
-                mesh(),
+                mesh().thisDb(),
                 IOobject::NO_READ,
                 IOobject::NO_WRITE
             ),
@@ -81,7 +81,7 @@ steadyStateD2dt2Scheme<Type>::fvcD2dt2
         (
             "d2dt2("+rho.name()+','+vf.name()+')',
             mesh().time().timeName(),
-            mesh(),
+            mesh().thisDb(),
             IOobject::NO_READ,
             IOobject::NO_WRITE
         ),
diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/CoEulerDdtScheme/CoEulerDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/CoEulerDdtScheme/CoEulerDdtScheme.C
index 0644e8866d..eff7207c0a 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/CoEulerDdtScheme/CoEulerDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/CoEulerDdtScheme/CoEulerDdtScheme.C
@@ -159,7 +159,7 @@ CoEulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+dt.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -209,7 +209,7 @@ CoEulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -261,7 +261,7 @@ CoEulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -313,7 +313,7 @@ CoEulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -369,7 +369,7 @@ CoEulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+alpha.name()+','+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -601,7 +601,7 @@ CoEulerDdtScheme<Type>::fvcDdtUfCorr
             (
                 "ddtCorr(" + U.name() + ',' + Uf.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), phiUf0, phiCorr)
            *rDeltaT*phiCorr
@@ -633,7 +633,7 @@ CoEulerDdtScheme<Type>::fvcDdtPhiCorr
             (
                 "ddtCorr(" + U.name() + ',' + phi.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), phi.oldTime(), phiCorr)
            *rDeltaT*phiCorr
@@ -676,7 +676,7 @@ CoEulerDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff(rhoU0, phiUf0, phiCorr, rho.oldTime())
                *rDeltaT*phiCorr
@@ -704,7 +704,7 @@ CoEulerDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -763,7 +763,7 @@ CoEulerDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -795,7 +795,7 @@ CoEulerDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -832,7 +832,7 @@ tmp<surfaceScalarField> CoEulerDdtScheme<Type>::meshPhi
             (
                 "meshPhi",
                 mesh().time().timeName(),
-                mesh(),
+                mesh().thisDb(),
                 IOobject::NO_READ,
                 IOobject::NO_WRITE,
                 false
diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicolsonDdtScheme/CrankNicolsonDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicolsonDdtScheme/CrankNicolsonDdtScheme.C
index 663a956f4f..207245e9d2 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicolsonDdtScheme/CrankNicolsonDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/CrankNicolsonDdtScheme/CrankNicolsonDdtScheme.C
@@ -157,7 +157,7 @@ CrankNicolsonDdtScheme<Type>::ddt0_
                     (
                         name,
                         mesh().time().timeName(),
-                        mesh(),
+                        mesh().thisDb(),
                         IOobject::NO_READ,
                         IOobject::AUTO_WRITE
                     ),
@@ -360,7 +360,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdt
     (
         "ddt(" + dt.name() + ')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     tmp<GeometricField<Type, fvPatchField, volMesh>> tdtdt
@@ -417,7 +417,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdt
     (
         "ddt(" + vf.name() + ')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     dimensionedScalar rDtCoef = rDtCoef_(ddt0);
@@ -505,7 +505,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdt
     (
         "ddt(" + rho.name() + ',' + vf.name() + ')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     dimensionedScalar rDtCoef = rDtCoef_(ddt0);
@@ -595,7 +595,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdt
     (
         "ddt(" + rho.name() + ',' + vf.name() + ')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     dimensionedScalar rDtCoef = rDtCoef_(ddt0);
@@ -696,7 +696,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdt
     (
         "ddt(" + alpha.name() + ',' + rho.name() + ',' + vf.name() + ')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     dimensionedScalar rDtCoef = rDtCoef_(ddt0);
@@ -1218,7 +1218,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdtUfCorr
             (
                 "ddtCorr(" + U.name() + ',' + Uf.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), mesh().Sf() & Uf.oldTime())
            *(
@@ -1280,7 +1280,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdtPhiCorr
             (
                 "ddtCorr(" + U.name() + ',' + phi.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), phi.oldTime())
            *(
@@ -1357,7 +1357,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -1421,7 +1421,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdtUfCorr
                 (
                     "ddtCorr(" + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -1514,7 +1514,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff(rhoU0, phi.oldTime(), rho.oldTime())
                *(
@@ -1574,7 +1574,7 @@ CrankNicolsonDdtScheme<Type>::fvcDdtPhiCorr
                 (
                     "ddtCorr(" + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff(U.oldTime(), phi.oldTime(), rho.oldTime())
                *(
@@ -1627,7 +1627,7 @@ tmp<surfaceScalarField> CrankNicolsonDdtScheme<Type>::meshPhi
             (
                 mesh().phi().name(),
                 mesh().time().timeName(),
-                mesh(),
+                mesh().thisDb(),
                 IOobject::NO_READ,
                 IOobject::NO_WRITE,
                 false
diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C
index 7dcf820781..584b9b8426 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C
@@ -55,7 +55,7 @@ EulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+dt.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -99,7 +99,7 @@ EulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -149,7 +149,7 @@ EulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -199,7 +199,7 @@ EulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -253,7 +253,7 @@ EulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+alpha.name()+','+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -317,7 +317,7 @@ EulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+sf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     return tmp<GeometricField<Type, fvsPatchField, surfaceMesh>>
@@ -509,7 +509,7 @@ EulerDdtScheme<Type>::fvcDdtUfCorr
             (
                 "ddtCorr(" + U.name() + ',' + Uf.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), phiUf0, phiCorr)
            *rDeltaT*phiCorr
@@ -541,7 +541,7 @@ EulerDdtScheme<Type>::fvcDdtPhiCorr
             (
                 "ddtCorr(" + U.name() + ',' + phi.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), phi.oldTime(), phiCorr)
            *rDeltaT*phiCorr
@@ -584,7 +584,7 @@ EulerDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff(rhoU0, phiUf0, phiCorr, rho.oldTime())
                *rDeltaT*phiCorr
@@ -612,7 +612,7 @@ EulerDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -671,7 +671,7 @@ EulerDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -703,7 +703,7 @@ EulerDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/SLTSDdtScheme/SLTSDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/SLTSDdtScheme/SLTSDdtScheme.C
index 5d441d8814..38ac5d0412 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/SLTSDdtScheme/SLTSDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/SLTSDdtScheme/SLTSDdtScheme.C
@@ -165,7 +165,7 @@ SLTSDdtScheme<Type>::fvcDdt
     (
         "ddt("+dt.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -214,7 +214,7 @@ SLTSDdtScheme<Type>::fvcDdt
     (
         "ddt("+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -266,7 +266,7 @@ SLTSDdtScheme<Type>::fvcDdt
     (
         "ddt("+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -318,7 +318,7 @@ SLTSDdtScheme<Type>::fvcDdt
     (
         "ddt("+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -374,7 +374,7 @@ SLTSDdtScheme<Type>::fvcDdt
     (
         "ddt("+alpha.name()+','+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     if (mesh().moving())
@@ -606,7 +606,7 @@ SLTSDdtScheme<Type>::fvcDdtUfCorr
             (
                 "ddtCorr(" + U.name() + ',' + Uf.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), phiUf0, phiCorr)
            *rDeltaT*phiCorr
@@ -638,7 +638,7 @@ SLTSDdtScheme<Type>::fvcDdtPhiCorr
             (
                 "ddtCorr(" + U.name() + ',' + phi.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), phi.oldTime(), phiCorr)
            *rDeltaT*phiCorr
@@ -681,7 +681,7 @@ SLTSDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff(rhoU0, phiUf0, phiCorr, rho.oldTime())
                *rDeltaT*phiCorr
@@ -709,7 +709,7 @@ SLTSDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -768,7 +768,7 @@ SLTSDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -800,7 +800,7 @@ SLTSDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -837,7 +837,7 @@ tmp<surfaceScalarField> SLTSDdtScheme<Type>::meshPhi
             (
                 "meshPhi",
                 mesh().time().timeName(),
-                mesh(),
+                mesh().thisDb(),
                 IOobject::NO_READ,
                 IOobject::NO_WRITE,
                 false
diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/backwardDdtScheme/backwardDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/backwardDdtScheme/backwardDdtScheme.C
index b81c180ca2..e470ec87c8 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/backwardDdtScheme/backwardDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/backwardDdtScheme/backwardDdtScheme.C
@@ -87,7 +87,7 @@ backwardDdtScheme<Type>::fvcDdt
     (
         "ddt("+dt.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     scalar deltaT = deltaT_();
@@ -140,7 +140,7 @@ backwardDdtScheme<Type>::fvcDdt
     (
         "ddt("+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     scalar deltaT = deltaT_();
@@ -212,7 +212,7 @@ backwardDdtScheme<Type>::fvcDdt
     (
         "ddt("+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     scalar deltaT = deltaT_();
@@ -284,7 +284,7 @@ backwardDdtScheme<Type>::fvcDdt
     (
         "ddt("+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     scalar deltaT = deltaT_();
@@ -360,7 +360,7 @@ backwardDdtScheme<Type>::fvcDdt
     (
         "ddt("+alpha.name()+','+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     scalar deltaT = deltaT_();
@@ -686,7 +686,7 @@ backwardDdtScheme<Type>::fvcDdtUfCorr
             (
                 "ddtCorr(" + U.name() + ',' + Uf.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), (mesh().Sf() & Uf.oldTime()))
            *rDeltaT
@@ -730,7 +730,7 @@ backwardDdtScheme<Type>::fvcDdtPhiCorr
             (
                 "ddtCorr(" + U.name() + ',' + phi.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), phi.oldTime())
            *rDeltaT
@@ -790,7 +790,7 @@ backwardDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -824,7 +824,7 @@ backwardDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -901,7 +901,7 @@ backwardDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff(rhoU0, phi.oldTime(), rho.oldTime())
                *rDeltaT
@@ -931,7 +931,7 @@ backwardDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff(U.oldTime(), phi.oldTime(), rho.oldTime())
                *rDeltaT
@@ -980,7 +980,7 @@ tmp<surfaceScalarField> backwardDdtScheme<Type>::meshPhi
             (
                 mesh().phi().name(),
                 mesh().time().timeName(),
-                mesh(),
+                mesh().thisDb(),
                 IOobject::NO_READ,
                 IOobject::NO_WRITE,
                 false
diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/ddtScheme/ddtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/ddtScheme/ddtScheme.C
index aeb6a14587..94a00c2a68 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/ddtScheme/ddtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/ddtScheme/ddtScheme.C
@@ -160,7 +160,7 @@ tmp<surfaceScalarField> ddtScheme<Type>::fvcDdtPhiCoeff
             (
                 "ddtCouplingCoeff",
                 U.mesh().time().timeName(),
-                U.mesh()
+                U.mesh().thisDb()
             ),
             U.mesh(),
             dimensionedScalar("one", dimless, 1.0)
@@ -234,7 +234,7 @@ tmp<surfaceScalarField> ddtScheme<Type>::fvcDdtPhiCoeffExperimental
             (
                 "ddtCouplingCoeff",
                 U.mesh().time().timeName(),
-                U.mesh()
+                U.mesh().thisDb()
             ),
             U.mesh(),
             dimensionedScalar("one", dimless, 1.0)
diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/localEulerDdtScheme/localEulerDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/localEulerDdtScheme/localEulerDdtScheme.C
index 539c59818b..bf3979fc21 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/localEulerDdtScheme/localEulerDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/localEulerDdtScheme/localEulerDdtScheme.C
@@ -66,7 +66,7 @@ localEulerDdtScheme<Type>::fvcDdt
     (
         "ddt(" + dt.name() + ')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     return tmp<GeometricField<Type, fvPatchField, volMesh>>::New
@@ -92,7 +92,7 @@ localEulerDdtScheme<Type>::fvcDdt
     (
         "ddt(" + vf.name() + ')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     return tmp<GeometricField<Type, fvPatchField, volMesh>>
@@ -120,7 +120,7 @@ localEulerDdtScheme<Type>::fvcDdt
     (
         "ddt(" + rho.name() + ',' + vf.name() + ')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     return tmp<GeometricField<Type, fvPatchField, volMesh>>
@@ -148,7 +148,7 @@ localEulerDdtScheme<Type>::fvcDdt
     (
         "ddt(" + rho.name() + ',' + vf.name() + ')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     return tmp<GeometricField<Type, fvPatchField, volMesh>>
@@ -177,7 +177,7 @@ localEulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+alpha.name()+','+rho.name()+','+vf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     return tmp<GeometricField<Type, fvPatchField, volMesh>>
@@ -208,7 +208,7 @@ localEulerDdtScheme<Type>::fvcDdt
     (
         "ddt("+sf.name()+')',
         mesh().time().timeName(),
-        mesh()
+        mesh().thisDb()
     );
 
     return tmp<GeometricField<Type, fvsPatchField, surfaceMesh>>
@@ -365,7 +365,7 @@ localEulerDdtScheme<Type>::fvcDdtUfCorr
             (
                 "ddtCorr(" + U.name() + ',' + Uf.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), phiUf0, phiCorr)
            *rDeltaT*phiCorr
@@ -397,7 +397,7 @@ localEulerDdtScheme<Type>::fvcDdtPhiCorr
             (
                 "ddtCorr(" + U.name() + ',' + phi.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             this->fvcDdtPhiCoeff(U.oldTime(), phi.oldTime(), phiCorr)
            *rDeltaT*phiCorr
@@ -440,7 +440,7 @@ localEulerDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff(rhoU0, phiUf0, phiCorr, rho.oldTime())
                *rDeltaT*phiCorr
@@ -468,7 +468,7 @@ localEulerDdtScheme<Type>::fvcDdtUfCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -527,7 +527,7 @@ localEulerDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -559,7 +559,7 @@ localEulerDdtScheme<Type>::fvcDdtPhiCorr
                     "ddtCorr("
                   + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
                     mesh().time().timeName(),
-                    mesh()
+                    mesh().thisDb()
                 ),
                 this->fvcDdtPhiCoeff
                 (
@@ -596,7 +596,7 @@ tmp<surfaceScalarField> localEulerDdtScheme<Type>::meshPhi
             (
                 "meshPhi",
                 mesh().time().timeName(),
-                mesh(),
+                mesh().thisDb(),
                 IOobject::NO_READ,
                 IOobject::NO_WRITE,
                 false
diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/steadyStateDdtScheme/steadyStateDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/steadyStateDdtScheme/steadyStateDdtScheme.C
index 9c44f21d1a..33fb99b030 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/steadyStateDdtScheme/steadyStateDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/steadyStateDdtScheme/steadyStateDdtScheme.C
@@ -54,7 +54,7 @@ steadyStateDdtScheme<Type>::fvcDdt
         (
             "ddt("+dt.name()+')',
             mesh().time().timeName(),
-            mesh()
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(dt.dimensions()/dimTime, Zero)
@@ -75,7 +75,7 @@ steadyStateDdtScheme<Type>::fvcDdt
         (
             "ddt("+vf.name()+')',
             mesh().time().timeName(),
-            mesh()
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(vf.dimensions()/dimTime, Zero)
@@ -97,7 +97,7 @@ steadyStateDdtScheme<Type>::fvcDdt
         (
             "ddt("+rho.name()+','+vf.name()+')',
             mesh().time().timeName(),
-            mesh()
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
@@ -119,7 +119,7 @@ steadyStateDdtScheme<Type>::fvcDdt
         (
             "ddt("+rho.name()+','+vf.name()+')',
             mesh().time().timeName(),
-            mesh()
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
@@ -142,7 +142,7 @@ steadyStateDdtScheme<Type>::fvcDdt
         (
             "ddt("+alpha.name()+','+rho.name()+','+vf.name()+')',
             mesh().time().timeName(),
-            mesh()
+            mesh().thisDb()
         ),
         mesh(),
         dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
@@ -250,7 +250,7 @@ steadyStateDdtScheme<Type>::fvcDdtUfCorr
             (
                 "ddtCorr(" + U.name() + ',' + Uf.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             mesh(),
             dimensioned<typename flux<Type>::type>
@@ -282,7 +282,7 @@ steadyStateDdtScheme<Type>::fvcDdtPhiCorr
             (
                 "ddtCorr(" + U.name() + ',' + phi.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             mesh(),
             dimensioned<typename flux<Type>::type>
@@ -317,7 +317,7 @@ steadyStateDdtScheme<Type>::fvcDdtUfCorr
               + rho.name()
               + ',' + U.name() + ',' + Uf.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             mesh(),
             dimensioned<typename flux<Type>::type>
@@ -352,7 +352,7 @@ steadyStateDdtScheme<Type>::fvcDdtPhiCorr
               + rho.name()
               + ',' + U.name() + ',' + phi.name() + ')',
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             mesh(),
             dimensioned<typename flux<Type>::type>
@@ -380,7 +380,7 @@ tmp<surfaceScalarField> steadyStateDdtScheme<Type>::meshPhi
         (
             "meshPhi",
             mesh().time().timeName(),
-            mesh(),
+            mesh().thisDb(),
             IOobject::NO_READ,
             IOobject::NO_WRITE,
             false
diff --git a/src/finiteVolume/finiteVolume/snGradSchemes/limitedSnGrad/limitedSnGrad.C b/src/finiteVolume/finiteVolume/snGradSchemes/limitedSnGrad/limitedSnGrad.C
index 6f0f95ee35..addf8e45e6 100644
--- a/src/finiteVolume/finiteVolume/snGradSchemes/limitedSnGrad/limitedSnGrad.C
+++ b/src/finiteVolume/finiteVolume/snGradSchemes/limitedSnGrad/limitedSnGrad.C
@@ -100,7 +100,7 @@ limitedSnGrad<Type>::correction
                 (
                     fieldName,
                     vf.mesh().time().timeName(),
-                    vf.mesh(),
+                    vf.mesh().thisDb(),
                     IOobject::NO_READ,
                     IOobject::NO_WRITE,
                     false
diff --git a/src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C b/src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C
index a82341e6ce..38dd6520a8 100644
--- a/src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C
+++ b/src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C
@@ -926,9 +926,10 @@ void Foam::fvMatrix<Type>::createOrUpdateLduPrimitiveAssembly()
     (
         lduAssemblyName_,
         psi_.mesh().time().timeName(),
-        psi_.mesh(),
+        psi_.mesh().thisDb(),
         IOobject::NO_READ,
-        IOobject::NO_WRITE
+        IOobject::NO_WRITE,
+        IOobject::REGISTER
     );
 
     UPtrList<lduMesh> uMeshPtr(nMatrices());
diff --git a/src/finiteVolume/fvMesh/wallDist/wallDist/wallDist.C b/src/finiteVolume/fvMesh/wallDist/wallDist/wallDist.C
index 6e417d37a4..be8b859bb7 100644
--- a/src/finiteVolume/fvMesh/wallDist/wallDist/wallDist.C
+++ b/src/finiteVolume/fvMesh/wallDist/wallDist/wallDist.C
@@ -49,7 +49,7 @@ void Foam::wallDist::constructn() const
             (
                 "n" & patchTypeName_,
                 mesh().time().timeName(),
-                mesh()
+                mesh().thisDb()
             ),
             mesh(),
             dimensionedVector(dimless, Zero),
@@ -115,7 +115,7 @@ Foam::wallDist::wallDist
         (
             "y" & patchTypeName_,
             mesh.time().timeName(),
-            mesh
+            mesh.thisDb()
         ),
         mesh,
         dimensionedScalar("y" & patchTypeName_, dimLength, SMALL),
diff --git a/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/blended/blended.H b/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/blended/blended.H
index 239d040f96..e99bf4889c 100644
--- a/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/blended/blended.H
+++ b/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/blended/blended.H
@@ -139,7 +139,7 @@ public:
                     (
                         vf.name() + "BlendingFactor",
                         this->mesh().time().timeName(),
-                        this->mesh(),
+                        this->mesh().thisDb(),
                         IOobject::NO_READ,
                         IOobject::NO_WRITE,
                         false
@@ -169,7 +169,7 @@ public:
                     (
                         "blendedLimiter",
                         this->mesh().time().timeName(),
-                        this->mesh(),
+                        this->mesh().thisDb(),
                         IOobject::NO_READ,
                         IOobject::NO_WRITE,
                         false
diff --git a/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/upwind/upwind.H b/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/upwind/upwind.H
index 2ba3bfb3f4..3831b19ccd 100644
--- a/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/upwind/upwind.H
+++ b/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/upwind/upwind.H
@@ -120,7 +120,8 @@ public:
                 (
                     "upwindLimiter",
                     this->mesh().time().timeName(),
-                    this->mesh(),
+                    this->mesh().time().timeName(),
+                    this->mesh().thisDb(),
                     IOobject::NO_READ,
                     IOobject::NO_WRITE,
                     false
diff --git a/src/finiteVolume/interpolation/surfaceInterpolation/schemes/midPoint/midPoint.H b/src/finiteVolume/interpolation/surfaceInterpolation/schemes/midPoint/midPoint.H
index 73e0350d9f..18d413df10 100644
--- a/src/finiteVolume/interpolation/surfaceInterpolation/schemes/midPoint/midPoint.H
+++ b/src/finiteVolume/interpolation/surfaceInterpolation/schemes/midPoint/midPoint.H
@@ -111,7 +111,7 @@ public:
                     (
                         "midPointWeights",
                         this->mesh().time().timeName(),
-                        this->mesh(),
+                        this->mesh().thisDb(),
                         IOobject::NO_READ,
                         IOobject::NO_WRITE,
                         false