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

Stability of results and time step ?

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

LinkBack Thread Tools Display Modes
Old   February 26, 2001, 09:10
Default Stability of results and time step ?
mnasri chokri
Posts: n/a
I am simulating equations in spectral variables , the general form of equations is d/dt+2.v.k^2)E(k,t)=T(k,t), the function T(k,t) contains linear and non-linear terms. I have problems with stability, I would like to know which terms may have influence on it, the time step for exemple , or k_min, k_max, k_pic(initial spectrum),the caractristic time, ...

thanks in advance.
  Reply With Quote

Old   February 26, 2001, 13:44
Default Re: Stability of results and time step ?
Posts: n/a
If you are using the EDQNM form of the transfer term, then you should have :

T(K) = int_{P,Q} (a_kpq E(P)E(Q) - b_kpq E(K)E(P)) dPdQ = T_nl - T_l

since the linear term T_l is negative, you can increase the stability of your scheme by writing it in an implicit form. This is possible if you are looking to a steady state solution - by forcing the small wave length for example.

If you are looking to a unsteady solution - turbulence decay for example - you should use a time scale at least as small as the smallest caracteristic time scale of your problem. In that case, I believe it is the turn over time of the smallest eddy - i.e. the turn over time associate to the Kolmogorov wave length K_kol. So if :

K_kol = (epsilon/nu^3)^{1/4}

then :

dt < 1/ ( 2 nu K_kol * K_kol )

where the dissipation rate epsilon is given by :

epsilon = int_{K} 2 nu K * K E(K) dK

hope that help
  Reply With Quote


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
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM Running, Solving & CFD 8 March 28, 2016 11:08
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24
DPM UDF particle position using the macro P_POS(p)[i] dm2747 FLUENT 0 April 17, 2009 01:29

All times are GMT -4. The time now is 15:08.