Roman Samulyak

An efficient numerical algorithm for laser driven multi-level ionization of high-atomic-number gases is proposed and implemented in an electromagnetic particle-in-cell code SPACE. The algorithm is based on analytical solutions to the system of differential equations describing ionization evolution. Using analytical solutions resolves the multiscale issue of ionization due to different characteristic time scales of ionization processes and the main code time step. Algorithm efficiency and memory requirements are significantly improved by using a locally reduced system of differential equations. The algorithm also assigns proper orbital quantum numbers and their projections to ionization states. The algorithm is verified and validated using experimental data.

Small (∼1 mm) neon pellet fragments are fired into DIII-D H-mode plasmas, and resulting trajectory-averaged photon efficiencies (neutral neon ionization events for every photon emitted) of S/XB≈85 are estimated for Ne-I 640 nm by dividing the estimated initial pellet fragment mass by the measured number of emitted Ne-I photons. The experiments are modeled by running the Lagrangian particle (LP) fluid/magneto-hydrodynamic pellet code to estimate axial ablation plume neon density profiles and temperature profiles at each pellet position. These solutions are then fed into the PrismSPECT collisional-radiative code, which calculates resulting neon charge states and photon emission rates, giving a profile-average of S/XB≈109. The burnthrough plasma minor radius predicted by LP (ρ≈0.63) is reasonably close to the experimental observation ρ≈0.6. The modeling indicates that local S/XB is not constant along the pellet trajectory but tends to increase with increasing ablation rate. Non-equilibrium kinetics are predicted to be very important, while line trapping is predicted to be relatively unimportant (for Ne-I 640 nm S/XB).