Next Article in Journal
Faces of 2-Dimensional Simplex of Order and Chain Polytopes
Next Article in Special Issue
Variational Methods for an Impulsive Fractional Differential Equations with Derivative Term
Previous Article in Journal
Set-Valued Interpolative Hardy–Rogers and Set-Valued Reich–Rus–Ćirić-Type Contractions in b-Metric Spaces
Previous Article in Special Issue
Numerical Solution of the Boundary Value Problems Arising in Magnetic Fields and Cylindrical Shells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Lyapunov-Type Inequality for a Laplacian System on a Rectangular Domain with Zero Dirichlet Boundary Conditions

Department of Mathematics, College of Science King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(9), 850; https://doi.org/10.3390/math7090850
Submission received: 1 July 2019 / Revised: 10 September 2019 / Accepted: 12 September 2019 / Published: 14 September 2019
(This article belongs to the Special Issue Mathematical Analysis and Boundary Value Problems)

Abstract

:
We consider a coupled system of partial differential equations involving Laplacian operator, on a rectangular domain with zero Dirichlet boundary conditions. A Lyapunov-type inequality related to this problem is derived. This inequality provides a necessary condition for the existence of nontrivial positive solutions.

1. Introduction

Consider the second order differential equation
u ( x ) = p ( x ) u ( x ) , a < x < b ,
subject to the Dirichlet boundary condition
u ( a ) = u ( b ) = 0 ,
where a , b R , a < b , and p C ( [ a , b ] ) . It is well known (see, e.g., [1,2]) that if u C 2 ( [ a , b ] ) is a nontrivial solution to Equations (1) and (2), then
a b | p ( x ) | d x > 4 b a .
The inequality in Equation (3) is known in the literature as Lyapunov inequality. It has many applications in the study of spectral properties of ODE (see, e.g., [3,4,5,6,7,8]).
In the multi-dimensional case, there are some important works dealing with Lyapunov-type inequalities for PDEs. In [9], the authors studied the Laplace equation
Δ u = q ( x ) u , x Ω ,
under Neumann boundary conditions, where Ω is an open bounded domain in R N ( N 2 ) and q : Ω R . In [10], the authors studied the p-Laplacian equation
Δ p u = q ( x ) | u | p 2 u , x Ω ,
under Dirichlet boundary conditions, where Ω is an open bounded domain in R N ( N 2 ), p > 1 and q L s ( Ω ) , for some s, which depends on p and N. In [11], the authors studied the fractional p–Laplacian equation
( Δ p ) s u = q ( x ) | u | p 2 u , x Ω ,
under the boundary conditions
u ( x ) = 0 , x R N \ Ω ,
where Ω is an open bounded domain in R N ( N 2 ), p > 1 , 0 < s < 1 and q L ( Ω ) . In [12], the authors extended the obtained results in [11] to a fractional p-Laplacian system. In [13], the authors studied the partial differential equation
G γ u ( x , y ) = q ( x ) u ( x , y ) , ( x , y ) Ω ,
under Dirichlet boundary conditions, where Ω = ] a , b [ × O , ( a , b ) R 2 , a < b , O is an open bounded subset in R N ( N 1 ), q C ( [ a , b ] ) and G γ , γ 0 , is the differential operator given by
G γ u ( x , y ) = 2 u x 2 + x 2 γ Δ y u ( x , y ) , ( x , y ) Ω .
The authors of [14] extended the obtained results in [13] to the differential operator
L γ , g u ( x , y ) = G γ u ( x , y ) + g ( x ) u x ( x , y ) , ( x , y ) Ω ,
where g C ( [ a , b ] ) and Δ y is the Laplacian operator with respect to the variable “y”.
Motivated by the above cited works, in this paper, we consider the two-dimensional coupled system of partial differential equations
Δ u ( x , y ) = p 11 ( x ) u ( x , y ) + p 12 ( x ) Γ ( α ) a x ( x z ) α 1 v ( z , y ) d z , ( x , y ) Ω , Δ v ( x , y ) = p 21 ( x ) Γ ( β ) a x ( x z ) β 1 u ( z , y ) d z + p 22 ( x ) v ( x , y ) , ( x , y ) Ω ,
under the Dirichlet boundary conditions
u ( x , y ) = v ( x , y ) = 0 , ( x , y ) Ω ,
where Ω = ] a , b [ × ] c , d [ , ( a , b ) , ( c , d ) R 2 , a < b , c < d , α > 0 , β > 0 , p i j C ( [ a , b ] ) , 1 i , j 2 , and Γ is the Gamma function. We derive a Lyapunov-type inequality, which provides a necessary condition for the existence of nontrivial positive solutions to Equations (4) and (5). Note that the used technique in this paper is different to that used in [11,12]. The approach used in this paper is based on an eigenvalue method from Kaplan [15]. Observe that the system in Equation (4) involves the nonlocal operators
1 Γ ( κ ) a x ( x z ) κ 1 f ( z ) d z ,
where κ { α , β } and f { u ( · , y ) , v ( · , y ) } . Such operators are known in the literature as Riemann–Liouville fractional integrals of order κ . For more details on fractional operators and their applications, see, for example, [16,17].
The rest of the paper is organized as follows. In Section 2, we recall and prove some results on matrices theory that are used in the proof of our main result. In Section 3, we establish a Lyapunov-type inequality for Equations (4) and (5), and we discuss some special cases of Equations (4) and (5).

2. Preliminaries

Let N 1 be a natural number. We denote by 0 N the zero vector in R N . Let · N be the Euclidean norm in R N , that is,
V N = V 1 2 + V 2 2 + + V N 2 , V = V 1 V 2 V N R N .
We define in R N the partial order N given by
U = U 1 U 2 U N N V = V 1 V 2 V N U i V i , i = 1 , 2 , , N .
It can be easily seen that
Lemma 1.
Let U and V be two vectors in R N such that
0 N N U N V .
then
U N V N .
We denote by M N ( R ) the set of square matrices of size N with coefficients in R , and by M N ( R + ) the subset of M N ( R ) with positive coefficients. We endow M N ( R ) with the subordinate matrix norm
M M N = sup X R N , X 0 N M X N X N , M M N ( R ) .
For a given matrix M M N ( R ) , let ρ ( M ) be its spectral radius, i.e.,
ρ ( M ) = max { | λ i ( M ) | : i = 1 , 2 , , N } ,
where λ i ( M ) , i = 1 , 2 , , N , are the (real or complex) eigenvalues of M.
The following result is standard in the theory of matrices.
Lemma 2.
Let M M N ( R ) . Then
ρ ( M ) < 1 lim n M n M N = 0 .
Lemma 3.
Let M M N ( R + ) and X R N , X 0 N . If
0 N N X N M X ,
then
ρ ( M ) 1 .
Proof. 
From Equation (6) and using the fact that M M N ( R + ) , for all natural number n 1 , we have
0 N N X N M n X .
Next, by Lemma 1, we obtain
X N M n X N M n M N X N , n 1 .
Since X 0 N , we get
M n M N 1 , n 1 .
Finally, using Lemma 2, Equation (7) follows. ☐
Lemma 4.
Let M = ( m i j ) 1 i , j 2 M 2 ( R + ) . Then,
ρ ( M ) = m 11 + m 22 + ( m 11 m 22 ) 2 + 4 m 21 m 12 2 .
Proof. 
Let P M be the characteristic polynomial of the matrix M, that is,
P M ( λ ) = λ 2 t r ( M ) λ + det ( M ) , λ C ,
where t r ( M ) is the trace of M and det ( M ) is its determinant. Then, the discriminant of P M is given by
Δ ( P M ) = [ t r ( M ) ] 2 4 det ( M ) ,
i.e.,
Δ ( P M ) = m 11 m 22 2 + 4 m 21 m 12 .
Hence, M admits two eigenvalues
λ 1 ( M ) = m 11 + m 22 + m 11 m 22 2 + 4 m 21 m 12 2
and
λ 2 ( M ) = m 11 + m 22 m 11 m 22 2 + 4 m 21 m 12 2 .
Observe that
λ 1 ( M ) λ 2 ( M ) .
Therefore, if λ 2 ( M ) 0 , then ρ ( M ) = λ 1 ( M ) . Further, if λ 2 ( M ) < 0 , we obtain
| λ 2 ( M ) | = λ 2 ( M ) = m 11 m 22 2 + 4 m 21 m 12 m 11 + m 22 2 λ 1 ( M ) .
Then, in all cases, we have ρ ( M ) = λ 1 ( M ) , which proves the desired result. ☐

3. Lyapunov-Type Inequalities

In this section, a Lyapunov-type inequality is derived for Equations (4) and (5) and some special cases are discussed.
Definition 1.
We say that ( u , v ) is a nontrivial positive solution to Equations (4) and (5), if and only if,
(i) 
( u , v ) C 2 ( Ω ¯ ) × C 2 ( Ω ¯ ) satisfies Equations (4) and (5).
(ii) 
For all ( x , y ) Ω , we have
u ( x , y ) v ( x , y ) 2 0 2 .
(iii) 
u 0 and v 0 .

3.1. From PDEs to ODEs

In this subsection, we reduce the study of Equation (4) to a coupled system of ordinary differential equations.
Suppose that ( u , v ) is a nontrivial positive solution to Equations (4) and (5). Let us introduce the functions
φ ( y ) = sin π ( y c ) d c , c y d ,
u ¯ ( x ) = c d u ( x , y ) φ ( y ) d y , a x b
and
v ¯ ( x ) = c d v ( x , y ) φ ( y ) d y , a x b .
Observe that due to the positivity of the function φ in [ c , d ] , (ii) and (iii), we have u ¯ 0 and v ¯ 0 . Further, multiplying the first equation in Equation (4) by φ ( y ) and integrating over ( c , d ) , we obtain
d 2 d x 2 c d u ( x , y ) φ ( y ) d y c d 2 u y 2 ( x , y ) φ ( y ) d y = p 11 ( x ) c d u ( x , y ) φ ( y ) d y + p 12 ( x ) Γ ( α ) c d a x ( x z ) α 1 v ( z , y ) d z φ ( y ) d y ,
for all a < x < b . On the other hand, using an integration by parts, we obtain
c d 2 u y 2 ( x , y ) φ ( y ) d y = u y ( x , y ) φ ( y ) y = c d π d c c d u y ( x , y ) cos π ( y c ) d c d y = π d c c d u y ( x , y ) cos π ( y c ) d c d y .
Again, using an integration by parts and Equation (5), we obtain
c d 2 u y 2 ( x , y ) φ ( y ) d y = π d c u ( x , y ) cos π ( y c ) d c y = c d π d c 2 c d u ( x , y ) φ ( y ) d y = π d c 2 c d u ( x , y ) φ ( y ) d y .
Next, using Fubini’s theorem, we have
c d a x ( x z ) α 1 v ( z , y ) d z φ ( y ) d y = a x ( x z ) α 1 c d v ( z , y ) φ ( y ) d y d z .
Combining Equations (10), (11) and (12), we obtain
u ¯ ( x ) = p 11 ( x ) π d c 2 u ¯ ( x ) + p 12 ( x ) Γ ( α ) a x ( x z ) α 1 v ¯ ( z ) d z ,
for a < x < b . Similarly, multiplying the second equation in Equation (4) by φ ( y ) and integrating over ( c , d ) , we obtain
v ¯ ( x ) = p 21 ( x ) Γ ( β ) a x ( x z ) β 1 u ¯ ( z ) d z + p 22 ( x ) π d c 2 v ¯ ( x ) ,
for a < x < b . Moreover, using the boundary conditions in Equation (5), we have
u ¯ ( a ) = u ¯ ( b ) = 0 = v ¯ ( a ) = v ¯ ( b ) .
Hence, we have the following result.
Proposition 1.
Let ( u , v ) be a nontrivial positive solution to Equations (4) and (5). Then, ( u ¯ , v ¯ ) is a nontrivial solution to the coupled system
u ¯ ( x ) = p 11 ( x ) π d c 2 u ¯ ( x ) + p 12 ( x ) Γ ( α ) a x ( x z ) α 1 v ¯ ( z ) d z , a < x < b , v ¯ ( x ) = p 21 ( x ) Γ ( β ) a x ( x z ) β 1 u ¯ ( z ) d z + p 22 ( x ) π d c 2 v ¯ ( x ) , a < x < b ,
subject to the boundary conditions
u ¯ ( a ) = u ¯ ( b ) = 0 = v ¯ ( a ) = v ¯ ( b ) ,
where u ¯ and v ¯ are given by Equations (8) and (9).

3.2. Main Result

Let us introduce the matrix M = ( m i j ) 1 i , j 2 M 2 ( R + ) given by
m i j = ( b a ) γ i j + 1 4 Γ ( γ i j + 1 ) a b p i j ( s ) π d c 2 δ i j d s , 1 i , j 2 ,
where
δ i j = 1 if i = j , 0 if i j
and
γ i j = 0 if i = j , α if ( i , j ) = ( 1 , 2 ) , β if ( i , j ) = ( 2 , 1 ) .
Now, we are able to state and prove our main result.
Theorem 1.
Let ( u , v ) be a nontrivial positive solution to Equations (4) and (5). Then,
m 11 + m 22 + ( m 11 m 22 ) 2 + 4 m 21 m 12 2 .
Proof. 
Let ( u , v ) be a nontrivial positive solution to Equations (4) and (5). By Proposition 1, we have
u ¯ ( x ) = f ( x , u ¯ , v ¯ ) , a < x < b , u ¯ ( a ) = u ¯ ( b ) = 0 ,
where u ¯ 0 and
f ( x , u ¯ , v ¯ ) = p 11 ( x ) π d c 2 u ¯ ( x ) + p 12 ( x ) Γ ( α ) a x ( x z ) α 1 v ¯ ( z ) d z , a < x < b .
Therefore, u ¯ is a solution to the integral equation
u ¯ ( x ) = a b G ( x , s ) f ( s , u ¯ , v ¯ ) d s , a x b ,
where
G ( x , s ) = ( x a ) ( b s ) b a , ( x , s ) [ a , b ] × [ a , b ] .
On the other hand, the arithmetic-geometric-harmonic mean inequality yields
max a x , s b | G ( x , s ) | b a 4 .
Further, let us estimate the term | f ( s , u ¯ , v ¯ ) | , a < s < b . We have
p 12 ( s ) Γ ( α ) a s ( s z ) α 1 v ¯ ( z ) d z | p 12 ( s ) | Γ ( α ) a s ( s z ) α 1 d s v ¯ = | p 12 ( s ) | Γ ( α + 1 ) ( s a ) α v ¯ ( b a ) α Γ ( α + 1 ) | p 12 ( s ) | v ¯ ,
where
h = max a t b | h ( t ) | , h C ( [ a , b ] ) .
Moreover, we have
p 11 ( s ) π d c 2 u ¯ ( s ) p 11 ( s ) π d c 2 u ¯ .
Next, combining Equations (18) and (19), we get
| f ( s , u ¯ , v ¯ ) | p 11 ( s ) π d c 2 u ¯ + ( b a ) α Γ ( α + 1 ) | p 12 ( s ) | v ¯ , a < s < b .
Now, combining Equations (16), (17) and (20), we deduce that
u ¯ m 11 u ¯ + m 12 v ¯ .
Similarly, by Proposition 1, we have
v ¯ ( x ) = g ( x , u ¯ , v ¯ ) , a < x < b , v ¯ ( a ) = v ¯ ( b ) = 0 ,
where v ¯ 0 and
g ( x , u ¯ , v ¯ ) = p 21 ( x ) Γ ( β ) a x ( x z ) β 1 u ¯ ( z ) d z + p 22 ( x ) π d c 2 v ¯ ( x ) , a < x < b .
Using a similar argument as above, we obtain
v ¯ m 21 u ¯ + m 22 v ¯ .
Using Equations (21) and (22), we deduce that
0 2 2 X 2 M X ,
where
X = u ¯ v ¯ 0 2 .
Hence, by Lemma 3, we deduce that
ρ ( M ) 1 .
Finally, using Lemma 4, Equation (15) follows. ☐

3.3. Particular Cases

In this subsection, we discuss some special cases following from Theorem 1.

3.3.1. The Case α = β

Let us consider the coupled system
Δ u ( x , y ) = p 11 ( x ) u ( x , y ) + p 12 ( x ) Γ ( θ ) a x ( x z ) θ 1 v ( z , y ) d z , ( x , y ) Ω , Δ v ( x , y ) = p 21 ( x ) Γ ( θ ) a x ( x z ) θ 1 u ( z , y ) d z + p 22 ( x ) v ( x , y ) , ( x , y ) Ω ,
where θ > 0 and p i j C ( [ a , b ] ) , 1 i , j 2 . Observe that Equation (23) is a special case of Equation (4) with α = β = θ . Using Theorem 1, we obtain the following Lyapunov-type inequality for Equations (23) and (5).
Corollary 1.
Let ( u , v ) be a nontrivial positive solution to Equations (23) and (5). Then,
a b i = 1 2 p i i ( s ) π d c 2 d s + a b i = 1 2 ( 1 ) i + 1 p i i ( s ) π d c 2 d s 2 + 4 ( b a ) 2 θ Γ ( θ + 1 ) a b | p 12 ( s ) | d s a b | p 21 ( s ) | d s 8 b a .

3.3.2. The Limit Case θ 0 +

In the limit case θ 0 + , Equation (23) reduces to
Δ u ( x , y ) = p 11 ( x ) u ( x , y ) + p 12 ( x ) v ( x , y ) , ( x , y ) Ω , Δ v ( x , y ) = p 21 ( x ) u ( x , y ) + p 22 ( x ) v ( x , y ) , ( x , y ) Ω .
Therefore, passing to the limit as θ 0 + in Equation (24), we deduce the following Lyapunov-type inequality for Equations (25) and (5).
Corollary 2.
Let ( u , v ) be a nontrivial positive solution to Equations (25) and (5). Then,
a b i = 1 2 p i i ( s ) π d c 2 d s + a b i = 1 2 ( 1 ) i + 1 p i i ( s ) π d c 2 d s 2 + 4 a b | p 12 ( s ) | d s a b | p 21 ( s ) | d s 8 b a .

3.3.3. The Case p 11 = p 22

Let us consider the system
Δ u ( x , y ) = p ( x ) u ( x , y ) + p 12 ( x ) v ( x , y ) , ( x , y ) Ω , Δ v ( x , y ) = p 21 ( x ) u ( x , y ) + p ( x ) v ( x , y ) , ( x , y ) Ω ,
where p 12 , p 21 , p C ( [ a , b ] ) . Observe that Equation (27) is a special case of Equation (25) with
p 11 = p 22 = p .
Therefore, taking in Equation (26) p 11 and p 22 as above, we deduce the following Lyapunov-type inequality for Equation (27) and (5).
Corollary 3.
Let ( u , v ) be a nontrivial positive solution to Equations (27) and (5). Then,
a b p ( s ) π d c 2 d s + a b | p 12 ( s ) | d s a b | p 21 ( s ) | d s 4 b a .

3.3.4. The Case of a Single Equation

Taking in Equation (27) p 12 = p 21 0 and u = v , we obtain the single equation
Δ u ( x , y ) = p ( x ) u ( x , y ) , ( x , y ) Ω .
Hence, by Equation (28), we deduce the following Lyapunov-type inequality (see [13]) for Equation (29) under the Dirichlet boundary conditions
u ( x , y ) = 0 , ( x , y ) Ω .
Corollary 4.
Let u be a nontrivial positive solution to Equations (29) and (30). Then,
a b p ( s ) π d c 2 d s 4 b a .
Remark 1.
In this paper, the system in Equation (4) is discussed under Dirichlet boundary conditions. It would be interesting to study other types of boundary conditions, as well as more general domains.

Author Contributions

M.J. contributed to Section 1 and Section 2. B.S. contributed to Section 3.

Funding

This research was funded by Deanship of Scientific Research at King Saud University (RGP-1435-034).

Acknowledgments

The authors would like to extend their sincere appreciation to the Deanship of Scientific Research at King Saud University for its funding of this research through the Research Group Project No RGP-1435-034.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Borg, G. On a Liapounoff criterion of stability. Am. J. Math. 1949, 71, 67–70. [Google Scholar] [CrossRef]
  2. Lyapunov, A. Problème Général de la Stabilité du Mouvement. Ann. Fac. Sci. Toulouse 1907, 9, 204–474. [Google Scholar]
  3. Das, K.M.; Vatsala, A.S. Green’s function for n-n boundary value problem and an analogue of Hartman’s result. J. Math. Anal. Appl. 1975, 51, 670–677. [Google Scholar] [CrossRef]
  4. Elbert, A. A half-linear second order differential equation. Colloq. Math. Soc. János Bolyai 1979, 30, 158–180. [Google Scholar] [CrossRef]
  5. Hartman, P.; Wintner, A. On an oscillation criterion of Liapunoff. Am. J. Math. 1951, 73, 885–890. [Google Scholar] [CrossRef]
  6. De Nápoli, P.L.; Pinasco, J.P. Estimates for eigenvalues of quasilinear elliptic systems. J. Differ. Equations 2006, 227, 102–115. [Google Scholar] [CrossRef] [Green Version]
  7. Nehari, Z. On the zeros of solutions of second-order linear differential equations. Am. J. Math. 1954, 76, 689–697. [Google Scholar] [CrossRef]
  8. Wintner, A. On the non-existence of conjugate points. Am. J. Math. 1951, 73, 368–380. [Google Scholar] [CrossRef]
  9. Cañada, A.; Montero, J.A.; Villegas, S. Lyapunov inequalities for partial differential equations. J. Funct. Anal. 2006, 237, 176–193. [Google Scholar] [CrossRef] [Green Version]
  10. De Nápoli, P.L.; Pinasco, J.P. Lyapunov-type inequalities for partial differential equations. J. Funct. Anal. 2016, 270, 1995–2018. [Google Scholar] [CrossRef]
  11. Jleli, M.; Kirane, M.; Samet, B. Lyapunov-type inequalities for fractional partial differential equations. Appl. Math. Lett. 2017, 66, 30–39. [Google Scholar] [CrossRef]
  12. Jleli, M.; Kirane, M.; Samet, B. Lyapunov-type inequalities for a fractional p-Laplacian system. Fract. Calc. Appl. Anal. 2017, 20, 1485–1506. [Google Scholar] [CrossRef]
  13. Jleli, M.; Kirane, M.; Samet, B. On Lyapunov-type inequalities for a certain class of partial differential equations. Appl. Anal. 2018. [Google Scholar] [CrossRef]
  14. Agarwal, R.P.; Jleli, M.; Samet, B. On De La Vallée Poussin-type inequalities in higher dimension and applications. Appl. Math. Lett. 2018, 86, 264–269. [Google Scholar] [CrossRef]
  15. Kaplan, S. On the growth of solutions of quasilinear parabolic equations. Comm. Pure Appl. Math. 1963, 16, 305–333. [Google Scholar] [CrossRef]
  16. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier Science Inc.: New York, NY, USA, 2006. [Google Scholar]
  17. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: Longhorne, PA, USA, 1993. [Google Scholar]

Share and Cite

MDPI and ACS Style

Jleli, M.; Samet, B. A Lyapunov-Type Inequality for a Laplacian System on a Rectangular Domain with Zero Dirichlet Boundary Conditions. Mathematics 2019, 7, 850. https://doi.org/10.3390/math7090850

AMA Style

Jleli M, Samet B. A Lyapunov-Type Inequality for a Laplacian System on a Rectangular Domain with Zero Dirichlet Boundary Conditions. Mathematics. 2019; 7(9):850. https://doi.org/10.3390/math7090850

Chicago/Turabian Style

Jleli, Mohamed, and Bessem Samet. 2019. "A Lyapunov-Type Inequality for a Laplacian System on a Rectangular Domain with Zero Dirichlet Boundary Conditions" Mathematics 7, no. 9: 850. https://doi.org/10.3390/math7090850

APA Style

Jleli, M., & Samet, B. (2019). A Lyapunov-Type Inequality for a Laplacian System on a Rectangular Domain with Zero Dirichlet Boundary Conditions. Mathematics, 7(9), 850. https://doi.org/10.3390/math7090850

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop