Numerical Periodic Normalization at Codim 1 Bifurcations of Limit Cycles in DDEs

TL;DR

This paper develops explicit computational formulas for analyzing bifurcations of limit cycles in delay differential equations, enabling precise classification of bifurcation types through numerical methods.

Contribution

It introduces a novel combination of periodic normalization and functional analysis to derive formulas for bifurcation coefficients in DDEs, with practical software implementation.

Findings

Formulas distinguish nondegenerate, sub- and supercritical bifurcations.

Characteristic operator enables robust numerical boundary-value algorithms.

Software implementation effectively analyzes various DDE models.

Abstract

Recent work in [53, 54] by the authors on periodic center manifolds and normal forms for bifurcations of limit cycles in delay differential equations (DDEs) motivates the derivation of explicit computational formulas for the critical normal form coefficients of all codimension one bifurcations of limit cycles. In this paper, we derive such formulas via an application of the periodic normalization method in combination with the functional analytic perturbation framework for dual semigroups (sun-star calculus). The explicit formulas allow us to distinguish between nondegenerate, sub- and supercritical bifurcations. To efficiently apply these formulas, we introduce the characteristic operator as this enables us to use robust numerical boundary-value algorithms based on orthogonal collocation. Although our theoretical results are proven in a more general setting, the software implementation and examples focus on discrete DDEs. The actual implementation is described in detail and its effectiveness is demonstrated on various models.