## A TOPOLOGICAL APPROACH TO PERIODIC OSCILLATIONS RELATED TO THE LIEBAU PHENOMENON

J. A. CID, G. INFANTE, M. TVRD ´ , AND M. ZIMA Y

Abstract. We give some sufficient conditions for existence, non-existence and localization of positive solutions for a periodic boundary value problem related to the Liebau phenomenon. Our approach is of topological nature and relies on the Krasnosel'ski˘ ı-Guo theorem on cone expansion and compression. Our results improve and complement earlier ones in the literature.

## 1. Introduction

In the 1950's the physician G. Liebau developed some experiments dealing with a valveless pumping phenomenon arising on blood circulation and that has been known for a long time: roughly speaking, Liebau showed experimentally that a periodic compression made on an asymmetric part of a fluid-mechanical model could produce the circulation of the fluid without the necessity of a valve to ensure a preferential direction of the flow [1, 12, 13]. After his pioneering work this effect has been known as the Liebau phenomenon.

In [14] G. Propst, with the aim of contributing to the theoretical understanding of the Liebau phenomenon, presented some differential equations modeling a periodically forced flow through different pipe-tank configurations. He was able to prove the presence of pumping effects in some of them, but the apparently simplest model, the 'one pipe-one tank' configuration, skipped his efforts due to a singularity in the corresponding differential equation model, namely

$$\overbrace { \overbrace { \infty } } ^ { \overbrace { \sim } } \quad ( 1 ) & \quad \begin{cases} \ u ^ { \prime \prime } ( t ) + a \, u ^ { \prime } ( t ) = \frac { 1 } { u } \left ( e ( t ) - b \left ( u ^ { \prime } ( t ) \right ) ^ { 2 } \right ) - c, & t \in [ 0, T ], \\ u ( 0 ) = u ( T ), & u ^ { \prime } ( 0 ) = u ^ { \prime } ( T ), \end{cases}$$

being a ≥ 0 , b &gt; 1 , c &gt; 0 and e t ( ) continuous and T -periodic on R .

The singular periodic problem (1) was studied in [4], where the authors gave general results for the existence and asymptotic stability of positive solutions by performing the change of variables u = x µ , where µ = 1 b +1 , which transforms the singular problem (1) into the regular

2010 Mathematics Subject Classification. Primary 34B18, secondary 34B27, 34B60.

Key words and phrases. Valveless pumping; Periodic boundary value problem; Krasnosel'ski˘ ı-Guo fixed point theorem; Green's function.

$$\begin{cases} \ x ^ { \prime \prime } ( t ) + a \, x ^ { \prime } ( t ) = \frac { e ( t ) } { \mu } \, x ^ { 1 - 2 \, \mu } ( t ) - \frac { c } { \mu } \, x ^ { 1 - \mu } ( t ), \quad t \in [ 0, T ], \\ x ( 0 ) = x ( T ), \quad x ^ { \prime } ( 0 ) = x ^ { \prime } ( T ). \end{cases}$$

Then the existence and stability of positive solutions for (2) were obtained by means of the lower and upper solution technique [5] and using tricks analogous to those used in [15].

In this paper we deal with the existence of positive solutions for the following generalization of the problem (2)

$$\begin{cases} \ x ^ { \prime \prime } ( t ) + a x ^ { \prime } ( t ) = r ( t ) x ^ { \alpha } ( t ) - s ( t ) x ^ { \beta } ( t ), & t \in [ 0, T ], \\ \ x ( 0 ) = x ( T ), & x ^ { \prime } ( 0 ) = x ^ { \prime } ( T ), \end{cases}$$

where we assume

$$( \text{H0} ) \ a \geq 0, \, r, s \in \mathcal { C } [ 0, T ], \, 0 < \alpha < \beta < 1.$$

Note that, by defining r t ( ) = e t ( ) µ , s t ( ) = c µ , α = 1 -2 µ and β = 1 -µ , the problem (2) fits within (3).

Our approach is essentially of topological nature: in Section 2 we rewrite problem (3) as an equivalent fixed point problem suitable to be treated by means of the Krasnosel'ski˘ ı-Guo cone expansion/compression fixed point theorem. A careful analysis of the related Green's function, necessary in our approach, is postponed to a final Appendix. In Section 3 we present our main results: existence, non-existence and localization criteria for positive solutions of the problem (3). Some corollaries with more ready-to-use results are also addressed. We point out that our results are valid not only for the more general problem (3), but also when applied to the singular model problem (1) we improve previous results of [4].

## 2. A fixed point formulation

First of all, by means of a shifting argument, we rewrite the problem (3) in the equivalent form

$$( 4 ) \quad \begin{cases} \ x ^ { \prime \prime } ( t ) + a x ^ { \prime } ( t ) + m ^ { 2 } x ( t ) = r ( t ) x ^ { \alpha } ( t ) - s ( t ) x ^ { \beta } ( t ) + m ^ { 2 } x ( t ) \coloneqq f _ { m } ( t, x ( t ) ), \ t \in [ 0, T ], \\ x ( 0 ) = x ( T ), \ \ x ^ { \prime } ( 0 ) = x ^ { \prime } ( T ), \end{cases}$$

with m ∈ R . A similar approach has been used, under a variety of boundary conditions in [7, 9, 16, 17].

