freegs.jtor.ConstrainBetapIp

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

Bases: Profile

Constrain poloidal Beta and plasma current

This is the constraint used in YoungMu Jeon arXiv:1503.03135

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

betap - Poloidal beta 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, betap, Ip, fvac[, alpha_m, ...])	betap - Poloidal beta 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