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
tprofilemethod -
Local electron density (\mathtt{dlocal}) is calculated using the
nprofilemethod -
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
cosangto 1.0 and calculateecgam1. - Setcosangto 0.5 and calculateecgam2. - Setcosangto -0.5 and calculateecgam3. - Setcosangto -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
ecgamas the average ofecgam1,ecgam2,ecgam3, andecgam4.
- Calculate the current drive efficiency by dividing
ecgamby(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. ↩