We say that problem (4) is non-resonant if zero is the unique solution of the homogeneous linear problem

$$\begin{cases} \ x ^ { \prime \prime } ( t ) + a x ^ { \prime } ( t ) + m ^ { 2 } x ( t ) = 0, \quad t \in [ 0, T ], \\ \ x ( 0 ) = x ( T ), \quad x ^ { \prime } ( 0 ) = x ^ { \prime } ( T ). \end{cases}$$

In this case the non-homogeneous linear problem

$$\begin{cases} \ x ^ { \prime \prime } ( t ) + a x ^ { \prime } ( t ) + m ^ { 2 } x ( t ) = h ( t ), & t \in [ 0, T ], \\ \ x ( 0 ) = x ( T ), & x ^ { \prime } ( 0 ) = x ^ { \prime } ( T ), \end{cases}$$

is also uniquely solvable and its unique solution is given by

$$\mathcal { K } h ( t ) = \int _ { 0 } ^ { T } G _ { m } ( t, s ) h ( s ) d s,$$

where G m ( t, s ) is the related Green's function (see [2, 3]).

In particular, if m&gt; 0 and m &lt; 2 ( π T ) 2 + ( a 2 ) 2 then the problem (4) is non-resonant and moreover G m ( t, s ) satisfies the following properties (see the Appendix for the details):

- i) G m ( t, s ) &gt; 0 , for all ( t, s ) ∈ [0 , T ] × [0 , T ].
- ii) ∫ T G m ( t, s ) ds = 1 m . 2
- 0 iii) There exists a constant c m ∈ (0 , 1) such that

$$G _ { m } ( t, s ) \geq G _ { m } ( s, s ) \geq c _ { m } G _ { m } ( t, s ), \text{ for all } ( t, s ) \in [ 0, T ] \times [ 0, T ].$$

A cone in a Banach space X is a closed, convex subset of X such that λx ∈ K for x ∈ K and λ ≥ 0 and K ∩ -( K ) = { } 0 . Here we work in the space X = C [0 , T ] endowed with the usual maximum norm ‖ x ‖ = max {| x t ( ) | : t ∈ [0 , T ] } and use the cone

$$P = \{ x \in X \, \colon x ( t ) \geq c _ { m } \| x \| \text{ on } [ 0, T ] \},$$

a type of cone firstly used by Krasnosel'ski˘ ı, see for example [10], and D. Guo, see for example [6]. The cone P is particularly useful when dealing with singular nonlinearities (see [11]), or for localizing the solutions (see for example [8]).

Let us define the operator F : P → X

$$F x ( t ) = \int _ { 0 } ^ { T } G _ { m } ( t, s ) f _ { m } ( s, x ( s ) ) d s.$$

Note that a fixed point of F in P is a non-negative solution of the problem (4). In order to get such a fixed point we employ the following well-known Krasnosel'ski˘ ı-Guo cone compression/expansion theorem.

Theorem 2.1. [6] Let P be a cone in X and suppose that Ω 1 and Ω 2 are bounded open sets in X such that 0 ∈ Ω 1 and Ω 1 ⊂ Ω 2 . Let F : P ∩ (Ω 2 \ Ω ) 1 → P be a completely continuous operator such that one of the following conditions holds:

- (i) ‖ Fx ‖ ≥ ‖ x ‖ for x ∈ P ∩ ∂ Ω 1 and ‖ Fx ‖ ≤ ‖ x ‖ for x ∈ P ∩ ∂ Ω 2 ,
- (ii) ‖ Fx ‖ ≤ ‖ x ‖ for x ∈ P ∩ ∂ Ω 1 and ‖ Fx ‖ ≥ ‖ x ‖ for x ∈ P ∩ ∂ Ω 2 .

Then F has a fixed point in the set P ∩ (Ω 2 \ Ω ) 1 .

In the sequel, for a given continuous function h : [0 , T ] → R we use the notation

$$h _ { * } = \min \{ h ( t ) \, \colon t \in [ 0, T ] \} \quad \text{and} \quad h ^ { * } = \max \{ h ( t ) \, \colon t \in [ 0, T ] \}.$$

## 3. Main results

Integrating on [0 , T ] the differential equation in the problem (3), and taking into account the boundary conditions, we get a necessary condition for the existence of positive solutions, namely

$$0 = \int _ { 0 } ^ { T } [ r ( t ) x ^ { \alpha } ( t ) - s ( t ) x ^ { \beta } ( t ) ] d t.$$

So it is easy to arrive to the following non-existence result.

Theorem 3.1. If one of the following conditions holds,

- (1) r ∗ ≥ 0 and s ∗ &lt; 0 ,
- (2) r ∗ &gt; 0 and s ∗ ≤ 0 ,
- (3) r ∗ ≤ 0 and s ∗ &gt; 0 ,
- (4) r ∗ &lt; 0 and s ∗ ≥ 0 ,

then problem (3) does not have positive solutions.

Our main result, concerning not only the existence but also the localization of positive solutions, is the following one.

Theorem 3.2. Assume that (H0) and the following condition hold:

$$( \text{H1} ) \text{ There exist } m > 0 \text{ and } 0 < R _ { 1 } < R _ { 2 } \text{ such that } m ^ { 2 } < \left ( \frac { \pi } { T } \right ) ^ { 2 } + \left ( \frac { a } { 2 } \right ) ^ { 2 },$$

