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

conjugate heat transfer in OpenFOAM

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

Like Tree28Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 25, 2016, 03:38
Default
  #81
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Thank you very much Paritosh and Bruno! I understand now the "*.OpenFOAM" file is an empty file for the paraView to know where to read the data. The "paraFoam -builtin" command combines all the "*.OpenFOAM" files into one file, and it works fine for the tutorial case "multiRegionHeater".

As Bruno as said, it is better to start with a simple case for which analytical solution is possible. I still can't get the files "planeWall2D{bottomAir}.OpenFOAM", "planeWall2D{topAir}.OpenFOAM" and "planeWall2D{wall}.OpenFOAM". I continue to compare the "planeWall2D" case with the "multiRegionHeater" to find out the reason.
sitajeje is offline   Reply With Quote

Old   September 26, 2016, 04:47
Default Try this setup
  #82
Senior Member
 
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 732
Rep Power: 20
vasava will become famous soon enough
Try this setup. I have made the following changes.
  1. The first thing I notice is that the blockMeshDict file is in the system folder. I added following line to Allrun file to copy the file to constant/polymesh folder.
    Code:
    cp system/blockMeshDict constant/polyMesh/blockMeshDict
  2. I combined Allrun and Allrun.pre files. This is personal preference, I like all the commands in one single file.
  3. Changed the region name from 'wall' to 'solidWall'. This is do avoid confusion with the boundary type and any other boundary named wall or consisting word 'wall'. I also made corresponding name changes in all the files.

I am using openFoam 2.3 and the case ran without any errors.

If there are errors, don't forget to attach the log files.

Edit:
If you want to learn how things work with Allrun-type files, I advise you to run the case in steps. You can start with Allrun file with one command. Run the case then read and understand the log file. If everything is alright, add the next command and so on.
Attached Files
File Type: gz planeWall2D_corr.tar.gz (6.6 KB, 13 views)

Last edited by vasava; September 26, 2016 at 04:55. Reason: edit
vasava is offline   Reply With Quote

Old   September 26, 2016, 11:45
Default
  #83
Senior Member
 
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 19
zfaraday will become famous soon enough
Quick answer: Check out the "log.<solverName>" file, there you will find the error that occurred!
__________________
Web site where I present my Master's Thesis: foamingtime.wordpress.com

The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in!
zfaraday is offline   Reply With Quote

Old   October 3, 2016, 10:46
Default
  #84
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Dear Paritosh,

Thank you very much for your help! It is very nice of you! Sorry that I couldn't reply earlier, it was too busy these days. I have several questions about fvSchemes in system/topAir.

What is div(Ji,Ii_h) please? I didn't figure it out with google. I understand div(((rho*nuEff)*dev2(T(grad(U))))) the same term as div((muEff*dev2(T(grad(U))))). Am I right? If so, why the both should be written? These two terms do not exist in Bruno's planeWall2D case.

I don't understand the divSchemes completely.
div(phi, U) is the advection of velocity,
what is div(phi, K) please?
div(phi, h) is the advection of enthalpy,
what is div(phi, k) please?
div(phi, K) is advection of the turbulent kinetic energy,
div(phi, epsilon) is advection of the dissipation rate
div(phi, R) is advection of the Reynolds stress,
div(R) is the divergence of the Reynolds stress.

Thank you very much in advance!
sitajeje is offline   Reply With Quote

Old   October 5, 2016, 10:51
Default
  #85
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Dear Paritosh,

Thank you very much for your help! I have run your "planeWall2D_corr", but it didn't work out completely. "log.blockMesh", "log.topoSet", "log.splitMeshRegions" reported no error. The error began in "log.changeDictionary.solidWall" as

--> FOAM FATAL IO ERROR
cannot find file
file: .../planeWall2D_corr/0/solidWall/dictionaryReplacement at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 72

FOAM exiting

I can run through the tutorial case "multiRegionHeater". I compared the two cases and found that in "mutiRegionHeater" the "changeDictionaryDict" in e.g. "system/bottomWater" begins without "dictionaryReplacement". Why must this function be in planeWall2D please?

I didn't figure out how to attach files and pictures in my Reply in the forum. Could you give me a hint please?

Thank you very much in advance! I really appreciate your help!
sitajeje is offline   Reply With Quote

Old   October 6, 2016, 10:43
Default
  #86
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Dear Paritosh,

It is quite strange for me that, I got the following error when I ran blockMesh in the original case files of Bruno,

fileName:: stripInvalid() called for invalid fileName Terminaloutput
For debug level (=2) > 1 this is considered fatal
Aborted (core dumped)

But the same blockMeshDict ran successfully in your planeWall_corr case files. I checked the directory, all the folder names are with letters and numbers without spaces. Did you have a similar experience?

I am new in OpenFOAM, I believe I will have quite a few strange mistakes. I would appreciate if you could help me on this. Thank you very much in advance!
sitajeje is offline   Reply With Quote

Old   October 6, 2016, 16:11
Default
  #87
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Dear Paritosh,

