The SIR model of epidemic spreading can be reduced to a nonlinear differential equation with an exponential nonlinearity. This differential equation can be approximated by a sequence of nonlinear differential equations with polynomial nonlinearities. The equations from the obtained sequence are treated by the Simple Equations Method (SEsM). This allows us to obtain exact solutions to some of these equations. We discuss several of these solutions. Some (but not all) of the obtained exact solutions can be used for the description of the evolution of epidemic waves. We discuss this connection. In addition, we use two of the obtained solutions to study the evolution of two of the COVID-19 epidemic waves in Bulgaria by a comparison of the solutions with the available data for the infected individuals.
Keywords: COVID-19; SIR model of epidemics; Simple Equations Method (SEsM); epidemic waves; exact solutions; nonlinear differential equations; real data for epidemic spreading of COVID-19.