$$( 7 ) \quad \ f _ { m } ( t, x ) = r ( t ) x ^ { \alpha } - s ( t ) x ^ { \beta } + m ^ { 2 } x \geq 0 \quad \text{for} \, t \in [ 0, T ] \text{ and } x \in [ c _ { m } R _ { 1 }, R _ { 2 } ],$$

$$f _ { m } ( t, x ) \geq m ^ { 2 } R _ { 1 } \ \text{ for } t \in [ 0, T ] \ a n d \ x \in [ c _ { m } R _ { 1 }, R _ { 1 } ]$$

and

$$f _ { m } ( t, x ) \leq m ^ { 2 } R _ { 2 } \ \ \text{ for } t \in [ 0, T ] \ \text{ and } x \in [ c _ { m } R _ { 2 }, R _ { 2 } ].$$

Then, the problem (3) has a positive solution x t ( ) such that

$$c _ { m } R _ { 1 } \leq x ( t ) \leq R _ { 2 }.$$

Proof. The problem (4) with the value of m given by (H1) is non-resonant and its Green's function G m ( t, s ) satisfies the properties i), ii) and iii) stated in the Introduction. Let us define

$$\Omega _ { j } = \{ x \in X \, \colon \| x \| < R _ { j } \} \ \ j = 1, 2.$$

Then by (7) we have f m ( t, x ( )) t ≥ 0 for x ∈ P ∩ (Ω 2 \ Ω ). Now, for every 1 x ∈ P ∩ (Ω 2 \ Ω ) 1 and for every t ∈ [0 , T ], we have that

$$\int _ { 0 } ^ { T } G _ { m } ( s, s ) f _ { m } ( s, x ( s ) ) d s \leq F x ( t ) = \int _ { 0 } ^ { T } G _ { m } ( t, s ) f _ { m } ( s, x ( s ) ) d s \leq \int _ { 0 } ^ { T } \frac { G _ { m } ( s, s ) } { c _ { m } } f _ { m } ( s, x ( s ) ) d s,$$

and therefore F P ( ∩ (Ω 2 \ Ω )) 1 ⊂ P . A standard argument shows that F is a completely continuous operator.

Now, we check that condition (i) in Theorem 2.1 is fulfilled.

Claim 1. ‖ Fx ‖ ≥ ‖ x ‖ for x ∈ P ∩ ∂ Ω 1 .

Let x ∈ P ∩ ∂ Ω , that is, 1 x ∈ P and ‖ x ‖ = R 1 . Then we have that c m R 1 ≤ x t ( ) ≤ R 1 for all t ∈ [0 , T ] and hence

$$\L ^ { \cdot } _ { \cdot } F x ( t ) \ & = \ \int _ { 0 } ^ { T } G _ { m } ( t, s ) f _ { m } ( s, x ( s ) ) d s \geq \int _ { 0 } ^ { T } G _ { m } ( t, s ) m ^ { 2 } R _ { 1 } d s \\ & = \ m ^ { 2 } R _ { 1 } \int _ { 0 } ^ { T } G _ { m } ( t, s ) d s = m ^ { 2 } R _ { 1 } \frac { 1 } { m ^ { 2 } } = R _ { 1 } = \| x \|. \\ \L _ { \cdot } \ \L _ { \cdot } \ \L _ { \cdot } \ \L _ { \cdot } \ \L _ { \cdot } \ \L _ { \cdot } \ \L _ { \cdot }$$

$$0$$

Claim 2. ‖ Fx ‖ ≤ ‖ x ‖ for x ∈ P ∩ ∂ Ω 2 .

Let x ∈ P ∩ ∂ Ω , that is, 2 x ∈ P and ‖ x ‖ = R 2 . Then c m R 2 ≤ x t ( ) ≤ R 2 for all t ∈ [0 , T ] and hence

$$F x ( t ) \ & = \ \int _ { 0 } ^ { T } G _ { m } ( t, s ) f _ { m } ( s, x ( s ) ) d s \leq \int _ { 0 } ^ { T } G _ { m } ( t, s ) m ^ { 2 } R _ { 2 } d s \\ & = \ m ^ { 2 } R _ { 2 } \int _ { 0 } ^ { T } G _ { m } ( t, s ) d s = m ^ { 2 } R _ { 2 } \frac { 1 } { m ^ { 2 } } = R _ { 2 } = \| x \|. \\ \cdots \cdots \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{array}$$

Thefore by Theorem 2.1 the existence of a solution of the problem (3) with the desired localization property immediately follows. glyph[square]

The following result is a consequence of Theorem 3.2.

Theorem 3.3. Assume that (H0) and the following conditions hold:

$$\text{(H2)} \ T h e r e \ e x i s t s \ m > 0 \ s u c h \ t h a t \ m ^ { 2 } < \left ( \frac { \pi } { T } \right ) ^ { 2 } + \left ( \frac { a } { 9 } \right ) ^ { 2 } \ a n d$$

$$\left ( \text{H2} \right ) \text{ Theme exists } m > 0 \text{ such that } m ^ { 2 } < \left ( \frac { n } { T } \right ) ^ { \kappa } + \left ( \frac { u } { 2 } \right ) ^ { \kappa } \text{ and }$$

