We develop a lattice Boltzmann (LB) model for immiscible two-phase flow simulations with central moments (CMs). This successfully combines a three-dimensional nonorthogonal CM-based LB scheme [De Rosis, Phys. Rev. E 95, 013310 (2017)] with our previous color-gradient LB model [Saito, Abe, and Koyama, Phys. Rev. E 96, 013317 (2017)]. Hydrodynamic melt-jet breakup simulations show that the proposed model is significantly more stable, even for flow with extremely high Reynolds numbers, up to O(106). This enables us to investigate the phenomena expected under actual reactor conditions
