suzuki.c File Reference

Suzuki beam-stopping coefficients. More...

#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include "ascot5.h"
#include "consts.h"
#include "error.h"
#include "suzuki.h"

Go to the source code of this file.

Functions
a5err	suzuki_sigmav (real *sigmav, real EperAmu, real vnorm, real ne, real te, integer nion, real *ni, const int *anum, const int *znum)
	Calculate beam-stopping cross-section according to Suzuki model.

Variables
real	A_highE [3][10]
	Fitting parameters for Aijk in equation 28.
real	A_lowE [3][10]
	Fitting parameters for Aijk in equation 28.
integer	Z_imp [] = { 2, 6, 6, 4, 8, 7, 3, 5, 26 }
	Charge number corresponding to the tabulated values of Bijk.
real	Zeffmin_imp [] = { 1.0, 1.0, 5.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 }
	Minimum valid Zeff corresponding to the tabulated values of Bijk.
real	Zeffmax_imp [] = { 2.1, 5.0, 6.0, 4.0, 5.0, 5.0, 3.0, 5.0, 5.0 }
	Maximum valid Zeff corresponding to the tabulated values of Bijk.
real	B_highE [9][12]
	Fitting parameters for Bijk in equation 28.
real	B_lowE [9][12]
	Fitting parameters for Bijk in equation 28.

Detailed Description

Suzuki beam-stopping coefficients.

The data and model is from S Suzuki et al 1998 Plasma Phys. Control. Fusion 40 2097 doi.org/10.1088/0741-3335/40/12/009

The model calculates beam-stopping coefficients for hydrogen beams in hydrogen plasma with impurities present.

Available impurities are: He, Li, Be, B, C, N, O, Fe.

Definition in file suzuki.c.

Function Documentation

suzuki_sigmav()

a5err suzuki_sigmav	(	real *	sigmav,
		real	EperAmu,
		real	vnorm,
		real	ne,
		real	te,
		integer	nion,
		real *	ni,
		const int *	anum,
		const int *	znum )

Calculate beam-stopping cross-section according to Suzuki model.

Parameters

sigmav	evaluated beam stopping cross section [m^2]
EperAmu	test particle energy divided by its atomic mass number [J]
vnorm	test particle velocity [m/s]
ne	electron particle density [m^-3]
te	electron temperature [J]
nion	number of ion species present in the plasma
ni	densities of the ion species with at least single hydrogen species
anum	background ion species' atomic mass number
znum	background ion species' charge number

Returns

zero if evaluation was succesfull

Definition at line 139 of file suzuki.c.

Variable Documentation

A_highE

real A_highE[3][10]

Initial value:

= {
{12.7, 1.25, 0.452, 0.0105, 0.547, -0.102,   0.360, -0.0298, -0.0959, 4.21e-3},
{14.1, 1.11, 0.408, 0.0105, 0.547, -0.0403,  0.345, -0.0288, -0.0971, 4.74e-3},
{12.7, 1.26, 0.449, 0.0105, 0.547, -0.00577, 0.336, -0.0282, -0.0974, 4.87e-3}
}

Fitting parameters for Aijk in equation 28.

Table 2a from Suzuki's paper for H, D, and T. Note the last column in the paper is for different magnetic field magnitude, whose effect was negligible and so that is not included here.

These values are valid in range 100 < E (keV/amu) < 10000.

Definition at line 31 of file suzuki.c.

A_lowE

real A_lowE[3][10]

Initial value:

= {
{-52.9, -1.36, 0.0719, 0.0137, 0.454, 0.403, -0.220, 0.0666, -0.0677, -1.48e-3},
{-67.9, -1.22, 0.0814, 0.0139, 0.454, 0.465, -0.273, 0.0751, -0.0630, -5.08e-4},
{-74.2, -1.18, 0.0843, 0.0139, 0.453, 0.491, -0.294, 0.0788, -0.0612, -1.85e-4}
}

Fitting parameters for Aijk in equation 28.

Table 3a from Suzuki's paper for H, D, and T. Note the last column in the paper is for different magnetic field magnitude, whose effect was negligible and so that is not included here.

These values are valid in range 10 < E (keV/amu) < 100.

Definition at line 46 of file suzuki.c.

Z_imp

integer Z_imp[] = { 2, 6, 6, 4, 8, 7, 3, 5, 26 }

Charge number corresponding to the tabulated values of Bijk.

Definition at line 55 of file suzuki.c.

Zeffmin_imp

real Zeffmin_imp[] = { 1.0, 1.0, 5.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 }

Minimum valid Zeff corresponding to the tabulated values of Bijk.

Definition at line 60 of file suzuki.c.

Zeffmax_imp

real Zeffmax_imp[] = { 2.1, 5.0, 6.0, 4.0, 5.0, 5.0, 3.0, 5.0, 5.0 }

Maximum valid Zeff corresponding to the tabulated values of Bijk.

Definition at line 65 of file suzuki.c.

B_highE

