current_drive_variables Module

Module containing global variables relating to the current drive system

References

-



Contents


Variables

TypeVisibility AttributesNameInitial
real(kind=dp), public :: beamwd

width of neutral beam duct where it passes between the TF coils (m) T Inoue et al, Design of neutral beam system for ITER-FEAT, Fusion Engineering and Design, Volumes 56-57, October 2001, Pages 517-521)

real(kind=dp), public :: bigq

Fusion gain; P_fusion / (P_injection + P_ohmic)

real(kind=dp), public :: bootstrap_current_fraction

bootstrap current fraction (enforced; see i_bootstrap_current)

real(kind=dp), public :: bootstrap_current_fraction_max

maximum fraction of plasma current from bootstrap; if bootstrap_current_fraction_max < 0, bootstrap fraction = abs(bootstrap_current_fraction_max)

real(kind=dp), public :: bscf_iter89

bootstrap current fraction, ITER 1989 model

real(kind=dp), public :: bscf_nevins

bootstrap current fraction, Nevins et al model

real(kind=dp), public :: bscf_sauter

bootstrap current fraction, Sauter et al model

real(kind=dp), public :: bscf_wilson

bootstrap current fraction, Wilson et al model

real(kind=dp), public :: bscf_sakai

Bootstrap current fraction, Sakai et al model

real(kind=dp), public :: cboot

bootstrap current fraction multiplier (i_bootstrap_current=1)

real(kind=dp), public :: cnbeam

neutral beam current (A)

real(kind=dp), public :: diacf_hender

diamagnetic current fraction, Hender fit

real(kind=dp), public :: diacf_scene

diamagnetic current fraction, SCENE fit

real(kind=dp), public :: diamagnetic_current_fraction

diamagnetic current fraction

real(kind=dp), public :: echpwr

ECH power (MW)

real(kind=dp), public :: echwpow

ECH wall plug power (MW)

real(kind=dp), public :: effcd

current drive efficiency (A/W)

real(kind=dp), public :: harnum

cyclotron harmonic frequency number, used in cut-off function

integer, public :: wave_mode

Switch for ECRH wave mode :

  • =0 O-mode
  • =1 X-mode
real(kind=dp), public :: enbeam

neutral beam energy (keV) (iteration variable 19)

real(kind=dp), public :: etacd

auxiliary power wall plug to injector efficiency

real(kind=dp), public :: etacdfix

secondary auxiliary power wall plug to injector efficiency

real(kind=dp), public :: etaech

ECH wall plug to injector efficiency

real(kind=dp), public :: etalh

lower hybrid wall plug to injector efficiency

real(kind=dp), public :: etanbi

neutral beam wall plug to injector efficiency

real(kind=dp), public :: fpion

fraction of beam energy to ions

real(kind=dp), public :: pnbitot

neutral beam power entering vacuum vessel

real(kind=dp), public :: pscf_scene

Pfirsch-Schlüter current fraction, SCENE fit

real(kind=dp), public :: nbshinemw

neutral beam shine-through power

real(kind=dp), public :: feffcd

current drive efficiency fudge factor (iteration variable 47)

real(kind=dp), public :: forbitloss

fraction of neutral beam power lost after ionisation but before thermalisation (orbit loss fraction)

real(kind=dp), public :: frbeam

R_tangential / R_major for neutral beam injection

real(kind=dp), public :: ftritbm

fraction of beam that is tritium

real(kind=dp), public :: gamcd

normalised current drive efficiency (1.0e20 A/(W m^2))

real(kind=dp), public :: gamma_ecrh

User input ECRH gamma (1.0e20 A/(W m^2))

real(kind=dp), public :: xi_ebw

User scaling input for EBW plasma heating. Default 0.43

integer, public :: iefrf

Switch for current drive efficiency model:

  • =1 Fenstermacher Lower Hybrid
  • =2 Ion Cyclotron current drive
  • =3 Fenstermacher ECH
  • =4 Ehst Lower Hybrid
  • =5 ITER Neutral Beam
  • =6 new Culham Lower Hybrid model
  • =7 new Culham ECCD model
  • =8 new Culham Neutral Beam model
  • =9 RFP option removed in PROCESS (issue #508)
  • =10 ECRH user input gamma
  • =11 ECRH "HARE" model (E. Poli, Physics of Plasmas 2019). Removed in #1811.
  • =12 EBW user scaling input. Scaling (S. Freethy)
integer, public :: iefrffix

Switch for 2nd current drive efficiency model:

  • =0 No fixed current drive
  • =1 Fenstermacher Lower Hybrid
  • =2 Ion Cyclotron current drive
  • =3 Fenstermacher ECH
  • =4 Ehst Lower Hybrid
  • =5 ITER Neutral Beam
  • =6 new Culham Lower Hybrid model
  • =7 new Culham ECCD model
  • =8 new Culham Neutral Beam model
  • =9 RFP option removed in PROCESS (issue #508)
  • =10 ECRH user input gamma
  • =11 ECRH "HARE" model (E. Poli, Physics of Plasmas 2019). Removed in #1811.
  • =12 EBW user scaling input. Scaling (S. Freethy)
integer, public :: irfcd

Switch for current drive calculation:

  • =0 turned off
  • =1 turned on
real(kind=dp), public :: nbshinef

neutral beam shine-through fraction

real(kind=dp), public :: nbshield

neutral beam duct shielding thickness (m)

real(kind=dp), public :: pheat

heating power not used for current drive (MW) (iteration variable 11)

real(kind=dp), public :: pheatfix

secondary fixed heating power not used for current drive (MW)

real(kind=dp), public :: pinjalw

maximum allowable value for injected power (MW) (constraint equation 30)

real(kind=dp), public :: pinjemw

auxiliary injected power to electrons (MW)

real(kind=dp), public :: pinjimw

auxiliary injected power to ions (MW)

real(kind=dp), public :: pinjmw

total auxiliary injected power (MW)

real(kind=dp), public :: pinjfixmw

secondary total fixed auxiliary injected power (MW)

real(kind=dp), public :: plasma_current_internal_fraction

plasma current fraction driven internally (Bootstrap + Diamagnetic + PS)

real(kind=dp), public :: plhybd

lower hybrid injection power (MW)

real(kind=dp), public :: pnbeam

neutral beam injection power (MW)

real(kind=dp), public :: porbitlossmw

neutral beam power lost after ionisation but before thermalisation (orbit loss power) (MW)

real(kind=dp), public :: ps_current_fraction

Pfirsch-Schlüter current fraction

real(kind=dp), public :: pwplh

lower hybrid wall plug power (MW)

real(kind=dp), public :: pwpnb

neutral beam wall plug power (MW)

real(kind=dp), public :: rtanbeam

neutral beam centreline tangency radius (m)

real(kind=dp), public :: rtanmax

maximum tangency radius for centreline of beam (m)

real(kind=dp), public :: taubeam

neutral beam e-decay lengths to plasma centre

real(kind=dp), public :: tbeamin

permitted neutral beam e-decay lengths to plasma centre


Subroutines

public subroutine init_current_drive_variables()

Initialise module variables

Arguments

None