Culham Neutral Beam Model | culnbi()
iefrf/iefrffix
= 8
This routine calculates Neutral Beam current drive parameters using the corrections outlined in AEA FUS 172 to the ITER method. The result cannot be guaranteed for devices with aspect ratios far from that of ITER (approx. 2.8).
Output | Description |
---|---|
\mathtt{effnbss} | Neutral beam current drive efficiency in Amperes per Watt |
\mathtt{fpion} | Fraction of NB power given to ions |
\mathtt{fshine} | Shine-through fraction of the beam |
Where R_{\text{tan}} is major radius at which the centre-line of the beam is tangential to the toroidal direction. This can be user defined
A quick sanity check is done to make sure no negative roots are formed when calculating \mathtt{dpath} this prevents setups where the NBI beam would miss the plasma
Beams topping cross section is calculated via \mathtt{sigbeam} found here. This produces \mathtt{sigstop}
Calculate number of decay lengths to centre
Calculate the shine through fraction of the beam
Deuterium and tritium beam densities
Power split to the ions and electrons is clauclated with the \mathtt{cfnbi()} method found here and outputs \mathtt{fpion}
Current drive efficiency | etanb2()
This routine calculates the current drive efficiency in A/W of a neutral beam system, based on the 1990 ITER model, plus correction terms outlined in Culham Report AEA FUS 172.
Input | Description |
---|---|
\mathtt{abeam} | beam ion mass (amu) |
\mathtt{alphan}, \alpha_n | density profile factor |
\mathtt{alphat}, \alpha_T | temperature profile factor |
\mathtt{aspect}, A | aspect ratio |
\mathtt{dene}, n_{\text{e}} | volume averaged electron density (\text{m}^{-3}) |
\mathtt{dnla}, n_{\text{e,0}} | line averaged electron density (\text{m}^{-3}) |
\mathtt{beam_energy} | neutral beam energy (\text{keV}) |
\mathtt{frbeam} | R_tangent / R_major for neutral beam injection |
\mathtt{fshine} | shine-through fraction of beam |
\mathtt{rmajor}, R | plasma major radius (\text{m}) |
\mathtt{rminor}, a | plasma minor radius (\text{m}) |
\mathtt{ten} | density weighted average electron temperature (\text{keV}) |
\mathtt{zeff}, Z_{\text{eff}} | plasma effective charge |
Charge of beam ions $$ \mathtt{zbeam} = 1.0 $$
Fitting factor (IPDG89)
Volume averaged electron density (10^{20} \text{m}^{-3})
Line averaged electron density (10^{20} \text{m}^{-3})
Critical energy (\text{MeV}) (power to electrons = power to ions) (IPDG89) N.B. ten is in keV
Beam energy in MeV
x and y coefficients of function J0(x,y) (IPDG89)
Fitting function J0(x,y)
Effective inverse aspect ratio, with a limit on its maximum value
Reduction in the reverse electron current due to neoclassical effects
Reduction in the net beam driven current due to the reverse electron current
Normalisation to allow results to be valid for non-ITER plasma size and density:
Line averaged electron density (10^{20} \text{m}^{-3}) normalised to ITER
Distance along beam to plasma centre
Distance along beam to plasma centre for ITER assuming a tangency radius equal to the major radius
Normalisation to beam energy (assumes a simplified formula for the beam stopping cross-section)
A_bd fitting coefficient, after normalisation with ebnorm
Normalised current drive efficiency (\text{A/W} \text{m}^{2}) (IPDG89)
Current drive efficiency (A/W)