A magnetohydrodynamic simulation of pellet ablation in the electrostatic approximation

A magnetohydrodynamic numerical model for hydrogenic pellet ablation in the electrostatic approximation has been developed based on the method of front tracking. The main features of the model are the explicit tracking of interfaces that separate the solid pellet from the ablated gas and the cold, dense and weakly ionized ablation cloud from the highly conducting fusion plasma, a surface ablation model, a kinetic model for the electron heat flux and an equation of state accounting for atomic processes in the ablation cloud. The interaction of the pellet ablation flow with the magnetic field including the J × B Lorentz force is studied here systematically for the first time. The model has also been validated through the comparison with the semi-analytic Transonic Flow model and previous purely hydrodynamic simulations. Contrary to prevailing expectations, the ablation rate is reduced only slightly when the geometry is changed from spherically symmetric to axially symmetric, in the case of purely hydrodynamic models. However, in the magnetohydrodynamic simulations the J × B force funnels the flow into an extended plasma shield, which intercepts the incident plasma heat flux and reduces the ablation rate, depending on the rise time of heat flux seen by the pellet. Shorter 'warm-up' times lead to narrower ablation channels, stronger shielding and reduced ablation rates. This new feature implies that pellets traversing strong plasma gradients, as in the edge pedestal region of the ITER plasma, could have significantly lower ablation rates if injected at higher velocity.