
[Sponsors] 
HeatSource BC to the whole region in chtMultiRegionHeater 

LinkBack  Thread Tools  Search this Thread  Display Modes 
November 3, 2016, 06:21 
HeatSource BC to the whole region in chtMultiRegionHeater

#1 
New Member
Xabi
Join Date: Sep 2016
Posts: 10
Rep Power: 6 
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; } 

November 4, 2016, 08:38 
Solved but problems now with temperature instability

#2 
New Member
Xabi
Join Date: Sep 2016
Posts: 10
Rep Power: 6 
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; } heater Code:
{ type patch; inGroups 1(wall); nFaces 608; startFace 1648; } heater_heater Code:
{ type zeroGradient; } heater_heater_slave { type zeroGradient; } 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; } 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.01709e005, 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.486511e005, 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.980892e005, No Iterations 3 GAMG: Solving for p_rgh, Initial residual = 0.0003848775, Final residual = 3.67114e005, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 5.09555e005, Final residual = 5.09555e005, 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.329981e005, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 3.627199e005, Final residual = 3.627199e005, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 3.627199e005, Final residual = 3.627199e005, 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.037174e005, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 3.1473e005, Final residual = 3.1473e005, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 3.1473e005, Final residual = 3.1473e005, 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.446879e005, 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.750125e005, 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.044495e005, No Iterations 3 GAMG: Solving for p_rgh, Initial residual = 0.0003264729, Final residual = 3.640505e005, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.76597e005, Final residual = 4.76597e005, 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.328459e005, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 3.673964e005, Final residual = 3.673964e005, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 3.673964e005, Final residual = 3.673964e005, 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.342866e005, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 1.898267e005, Final residual = 1.898267e005, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 1.898267e005, Final residual = 1.898267e005, 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.695403e005, No Iterations 7 DICPCG: Solving for h, Initial residual = 0.01523796, Final residual = 2.769823e005, No Iterations 4 DICPCG: Solving for h, Initial residual = 0.0003536687, Final residual = 5.470599e005, 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.178876e005, No Iterations 6 DICPCG: Solving for h, Initial residual = 0.01340643, Final residual = 3.755647e005, No Iterations 3 DICPCG: Solving for h, Initial residual = 0.0001187311, Final residual = 1.405485e005, 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.471222e005, No Iterations 7 DICPCG: Solving for h, Initial residual = 0.01153215, Final residual = 3.378124e005, No Iterations 3 DICPCG: Solving for h, Initial residual = 0.0001115365, Final residual = 1.417605e005, 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.385133e009 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.52002e005, No Iterations 2 Min/max T:4.170379e+007 4.004755e+007 GAMG: Solving for p_rgh, Initial residual = 0.9937021, Final residual = 5.525147e018, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.745797e009, Final residual = 4.745797e009, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.745797e009, Final residual = 4.745797e009, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (bottomAir): sum local = 1.663656e008, global = 7.308624e018, cumulative = 0.001561362 GAMG: Solving for p_rgh, Initial residual = 0.0002569213, Final residual = 4.207221e014, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 8.978004e010, Final residual = 8.978004e010, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 8.978004e010, Final residual = 8.978004e010, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (bottomAir): sum local = 5.050287e009, global = 1.075352e013, cumulative = 0.001561362 GAMG: Solving for p_rgh, Initial residual = 2.600384e007, Final residual = 2.600384e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 2.600384e007, Final residual = 2.600384e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 2.600384e007, Final residual = 2.600384e007, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (bottomAir): sum local = 1.462762e006, global = 1.07533e013, cumulative = 0.001561362 GAMG: Solving for p_rgh, Initial residual = 2.600384e007, Final residual = 2.600384e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 2.600384e007, Final residual = 2.600384e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 2.600384e007, Final residual = 2.600384e007, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (bottomAir): sum local = 1.462762e006, global = 1.075335e013, 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.881394e005, No Iterations 2 Min/max T:3.161487e+007 3.228481e+007 GAMG: Solving for p_rgh, Initial residual = 0.9917856, Final residual = 4.976222e018, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.345829e009, Final residual = 4.345829e009, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.345829e009, Final residual = 4.345829e009, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (topAir): sum local = 7.226967e009, global = 2.21906e019, cumulative = 0.0008474586 GAMG: Solving for p_rgh, Initial residual = 0.0004717654, Final residual = 4.723399e014, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 1.770327e009, Final residual = 1.770327e009, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 1.770327e009, Final residual = 1.770327e009, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (topAir): sum local = 5.271096e009, global = 4.341831e014, cumulative = 0.0008474586 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (topAir): sum local = 1.369762e006, global = 4.342164e014, cumulative = 0.0008474586 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (topAir): sum local = 1.369762e006, global = 4.342164e014, cumulative = 0.0008474586 Solving for solid region heater DICPCG: Solving for h, Initial residual = 4.345812e009, Final residual = 4.345812e009, No Iterations 0 DICPCG: Solving for h, Initial residual = 4.352342e009, Final residual = 4.352342e009, No Iterations 0 DICPCG: Solving for h, Initial residual = 4.352342e009, Final residual = 4.352342e009, No Iterations 0 Min/max T:174.5245 318.7215 Solving for solid region leftSolid DICPCG: Solving for h, Initial residual = 5.845175e008, Final residual = 5.845175e008, No Iterations 0 DICPCG: Solving for h, Initial residual = 5.849886e008, Final residual = 5.849886e008, No Iterations 0 DICPCG: Solving for h, Initial residual = 5.849886e008, Final residual = 5.849886e008, No Iterations 0 Min/max T:215.0183 888.1072 Solving for solid region rightSolid DICPCG: Solving for h, Initial residual = 2.489504e008, Final residual = 2.489504e008, No Iterations 0 DICPCG: Solving for h, Initial residual = 2.488219e008, Final residual = 2.488219e008, No Iterations 0 DICPCG: Solving for h, Initial residual = 2.488219e008, Final residual = 2.488219e008, 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.323624e011 max: 1.656567e011 Region: leftSolid Diffusion Number mean: 1.32393e011 max: 1.640852e011 Region: rightSolid Diffusion Number mean: 1.323903e011 max: 1.63954e011 deltaT = 2.862159e009 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.052903e006, 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/OpenFOAM3.0.x/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66. FOAM aborting I can supply the case folder if needed. 

