This work based on the work of hybridCentralSolvers in https://github.com/unicfdlab/hybridCentralSolvers.

The way of the blending method for convective fluxes and the kappa function used are updated to supress the unphysical oscillation at two-phase buffer regions.