The computational simulation of the Bray-Liebhafsky (BL) oscillating chemical reaction by differential kinetic methodology is carried out in this work. According to the mechanism of Treindl and Noyes involving 10 reaction steps, the changes of the concentrations of I 2 and O 2 in solution are simulated. When the control parameters are α = 0.55, β = 0.2882 and δ < 0.6, the differential equations present periodic solution, and the oscillations can be observed in 150 min. If α , β and δ are taken as the control parameters, respectively, the bifurcation points would be observed in the processes of control parameters, changing successively with the critical values of α = 0.55,β = 0.2882, and δ = 0.6. The acidity of solution on the nonlinear phenomena is also investigated in detail.