Two-dimensional simulation of pellet ablation with atomic processes

A two-dimensional hydrodynamic simulation code CAP has been developed in order to investigate the dynamics of hydrogenic pellet ablation in magnetized plasmas throughout their temporal evolution. One of the properties of the code is that it treats the solid-to-gas phase change at the pellet surface without imposing artificial boundary conditions there, as done in previous ablation models. The simulation includes multispecies atomic processes, mainly molecular dissociation and thermal ionization in the ablation flow beyond the pellet, with a kinetic heat flux model. It was found that ionization causes the formation of a quasistationary shock front in the supersonic region of the ablation flow, followed by a "second" sonic surface farther out. Anisotropic heating, due to the directionality of the magnetic field, contributes to a nonuniform ablation (recoil) pressure distribution over the pellet surface. Since the shear stress can exceed the yield strength of the solid for a sufficiently high heat flux, the solid pellet can be fluidized and flattened into a "pancake" shape while the pellet is ablating and losing mass. The effect of pellet deformation can shorten the pellet lifetime almost 3x from that assuming the pellet remains rigid and stationary during ablation. (C) 2004 American Institute of Physics.