IJANAL International Journal of Analysis 2314-4998 2314-498X Hindawi Publishing Corporation 467831 10.1155/2014/467831 467831 Research Article Unexpected Solutions of the Nehari Problem Gritsans A. 1 Sadyrbaev F. 2 Molinet Luc 1 University of Daugavpils Parades Street 1 Daugavpils LV-5400 Latvia du.lv 2 Institute of Mathematics and Computer Science University of Latvia Raynis Boulevard 29 Riga LV-1459 Latvia lu.lv 2014 1042014 2014 28 11 2013 19 03 2014 10 4 2014 2014 Copyright © 2014 A. Gritsans and F. Sadyrbaev. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The Nehari characteristic numbers λ n ( a , b ) are the minimal values of an integral functional associated with a boundary value problem (BVP) for nonlinear ordinary differential equation. In case of multiple solutions of the BVP, the problem of identifying of minimizers arises. It was observed earlier that for nonoscillatory (positive) solutions of BVP those with asymmetric shape can provide the minimal value to a functional. At the same time, an even solution with regular shape is not a minimizer. We show by constructing the example that the same phenomenon can be observed in the Nehari problem for the fifth characteristic number λ n ( a , b ) which is associated with oscillatory solutions of BVP (namely, with those having exactly four zeros in ( a , b )) .

1. Introduction

The variational theory of eigenvalues in Sturm-Liouville problems for linear ordinary differential equations provides variational interpretation of eigenvalues which emerge as minima of some quadratic functionals being considered with certain restrictions .

As to nonlinear boundary value problems for ordinary differential equations, the Nehari theory of characteristic values provides some analogue of the linear theory. The Nehari theory deals in particular with superlinear differential equations of the form (1) x ′′ = - q ( t ) | x | 2 ε x , ε > 0 . The Nehari numbers λ n ( a , b ) , by definition, are minimal values of the functional (2) H ( x ) = a b [ x 2 ( t ) - ( 1 + ε ) - 1 q ( t ) x 2 + 2 ε ( t ) ] d t over the set Γ n of all functions x ( t ) , which are ( 1 ) continuous and piecewise continuously differentiable in [ a , b ] ; ( 2 ) there exist numbers a ν such that a = a 1 < < a n - 1 = b and x ( a ν ) = 0 in any a ν ; ( 3 ) in any [ a ν - 1 , a ν ] , x ( t ) 0 and (3) a ν - 1 a ν x 2 ( t ) d t = a ν - 1 a ν q ( t ) x 2 | x | 2 ε d t . It was proved in  (see also ) that minimizers in the above variational problem are C 2 -solutions of the boundary value problem (4) x ′′ = - q ( t ) | x | 2 ε x , x ( a ) = x ( b ) = 0 , x ( t )    has    exactly    n - 1    zeros    in    ( a , b ) . Putting (3) into (2) one gets (5) λ n ( a , b ) = min x Γ n H ( x ) = ε 1 + ε a b q ( t ) x 2 + 2 ε ( t ) d t = ε 1 + ε a b x 2 ( t ) d t , where x ( t ) is an appropriate solution of the BVP (4). The BVP (4) may have multiple solutions but not all of them are minimizers.

It appears that in order to detect λ n ( a , b ) it is sufficient to consider solutions of the boundary value problem (4).

Z. Nehari posed the question is it true that there is only one minimizer associated with λ n ( a , b ) . It was shown implicitly in  that there may be multiple minimizers associated with the first characteristic number λ 1 ( a , b ) . These minimizers are positive solutions of the problem (4) ( n = 1 ).

Later in the work by the authors  the example was constructed showing three solutions of the BVP (4). They are depicted in Figure 1.

Three solutions of the BVP from .

Two of these solutions are asymmetric and one is an even function. Surprisingly, two asymmetric solutions are the minimizers.

The same phenomenon was observed later by Kajikiya  who studied “the Emden-Fowler equation whose coefficient is even in the interval ( - 1 , 1 ) under the Dirichlet boundary condition.” It was proved that if the density of the coefficient function is thin in the interior of ( - 1,1 ) and thick on the boundary, then a least energy solution is not even. Therefore, the equation has at least three positive solutions: the first one is even, the second one is a non-even least energy solution u ( t ) , and the third one is the reflection u ( - t ) . Similar phenomena were discussed in [7, 8].

In this note, we confirm this phenomenon for the characteristic value λ 5 . Solutions for this characteristic value have exactly four zeros in an open interval ( a , b ) . We construct the example and provide all the calculations and visualizations.

2. Preliminaries 2.1. Conventions and Definitions

The problem of finding characteristic values λ n ( a , b ) is called the Nehari problem. A function x ( t ) that supplies minimal value in the Nehari problem will be called the Nehari solution. Nehari solutions associated with the equation, the interval ( a , b ) , and the number n , all are solutions of the respective boundary value problem (4).

