GNGA for general regions: Semilinear elliptic PDE and crossing eigenvalues

We consider the semilinear elliptic PDE Δu + f(λ, u) = 0 with the zero-Dirichlet boundary condition on a family of regions, namely stadions. Linear problems on such regions have been widely studied in the past. We seek to observe the corresponding phenomena in our nonlinear setting. Using the Gradient Newton Galerkin Algorithm (GNGA) of Neuberger and Swift, we document bifurcation, nodal structure, and symmetry of solutions. This paper provides the first published instance where the GNGA is applied to general regions. Our investigation involves both the dimension of the stadions and the value λ as parameters. We find that the so-called crossings and avoided crossings of eigenvalues as the dimension of the stadions vary influences the symmetry and variational structure of nonlinear solutions in a natural way.