Accelerating Simulation of Two-Phase Flows with Neural PDE Surrogates