An accurate computation of the global solution curve for the Gelfand problem through a two point approximation

Authors:
Philip Korman
Affiliations:
Department of Mathematical Sciences, University of Cincinnati, Cincinnati, OH
Venue:
Applied Mathematics and Computation
Year:
2003

Citing 2
Cited 0

Numerical continuation and the Gelfand problem

Applied Mathematics and Computation - Special issue on differential equations and computational simulations II
Numerical methods for bifurcations of dynamical equilibria

Numerical methods for bifurcations of dynamical equilibria

Quantified Score

Hi-index	0.48

Visualization

Abstract

When one solves numerically the Gelfand boundary value problem in two dimensions (with u = u(x,y), r = √x2 + y2 and λ a positive parameter) Δu + λeu = 0 for r = 0 when r = R. with 10 mesh points, two spurious turns appear on the solution curve. When one increases the number of mesh points, the spurious turns persist, occurring at smaller values of λ. It turns out that spurious turns are avoided in the other direction, by using two to six mesh points. We prove this in case of two mesh points. This gives us the correct form of the solution curve, and the accuracy can then be improved by using a linear search on u(0), combined with shooting.