Python program to create inputs for two dimensional full-wave COMSOL-RF simulations based on Miller flux surfaces. COMSOL simulations require the plasma density (ne) and the magnetic field components (Bx, By, Bz) over an R, Z grid.
Toy models of plasma turbulence can also be included and written to a file for use in COMSOL.
The program can also generate inputs for the 3D Scotty beam-tracing code (ne.dat and topfile.json).
The user describes the plasma with,
- Miller parameters (
R0(r), delta(r), kappa(r)) - Toroidal field function (
I(r) = R*B_tor) - Saftey factor profile (
q(r)) - Density profile (
ne(r))
The user also specifies the rectangular domain of the COMSOL simulation. COMSOL inputs are interpolated onto the rectangular domain with a higher resolution.
The resulting flux surfaces (and density profile) are not necessarily solutions to the Grad-Shafranov equation.
This plot illustrates the flux surfaces and calculated poloidal magnetic field given the profiles input on the right side of the figure. The red rectangle indicates an example COMSOL domain.
This plot illustrates the mapped density profile as well as calculations of the midplane plasma cutoff frequencies. In this case we use a simple piecewise-continuous linear density profile.
This plot shows an example with include_turb=True. We use a toy model density perturbation with k_theta = 0.3 [rad/cm] and Gaussian localization in the r-direction.
