From 9bed548227cf35517c82cbb81b7456e212682f8b Mon Sep 17 00:00:00 2001
From: sergio <s.ferraris@opencfd.co.uk>
Date: Thu, 14 Oct 2021 14:40:02 -0700
Subject: [PATCH] ENH: Adding coupled matrix option for solid-only cht cases

---
 .../include/createCoupledRegions.H             | 18 ++++++++++++++++++
 1 file changed, 18 insertions(+)

diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/include/createCoupledRegions.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/include/createCoupledRegions.H
index 8503be210a..2d62b08dd4 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/include/createCoupledRegions.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/include/createCoupledRegions.H
@@ -49,3 +49,21 @@ forAll(fluidRegions, i)
         break;
     }
 }
+
+forAll(solidRegions, i)
+{
+    if (coupled && fvMatrixAssemblyPtr.empty())
+    {
+        solidThermo& thermo = thermos[i];
+        Info << "Create fvMatrixAssembly." << endl;
+        fvMatrixAssemblyPtr.reset
+        (
+            new fvMatrix<scalar>
+            (
+                thermo.he(),
+                dimEnergy/dimTime
+            )
+        );
+        break;
+    }
+}