“…To avoid division by zero, a tiny value γ = 10 –10 is utilized here. φ and ϕ are employed to maintain the scheme’s local monotonicity (positivity), which is defined as , φ ( α l , m ) = max ( 0 , min true( 2 α normall , m , min ( 1 3 + 2 3 α l , m , 2 ) true) ) ϕ ( β m , n ) = max ( 0 , min true( 2 β m , n , min ( 1 3 + 2 3 β m , n , 2 ) true) ) The proposed method is inapplicable at the border intervals. Thus, the backward approximation of first order is utilized for the boundary of cell fluxes, whereas the high-resolution approach is employed for all internal cell interface fluxes.…”