Comments written by Greg Hammett May 27, 1997, with input from several others...

These are comments of work in progress for discussion among those involved. I am not completely confident that this is all correct...To summarize the following: 1. At least for some initial conditions, the gyrofluid code appears to incorrectly damp the ExB flow to zero, instead of the residual value given by the Rosenbluth-Hinton formula and reproduced by Dimits' gyrokinetic code. However, for other initial conditions, the gyrofluid code does give a significant non-zero residual phi, but the scaling and dependence of this residual level is not yet clear. 2. The gyrofluid equations probably do have too much long-term damping of higher k_r rotation modes, due to an excessive neoclassical viscosity. I'm guessing the resulting damping rate scales like (v_t/R)*(k_r*rho)**2 (give or take factors of q or r/R perhaps). We have assumed that the short-term damping (of order v_t/R) and/or the turbulent viscosity exceeds this level, so that it doesn't matter much (see point 4 below). However, there could be regimes close to marginal stability where the turbulent viscosity isn't large enough to be dominant, and the excessive long-term neoclassical viscosity might reduce the rotation modes and cause us to overestimate chi. We could eventually try to quantify this criterion so we can quantify how close to marginal stability this occurs. 3. There are several approaches to fix these problems which we are looking into. One might be to adjust the existing Landau-fluid closure coefficients in various ways. [There is some freedom in the way these are chosen. The present coefficients were chosen to give a best fit to the toroidal response functions, but there are other possible coefficients which do not change the goodness of fit very much.] Another would be to explicitly add in some terms in the closures so that the Landau-damping causes the parallel heat fluxes to relax to their neoclassical values (a Pfirsch-Schluter heat flux analogous to the Pfirsch-Schluter flow). We are double checking the form of the Landau-damping operator (should |k_par| operate on q_par or on q_par/B, etc.)? Another approach would be to separate the ion gyrofluid eqs. into trapped and passing components and treat them separately (as Waltz has suggested). And another would be to treat the m=n=0 component with a fully-kinetic Vlasov treatment. 4. Even if the gyrofluid equations are getting the long-term behaviour of the m=n=0 rotation mode incorrectly, there are some reasons why the nonlinear chi might not change much. Nevertheless, we are taking this as a serious issue, because there is some chance that this causes the gyrofluid equations to overestimate chi. The planned GF/GK comparisons and additional diagnostics, which we are all pursuing, will be important. There are several ways to look at point 4. Before going on with the details of point 1 above, let me go into more detail on point 4: In the DIA treatment of turbulence, the turbulence generates both "nonlinear noise" which drives the amplitude of various modes and "nonlinear damping" which takes energy out of various modes. (Note that this is a "real" noise, i.e., treating part of the turbulence as a random source, as opposed to "particle noise", which is another topic...) For example, in regular Navier-Stokes turbulence, there is a wide "inertial range" of undamped modes, but their amplitude doesn't grow to infinity because the the nonlinear driving of them is balanced by the nonlinear damping of them. Rosenbluth and Hinton focussed on the turbulent "nonlinear noise" driving the linearly undamped m=n=0 modes, but there should also be nonlinear damping of these modes. Physically, this correponds to a turbulent viscosity which mixes regions of oppositely directed poloidal flow. I believe there is some evidence for this in Dimits' gyrokinetic paper at the 1994 IAEA meeting (Vol. III, p. 457). In his Fig. 2 is a scan of chi vs. r/R. If Rosenbluth and Hinton's mechanism were dominant (that the turbulent viscosity is small and that larger r/R allows larger levels of the rotation mode to build up), then one would expect that chi would drop as r/R is increased from 0 to 0.2. However, Dimits found that in fact chi_hat increased from about 0.75 to 1 as r/R was increased from 0 to 0.2. This would be consistent with our general experience as well, that increasing r/R increases the trapped-ion drive of the turbulence some. This suggests that r/R=0 is a useful GF/GK benchmark point, since then the rotation mode should damp completely to zero, eliminating any concerns about its possible long-term build up. Does the r/R scan from Dimits' 1994 IAEA paper need to be redone with the present version of his GK code before the GF code compares with it?