Skip to main content

Periodic solutions to the Liénard type equations with phase attractive singularities

Abstract

Sufficient conditions are established guaranteeing the existence of a positive ω-periodic solution to the equation

u +f(u) u +g(u)=h(t,u),

where f,g:(0,+)R are continuous functions with possible singularities at zero and h:[0,ω]×RR is a Carathéodory function. The results obtained are rewritten for the equation of the type

u + c u u μ + g 1 u ν g 2 u γ = h 0 (t) u δ ,

where g 1 , g 2 , δ are non-negative constants, c, μ, ν, γ are real numbers, and h 0 L([0,ω];R). The last equation also covers the so-called Rayleigh-Plesset equation, frequently used in fluid mechanics to model the bubble dynamics in liquid. In the paper, the case when ν>γ, i.e., the case which covers the attractive singularity of the function g, is studied. The results obtained assure that there exists a positive ω-periodic solution to the above-mentioned equation if the power μ or ν is sufficiently large.

MSC:34C25, 34B16, 34B18, 76N15.

1 Introduction

The topic of singular boundary value problems has been of substantial and rapidly growing interest for many scientists and engineers. The importance of such investigation is emphasized by the fact that numerical simulations of solutions to such problems usually break down near singular points.

On the other hand, problems of this type arise frequently in applied science. Namely, in fluid mechanics, since 1917 the physicists have used the Rayleigh equation,

ρ [ R R ¨ + 3 2 R ˙ 2 ] =p(R) p ,

to model the bubble dynamics in liquid, where R(t) is the ratio of the bubble at the time t, ρ is the liquid density, p is the pressure in the liquid at a large distance from the bubble, and p(R) is the pressure in the liquid at the bubble boundary. In 1949, Plesset proposed to use a more exact equation involving the surface-tension constant S and the coefficient of the liquid viscosity μ , which was finally improved by adding a term with polytropic coefficient k in 1977, nowadays known as a Rayleigh-Plesset equation (see [1])

ρ [ R R ¨ + 3 2 R ˙ 2 ] = [ P v p ( t ) ] + P g 0 ( R 0 R ) 3 k 2 S R 4 μ R ˙ R .

The transformation R= u 2 5 in the previous equation leads to the equation

u ¨ = 5 [ P v p ( t ) ] 2 ρ u 1 5 + ( 5 P g 0 R 0 3 k 2 ρ ) 1 u 6 k 1 5 5 S u 1 5 4 μ u ˙ u 4 5 .

Consequently, the class of equations

u (t)+ c u ( t ) u μ ( t ) + g 1 u ν ( t ) g 2 u γ ( t ) = h 0 (t) u δ (t)for a.e. t[0,ω],
(1.1)

with non-negative constants g 1 , g 2 , δ, real numbers c, μ, ν, γ, and h 0 L([0,ω];R), plays an important role in fluid mechanics. Therefore, the equation

u (t)+f ( u ( t ) ) u (t)+g ( u ( t ) ) =h ( t , u ( t ) ) for a.e. t[0,ω],
(1.2)

subjected to the periodic conditions

u(0)=u(ω), u (0)= u (ω),
(1.3)

is investigated in the presented paper. Here, f,g:(0,+)R are continuous, having possible singularities at zero, and h:[0,ω]×RR is a Carathéodory function, i.e., h(,x):[0,ω]R is measurable for all xR, h(t,):RR is continuous for a.e. t[0,ω], and for every r>0, there exists a non-negative function q r L([0,ω];R) such that |h(t,x)| q r (t) for a.e. t[0,ω] and all |x|r. By a solution to (1.2), (1.3) we understand a function u:[0,ω]R which is positive, absolutely continuous together with its first derivative, satisfies (1.2) almost everywhere on [0,ω], and verifies (1.3). In spite of the fact that the problem (1.2), (1.3) was investigated by many mathematicians (see, e.g., [227]), most of the mentioned works deal with the repulsive case and/or when f has no singularity. However, the physical model, covered by equation (1.1), justifies considering the types of equations with a singular friction-like term.

A particular case of (1.1) is the equation

u + 1 u ν = h 0 (t)
(1.4)

studied by Lazer and Solimini. Their results were published in 1987 (see [11]) and they proved, among others, that (1.4), (1.3) has at least one solution if and only if h ¯ 0 >0, provided h 0 is bounded. Recently, we have proved (see [28]) that this result cannot be extended to the case when h 0 is a general integrable (and so unbounded) function unless some additional conditions are introduced. In particular, (1.4), (1.3) is solvable for any h 0 L([0,ω];R) with h ¯ 0 >0 if ν1; and, moreover, for any ν(0,1), there exists a function h 0 L([0,ω];R) with h ¯ 0 >0 such that (1.4), (1.3) has no solution. At this point, we would like to emphasize the important fact that the condition ν1 can be weakened if (1.4) is generalized to equation (1.1), see Remark 2.2 below.

The structure of the paper is as follows. After the introduction and basic notation, we recall the definition of lower and upper functions to the problem (1.2), (1.3), and we formulate the classical theorem on the existence of a solution to (1.2), (1.3) in the case when there exists a couple of well-ordered lower and upper functions. In Section 2, we establish our main results and their consequences. Sections 3 and 4 are devoted to auxiliary propositions and proofs of the main results, respectively.

For convenience, we finish the introduction with a list of notations which are used throughout the paper:

is the set of all natural numbers, is the set of all real numbers, R + =(0,+), [ x ] + =max{x,0}, [ x ] =max{x,0}.

C([0,ω];R) is the Banach space of continuous functions u:[0,ω]R with the norm

u =max { | u ( t ) | : t [ 0 , ω ] } .

C( R + ;R), resp. C( R + ; R + ), is the set of continuous functions u: R + R, resp. u: R + R + .

C 1 ( R + ; R + ) is the set of functions u: R + R + which are continuous together with their first derivative.

A C 1 ([0,ω];R) is a set of all functions u:[0,ω]R such that u and u are absolutely continuous.

L([0,ω];R) is the Banach space of the Lebesgue integrable functions p:[0,ω]R endowed with the norm

p 1 = 0 ω | p ( s ) | ds.

For a given pL([0,ω];R), its mean value is defined by

p ¯ = 1 ω 0 ω p(s)ds.

Given φ,ψL([0,ω];R), then