No matter how I rename the folder of Bruno's original case, blockMesh never worked. However, I created a new folder with a valid name as "planeWallTest", and copied all the files of Bruno to this new folder, the blockMesh functioned! Do you have an idea why it is so folder dependent? I suspect whether I have dong something wrong with the installation.
sitajeje is offline   Reply With Quote

Old   December 31, 2016, 11:35
Default
  #88
New Member
 
Sandeep Pandey
Join Date: Jul 2016
Location: Germany
Posts: 10
Rep Power: 7
sandeep.pandey is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Greetings skuznet and welcome to the forum!

Have a look into this (unofficial) tutorial: http://openfoamwiki.net/index.php/Ge..._-_planeWall2D - and feel free to add information to it, when you've understood better how things work!

As for Abaqus, I have no idea if there is a direct way. You can try looking here in the forum. You can even use Google:
Code:
site:www.cfd-online.com/Forums/ Abaqus OpenFOAM
If there is no such direct conversion, have a look into this wiki page and the links therein: http://openfoamwiki.net/index.php/SnappyHexMesh

Best regards,
Bruno
Hey Bruno,
Do you have this CHT tutorial for OpenFOAM v4.1? Or can you suggest what changes is to be made?

Thanks and HNY,
Sandeep
sandeep.pandey is offline   Reply With Quote

Old   January 3, 2017, 13:55
Default
  #89
Member
 
Jaydeep
Join Date: Jun 2015
Posts: 34
Rep Power: 8
jaydeep is on a distinguished road
@sandeep.pandey

I am making a short tutorial video for the Planewall2D case with OF 4.1. I will post a link in this thread in a couple of days.
LeeRuns likes this.
jaydeep is offline   Reply With Quote

Old   January 3, 2017, 15:38
Default
  #90
New Member
 
Sandeep Pandey
Join Date: Jul 2016
Location: Germany
Posts: 10
Rep Power: 7
sandeep.pandey is on a distinguished road
Quote:
Originally Posted by jaydeep View Post
@sandeep.pandey

I am making a short tutorial video for the Planewall2D case with OF 4.1. I will post a link in this thread in a couple of days.
Hey Jaydeep,
Thanks for taking this initiative, eagerly waiting!
sandeep.pandey is offline   Reply With Quote

Old   January 5, 2017, 06:59
Default planeWall2D case for OpenFOAM4.0
  #91
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Happy New Year everybody,

I have made a more simplified version of the planeWall2D. I would appreciate if anyone would like to comment on the resulting temperature and velocity fields. Are they correct?

I wonder how to display the velocity in the solid domain as 0 instead of red in my image?

When I plot the velocity vector field with Glyph in paraView, there are always two arrows in one cell. I attached my state file, any suggestions would be very helpful for me.

I wonder how can I monitor the residuals during computing. I tried
Code:
functions
{
#includeFunc residuals
}
and add
Code:
foamMonitor -l postProcessing/residuals/0/residuals.dat
below the "runApplication 'getApplication' " in "Allrun", but it was unsuccessful with the following text in the terminal:
Code:
File postProcessing/residuals/0/residuals.dat does not exit
Does anyone know how to manage it?

Thank you very much in advance,

All the best,

Sitajeje
Attached Images
File Type: jpg boundaryCondition.jpg (51.2 KB, 57 views)
File Type: jpg velocity.jpg (53.2 KB, 48 views)
File Type: jpg temperature.jpg (43.8 KB, 42 views)
Attached Files
File Type: gz planeWall2DSitajeje.tar.gz (5.8 KB, 16 views)
File Type: zip state file paraView.zip (16.2 KB, 4 views)
sitajeje is offline   Reply With Quote

Old   January 11, 2017, 16:24
Default
  #92
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Happy New Year Alex,

your master thesis is really impressive! I learned a lot from your descriptions and case files, especially about the fvOption! Thank you very much for sharing your master thesis with the public!

Best,

Sitajeje
sitajeje is offline   Reply With Quote

Old   January 12, 2017, 04:33
Default
  #93
Senior Member
 
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 19
zfaraday will become famous soon enough
Happy New Year,

I'm glad you find my site useful, Sitajeje, and I really appreciate your thoughts about it. I hope more people can find a little help in my web site.

Best regards,

Alex
__________________
Web site where I present my Master's Thesis: foamingtime.wordpress.com

The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in!
zfaraday is offline   Reply With Quote

Old   January 14, 2017, 03:20
Default
  #94
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Hallo Alex,

I am struggling with adding a heat source to chtMultiRegionSimpleFoam with fvOption. I divided my computation domain into 2 blocks and gave name to them. I found that topoSet is then not necessary. Am I right?

I copy your fvOption file into constant/heater and set cellZone to heater. blockMesh, splitMeshRegions, changeDictionary ran through till chtMultiRegionSimpleFoam. The error appears as

keyword selectionMode ist undefined in ..../fvOption.volumetricHeatSource.scalarSemiImplicitSo urceCoeffs

The same error occurred when I added selectionMode into the dictionary scalarSemiImplicitSourceCoeffs.

I found quite a few your comments and questions on chutMultiRegionSimpleFoam. You mentioned in a thread in 2014 that the source term is of no effect in the new chtMultiRegionSimpleFoam solver. Is it true? I would appreciate greatly if you would share your experience with me.

