In this manuscript, the fractional-in-time NoyesField Casein Protein model for Belousov-Zhabotinsky reaction transport is considered with a novel numerical technique, which was used to approximate the Atangana-Baleanu (ABC) operator which models the subdiffusion partial derivative in time.The effect of the ABC operator is observed and Oven Glass Spacer captured more interesting physical behavior of some real-life phenomena.Applicability and suitability of the adopted method were carried out on some cases of nonlinear Belousov-Zhabotinsky sub-reaction-diffusion models, their dynamic behaviors with respect to fractional-order parameters were displayed in figures.