Enhanced schemes for resolution of the continuity equation in projection-based SPH

Abstract

To the best knowledge of the authors, the numerical resolution of the continuity equation, i.e., the simultaneous achievement of accurate and physically consistent velocity divergence and density fields, has not been investigated in the context of projection-based particle methods. In this study, based on the Velocity-divergence Error Mitigating (VEM)/Volume Conservation Shifting (VCS) schemes developed for the WCSPH method, which improve the numerical resolution of the continuity equation, VEM/VCS schemes are proposed for projection-based particle methods. Furthermore, two other novel schemes are introduced. First, a Dynamic VEM (DVEM) scheme is proposed, where VEM is dynamically activated only when the velocity divergence instantaneously diverge from the reference solution at individual particles. The second scheme is referred to as Enhanced VCS by reference density (EVCS), where precise imposition of dynamic free surface boundary condition is considered. Several benchmark tests with different features have been conducted, and it is shown that the proposed schemes enhance the numerical resolution of the continuity equation for incompressible flows, i.e., the divergence free velocity and the constant density conditions. Application of DVEM and EVCS has provided enhanced resolution of the continuity equation along with clear reductions in the numerical energy dissipation associated with the approximated additional numerical terms.