[reoverstreet]

I'm working on a toy model of a simple tank with a band heater and don't see any heat transfer after meshing and running chtMultiRegionFoam in a 1 hour long simulation. The band heater should be providing 1000 watts of heat to the tank however, the temperature distribution on the tank/heater remains constant throughout the simulation. Right now I'm not sure what is wrong with the simulation configuration. The attached figure shows the simulation result at 3600 seconds. I am using OpenFoam version 2312.

Ultimately I would like to simulate water re-circulation from the inlet to the outlet in conjunction with heat transfer from the band heater.

The full case is attached in a zip file.

My work flow for meshing is:

surfaceFeatureExtract
blockMesh
decomposePar -force
mpirun -np 16 snappyHexMesh -parallel -overwrite
reconstructParMesh -latestTime -constant
splitMeshRegions -cellZones

system/snappyHexMeshDict:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      snappyHexMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

castellatedMesh true;
snap            true;
addLayers       false;

geometry
{
  tank.stl
      {
    type triSurfaceMesh;
    name tank;
      }
};

castellatedMeshControls
{
    maxLocalCells 1000;
    maxGlobalCells 100000;
    minRefinementCells 10;
    nCellsBetweenLevels 3;
    features
    (
     {
       file    "tank.eMesh";
       level 5;
     }
    );
    refinementSurfaces
    {
      tank
        {
      level (2 5);
        }
    }
    resolveFeatureAngle 30;
    refinementRegions
    {
    }
    locationsInMesh         
    (
     (( 50 75 130 ) water)
     (( 50 75 80 ) tank)
     (( 5 75 130 ) heater)
     );
    allowFreeStandingZoneFaces false;
}

snapControls
{
    nSmoothPatch 3;
    tolerance 1.0;
    nSolveIter 30;
    nRelaxIter 5;
    nFeatureSnapIter 10;
    implicitFeatureSnap true;
    explicitFeatureSnap true;
    multiRegionFeatureSnap false;
}

addLayersControls
{
    relativeSizes false;
    expansionRatio 0.0001;
    finalLayerThickness 0.0001;
    minThickness 0.0001;
    nGrow 0;
    featureAngle 10;
    nRelaxIter 5;
    nSmoothSurfaceNormals 1;
    nSmoothNormals 3;
    nSmoothThickness 10;
    maxFaceThicknessRatio 0.0005;
    maxThicknessToMedialRatio 0.0003;
    minMedianAxisAngle 90;
    nBufferCellsNoExtrude 0;
    nLayerIter 50;
    nRelaxedIter 20;
}

meshQualityControls
{
    #include "meshQualityDict"
    nSmoothScale 4;
    errorReduction 0.75;
}
writeFlags
(
    scalarLevels
    layerSets   
    layerFields    
);

mergeTolerance 1E-6;

// ************************************************************************* //

Boundary for the tank:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2312                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      binary;
    arch        "LSB;label=32;scalar=64";
    class       polyBoundaryMesh;
    location    "1/tank/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

3
(
    tank
    {
        type            mappedWall;
        inGroups        1(wall);
        nFaces          32450;
        startFace       444486;
        sampleMode      nearestPatchFace;
        sampleRegion    tank;
        samplePatch     tank_to_water;
    }
    tank_to_water
    {
        type            mappedWall;
        inGroups        1(wall);
        nFaces          38375;
        startFace       476936;
        sampleMode      nearestPatchFace;
        sampleRegion    tank;
        samplePatch     tank;
    }
    tank_to_heater
    {
        type            mappedWall;
        inGroups        1(wall);
        nFaces          6968;
        startFace       515311;
        sampleMode      nearestPatchFace;
        sampleRegion    heater;
        samplePatch     heater_to_tank;
    }
)

// ************************************************************************* //

Boundary for the heater:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2312                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      binary;
    arch        "LSB;label=32;scalar=64";
    class       polyBoundaryMesh;
    location    "1/heater/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

2
(
    heater
    {
        type            wall;
        inGroups        1(wall);
        nFaces          9886;
        startFace       102389;
    }
    heater_to_tank
    {
        type            mappedWall;
        inGroups        1(wall);
        nFaces          6968;
        startFace       112275;
        sampleMode      nearestPatchFace;
        sampleRegion    tank;
        samplePatch     tank_to_heater;
    }
)

// ************************************************************************* //

Initial conditions 0/tank/T

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2306                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 300;

boundaryField
{
    tank_to_water
    {
      type            compressible::turbulentTemperatureRadCoupledMixed;
      value           $internalField;
      Tnbr            T;
      kappaMethod     solidThermo;
      qrNbr           none;
      qr              none;
    }
    tank_to_heater
    {
      type            compressible::turbulentTemperatureRadCoupledMixed;
      value           $internalField;
      Tnbr            T;
      kappaMethod     solidThermo;
      qrNbr           none;
      qr              none;
    }
   tank
    {
      type            compressible::turbulentTemperatureRadCoupledMixed;
      value           $internalField;
      Tnbr            T;
      kappaMethod     solidThermo;
      qrNbr           none;
      qr              none;
    }

}

// ************************************************************************* //

Initial conditions 0/heater/T

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2306                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    arch        "LSB;label=32;scalar=64";
    class       volScalarField;
    location    "0/heater";
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [ 0 0 0 1 0 0 0 ];

internalField   uniform 338;

boundaryField
{
heater
  {
    type            externalWallHeatFluxTemperature;  
    mode            power;
    Q               uniform 1000;
    kappaMethod     solidThermo;
    value           $internalField;
  }
heater_to_tank
    {
      type            compressible::turbulentTemperatureRadCoupledMixed;
      qrNbr           none;
      qr              none;
      Tnbr            T;
      kappaMethod     solidThermo;
      useImplicit     false;
      value           $internalField;
    }

}


// ************************************************************************* //

fvSolution for both the tank and heater:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2306                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    h
    {
        solver           PCG;
        preconditioner   DIC;
        tolerance        1e-06;
        relTol           0.1;
    }

    hFinal
    {
        $h;
        tolerance        1e-06;
        relTol           0;
    }
}

PIMPLE
{
    nNonOrthogonalCorrectors 0;
}


// ************************************************************************* //