Normalisations

Introduction

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.

General Notes

  • All perturbed quantities are scaled up by , the gyrokinetic epsilon
  • Unless otherwise stated, the subscript 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.
  • The flux surface label should not be confused with the species gyroradii:

Definitions

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.

  • - All lengths are normalized to the quantity , which is chosen to be half the diameter of the last closed flux surface, measured on the midplane (the level of the magnetic axis). You can relate a to using the input parameter epsl, which is equal to
  • - is the average of the maximum and minimum of R on the flux surface. With circular flux surfaces this is obviously the major radius at the center of the flux surface. GS2 users often use and interchangeably.
  • - 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.

Coordinates

See main article GS2 Coordinates

Flux Surface Label

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.

Lengths

All lengths are normalized to the quantity a.

Time

Times are normalised to

Geometrical Factors

Safety Factor

Magnetic Shear

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,

Miscellaneous Parameters

Flow Shear

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

Fields

Magnetic Field

The equilibrium field is defined as:

The magnetic field is normalized to .

Magnetic Flux

The dimensions of are length * length * magnetic field so

The Perturbed Field

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

The Electric Potential

See #The Perturbed Field.

Fluxes

Heat Flux

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:

Toroidal Angular Momentum Flux

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:

Collisions

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.

Collision Frequency

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.

GS2/AstroGK

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 .

Abel et al.

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.

Rosenbluth Hazeltine & Hinton

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.

Hinton & Wong

This is taken from Braginskii.

There is a factor of between this and the GS2 frequency.