ENH: Improving settings for multiRegionHeater, adding example of contact resistance

and coded source fvOption in heater/fvOptions

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/bottomAir/fvSchemes

@@ -40,7 +40,7 @@ divSchemes
 
 laplacianSchemes
 {
-    default         Gauss linear limited corrected 0.333;
+    default         Gauss linear corrected;
 }
 
 interpolationSchemes
@@ -50,7 +50,7 @@ interpolationSchemes
 
 snGradSchemes
 {
-    default         limited corrected 0.333;
+    default         corrected;
 }
 
 fluxRequired

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/bottomAir/fvSolution

@@ -20,21 +20,21 @@ solvers
     {
         solver          PCG;
         preconditioner  DIC;
-        tolerance       1e-7;
+        tolerance       1e-6;
         relTol          0.1;
     }
 
     rhoFinal
     {
         $rho;
-        tolerance       1e-7;
+        tolerance       1e-6;
         relTol          0;
     }
 
     p_rgh
     {
         solver           GAMG;
-        tolerance        1e-7;
+        tolerance        1e-6;
         relTol           0.01;
 
         smoother         GaussSeidel;
@@ -43,12 +43,12 @@ solvers
         nCellsInCoarsestLevel 10;
         agglomerator     faceAreaPair;
         mergeLevels      1;
+        maxIter          100;
     }
 
     p_rghFinal
     {
         $p_rgh;
-        tolerance        1e-7;
         relTol           0;
     }
 
@@ -56,14 +56,14 @@ solvers
     {
         solver           PBiCG;
         preconditioner   DILU;
-        tolerance        1e-7;
+        tolerance        1e-6;
         relTol           0.1;
     }
 
     "(U|h|k|epsilon|R)Final"
     {
         $U;
-        tolerance        1e-7;
+        tolerance        1e-6;
         relTol           0;
     }
 }

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/controlDict

@@ -39,7 +39,7 @@ writeInterval   10;
 
 purgeWrite      0;
 
-writeFormat     binary;
+writeFormat     ascii;
 
 writePrecision  8;
 
@@ -51,7 +51,7 @@ timePrecision   6;
 
 runTimeModifiable yes;
 
-maxCo           0.3;
+maxCo           0.6;
 
 // Maximum diffusion number
 maxDi           10.0;

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/heater/changeDictionaryDict

@@ -51,6 +51,18 @@ dictionaryReplacement
                 kappaName       none;
                 value           uniform 300;
             }
+
+            heater_to_leftSolid
+            {
+                type            compressible::turbulentTemperatureCoupledBaffleMixed;
+                Tnbr            T;
+                kappa           solidThermo;
+                kappaName       none;
+                thicknessLayers (1e-3);
+                kappaLayers     (5e-4);
+                value           uniform 300;
+            }
+
             minY
             {
                 type            fixedValue;

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/heater/fvOptions

+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvOptions;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+energySource
+{
+    type            scalarCodedSource;
+
+    active          true;
+    selectionMode   all;
+
+    scalarCodedSourceCoeffs
+    {
+        fieldNames      (h);
+        redirectType    sourceTime;
+
+        codeInclude
+        #{
+
+        #};
+
+        codeCorrect
+        #{
+            //Pout<< "**codeCorrect**" << endl;
+        #};
+
+        codeAddSup
+        #{
+            const Time& time = mesh().time();
+            const scalarField& V = mesh_.V();
+            scalarField& heSource = eqn.source();
+            heSource -= 1e4*sqr(time.value())*V; // J/s/m3
+            Pout << "heat Source :" << -gSum(heSource) << endl;
+        #};
+
+        codeSetValue
+        #{
+           // Pout<< "**codeSetValue**" << endl;
+        #};
+
+        // Dummy entry. Make dependent on above to trigger recompilation
+        code
+        #{
+            $codeCorrect
+            $codeAddSup
+            $codeSetValue
+        #};
+    }
+
+    sourceTimeCoeffs
+    {
+        // Dummy entry
+    }
+ }
+
+
+
+
+
+// ************************************************************************* //

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/heater/fvSchemes

@@ -32,7 +32,7 @@ divSchemes
 laplacianSchemes
 {
     default             none;
-    laplacian(alpha,h)  Gauss linear limited corrected 0.333;
+    laplacian(alpha,h)  Gauss linear corrected;
 }
 
 interpolationSchemes
@@ -42,7 +42,7 @@ interpolationSchemes
 
 snGradSchemes
 {
-    default         limited corrected 0.333;
+    default         corrected;
 }
 
 fluxRequired

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/heater/fvSolution

@@ -34,7 +34,7 @@ solvers
 
 PIMPLE
 {
-    nNonOrthogonalCorrectors 1;
+    nNonOrthogonalCorrectors 0;
 }
 
 // ************************************************************************* //

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/leftSolid/changeDictionaryDict

@@ -47,6 +47,17 @@ dictionaryReplacement
                 kappaName       none;
                 value           uniform 300;
             }
