Unspecific peroxygenases have attracted interest due to their ability to catalyze the oxygenation of various types of C-H bonds using only hydrogen peroxide as a cosubstrate. Due to the instability of these enzymes at even low hydrogen peroxide concentrations, careful fed-batch addition of the cosubstrate or ideally in situ production is required. While various approaches for hydrogen peroxide addition have been qualitatively assessed, only limited kinetic data concerning enzyme inactivation and peroxide accumulation has been reported so far. To obtain quantitative insights into the kinetics of such a process, a detailed data set for a peroxygenase-catalyzed benzylic hydroxylation coupled with electrochemical hydrogen peroxide production is presented. Based on this data set, we set out to model such an electroenzymatic process. For this, initial velocity data for the benzylic hydroxylation is collected and an extended Ping-Pong-Bi-Bi type rate equation is established, which sufficiently describes the enzyme kinetic. Moreover, we propose an empirical inactivation term based on the collected data set. Finally, we show that the full model does not only describe the process with sufficient accuracy, but can also be used predictively to control hydrogen peroxide feeding rates To limit the concentration of this critical cosubstrate in the system.