|
|
|
[[_TOC_]]
|
|
|
|
|
|
|
|
# Introduction
|
|
|
|
|
|
|
|
The beam tracing code [GRAY](https://dx.doi.org/10.13182/fst07-a1494) performs the computation of the quasi-optical propagation of a Gaussian beam of electron cyclotron waves in a general tokamak equilibrium, and of the power absorption and driven current.
|
|
|
|
|
|
|
|
The propagation of a general astigmatic Gaussian beam is described within the framework of the complex eikonal approach in terms of a set of "extended" rays that allow for diffraction effects.
|
|
|
|
|
|
|
|
The absorbed power and the driven current density are computed along each ray solving the fully relativistic dispersion relation for electron cyclotron wave and by means of the neoclassical response
|
|
|
|
function for the current.
|
|
|
|
|
|
|
|
The equations and models used in GRAY to compute Gaussian beam propagation, absorption and current drive are summarized in the [Physics](Physics-model) section.
|
|
|
|
The absorbed power and the driven current density are computed along each ray solving the fully relativistic dispersion relation for electron cyclotron wave and by means of the neoclassical response function for the current.
|
|
|
|
|
|
|
|
The main development version of the source code is currently managed under a Git [repository](https://maxwell.ydns.eu/git/istp/gray), and it is available as a stand-alone executable or as a library to be called from within the [JINTRAC](https://git.ccfe.ac.uk/jintrac/jetto-sanco) suite.
|
|
|
|
|
|
|
|
For such version of GRAY, we provide here:
|
|
|
|
- Instructions to [obtain and install](Installation) the code
|
|
|
|
- Description of the [interface](Interface)
|
|
|
|
- Description of the [code parameters](Code-parameters), [input](Input-data) and [output](Output-data) data.
|
|
|
|
|
|
|
|
Two forks of the code happened in the past, in the process of developing the interface for the EUROfusion [Integrated Modelling](https://wpcd-workflows.github.io) workflows, and they are currently maintained in different repositories:
|
|
|
|
|
|
|
|
- [IMAS version at EUROfusion](https://gforge-next.eufus.eu/svn/gray)
|
|
|
|
- [IMAS version at ITER](https://git.iter.org/projects/HEAT/repos/gray/) |
|
|
\ No newline at end of file |
|
|
|
- [IMAS version at ITER](https://git.iter.org/projects/HEAT/repos/gray/)
|
|
|
|
|
|
|
|
# Coordinate systems and conventions
|
|
|
|
|
|
|
|
## Coordinate systems
|
|
|
|
|
|
|
|
## Axisymmetric magnetic field
|
|
|
|
|
|
|
|
## Flux labels
|
|
|
|
|
|
|
|
## Polarisation
|
|
|
|
|
|
|
|
# Physics models
|
|
|
|
|
|
|
|
## Beam propagation
|
|
|
|
|
|
|
|
### Quasi-optical approximation
|
|
|
|
|
|
|
|
### Beam-tracing equations
|
|
|
|
|
|
|
|
### Ray initial conditions
|
|
|
|
|
|
|
|
### Launch point and wavevector
|
|
|
|
|
|
|
|
## Wave absorption
|
|
|
|
|
|
|
|
### Weakly relativistic model
|
|
|
|
|
|
|
|
### Fully relativistic model
|
|
|
|
|
|
|
|
## Current drive
|
|
|
|
|
|
|
|
### High-speed-limit model (Lin-liu)
|
|
|
|
|
|
|
|
### Momentum-conserving model (Marushchenko)
|
|
|
|
|
|
|
|
## Multi-pass absorption
|
|
|
|
|
|
|
|
## Appendix: Gaussian beams
|
|
|
|
|
|
|
|
# Installation
|
|
|
|
|
|
|
|
# Code interface
|
|
|
|
|
|
|
|
## Running GRAY
|
|
|
|
|
|
|
|
## Code parameters (2, 602)
|
|
|
|
|
|
|
|
## Input data
|
|
|
|
|
|
|
|
### Beam characteristics (97, 603)
|
|
|
|
|
|
|
|
### Kinetic profiles (98, 644)
|
|
|
|
|
|
|
|
### Magnetic equilibrium (99, 644)
|
|
|
|
|
|
|
|
### First wall shape (99, 644)
|
|
|
|
|
|
|
|
## Output data
|
|
|
|
|
|
|
|
### Global quantities (7, 607)
|
|
|
|
|
|
|
|
### Profiles in plasma (48, 648)
|
|
|
|
|
|
|
|
### Beam path (4-33, 604-633)
|
|
|
|
|
|
|
|
### Beam shape (8-12, 608-612)
|
|
|
|
|
|
|
|
### (DEBUG) Dispersion relation (17, 617)
|
|
|
|
|
|
|
|
### (DEBUG) Input profiles (55, 645)
|
|
|
|
|
|
|
|
### (DEBUG) Flux averages (56, 646)
|
|
|
|
|
|
|
|
### (DEBUG) Resonances (70, 630)
|
|
|
|
|
|
|
|
### (DEBUG) Flux surfaces (71, 631)
|
|
|
|
|
|
|
|
### (DEBUG) Poloidal plane maps (72, 632) |
|
|
\ No newline at end of file |
