A non-oscillatory method for spallation studies

This paper introduces a non-oscillatory method, the finite element flux-corrected transport (FE-FCT) method for spallation studies. This method includes the implementation of a one-dimensional FCT algorithm into a total Lagrangian finite element method. Consequently, the FE-FCT method can efficiently eliminate fluctuations behind shock wave fronts without smearing them. In multidimensional simulations, the one-dimensional FCT algorithm is used on each grid line of the structured meshes to correct the corresponding component of nodal velocities separately. The requirement of structured meshes is satisfied by using an implicit function so that arbitrary boundaries of the simulated object can be described. In this paper, the proposed FE-FCT method is applied in spallation studies. Oneand two-dimensional examples show this non-oscillatory method could be one of the candidates to accurately predict spallation and the spall thickness