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

Problem with Nonlinear k-epsilon model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 11, 2003, 08:15
Default Problem with Nonlinear k-epsilon model
  #1
bylin
Guest
 
Posts: n/a
I'm trying to compute a channel flow with nonlinear k-eplison model ,which is added as a source term to the momentum equations of the standard k-epsilon model. The source term involves derivatives of Reynolds stress,and which are derived by the nonlinear k-epsilon model.I used Gauss theorum to calculate the derivatives,and assumed zero Reynolds stress at the walls and Neumann boundary condition at the inlet and outlet. I started with standard k-epsilon model,but when I switched to the nonlinear model,divergence occurred to the x-monmentum equation immediately. What's the problem?
  Reply With Quote

Old   November 11, 2003, 10:19
Default Re: Problem with Nonlinear k-epsilon model
  #2
ap
Guest
 
Posts: n/a
Try to reduce the under relazation factor of the momentum equation and of k and eps equations. Try values around 0.4.

If this doesn't work, check boundary conditions.

Hi

ap
  Reply With Quote

Old   November 11, 2003, 21:43
Default Re: Problem with Nonlinear k-epsilon model
  #3
bylin
Guest
 
Posts: n/a
I used an outflow profile of a sufficiently long channel flow computed by standard k-epsilon model as the inlet boundary condition,and "outflow" as the outlet boundary condition. The same problem occurred even when the relaxation factor is reduced to 0.01.
  Reply With Quote

Old   November 23, 2003, 13:28
Default Re: Problem with Nonlinear k-epsilon model
  #4
ap
Guest
 
Posts: n/a
Check the source term implementation.

Be sure to explicitly define the term dS in the DEFINE_SOURCE macro, which is used by FLUENT to stabilize the solution.

Try to iterate using first order discretization scheme for turbulent variables, at least for the first iterations.

P.S. An underrelaxation factor of 0.01 is too low, in my opinion. You risk to obtain unphysical results using such a low value. If divergence problem persist even with low values of the underrelaxation factor, usually there is some problem in the implementation of the source term or in boundary conditions.

Hi and sorry for my late answer

ap
  Reply With Quote

Old   November 23, 2003, 22:45
Default Re: Problem with Nonlinear k-epsilon model
  #5
bylin
Guest
 
Posts: n/a
I'm sorry to say that my problem was due to my carelessness in programming the source code.I did not not notice an extremly small number in some denominators. Anyway,thank you so much!I have learned a lot from you. Now the problem is that the results oscillate wildly. Another question:what's the accuracy of Gauss theorum for calculating derivatives?I think it's 2nd order accuracy.Is it right?
  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
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 13 November 4, 2013 15:13
Polyflow Simulation with ViscoElastic Model memory problem sid2909 FLUENT 0 May 10, 2010 05:29
Polyflow Simulation with ViscoElastic Model memory problem sid2909 ANSYS 0 May 10, 2010 05:27
Problem with LaunderSharma compressible turbulence model wrong formulation luca OpenFOAM Running, Solving & CFD 4 July 3, 2006 06:24
How to define Epsilon in RSM model? David FLUENT 0 March 17, 2004 10:39


All times are GMT -4. The time now is 05:44.