Abstract
In this note, we study the oscillating global continua of the differential inclusion of the form
where F is a "setvalued representation" of a function with jump discontinuities along the line segment [0, 1] × {0}, and λ ∈ [0, ∞) is a parameter. The proof of our main result relies on an approximation procedure.
Mathematics Subject Classification 2000: 34B16; 34B18.
Keywords:
climate model; differential inclusion; eigenvalue; positive solutions1 Introduction
In recent years, nonsmooth analysis has come to play an important role in functional analysis [1], dynamical systems [2], control theory [3], optimization [4], mechanical systems [5], differential equation [6,7] etc. Since many mathematical and physical problems may be reduced to ODES or PDES with discontinuous nonlinearities, the existence of multiple solutions for differential inclusion problems has been widely investigated [819].
In this article, we are concerned with the following differential inclusion problem which raises from a BudykoNorth type energy balance climate models:
(see [2025] and the references therein). In particular, the setvalued right hand side arise from a jump discontinuity of the albedo at the iceedge in these models. By filling such a gap, one arrives at the setvalued problem (1.1). As in [25], we are here interested in a considerably simplify version as compared to the situation from climate modeling, e.g. a onedimensional regular SturmLiouville differential operator substitutes for a twodimensional LaplaceBeltrami operator or a singular Legendretype operator, and the jump discontinuity is transformed to u = 0 in a way, which resembles only locally the climatological problem.
We are concerned with the setvalued problem (1.1) under the following assumptions
(H1) q∈C([0, 1],(0,+∞));
(H2) f^{+}∈ C ([0, 1] × [0,+∞), (0,+∞)), .
Let the setvalued function F in (1.1) is given by
Notice that if f^{+}(x, 0) ≡ 0, x ∈ [0, 1], then the differential inclusion problem (1.1) reduces to the BVP of differential equation
In the last 20 years, the positive solutions of (1.2) have been studied by several authors, see Jiang and Liu [26], Chu et al. [27] and Sun et al. [28].
The purpose of this article is to investigate the oscillating global continua of positive solutions of the differential inclusion problem (1.1). The proof of our main result relies on an approximation procedure. The rest of the article is organized as follows. In Section 2, we state some notations and prove some preliminary results. In Section 3, we state and prove our main result. In Section 4, an example is given to illustrate the application of our main result.
2 Notations and preliminaries
Recall Kuratowski's notion of lower and upper limits of sequence of sets.
Definition 2.1. [29]Let X be a metric space and {Z_{l}}_{l∈ℕ}be a sequence of subsets of X. The set
is called the upper limit of the sequence {Z_{l}}, whereas
is called the lower limit of the sequence {Z_{l}}.
Definition 2.2. [29]A component of a set M is meant a maximal connected subset of M.
Lemma 2.1. [29]Suppose that Y is a compact metric space, A and B are nonintersecting closed subsets of Y, and no component of Y intersects both A and B. Then there exist two disjoint compact subsets Y_{A}and Y_{B}, such that Y = Y_{A}∪Y_{B}, A ⊂ Y_{A}, B ⊂ Y_{B}.
Using the above Whyburn Lemma, Ma and An [30] proved the following
Lemma 2.2. [30, Lemma 2.1] Let Z be a Banach space and let {A_{n}} be a family of closed connected subsets of Z. Assume that
(i) there exist z_{n} ∈ A_{n}, n = 1, 2, ..., and z* ∈ Z, such that z_{n} → z*;
(ii) r_{n} = sup {∥x∥  x ∈ A_{n}} = ∞;
(iii) for every R > 0, is a relatively compact set of Z, where B_{R}= {x ∈ Z  ∥x∥ ≤ R}. Then there exists an unbounded component in and .
Remark 2.1. The limiting processes for sets go back at least to the work of Kuratowski [31]. Lemma 2.2 is a slight generalization of the following wellknow result due to Whyburn [29]:
Proposition 2.1. (Whyburn [29, p. 12]) Let Z be a Banach space and {A_{n}} be a family of closed connected subsets of Z. Let and ∪_{l∈ℕ}A_{l} is relatively compact. Then is nonempty, compact and connected.
Next, we introduce the result of global solution behavior of the bifurcation branches of the equation
to wit the following lemma.
Lemma 2.3. [32] (Dancer (1974)) Assume that
(C1) The operators L, N: X → X are compact on the Banach space X over R. Furthermore, L is linear and ∥Nx∥/∥x∥ → 0 as ∥x∥ → 0;
(C2) The real number μ_{0}is a characteristic number of L of odd algebraic multiplicity;
(C1+) The real Banach space X has an order cone K with X = KK, i.e., every x ∈ X can be represented as x = x_{1}  x_{2}, where x_{1}, x_{2} ∈ K. Furthermore, L + N is positive, i.e., L + N maps K into K;
(C2+) The spectral radius r(L) of L is positive. We set μ_{0} = (r(L))^{1}.
Then (μ_{0}, 0) is a bifurcation point of equation (2.1) and
contains an unbounded solution component which passes through (μ_{0}, 0).
If additionally
(C3+) The linear operator L is strongly positive, then and μ ≠ μ_{0}always implies x > 0 and μ > 0.
Remark 2.2. This result is often called the nonlinear KreinRutman theorem. It will play an important role in the proof of our main result.
Let φ and ψ be the unique solution of the problems
and
respectively. Then it is easy to check φ(·) is nondecreasing on (0,1), ψ(·) is nonincreasing on (0,1), and the Green's function G(x, s) of
is explicitly given by
Moreover, we have that
3 The main result
Let Σ be the closure of the set of positive solutions of (1.1) in [0, ∞) × C^{1}[0, 1], and ℕ* := {1, 2,..., N}. The main result of this article is the following theorem.
Theorem 3.1. Assume that (H1)(H2) hold. If
(H3) there is an increasing sequence of positive numbers and a small enough constant δ such that ξ_{1} < σ(ξ_{2} δ) and
where
then there exits an unbounded component in Σ with . Moreover,
(i) with ∥u∥_{∞} = ξ_{2j1}for some j ∈ ℕ* implies that λ ≥ 2;
(ii) with ∥u∥_{∞} = ξ_{2j}for some j ∈ ℕ* implies that .
Actually, such continua can be obtained as upper limits in the sense of Kuratowski of sequence of solution continua from associated continuous problems. To this end one sets
fixes l_{0} ∈ ℕ such that , and selects an approximation sequence {f_{l}} ⊂ C ([0, 1] × ℝ, ℝ) (l > l_{0}) of F satisfying:
(A1) f_{l} (x, y) = ly for x ∈ [0, 1] and ;
(A3) f_{l}(x,y) = f^{+}(x, y) for x ∈ [0, 1] and ;
(A4) {f_{l} (x, y)}_{l∈ℕ} is nondecreasing in l for (x, y) ∈ [0, 1] × (0,∞).
Next, we show that the continuous problem
has an unbounded closed subsets of the positive solutions set of (3.2_{l}) with
(a) is the bifurcation point contained in ;
(b) If and ϑ ≢ 0, then ϑ is positive on (0,1).
It is easy to see that (3.2_{l}) equivalent to
Let
Then according to (3.3), (3.2_{l}) can be written as the following operator equation
Clearly, the operators L, N : C[0, 1] → C[0, 1] are compact on the Banach space
C[0, 1]. Furthermore, L is linear and thanks to (2.3)(A1) that
which implies that the condition (C1) of Lemma 2.3 is satisfied.
Denote the principal eigenvalue of
by λ_{1}, then we know that λ_{1}> 0 (see [33]). Since (3.4) is equivalent to operator equation
we have that . Therefore, the conditions (C2)(C2+) of Lemma 2.3 are satisfied.
Let the cone K in C[0, 1] is given by
It is easy to see thanks to (A1)(A4) and (2.3) that the (C1+)(C3+) conditions of Lemma 2.3 are satisfied.
According to Lemma 2.3, we obtain that is a bifurcation point of the positive solutions set of (3.2_{l}) for every l ∈ {l_{0} + 1, l_{0} + 2, ...} =: ℕ_{0}, and for each l ∈ ℕ_{0} there exits an unbounded closed subsets of the positive solutions set of (3.2_{l}) with (a) and (b).
Combining the above with the fact
and utilizing Lemma 2.2, it concludes that there exits an unbounded component with
and
Denote the cone P in C[0, 1] by
Define an operator T_{λ} : P → C[0, 1] by
It is easy to get the following lemma.
Lemma 3.1. Assume that (H1), (H2) and (A1)(A4) hold. Then T_{λ}: P → P is completely continuous.
Lemma 3.2. Assume that (H1), (H2) and (A1)(A4) hold. If 0 ≤ u(x) ≤ r, r > 0, for x ∈ [0, 1], then
Proof. Since f_{l}(x, u(x)) ≤ M_{r} for x ∈ [0, 1], it follows from (2.3) that
Lemma 3.3. Assume that (H1), (H2) and (A1)(A4) hold. If σ(r  δ) ≤ u(x) ≤ r + δ, r > δ, for x ∈ [0, 1], then
Proof. Since f_{l}(x, u(x)) ≥ m_{r} for x ∈ [0, 1], it follows that
Lemma 3.4. Assume that (H1), (H2), (H3) and (A1)(A4) hold. then
(i) with ∥u∥_{∞} ∈ (ξ_{2j1}  δ,ξ_{2j1} + δ) for some j ∈ ℕ* implies that λ > 2;
(ii) , with ∥u∥_{∞} ∈ ( ξ_{2j} δ,ξ_{2j} + δ) for some j ∈ ℕ* implies that .
Proof. (i) Assume that with ∥u∥_{∞} ∈ (ξ_{2j1}  δ, ξ_{2j1} + δ) for some j ∈ ℕ*, then u = T_{λ}u and
By Lemma 3.2 and (H3), it follows that
Thus λ > 2.
(ii) Assume that with ∥u∥_{∞} ∈ (ξ_{2j}  δ, ξ_{2j} + δ) for some j ∈ ℕ*, then u = T_{λ}u and
By Lemma 3.3 and the assumption (H3), it follows that
Lemma 3.5. If , then (λ, u) is a solution of (1.1) and u ∈ W^{2,∞}(0, 1).
Proof. Let . By the definition of there exists a sequence {l_{k}} ∈ ℕ_{0} strictly increasing, and with for k ∈ ℕ and
Since is uniformly bounded, i.e.
we can assume after passing to a subsequence, if necessary, that it converges weekly in L^{2}(0, 1) to some ϕ. We claim that ϕ(x) ∈ F(x, u(x)) a.e. on (0, 1).
Let x_{0} ∈ (0, 1) with u(x_{0}) > 0. Then there exist ρ > 0 and τ ∈ (0, min{x_{0}, 1x_{0}}) with u(x) > ρ for all x ∈ (x_{0}  τ, x_{0} + τ), hence there is a k_{0} ∈ ℕ with for all k > k_{0} and x ∈ (x_{0}  τ, x_{0} + τ). Choose k_{1} > k_{0} with . Then for all k ≥ k_{1} and x ∈ (x_{0}  τ, x_{0} + τ), which yields ϕ(x) = f^{+}(x, u(x)) for x ∈ (x_{0}  τ, x_{0} + τ) a.e.
Next, if u ≡ 0, let K: = {x ∈ (0, 1) : ϕ(x) > f^{+}(x, 0)}. We claim that meas(K) = 0. Suppose that meas(K) > 0. Then ε := ∫_{K} [ϕ(x)  f^{+}(x, 0)] dx > 0, and one finds η ∈ (0, ∞) with for x ∈ [0, 1] and y ∈ [0, η]. Choosing k_{2} ∈ ℕ with for k ≥ k_{2}. One obtains for k ≥ k_{2}:
which contradicts . Thus, meas(K) = 0.
Now, let A be the closed linear operator in L^{2}(0, 1) defined by
and Aφ := φ" + qφ. Clearly,
hence and the fact that A is weakly closed yields
i.e.
Finally, we show that u ∈ W^{2,∞}(0,1). In fact, from (3.9) we have
According to (H1) and the boundedness of u we have
We claim that ϕ ∈ L^{∞}(0,1). Suppose on the contrary that there exists a set E ⊂ [0, 1], meas(E) > 0 such that ϕ is unbounded on E. Without loss of generality, we assume that
where M is given by (3.7) and w ∈ L^{2}(0,1). On the one hand, for k larger enough from (3.7), (3.8) and (H2) we have
On the other hand, from (3.12) we have
which contradicts (3.13). Thus,
Therefore, from (3.10), (3.11) and (3.14) we obtain u ∈ W^{2,∞}(0,1).
Now we are in the position to prove Theorem 3.1.
Proof of Theorem 3.1.
Assume that . We divide the proof into two cases.
Case l. If ∥u∥_{∞} = ξ_{2j1} for some j ∈ ℕ*, then λ ≥ 2.
Since , there exists a sequence , such that
Hence, for δ > 0 there exists i_{0} ∈ ℕ, such that
i.e.
By using Lemma 3.4, we obtain that
Hence, we get
Case 2. If ∥u∥_{∞} = ξ_{2j} for some j ∈ ℕ*, then .
Since , there exists a sequence , such that
Hence, for δ > 0 there exists i_{0} ∈ ℕ, such that
i.e.
By using lemma 3.4, we obtain that
Hence, we get
Corollary 3.1. Assume that (H1)(H3) hold. Then
(i) for each , (1.1) has at least one positive solution: ;
(ii) for each , (1.1) has N positive solutions:
Proof. According to Theorem 3.1, the boundary value problem (1.1) has an unbounded component in Σ with . Moreover,
with ∥u∥_{∞} = ξ_{2j1} for some j ∈ ℕ* implies that λ ≥ 2;
with ∥u∥_{∞} = ξ_{2j} for some j ∈ ℕ* implies that .
From the facts , and with ∥u∥_{∞}= ξ_{1} implies that λ ≥ 2 and the connectivity of , we obtain
which implies for each , (1.1) has at least one positive solution: .
Let
where ξ_{0} = 0, ξ_{k} (k = 1, 2,..., N) is given by (H3). Then according to and the connectivity of , we obtain
which implies for each , (1.1) has N positive solutions:
4 Example
In this section, an example is given to illustrate the application of our main result (Theorem 3.1). Consider second order Neumann differential inclusion problem
where the setvalued function F in (4.1) is given by
Obviously, (H1), (H2) conditions of Theorem 3.1 are satisfied. Moreover, Green's function of the associated linear problem
can be explicitly expressed by
By calculation we can get and .
Let ξ_{1} = 3, ξ_{2} = 11, δ = 1, then we can check that ξ_{1} = 3 < 5 < σ(ξ_{2} δ), and
So that (H3) condition of Theorem 3.1 is satisfied. Therefore, according to Theorem 3.1 the differential inclusion problem (4.1) has an unbounded component in Σ with . Moreover,
Competing interests
The authors declare that they have no competing interests.
Acknowledgements
The authors express their gratitude to Professors Ma Tian and Ma Ruyun for their guidance and encouragement, also to an anonymous referee for a number of valuable comments and suggestions.
References

