The perennial question: Given a normalised quantity in gs2, what factor do I need to multiply it by to get a real quantity?
The answer should be answered somewhat by individual variable pages, but this page gives an overview of the gs2 normalisation philosophy. Be sure to read the section Definitions, containing definitions of some key normalizing quantities, and also the General Notes, before looking at the rest of the page. Once you have read these sections, each definition of a normalization should be self-contained.
An extensive discussion of normalisations is in Chapter 3 of this thesis.
See also Chapter 26 of the AstroGK manual. Beware: AstroGK does not have variation of the magnetic field magnetic along the field line, but GS2 does. Important differences exist as a result.
There is also a good paper on G. Hammett's website introducing the GS2 geometry module, from which a lot of this page is taken. However, be aware that there is a crucial sign error in the definition of the equilibrium magnetic field (which also propagates into other work). Use the form given on this page.
s
is used to denote the species of interest,
and the subscript 1
is used to denote the reference species, species 1. Note that
species 1 in the GS2 input file does not necessarily set the reference mass, density,
etc. These reference quantities are determined by the definition of mass, density,
etc. in GS2 to be the ratio of the species quantity to the reference quantity. So, if
dens for a given species is set to 1, then the reference density is equal to the
density for that species. The reference mass, temperature, etc. are defined in the
same way.Here we present the definitions of some key normalising quantities. Note that many of these quantities are defined on a given flux surface: this means normalizations are in general dependent on the flux surface label.
- the toroidal field on the given flux surface evaluated at :
- the gyrofrequency of the reference species on the given flux surface evaluated at :
- the gyroradius of the reference species on the given flux surface evaluated at :
- the thermal velocity of species 1 (usually ions). Note that the definition of velocity changes according to the value of norm_option.
See main article GS2 Coordinates
The exact definition of the flux surface label depends on the value of irho. However, irho is almost always set to 2, in which case
where r is the half-diameter of the current flux surface.
All lengths are normalized to the quantity a.
Times are normalised to
Note that if irho = 2 (the default), then this is simply
where is the half-diameter of the flux surface.
In the high aspect ratio, zero circular flux surface limit,
The flow shear is defined as
where is the flux surface label and is the toroidal angular velocity. For historical, numerical and ordering reasons, in GS2 flow shear is specified rather than the more intuitive (when higher flow terms such as Coriolis and centrifugal drifts were all absent, flow shear was the only term that appeared in the equation).
If irho = 2, then and
Now the flow shear has units of inverse time, so it is normalised to
So
The equilibrium field is defined as:
The magnetic field is normalized to .
The dimensions of are length * length * magnetic field so
The perturbed field is normalised to and is also scaled up by , i.e. gyrokinetic epsilon. The perturbed generalised electromagnetic potential is defined as:
and with the normalizations we find that
See #The Perturbed Field.
The radial heat flux for species s is defined (Abel et al., forthcoming) as
In GS2 it is output in spectral components as well as in aggregate.
It is normalised to the GyroBohm diffusion estimate (Cowley et al, 1991):
Replacing by , the pressure of the reference species, species 1, and replacing by , dropping the temperature ratio and scaling up by as is normal for perturbed quantities, we arrive at the normalisation:
The radial component of the toroidal angular momentum flux for species is defined (Abel et al., forthcoming) as
In GS2 it is output in spectral components as well as in aggregate.
It is normalised in an equivalent way to the heat flux:
The implementation of collisions in GS2/AstroGK is described in these two papers by Abel & Barnes:
Abel et. al - Theory
Barnes et. al - Numerical implementation.
All collision frequencies are usually defined in terms of the like particle collision frequency .
The like particle frequency for each species is specified by the input parameter vnewk_s
where is the species index.
There are many definitions of the like particle frequency floating around. Here we display some of them with comments.
It is very important to note that choosing a collision frequency specifies the temperature and density you are considering. In a flux tube calculation, temperature and density are scaled out, except in the collision frequency.
The electron and ion collision frequencies used in GS2 and AstroGK are:
vnewk_e =
vnewk_i =
Note that in the second equation where you might expect to see you see . These formulate were taken from a note of Greg Hammett's on the collision frequency, assuming that here is different from the symbols defined earlier on this page (which depend on the choice of \"norm_option\") and is instead equivalent to as used in Hammett\'s notes. Note that changing the choice of the \"norm_option\" variable (above) changes the normalizations of most of the inputs and outputs involving a timescale by a factor of , but I think that the definition of vnewk_e and vnewk_i are unchanged because they are the collision frequencies normalized by GS2\'s internal time scale .
The collision frequency in Abel et al, 2009, which presents the theory for the GS2 collision operator, is also the same as that of Trebnetov, and also Rosenbluth, Macdonald & Judd
Thus we see that there is no difference between that of Abel and that in GS2/AstroGK, except the normalisation and the treatment of species charge and density.
The collision frequency in Rosenbluth Hazeltine & Hinton 1972 which derive neoclassical transport coefficients is
.
There is a factor of between this and the GS2 frequency.
This is taken from Braginskii.
There is a factor of between this and the GS2 frequency.