Numerical research on the ion-beam-driven hydrodynamic motion of fissile targets for nuclear safety studies

As a method to evaluate high-temperature equation of state (EOS) data of fissile materials precisely and safely, we numerically examined an experimental setup based on a sub-range fissile target and a high-intensity short-pulsed heavy-ion beam. As an example, we calculated one-dimensional hydrodynamic motion of a uranium target with ρ = 0.03ρsolid (ρsolid ≡ solid density = 19.05 g/cm3) induced by a pulsed 23Na+ beam with a duration of 2 ns and a peak power of 5 GW/mm2. The projectile stopping power was calculated using a density- and temperature-dependent dielectric response function. To heat the target uniformly, we optimized the experimental condition so that the energy deposition could occur almost at the top of the Bragg peak. The energy deposition inhomogeneity could be reduced to ±5% by adjusting the incident energy and the target thickness to be 2.02 MeV/u and 180 μm, respectively. The target could be heated homogeneously up to kT =7 eV well before the arrival of the rarefaction waves at the center of the target. In principle, the EOS data can be evaluated by iteratively adjusting the data embedded in the hydro code until the measured hydrodynamic motion is reproduced by the calculation. This method is consistent with the conditions of nuclear nonproliferation, because a very small amount of fissile material is enough to perform the experiment, and no shock compression occurs in the target.