diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousSolid/setPorousRegionSolidFields.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousSolid/setPorousRegionSolidFields.H
index beef8af1f5..102ba89a09 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousSolid/setPorousRegionSolidFields.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousSolid/setPorousRegionSolidFields.H
@@ -9,7 +9,7 @@
     tmp<volScalarField> tcp = thermo.Cp();
     const volScalarField& cp = tcp();
 
-    tmp<volScalarField> tkappa = betav*thermo.K();
+    tmp<volScalarField> tkappa = thermo.K();
     //tmp<volSymmTensorField> tkappa = thermo.directionalK()*betav;
 
     const volScalarField& kappa = tkappa();
diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousSolid/solvePorousSolid.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousSolid/solvePorousSolid.H
index 75d8d1347d..b4346275e4 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousSolid/solvePorousSolid.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousSolid/solvePorousSolid.H
@@ -3,7 +3,7 @@
     {
         tmp<fvScalarMatrix> TEqn
         (
-          - fvm::laplacian(kappa, T, "laplacian(K,T)")
+          - fvm::laplacian(betav*kappa, T, "laplacian(K,T)")
           + sources(rho, T)
         );
 
diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/porousSolid/setPorousRegionSolidFields.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/porousSolid/setPorousRegionSolidFields.H
index 97df0171dc..f4c7081e9b 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/porousSolid/setPorousRegionSolidFields.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/porousSolid/setPorousRegionSolidFields.H
@@ -9,7 +9,7 @@
     tmp<volScalarField> tcp = thermo.Cp();
     const volScalarField& cp = tcp();
 
-    tmp<volScalarField> tkappa = betav*thermo.K();
+    tmp<volScalarField> tkappa = thermo.K();
     //tmp<volSymmTensorField> tkappa = thermo.directionalK()*betav;
 
     const volScalarField& kappa = tkappa();
diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/porousSolid/solvePorousSolid.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/porousSolid/solvePorousSolid.H
index f476df479d..e8941cf2ad 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/porousSolid/solvePorousSolid.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/porousSolid/solvePorousSolid.H
@@ -9,7 +9,7 @@ if (finalIter)
         tmp<fvScalarMatrix> TEqn
         (
             fvm::ddt(betav*rho*cp, T)
-          - fvm::laplacian(kappa, T, "laplacian(K,T)")
+          - fvm::laplacian(betav*kappa, T, "laplacian(K,T)")
           + sources(rho, T)
         );
 
diff --git a/src/fieldSources/basicSource/interRegionHeatTransferModel/interRegionHeatTransferModel/interRegionHeatTransferModel.C b/src/fieldSources/basicSource/interRegionHeatTransferModel/interRegionHeatTransferModel/interRegionHeatTransferModel.C
index 497ff86ef5..24c0f612dd 100644
--- a/src/fieldSources/basicSource/interRegionHeatTransferModel/interRegionHeatTransferModel/interRegionHeatTransferModel.C
+++ b/src/fieldSources/basicSource/interRegionHeatTransferModel/interRegionHeatTransferModel/interRegionHeatTransferModel.C
@@ -167,16 +167,13 @@ void Foam::interRegionHeatTransferModel::addSup
     const fvMesh& secondaryMesh =
         mesh_.time().lookupObject<fvMesh>(mapRegionName_);
 
-    basicThermo&  secondaryThermo =
-        const_cast<basicThermo&>
-        (
-            secondaryMesh.lookupObject<basicThermo>("thermophysicalProperties")
-        );
+    const volScalarField& Tsecondary =
+        secondaryMesh.lookupObject<volScalarField>("T");
 
     secondaryToPrimaryInterpPtr_->interpolateInternalField
     (
         Tmapped,
-        secondaryThermo.T(),
+        Tsecondary,
         meshToMesh::MAP,
         eqOp<scalar>()
     );