pub fn calculate_forces(
psi: &Array2<f64>,
grid: &Grid2D,
r: f64,
z: f64,
i_plasma_ma: f64,
) -> FusionResult<(f64, f64)>Expand description
Calculate forces acting on a current-carrying plasma ring.
Returns (F_radial [N], F_vertical [N]).
F_hoop = (μ₀ I²/2) · [ln(8R/a) + β_p + l_i/2 - 1.5] / R F_Lorentz_R = I · B_Z · 2πR F_Lorentz_Z = -I · B_R · 2πR