real B_highE[9][12]

Initial value:

= {
    {  0.231,     0.343,    -0.185,    -0.162e-1,  0.105,    -0.703e-1,
       0.531e-1,  0.342e-2, -0.838e-2,  0.415e-2, -0.335e-2, -0.221e-3 },
    { -0.101e1,  -0.865e-2, -0.124,    -0.145e-1,  0.391,     0.161e-1,
       0.298e-1,  0.332e-2, -0.248e-1, -0.104e-2, -0.152e-2, -0.189e-3 },
    { -0.100e1,  -0.255e-1, -0.125,    -0.142e-1,  0.388,     0.206e-1,
       0.297e-1,  0.326e-2, -0.246e-1, -0.131e-2, -0.148e-2, -0.180e-3 },
    { -0.613,     0.552e-1, -0.167,    -0.159e-1,  0.304,     0.154e-2,
       0.436e-1,  0.378e-2, -0.201e-1, -0.216e-3, -0.251e-2, -0.227e-3 },
    { -0.102e1,  -0.148e-1, -0.674e-1, -0.917e-2,  0.359,     0.143e-1,
       0.139e-1,  0.184e-2, -0.209e-1, -0.732e-3, -0.502e-3, -0.949e-4 },
    { -0.102e1,  -0.139e-1, -0.979e-1, -0.117e-1,  0.375,     0.156e-1,
       0.224e-1,  0.254e-2, -0.226e-1, -0.889e-3, -0.104e-2, -0.139e-3 },
    { -0.441,     0.129,    -0.170,    -0.162e-1,  0.277,    -0.156e-1,
       0.466e-1,  0.379e-2, -0.193e-1,  0.753e-3, -0.286e-2, -0.239e-3 },
    { -0.732,     0.183e-1, -0.155,    -0.172e-1,  0.321,     0.946e-2,
       0.397e-1,  0.420e-2, -0.204e-1, -0.619e-3, -0.224e-2, -0.254e-3 },
    { -0.820,    -0.636e-2,  0.542e-1,  0.395e-2,  0.202,     0.806e-3,
      -0.200e-2, -0.178e-2, -0.610e-2,  0.651e-3,  0.175e-2,  0.146e-3 }
}

Fitting parameters for Bijk in equation 28.

Table 2bc (partially) from Suzuki's paper for H, D, and T.

These values are valid in range 100 < E (keV/amu) < 10000.

Definition at line 74 of file suzuki.c.

B_lowE

real B_lowE[9][12]

Initial value:

= {
    { -0.792,     0.420e-1,  0.530e-1,  -0.139e-1,  0.301,    -0.264e-1,
      -0.299e-1,  0.607e-2,  0.272e-3,   0.611e-2,  0.347e-2, -0.919e-3 },
    {  0.161,     0.598e-1, -0.336e-2,  -0.426e-2, -0.157,    -0.396e-1,
       0.460e-2,  0.219e-2,  0.391e-1,   0.711e-2, -0.144e-2, -0.385e-3 },
    {  0.158,     0.554e-1, -0.431e-2,  -0.335e-2, -0.155,    -0.374e-1,
       0.537e-2,  0.174e-2,  0.388e-1,   0.683e-2, -0.160e-2, -0.322e-3 },
    {  0.112,     0.495e-1,  0.116e-1,  -0.286e-2, -0.149,    -0.331e-1,
      -0.426e-2,  0.980e-3,  0.447e-1,   0.652e-2, -0.356e-3, -0.203e-3 },
    {  0.111,     0.541e-1, -0.346e-3,  -0.368e-2, -0.108,    -0.347e-1,
       0.193e-2,  0.181e-2,  0.280e-1,   0.604e-2, -0.841e-3, -0.317e-3 },
    {  0.139,     0.606e-1, -0.306e-2,  -0.455e-2, -0.133,    -0.394e-1,
       0.399e-2,  0.236e-2,  0.335e-1,   0.690e-2, -0.124e-2, -0.405e-3 },
    {  0.112,     0.495e-1,  0.116e-1,  -0.286e-2, -0.149,    -0.331e-1,
      -0.426e-2,  0.980e-3,  0.447e-1,   0.652e-2, -0.356e-3, -0.203e-3 },
    {  0.122,     0.527e-1, -0.430e-3,  -0.318e-2, -0.151,    -0.364e-1,
       0.343e-2,  0.151e-2,  0.420e-1,   0.692e-2, -0.141e-2, -0.290e-3 },
    { -0.110e-1,  0.202e-1,  0.946e-3,  -0.409e-2, -0.666e-2, -0.117e-1,
      -0.236e-3,  0.202e-2,  0.408e-2,   0.185e-2, -0.648e-4, -0.313e-3 }
}

Fitting parameters for Bijk in equation 28.

Table 3bc (partially) from Suzuki's paper for H, D, and T.

These values are valid in range 10 < E (keV/amu) < 100.

Definition at line 102 of file suzuki.c.