CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Alter turbulence parameters

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 28, 2006, 07:20
Default Hi all, Is it theoretically
  #1
newbee
Guest
 
Posts: n/a
Hi all,

Is it theoretically correct to alter the turbulece parameters in order to get better resemlence to experemental data and if so what guidelines are there for doing this.

I doent seem to get results close enough to real values.

/Erik
  Reply With Quote

Old   August 28, 2006, 08:12
Default No, it's not correct! Of cours
  #2
Member
 
Ola Widlund
Join Date: Mar 2009
Location: Sweden
Posts: 87
Rep Power: 7
olwi is on a distinguished road
No, it's not correct! Of course, it depends on your application, you may of course stumble on a very special case for which your model is not appropriate. If you are still talking about a simple channel flow, however, then it's a very bad idea to change constants. Standard turbulence ARE tuned to do that very well. It may give you a "better" result, but for the wrong reason. Two wrongs do not make one right. You should try to focus on what is so special with your case, if there are things you didn't think of:
- Compressible?
- Do you have buoyancy effects?
- Is it single phase?
- Chemical reactions/combustion?
- Swirling flows?
- Permeable walls?
- Wall jets?
- Is it really steady-state? If you think it is, do you get the same solution at the end of a transient run?
- Are you using symmetries where you shouldn't?
- Are your measurements really that good? Do they agree with other results in the literature?
- Are your inlet and outlet conditions placed in a good position upstream/downstream? Are they placed well in the experiment?!

If you want to read more about how the model constants are chosen you could check "Turbulence modeling for CFD" by David C Wilcox, or "Turbulent Flows" by Stephen Pope. Wilcox is good for this question, I think.

/Ola


/Ola
olwi is offline   Reply With Quote

Old   August 28, 2006, 08:35
Default Thank you Ola for the thorough
  #3
newbee
Guest
 
Posts: n/a
Thank you Ola for the thorough answer. It got me thinking of some new things i shall try. I am modelling a air flow with a incompressible because im hopeful to use the model for water flow aswell. will also try a compressible solver

Thanks again
/Erik
  Reply With Quote

Old   August 28, 2006, 08:51
Default Hi again, Air will indeed b
  #4
Member
 
Ola Widlund
Join Date: Mar 2009
Location: Sweden
Posts: 87
Rep Power: 7
olwi is on a distinguished road
Hi again,

Air will indeed be incompressible, as long as you have a Mach number of less than 0.2-0.3. That means it's ok for speeds below 100 m/s, if you're at ambient pressure. Feel free to tell us a bit more about your case, and about the expriment you compare with. Maybe that will give us some ideas...

/Ola
olwi is offline   Reply With Quote

Old   August 28, 2006, 11:24
Default Thanks a lot! This is the g
  #5
newbee
Guest
 
Posts: n/a
Thanks a lot!

This is the geometry of the case that I run my simulations on:



It is 1/12 of the circumference of a rod. the three other sides are symmetry sides. I have used wallFunctions for the rod boundary, inlet for the inlet and inletOutlet for the outlet. The channel is 11.5 meter long and has a inside air flow velocity of 20.57 m/s. The Re number is 64590.
I count the air flow as incompressible and steady state.


this attached papper explains the experiment that I attempt to reconstruct:



As the papper mentions the next step is to heat the flow through the wall with the following implemented code:
solve
(
fvm::div(phi, T)
- fvm::laplacian(alphaEff, T)
- S*q/(rho*Cp)
- (mu*(gradU + gradU.T()) && gradU)/(rho*Cp)
- turbulence().epsilon()()/Cp
);

where S is the cell wallarea divided by the cell volume. q is the applied heat/area. The last two term hopefully models vilcous dissipation heating.

I hope I havent shoveled too much uninteresting information on you.

Thankful
/Erik
  Reply With Quote

Old   August 28, 2006, 11:47
Default the papper was to big to send.
  #6
newbee
Guest
 
Posts: n/a
the papper was to big to send. I will also lay up my graph results once iv'e finnished my last turbmodel.
  Reply With Quote

Old   August 28, 2006, 12:38
Default The flow is oriented thou ward
  #7
newbee
Guest
 
Posts: n/a
The flow is oriented thou wards the screen in the colored picture above.

here are my graph results:

sampling on the narrow symmetry side to the left

sampling on the wide symmetry side to the right


where the red crosses are from the experement and the other signs are from the models

Should I expect better coherence?
  Reply With Quote

Old   August 29, 2006, 03:31
Default I would expect velocity curves
  #8
Member
 
Ola Widlund
Join Date: Mar 2009
Location: Sweden
Posts: 87
Rep Power: 7
olwi is on a distinguished road
I would expect velocity curves to be a lot more smooth... You sure it's properly converged? I wouldn't be surprised to see a difference between the lines, but their shape should be even and smooth, as in the experiment.

Do you model the whole 11.5 meters of the rod bundle? Do you make the grid using blockMesh? You can make cells very long in the streamwise direction (200-300 mm, maybe even longer), that would save speed and improve convergence. Feel free to send the blockMeshDict to my private mail address, maybe I can have a look: olwix@yahoo.com

/Ola
olwi 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
parameters of turbulence models george CD-adapco 2 January 16, 2007 10:00
Turbulence model parameters Vidya Raja FLUENT 2 October 6, 2006 14:38
Turbulence Parameters Bob FLUENT 0 July 28, 2004 16:32
turbulence parameters Lio Main CFD Forum 0 July 12, 2004 05:38
Turbulence parameters Mikael Ersson Main CFD Forum 0 May 13, 2004 22:28


All times are GMT -4. The time now is 19:46.