One of the significant models in chemical reactions with oscillations is the Brusselator model. This model essentially describes a nonlinear reaction-diffusion equation. Brusselator system arises in applications of many physical and chemical models. In this study, the Brusselator model is solved numerically with the help of a time-splitting method. Consistency and stability of the method are proved with the help of auxiliary lemmas. Additionally, the positivity preservation of the method is analyzed. The accuracy of the presented method is also tested on numerical examples and all theoretical results are supported by the tables and figures.