CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Pre-Processing

Channel Flow with Heat Flux

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 18, 2009, 15:00
Default Channel Flow with Heat Flux
  #1
Member
 
Sven Schweikert
Join Date: Jun 2009
Posts: 38
Rep Power: 8
svens is on a distinguished road
Hello

I am an absolutely OF beginner and to learn the ropes I'm trying to create and solve a simple 2D laminar channel flow with heated walls.

-What solver fits best for my problem? I successfully calculated laminar flow in same channel geometry without heat flux with the icoFoam solver. Is there a way to modify icoFoam to use it also for my new case with heat transfer?

-I tried to solve my problem with the buoyantFoam solver. At some point the Courant-Number exploded which led to an abort of calculation. Is there a problem using buoyantFoam for a laminar case?

-Is there a case comparable to mine which I can use as assistance to get a little bit easier in the problem?

Thanks so much - svens
svens is offline   Reply With Quote

Old   June 19, 2009, 08:54
Default
  #2
New Member
 
Quentin
Join Date: May 2009
Posts: 22
Rep Power: 8
Bedotto is on a distinguished road
Hi Svens

For your first point, you should go on the OpenFoam Wiki. There is an example about how to implement the heat equation in a solver with icoFoam:
http://openfoamwiki.net/index.php/Ho...ure_to_icoFoam

For the problem concerning your courant number, sorry I've no idea (I'm a beginner like you )

However if you want to compare your results with a DNS , you should have a look on this site:
http://murasun.me.noda.tus.ac.jp/turbulence/

Regards

Quentin
Bedotto is offline   Reply With Quote

Old   June 19, 2009, 12:19
Default
  #3
Member
 
Sven Winkler
Join Date: May 2009
Posts: 70
Rep Power: 8
sven is on a distinguished road
I also faced the same problem as svens. I used the manual from the link posted. It worked, but I still have some questions on this manual:

1. What is the equation that is implemented in step 4? This is not the normal energy equation, is it?

2. What is the DT specified in the transport properties (step 5)?

Thanks for your help.
sven is offline   Reply With Quote

Old   June 19, 2009, 13:23
Default
  #4
Member
 
Sven Schweikert
Join Date: Jun 2009
Posts: 38
Rep Power: 8
svens is on a distinguished road
Thanks Quentin for your reply - it seems to work.

Unfortunately I've got the same questions like sven to get a fully understanding what I did. I would be really grateful about some further explanations.

Another question which occurred is the lack of an implementation for the Thermophysical Model which is performed in Chapter7 of the OF Users Guide.

I guess it's necessary to create such a model for the complete definition of this case, isn't it?

Thanks a lot for your help & greetings.
svens

couple of hours later:

Now I am pretty sure that DT labels the thermal diffusion also know as alpha [=k/(roh*c_p)].

With DT my fluid properties are completely defined and I guess the above mentioned thermophysical model isn't necessary anymore.

But I still have some huge problems to understand how the energy equation is implemented in my new icoFoam solver.

fvScalarMatrix TEqn
(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(DT, T)
);
TEqn.solve()
;

About some further explanatinons of the implemtation I would be really grateful.

Thanks a lot




Last edited by svens; June 19, 2009 at 20:48. Reason: new results
svens is offline   Reply With Quote

Old   June 22, 2009, 06:04
Default
  #5
New Member
 
Quentin
Join Date: May 2009
Posts: 22
Rep Power: 8
Bedotto is on a distinguished road
fvScalarMatrix TEqn
(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(DT, T)
);
TEqn.solve()
;

In the tutorial these lines are written in the file .C, more precisely in the time loop. The term fvm means a numerical scheme used to discretize the derivative operation. There's a good example on the wiki which shows the 'traduction' of a equation written in OpenFoam. Here's the link:
http://en.wikipedia.org/wiki/OpenFOAM

Tell if I'm wrong.

Regards
Quentin
Bedotto is offline   Reply With Quote

Old   June 22, 2009, 10:40
Default
  #6
New Member
 
Quentin
Join Date: May 2009
Posts: 22
Rep Power: 8
Bedotto is on a distinguished road
By the way, I would like to add to this equation the term nu_t/Pr_t where nu_t is the SGS model term and Pr_t is the turbulent Prandt. But I don't know how to "call" the nu_t variable. Does someone know how to implement this term?

Regards

Quentin
Bedotto is offline   Reply With Quote

Old   June 22, 2009, 11:01
Default
  #7
Senior Member
 
santos's Avatar
 
Jose Luis Santos
Join Date: Mar 2009
Location: Portugal
Posts: 213
Rep Power: 9
santos is on a distinguished road
Send a message via Skype™ to santos
Hi,

I think you may call it with:

Code:
sgsModel->nuSgs()
Regards,
Jose Santos
santos is offline   Reply With Quote

Old   June 23, 2009, 10:13
Default
  #8
New Member
 
Quentin
Join Date: May 2009
Posts: 22
Rep Power: 8
Bedotto is on a distinguished road
I've checked it, it works.
Thanks Santos!
Bedotto is offline   Reply With Quote

Reply

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
Yacht in Open Channel Flow andreimour FLUENT 1 October 14, 2010 23:54
Heat Flux Wall Boundary Confusion. Joee FLUENT 1 August 21, 2010 12:20
Wall Heat Flux not a listed variable mike CFX 8 February 5, 2009 14:32
Concentric tube heat exchanger (Air-Water) Young CFX 5 October 6, 2008 23:17
compressible channel flow.. R.D.Prabhu Main CFD Forum 0 July 17, 1998 17:23


All times are GMT -4. The time now is 11:15.