Text this: A differential equations model-fitting analysis of COVID-19 epidemiological data to explain multi-wave dynamics