CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   HeatSource BC to the whole region in chtMultiRegionHeater (https://www.cfd-online.com/Forums/openfoam-solving/179600-heatsource-bc-whole-region-chtmultiregionheater.html)

xsa November 3, 2016 07:21

HeatSource BC to the whole region in chtMultiRegionHeater
 
1 Attachment(s)
Hello,

I am having some troubles in implementing BC to a patch in chtMultiRegionFoam.
I have taken the tutorial "snappyMultiRegionHeater" and implemented in the region “heater” in the boundary “minY” a heat source of 500 W with the command “turbulentHeatFluxTemperature” instead of having the fixed 500 K temperature. So far the simulation goes well (image attached).


But my intention is to set that BC to the whole heater. As it is set now, the heat source is just set in the “minT” boundary but I want the complete heater to be the heat source. The problem is that with “splitMeshRegions”, I can set BCs just to the interface between different regions, but not to the whole region. I have tried using “topoSetDict” and “createPatch” without luck. So my question is:
Is it possible to fix the region “heater” as a boundary, maintaining the interface boundaries? For example having this:


Code:

heater_to_bottomAir
      {
          type            compressible::turbulentTemperatureCoupledBaffleMixed;
          value          uniform 300;
          Tnbr            T;
          kappa          solidThermo;
          kappaName      none;
      }
      heater
      {
          type            compressible::turbulentHeatFluxTemperature;
          value          300;
          heatSource      power;
          q              uniform 500;
          kappa          solidThermo;
          kappaName      none;
          Qr              none;
      }

I have tried just adding this but OpenFoam does not recognize “heater”. So does somebody know how to set "heater" for BC as well?

xsa November 4, 2016 09:38

Solved but problems now with temperature instability
 
1 Attachment(s)
Hi all,
After a lot of trial and errors i have managed to set a BC to the whole heater, even though the process is not very elegant.
After creating the mesh and splitting it, I had to edit different files manually.
I entered the boundary file in polyMesh and add these two lines in the end plus increment with two the boundaries.


Code:


heater_heater
      {
          type            wall;
          inGroups        1(wall);
          nFaces          0;
          startFace      49094;
      }
      heater_heater_slave
      {
          type            wall;
          inGroups        1(wall);
          nFaces          0;
          startFace      49094;
      }

After that the same in the boundary file in heater/polyMesh:
heater
Code:


  {
          type            patch;
          inGroups        1(wall);
          nFaces          608;
          startFace      1648;
 
      }

The edit the cellToRegion file in 0 folder and add
heater_heater
Code:


 {
          type            zeroGradient;
      }
      heater_heater_slave
      {
          type            zeroGradient;
      }

And finally add the heater BC in all files inside 0/heater.
I edited the changeDirectoryDict file to insert the heatsource:
Code:


  heater
              {
          type            compressible::turbulentHeatFluxTemperature;
                  heatSource        power;
                  q                            uniform 500;
                  kappa                  solidThermo;
                  kappaName      none;
                  Qr                          none;
                  value                    300;
 
              }

So I got the heat source in the whole heater (see attached image)

But now, I wanted to try the same but making the blockMesh bigger so that I have to erase the boundaries and isolate the mesh, cause later I will have to do that for simulating a motor. So i made the mesh double and double cells so that the final mesh is the same for the desired region. Then I used batch.setSet and “subsetMesh isolate” to isolate the mesh. After creating the mesh I had to do some editing but similar to the previous simulation. But when I start simulating, the temperature starts getting unrealistic values and then it crashes, even though I have not changed any boundary condition.
Code:

Solving for fluid region bottomAir
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 5.01709e-005, No Iterations 7
Min/max T:291.332 360.629
GAMG:  Solving for p_rgh, Initial residual = 0.9999684, Final residual = 0.004615176, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.05179337, Final residual = 0.0002363303, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.003098775, Final residual = 7.486511e-005, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (bottomAir): sum local = 329.3255, global = 0.0003897242, cumulative = 0.0003897242
GAMG:  Solving for p_rgh, Initial residual = 0.005565491, Final residual = 5.980892e-005, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.0003848775, Final residual = 3.67114e-005, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 5.09555e-005, Final residual = 5.09555e-005, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (bottomAir): sum local = 232.87, global = 0.0003904275, cumulative = 0.0007801518
GAMG:  Solving for p_rgh, Initial residual = 0.0002054582, Final residual = 2.329981e-005, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.627199e-005, Final residual = 3.627199e-005, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 3.627199e-005, Final residual = 3.627199e-005, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (bottomAir): sum local = 165.7727, global = 0.0003905498, cumulative = 0.001170702
GAMG:  Solving for p_rgh, Initial residual = 0.0001924865, Final residual = 2.037174e-005, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.1473e-005, Final residual = 3.1473e-005, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 3.1473e-005, Final residual = 3.1473e-005, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (bottomAir): sum local = 143.8902, global = 0.0003906605, cumulative = 0.001561362

Solving for fluid region topAir
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 9.446879e-005, No Iterations 6
Min/max T:290.8782 359.4351
GAMG:  Solving for p_rgh, Initial residual = 0.9999777, Final residual = 0.003114044, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.04809523, Final residual = 0.0001513841, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.00286152, Final residual = 4.750125e-005, No Iterations 2
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (topAir): sum local = 227.8107, global = 0.0002112262, cumulative = 0.0002112262
GAMG:  Solving for p_rgh, Initial residual = 0.005279227, Final residual = 6.044495e-005, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.0003264729, Final residual = 3.640505e-005, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 4.76597e-005, Final residual = 4.76597e-005, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (topAir): sum local = 238.8201, global = 0.0002117873, cumulative = 0.0004230135
GAMG:  Solving for p_rgh, Initial residual = 0.0001569056, Final residual = 2.328459e-005, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 3.673964e-005, Final residual = 3.673964e-005, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 3.673964e-005, Final residual = 3.673964e-005, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (topAir): sum local = 184.0833, global = 0.0002121023, cumulative = 0.0006351158
GAMG:  Solving for p_rgh, Initial residual = 0.0001085723, Final residual = 1.342866e-005, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.898267e-005, Final residual = 1.898267e-005, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 1.898267e-005, Final residual = 1.898267e-005, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (topAir): sum local = 95.12819, global = 0.0002123428, cumulative = 0.0008474586

Solving for solid region heater
DICPCG:  Solving for h, Initial residual = 1, Final residual = 3.695403e-005, No Iterations 7
DICPCG:  Solving for h, Initial residual = 0.01523796, Final residual = 2.769823e-005, No Iterations 4
DICPCG:  Solving for h, Initial residual = 0.0003536687, Final residual = 5.470599e-005, No Iterations 1
Min/max T:286.3776 318.7215

Solving for solid region leftSolid
DICPCG:  Solving for h, Initial residual = 1, Final residual = 8.178876e-005, No Iterations 6
DICPCG:  Solving for h, Initial residual = 0.01340643, Final residual = 3.755647e-005, No Iterations 3
DICPCG:  Solving for h, Initial residual = 0.0001187311, Final residual = 1.405485e-005, No Iterations 1
Min/max T:294.1181 325.7575

Solving for solid region rightSolid
DICPCG:  Solving for h, Initial residual = 1, Final residual = 2.471222e-005, No Iterations 7
DICPCG:  Solving for h, Initial residual = 0.01153215, Final residual = 3.378124e-005, No Iterations 3
DICPCG:  Solving for h, Initial residual = 0.0001115365, Final residual = 1.417605e-005, No Iterations 1
Min/max T:293.677 328.5112
ExecutionTime = 0.905 s  ClockTime = 1 s

Region: bottomAir Courant Number mean: 2461431 max: 1.04816e+008
Region: topAir Courant Number mean: 2125396 max: 9.421526e+007
Region: heater Diffusion Number mean: 0.004624564 max: 0.005787824
Region: leftSolid Diffusion Number mean: 0.004625633 max: 0.005732918
Region: rightSolid Diffusion Number mean: 0.00462554 max: 0.005728334
deltaT = 2.385133e-009
Time = 0.833333


Solving for fluid region bottomAir
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 0.9999842, Final residual = 2.52002e-005, No Iterations 2
Min/max T:-4.170379e+007 4.004755e+007
GAMG:  Solving for p_rgh, Initial residual = 0.9937021, Final residual = 5.525147e-018, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 4.745797e-009, Final residual = 4.745797e-009, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 4.745797e-009, Final residual = 4.745797e-009, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (bottomAir): sum local = 1.663656e-008, global = -7.308624e-018, cumulative = 0.001561362
GAMG:  Solving for p_rgh, Initial residual = 0.0002569213, Final residual = 4.207221e-014, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 8.978004e-010, Final residual = 8.978004e-010, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 8.978004e-010, Final residual = 8.978004e-010, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (bottomAir): sum local = 5.050287e-009, global = -1.075352e-013, cumulative = 0.001561362
GAMG:  Solving for p_rgh, Initial residual = 2.600384e-007, Final residual = 2.600384e-007, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 2.600384e-007, Final residual = 2.600384e-007, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 2.600384e-007, Final residual = 2.600384e-007, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (bottomAir): sum local = 1.462762e-006, global = -1.07533e-013, cumulative = 0.001561362
GAMG:  Solving for p_rgh, Initial residual = 2.600384e-007, Final residual = 2.600384e-007, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 2.600384e-007, Final residual = 2.600384e-007, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 2.600384e-007, Final residual = 2.600384e-007, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (bottomAir): sum local = 1.462762e-006, global = -1.075335e-013, cumulative = 0.001561362

Solving for fluid region topAir
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 0.9999844, Final residual = 5.881394e-005, No Iterations 2
Min/max T:-3.161487e+007 3.228481e+007
GAMG:  Solving for p_rgh, Initial residual = 0.9917856, Final residual = 4.976222e-018, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 4.345829e-009, Final residual = 4.345829e-009, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 4.345829e-009, Final residual = 4.345829e-009, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (topAir): sum local = 7.226967e-009, global = -2.21906e-019, cumulative = 0.0008474586
GAMG:  Solving for p_rgh, Initial residual = 0.0004717654, Final residual = 4.723399e-014, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 1.770327e-009, Final residual = 1.770327e-009, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 1.770327e-009, Final residual = 1.770327e-009, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (topAir): sum local = 5.271096e-009, global = -4.341831e-014, cumulative = 0.0008474586
GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (topAir): sum local = 1.369762e-006, global = -4.342164e-014, cumulative = 0.0008474586
GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (topAir): sum local = 1.369762e-006, global = -4.342164e-014, cumulative = 0.0008474586

Solving for solid region heater
DICPCG:  Solving for h, Initial residual = 4.345812e-009, Final residual = 4.345812e-009, No Iterations 0
DICPCG:  Solving for h, Initial residual = 4.352342e-009, Final residual = 4.352342e-009, No Iterations 0
DICPCG:  Solving for h, Initial residual = 4.352342e-009, Final residual = 4.352342e-009, No Iterations 0
Min/max T:174.5245 318.7215

Solving for solid region leftSolid
DICPCG:  Solving for h, Initial residual = 5.845175e-008, Final residual = 5.845175e-008, No Iterations 0
DICPCG:  Solving for h, Initial residual = 5.849886e-008, Final residual = 5.849886e-008, No Iterations 0
DICPCG:  Solving for h, Initial residual = 5.849886e-008, Final residual = 5.849886e-008, No Iterations 0
Min/max T:-215.0183 888.1072

Solving for solid region rightSolid
DICPCG:  Solving for h, Initial residual = 2.489504e-008, Final residual = 2.489504e-008, No Iterations 0
DICPCG:  Solving for h, Initial residual = 2.488219e-008, Final residual = 2.488219e-008, No Iterations 0
DICPCG:  Solving for h, Initial residual = 2.488219e-008, Final residual = 2.488219e-008, No Iterations 0
Min/max T:-239.0443 767.7877
ExecutionTime = 1.185 s  ClockTime = 1 s

Region: bottomAir Courant Number mean: 0.0007467967 max: 0.02241564
Region: topAir Courant Number mean: 0.000683475 max: 0.02254898
Region: heater Diffusion Number mean: 1.323624e-011 max: 1.656567e-011
Region: leftSolid Diffusion Number mean: 1.32393e-011 max: 1.640852e-011
Region: rightSolid Diffusion Number mean: 1.323903e-011 max: 1.63954e-011
deltaT = 2.862159e-009
Time = 0.833333


Solving for fluid region bottomAir
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 0.1387004, Final residual = 6.052903e-006, No Iterations 1


--> FOAM FATAL ERROR:
Maximum number of iterations exceeded

    From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>]
    in file /opt/OpenFOAM/OpenFOAM-3.0.x/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.

