The study of deterministic models of Chemical Reaction Networks has been a central topic in Chemical Physics since the 1970s thanks to the work of M. Feinberg, F. Horn, and R. Jackson. The stochastic counterpart was formulated in the same period and investigated mostly through simulations. Stochastic simulation algorithms have been developed, the most known one being the exact simulation method usually associated to the name of Gillespie who made it popular in the field of Chemical Physics. A mathematical theory of stochastic models of Chemical Reaction Networks has been developed more recently, but efficient simulation methods are still needed to address most models that remains intractable with analytical tools.