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

calculation of k, epsilon and omega

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

Like Tree8Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   January 5, 2013, 05:13
Default
  #21
Member
 
sqing
Join Date: Sep 2012
Location: Dalian
Posts: 76
Rep Power: 4
Sunxing is on a distinguished road
Quote:
Originally Posted by umer.chaudrey View Post
Hey everyone,

I have a query regarding K-epsilon values too.

I am using a transient compressible rhoPimpleFoam solver, RAS turbulence model but without any wall functions for k and epsilon.

Whatever I read in various threads on this forum, I understand that initial values of internal fields of k & epsilon are independent of the solution.

But however in my case, if I use different values my solution behaviour and values change. Can anyone please comment on that?

I had written on another thread too, but no one replied regarding specific k & epsilon intial values for internal field. I am not using a velocity boundary condition, rather pressure difference.

Thanks, I will really appreciate some advice.
Hi umer.chaudrey

I am also simulating a case using the rhoPimpleFoam solver, RAS turbulence model. In my case--wall-jet flow, I set the inlet velocity as fixed value. However, I don't know how to set the the outlet condition for velocity and inlet/outlet condition for pressure. I wander that how do you set these condition in your case?

Regards,
Sunxing
Sunxing is offline   Reply With Quote

Old   January 11, 2013, 13:01
Default
  #22
Member
 
Suranga Dharmarathne
Join Date: Jan 2011
Location: TX, USA
Posts: 39
Rep Power: 5
sdharmar is on a distinguished road
Hi Sunxing,

You can use zeroGradient BC or inletOutlet BC for U at the outlet.
For p you can use zeroGradient at the inlet and fixedvalue at the outlet. I am also a new member so please double check my suggestions.

If any senior members can support/oppose my comment that would be really appreciated.

BR,
Suranga.
sdharmar is offline   Reply With Quote

Old   January 11, 2013, 21:09
Default
  #23
Member
 
sqing
Join Date: Sep 2012
Location: Dalian
Posts: 76
Rep Power: 4
Sunxing is on a distinguished road
Hi Suranga,

Thanks for your reply. I have already tested this condition that you have suggested, however it has turned out to be a bad setting.

regards,
Sunxing
Sunxing is offline   Reply With Quote

Old   May 26, 2013, 03:55
Default
  #24
Member
 
Join Date: Aug 2011
Posts: 63
Rep Power: 5
idefix is on a distinguished road
Hello,

at the moment I have some problems with one simulation.
k and epsilon have very high values at the wall. I think that there is something wrong with the wall function.
But I am not familiar with this topic.
I read this sentence here in the discussion.

Quote:
Yes. If you have a look in the source code, you will see that the wall function for k simply impose a zeroGradient condition at the wall, while the while function for epsilon calculate (for each timestep) the first grid point value by using an algebraic expression derived from the classical logarithmic law-of the wall approach.
Where can I find this mentioned source code?
Are there any literature explaining the wall function model used in OF in k epsilon-model?

Sorry for this simple question, but I am very confused at the moment and I am not very good in C++.

Thanks for your help
Idefix
idefix is offline   Reply With Quote

Old   June 6, 2013, 10:32
Default
  #25
New Member
 
ali
Join Date: Feb 2012
Posts: 8
Rep Power: 4
cfdman_aero is on a distinguished road
hi, I have a difficulty with discretization of source terms of k and epsilon. for example we have (du/dy)**2 in the k source term. my code is on FVM and explicit. Please guide me through.
cfdman_aero is offline   Reply With Quote

Old   June 7, 2013, 03:15
Default
  #26
Member
 
Join Date: Mar 2010
Posts: 77
Rep Power: 6
cutter is on a distinguished road
Please don't hijack other threads! What have you done so far and which discretization scheme are your trying to employ? What's your intention - writing an own solver?

Have a look at the book "Essential Computational Fluid Dynamics" by Oleg Zikanov. It covers the basics of FVM discretization pretty thoroughly. Many other CFD books do as well.

cutter
cutter is offline   Reply With Quote

Old   June 7, 2013, 03:24
Default
  #27
New Member
 
ali
Join Date: Feb 2012
Posts: 8
Rep Power: 4
cfdman_aero is on a distinguished road
Quote:
Originally Posted by cutter View Post
Please don't hijack other threads! What have you done so far and which discretization scheme are your trying to employ? What's your intention - writing an own solver?

Have a look at the book "Essential Computational Fluid Dynamics" by Oleg Zikanov. It covers the basics of FVM discretization pretty thoroughly. Many other CFD books do as well.

cutter
Thanks for reply. I am writing FVM for incompressible turbulent flow utilizing K-Epsilon. The scheme is explicit and I do not know how to discretize the nonlinear production in wich we have i.e. (du/dy)**2. I use some methods but the code doses not run and exit with "floating point error: invalid". I wonder if there is any method to linearized form of those terms.
cfdman_aero is offline   Reply With Quote

Old   August 20, 2013, 13:05
Default regarding boundary conditions for omega
  #28
New Member
 
payal
Join Date: Aug 2013
Location: banglore
Posts: 14
Rep Power: 3
payal05 is on a distinguished road
helo
I am relatively new in OF and would be grateful if some one can help me please. I am trying to set up a case using rhopimplefoam solver and the kOmegaSST model but I have difficulties defining the boundary conditions for omega. frankly speaking i dont have any clue about omega BC"S.

in my case i have inlet , axis , upperwall, outlet so hw to define BC's for that ?
payal05 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
Help with k epsilon values of turbulence vw.cfd OpenFOAM 17 December 12, 2013 05:57
Bounding epsilon or bounding omega Stylianos OpenFOAM 5 June 9, 2011 00:40
boundary conditions for simpleFoam calculation foam_noob OpenFOAM Running, Solving & CFD 7 March 3, 2010 20:17
Use of k omega turbulence model john_w OpenFOAM Running, Solving & CFD 2 September 22, 2009 05:15
Values of k & epsilon in the moving wall ghlee Main CFD Forum 1 August 16, 1999 12:42


All times are GMT -4. The time now is 04:30.