CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

HeatSource BC to the whole region in chtMultiRegionHeater

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By Zeppo

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 3, 2016, 07:21
Default HeatSource BC to the whole region in chtMultiRegionHeater
  #1
xsa
New Member
 
Xabi
Join Date: Sep 2016
Posts: 10
Rep Power: 10
xsa is on a distinguished road
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?
Attached Images
File Type: bmp heater.bmp (70.9 KB, 8 views)
xsa is offline   Reply With Quote

Old   November 4, 2016, 09:38
Default Solved but problems now with temperature instability
  #2
xsa
New Member
 
Xabi
Join Date: Sep 2016
Posts: 10
Rep Power: 10
xsa is on a distinguished road
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.
Attached Images
File Type: bmp heater.bmp (149.4 KB, 2 views)
xsa is offline   Reply With Quote

Old   November 4, 2016, 19:01
Default
  #3
Senior Member
 
Zeppo's Avatar
 
Sergei
Join Date: Dec 2009
Posts: 261
Rep Power: 21
Zeppo will become famous soon enough
fvOptions can help you with defining volumetric heat source
huadian001 likes this.
Zeppo is offline   Reply With Quote

Old   November 7, 2016, 06:07
Default
  #4
xsa
New Member
 
Xabi
Join Date: Sep 2016
Posts: 10
Rep Power: 10
xsa is on a distinguished road
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
Attached Images
File Type: bmp image1.bmp (94.6 KB, 9 views)
File Type: bmp image2.bmp (122.7 KB, 8 views)
xsa is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[mesh manipulation] Importing Multiple Meshes thomasnwalshiii OpenFOAM Meshing & Mesh Conversion 18 December 19, 2015 19:57
Unstabil Simulation with chtMultiRegionFoam mbay101 OpenFOAM Running, Solving & CFD 13 December 28, 2013 14:12
[Commercial meshers] Using starToFoam clo OpenFOAM Meshing & Mesh Conversion 33 September 26, 2012 05:04
[Other] StarToFoam error Kart OpenFOAM Meshing & Mesh Conversion 1 February 4, 2010 05:38
[Gmsh] Import gmsh msh to Foam adorean OpenFOAM Meshing & Mesh Conversion 24 April 27, 2005 09:19


All times are GMT -4. The time now is 20:11.