Φ + = 0 ω [ φ ( s ) ] + d s , Φ = 0 ω [ φ ( s ) ] d s , Ψ + = 0 ω [ ψ ( s ) ] + d s , Ψ = 0 ω [ ψ ( s ) ] d s .

The following definitions of lower and upper functions are suitable for us. For more general definitions, one can see, e.g., [[12], Definition 8.2].

Definition 1.1 A function αA C 1 ([0,ω];R) is called a lower function to the problem (1.2), (1.3) if α(t)>0 for every t[0,ω] and

α ( t ) + f ( α ( t ) ) α ( t ) + g ( α ( t ) ) h ( t , α ( t ) ) for a.e.  t [ 0 , ω ] , α ( 0 ) = α ( ω ) , α ( 0 ) α ( ω ) .

Definition 1.2 A function βA C 1 ([0,ω];R) is called an upper function to the problem (1.2), (1.3) if β(t)>0 for every t[0,ω] and

β ( t ) + f ( β ( t ) ) β ( t ) + g ( β ( t ) ) h ( t , β ( t ) ) for a.e.  t [ 0 , ω ] , β ( 0 ) = β ( ω ) , β ( 0 ) β ( ω ) .

The following theorem is well known in the theory of differential equations (see, e.g., [[12], Theorem 8.12]).

Theorem 1.1 Let α and β be lower and upper functions to the problem (1.2), (1.3) such that

α(t)β(t) for t[0,ω].
(1.5)

Then there exists a solution u to the problem (1.2), (1.3) such that

α(t)u(t)β(t) for t[0,ω].

2 Main results

Theorem 2.1 Let ρ 0 C 1 ( R + ; R + ) and ρ 1 C( R + ; R + ) be non-decreasing functions, h 0 , h 1 L([0,ω];R), and x 0 >0 be such that

h 1 (t) ρ 1 (x)h(t,x) h 0 (t) ρ 0 (x) for a.e. t[0,ω],x x 0 ,
(2.1)

and let there exist c 0 , c 1 R such that

g ( x ) ρ 0 ( x ) c 0 < h ¯ 0 for x x 0 ,
(2.2)
g ( x ) ρ 1 ( x ) c 1 h ¯ 1 for x x 0 .
(2.3)

Let, moreover, there exist λ[0,1] such that

0 1 d s ρ 0 λ ( s ) <+,
(2.4)
lim x 0 + g ( x ) ρ 0 λ ( x ) =+,
(2.5)

and let either

0 1 ( [ f ( s ) ] + ρ 0 λ ( s ) + [ g ( s ) ] + ρ 0 2 λ ( s ) ) ds=+, 0 1 [ f ( s ) ] ρ 0 λ ( s ) ds<+
(2.6)

or

0 1 ( [ f ( s ) ] ρ 0 λ ( s ) + [ g ( s ) ] + ρ 0 2 λ ( s ) ) ds=+, 0 1 [ f ( s ) ] + ρ 0 λ ( s ) ds<+.
(2.7)

Furthermore, let us suppose that ρ 0 fulfills at least one of the following conditions:

  1. (a)

    there exists a sequence { y n } n = 1 + of positive numbers such that

    lim n + y n =+, lim n + ρ 0 1 λ ( y n ) σ ( y n ) =0,
    (2.8)

and there exist ε 0 >0, ε 1 (0, ε 0 ], and n 0 N such that

ρ 0 1 λ ( ( 1 + ε 0 ) y n ) ρ 0 1 λ ( y n ) Φ Φ + ε 0 for n n 0 ,
(2.9)
(1+ ε 1 )σ( y n )σ ( ( 1 + ε 0 ) y n ) for n n 0 ,
(2.10)

where φ(t)= h 0 (t) c 0 for almost every t[0,ω] and

σ(x)= 0 x d s ρ 0 λ ( s ) ;
(2.11)
  1. (b)

    the function ρ 0 1 λ ( x ) σ ( x ) is non-increasing and

    ω 4 Φ + Φ ρ 0 1 λ ( x 0 ) σ ( x 0 ) < Φ + Φ ,
    (2.12)

where φ(t)= h 0 (t) c 0 for almost every t[0,ω] and σ is given by (2.11).

Besides, let us suppose that ρ 1 fulfills at least one of the following conditions:

  1. (c)

    there exists a sequence { z n } n = 1 + of positive numbers such that

    lim n + z n =+, lim n + ρ 1 ( z n ) z n =0,

and there exist ε 2 >0 and n 1 N such that

ρ 1 ( z n ( 1 + ε 2 ) ) ρ 1 ( z n ) Ψ Ψ + for n n 1 ,

where ψ(t)= h 1 (t) c 1 for almost every t[0,ω];

  1. (d)

    the function ρ 1 ( x ) x is non-increasing and

    ω 4 Ψ + Ψ ρ 1 ( x 0 ) x 0 Ψ + Ψ ,

where ψ(t)= h 1 (t) c 1 for almost every t[0,ω].

Then there exists at least one solution to the problem (1.2), (1.3).

Remark 2.1 Note that there exists a suitable ε 1 such that (2.10) holds, e.g., if

lim sup x + ρ 0 λ ( ( 1 + ε 0 ) x ) ρ 0 λ ( x ) <1+ ε 0 .

For equation (1.1), from Theorem 2.1 we get the following assertion.

Corollary 2.1 Let g 1 >0, g 2 0, 0δ<1, ν>γ, ν+δ>0 and

either(μ+δ)sgn|c|1orν+2δ1.

If

h ¯ 0 > lim x + g 2 x γ + δ ,

then (1.1), (1.3) has at least one solution.

Remark 2.2 In [28], it is proved, among others, that the equation

u + 1 u ν = h 0 (t),
(2.13)

with h 0 L([0,ω];R) and h ¯ 0 >0, has a positive ω-periodic solution if ν1. Moreover, there is also an example introduced showing that for any ν(0,1), there exists h 0 L([0,ω];R) with h ¯ 0 >0 such that (2.13), (1.3) has no positive solution.

Corollary 2.1 says that if a friction-like term or sub-linear term are added to (2.13), the condition ν1 can be weakened. For example,

u + u u μ + 1 u ν = h 0 (t)

has a positive solution satisfying (1.3) for any ν>0 if μ1, provided h ¯ 0 >0. Also, the equation

