diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/BlendedInterfacialModel/BlendedInterfacialModel.C b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/BlendedInterfacialModel/BlendedInterfacialModel.C
index 194a224d3e9df4661fa1c33a0e3e95aea1e3896c..e02f2c8e23e5dd9ff6a40a181278d316971317ee 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/BlendedInterfacialModel/BlendedInterfacialModel.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/BlendedInterfacialModel/BlendedInterfacialModel.C
@@ -270,22 +270,19 @@ Foam::BlendedInterfacialModel<modelType>::F() const
         f2 = blending_.f2(pair1In2_.dispersed(), pair2In1_.dispersed());
     }
 
-    tmp<GeometricField<Type, fvPatchField, volMesh>> x
+    auto x = tmp<GeometricField<Type, fvPatchField, volMesh>>::New
     (
-        new GeometricField<Type, fvPatchField, volMesh>
+        IOobject
         (
-            IOobject
-            (
-                modelType::typeName + ":F",
-                pair_.phase1().mesh().time().timeName(),
-                pair_.phase1().mesh(),
-                IOobject::NO_READ,
-                IOobject::NO_WRITE,
-                false
-            ),
+            modelType::typeName + ":F",
+            pair_.phase1().mesh().time().timeName(),
             pair_.phase1().mesh(),
-            dimensioned<Type>(modelType::dimF, Zero)
-        )
+            IOobject::NO_READ,
+            IOobject::NO_WRITE,
+            false
+        ),
+        pair_.phase1().mesh(),
+        dimensioned<Type>(modelType::dimF, Zero)
     );
 
     if (model_.valid())
@@ -338,24 +335,23 @@ Foam::BlendedInterfacialModel<modelType>::Ff() const
         );
     }
 
-    tmp<surfaceScalarField> x
+    auto x = tmp<surfaceScalarField>::New
     (
-        new surfaceScalarField
+        IOobject
         (
-            IOobject
-            (
-                modelType::typeName + ":Ff",
-                pair_.phase1().mesh().time().timeName(),
-                pair_.phase1().mesh(),
-                IOobject::NO_READ,
-                IOobject::NO_WRITE,
-                false
-            ),
+            modelType::typeName + ":Ff",
+            pair_.phase1().mesh().time().timeName(),
             pair_.phase1().mesh(),
-            dimensionedScalar(modelType::dimF*dimArea, Zero)
-        )
+            IOobject::NO_READ,
+            IOobject::NO_WRITE,
+            false
+        ),
+        pair_.phase1().mesh(),
+        dimensionedScalar(modelType::dimF*dimArea, Zero)
     );
 
+    x.ref().setOriented();
+
     if (model_.valid())
     {
         x.ref() += model_->Ff()*(f1() - f2());