
[Sponsors] 
How to define the turbulence intensity and mixing length at an outlet (for keps)? 

LinkBack  Thread Tools  Display Modes 
January 27, 2010, 14:14 
How to define the turbulence intensity and mixing length at an outlet (for keps)?

#1 
New Member
Join Date: Jan 2010
Posts: 7
Rep Power: 8 
Hello,
I am new in the OpenFOAM world so please apology if my question is simple. I have a model with an inlet defined as total pressure inlet boundary condition and an outlet defined as pressure outlet boundary condition. The aim of the simulation is to determine the massflow in a pipe. For the inlet boundary, I have used the type "turbulentIntensityKineticEnergyInlet" for 0/k and "turbulentMixingLengthDissipationRateInlet" for 0/epsilon in order to define the turbulence quantities at the boundary. Thanks to this boundary condition type, the turbulence quantities depend on the velocity at the inlet (massflow). At the outlet, I have backflow. I already have added an extrusion to limit the backflow but there are still recirculation at the oulet boundary. For the moment, I use a "inletOutlet" to define the turbulence for the backflow but k and epsilon are constant. So, for the outlet, I would like to define 0/k and 0/epsilon in a similar way as at the inlet (by using the turbulent intensity and the mixing length) to make the backflow turbulence depending on the backflowmassflow but for the moment, I have not found how to define the turbulence at the outlet. How should I define 0/k and 0/epsilon for the outlet in order to specify a turbulent intensity and a mixing length for the k and epsilon calculation ? Is it possible in OpenFOAM? It would be nice if someone could help me. Thanks, David 

January 28, 2010, 04:25 

#2  
Senior Member
Mark Olesen
Join Date: Mar 2009
Location: http://olesenm.github.io/
Posts: 780
Rep Power: 19 
Quote:
This does not work directly out of the box. You'd need to roll your own boundary conditions that would look like "inletOutlet", but calculate k/epsilon along the same lines as your inlet condition. But the fundamental question: if you have backflow, how can you know a fixed value for the turbulent intensity and mixing length a priori? 

January 28, 2010, 13:49 

#3 
New Member
Join Date: Jan 2010
Posts: 7
Rep Power: 8 
Thank you Mark for your answer.
I agree with your remark : we cannot know the turbulent intensity and the mixing length a priori. So the question is to know what is the "best" pratice : 1) To use constant values for k and epsilon ? or 2) To use an turbulent intensity and a mixing length to calculate k and espilon depending on the velocity of the backflow ? I will make some tests to investigate the influence of constant values for k and epsilon on the solution results and on the convergence of the simulation. Thanks. David 

September 10, 2010, 02:01 

#4 
New Member
Join Date: Jul 2010
Posts: 20
Rep Power: 8 
Hello David,
have you made any progress on this problem? I guess I'm facing a similar issue. I'm having convergence problems and I suspect that they are due to wron/inappropriate boundary conditions. As U and p should be fine, I'm guessing it has something to do with the turbulent boundary conditions. Thanks, bb 

September 10, 2010, 05:10 

#5 
Senior Member
Nilesh Rane
Join Date: Apr 2010
Posts: 122
Rep Power: 8 
Why do you want to be so specific in giving outlet k and epsilon BCs??? A crude estimate is all what is needed, as the amount of reversed flow is not so significant. A good guess would be: specify k and epsilon as "$internalField" value. This is generally good enough. If you have large reversed flow which affecting your solution then you better rethink on your domain size.
__________________
Imagination is more important than knowledge..


October 13, 2010, 07:01 

#6 
New Member
Join Date: Jul 2010
Posts: 20
Rep Power: 8 
Thanks for your reply. You're right, the workarounds you posted are probably the better ways of dealing with this kind of problem.
For my case, I think the problem was something different, not the turbulent boundary conditions ... now I got everything smoothly running, without changing the domains boundaries. Thanks anyway for your help! 

February 2, 2011, 04:01 

#7  
Member
Join Date: Sep 2010
Posts: 36
Rep Power: 9 
Quote:
I am new to this openfoam. Even i am facing the same situation. I am trying to simulate diesel flow in the nozzle of the injector. At the outlet i defined inletOutlet boundary condition for velocity. I am using cavitatingFoam solver. I am using komega model. for k and omega at the outlet i have defined zerogradient. I have backflow at the outlet. I tried extruding the outlet but its even worse now. So could you tell me what could be the other reason that you had as you mentioned.or how i can improve my solution?? 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
mixing length  lili  CDadapco  1  February 5, 2007 20:22 
Turbulence boundary values  lego  CFX  9  October 25, 2002 11:55 
Turbulence boundary values  lego  Main CFD Forum  0  October 24, 2002 13:47 
Mixing length  Sabatini  CDadapco  3  December 31, 2000 06:08 