Struct FusionKernel 

pub struct FusionKernel { /* private fields */ }

The Grad-Shafranov equilibrium solver.

Implementations

impl FusionKernel

pub fn new(config: ReactorConfig) -> Self

Create a new kernel from configuration.

pub fn from_file(path: &str) -> FusionResult<Self>

Create a new kernel from a JSON config file.

pub fn solve_equilibrium(&mut self) -> FusionResult<EquilibriumResult>

Main equilibrium solver. Port of solve_equilibrium().

Algorithm:

1. Compute vacuum field from coils
2. Seed Gaussian current distribution
3. Initial Jacobi solve
4. Picard iteration loop: a. Find O-point (axis) and X-point b. Update nonlinear source term c. Jacobi elliptic solve d. Apply vacuum boundary conditions e. Relax: Ψ = (1-α)Ψ_old + α Ψ_new f. Check convergence

pub fn psi(&self) -> &Array2<f64>

Read-only access to flux.

pub fn j_phi(&self) -> &Array2<f64>

Read-only access to current density.

pub fn state(&self) -> &PlasmaState

Read-only access to the plasma state.

pub fn grid(&self) -> &Grid2D

Read-only access to the grid.

pub fn set_solver_method(&mut self, method: SolverMethod)

Set the inner linear solver method.

pub fn solver_method(&self) -> SolverMethod

Get the current inner linear solver method.

pub fn config(&self) -> &ReactorConfig

Read-only access to the configuration.

pub fn set_external_profiles( &mut self, params_p: ProfileParams, params_ff: ProfileParams, )

Enable externally specified mTanh profiles for source update.

pub fn clear_external_profiles(&mut self)

Disable external profile mode and revert to default nonlinear source update.

pub fn set_particle_current_feedback( &mut self, particle_j_phi: Array2<f64>, coupling: f64, ) -> FusionResult<()>

Inject a static particle-current map blended into fluid J_phi each Picard step.

pub fn clear_particle_current_feedback(&mut self)

Disable particle-current feedback for subsequent solves.

pub fn set_particle_feedback_from_population( &mut self, particles: &[ChargedParticle], coupling: f64, runaway_threshold_mev: f64, ) -> FusionResult<ParticlePopulationSummary>

Build particle feedback map from macro-particle population and enable coupling.

pub fn solve_equilibrium_with_profiles( &mut self, params_p: ProfileParams, params_ff: ProfileParams, ) -> FusionResult<EquilibriumResult>

Convenience wrapper to solve with temporary external profile parameters.

pub fn sample_psi_at(&self, r: f64, z: f64) -> FusionResult<f64>

Sample solved flux at nearest grid point to (R, Z).

pub fn sample_psi_at_probes( &self, probes: &[(f64, f64)], ) -> FusionResult<Vec<f64>>

Sample solved flux at multiple probe coordinates.