diff --git a/applications/solvers/modules/compressibleVoF/compressibleInterPhaseThermophysicalTransportModel/compressibleInterPhaseThermophysicalTransportModel.C b/applications/solvers/modules/compressibleVoF/compressibleInterPhaseThermophysicalTransportModel/compressibleInterPhaseThermophysicalTransportModel.C
index 7936f8eef9..5cc6633faf 100644
--- a/applications/solvers/modules/compressibleVoF/compressibleInterPhaseThermophysicalTransportModel/compressibleInterPhaseThermophysicalTransportModel.C
+++ b/applications/solvers/modules/compressibleVoF/compressibleInterPhaseThermophysicalTransportModel/compressibleInterPhaseThermophysicalTransportModel.C
@@ -182,9 +182,12 @@ Foam::compressibleInterPhaseThermophysicalTransportModel::alphaEff() const
 Foam::tmp<Foam::surfaceScalarField>
 Foam::compressibleInterPhaseThermophysicalTransportModel::q() const
 {
-    NotImplemented;
-
-    return tmp<surfaceScalarField>(nullptr);
+    return surfaceScalarField::New
+    (
+        "q",
+        -fvc::interpolate(kappaEff())
+        *fvc::snGrad(momentumTransport_.mixture_.T())
+    );
 }
 
 
diff --git a/applications/solvers/modules/compressibleVoF/compressibleInterPhaseThermophysicalTransportModel/compressibleInterPhaseThermophysicalTransportModel.H b/applications/solvers/modules/compressibleVoF/compressibleInterPhaseThermophysicalTransportModel/compressibleInterPhaseThermophysicalTransportModel.H
index 9c240b0aa5..45aa0d4eac 100644
--- a/applications/solvers/modules/compressibleVoF/compressibleInterPhaseThermophysicalTransportModel/compressibleInterPhaseThermophysicalTransportModel.H
+++ b/applications/solvers/modules/compressibleVoF/compressibleInterPhaseThermophysicalTransportModel/compressibleInterPhaseThermophysicalTransportModel.H
@@ -25,6 +25,10 @@ Class
     Foam::compressibleInterPhaseThermophysicalTransportModel
 
 Description
+    Compressible two-phase VoF thermophysicalTransportModel
+
+    Provides the mixture thermal conductivity and diffusivity for
+    Fourier's temperature gradient heat flux model.
 
 SourceFiles
     compressibleInterPhaseThermophysicalTransportModel.C
diff --git a/src/twoPhaseModels/compressibleTwoPhaseMixture/compressibleTwoPhaseMixture.H b/src/twoPhaseModels/compressibleTwoPhaseMixture/compressibleTwoPhaseMixture.H
index ba064a6d44..ba313639e3 100644
--- a/src/twoPhaseModels/compressibleTwoPhaseMixture/compressibleTwoPhaseMixture.H
+++ b/src/twoPhaseModels/compressibleTwoPhaseMixture/compressibleTwoPhaseMixture.H
@@ -119,6 +119,12 @@ public:
             return T_;
         }
 
+        //- Return mixture temperature [K]
+        const volScalarField& T() const
+        {
+            return T_;
+        }
+
         //- Return the thermo for phase 1
         const rhoThermo& thermo1() const
         {