Aubin, JP, Cellina, A: Differential Inclusion, vol. 264. Grundlehren Math Wiss, SpringerVerlag, Berlin (1984)

Kunze, M: Nonsmooth dynamical systems. Lecture Notes in Mathematics, vol. 1744. SpringerVerlag, Berlin (2000)

Clarke, FH, Ldeyaev, YS, Stern, RJ, Wolenski, PR: Nonsmooth Analysis and Control Theory. SpringerVerlag, New York (1998)

Clarke, FH: Optimization and Nonsmooth Analysis. SIAM, Philadelphia (1990)

Leine, RI, Nijjmeijer, H: Dynamics and bifurcation of nonsmooth mechanical systems. Lecture Notes in Applied and Computational Mechanics, vol. 18. SpringerVerlag, Berlin (2004)

Gasiésdi, L, Papageorgion, NS: Nonsmooth Critical Point Theory and Nonlinear Boundary Value Problems. Chapman and Hall/CRC, Boca Raton (2005)

Deimling, K: Multivalued Differential Equation. SpringerVerlag, Berlin (1985)

Kowalczyk, P, Piiroinen, PT: Twoparameter sliding bifurcation of periodic solutions in a dryfriction oscillator. Physica D. 237, 1053–1073 (2008). Publisher Full Text