FOAM aborting

Does somebody know where could be the problem?
I can supply the case folder if needed.

Zeppo November 4, 2016 19:01

fvOptions can help you with defining volumetric heat source

xsa November 7, 2016 06:07

2 Attachment(s)
Thanks for answering, but the issue of setting the heat source is solved, i generate it with ""compressible::turbulentHeatFluxTemperature" where i fix 500 W. The problem was that I did not have the patch "heater", inside the region heater. The only patches where the intersections between the heater and the adjacent regions (heater_to_topAir, heater_to_leftSolid, etc). But I solved it by inserting the patch manually in the “boundary” file.
My problem now is concerned with the mesh I guess, cause I have been able to simulate it in the “snappyMultiRegionHeater” tutorial. Now, I want to create a bigger mesh with blockmesh and then eliminate the outside regions (as I will have to do that for my project). So I created the box double the size as seen in the image1.


But meshing this way, a new domain called “domain0” appears as the outside boundary. I searched in the forums and the way to eliminate this is by using the commands setSet and subsetMesh in “allrun”:
Code:

  runApplication blockMesh
  runApplication surfaceFeatureExtract
  runApplication snappyHexMesh -overwrite
  runApplication setSet -batch batch.setSet
  runApplication subsetMesh -overwrite isolation
  runApplication splitMeshRegions -cellZones –overwrite

