Foam::divide, transient conjugate heat transfer
 

Hello foamers, I am trying to simulate transient conjugate heat transfer inside a cylinder and to the air. When I type chtMultiRegionFoam on the command window I get this error. Dictionaries are uploaded. Any help will be appreciated. Thank you very much.

--> FOAM Warning :
From function virtual void Foam::probes::findElements(const Foam::fvMesh&)
in file probes/probes.C at line 122
Did not find location (0.7 0.8 0.25) in any cell. Skipping location.
Region: fluid Courant Number mean: 5.294362 max: 509.8198
Region: solid Diffusion Number mean: 401.0293 max: 16244.76
Time = 1


Solving for fluid region fluid

Program received signal SIGFPE, Arithmetic exception.
0x00007ffff24ca208 in Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) () from /home/user2/OpenFOAM-v1706/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
(gdb) backtrace
#0 0x00007ffff24ca208 in Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) () from /home/user2/OpenFOAM-v1706/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#1 0x00007ffff24cbf64 in Foam::operator/(Foam::UList<double> const&, Foam::UList<double> const&) () from /home/user2/OpenFOAM-v1706/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 0x00007ffff2318e84 in Foam::diagonalSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const () from /home/user2/OpenFOAM-v1706/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#3 0x00007ffff480e41e in Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) () from /home/user2/OpenFOAM-v1706/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#4 0x00007ffff434b467 in Foam::fvMatrix<double>::solveSegregatedOrCoupled(F oam::dictionary const&) () from /home/user2/OpenFOAM-v1706/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#5 0x00007ffff43116af in Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const () from /home/user2/OpenFOAM-v1706/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#6 0x000000000046ebdf in Foam::fvMatrix<double>::solve() ()
#7 0x000000000043f716 in main ()

Hello Cagatayemre!

Where is your heat source?

Why you are using chtMultiRegionFoam? There is just one region in the example you uploaded!

Usually you need at least one region as a fluid and one as solid to use chtMultiRegion(Simple)Foam.

Please tell The heat source or which surface has which temperature...

Why you are using the solver potetialFoam in your example? As I understood you want to use chtMultiRegionFoam! You need to change the solver (application) in controlDict to chtMultiRegionFoam!

There is no T dictionary in your 0.orig folder... This needs to be defined if you are running a conjugate heat transfer!

Here is an example for chtMultiRegionSimpleFoam with heat source in region LED!

Regards

Peter

Here is an example like your case!

Heat source is in SOURCE! fvOptions!

Regards

Peter

Sorry Peter, I just uploaded a wrong file. This is the correct file. And thank you for the test file. I am new to OpenFOAM. By the way, I did look at these topics and I just changed the diameter of the cylinder and add PIMPLE solver. I used Sandymech's dictionaries. I didnt use yours because you did with snappyHexMesh. but I need to change the diameter and length of the cylinder.


p { margin-bottom: 0.1in; line-height: 120%; }a:link { } https://www.cfd-online.com/Forums/op...t-transfe.html

https://www.cfd-online.com/Forums/openfoam-solving/191550-conjugate-heat-transfer-heat-generating-cylinder-natural-convection.html

Hello Cagatayemre!

In your case (the last one you sent!) you need to change the file fvSolutions in your system/fluid with one that solve a simple instead of pimple!

Take simply the file in my example LUM from system/INNEN/fvSolutions and replace it in your case.

Then the simulation (yours) will run!

Anyway, you will have a divergence after 3 itterations :)

Somthing still going wrong here... The temperature field is exploding...

Decreasing the relaxationFaktors of h to 0.7 (for fluid) will let your case works fine!

Regards

Peter

Thank you very much Peter but now I have following error. I am trying to solve it with chtMultiRegionFoam. Is it correct ? Which command should I type? I use
blockMesh -region solid
blockMesh -region fluid
changeDictionary - region solid
changeDictionary - region fluid
chtMultiRegionFoam
keyword PIMPLE is undefined in dictionary "/home/user2/OpenFOAM-v1706/tutorials/OpenFoamRunFiles/flowAroundHeatedCylinder/cylnatural_send4/system/fluid/fvSolution"

file: /home/user2/OpenFOAM-v1706/tutorials/OpenFoamRunFiles/flowAroundHeatedCylinder/cylnatural_send4/system/fluid/fvSolution from line 21 to line 60.

Sure it will not work! All setups are made for chtMultiRegionSimpleFoam!

To do that you need to do the following:

- Replace your

controlDict

with the one in:

../tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/


- Replace:

fvSolution

and

fvSchemes

for both solid and fluid with those in:

../tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/heater/

and

,,/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/bottomWater/


Then it will work :)


Regards

Peter

Thank you very much Peter. Now it iterates without any error. I dont know how can I give a reputation point for you. Where is the button do you know?