Thank you very much in advance!

Best regards,

sitajeje
sitajeje is offline   Reply With Quote

Old   January 14, 2017, 10:17
Default adding heat source to chtMultiRegionSimpleFoam with fvOptions
  #95
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Dear Foamers,

I want to add heat source to the solid domain while calculating with the chtMultiRegionSimpleFoam solver. I followed the following link which is a discussion of Bruno with another problem reporter on adding heat source with fvOptions.

https://bugs.openfoam.org/view.php?id=1266

The test case is the same as my post #91, except that fvOptions is added in system/heater and the bottom face of the heater is changed to zeroGradient instead of fixedValue in temperature.

I met a strange error as

Code:
keyword selectionMode ist undefined in ..../fvOption.heatSource.scalarSemiImplicitSourceCoeffs
I would appreciate greatly if someone can help me!

Best regards,

sitajeje
Attached Files
File Type: gz planeWall2Dsitajeje.tar.gz (50.3 KB, 5 views)
sitajeje is offline   Reply With Quote

Old   January 14, 2017, 11:10
Default
  #96
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
got it done. It was a version problem. The case from Bruno is excellent!
sitajeje is offline   Reply With Quote

Old   January 19, 2017, 18:02
Default
  #97
mom
New Member
 
Join Date: Jan 2017
Posts: 5
Rep Power: 6
mom is on a distinguished road
Quote:
Originally Posted by sitajeje View Post
got it done. It was a version problem. The case from Bruno is excellent!
Hi sitajeje,
I want to use chtMultiRegionFoam for a problem. But I thought it is better to start with chtMultiRegionSimpleFoam. Since the old files are not compatible, I don't know how to modify them. Did you apply it on OF v 4.0? Do you know whether it is compatible with v 4.1 or not?
Is it possible for you to share your file?

mom
mom is offline   Reply With Quote

Old   January 23, 2017, 06:17
Default
  #98
mom
New Member
 
Join Date: Jan 2017
Posts: 5
Rep Power: 6
mom is on a distinguished road
Quote:
Originally Posted by jaydeep View Post
@sandeep.pandey

I am making a short tutorial video for the Planewall2D case with OF 4.1. I will post a link in this thread in a couple of days.
Dear Jaydeep,
I am eagerly waiting for your tutorial video
Thank you all for allocating time and sharing your knowledge

mom
mom is offline   Reply With Quote

Old   March 14, 2017, 16:08
Default reversal flow at the outlet in a free convection
  #99
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Dear OpenFOAMers,

I met a problem with chtMultiRegionFoam while simulation natural convection. I believe I am not the only person who met this difficulty. I would appreciate if anyone can help me. It is going to be helpful for our community.

Please find my computation domain in the attached image. It is a 2D simulation. The domain is 25mm wide and 135mm long. The Grashof number is roughly 1e6 < 1e9, so the flow should stay laminar. The left border is an insulated wall. A rectangular solid with a heat source of 1W locates next to the wall. The other three sides are atmosphere with totalPressure, p0=1e5 for p_rgh and pressureInletOutletVelocity, value uniform (0 0 0) for u.

I would expect a thermal boundary layer next to the wall, which extends out of the domain with a weak scattering after passing the heating solid. However the boundary layer decreases its speed pretty far away from the outlet and becomes reversal flow next to the right border. I doubled the width of the computation domain, and the same result occurs: the flow reverses and turn down next to the right border. This indicates that the simulation is wrong because of the boundary condition.

I wonder whether anyone has experience with setting the boundary conditions of atmosphere for free convection? Any suggestions, comments would be very helpful for me. Thank you very much in advance!

Best regards,

Sitajeje
Attached Images
File Type: jpg Screenshot from 2017-03-14 21_38_09.jpg (154.3 KB, 44 views)
File Type: jpg Screenshot from 2017-03-14 21_35_31.jpg (103.5 KB, 55 views)
Attached Files
File Type: gz HSPlaminar2.tar.gz (6.9 KB, 14 views)
sitajeje is offline   Reply With Quote

Old   March 16, 2017, 05:07
Default
  #100
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 7
sitajeje is on a distinguished road
Dear Foamers,

Problem is solved. The temperature should be zeroGradient at the upper border instead of inletOutlet.

Thank you all the same!

Best regards,

Sitajeje
LeeRuns likes this.
sitajeje is offline   Reply With Quote

Reply

Tags
cht, solid-fluid interface

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
Conjugate Heat Transfer: Wall Heat Flux at Coupled Walls? MaxHeat FLUENT 4 September 14, 2017 10:44
openfoam for heat transfer kirankarki OpenFOAM Running, Solving & CFD 29 February 12, 2015 18:46
Conjugate Heat transfer in CFX ksp1717 CFX 11 December 10, 2010 22:07
Conjugate Heat Transfer of Motorized EGR enr_venkat CFX 1 October 12, 2010 18:17
best mesh generator for conjugate heat transfer? phsieh2005 Main CFD Forum 1 June 1, 2007 17:35


All times are GMT -4. The time now is 05:03.