2.2. Auxiliary Functions

In our constructions below we use the auxiliary functions called the lemniscatic sine sl t and cosine cl t . These functions can be introduced as solutions of the Cauchy problems, respectively, as follows: (6) x ( 0 ) = 0 , x ( 0 ) = 1 , x ( 0 ) = 1 , x ( 0 ) = 0 for the equation x ′′ = - 2 x 3 . These functions are much like the usual sine and cosine functions, but they are not the derivatives of each other. Instead the following holds: (7) sl t = cl t ( 1 + sl 2 t ) , cl t = - sl t ( 1 + cl 2 t ) . A complete list of formulas for the lemniscatic functions can be found in . The lemniscatic functions can be handled symbolically by the Wolfram Mathematica program using the representation in terms of the built-in Jacobi functions.

3. Construction of the Equation

Consider the interval [ - 1,1 ] . Define the piecewise linear function (8) ξ ( t ) = { h t + η , t [ - 1,0 ] , - h t + η , t [ 0,1 ] , where η = h + 1 and h > 0 is a selected number. Define q ( t ) = 2 / ξ ( t ) 6 . The function q ( t ) (depicted in Figure 2) is U-shaped function “thin” in the middle of the interval [ - 1,1 ] .

The function q ( t ) in (9) for h = 10 and η = 11 .

Consider equation (9) x ′′ = - q ( t ) x 3 , t ( - 1,1 ) , together with the boundary conditions (10) x ( - 1 ) = 0 ,          x ( 1 ) = 0 .

Consider the Cauchy problems (11) x 1 ′′ = - k ( h t + η ) 6    x 1 3 , x 1 ( - 1 ) = 0 , x 1 ( - 1 ) = β , t ( - 1,0 ) ; (12) x 2 ′′ = - k ( - h t + η ) 6    x 2 3 , x 2 ( 1 ) = 0 , x 2 ( 1 ) = - γ , t ( 0,1 ) . Let

x 1 ( t , β ) be a solution of (11) in [ - 1,0 ] ,

x 2 ( t , γ ) a solution of (12) in [ 0,1 ] .

Then, the function (13) x ( t ) = { x 1 ( t , β ) , if    - 1 t 0 , x 2 ( t , γ ) , if    0 t 1 is a C 2 -solution of the problem (9) and the problem (10) if the continuity and smoothness conditions (14) x 1 ( 0 , β ) = x 2 ( 0 , γ ) , x 1 ( 0 , β ) = x 2 ( 0 , γ ) are satisfied. The problems (11) and (12) can be solved explicitly as (15) x 1 ( t , β ) = β 1 / 2 ( h t + η ) · sl ( β 1 / 2 t + 1 h t + η ) , t [ - 1,0 ] , (16) x 2 ( t , γ ) = - γ 1 / 2 ( - h t + η ) · sl ( γ 1 / 2 t - 1 - h t + η ) , t [ 0,1 ] , where (17) β = x 1 ( - 1 ) > 0 , - γ = x 2 ( 1 ) . The derivatives can be computed as (18) x 1 ( t ; β ) = β 1 / 2 h · sl ( β 1 / 2 t + 1 h t + η ) + β - h + η h t + η · sl ( β 1 / 2 t + 1 h t + η ) . In order to get an explicit formula for a solution of the BVP (9) and (10), one has to solve a system of two equations with respect to ( β , γ ) : (19) x 1 ( 0 ; β ) = x 2 ( 0 ; γ ) , x 1 ( 0 ; β ) = x 2 ( 0 ; γ ) . This system after replacements and simplifications is (20) β 1 / 2 · sl ( β 1 / 2 η ) = γ 1 / 2 · sl ( γ 1 / 2 η ) , β 1 / 2 h · sl ( β 1 / 2 η ) + β η · sl ( β 1 / 2 η ) = - γ 1 / 2 h · sl ( γ 1 / 2 η ) -    γ η · sl ( γ 1 / 2 η ) . In new variables u : = β 1 / 2 / η and v : = γ 1 / 2 / η , the system takes the form (21) u sl u = v sl v , h u sl u + u 2 sl u = - h v sl v - v 2 sl v , h > 0 . Notice that if a solution ( u 0 , v 0 ) of the system (21) is known, then a solution x ( t ) of the BVP (9) and (10) can be constructed such that (22) x ( - 1 ) = β = u 0 2 η 2 , x ( 1 ) = - γ = - v 0 2 η 2 . If we introduce the functions (23) Φ ( u , v ) = u sl u - v sl v , Ψ ( u , v ) = h u sl u + u 2 sl u + h v sl v + v 2 sl v , the system (21) can be rewritten as (24) Φ ( u , v ) = 0 , Ψ ( u , v ) = 0 .

Zeros of the functions Φ ( u , v ) and Ψ ( u , v ) in the square Q = { ( u , v ) : 0 u , v 9 } are depicted in Figure 3. Notice that a set of zeros of Φ consists of the diagonal u = v and the “wings.”

Zeros of Φ ( u , v ) (solid line) and Ψ ( u , v ) (dashed line), h = 10 .

The intersection points of smaller “hat” with the zero set of Φ ( u , v ) reflect three solutions of the problem (9) and (10) depicted in Figure 1. The cross-point on the bisectrix relates to the even solution, and two symmetric cross-points on the “wings” relate to the remaining two solutions of asymmetric shape. The latter two solutions are “unexpected” minimizers (or, as in , “non-even least energy solutions”).

It was proved in [5, Proposition 1] that for h sufficiently large there are exactly three cross-points on a smaller “hat” (probably, for any h > 1 ). The similar proof can be conducted for the middle “hat” in Figure 3. There are exactly three cross-points (and exactly three solutions of the system (24)) that give rise to solutions of the problem (9) and the problem (10) that have exactly two zeros in the interval ( - 1,1 ) . The respective values of the Nehari functional H ( x ) (5) were calculated and the result is the same: the even solution supplies H ( x even ) = 505549 , the two solutions of asymmetric shape supply the value H ( x asym ) = 332861 .

Therefore, we confirm the phenomenon observed in [5, 6] also for oscillatory (with exactly two zeros in ( - 1,1 ) ) solutions.

3.1. Nehari Solutions with Four Internal Zeros

We study in details the case of the Nehari characteristic number λ 5 ( - 1,1 ) . Related solutions of the boundary value problem have exactly four zeros in the interval ( - 1,1 ) . Solving the system (24) on the third (counting from the origin) “hat” provides us with values (25) u 1 = 3666.804 , v 1 = 7135.523 ; u 2 = v 2 = 6275.287 ; u 3 = v 1 , v 3 = u 1 . The respective solutions are known analytically through (15) and (16) and can be computed numerically. The second way yields the three graphs depicted in Figures 4, 5 and 6.

x ( - 1 ) = 3666.80 .

x ( - 1 ) = 6275.29 .

x ( - 1 ) = 7135.52 .

In order to detect the Nehari solutions, we compute the expression (26) ε 1 + ε a b q ( t ) x 2 + 2 ε ( t ) d t in (5) for the three solutions depicted in Figure 4 to Figure 6. Recall that ε = 1 ,   a = - 1 , and   b = 1 .

Let I i be the numerical value of the above expression for solutions defined by the initial data x ( - 1 ) = 0 ,   x ( - 1 ) = ( u i η ) 2 ,   i = 1,2 , 3 .

One gets that I 1 = I 3 = 1968611.835 . The integral over the even solution is I 2 = 2582181.527 . The Nehari solutions are those of asymmetric shape.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Courant R. Hilbert D. Methods of Mathematical Physics 1953 1 New York, NY, USA Interscience MR0065391 Nehari Z. Characteristic values associated with a class of non-linear second-order differential equations Acta Mathematica 1961 105 141 175 MR0123775 10.1007/BF02559588 Nehari Z. On a class of nonlinear second-order differential equations Transactions of the American Mathematical Society 1960 95 101 123 MR0111898 10.1090/S0002-9947-1960-0111898-8 ZBL0097.29501 Sadyrbaev F. Solutions of an equation of Emden-Fowler type Differential Equations 1989 25 5 560 565 MR1003036 ZBL0792.34032 Gritsans A. Sadyrbaev F. Characteristic numbers of non-autonomous Emden-Fowler type equations Mathematical Modelling and Analysis 2006 11 3 243 252 MR2268126 ZBL1107.34024 Kajikiya R. Non-even least energy solutions of the Emden-Fowler equation Proceedings of the American Mathematical Society 2012 140 4 1353 1362 10.1090/S0002-9939-2011-11172-9 MR2869119 ZBL1248.34023 Kajikiya R. Least energy solutions of the generalized Hénon equation in reflectionally symmetric or point symmetric domains Journal of Differential Equations 2012 253 5 1621 1646 10.1016/j.jde.2012.05.004 MR2927393 ZBL1248.35062 Kajikiya R. Non-even positive solutions of the one-dimensional p -Laplace Emden-Fowler equation Applied Mathematics Letters 2012 25 11 1891 1895 10.1016/j.aml.2012.02.057 MR2957774 Gritsans A. Sadyrbaev F. Lemniscatic functions in the theory of the Emden-Fowler diferential equation Mathematics. Differential Equations 2003 3 5 27