$$f _ { m } ( t, x ) = r ( t ) x ^ { \alpha } - s ( t ) x ^ { \beta } + m ^ { 2 } x \geq 0 \ \ \text{ for } t \in [ 0, T ] \text{ and } x \geq 0.$$

(H3) r ∗ &gt; 0 and s ∗ &gt; 0 .

Then the problem (3) has a positive solution.

Proof. We are going to prove that assumption (H1) in Theorem 3.2 is satisfied. By (H0), (H3) and since c m ∈ (0 , 1) we can choose 0 &lt; R 1 &lt; R 2 such that

$$( 1 - c _ { m } ) m ^ { 2 } R _ { 1 } ^ { 1 - \alpha } + s ^ { * } R _ { 1 } ^ { \beta - \alpha } \leq r _ { * } c _ { m } ^ { \alpha }$$

and

$$r ^ { * } \leq s _ { * } c _ { m } ^ { \beta } R _ { 2 } ^ { \beta - \alpha }.$$

Note that the inequality (7) follows from (H2) so it only remains to check the inequalities (8) and (9).

$$\text{If } t \in [ 0, T ] \text{ and } x \in [ c _ { m } R _ { 1 }, R _ { 1 } ] \text{ then, taking into account (H3) and (10), we have that} \\ f _ { m } ( t, x ) \geq r _ { * } ( c _ { m } R _ { 1 } ) ^ { \alpha } - s ^ { * } R _ { 1 } ^ { \beta } + m ^ { 2 } c _ { m } R _ { 1 } \geq m ^ { 2 } R _ { 1 },$$

so the inequality (8) is fulfilled.

$$\text{If } t \in [ 0, T ] \text{ and } x \in [ c _ { m } R _ { 2 }, R _ { 2 } ] \text{ then, taking into account (H3) and (11), we have that } \\ f _ { m } ( t, x ) \leq r ^ { * } R _ { 2 } ^ { \alpha } - s _ { * } ( c _ { m } R _ { 2 } ) ^ { \beta } + m ^ { 2 } R _ { 2 } \leq m ^ { 2 } R _ { 2 },$$

so the inequality (9) is also fulfilled.

glyph[square]

Remark 3.4. Note that (H0) and (H2) imply that r ∗ &gt; 0, thus the first part of condition (H3) is redundant but we have included it for the sake of clarity. On the other hand, notice that condition (7) in Theorem 3.2 could be satisfied even if r t ( ) assumes negative values.

Next, we present an explicit condition sufficient in order to get the condition (H2) in Theorem 3.3.

Corollary 3.5. Assume that (H0) and (H3) hold and, moreover, that

$$s ^ { * } < \min \{ \left ( \frac { \pi } { T } \right ) ^ { 2 } + \left ( \frac { a } { 2 } \right ) ^ { 2 }, r _ { * } \}.$$

Then the problem (3) has a positive solution.

Proof. It is enough to show that condition (12) implies (H2). Indeed, by (12) we can choose m&gt; 0 such that

$$s ^ { * } < m ^ { 2 } < \left ( \frac { \pi } { T } \right ) ^ { 2 } + \left ( \frac { a } { 2 } \right ) ^ { 2 }.$$

Now, we are going to prove that for such m we have that f m ( t, x ) ≥ 0 for all t ∈ [0 , T ] and all x ≥ 0.

Step 1. We show that f m ( t, x ) ≥ 0 for all t ∈ [0 , T ] and all x ∈ [0 , 1] .

Since 0 &lt; α &lt; β &lt; 1 we have x β ≤ x α for 0 ≤ x ≤ 1. Moreover s ∗ &lt; r ∗ by (12) and then

$$f _ { m } ( t, x ) \ & = \ r ( t ) x ^ { \alpha } - s ( t ) x ^ { \beta } + m ^ { 2 } x \geq r _ { * } x ^ { \alpha } - s ^ { * } x ^ { \beta } + m ^ { 2 } x \\ & \geq \ r _ { * } x ^ { \beta } - s ^ { * } x ^ { \beta } + m ^ { 2 } x = ( r _ { * } - s ^ { * } ) x ^ { \beta } + m ^ { 2 } x > 0.$$

Step 2. We show that f m ( t, x ) ≥ 0 for all t ∈ [0 , T ] and all x ≥ 1 . Since 0 &lt; β &lt; 1 we have x β ≤ x for x ≥ 1 and therefore

$$f _ { m } ( t, x ) \ & = \ r ( t ) x ^ { \alpha } - s ( t ) x ^ { \beta } + m ^ { 2 } x \geq r _ { * } x ^ { \alpha } - s ^ { * } x ^ { \beta } + m ^ { 2 } x \\ & \geq \ r _ { * } x ^ { \alpha } - s ^ { * } x ^ { \beta } + m _ { \substack { \epsilon \\ \epsilon } } ^ { 2 } x ^ { \beta } = r _ { * } x ^ { \alpha } + ( m ^ { 2 } - s ^ { * } ) x ^ { \beta } \geq 0.$$

glyph[square]

In the particular case of problem (2) we recover [4, Theorem 1.8].

Corollary 3.6. Assume that a ≥ 0 , b &gt; 1 , c &gt; 0 and e ∗ &gt; . 0 Then problem (2) has a positive solution provided that