Deimling, K: Resonance and clulomb friction. Diff Integr Equ. 7(3), 759–765 (1994)

Ma, R: Existence of periodic solutions of a generalized friction oscillator. Nonlinear Anal Real World Appl. 11, 3316–3322 (2010). Publisher Full Text

Chang, KC: Variational methods for nondifferentable functionals and their applications to partial differential equations. J Math Anal Appl. 80, 102–112 (1981). Publisher Full Text

Kourogenis, NC, Papageorgiou, NS: Nonsmooth critical point theory and nonlinear elliptic equation at resonance. Kodai Math J. 23, 128–135 (2000)

Zykov, PS: On twopoint boundary value problems for secondorder differential inclusions on manifolds. Appl Anal. 88(6), 895–902 (2009). Publisher Full Text

Hannelore, L, Csaba, V: Multiple solutions for a differential inclusion problem with nonhomogeneous boundary conditions. Numer Funct Anal Optim. 30(56), 566–581 (2009). Publisher Full Text

Ntouyas, SK, O' Regan, D: Existence results for semilinear neutral differential inclusions with nonlocal conditions. Diff Equ Appl. 1(1), 41–65 (2009)

Papageorgion, NS, Staicu, V: The method of upperlower solutions for nonlinear second order differential inclusions. Nonlinear Anal. 67(3), 708–726 (2007). Publisher Full Text

