diff --git a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
index a0c76f5abe..e6a7388f4a 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
@@ -106,7 +106,8 @@ int main(int argc, char *argv[])
 
                 if (correctPhi)
                 {
-                    divU = volScalarField::New
+                    // Construct and register divU for mapping
+                    divU = new volScalarField
                     (
                         "divU0",
                         fvc::div(fvc::absolute(phi, U))
diff --git a/applications/solvers/multiphase/interFoam/interFoam.C b/applications/solvers/multiphase/interFoam/interFoam.C
index 264c3bd578..ad4166825f 100644
--- a/applications/solvers/multiphase/interFoam/interFoam.C
+++ b/applications/solvers/multiphase/interFoam/interFoam.C
@@ -111,7 +111,8 @@ int main(int argc, char *argv[])
                  && !isType<twoPhaseChangeModels::noPhaseChange>(phaseChange)
                 )
                 {
-                    divU = volScalarField::New
+                    // Construct and register divU for mapping
+                    divU = new volScalarField
                     (
                         "divU0",
                         fvc::div(fvc::absolute(phi, U))
diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam/multiphaseEulerFoam.C b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam/multiphaseEulerFoam.C
index 15ce43cb3b..c2c48e9a50 100644
--- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam/multiphaseEulerFoam.C
+++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam/multiphaseEulerFoam.C
@@ -136,7 +136,8 @@ int main(int argc, char *argv[])
                         correctPhi
                     )
                     {
-                        divU = volScalarField::New
+                        // Construct and register divU for mapping
+                        divU = new volScalarField
                         (
                             "divU0",
                             fvc::div