November 4, 2016, 18:01 

#3 
Senior Member
Sergei
Join Date: Dec 2009
Posts: 253
Rep Power: 17 
fvOptions can help you with defining volumetric heat source


November 7, 2016, 05:07 

#4 
New Member
Xabi
Join Date: Sep 2016
Posts: 10
Rep Power: 6 
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 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 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; } 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; } Code:
… oldInternalFaces { type empty; } heater { type zeroGradient; } heater_slave { type zeroGradient; } 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.881394e005, No Iterations 2 Min/max T:3.161487e+007 3.228481e+007 GAMG: Solving for p_rgh, Initial residual = 0.9917856, Final residual = 4.976222e018, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.345829e009, Final residual = 4.345829e009, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.345829e009, Final residual = 4.345829e009, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (topAir): sum local = 7.226967e009, global = 2.21906e019, cumulative = 0.0008474586 GAMG: Solving for p_rgh, Initial residual = 0.0004717654, Final residual = 4.723399e014, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 1.770327e009, Final residual = 1.770327e009, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 1.770327e009, Final residual = 1.770327e009, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (topAir): sum local = 5.271096e009, global = 4.341831e014, cumulative = 0.0008474586 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (topAir): sum local = 1.369762e006, global = 4.342164e014, cumulative = 0.0008474586 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 GAMG: Solving for p_rgh, Initial residual = 4.600421e007, Final residual = 4.600421e007, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (topAir): sum local = 1.369762e006, global = 4.342164e014, cumulative = 0.0008474586 Solving for solid region heater DICPCG: Solving for h, Initial residual = 4.345812e009, Final residual = 4.345812e009, No Iterations 0 DICPCG: Solving for h, Initial residual = 4.352342e009, Final residual = 4.352342e009, No Iterations 0 DICPCG: Solving for h, Initial residual = 4.352342e009, Final residual = 4.352342e009, No Iterations 0 Min/max T:174.5245 318.7215 Solving for solid region leftSolid DICPCG: Solving for h, Initial residual = 5.845175e008, Final residual = 5.845175e008, No Iterations 0 DICPCG: Solving for h, Initial residual = 5.849886e008, Final residual = 5.849886e008, No Iterations 0 DICPCG: Solving for h, Initial residual = 5.849886e008, Final residual = 5.849886e008, No Iterations 0 Min/max T:215.0183 888.1072 Solving for solid region rightSolid DICPCG: Solving for h, Initial residual = 2.489504e008, Final residual = 2.489504e008, No Iterations 0 DICPCG: Solving for h, Initial residual = 2.488219e008, Final residual = 2.488219e008, No Iterations 0 DICPCG: Solving for h, Initial residual = 2.488219e008, Final residual = 2.488219e008, 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.323624e011 max: 1.656567e011 Region: leftSolid Diffusion Number mean: 1.32393e011 max: 1.640852e011 Region: rightSolid Diffusion Number mean: 1.323903e011 max: 1.63954e011 deltaT = 2.862159e009 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.052903e006, 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/OpenFOAM3.0.x/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66. FOAM aborting So does somebody have any idea of where the problem could be? Regards 

Thread Tools  Search this Thread 
Display Modes  


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