Kyritsi, S, Matzakos, N, Papageorgion, NS: Periodic problems for strongly nonlinear secondorder differential inclusions. J Diff Equ. 183, 279–302 (1982)

Dhage, BC, Ntouyas, SK, Cho, YJ: On the second order discontinuous differential inclusions. J Appl Funt Anal. 1(4), 469–476 (2006)

Benchchra, M, Graef, JR, Ouahab, A: Oscillatory and nonoscillatory solutions of multivalued differential inclusions. Comput Math Appl. 49(910), 1347–1354 (2005). Publisher Full Text

Budyko, MI: The effect of solar radiation variations on the climate of the earth. Tellus. 21, 611–619 (1969). Publisher Full Text

Diaz, JI: Mathematical Analysis of Some Diffusive Energy Balance Models. Math Climate Environ, Mason, Paris (1993)

Diaz, JI: The Mathematics of Models for Climatology and Environment, NATO ASI Series I: Global Environmental Changes, vol. 48. SpringerVerlag, New York (1997)

Diaz, JI, Hernandez, J, Tello, L: On the multiplicity of equilibrium solutions to a nonlinear diffusion equation on a manifold arising in climatology. J Math Anal Appl. 216, 593–613 (1997). Publisher Full Text

HendersonSellers, A, McGuffie, KA: A Climate Modeling Primer. Wiley, Chichester (1987)