Where the file “batch.setSet” contains:

Code:

  cellSet isolation new zoneToCell bottomAir
  cellSet isolation add zoneToCell heater
  cellSet isolation add zoneToCell leftSolid
  cellSet isolation add zoneToCell rightSolid
  cellSet isolation add zoneToCell topAir
  cellSet isolation subset

Si if I run this way, my mesh looks as seen in image2, without the outside box


A new folder appears in the case directory called “VTK” and a new patch in the “boundary file” called “oldInternalFaces”.
Then I follow the steps I did with the successful simulation.
I add two patches for "heater" in “boundary” file and increase the number +2
Code:

    …
  heater_heater_to_bottomAir_slave
      {
          type            wall;
          inGroups        1(wall);
          nFaces          0;
          startFace      49094;
      }
      oldInternalFaces
      {
          type            empty;
          inGroups        1(empty);
          nFaces          5600;
          startFace      49094;
      }
      heater_heater
      {
          type            wall;
          inGroups        1(wall);
          nFaces          0;
          startFace      49094;
      }
      heater_heater_slave
      {
          type            wall;
          inGroups        1(wall);
          nFaces          0;
          startFace      49094;
      }

Then add a patch "heater" to the heater/polyMesh/boundary file and increase +1

Code:

  …
      heater_to_bottomAir
      {
          type            mappedWall;
          inGroups        1(wall);
          nFaces          368;
          startFace      1888;
          sampleMode      nearestPatchFace;
          sampleRegion    bottomAir;
          samplePatch    bottomAir_to_heater;
      }
      heater
      {
          type            patch;
          inGroups        1(wall);
          nFaces          608;
          startFace      1648;
      }

