Scalable, conservative, and converging methods for high performance computations of turbulent liquid atomization