freegs.jtor.ConstrainPaxisIp

class freegs.jtor.ConstrainPaxisIp(eq, paxis, Ip, fvac, alpha_m=1.0, alpha_n=2.0, Raxis=1.0)

Bases: Profile

Constrain pressure on axis and plasma current

__init__(eq, paxis, Ip, fvac, alpha_m=1.0, alpha_n=2.0, Raxis=1.0)

paxis - Pressure at magnetic axis [Pa] Ip - Plasma current [Amps] fvac - Vacuum f = R*Bt

Raxis - R used in p’ and ff’ components

Methods

Jtor(R, Z, psi[, psi_bndry])	Calculate toroidal plasma current
__init__(eq, paxis, Ip, fvac[, alpha_m, ...])	paxis - Pressure at magnetic axis [Pa] Ip - Plasma current [Amps] fvac - Vacuum f = R*Bt
ffprime(pn)	f * df/dpsi as a function of normalised psi.
fpol(psinorm[, out])	Return f as a function of normalised psi
fvac()	Return f = R*Bt in vacuum
pprime(pn)	dp/dpsi as a function of normalised psi.
pressure(psinorm[, out])	Return p as a function of normalised psi by integrating pprime

Jtor(R, Z, psi, psi_bndry=None)

Calculate toroidal plasma current

Jtor = L * (Beta0*R/Raxis + (1-Beta0)*Raxis/R)*jtorshape

where jtorshape is a shape function L and Beta0 are parameters which are set by constraints

ffprime(pn)

f * df/dpsi as a function of normalised psi. 0 outside core. Calculate ffprimeshape inside the core only.

fvac()

Return f = R*Bt in vacuum

pprime(pn)

dp/dpsi as a function of normalised psi. 0 outside core Calculate pprimeshape inside the core only