$$\frac { \left ( b + 1 \right ) c ^ { 2 } } { 4 \, e _ { * } } < \left ( \frac { \pi } { T } \right ) ^ { 2 } + \frac { a ^ { 2 } } { 4 } \,.$$

Proof. We recall that the problem (2) is of the form (3) with

$$\mu = \frac { 1 } { b + 1 }, \, r ( t ) = \frac { e ( t ) } { \mu }, \, s ( t ) = \frac { c } { \mu }, \, \alpha = 1 - 2 \, \mu, \ \beta = 1 - \mu.$$

From our assumptions it follows immediately that the conditions (H0) and (H3) of Theorem 3.3 hold. Thus it only remains to check the condition (H2). By (13) we can choose m &gt; 0 such that

$$\frac { \left ( b + 1 \right ) c ^ { 2 } } { 4 \, e _ { * } } < m ^ { 2 } < \left ( \frac { \pi } { T } \right ) ^ { 2 } + \frac { a ^ { 2 } } { 4 } \,.$$

Then for all t ∈ [0 , T ] and x ≥ 0 we have

$$f _ { m } ( t, x ) \ & = \ \frac { e ( t ) } { \mu } \, x ^ { 1 - 2 \mu } - \frac { c } { \mu } \, x ^ { 1 - \mu } + m ^ { 2 } x \geq \frac { e _ { * } } { \mu } \, x ^ { 1 - 2 \mu } - \frac { c } { \mu } \, x ^ { 1 - \mu } + m ^ { 2 } x \\ & = \ \frac { x ^ { 1 - 2 \mu } } { \mu } \left ( e _ { * } - c x ^ { \mu } + m ^ { 2 } \, \mu \, x ^ { 2 \mu } \right ) \geq 0, \\. \. \. \. \ \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \partial _ { \lambda } \quad \odot.$$

because e ∗ -cx µ + m µx 2 2 µ &gt; 0 (it is a parabola in the variable x µ which has by (14) negative discriminant c 2 -4 m µe 2 ∗ &lt; 0, so it has constant sign, and positive independent term e ∗ &gt; 0). glyph[square]

The following example shows that Theorem 3.2 is in fact more general than Corollary 3.6. Therefore, even in the case of the model problem (2), Theorem 3.2 is a true extension of [4, Theorem 1.8]. Moreover, notice that we also obtain information about the localization of the solution which was not provided in [4]. Some computations here were made with MAPLE.

Example 3.7. Consider the problem (2), that is,

$$\begin{cases} \ x ^ { \prime \prime } ( t ) + a \, x ^ { \prime } ( t ) = \frac { e ( t ) } { \mu } \, x ^ { 1 - 2 \, \mu } - \frac { c } { \mu } \, x ^ { 1 - \mu }, \quad t \in [ 0, T ], \\ \ x ( 0 ) = x ( T ), \quad x ^ { \prime } ( 0 ) = x ^ { \prime } ( T ), \end{cases}$$

with the parameter values a = 1 6, . e t ( ) ≡ 1 54, . µ = 0 01, . c = 1 49 and . T = 1.

Therefore we have

$$b = 9 9 \, ( \text{recall that } \mu \coloneqq \frac { 1 } { b + 1 } ), \, r ( t ) \coloneqq \frac { e ( t ) } { \mu } \equiv 1 5 4, \, s ( t ) \coloneqq \frac { c } { \mu } \equiv 1 4 9, \\ \alpha \coloneqq 1 - 2 \mu = 0. 9 8 \quad \text{and} _ { 7 } \quad \beta \coloneqq 1 - \mu = 0. 9 9.$$

Thus, condition (13) in Corollary 3.6 is not satisfied because

$$\frac { \left ( b + 1 \right ) c ^ { 2 } } { 4 \, e _ { * } } \approx 3 6. 0 4 0 6 > 1 0. 5 0 9 6 \approx \left ( \frac { \pi } { T } \right ) ^ { 2 } + \frac { a ^ { 2 } } { 4 }.$$

Notice that neither condition (12) in Corollary 3.5 is satisfied since

$$s ^ { * } = 1 4 9 > 1 0. 5 0 9 6 \approx \left ( \frac { \pi } { T } \right ) ^ { 2 } + \frac { a ^ { 2 } } { 4 }.$$

However, the condition (H0) in Theorem 3.2 is clearly satisfied and if we define m := 0 7, . R 1 := 27 and R 2 := 29 then (H1) is also fulfilled (notice that c m ≈ 0 9414 as it is given by . (20)).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0130v1-figure-8ddef821-721a-41e2-b71c-c516c4c82044.jpg)

Thus, Theorem 3.2 ensures the existence of a solution x t ( ) of this problem such that

$$2 5. 4 1 8 9 \approx c _ { m } R _ { 1 } \leq x ( t ) \leq R _ { 2 } = 2 9.$$

Observe that f m is not positive for all x ≥ 0.

Remark 3.8. From a careful reading of the proof of [4, Theorem 1.8] it follows that the existence of a positive periodic solution to problem (2) is ensured if condition (13) in Corollary 3.6 is replaced by the more general one:

