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

2D counter flow heat exchanger analysis using chtMultiRegionSimpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 31, 2013, 05:31
Default 2D counter flow heat exchanger analysis using chtMultiRegionSimpleFoam
  #1
New Member
 
Join Date: Jun 2013
Posts: 5
Rep Power: 4
dprasad_p is on a distinguished road
Dear Foamers,

I am doing Maters project in CFD using openFOAM. As a part of a project, I am working on chtMultiRegionSimpleFoam solver using the test case "2D counter flow heat exchanger" from the literature "A note on the solution of conjugate heat transfer problems using SIMPLE-like algorithms. Xi Chen *, Peng Han"

I am getting the fluid part temperatures almost correct but a difference of 50K is observed in solid part. The difference is constant thought the solid part. I have checked BCs and input values and they are ok. I am clueless where the problem is.

I have attached the case and the comparison graph. can anyone please go through the case and suggest?

Thanks in advance,
Prasad
Attached Images
File Type: jpg openFOAM-Results.jpg (43.7 KB, 22 views)
Attached Files
File Type: zip heatExCounterWater.zip (67.6 KB, 19 views)
dprasad_p is offline   Reply With Quote

Old   September 7, 2013, 08:46
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings Prasad and welcome to the forum!

It looks like:
  1. You're using OpenFOAM 2.1.x.
  2. The Allrun script seems to have been edited on Windows.
  3. You are using blockMeshDG for generating the mesh: http://openfoamwiki.net/index.php/Contrib_blockMeshDG
This kind of information is crucial for anyone to help you!

About your case:
  • It looks similar to this one: http://openfoamwiki.net/index.php/Ge..._-_planeWall2D - this is an incomplete example. But you might want to have a look into the chapter "4 References" on that page.
  • I have ran your case with openFOAM 2.1.s, but I don't know where you are measuring the Temperature.
  • Therefore, all I can do is guess:
    1. The settings in "system/*/fvSchemes" for the air regions might not be well done, because most of the divergence schemes are set to 1st order, namely "upwind". Switching to second order might help.
    2. You might not have enough mesh resolution near the walls of the solid.
Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   September 8, 2013, 06:25
Default
  #3
New Member
 
Join Date: Jun 2013
Posts: 5
Rep Power: 4
dprasad_p is on a distinguished road
Hi Bruno,

Thanks for your time. I am using blueCFD2.1.1. I have taken the values across the HE at L = 0.5.

I have tried couple of variations in numerical schemes and got mixed results. First I tried with SFCD, linearUpwind and both had no improvement in the results. With QUICK scheme, the results were improved and the deviation this time was around 10K. :-) The mesh refinement in the fluid and solid regions did not help much.

Revised comparison graph is attached.

I have a quick question. Interpolation scheme 'linearUpwind' is described as 'first/second order, bounded' and how do we ensure to enforce first or second order specifically? Or on what basis program will consider the specific order?

Thanks,
Prasad.
Attached Images
File Type: jpg Revised_Results.jpg (20.3 KB, 8 views)
dprasad_p is offline   Reply With Quote

Old   September 8, 2013, 09:48
Default
  #4
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Dear Prasad,

I don't have access to the paper, so I don't know the specifics of how the paper generated the data, in order to estimate what might be the best approach. It could also be a problem with the turbulence model or fluid properties.

As for "linearUpwind", a quick search lead me to this thread, although I did not fully read it: A few simple questions about linearUpwind and limitedLinear

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   September 14, 2013, 04:16
Default
  #5
New Member
 
Join Date: Jun 2013
Posts: 5
Rep Power: 4
dprasad_p is on a distinguished road
Dear Bruno,

Sorry for the late reply. I got tied up with some other work.

The case is based on laminar flow and fluid properties are given as stated in the literature. I could not attach the literature document because of large file size. Instead please download it from the below dropbox link.

http://db.tt/6qrj6vRH

I presume chtMultiRegionSimpleFoam can be used for incompressible fluids such as water. I searched in forum and found no concluded answer or did miss the right post?

I am also trying to work out on incompressible version of cht...foam chtIcoMultiRegionSimpleFoam which is available in forum. I am finding it difficult to set up a test case. Is there any test case already run for this solver? I could not locate one on forum.

Thanks,
Prasad
dprasad_p is offline   Reply With Quote

Old   September 14, 2013, 14:00
Default
  #6
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Prasad,

Sorry, I don't have much patience to look at this with much more attention (it's been a very long week), but the following properties look suspicious to me, at least on the solid side:
Code:
constSolidThermoCoeffs
{
    //- thermo properties
    rho rho [1 -3  0  0 0 0 0] 8000;// kg / m^3
    Cp  Cp  [0  2 -2 -1 0 0 0] 500; // m^2 / ( sē . K )
    K   K   [1  1 -3 -1 0 0 0] 50;  // kg . m / ( s^3 . K )

   //[omitted the other entries]
}
The comments about the units are based on this: http://www.openfoam.org/docs/user/ba...18-1000004.2.6

The paper stipulates that the solid part has the following properties:
Quote:
metallic plate properties
qs = 8000 kg/m3
ks = 50 W / m.K
Cps = 500 J / (kg . K)
I'm too tired to double-check the units, but my suspicion is that the unit conversion is not direct for all 3 properties.


As for the other solver you mentioned: I'm not familiar with it.

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Reply

Tags
chtmultiregionsimplefoam, conjugate heat transfer

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
two-phase flow in compact heat exchanger kakelliravi FLUENT 3 April 3, 2010 12:46
Heat Exchanger Simulation tightaznbreaka FLUENT 0 March 18, 2010 16:45
Concentric tube heat exchanger (Air-Water) Young CFX 5 October 6, 2008 23:17
Counter flow heat exchanger with two different liq Ram Dayal CFX 3 August 6, 2006 20:42
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 16:55


All times are GMT -4. The time now is 23:37.