WWW Public <fusion.gat.com> Only All <fusion.gat.com> 

Jump to: navigation, search
Documentation
Install/Build
Geometry
INPUT_profiles
Platform files
GYRO Basics
GYRO Input
GYRO Output
NEO Basics
NEO Input
NEO Output
TGYRO Basics
TGYRO Input
TGYRO Output

Contents

INPUT_profiles

Example

A sample INPUT_profiles file is given here.

Overview

The file INPUT_profiles contains the entire dataset required for specification of experimental profiles. All such profiles are specified on an n_grid_exp-point grid. The information included is sufficient to use the Miller local equilibrium model. The format of each data-block is fixed, but between data blocks any number of comment lines (starting with # in the first column) can be added for convenience. These comment lines are ignored by the parser. Unlike INPUT, no inputs are supplied by default; a crash will result if data is missing. Below, we list the individual scalar and vector data names (which coincide with the GYRO internal names). In these definitions, recall that n_spec is the number of species.

Radial coordinates

The connection between radial variables can be confusing. We introduce two primitive radial variables:

  • the midplane minor radius r\,\! ,
  • the areal coordinate \rho\,\! .

The areal coordinate serves to parameterize the toroidal flux:

\Psi_t = \iint {\mathbf B} \cdot d{\mathbf S} \longrightarrow \frac{\Psi_t}{2\pi} = B_0 \, \frac{\rho(r)^2}{2} \,\!

The approach and definitions above are not strictly standard. However, the general approach of parameterizing the toroidal flux in terms of a new radial coordinate is ubiquitous. One source of ambiguity is the choice of reference magnetic field B_0\,\!. This is normally the vacuum toroidal field, or close to it. However, it is important to keep in mind that the choice is ultimately arbitrary. What is required is that by knowing B_0\,\! and \rho\,\! one can exactly calculate \Psi_t\,\!. Note that in terms of \rho\,\!, the area enclosed by a flux surface is approximately given by A \simeq \pi\rho^2\,\!.

The effective magnetic field strength

The information connecting r\,\! with \rho\,\! is defined through the effective magnetic field strength, B_{\rm unit}\,\!, where

B_{\rm unit}(r) = B_0 \, \frac{\rho \, d\rho}{r \, dr} = \frac{\Psi_t(a)}{\pi} \,\frac{\hat\rho \,d\hat\rho}{r\,dr},

Above, we have introduced the root of the nomalized toroidal flux according to

{\hat\rho}(r) = \frac{\rho(r)}{\rho(a)} = \sqrt{\frac{\Psi_t(r)}{\Psi_t(a)}} .

File Structure

Scalar Data

Scalar data comprise the first three non-comment lines:

  1. n_grid_exp
    • Number of experimental data gridpoints.
    • This is the dimension of the vectors which follow, for example, size(rhogrid_exp(:))=n_grad_exp.
  2. bt_exp
    • On-axis magnetic field, B_0\,\!, in Tesla.
    • The sign is taken to be positive if the direction of B points in the clockwise direction, viewed from above. For example, B_0 \simeq 2.1 in DIII-D and B_0 \simeq -3.0 in JET.
  3. arho_exp
    • \rho(a)\,\! (i.e., \rho\,\! at the separatrix) in units of m\,\!.
    • This will serve to denormalize the vector \hat\rho\,\! defined in the next section.

Vector data, block 1 (1-5)

  1. rhogrid_exp(:)
    • The dimensionless areal variable, \hat\rho = \rho(r)/\rho(a).
    • This should span the region {\hat\rho} \in [0,1] using n_grid_exp values.
    • The grid-spacing can be nonuniform.
  2. rmin_exp(:)
    • The midplane minor radius, r\,\!, in units of m\,\!.
  3. rmaj_exp(:)
    • The flux-surface-center major radius, R_0(r)\,\!, in units of m\,\!.
  4. q_exp(:)
    • The dimensionless safety factor, q\,\!.
  5. kappa_exp(:)
    • The dimensionless plasma elongation, \kappa\,\!.

Vector data, block 2 (6-10)

  1. delta_exp(:)
    • The dimensionless plasma triangularity, \delta\,\!.
  2. tem_exp(n_spec,:)
    • The electron temperature, T_e\,\!, in units of keV.
  3. den_exp(n_spec,:)
    • The electron density, n_e\,\!, in units of 10^{19}/{\rm m}^3\,\!.
  4. z_eff_exp(:)
    • The (dimensionless) effective ion charge, z_{\rm eff}\,\!.
  5. er_exp(:)
    • Radial electric field, in units of kV/m.

Vector data, block 3 (11-15)

  1. flow_mom_exp(:)
    • The convected and conducted angular momentum flow [Nt-m].
    • The rate of change of angular momentum within a flux surface has the units of a torque.
    • The transp variable is VOLINT(MVISC+MCOND).
  2. pow_e_exp(:)
    • The convected and conducted electron energy flow [MW].
    • This is the electron heating power, corrected for the ion-electron energy transfer, the radiated power, and the rate of change of the energy stored by the electrons.
    • The transp variable is VOLINT(PCONV+PCOND)*E-6
  3. pow_i_exp(:)
    • The convected and conducted electron energy flow [MW].
    • This is the ion heating power, corrected for the ion-electron energy transfer, the net charge-exchange loss power, and the rate of change of the energy stored by the ions.
    • The transp variable is VOLINT(PCNVE+PCNDE)*E-6
  4. pow_ei_exp(:)
    • The volume-integral of the electron-ion energy transfer rate [MW]
    • The transp variable is VOLINT(QIE)*E-6
  5. delta1
    • Asymmetry of triangularity (see tech. manual).

Vector data, block 4 (16-20)

  1. flow_beam_exp(:)
    • DEFINE
  2. flow_wall_exp(:)
    • DEFINE
  3. z_mag_axis(:)
    • Flux-surface elevation, in units of m\,\!.
  4. (empty)
    • DEFINE
  5. (empty)
    • DEFINE

Vector data, block 5 (21-25)

  1. den_exp(1,:)
    • Ion density, n_{i1}\,\!, in units of 10^{19}/{\rm m}^3\,\!.
  2. den_exp(2,:)
    • Ion density, n_{i2}\,\!, in units of 10^{19}/{\rm m}^3\,\!.
  3. den_exp(3,:)
    • Ion density, n_{i3}\,\!, in units of 10^{19}/{\rm m}^3\,\!.
  4. den_exp(4,:)
    • Ion density, n_{i4}\,\!, in units of 10^{19}/{\rm m}^3\,\!.
  5. den_exp(5,:)
    • Ion density, n_{i5}\,\!, in units of 10^{19}/{\rm m}^3\,\!.

Vector data, block 6 (26-30)

  1. tem_exp(1,:)
    • Ion temperature, T_{i1}\,\!, in units of keV.
  2. tem_exp(2,:)
    • Ion temperature, T_{i2}\,\!, in units of keV.
  3. tem_exp(3,:)
    • Ion temperature, T_{i3}\,\!, in units of keV.
  4. tem_exp(4,:)
    • Ion temperature, T_{i4}\,\!, in units of keV.
  5. tem_exp(5,:)
    • Ion temperature, T_{i5}\,\!, in units of keV.

Vector data, block 7 (31-35)

  1. vphi_exp(1,:)
    • Ion toroidal velocity, v_{\phi,i1}\,\!, in units of m/s\,\!.
  2. vphi_exp(2,:)
    • Ion toroidal velocity, v_{\phi,i2}\,\!, in units of m/s\,\!.
  3. vphi_exp(3,:)
    • Ion toroidal velocity, v_{\phi,i3}\,\!, in units of m/s\,\!.
  4. vphi_exp(4,:)
    • Ion toroidal velocity, v_{\phi,i4}\,\!, in units of m/s\,\!.
  5. vphi_exp(5,:)
    • Ion toroidal velocity, v_{\phi,i5}\,\!, in units of m/s\,\!.

Vector data, block 8 (36-40)

  1. vpol_exp(1,:)
    • Ion poloidal velocity, v_{\theta,i1}\,\!, in units of m/s\,\!.
  2. vpol_exp(2,:)
    • Ion poloidal velocity, v_{\theta,i2}\,\!, in units of m/s\,\!.
  3. vpol_exp(3,:)
    • Ion poloidal velocity, v_{\theta,i3}\,\!, in units of m/s\,\!.
  4. vpol_exp(4,:)
    • Ion poloidal velocity, v_{\theta,i4}\,\!, in units of m/s\,\!.
  5. vpol_exp(5,:)
    • Ion toroidal velocity, v_{\theta,i5}\,\!, in units of m/s\,\!.
Retrieved from "http://fusion.gat.com/theory/INPUT_profiles"