Then add the patches in 0/cellToRegion file
Code:

  …   
  oldInternalFaces
      {
          type            empty;
      }
      heater
      {
          type            zeroGradient;
      }
      heater_slave
      {
          type            zeroGradient;
      }

And then add the patch to all files inside 0/heater, which are cellToRegion, p, rho, T.

After that I use the command changeDictionary and everything seems fine, as the only thing I have done is erase the BCs of maxX, minY, etc. and add empty BC to “oldInternalFaces” in the required regions. But when I simulate, the next lines appear.
Code:

  Solving for fluid region topAir
  diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
  DILUPBiCG:  Solving for h, Initial residual = 0.9999844, Final residual = 5.881394e-005, No Iterations 2
  Min/max T:-3.161487e+007 3.228481e+007
  GAMG:  Solving for p_rgh, Initial residual = 0.9917856, Final residual = 4.976222e-018, No Iterations 1
  GAMG:  Solving for p_rgh, Initial residual = 4.345829e-009, Final residual = 4.345829e-009, No Iterations 0
  GAMG:  Solving for p_rgh, Initial residual = 4.345829e-009, Final residual = 4.345829e-009, No Iterations 0
  diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
  time step continuity errors (topAir): sum local = 7.226967e-009, global = -2.21906e-019, cumulative = 0.0008474586
  GAMG:  Solving for p_rgh, Initial residual = 0.0004717654, Final residual = 4.723399e-014, No Iterations 1
  GAMG:  Solving for p_rgh, Initial residual = 1.770327e-009, Final residual = 1.770327e-009, No Iterations 0
  GAMG:  Solving for p_rgh, Initial residual = 1.770327e-009, Final residual = 1.770327e-009, No Iterations 0
  diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
  time step continuity errors (topAir): sum local = 5.271096e-009, global = -4.341831e-014, cumulative = 0.0008474586
  GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
  GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
  GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
  diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
  time step continuity errors (topAir): sum local = 1.369762e-006, global = -4.342164e-014, cumulative = 0.0008474586
  GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
  GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
  GAMG:  Solving for p_rgh, Initial residual = 4.600421e-007, Final residual = 4.600421e-007, No Iterations 0
  diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
  time step continuity errors (topAir): sum local = 1.369762e-006, global = -4.342164e-014, cumulative = 0.0008474586
 
  Solving for solid region heater
  DICPCG:  Solving for h, Initial residual = 4.345812e-009, Final residual = 4.345812e-009, No Iterations 0
  DICPCG:  Solving for h, Initial residual = 4.352342e-009, Final residual = 4.352342e-009, No Iterations 0
  DICPCG:  Solving for h, Initial residual = 4.352342e-009, Final residual = 4.352342e-009, No Iterations 0
  Min/max T:174.5245 318.7215
 
  Solving for solid region leftSolid
  DICPCG:  Solving for h, Initial residual = 5.845175e-008, Final residual = 5.845175e-008, No Iterations 0
  DICPCG:  Solving for h, Initial residual = 5.849886e-008, Final residual = 5.849886e-008, No Iterations 0
  DICPCG:  Solving for h, Initial residual = 5.849886e-008, Final residual = 5.849886e-008, No Iterations 0
  Min/max T:-215.0183 888.1072
 
  Solving for solid region rightSolid
  DICPCG:  Solving for h, Initial residual = 2.489504e-008, Final residual = 2.489504e-008, No Iterations 0
  DICPCG:  Solving for h, Initial residual = 2.488219e-008, Final residual = 2.488219e-008, No Iterations 0
  DICPCG:  Solving for h, Initial residual = 2.488219e-008, Final residual = 2.488219e-008, No Iterations 0
  Min/max T:-239.0443 767.7877
  ExecutionTime = 1.17 s  ClockTime = 1 s
 
  Region: bottomAir Courant Number mean: 0.0007467967 max: 0.02241564
  Region: topAir Courant Number mean: 0.000683475 max: 0.02254898
  Region: heater Diffusion Number mean: 1.323624e-011 max: 1.656567e-011
  Region: leftSolid Diffusion Number mean: 1.32393e-011 max: 1.640852e-011
  Region: rightSolid Diffusion Number mean: 1.323903e-011 max: 1.63954e-011
  deltaT = 2.862159e-009
  Time = 0.833333
 
 
  Solving for fluid region bottomAir
  diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
  DILUPBiCG:  Solving for h, Initial residual = 0.1387004, Final residual = 6.052903e-006, No Iterations 1
 
 
  --> FOAM FATAL ERROR:
  Maximum number of iterations exceeded
 
      From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>]
      in file /opt/OpenFOAM/OpenFOAM-3.0.x/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.
 
  FOAM aborting

As you can see, the temperature gets unrealistic values.I hope I have explained well. My project is to simulate the cooling of an induction motor where mesh is quite complex and I will have to define a bigger box for the geometry as here, so it is important to me to fix this issue.
So does somebody have any idea of where the problem could be?
Regards


All times are GMT -4. The time now is 14:43.