|
[Sponsors] |
February 26, 2001, 08:10 |
Stability of results and time step ?
|
#1 |
Guest
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. |
|
February 26, 2001, 12:44 |
Re: Stability of results and time step ?
|
#2 |
Guest
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 |
|
|
|
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 05: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 |