In this paper, we apply the particle method to solve the numerical solution of a family of non-li-near Evolutionary Partial Differential Equations. It is called b-equation because of its bi-Hamiltonian structure. We introduce the particle method as an approximation of these equations in Lagrangian representation for simulating collisions between wave fronts. Several numerical examples will be set to illustrate the feasibility of the particle method.