Hostname: page-component-cd9895bd7-p9bg8 Total loading time: 0 Render date: 2024-12-23T12:17:21.439Z Has data issue: false hasContentIssue false

Bifurcation of positive solutions for a Neumann boundary value problem

Published online by Cambridge University Press:  17 February 2009

E. L. Montagu
Affiliation:
Mathematical Institute, 24/29 St Giles, Oxford OX1 3LB, UK.
John Norbury
Affiliation:
Mathematical Institute, 24/29 St Giles, Oxford OX1 3LB, UK.
Rights & Permissions [Opens in a new window]

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the ‘Save PDF’ action button.

Analytical, approximate and numerical methods are used to study the Neumann boundary value problem

uxx + q2u = u2(1 + sin x), for 0 < x < π,

subject to ux(0) = 0, ux(π) = 0,

for q2 ∈ (0,∞). Asymptotic approximations to (1) are found for q2 small and q2 large. In the case where q2 is large u(x) ≈ 3(xπ/2). When q2 = 0 we show that the only possible solution is u ≡ 0. However, there exist non-zero solutions for q2 > 0 as well as the trivial solution u ≡ 0. To O(q4) in the q2 small case u(x) = q2π(π + 2)−1, so that bifurcation occurs about the trivial solution branch u ≡ 0 at the first eigenvalue λ0 = 0 and in the direction of the first eigenfunction ξ0 = constant.

We obtain a bifurcation diagram for (1), which confirms that there exists a positive solution for q2 ∈ (0, 10). Symmetry-breaking bifurcations and blow-up behaviour occur on certain regions of the diagram. We show that all non-trival solutions to the problem must be positive.

The formal outer solution u = q2û appears to satisfy û = û2(1 + sin x), so that û ≡ 0 and û = (1 + sin x)−1 are possible limit solutions. However, in the non-trivial case ûx(0) = −1 and ûx(π) = 1; this means that û does not satisfy the boundary conditions required for a solution of (1). This behaviour usually implies that for q2 large a boundary layer exists near x = 0 (and one near x = π), which corrects the slope. However, we find no evidence for such a solution structure, and only find perturbations in the direction of a delta function about u ≡ 0. We show using the monotone convergence theorem for quadratic forms that the inverse of the operator on the left-hand side of (1) is strongly convergent as q2 → ∞. We show that strong convergence of the operator is sufficient to stop outer-layer behaviour occurring.

Type
Research Article
Copyright
Copyright © Australian Mathematical Society 2001

References

[1] Benjamin, T. B., “A new kind of solitary wave”, J. Fluid Mech 245 (1992) 401411.CrossRefGoogle Scholar
[2] Benjamin, T. B., “Solitary and periodic waves of a new kind”, Phil. Trans. R. Soc. Lond A354 (1996) 17751806.Google Scholar
[3] Forsythe, G. E., Malcolm, M. A. and Moler, C. B., Computer Methods for Mathematical Computations (Prentice-Hall, New York, 1977).Google Scholar
[4] Kato, T., Perturbation theory for linear operators, 2nd (corrected) ed. (Springer-Verlag, Berlin, 1966).Google Scholar
[5] Kaye, G. W. C. and Laby, T. H., Tables of Physical and Chemical Constants, 13th ed. (Longmans, London, 1966).Google Scholar
[6] Mays, L. J. H. O. and Norbury, J., “Existence of solutions for nonautonomous nonlocal elliptic equations with Neumann boundary conditions”, SIAM J. Appl. Anal. (1999), submitted for publication.Google Scholar
[7] Mays, L. J. H. O. and Norbury, J., “Existence of solution of nonlinear problems using a new method based on positive square root operators”, Proc. R. Soc. London (2000), to be resubmitted with corrections.Google Scholar
[8] Pro-Matlab, , User's Guide (The Maths Works, South Natic, Massachusetts, 1989).Google Scholar
[9] Reid, W. T., Ordinary differential equations (Wiley, New York, 1989).Google Scholar