Culham Electron Cyclotron Model | culecd()
iefrf/iefrffix
= 7
This routine calculates the current drive parameters for a electron cyclotron system, based on the AEA FUS 172 model1
- Local electron temperature (\mathtt{tlocal}) is calculated using the
tprofile
method -
Local electron density (\mathtt{dlocal}) is calculated using the
nprofile
method -
Calculate the inverse aspect ratio
epsloc
.
- Calculate the Coulomb logarithm for ion-electron collisions
coulog
.2
Calculate normalised current drive efficiency at four different poloidal angles, and average. cosang = cosine of the poloidal angle at which ECCD takes place = +1 outside, -1 inside.
Normalised current drive efficiency
Uses the eccdef
model found here
-
Calculate the normalised current drive efficiency at four different poloidal angles, and average. - Set
cosang
to 1.0 and calculateecgam1
. - Setcosang
to 0.5 and calculateecgam2
. - Setcosang
to -0.5 and calculateecgam3
. - Setcosang
to -1.0 and calculateecgam4
.cosang = 1.0e0 ecgam1 = self.eccdef(tlocal, epsloc, zlocal, cosang, coulog) cosang = 0.5e0 ecgam2 = self.eccdef(tlocal, epsloc, zlocal, cosang, coulog) cosang = -0.5e0 ecgam3 = self.eccdef(tlocal, epsloc, zlocal, cosang, coulog) cosang = -1.0e0 ecgam4 = self.eccdef(tlocal, epsloc, zlocal, cosang, coulog)
-
Calculate the normalised current drive efficiency
ecgam
as the average ofecgam1
,ecgam2
,ecgam3
, andecgam4
.
- Calculate the current drive efficiency by dividing
ecgam
by(dlocal * physics_variables.rmajor)
.
Note: The eccdef
method is called to calculate the current drive efficiency at each poloidal angle.
-
Hender, T.C., Bevir, M.K., Cox, M., Hastie, R.J., Knight, P.J., Lashmore-Davies, C.N., Lloyd, B., Maddison, G.P., Morris, A.W., O’Brien, M.R. and Turner, M.F., 1992. "Physics assessment for the European reactor study." AEA FUS, 172. ↩
-
Wesson, J. and Campbell, D.J., 2011. "Tokamaks" (Vol. 149). Oxford university press. ↩