CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

Simple test case for k-epilon model

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

LinkBack Thread Tools Search this Thread Display Modes
Old   September 27, 2012, 18:05
Default Simple test case for k-epilon model
New Member
Join Date: Sep 2011
Location: Germany
Posts: 7
Rep Power: 14
mathletic is on a distinguished road
I wrote a CFD code using the k-epsilon turbulence model, currently only with Dirichlet and outflow boundary conditions. I tested it with the lid-driven cavity flow but I get false solution with k and epsilon dropping to zero everywhere except for the Dirichlet boundary.

So I am looking for a simple test case that helps me to find my mistake. Preferably it would be 1D, so I can examine the complete stiffness matrix. And the example gets extra points, if an analytic solution exists. Does anyone know such a test case or a paper / a book describing it?

I scanned many CFD books and found only the "Grid turbulance" example in [1]. It is 1D, steady state and give analytic solutions for k and epsilon. But still the informations given seems incomplete. Left are Dirichlet boundary conditions, what do I use on the right? What are the solutions for velocity and pressure? What is the turbulent viscosity nu_t and how large is the production term? Where is the example from, why is no source given? How is the analytic solution obtained? Has someone done the calculation for this example?

The lack of good test cases for the k-epsilon model with reference data is discouraging. Everybody seems happy with the back-ward facing step...

[1] Mohammadi, Pironneau: "Analysis of the K-Epsilon Turbulence Model", V.5.1, page 74
mathletic is offline   Reply With Quote

Old   September 28, 2012, 02:12
Senior Member
Oliver Gloth
Join Date: Mar 2009
Location: Todtnau, Germany
Posts: 121
Rep Power: 17
ogloth is on a distinguished road
Good morning,

Dirichlet boundary conditions alone will not enable you to compute wall bounded flows, because for the standard k-epsilon model you need to implement a wall function. You could, however, test your code with free shear flows (e.g. jets). I believe Wilcox's famous book on "Turbulence Modelling for CFD" should have a few suitable test cases.

Once you have a decent wall boundary condition you can solve a flat plate boundary layer.

Since I am not at all a turbulence expert -- maybe someone else has some better proposals.

ogloth is offline   Reply With Quote

Old   September 28, 2012, 03:04
New Member
Join Date: Sep 2011
Location: Germany
Posts: 7
Rep Power: 14
mathletic is on a distinguished road
Godd morning Oliver,
thanks for your quick reply. I hesitate to implement wall functions and add complexity to my program because I want to fix my bug-ridden k-espilon code first.

Wolcox's book has three 2d examples but they are described only vague and I'd like to have a complete description - idiot proof :-)

I am still hoping to get a 1D example or more details about the example from Mohammadi and Pironneau.
mathletic is offline   Reply With Quote


analytic solution, k-epsilon, test case

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
Planing Hull Test Case lava12005 Main CFD Forum 7 August 9, 2013 16:15
What model of turbulence choose to study an external aerodynamics case raffale OpenFOAM 0 August 23, 2012 05:45
[stressAnalysis] Problems with simple test case matthiass OpenFOAM CC Toolkits for Fluid-Structure Interaction 19 April 20, 2012 12:44
Interfoam Droplet under shear test case adona058 OpenFOAM Running, Solving & CFD 3 May 3, 2010 18:46
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24

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