Physical constraints must be taken into account in solving partial differential equations (PDEs) in modeling physical phenomenon time evolution of chemical or biological species. In other words, numerical schemes ought to ...