Hetzer, G: A bifurcation result for SturmLiouville problem with a setvalued term. Mississippi State University (1997) In Proceedings of the Third Mississippi State Conference on Difference Equations and Computational Simulations: 1617 May 1997

Jiang, D, Liu, H: Existence of positive solutions to second order Neumann boundary value problem. J Math Res Expo. 20, 360–364 (2000)

Chu, JF, Sun, YG, Chen, H: Positive solutions of Neumann problems with singularities. J Math Anal Appl. 337, 1267–1272 (2008). Publisher Full Text

Sun, Y, Cho, YJ, O' Regan, D: Positive solution for singular second order Neumann boundary value problems via a cone fixed point theorem. Appl Math Comput. 210, 80–86 (2009). Publisher Full Text

Whyburn, GT: Topological Analysis. Princeton University Press, Princeton, NJ (1964)

Ma, R, An, Y: Global structure of positive solutions for nonlocal boundary value problems involving integral conditions. Nonlinear Anal TMA. 71(10), 4364–4376 (2009). Publisher Full Text

Zeidler, E: Nonlinear Functional Analysis and its Applications I (Fixedpoint theorems). SpringerVerlag, New York (1986)

Mavinga, N, Nkashama, MN: SteklovNeumann eigenproblems and nonlinear elliptic equations with nonlinear boundary conditions. J Diff Equ. 248, 1212–1229 (2010). Publisher Full Text