From e3298597e1a72ec2ed2f8698bfbb1e7e1cdffb75 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Sat, 4 Jun 2022 11:14:26 +0100
Subject: [PATCH] compressibleInterFoam: Register dgdt to ensure it is mapped
 after mesh topology change

---
 .../multiphase/compressibleInterFoam/createFields.H | 13 ++++++++++++-
 .../phaseModel/phaseModel.H                         |  8 ++++----
 2 files changed, 16 insertions(+), 5 deletions(-)

diff --git a/applications/solvers/multiphase/compressibleInterFoam/createFields.H b/applications/solvers/multiphase/compressibleInterFoam/createFields.H
index cc4415c4a4..9a2cd67713 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/createFields.H
+++ b/applications/solvers/multiphase/compressibleInterFoam/createFields.H
@@ -76,7 +76,18 @@ surfaceScalarField rhoPhi
     fvc::interpolate(rho)*phi
 );
 
-volScalarField dgdt(alpha1*fvc::div(phi));
+volScalarField::Internal dgdt
+(
+    IOobject
+    (
+        "dgdt",
+        runTime.timeName(),
+        mesh,
+        IOobject::READ_IF_PRESENT,
+        IOobject::AUTO_WRITE
+    ),
+    alpha1*fvc::div(phi)
+);
 
 #include "createAlphaFluxes.H"
 
diff --git a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/compressibleMultiphaseMixture/phaseModel/phaseModel.H b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/compressibleMultiphaseMixture/phaseModel/phaseModel.H
index aa0389dfd9..20f1bb9a2f 100644
--- a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/compressibleMultiphaseMixture/phaseModel/phaseModel.H
+++ b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/compressibleMultiphaseMixture/phaseModel/phaseModel.H
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2013-2021 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2013-2022 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -62,7 +62,7 @@ class phaseModel
         const volScalarField& T_;
         autoPtr<rhoThermo> thermo_;
         volScalarField Alpha_;
-        volScalarField dgdt_;
+        volScalarField::Internal dgdt_;
 
 
 public:
@@ -143,13 +143,13 @@ public:
         }
 
         //- Return const-access to phase divergence
-        const volScalarField& dgdt() const
+        const volScalarField::Internal& dgdt() const
         {
             return dgdt_;
         }
 
         //- Return access to phase divergence
-        volScalarField& dgdt()
+        volScalarField::Internal& dgdt()
         {
             return dgdt_;
         }