u + 1 u ν = h 0 (t) u δ

subjected to the boundary conditions (1.3) is solvable for any ν>0 if δ[1/2,1), provided h ¯ 0 >0.

Example 2.1 As it was mentioned in the introduction, the particular case of (1.1) is the so-called Rayleigh-Plesset equation frequently used in fluid mechanics. This equation has the following form:

u + c u u 4 5 + g 1 u 1 5 g 2 u γ = h 0 (t) u 1 5 ,
(2.14)

where h 0 L([0,ω];R), c, g 1 , g 2 are positive constants and γR (see [9, 10]).

The results dealing with the existence of positive ω-periodic solutions of (2.14) were established in [10] provided h 0 is bounded from above (see [[10], Theorems 4.4, 4.6, 4.7]). However, Corollary 2.1 says that in the case when γ<1/5, the problem (2.14), (1.3) is solvable if one of the following items is fulfilled:

  1. 1.

    γ>1/5 and h ¯ 0 >0;

  2. 2.

    γ=1/5 and h ¯ 0 > g 2 ;

  3. 3.

    γ<1/5.

Thus, Corollary 2.1 assures that the boundedness of h 0 is not necessary.

Corollary 2.2 Let g 1 >0, g 2 0, ν>γ, ν+1>0. Let, moreover, either g = or

h ¯ 0 > g >, ω 4 Φ + Φ < Φ + Φ ,
(2.15)

where φ(t)= h 0 (t) g for almost every t[0,ω], and

g = lim x + g 2 x γ + 1 .

Then the problem (1.1), (1.3) with δ=1 has at least one solution.

Remark 2.3 According to [29] and Theorem 1.1, it can be easily verified that the problem

u + g 1 u ν = h 0 (t)u;u(0)=u(ω), u (0)= u (ω),
(2.16)

with g 1 >0 and ν>0, has a positive solution if and only if the inclusion L[0, h 0 ] V holds (see notation in [29]).

Indeed, according to [[29], Definition 1.1], the inclusion L[0, h 0 ] V implies the existence of a positive solution v to the problem

v = h 0 (t)v g 1 ;u(0)=u(ω), u (0)= u (ω).

Therefore there exist constants x>0 and y>0 such that x 1 + ν v(t) y 1 + ν for t[0,ω]. By setting

α(t) = def v ( t ) y ν ,β(t) = def v ( t ) x ν for t[0,ω],

one can easily realize that α and β are, respectively, lower and upper functions to (2.16) satisfying (1.5). Now, the existence of a positive solution to (2.16) follows from Theorem 1.1.

On the other hand, the existence of a positive solution to (2.16) implies the inclusion L[0, h 0 ] V (see [[29], Theorem 2.1]).

However, one of the optimal effective conditions guaranteeing such an inclusion is h 0 0 and

ω 4 0 ω [ h 0 ( s ) ] + ds 0 ω [ h 0 ( s ) ] ds 0 ω [ h 0 ( s ) ] + ds 0 ω [ h 0 ( s ) ] ds

(see [[29], Corollary 2.5]). Therefore, the condition (2.15) is natural in a certain sense.

When the right-hand side of equation (1.3) does not depend on u, i.e., when h(t,x) h 0 (t), then (1.3) has the form

u (t)+f ( u ( t ) ) u (t)+g ( u ( t ) ) = h 0 (t)for a.e. t[0,ω].
(2.17)

From Theorem 2.1, for equation (2.17) we get the following assertion.

Corollary 2.3 Let there exist x 0 >0 and c 0 R such that

g(x) c 0 < h ¯ 0 for x x 0

and let

lim x 0 + g(x)=+.

Let, moreover, either

0 1 ( [ f ( s ) ] + + [ g ( s ) ] + ) ds=+, 0 1 [ f ( s ) ] ds<+

or

0 1 ( [ f ( s ) ] + [ g ( s ) ] + ) ds=+, 0 1 [ f ( s ) ] + ds<+.

Then there exists at least one solution to the problem (2.17), (1.3).

In the following result, the assumptions do not depend on the friction-like term. On the other hand, a certain smallness of oscillation of the primitive to h 0 is supposed. Clearly, Theorems 2.1 and 2.2 are independent.

Theorem 2.2 Let ρ 0 C 1 ( R + ; R + ) and ρ 1 C( R + ; R + ) be non-decreasing functions, h 0 , h 1 L([0,ω];R), and 0< x 0 x 1 <+ be such that

h(t,x) h 0 (t) ρ 0 (x) for a.e. t[0,ω],0<x x 0 ,
(2.18)
h(t,x) h 1 (t) ρ 1 (x) for a.e. t[0,ω],x x 1 .
(2.19)

Let, moreover,

ω 8 h 0 h ¯ 0 1 < 0 x 0 d s ρ 0 ( s ) <+,
(2.20)
g ( x ) ρ 0 ( x ) h ¯ 0 for 0<x x 0 ,
(2.21)

and let there exist c 1 R such that

g ( x ) ρ 1 ( x ) c 1 h ¯ 1 for x x 1 .
(2.22)

Besides, let us suppose that ρ 1 fulfills at least one of the conditions (c) or (d) of Theorem  2.1. Then there exists at least one solution to the problem (1.2), (1.3).

In the particular case, when equation (1.2) has the form (1.1), the following assertion immediately follows from Theorem 2.2.

Corollary 2.4 Let 0δ<1, and let 0< x 0 x 1 <+ be such that

(1δ) ω 8 h 0 h ¯ 0 1 < x 0 1 δ .

Let, moreover,

g 1 x ν + δ g 2 x γ + δ h ¯ 0 if 0 < x x 0 , g 1 x ν + δ g 2 x γ + δ h ¯ 0 if x x 1 .

Then the problem (1.1), (1.3) has at least one solution.

Remark 2.4 The consequence of Theorem 2.2 for the problem (2.17), (1.3) coincides with the result obtained in [[10], Theorem 3.6].

3 Auxiliary propositions

In what follows, we will show the existence of a solution to the equation

u (t)+f ( u ( t ) ) u (t)+g ( u ( t ) ) = h 0 (t) ρ 0 ( u ( t ) ) for a.e. t[0,ω],
(3.1)

