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

Foam::divide, transient conjugate heat transfer

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By peterhess

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 9, 2018, 08:16
Default Foam::divide, transient conjugate heat transfer
  #1
Member
 
Çağatay Emre Ayhan
Join Date: Sep 2017
Location: Istanbul, Turkey
Posts: 31
Rep Power: 8
Cagatayemre is on a distinguished road
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:robes::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:perator/(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 ()
Attached Files
File Type: gz cylinder.tar.gz (9.7 KB, 3 views)

Last edited by Cagatayemre; January 9, 2018 at 09:56.
Cagatayemre is offline   Reply With Quote

Old   January 9, 2018, 09:04
Default
  #2
Senior Member
 
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17
peterhess is on a distinguished road
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
Attached Files
File Type: gz LUM.tar.gz (10.4 KB, 5 views)
peterhess is offline   Reply With Quote

Old   January 9, 2018, 09:16
Default
  #3
Senior Member
 
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17
peterhess is on a distinguished road
Here is an example like your case!

Heat source is in SOURCE! fvOptions!

Regards

Peter
Attached Files
File Type: gz Test.tar.gz (8.1 KB, 5 views)
peterhess is offline   Reply With Quote

Old   January 9, 2018, 09:49
Default
  #4
Member
 
Çağatay Emre Ayhan
Join Date: Sep 2017
Location: Istanbul, Turkey
Posts: 31
Rep Power: 8
Cagatayemre is on a distinguished road
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 { } Volumetric Heat generation in cylinder with natural convection conjugate heat transfe

https://www.cfd-online.com/Forums/openfoam-solving/191550-conjugate-heat-transfer-heat-generating-cylinder-natural-convection.html
Attached Files
File Type: gz cylinder.tar.gz (9.7 KB, 6 views)
Cagatayemre is offline   Reply With Quote

Old   January 9, 2018, 10:13
Default
  #5
Senior Member
 
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17
peterhess is on a distinguished road
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
Attached Files
File Type: gz cylinder_Works.tar.gz (8.7 KB, 7 views)
peterhess is offline   Reply With Quote

Old   January 9, 2018, 10:50
Default
  #6
Member
 
Çağatay Emre Ayhan
Join Date: Sep 2017
Location: Istanbul, Turkey
Posts: 31
Rep Power: 8
Cagatayemre is on a distinguished road
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.

Quote:
Originally Posted by peterhess View Post
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
Cagatayemre is offline   Reply With Quote

Old   January 9, 2018, 11:14
Default
  #7
Senior Member
 
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17
peterhess is on a distinguished road
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
Attached Files
File Type: gz cylinder_Pimple_Works.tar.gz (8.6 KB, 8 views)
Cagatayemre likes this.
peterhess is offline   Reply With Quote

Old   January 9, 2018, 12:13
Default
  #8
Member
 
Çağatay Emre Ayhan
Join Date: Sep 2017
Location: Istanbul, Turkey
Posts: 31
Rep Power: 8
Cagatayemre is on a distinguished road
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?
Quote:
Originally Posted by peterhess View Post
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
Cagatayemre is offline   Reply With Quote

Reply


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
Thermal non-equilibrium porous media model with conjugate heat transfer Hexahedron FLUENT 9 February 22, 2023 02:55
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 07:30
Inverse and Transient Heat Transfer Problem on commercial software: is it possible? rogbrito Main CFD Forum 1 February 19, 2019 02:11
Convective / Conductive Heat Transfer in Hypersonic flows enigma Main CFD Forum 2 November 1, 2009 22:53
Conjugate heat transfer problem hvem10 FLUENT 2 October 29, 2009 17:31


All times are GMT -4. The time now is 00:09.