- (H) There is δ 0 &gt; 0 such that f ( t, x ) &lt; 0 for all t ∈ [0 , T ] and all x ∈ (0 , δ 0 ) and λ x ∗ -f ( t, x ) &gt; 0 for all t ∈ [0 , T ] and all x &gt; 0, where λ ∗ = ( π T ) 2 + a 2 4 .

However, condition (H) is neither satisfied for the problem presented in Example 3.7 because in that case

$$\lambda ^ { * } x - f ( t, x ) \approx 1 0. 5 0 9 6 x - 1 4 9 x ^ { 0. 9 9 } + 1 5 4 x ^ { 0. 9 8 },$$

which is a sing-changing nonlinearity on the positive real line. Indeed, dividing by x 0 98 . and making the change y = x 0 01 . , it is clear that λ x ∗ -f ( t, x ) is negative on ( r , r 1 2 ) and positive on (0 , r 1 ) ∪ ( r , 2 ∞ ), being r 1 ≈ 103620 4135 and . r 2 ≈ 3 7845 . ∗ 10 111 .

$$4. \ A P P E N D I X \colon \text{THE GREEN'S FUNCTION}$$

4.1. The case 0 &lt; m &lt; a 2 . The Green's function related to problem (5) is given by

$$G _ { m } ( t, s ) = \frac { 1 } { \lambda _ { 2 } - \lambda _ { 1 } } \begin{cases} \frac { e ^ { \lambda _ { 2 } ( t - s ) } } { 1 - e ^ { \lambda _ { 2 } T } } - \frac { e ^ { \lambda _ { 1 } ( t - s ) } } { 1 - e ^ { \lambda _ { 1 } T } }, \quad & 0 \leq s \leq t \leq T, \\ \frac { e ^ { \lambda _ { 2 } ( t - s + T ) } } { 1 - e ^ { \lambda _ { 2 } T } } - \frac { e ^ { \lambda _ { 1 } ( t - s + T ) } } { 1 - e ^ { \lambda _ { 1 } T } }, \quad 0 \leq t \leq s \leq T, \end{cases} \\ \text{where $\lambda$} \text{ and $\lambda$} \text{ are the characteristic roots of the homogeneous equation}$$

where λ 1 and λ 2 are the characteristic roots of the homogeneous equation

$$x ^ { \prime \prime } ( t ) + a x ^ { \prime } ( t ) + m ^ { 2 } x ( t ) = 0,$$

that is,

$$\lambda _ { 1 } = \frac { - a - \sqrt { a ^ { 2 } - 4 m ^ { 2 } } } { 2 }, \ \lambda _ { 2 } = \frac { - a + \sqrt { a ^ { 2 } - 4 m ^ { 2 } } } { 2 }.$$

Note that λ 1 &lt; λ 2 &lt; 0.

Remark 4.1. G m ( s, s ) is constant for s ∈ [0 , T ],

$$G _ { m } ( s, s ) = \frac { 1 } { \lambda _ { 2 } - \lambda _ { 1 } } \left ( \frac { 1 } { 1 - e ^ { \lambda _ { 2 } T } } - \frac { 1 } { 1 - e ^ { \lambda _ { 1 } T } } \right ),$$

and

$$\int _ { 0 } ^ { T } G _ { m } ( t, s ) d s = \frac { 1 } { \lambda _ { 1 } \lambda _ { 2 } } = \frac { 1 } { m ^ { 2 } }.$$

The Green's function (15) satisfies

$$G _ { m } ( t, s ) \geq G _ { m } ( s, s ) > 0,$$

and

$$G _ { m } ( s, s ) \geq \underset { 9 } { c } _ { m } G _ { m } ( t, s ),$$

for all ( t, s ) ∈ [0 , T ] × [0 , T ], where

$$c _ { m } = \frac { G _ { m } ( s, s ) } { \max \{ G _ { m } ( t, s ) \, \colon t, s \in [ 0, T ] \} },$$

that is,

$$c _ { m } = - \frac { \lambda _ { 2 } } { \lambda _ { 2 } - \lambda _ { 1 } } \frac { e ^ { \lambda _ { 2 } T } - e ^ { \lambda _ { 1 } T } } { ( 1 - e ^ { \lambda _ { 2 } T } ) \left ( \frac { ( 1 - e ^ { \lambda _ { 2 } T } ) \lambda _ { 1 } } { ( 1 - e ^ { \lambda _ { 1 } T } ) \lambda _ { 2 } } \right ) ^ { \frac { \lambda _ { 1 } } { \lambda _ { 2 } - \lambda _ { 1 } } } }.$$

4.2. The case m = a 2 . The Green's function related to the problem (5) is given by

$$G _ { m } ( t, s ) & = \frac { 1 } { e ^ { m T } - 1 } \begin{cases} e ^ { m ( t - s ) } \left [ \frac { T e ^ { m T } } { e ^ { m T } - 1 } + s - t \right ], & 0 \leq s \leq t \leq T, \\ e ^ { m ( T - s + t ) } \left [ \frac { T } { e ^ { m T } - 1 } + s - t \right ], & 0 \leq t \leq s \leq T. \end{cases}$$

Then we have

$$\int _ { 0 } ^ { T } G _ { m } ( t, s ) d s = \frac { 1 } { m ^ { 2 } }, \\ G _ { m } ( t, s ) \geq G _ { m } ( s, s ) > 0,$$

and

$$G _ { m } ( s, s ) \geq c _ { m } \, G _ { m } ( t, s ),$$

for all ( t, s ) ∈ [0 , T ] × [0 , T ], where

$$c _ { m } = \frac { G _ { m } ( s, s ) } { \max \{ G _ { m } ( t, s ) \colon t, s \in [ 0, T ] \} },$$

that is,

$$c _ { m } = \kappa e ^ { 1 - \kappa },$$

with

$$\kappa = \frac { m T } { e ^ { m T } - 1 }.$$

4.3. The case m&gt; a 2 and m &lt; 2 ( π T ) 2 + ( a 2 ) 2 . We use the following notation:

$$& 2 & \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad \quad \quad$$

The Green's function related to the problem (5) is given by

(22)

$$G _ { m } ( t, s ) = \frac { 1 } { \delta D } \begin{cases} e ^ { \gamma ( T + t - s ) } \sin ( \delta ( T + s - t ) ) + e ^ { \gamma ( t - s ) } \sin ( \delta ( t - s ) ), \quad & 0 \leq s \leq t \leq T, \\ e ^ { \gamma ( T + t - s ) } \sin ( \delta ( T + t - s ) ) + e ^ { \gamma ( 2 T + t - s ) } \sin ( \delta ( s - t ) ), \quad 0 \leq t \leq s \leq T. \end{cases}$$

Then we have and

$$G _ { m } ( s, s ) \geq c _ { m } \, G _ { m } ( t, s ),$$

for all ( t, s ) ∈ [0 , T ] × [0 , T ], where

$$c _ { m } = \frac { G _ { m } ( s, s ) } { \max \{ G _ { m } ( t, s ) \, \colon t, s \in [ 0, T ] \} } = \frac { \frac { 1 } { \delta D } e ^ { \gamma T } \sin ( \delta T ) } { \max \{ G _ { m } ( t, s ) \, \colon t, s \in [ 0, T ] \} }.$$

If a = 0 then γ = 0 , δ = m, D = 2(1 -cos mT ) and it is known that max { G m ( t, s ) : t, s ∈ [0 , T ] } = 1 2 m sin mT 2 , see [16]. Therefore we have c m = cos mT 2 .

On the other hand, if a &gt; 0 it is difficult to find exactly the maximum of G m ( t, s ), but we are able to get an upper bound in the following way: since 0 &lt; δT &lt; π we have for 0 ≤ s ≤ t ≤ T (the case 0 ≤ t ≤ s ≤ T is analogous)

$$e ^ { \gamma ( T + t - s ) } \sin ( \delta ( T + s - t ) ) & + e ^ { \gamma ( t - s ) } \sin ( \delta ( t - s ) ) \leq \sin ( \delta ( T + s - t ) ) + \sin ( \delta ( t - s ) ) \\ & = 2 \sin \frac { \delta T } { 2 } \cos \left ( \delta ( \frac { T } { 2 } - ( t - s ) ) \right ) \leq 2 \sin \frac { \delta T } { 2 }.$$

So, max { G m ( t, s ) : t, s ∈ [0 , T ] } ≤ 2 δD sin δT 2 which yields c m ≥ e γT cos δT 2 .

$$\text{Remark } 4. 2. \text{ If } m ^ { 2 } = \left ( \frac { \pi } { T } \right ) ^ { 2 } + \left ( \frac { a } { 2 } \right ) ^ { 2 } \text{ then } G _ { m } ( s, s ) = 0.$$

## Acknowledgements

J. A. Cid was partially supported by Ministerio de Educaci´n y Ciencia, o Spain, and FEDER, Project MTM2010-15314, G. Infante was partially supported by G.N.A.M.P.A. -INdAM (Italy), M. Tvrd´ was supported by GA CR Grant P201/14-06958S and RVO: y ˇ 67985840 and M. Zima was partially supported by the Centre for Innovation and Transfer of Natural Science and Engineering Knowledge of University of Rzesz´w. o

## References

[1] A. Borz` ı and G. Propst, Numerical investigation of the Liebau phenomenon, Z. Angew. Math. Phys. 54 (2003) 1050-1072.

[2] A. Cabada, Green's Functions in the Theory of Ordinary Differential Equations, SpringerBriefs in Mathematics, 2014.

[3] A. Cabada, J. A. Cid and B. M´quez-Villamar´ ın, Computation of Green's functions for boundary value a problems with Mathematica, Appl. Math. Comput. 219 (2012), 1919-1936.

[4] J. A. Cid, G. Propst and M. Tvrd´, On the pumping effect in a pipe/tank flow configuration with y friction, Phys. D 273-274 (2014) 28-33.

$$\int _ { 0 } ^ { T } G _ { m } ( t, s ) d s = \frac { 1 } { m ^ { 2 } }, \\ G _ { m } ( t, s ) \geq G _ { m } ( s, s ) > 0,$$

$$\dots$$

| [5]                                                                                                                                | C. de Coster and P. Habets, Two-point Boundary Value Problems. Lower and Upper Solutions , Elsevier, Amsterdam (2006).                                                                                                       |
|------------------------------------------------------------------------------------------------------------------------------------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [6]                                                                                                                                | D. Guo and V. Lakshmikantham, Nonlinear problems in abstract cones , Academic Press, Boston, 1988.                                                                                                                           |
| [7]                                                                                                                                | X. Han, Positive solutions for a three-point boundary value problems at resonance, J. Math. Anal. Appl. 336 (2007), 556-568.                                                                                                 |
| [8]                                                                                                                                | G. Infante, P. Pietramala and M. Tenuta, Existence and localization of positive solutions for a nonlocal BVP arising in chemical reactor theory, Commun. Nonlinear Sci. Numer. Simulat. , 19 (2014), 2245- 2251.             |
| [9]                                                                                                                                | G. Infante, P. Pietramala and F. A. F. Tojo, Nontrivial solutions of local and nonlocal Neumann boundary value problems, arXiv:1404.1390 (2014).                                                                             |
| [10]                                                                                                                               | M. A. Krasnosel'skii ˘ and P. P. Zabreiko, ˘ Geometrical methods of nonlinear analysis , Springer-Verlag, Berlin, (1984).                                                                                                    |
| [11]                                                                                                                               | K. Q. Lan, Multiple positive solutions of Hammerstein integral equations and applications to periodic boundary value problems, Appl. Math. Comput. , 154 (2004), 531-542.                                                    |
| [12]                                                                                                                               | G. Liebau, ¨ Uber ein ventilloses Pumpprinzip, Naturwissenschaften 41 (1954), 327.                                                                                                                                           |
| [13]                                                                                                                               | M. Moser, J. W. Huang, G. S. Schwarz, T. Kenner and A. Noordergraaf, Impedance defined flow. Generalisation of William Harvey's concept of the circulation-370 years later, Int. J. Cardiovasc. Med. Sci. 1 (1998), 205-211. |
| [14]                                                                                                                               | G. Propst, Pumping effects in models of periodically forced flow configurations, Phys. D 217 (2006), 193-201.                                                                                                                |
| [15]                                                                                                                               | I. Rach˚nkov´, u a M Tvrd´ y and I. Vrkoˇ, c Existence of nonnegative and nonpositive solutions for second order periodic boundary value problems, J. Differential Equations 176 (2001), 445-469.                            |
| [16]                                                                                                                               | P. J. Torres, Existence of one-signed periodic solutions of some second-order differential equations via a Krasnoselskii fixed point theorem, J. Differential Equations , 190 (2003), 643-662.                               |
| [17]                                                                                                                               | J. R. L. Webb and M. Zima, Multiple positive solutions of resonant and non-resonant nonlocal boundary value problems, Nonlinear Anal. 71 (2009), 1369-1378.                                                                  |
| Jos´ e ´ Angel Cid, Departamento de Matem´ticas, a Universidade de Vigo, 32004, Pabell´n o 3, Campus de Ourense, Spain             | Jos´ e ´ Angel Cid, Departamento de Matem´ticas, a Universidade de Vigo, 32004, Pabell´n o 3, Campus de Ourense, Spain                                                                                                       |
| Gennaro Infante, Dipartimento di Matematica e Informatica, Universit` a della Calabria, 87036 Arcavacata di Rende, Cosenza, Italy  | Gennaro Infante, Dipartimento di Matematica e Informatica, Universit` a della Calabria, 87036 Arcavacata di Rende, Cosenza, Italy                                                                                            |
| Milan Tvrd´, y Mathematical Institute, Academy of Sciences of the Czech Republic, CZ 115                                           | Milan Tvrd´, y Mathematical Institute, Academy of Sciences of the Czech Republic, CZ 115                                                                                                                                     |
| E-mail address : tvrdy@math.cas.cz Mirosglyph[suppress] lawa Zima, Department of Differential Equations and Statistics, Faculty of | E-mail address : tvrdy@math.cas.cz Mirosglyph[suppress] lawa Zima, Department of Differential Equations and Statistics, Faculty of                                                                                           |
| and Natural Sciences, University of Rzesz´w, o Pigonia 1, 35-959 Rzesz´w, o Poland.                                                | and Natural Sciences, University of Rzesz´w, o Pigonia 1, 35-959 Rzesz´w, o Poland.                                                                                                                                          |
| E-mail address : mzima@ur.edu.pl                                                                                                   | E-mail address : mzima@ur.edu.pl                                                                                                                                                                                             |
| ematics                                                                                                                            | ematics                                                                                                                                                                                                                      |
| Praha 1, a 25, Czech Republic                                                                                                      | Praha 1, a 25, Czech Republic                                                                                                                                                                                                |
| Zitn´                                                                                                                              | Zitn´                                                                                                                                                                                                                        |
| ˇ                                                                                                                                  | ˇ                                                                                                                                                                                                                            |
| E-mail address : gennaro.infante@unical.it                                                                                         | E-mail address : gennaro.infante@unical.it                                                                                                                                                                                   |
|                                                                                                                                    | Math-                                                                                                                                                                                                                        |
José Ángel Cid Gennaro Infante Milan Tvrdý Mirosława Zima
A topological approach to periodic oscillations related to the Liebau phenomenon
We give some sufficient conditions for existence, non-existence and
localization of positive solutions for a periodic boundary value problem
related to the Liebau phenomenon. Our approach is of topological nature and
relies on the Krasnosel'ski\u\i{}-Guo theorem on cone expansion and
compression. Our results improve and complement earlier ones in the literature.