+
+            leftSolid_to_heater
+            {
+                type            compressible::turbulentTemperatureCoupledBaffleMixed;
+                Tnbr            T;
+                kappa           solidThermo;
+                kappaName       none;
+                thicknessLayers (1e-3);
+                kappaLayers     (5e-4);
+                value           uniform 300;
+            }
         }
     }
 }

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/leftSolid/fvSchemes

@@ -32,7 +32,7 @@ divSchemes
 laplacianSchemes
 {
     default                 none;
-    laplacian(alpha,h)      Gauss linear limited corrected 0.333;
+    laplacian(alpha,h)      Gauss linear corrected;
 }
 
 interpolationSchemes
@@ -42,7 +42,7 @@ interpolationSchemes
 
 snGradSchemes
 {
-    default         limited corrected 0.333;
+    default         corrected;
 }
 
 fluxRequired

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/leftSolid/fvSolution

@@ -34,7 +34,7 @@ solvers
 
 PIMPLE
 {
-    nNonOrthogonalCorrectors 1;
+    nNonOrthogonalCorrectors 0;
 }
 
 // ************************************************************************* //

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/rightSolid/fvSchemes

@@ -32,7 +32,7 @@ divSchemes
 laplacianSchemes
 {
     default              none;
-    laplacian(alpha,h)   Gauss linear limited corrected 0.333;
+    laplacian(alpha,h)   Gauss linear corrected;
 }
 
 interpolationSchemes
@@ -42,7 +42,7 @@ interpolationSchemes
 
 snGradSchemes
 {
-    default         limited corrected 0.333;
+    default         corrected;
 }
 
 fluxRequired

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/rightSolid/fvSolution

@@ -34,7 +34,7 @@ solvers
 
 PIMPLE
 {
-    nNonOrthogonalCorrectors 1;
+    nNonOrthogonalCorrectors 0;
 }
 
 // ************************************************************************* //

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/topAir/fvSchemes

@@ -40,7 +40,7 @@ divSchemes
 
 laplacianSchemes
 {
-    default         Gauss linear limited corrected 0.333;
+    default         Gauss linear corrected;
 }
 
 interpolationSchemes
@@ -50,7 +50,7 @@ interpolationSchemes
 
 snGradSchemes
 {
-    default         limited corrected 0.333;
+    default         corrected;
 }
 
 fluxRequired

tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/topAir/fvSolution

@@ -20,21 +20,21 @@ solvers
     {
         solver          PCG;
         preconditioner  DIC;
-        tolerance       1e-7;
+        tolerance       1e-6;
         relTol          0.1;
     }
 
     rhoFinal
     {
         $rho;
-        tolerance       1e-7;
+        tolerance       1e-6;
         relTol          0;
     }
 
     p_rgh
     {
         solver           GAMG;
-        tolerance        1e-7;
+        tolerance        1e-6;
         relTol           0.01;
 
         smoother         GaussSeidel;
@@ -43,14 +43,13 @@ solvers
         nCellsInCoarsestLevel 10;
         agglomerator     faceAreaPair;
         mergeLevels      1;
-
         maxIter          100;
     }
 
     p_rghFinal
     {
         $p_rgh;
-        tolerance        1e-7;
+        tolerance        1e-6;
         relTol           0;
     }
 
@@ -58,14 +57,13 @@ solvers
     {
         solver           PBiCG;
         preconditioner   DILU;
-        tolerance        1e-7;
+        tolerance        1e-6;
         relTol           0.1;
     }
 
     "(U|h|k|epsilon|R)Final"
     {
         $U;
-        tolerance        1e-7;
         relTol           0;
     }
 }