Ferroresonance study using Galerkin method with pseudo-arclength continuation method

The ferroresonance phenomenon shows at least two different steady states for the same applied voltage, depending on initial conditions. The Galerkin method is well suited for computing periodic steady states, but often shows difficulties in terms of initialization or convergence. To overcome these difficulties, a continuation method solution is proposed: the pseudo-arclength method. This method has been applied to draw certain characteristic curves (flux magnitude versus applied voltage) for the fundamental, third subharmonic, and third harmonic series ferroresonance case. A number of solution branches with the determination of limit points (turning and bifurcation points) are presented