Skip to content
Projects

Neutralizer

The equations of beamtracer.neutralizer, declared on the code that implements them and checked by the tests.

NeutralizerChannel

classsource

One neutralizer channel traced in (x, y); see the module docstring.

beam
J=IwbH,T⊥=2Eθ2J = \frac{I}{w_b H}, \qquad T_\perp = 2 E \theta^2
path fractions
fs=∑kNk∫Fs,k dsN˙L,N˙=IeHf_s = \frac{\sum_k N_k \int F_{s,k}\, ds}{\dot N L}, \qquad \dot N = \frac{I}{e H}
relaxation
f←(1−ω) f+ω ftraced,f←ftraced in the first iterationf \leftarrow (1 - \omega)\, f + \omega\, f_\mathrm{traced}, \qquad f \leftarrow f_\mathrm{traced} \text{ in the first iteration}
convergence
max⁡s∣ftraced,s−fs∣<toleranceand the plasma converged\max_s |f_{\mathrm{traced}, s} - f_s| < \text{tolerance} \quad \text{and the plasma converged}
exit
Fsexit=Is(xexit)eN˙F_s^\mathrm{exit} = \frac{I_s(x_\mathrm{exit})}{e \dot N}
Assumes
  • the channel's section in (x,y)(x, y), planar per metre of the height HH; grounded walls at y=±w/2y = \pm w/2; the beam of width wbw_b enters as negative ions, centred
  • rms divergence θ=σ⊥/v0\theta = \sigma_\perp / v_0, hence T⊥=2Eθ2T_\perp = 2 E \theta^2; the space charge reduced by the plasma's compensation
  • the gas target at the plasma's gas target density, the plasma target at its electron density; the plasma model (wellbalanced) takes the path fractions in place of fractions averaged along the axis
  • ω\omega = relaxation (0.7); at most max_iterations, the last state returned with converged=False and no warning
References
Built from