Plasma simulation

The complexity of fusion-grade plasmas and the increased computational power that has become available in recent years has made the simulation of plasmas a prime object of study in the field of fusion research. Although the basic equations governing the behaviour of magnetised plasmas are known, approximations are always necessary in any code of practical interest; e.g. the extreme disparity of the transport timescales (seconds) and turbulent timescales (microseconds) make it hard to perform detailed turbulence simulations for the whole three-dimensional plasma volume and for several transport timescales.

Codes

Codes can either be interpretative (taking some input from experiment) or predictive. They can be full-tokamak (or full-stellarator), or simulate only a small portion of plasma (a flux tube, the edge, or the Scrape-Off Layer). They can be fluid models for one (electrons), two (electrons + ions) or more (impurities) fluid species, Monte Carlo type (particle tracing) codes, or gyro-kinetic codes. The latter are again subdivided into full-f or delta-f codes (delta-f referring to the fact that only the deviation from a background Maxwellian particle velocity distribution function is simulated).

Recent years have seen an increased effort in the field of cross code benchmarking.

Fluid codes

• CUTIE (predictive, 3-D, full-tokamak)
• PRETOR
• PROCTR (1-D)
• TRANSP
• JETTO
• MMM95
• EDGE2D-NIMBUS (edge)
• UEDGE
• SOLPS

Monte Carlo

• MOCA
• EIRENE (edge)

Gyrokinetic codes

The gyrokinetic treatment simplifies the Vlasov equation for the evolution of the single-particle distribution function $f({\vec {q}},{\vec {p}},t)$ by averaging over the gyration angle, resulting in an evolution equation for the particle guiding centre.

• GYRO
• GS2
• GENE
• GEM (delta f)
• EUTERPE
• SUMMIT/PG3EQ_NC