CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Main CFD Forum (http://www.cfd-online.com/Forums/main/)
-   -   k-epsilon problem: non-physical turbulent viscosity production in coaxial free-jet (http://www.cfd-online.com/Forums/main/78569-k-epsilon-problem-non-physical-turbulent-viscosity-production-coaxial-free-jet.html)

 Felipeb July 26, 2010 09:45

k-epsilon problem: non-physical turbulent viscosity production in coaxial free-jet

What are some unexpected producers of turbulent viscosity when using k-epsilon?

In our simulation, concentrically aligned free-jets are simulated for comparison with experimental data.

In this figure, turbulent viscosity measured on the domain's centerline is shown in red, the axial velocity measured on the domain's centerline is shown in green, and the experimental axial velocity data which we seek to match is shown in orange. These curves are overlaid on top of the domain which is shown in blue.

http://lh4.ggpht.com/_PqzKw7xCW0A/TE...c0/s800/tv.jpg

Why would k-epsilon produce the large amount of turbulent viscosity starting at about 40% of the length of the domain?

Grid dependency in the direction of flow has been eliminated.

Further simulation details:
1. Fluent
2. 2D Axisymmetric
3. Incompressible
4. Species Transport model Hydrogen-Air
5. Concentric velocity-inlets: inner: Hydrogen, turbulence intensity 4%
6. Concentric velocity-inlets: outer: Air, turbulence intensity of 3%
7. Pressure inlet: Air, 1 atm
8. Pressure outlet: 1 atm
Very grateful for any advice or suggestions on how to remove/remedy this non-physical turbulent viscosity production.

 Felipeb July 26, 2010 17:02

k, epsilon, rho plot

Below are plotted the constituent quantities comprising the turbulent viscosity in our simulation.
The turbulent viscosity is defined as follows:

http://www.cfd-online.com/W/images/m...7849438c8e.png

From initial inspection of the plots, it looks like k is at fault.

Why would k suddenly jump at X/D = 8 ? Any suggestions?
http://lh4.ggpht.com/_PqzKw7xCW0A/TE...lots_bitte.jpg

 All times are GMT -4. The time now is 01:45.