CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > Siemens

Modelling a turbulent jet and k-epsilon constants

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 24, 2005, 06:22
Default Modelling a turbulent jet and k-epsilon constants
  #1
Ant
Guest
 
Posts: n/a
Hi,

I am trying to model a turbulent jet using STAR CD. Initially I have used the high Re number turbulence model and this has given reasonably good results, but not good enough. As one might expect, the spreading rate is too high, and I want to implement the Launder, Morse, Rodi and Spalding model which varies the constants C_{mu} and C_{epsilon 3} as:

C_{mu} = 0.09 - 0.04f

C_{epsilon 3} = 1.92 - 0.0667f

where f = | (R_{mu}/1.6u_{0} |du_{cl}/dx - |du_{cl}/dx|||

R_{mu} is the radius of the mixing region u_{cl} is the centreline velocity u_{0} is the jet velocity.

More details are in

Launder, B.E. and Spalding, D.B., " The numerical computation of turbulent flows", Computer Methods in Applied Mechanics and Engineering 3 (1974) pp 269-289

available online.

I need to use this model rather than a more up to date one to compare my results with those of my supervisor. How can I do this using STAR CD? I have seen that you can make the Prandtl number variable, and also you can write a user routine to supply turbulent viscosity, but is there any way so that one can write and use one's own turbulence model, or simply let the constants vary?

Thanks in advance for your responses,

Ant
  Reply With Quote

Old   January 24, 2005, 07:25
Default Re: Modelling a turbulent jet and k-epsilon consta
  #2
Richard
Guest
 
Posts: n/a
It's possible ...

For C_mu, you could use the vistur user subroutine. The value of vist which is passed into this subroutine will be the one using the standard C_mu, so all you would need to do would be to multiply this by C_mu_new/C_mu.

For C_eps3 you would need to re-implement the source term through sorkep. The sample coding in that subroutine should give you an idea how to do it.

If you haven't thought of this already: make sure you test your subroutines by implementing the standard k-epsilon model through them first.
  Reply With Quote

Old   January 24, 2005, 12:21
Default Re: Modelling a turbulent jet and k-epsilon consta
  #3
Ant
Guest
 
Posts: n/a
Thanks for your help. I think I'll be able to do it now. I'll presumably have to use POSDAT in order to get the velocity data required to calculate f, and then use VISTUR to deal with the C_mu term and SORKEP to deal with the rest.

Seems like a lot of work for something which is conceptually straightforward, do you think this method is overcomplicated?

Thanks again for your help and suggestions,

Ant
  Reply With Quote

Old   January 24, 2005, 15:56
Default Re: Modelling a turbulent jet and k-epsilon consta
  #4
Richard
Guest
 
Posts: n/a
One thing which you could do to simplify the coding would be to contact star support and ask them for the common block which contains the velocity field. Then you could calculate f locally (and separately) in vistur and sorkep, without relying on posdat. That's probably about as simples as it could be.
  Reply With Quote

Reply

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
Model constants in k-epsilon model Sachin T Main CFD Forum 17 July 18, 2016 14:35
k epsilon model constants for specific application Flavio Main CFD Forum 0 February 5, 2009 04:43
turbulent model k-epsilon yann FLUENT 6 August 7, 2004 08:18
How to adjust turbulent transport constants? Wang FLUENT 0 December 6, 2000 07:57
Wanted buoyancy constants in k-epsilon model grig Main CFD Forum 1 February 5, 2000 18:21


All times are GMT -4. The time now is 07:38.