satisfying the boundary conditions (1.3). Here, ρ 0 C( R + ; R + ) is a non-decreasing function, h 0 L([0,ω];R), and f,gC( R + ;R). Together with (3.1), for every kN, consider the auxiliary equation

u (t)+f ( u ( t ) ) u (t)+g ( u ( t ) ) = h 0 k (t) ρ 0 ( u ( t ) ) for a.e. t[0,ω],
(3.2)

where

h 0 k (t)={ k if  h 0 ( t ) > k , h 0 ( t ) if  h 0 ( t ) k for a.e. t[0,ω],kN.
(3.3)

Obviously,

h 0 k (t) h 0 m (t) h 0 (t)for a.e. t[0,ω],km,
(3.4)

and

lim k + h ¯ 0 k = h ¯ 0 .
(3.5)

The following three results can be found in [10].

Lemma 3.1 (see [[10], Corollary 2.17])

Let x 0 >0 and cR be such that

g ( x ) ρ 0 ( x ) c h ¯ 0 for x x 0 .
(3.6)

Let, moreover, there exist a sequence { y n } n = 1 + of positive numbers such that

lim n + y n =+, lim n + ρ 0 ( y n ) y n =0,
(3.7)

and let there exist ε>0 and n 0 N such that

ρ 0 ( ( 1 + ε ) y n ) ρ 0 ( y n ) Φ Φ + for n n 0 ,

where φ(t)= h 0 (t)c for almost every t[0,ω]. Then there exists an upper function β to the problem (3.1), (1.3) satisfying

β(t) x 0 for t[0,ω].
(3.8)

Lemma 3.2 (see [[10], Corollary 2.18])

Let x 0 >0 and cR be such that (3.6) holds. If ρ 0 ( x ) x is a non-increasing function such that

ω 4 Φ + Φ ρ 0 ( x 0 ) x 0 Φ + Φ ,

where φ(t)= h 0 (t)c for almost every t[0,ω], then there exists an upper function β to the problem (3.1), (1.3) satisfying (3.8).

Lemma 3.3 (see [[10], Corollary 2.11])

Let x 0 > ω 8 h 0 h ¯ 0 1 be such that

g(x) h ¯ 0 for 0<x x 0 .

Then there exists a lower function α to the problem (2.17), (1.3) with

0<α(t) x 0 for t[0,ω].
(3.9)

Lemma 3.4 Let x 0 >0 and c 0 R be such that (2.2) holds. Let, moreover, there exist a sequence { y n } n = 1 + of positive numbers such that (3.7) is fulfilled, and let there exist ε>0 and n 0 N such that

ρ 0 ( ( 1 + ε ) y n ) ρ 0 ( y n ) Φ Φ + ε for n n 0 ,
(3.10)

where φ(t)= h 0 (t) c 0 for almost every t[0,ω]. Then there exist k 0 N and an upper function β to the problems (3.2), (1.3) for k k 0 satisfying (3.8).

Proof

Put

φ k (t)= h 0 k (t) c 0 for a.e. t[0,ω],
(3.11)
Φ k + = 0 ω [ φ k ( s ) ] + ds, Φ k = 0 ω [ φ k ( s ) ] ds.
(3.12)

Then, obviously, in view of (3.3), we have

lim k + Φ k + = Φ + , lim k + Φ k = Φ
(3.13)

and, consequently, on account of (2.2), (3.5), (3.10), and (3.13), there exists k 0 N such that

g ( x ) ρ 0 ( x ) c 0 h ¯ 0 k 0 h ¯ 0 for x x 0 ,
(3.14)
ρ 0 ( ( 1 + ε ) y n ) ρ 0 ( y n ) Φ k 0 Φ k 0 + for n n 0 .
(3.15)

Therefore, according to Lemma 3.1, there exists an upper function β to (3.2), (1.3) with k= k 0 satisfying (3.8). Obviously, in view of (3.4) and the non-negativity of ρ 0 , it follows that β is also an upper function to (3.2), (1.3) for k k 0 . □

Lemma 3.5 Let x 0 >0 and c 0 R be such that (2.2) holds. If ρ 0 ( x ) x is a non-increasing function such that

ω 4 Φ + Φ ρ 0 ( x 0 ) x 0 < Φ + Φ ,
(3.16)

where φ(t)= h 0 (t) c 0 for almost every t[0,ω], then there exist k 0 N and an upper function β to the problems (3.2), (1.3) for k k 0 satisfying (3.8).

Proof Define φ k , Φ k + , and Φ k by (3.11) and (3.12). Then, obviously, in view of (3.3), we have that (3.13) holds and, consequently, on account of (2.2), (3.5), (3.13), and (3.16), there exists k 0 N such that (3.14) is valid and

ω 4 Φ k 0 + Φ k 0 ρ 0 ( x 0 ) x 0 Φ k 0 + Φ k 0 .
(3.17)

Therefore, according to Lemma 3.2, there exists an upper function β to (3.2), (1.3) with k= k 0 satisfying (3.8). Obviously, in view of (3.4) and the non-negativity of ρ 0 , it follows that β is also an upper function to (3.2), (1.3) for k k 0 . □

Lemma 3.6 Let

lim inf x 0 + g(x)>,
(3.18)

and let either

0 1 [ f ( s ) ] + ds<+
(3.19)

or

0 1 [ f ( s ) ] ds<+.
(3.20)

Then, for every K>0, there exists a constant K 1 >0 such that for any kN and any positive solution u of (3.2), (1.3) with

u K,
(3.21)

we have the estimate

u K 1 .
(3.22)

Proof Assume that (3.20) is fulfilled. Let u be a positive solution to (3.2), (1.3) satisfying (3.21). Then there exist t 0 , t 1 [0,ω] such that

u( t 0 )=min { u ( t ) : t [ 0 , ω ] } ,u( t 1 )=max { u ( t ) : t [ 0 , ω ] } .
(3.23)

Define the operator ϑ of ω-periodic prolongation by

ϑ(v)(t)={ v ( t ) if  t [ 0 , ω ] , v ( t ω ) if  t ( ω , 2 ω ] .
(3.24)

Then, obviously, from (3.2) and (1.3) it follows that

ϑ ( u ) ( t ) + f ( ϑ ( u ) ( t ) ) ϑ ( u ) ( t ) + g ( ϑ ( u ) ( t ) ) = ϑ ( h 0 k ) ( t ) ρ 0 ( ϑ ( u ) ( t ) ) for a.e.  t [ 0 , 2 ω ] .
(3.25)

The integration of (3.25) from t 0 to t, on account of (3.23), yields

ϑ ( u ) ( t ) = t 0 t f ( ϑ ( u ) ( s ) ) ϑ ( u ) ( s ) d s t 0 t g ( ϑ ( u ) ( s ) ) d s + t 0 t ϑ ( h 0 k ) ( s ) ρ 0 ( ϑ ( u ) ( s ) ) d s for  t [ t 0 , t 0 + ω ] .
(3.26)

From (3.21), (3.23), and (3.24) it follows that

0<ϑ(u)( t 0 )ϑ(u)(t)Kfor t[ t 0 , t 0 +ω].
(3.27)

Put

μ=sup { [ g ( s ) ] : s ( 0 , K ] } .
(3.28)

According to (3.18), we have

0μ<+.
(3.29)

Thus, using (3.3), (3.20), (3.21), and (3.27)-(3.29) in (3.26), we arrive at

ϑ ( u ) (t) 0 K [ f ( s ) ] ds+ωμ+ h 0 1 ρ 0 (K)for t[ t 0 , t 0 +ω].
(3.30)

Put

K 1 = 0 K [ f ( s ) ] ds+ωμ+ h 0 1 ρ 0 (K).

Then, on account of (3.24) and (3.30), we have

u (t) K 1 for t[0,ω].
(3.31)

On the other hand, the integration of (3.25) from t to t 1 +ω, with respect to (3.23), results in

ϑ ( u ) ( t ) = t t 1 + ω f ( ϑ ( u ) ( s ) ) ϑ ( u ) ( s ) d s + t t 1 + ω g ( ϑ ( u ) ( s ) ) d s t t 1 + ω ϑ ( h 0 k ) ( s ) ρ 0 ( ϑ ( u ) ( s ) ) d s for  t [ t 1 , t 1 + ω ] .
(3.32)

Now, using (3.3), (3.20), (3.21), and (3.27)-(3.29) in (3.32), we obtain

ϑ ( u ) (t) K 1 for t[ t 1 , t 1 +ω].
(3.33)

Therefore, in view of (3.24), from (3.33) we get

u (t) K 1 for t[0,ω].
(3.34)

Consequently, (3.31) and (3.34) result in (3.22).

Now, suppose that (3.19) is fulfilled. Put

v(t)=u(ωt)for t[0,ω].
(3.35)

Then, according to (3.2), we have

v (t)f ( v ( t ) ) v (t)+g ( v ( t ) ) = h ˜ 0 k (t) ρ 0 ( v ( t ) ) for a.e. t[0,ω],
(3.36)

where

h ˜ 0 k (t)= h 0 k (ωt)for a.e. t[0,ω].

Analogously to the above-proved, using (3.19) instead of (3.20), we obtain

v K 1
(3.37)

with

K 1 = 0 K [ f ( s ) ] + ds+ωμ+ h 0 1 ρ 0 (K).

Thus, (3.35) and (3.37) yield (3.22). □

Lemma 3.7 Let

lim x 0 + g(x)=+
(3.38)

and let either

0 1 ( [ f ( s ) ] + + [ g ( s ) ] + ) ds=+, 0 1 [ f ( s ) ] ds<+
(3.39)

or

0 1 ( [ f ( s ) ] + [ g ( s ) ] + ) ds=+, 0 1 [ f ( s ) ] + ds<+.
(3.40)

Then, for every K>0, there exists a constant a>0 such that for any kN and any positive solution u of (3.2), (1.3) satisfying (3.21), we have the estimate

au(t) for t[0,ω].
(3.41)

Proof Let u be a positive solution to (3.2), (1.3) satisfying (3.21). Thus, the integration of (3.2) from 0 to ω, in view of (1.3) and (3.4), yields

0 ω g ( u ( s ) ) ds h 0 1 ρ 0 (K).
(3.42)

On the other hand, (3.38) implies the existence of x 0 (0,+) such that

g(x)> h 0 1 ρ 0 ( K ) ω 0for x(0, x 0 ).
(3.43)

Let t m [0,ω] be such that

u( t m )=min { u ( t ) : t [ 0 , ω ] } .
(3.44)

Obviously, either

u( t m ) x 0

or

u( t m )< x 0 .
(3.45)

Obviously, it is sufficient to show the estimate (3.41) is valid just in the case when (3.45) is fulfilled. Let, therefore, (3.45) hold.

If u(t)< x 0 for t[0,ω], then applying (3.43) in (3.42) we obtain a contradiction. Thus, there exist points t 1 , t 2 ( t m , t m +ω) such that

ϑ(u)(t)< x 0 for t[ t m , t 1 ),ϑ(u)( t 1 )= x 0 ,
(3.46)
ϑ(u)(t)< x 0 for t( t 2 , t m +ω],ϑ(u)( t 2 )= x 0 ,
(3.47)

where ϑ is the operator defined by (3.24). Obviously, (3.25) holds.

Assume that (3.39) holds. Then, according to Lemma 3.6, there exists K 1 >0 such that (3.22) holds. The integration of (3.25) from t m to t 1 , in view of (3.4), (3.21), (3.22), (3.39), (3.43), (3.44), and (3.46), results in

ϑ ( u ) ( t 1 ) + ϑ ( u ) ( t m ) x 0 [ f ( s ) ] + d s + 1 K 1 ϑ ( u ) ( t m ) x 0 [ g ( s ) ] + d s 0 x 0 [ f ( s ) ] d s + h 0 1 ρ 0 ( K ) .
(3.48)

Note that in view of (3.46), we have ϑ ( u ) ( t 1 )0. Consequently, from (3.48) we obtain

ϑ ( u ) ( t m ) x 0 ( [ f ( s ) ] + + [ g ( s ) ] + ) ds K 2 ,
(3.49)

where

K 2 =( K 1 +1) ( 0 x 0 [ f ( s ) ] d s + h 0 1 ρ 0 ( K ) ) .

Note that K 2 does not depend on k. Therefore, if we apply (3.39) in (3.49), it can be easily seen, with respect to (3.44), that there exists a constant a>0 such that (3.41) holds.

If (3.40) holds, we integrate (3.25) from t 2 to t m +ω and apply similar steps as above, just using (3.47) instead of (3.46). Finally, we arrive at

ϑ ( u ) ( t m + ω ) x 0 ( [ f ( s ) ] + [ g ( s ) ] + ) ds K 2

with

K 2 =( K 1 +1) ( 0 x 0 [ f ( s ) ] + d s + h 0 1 ρ 0 ( K ) ) .

Therefore, also in this case, there exists a constant a>0 such that (3.41) holds. □

Lemma 3.8 Let x 0 >0 and c 0 R be such that (2.2) holds. Let, moreover, (3.38) be fulfilled, and let either (3.39) or (3.40) be valid. Let, in addition, there exist a sequence { y n } n = 1 + of positive numbers such that (3.7) holds, and let there exist ε>0 and n 0 N such that (3.10) is fulfilled, where φ(t)= h 0 (t) c 0 for almost every t[0,ω]. Then there exists a positive solution u to (3.1), (1.3).

Proof According to Lemma 3.4, there exist k 0 N and an upper function β to the problems (3.2), (1.3) for k k 0 satisfying (3.8). On the other hand, in view of (3.4) and (3.38), there exists x k (0, x 0 ] for k k 0 such that

g( x k ) h 0 k (t) ρ 0 ( x k )for a.e. t[0,ω].

Thus, if we put α k (t)= x k for t[0,ω], according to Theorem 1.1, there exists a solution u k to (3.2), (1.3) for k k 0 satisfying

0< α k (t) u k (t)β(t)for t[0,ω].
(3.50)

Moreover, according to Lemmas 3.6 and 3.7, in view of (3.50), there exist constants K>0, K 1 >0, and a>0, not depending on k, such that

u k K, u k K 1 for k k 0 ,
(3.51)
a u k (t)for t[0,ω],k k 0 ,
(3.52)
| u k ( t ) | f 0 K 1 + g 0 + | h 0 ( t ) | ρ 0 (K)for a.e. t[0,ω],k k 0 ,
(3.53)

where

f 0 =max { | f ( x ) | : x [ a , K ] } , g 0 =max { | g ( x ) | : x [ a , K ] } .

Therefore, according to the Arzelà-Ascoli theorem, there exist u 0 C([0,ω];R) and v 0 C([0,ω];R) such that

lim k + u k u 0 =0, lim k + u k v 0 =0.
(3.54)

Moreover, since u k are solutions to (3.2), (1.3), in view of (3.3), (3.52), and (3.54), we have u 0 A C 1 ([0,ω];R), u 0 v 0 , and u 0 is a positive solution to (3.1), (1.3). □

The following assertion can be proved analogously to Lemma 3.8, just Lemma 3.5 is used instead of Lemma 3.4.

Lemma 3.9 Let x 0 >0 and c 0 R be such that (2.2) holds. Let, moreover, (3.38) be fulfilled, and let either (3.39) or (3.40) be valid. Let, in addition, ρ 0 ( x ) x be a non-increasing function and let (3.16) be fulfilled, where φ(t)= h 0 (t) c 0 for almost every t[0,ω]. Then there exists a positive solution u to (3.1), (1.3).

Lemma 3.10 Let ρ 0 C 1 ( R + ; R + ) be non-decreasing, x 0 >0, and c 0 R be such that (2.2) holds. Let, moreover, there exist λ[0,1] such that (2.4) and (2.5) are valid, and let either (2.6) or (2.7) be fulfilled. Let, in addition, there exist a sequence { y n } n = 1 + of positive numbers such that (2.8) holds and let there exist ε 0 >0, ε 1 (0, ε 0 ], and n 0 N such that (2.9) and (2.10) are fulfilled, where φ(t)= h 0 (t) c 0 for almost every t[0,ω] and σ is given by (2.11). Then there exists a lower function α to the problem (3.1), (1.3).

Proof Because ρ 0 is a positive function, from (2.4) and (2.11) we obtain that σ is a positive increasing function. Therefore, there exists an inverse function σ 1 to σ which is also increasing. Moreover, in view of (2.4) and (2.11), it follows that

lim x 0 + σ ( x ) = 0 , lim x 0 + σ 1 ( x ) = 0 , lim x + σ ( x ) = + , lim x + σ 1 ( x ) = + .
(3.55)

Consider the auxiliary equation

u ( t ) + f ( σ 1 ( u ( t ) ) ) u ( t ) + g ( σ 1 ( u ( t ) ) ) ρ 0 λ ( σ 1 ( u ( t ) ) ) = h 0 ( t ) ρ 0 1 λ ( σ 1 ( u ( t ) ) ) for a.e.  t [ 0 , ω ] .
(3.56)

Put z=σ(x), z 0 =σ( x 0 ). Then from (2.2) we get

g ( σ 1 ( z ) ) ρ 0 ( σ 1 ( z ) ) c 0 < h ¯ 0 for z z 0
(3.57)

and, in view of (3.55), from (2.5) we have

lim z 0 + g ( σ 1 ( z ) ) ρ 0 λ ( σ 1 ( z ) ) =+.
(3.58)

Furthermore, the substitution r=σ(s) in (2.6), resp (2.7), with respect to (2.11) yields

0 1 ( [ f ( σ 1 ( r ) ) ] + + [ g ( σ 1 ( r ) ) ] + ρ 0 λ ( σ 1 ( r ) ) ) dr=+, 0 1 [ f ( σ 1 ( r ) ) ] dr<+,
(3.59)

resp.

0 1 ( [ f ( σ 1 ( r ) ) ] + [ g ( σ 1 ( r ) ) ] + ρ 0 λ ( σ 1 ( r ) ) ) dr=+, 0 1 [ f ( σ 1 ( r ) ) ] + dr<+.
(3.60)

Moreover, put z n =σ( y n ) for nN. Then from (2.8), in view of (3.55), we get

lim n + z n =+, lim n + ρ 0 1 λ ( σ 1 ( z n ) ) z n =0.
(3.61)

Finally, (2.10) results in

σ 1 ( ( 1 + ε 1 ) z n ) (1+ ε 0 ) y n for n n 0 ,

and so, since ρ 0 is a non-decreasing function, from (2.9) we obtain

ρ 0 1 λ ( σ 1 ( ( 1 + ε 1 ) z n ) ) ρ 0 1 λ ( σ 1 ( z n ) ) Φ Φ + ε 1 for n n 0 .
(3.62)

Therefore, applying Lemma 3.8, according to (3.57)-(3.62), there exists a positive solution u to the problem (3.56), (1.3).

Now, we put α(t)= σ 1 (u(t)) for t[0,ω], i.e., in view of (2.11),

u(t)= 0 α ( t ) d s ρ 0 λ ( s ) for t[0,ω].

Obviously, αA C 1 ([0,ω];R) is a positive function and

u ( t ) = α ( t ) ρ 0 λ ( α ( t ) ) for  t [ 0 , ω ] , u ( t ) = α ( t ) ρ 0 λ ( α ( t ) ) λ α 2 ( t ) ρ 0 ( α ( t ) ) ρ 0 1 + λ ( α ( t ) ) α ( t ) ρ 0 λ ( α ( t ) ) for a.e.  t [ 0 , ω ] .

Thus, it can be easily seen that α is a lower function to the problem (3.1), (1.3). □

Analogously to the proof of Lemma 3.10, one can prove the following assertion applying Lemma 3.9 instead of Lemma 3.8.

Lemma 3.11 Let ρ 0 C 1 ( R + ; R + ) be non-decreasing, x 0 >0, and c 0 R be such that (2.2) holds. Let, moreover, there exist λ[0,1] such that (2.4) and (2.5) are valid, and let either (2.6) or (2.7) be fulfilled. Let, in addition, ρ 0 1 λ ( x ) σ ( x ) be a non-increasing function and let (2.12) be fulfilled, where φ(t)= h 0 (t) c 0 for almost every t[0,ω] and σ is given by (2.11). Then there exists a lower function α to the problem (3.1), (1.3).

4 Proofs of the main results

Proof of Theorem 2.1 According to Lemmas 3.1, 3.2, 3.10, and 3.11, the conditions of the theorem guarantee a well-ordered couple of lower and upper functions, therefore the result is a direct consequence of Theorem 1.1. □

Proof of Corollary 2.1 It follows from Theorem 2.1 with h 1 h 0 , ρ 0 (x)= ρ 1 (x)= x δ , λ=1, and c 0 = c 1 such that

h ¯ 0 > c 0 > lim x + g 2 x γ + δ .

Then items (a) and (c) of Theorem 2.1 are fulfilled. □

Proof of Corollary 2.2 It follows from Theorem 2.1 with h 1 h 0 , ρ 0 (x)= ρ 1 (x)=x, and λ<1 such that ν+λ>0, ν+2λ1. Then items (b) and (d) of Theorem 2.1 are fulfilled. □

Proof of Corollary 2.3 It immediately follows from Theorem 2.1 with h 1 h 0 , ρ i (x)1 (i=0,1). □

Proof of Theorem 2.2

Put

σ(x)= 0 x d s ρ 0 ( s ) for x0.
(4.1)

Because ρ 0 is a positive function, from (2.20) and (4.1) we obtain that σ is an increasing function. Therefore, there exists an inverse function σ 1 to σ which is also increasing.

Consider the auxiliary equation

u (t)+f ( σ 1 ( u ( t ) ) ) u (t)+ g ( σ 1 ( u ( t ) ) ) ρ 0 ( σ 1 ( u ( t ) ) ) = h 0 (t)for a.e. t[0,ω].
(4.2)

Put z=σ(x), z 0 =σ( x 0 ). Then from (2.20) and (2.21), in view of (4.1), we get

ω 8 h 0 h ¯ 0 1 < z 0 , g ( σ 1 ( z ) ) ρ 0 ( σ 1 ( z ) ) h ¯ 0 for  0 < z z 0 .

Therefore, according to Lemma 3.3, there exists a lower function w to the problem (4.2), (1.3) satisfying

0<w(t)σ( x 0 )for t[0,ω].
(4.3)

Now, we put α= σ 1 (w(t)) for t[0,ω], i.e., in view of (4.1),

w(t)= 0 α ( t ) d s ρ 0 ( s ) for t[0,ω].

Obviously, with respect to (4.3), αA C 1 ([0,ω];R) is a positive function satisfying (3.9), and

w ( t ) = α ( t ) ρ 0 ( α ( t ) ) for  t [ 0 , ω ] , w ( t ) = α ( t ) ρ 0 ( α ( t ) ) α 2 ( t ) ρ 0 ( α ( t ) ) ρ 0 2 ( α ( t ) ) α ( t ) ρ 0 ( α ( t ) ) for a.e.  t [ 0 , ω ] .

Thus, on account of (2.18), (3.9), and (4.2), it can be easily seen that α is a lower function to the problem (1.2), (1.3).

The existence of an upper function β to (1.2), (1.3) satisfying

β(t) x 1 for t[0,ω]
(4.4)

follows from (2.19) and Lemma 3.1, resp. 3.2.

Obviously, in view of (3.9) and (4.4), we have that (1.5) holds. Thus the theorem follows from Theorem 1.1. □

Proof of Corollary 2.4 It follows from Theorem 2.2 with h 1 h 0 , ρ 0 (x)= ρ 1 (x)= x δ , and c 1 = h ¯ 0 . □

References

  1. Plesset MS, Prosperetti A: Bubble dynamics and cavitation. Annu. Rev. Fluid Mech. 1977, 9: 145-185. 10.1146/annurev.fl.09.010177.001045

    Article  Google Scholar 

  2. Habets P, Sanchez L: Periodic solutions of some Liénard equations with singularities. Proc. Am. Math. Soc. 1990, 109: 1135-1144. 10.1090/S0002-9939-1990-1009992-7

    Article  MathSciNet  Google Scholar 

  3. Bonheure D, Fabry C, Smets D: Periodic solutions of forced isochronous oscillators at resonance. Discrete Contin. Dyn. Syst. 2002, 8(4):907-930.

    Article  MathSciNet  Google Scholar 

  4. Bonheure D, De Coster C: Forced singular oscillators and the method of upper and lower solutions. Topol. Methods Nonlinear Anal. 2003, 22: 297-317.

    MathSciNet  Google Scholar 

  5. Martins RF: Existence of periodic solutions for second-order differential equations with singularities and the strong force condition. J. Math. Anal. Appl. 2006, 317: 1-13. 10.1016/j.jmaa.2004.07.016

    Article  MathSciNet  Google Scholar 

  6. Mawhin J: Topological degree and boundary value problems for nonlinear differential equations. Lecture Notes in Mathematics 1537. In Topological Methods for Ordinary Differential Equations. Edited by: Furi M, Zecca P. Springer, Berlin; 1993:74-142.

    Chapter  Google Scholar 

  7. Omari P, Ye W: Necessary and sufficient conditions for the existence of periodic solutions of second-order ordinary differential equations with singular nonlinearities. Differ. Integral Equ. 1995, 8: 1843-1858.

    MathSciNet  Google Scholar 

  8. Yuan R, Zhang Z: Existence of positive periodic solutions for the Liénard differential equations with weakly repulsive singularity. Acta Appl. Math. 2010, 111: 171-178. 10.1007/s10440-009-9538-x

    Article  MathSciNet  Google Scholar 

  9. Hakl R, Torres P, Zamora M: Periodic solutions of singular second order differential equations: the repulsive case. Topol. Methods Nonlinear Anal. 2012, 39: 199-220.

    MathSciNet  Google Scholar 

  10. Hakl R, Torres P, Zamora M: Periodic solutions of singular second order differential equations: upper and lower functions. Nonlinear Anal. 2011, 74: 7078-7093. 10.1016/j.na.2011.07.029

    Article  MathSciNet  Google Scholar 

  11. Lazer AC, Solimini S: On periodic solutions of nonlinear differential equations with singularities. Proc. Am. Math. Soc. 1987, 99: 109-114. 10.1090/S0002-9939-1987-0866438-7

    Article  MathSciNet  Google Scholar 

  12. Rachůnková I, Staněk S, Tvrdý M Contemporary Mathematics and Its Applications 5. In Solvability of Nonlinear Singular Problems for Ordinary Differential Equations. Hindawi Publishing Corporation, New York; 2008.

    Google Scholar 

  13. Zhang M: Periodic solutions of Liénard equations with singular forces of repulsive type. J. Math. Anal. Appl. 1996, 203: 254-269. 10.1006/jmaa.1996.0378

    Article  MathSciNet  Google Scholar 

  14. Chu J, Li M: Positive periodic solutions of Hill’s equation with singular nonlinear perturbations. Nonlinear Anal. 2008, 69: 276-286. 10.1016/j.na.2007.05.016

    Article  MathSciNet  Google Scholar 

  15. Chu J, Nieto JJ: Recent existence results for second-order singular periodic differential equations. Bound. Value Probl. 2009., 2009: Article ID 540863

    Google Scholar 

  16. Chu J, Torres PJ: Applications of Schauder’s fixed point theorem to singular differential equations. Bull. Lond. Math. Soc. 2007, 39: 653-660. 10.1112/blms/bdm040

    Article  MathSciNet  Google Scholar 

  17. Chu J, Torres PJ, Zhang M: Periodic solutions of second order non-autonomous singular dynamical systems. J. Differ. Equ. 2007, 239: 196-212. 10.1016/j.jde.2007.05.007

    Article  MathSciNet  Google Scholar 

  18. Franco D, Torres PJ: Periodic solutions of singular systems without the strong force condition. Proc. Am. Math. Soc. 2008, 136: 1229-1236.

    Article  MathSciNet  Google Scholar 

  19. Franco D, Webb JRL: Collisionless orbits of singular and nonsingular dynamical systems. Discrete Contin. Dyn. Syst. 2006, 15: 747-757.

    Article  MathSciNet  Google Scholar 

  20. Hakl R, Torres P: On periodic solutions of second-order differential equations with attractive-repulsive singularities. J. Differ. Equ. 2010, 248: 111-126. 10.1016/j.jde.2009.07.008

    Article  MathSciNet  Google Scholar 

  21. Rachůnková I, Tvrdý M, Vrkoč I: Existence of nonnegative and nonpositive solutions for second order periodic boundary value problems. J. Differ. Equ. 2001, 176: 445-469. 10.1006/jdeq.2000.3995

    Article  Google Scholar 

  22. Torres PJ: Existence of one-signed periodic solutions of some second order differential equations via a Kranoselskii fixed point theorem. J. Differ. Equ. 2003, 190: 643-662. 10.1016/S0022-0396(02)00152-3

    Article  Google Scholar 

  23. Torres PJ: Non-collision periodic solutions of forced dynamical systems with weak singularities. Discrete Contin. Dyn. Syst. 2004, 11: 693-698.

    Article  MathSciNet  Google Scholar 

  24. Torres PJ: Weak singularities may help periodic solutions to exists. J. Differ. Equ. 2007, 232: 277-284. 10.1016/j.jde.2006.08.006

    Article  Google Scholar 

  25. Torres PJ: Existence and stability of periodic solutions for second order semilinear differential equations with a singular nonlinearity. Proc. R. Soc. Edinb., Sect. A, Math. 2007, 137: 195-201.

    Article  Google Scholar 

  26. Yan P, Zhang M: Higher order nonresonance for differential equations with singularities. Math. Methods Appl. Sci. 2003, 26: 1067-1074. 10.1002/mma.413

    Article  MathSciNet  Google Scholar 

  27. Zhang M: Periodic solutions of equations of Emarkov-Pinney type. Adv. Nonlinear Stud. 2006, 6: 57-67.

    MathSciNet  Google Scholar 

  28. Hakl, R, Zamora, M: On the open problems connected to the results of Lazer and Solimini. Proc. R. Soc. Edinb., Sect. A, Math. (to appear) http://www.math.cas.cz/fichier/preprints/IM_20120709132356_18.pdf

  29. Hakl R, Torres P: Maximum and antimaximum principles for a second order differential operator with variable coefficients of indefinite sign. Appl. Math. Comput. 2011, 217: 7599-7611. 10.1016/j.amc.2011.02.053

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

The first author was supported by RVO: 67985840; the second author was supported by Ministerio de Educación y Ciencia, Spain, project MTM2011-23652.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Robert Hakl.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RH and MZ obtained the results in a joint research. Both authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Hakl, R., Zamora, M. Periodic solutions to the Liénard type equations with phase attractive singularities. Bound Value Probl 2013, 47 (2013). https://doi.org/10.1186/1687-2770-2013-47

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-47

Keywords