Periodic orbits in glycolytic oscillators: from elliptic orbits to relaxation oscillations

We consider the Sel'kov model of glycolytic oscillator for a quantitative study of the limit cycle oscillations in the system. We identify a region of parameter space where perturbation theory holds and use both Linstedt Poincaré technique and harmonic balance to obtain the shape and frequency of the limit cycle. The agreement with the numerically obtained result is excellent. We also find a different extreme, where the limit cycle is of the relaxation oscillator variety, has a large time period and it is seen that, as a particular parameter in the model is varied, the time period increases indefinitely. We characterize this divergence numerically. A calculational method is devised to capture the divergence approximately.