Citation: Murali Ramdoss, Ponmana Selvan-Arumugam, Choonkil Park. Ulam stability of linear differential equations using Fourier transform[J]. AIMS Mathematics, 2020, 5(2): 766-780. doi: 10.3934/math.2020052
[1] | Hui Sun, Zhongyang Sun, Ya Huang . Equilibrium investment and risk control for an insurer with non-Markovian regime-switching and no-shorting constraints. AIMS Mathematics, 2020, 5(6): 6996-7013. doi: 10.3934/math.2020449 |
[2] | Gideon Simpson, Daniel Watkins . Relative entropy minimization over Hilbert spaces via Robbins-Monro. AIMS Mathematics, 2019, 4(3): 359-383. doi: 10.3934/math.2019.3.359 |
[3] | Abdullah Ali H. Ahmadini, Amal S. Hassan, Ahmed N. Zaky, Shokrya S. Alshqaq . Bayesian inference of dynamic cumulative residual entropy from Pareto Ⅱ distribution with application to COVID-19. AIMS Mathematics, 2021, 6(3): 2196-2216. doi: 10.3934/math.2021133 |
[4] | I. A. Husseiny, M. Nagy, A. H. Mansi, M. A. Alawady . Some Tsallis entropy measures in concomitants of generalized order statistics under iterated FGM bivariate distribution. AIMS Mathematics, 2024, 9(9): 23268-23290. doi: 10.3934/math.20241131 |
[5] | Hanan Haj Ahmad, Osama E. Abo-Kasem, Ahmed Rabaiah, Ahmed Elshahhat . Survival analysis of newly extended Weibull data via adaptive progressive Type-Ⅱ censoring and its modeling to Carbon fiber and electromigration. AIMS Mathematics, 2025, 10(4): 10228-10262. doi: 10.3934/math.2025466 |
[6] | Huifang Yuan, Tao Jiang, Min Xiao . The ruin probability of a discrete risk model with unilateral linear dependent claims. AIMS Mathematics, 2024, 9(4): 9785-9807. doi: 10.3934/math.2024479 |
[7] | Refah Alotaibi, Mazen Nassar, Zareen A. Khan, Wejdan Ali Alajlan, Ahmed Elshahhat . Entropy evaluation in inverse Weibull unified hybrid censored data with application to mechanical components and head-neck cancer patients. AIMS Mathematics, 2025, 10(1): 1085-1115. doi: 10.3934/math.2025052 |
[8] | Amal S. Hassan, Najwan Alsadat, Oluwafemi Samson Balogun, Baria A. Helmy . Bayesian and non-Bayesian estimation of some entropy measures for a Weibull distribution. AIMS Mathematics, 2024, 9(11): 32646-32673. doi: 10.3934/math.20241563 |
[9] | M. G. M. Ghazal . Modified Chen distribution: Properties, estimation, and applications in reliability analysis. AIMS Mathematics, 2024, 9(12): 34906-34946. doi: 10.3934/math.20241662 |
[10] | Lin Xu, Linlin Wang, Hao Wang, Liming Zhang . Optimal investment game for two regulated players with regime switching. AIMS Mathematics, 2024, 9(12): 34674-34704. doi: 10.3934/math.20241651 |
This article studies the questions of existence and nonexistence of weak solutions to the system of polyharmonic wave inequalities
{utt+(−Δ)mu≥|x|a|v|p,(t,x)∈(0,∞)×RN∖¯B1,vtt+(−Δ)mv≥|x|b|u|q,(t,x)∈(0,∞)×RN∖¯B1. | (1.1) |
Here, (u,v)=(u(t,x),v(t,x)), N≥2, B1 is the open unit ball of RN, m≥1 is an integer, a,b≥−2m, (a,b)≠(−2m,−2m), and p,q>1. We will investigate (1.1) under the Navier-type boundary conditions
{(−Δ)iu≥fi(x),i=0,⋯,m−1,(t,x)∈(0,∞)×∂B1,(−Δ)iv≥gi(x),i=0,⋯,m−1,(t,x)∈(0,∞)×∂B1, | (1.2) |
where fi,gi∈L1(∂B1) and (−Δ)0 is the identity operator. Notice that no restriction on the signs of fi or gi is imposed.
The study of semilinear wave inequalities in RN was firstly considered by Kato [1] and Pohozaev & Véron [2]. It was shown that the problem
utt−Δu≥|u|p,(t,x)∈(0,∞)×RN | (1.3) |
possesses a critical exponent pK=N+1N−1 in the following sense:
(ⅰ) If N≥2 and 1<p≤pK, then (1.3) possesses no global weak solution, provided
∫RNut(0,x)dx>0. | (1.4) |
(ⅱ) If p>pK, there are global positive solutions satisfying (1.4).
Caristi [3] studied the higher-order evolution polyharmonic inequality
∂ju∂tj−|x|αΔmu≥|u|p,(t,x)∈(0,∞)×RN, | (1.5) |
where α≤2m. Caristi discussed separately the cases α=2m and α<2m. For instance, when j=2 and α=0, it was shown that, if N≥m+1 and 1<p≤N+mN−m, then (1.5) possesses no global weak solution, provided (1.4) holds. Other existence and nonexistence results for evolution inequalities involving the polyharmonic operator in the whole space can be found in [4,5,6].
The study of the blow-up for semilinear wave equations in exterior domains was firstly considered by Zhang [7]. Namely, among many other problems, Zhang investigated the equation
utt−Δu=|x|a|u|p,(t,x)∈(0,∞)×RN∖D, | (1.6) |
where N≥3, a>−2, and D is a smooth bounded subset of RN. It was shown that (1.6) under the Neumann boundary condition
∂u∂ν=f(x)≥0,(t,x)∈(0,∞)×∂D, |
admits a critical exponent N+aN−2 in the following sense:
(ⅰ) If 1<p<N+aN−2, then (1.6) admits no global solution, provided f≢0.
(ⅱ) If p>N+aN−2, then (1.6) admits global solutions for some f>0.
In [8,9], it was shown that the critical value p=N+aN−2 belongs to case (ⅰ). Furthermore, the same result holds true, if (1.6) is considered under the Dirichlet boundary condition
u=f(x)≥0,(t,x)∈(0,∞)×∂D, |
where D=¯B1.
In [10], the authors considered the system of wave inequalities (1.1) in the case m=1. The system was studied under different types of inhomogeneous boundary conditions. In particular, under the boundary conditions (1.2) with m=1 (Dirichlet-type boundary conditions), the authors obtained the following result: Assume that a,b≥−2, (a,b)≠(−2,−2), If0:=∫∂B1f0dSx≥0, Ig0:=∫∂B1g0dSx≥0, (If0,Ig0)≠(0,0), and p,q>1. If N=2; or N≥3 and
N<max{sign(If0)2p(q+1)+pb+apq−1,sign(Ig0)2q(p+1)+qa+bpq−1}, |
then (1.1)-(1.2) (with m=1) admits no weak solution. Moreover, the authors pointed out the sharpness of the above condition.
In the case m=2, the system (1.1) was recently studied in [11] under different types of boundary conditions. In particular, under the boundary conditions (1.2) with f0≡0 and g0≡0, i.e.,
{u≥0,−Δu≥f1(x),(t,x)∈(0,∞)×∂B1,v≥0,−Δv≥g1(x),(t,x)∈(0,∞)×∂B1. | (1.7) |
Namely, the following result was obtained: Let N≥2, a,b≥−4, (a,b)≠(−4,−4), ∫∂B1f1dSx>0, ∫∂B1g1dSx>0, and p,q>1. If N∈{2,3,4}; or
N≥5,N<max{4p(q+1)+pb+apq−1,4q(p+1)+qa+bpq−1}, |
then (1.1) (with m=2) under the boundary conditions (1.7) admits no weak solution. Moreover, it was shown that the above condition is sharp.
Further results related to the existence and nonexistence of solutions for evolution problems in exterior domains can be found in [12,13,14,15,16,17].
The present work aims to extend the obtained results in [10,11] from m∈{1,2} to an arbitrary m≥1. Before presenting our main results, we need to define weak solutions to the considered problem.
Let
Q=(0,∞)×RN∖B1,ΣQ=(0,∞)×∂B1. |
Notice that ΣQ⊂Q.
Definition 1.1. We say that φ is an admissible test function, if
(i) φ∈C2,2mt,x(Q);
(ii) supp(φ)⊂⊂Q (φ is compactly supported in Q);
(iii) φ≥0;
(iv) For all j=0,1,⋯,m−1,
Δjφ|ΣQ=0,(−1)j∂(Δjφ)∂ν|ΣQ≤0, |
where ν denotes the outward unit normal vector on ∂B1, relative to RN∖B1.
The set of all admissible test functions is denoted by Φ.
Definition 1.2. We say that the pair (u,v) is a weak solution to (1.1)-(1.2), if
(u,v)∈Lqloc(Q)×Lploc(Q),∫Q|x|a|v|pφdxdt−m−1∑i=0∫ΣQfi(x)∂((−Δ)m−1−iφ)∂νdσdt≤∫Qu(−Δ)mφdxdt+∫Quφttdxdt | (1.8) |
and
∫Q|x|b|u|qφdxdt−m−1∑i=0∫ΣQgi∂((−Δ)m−1−iφ)∂νdσdt≤∫Qv(−Δ)mφdxdt+∫Qvφttdxdt | (1.9) |
for every φ∈Φ.
Notice that, if (u,v) is a regular solution to (1.1)-(1.2), then (u,v) is a weak solution in the sense of Definition 1.2.
For every function f∈L1(∂B1), we set
If=∫∂B1f(x)dσ. |
Our first main result is stated in the following theorem.
Theorem 1.1. Let p,q>1, N≥2, and a,b≥−2m with (a,b)≠(−2m,−2m). Let fi,gi∈L1(∂B1) for every i=0,⋯,m−1. Assume that Ifm−1,Igm−1≥0 and (Ifm−1,Igm−1)≠(0,0). If N≤2m; or N≥2m+1 and
N<max{sign(Ifm−1)×2mp(q+1)+pb+apq−1,sign(Igm−1)×2mq(p+1)+qa+bpq−1}, | (1.10) |
then (1.1)-(1.2) possesses no weak solution.
Remark 1.1. Notice that (1.10) is equivalent to
N−2m<α,Ifm−1>0; orN−2m<β,Igm−1>0, | (1.11) |
where
α=a+2m+p(b+2m)pq−1 | (1.12) |
and
β=b+2m+q(a+2m)pq−1. | (1.13) |
On the other hand, due to the condition a,b≥−2m and (a,b)≠(−2m,−2m), we have α,β>0, which shows that, if N≤2m, then (1.10) is always satisfied.
The proof of Theorem 1.1 is based on the construction of a suitable admissible test function and integral estimates. The construction of the admissible test function is specifically adapted to the polyharmonic operator (−Δ)m, the geometry of the domain, and the Navier-type boundary conditions (1.2).
Remark 1.2. By Theorem 1.1, we recover the nonexistence result obtained in [10] in the case m=1. We also recover the nonexistence result obtained in [11] in the case m=2.
Next, we are concerned with the existence of solutions to (1.1)-(1.2). Our second main result shows the sharpness of condition (1.10).
Theorem 1.2. Let p,q>1 and a,b≥−2m with (a,b)≠(−2m,−2m). If
N−2m>max{α,β}, | (1.14) |
where α and β are given by (1.12) and (1.13), then (1.1)-(1.2) admits stationary solutions for some fi,gi∈L1(∂B1) (i=0,⋯,m−1) with Ifm−1,Igm−1>0.
Theorem 1.2 will be proved by the construction of explicit stationary solutions to (1.1)-(1.2).
Remark 1.3. At this moment, we don't know whether there is existence or nonexistence in the critical case N≥2m+1,
N=max{sign(Ifm−1)×2mp(q+1)+pb+apq−1,sign(Igm−1)×2mq(p+1)+qa+bpq−1}. |
This question is left open.
From Theorem 1.1, we deduce the following nonexistence result for the corresponding stationary polyharmonic system
{(−Δ)mu≥|x|a|v|p,x∈RN∖¯B1,(−Δ)mv≥|x|b|u|q,x∈RN∖¯B1, | (1.15) |
under the Navier-type boundary conditions
{(−Δ)iu≥fi(x),i=0,⋯,m−1,x∈∂B1,(−Δ)iv≥gi(x),i=0,⋯,m−1,x∈∂B1. | (1.16) |
Corollary 1.1. Let p,q>1, N≥2, and a,b≥−2m with (a,b)≠(−2m,−2m). Let fi,gi∈L1(∂B1) for every i=0,⋯,m−1. Assume that Ifm−1,Igm−1≥0 and (Ifm−1,Igm−1)≠(0,0). If N≤2m; or N≥2m+1 and (1.10) holds, then (1.15)-(1.16) possesses no weak solution.
The rest of this manuscript is organized as follows: Section 2 is devoted to some auxiliary results. Namely, we first construct an admissible test function in the sense of Definition 1.1. Next, we establish some useful integral estimates involving the constructed test function. The proofs of Theorems 1.1 and 1.2 are provided in Section 3.
Throughout this paper, the letter C denotes a positive constant that is independent of the scaling parameters T, τ, and the solution (u,v). The value of C is not necessarily the same from one line to another.
In this section, we establish some auxiliary results that will be used later in the proof of our main result.
Let us introduce the radial function H defined in RN∖B1 by
H(x)={ln|x|ifN=2,1−|x|2−NifN≥3. | (2.1) |
We collect below some useful properties of the function H.
Lemma 2.1. The function H satisfies the following properties:
(i) H≥0;
(ii) H∈C2m(RN∖B1);
(iii) H|∂B1=0;
(iv) ΔH=0 in RN∖B1;
(v) For all j≥1,
ΔjH|∂B1=∂(ΔjH)∂ν|∂B1=0; |
(vi) ∂H∂ν|∂B1=−C.
Proof. (ⅰ)–(ⅴ) follow immediately from (2.1). On the other hand, we have
∂H∂ν|∂B1={−1ifN=2,−(N−2)ifN≥3, |
which proves (ⅵ).
We next consider a cut-off function ξ∈C∞(R) satisfying the following properties:
0≤ξ≤1,ξ(s)=1 if |s|≤1,ξ(s)=0 if |s|≥2. | (2.2) |
For all τ≫1, let
ξτ(x)=ξ(|x|τ),x∈RN∖B1, |
that is (from (2.2)),
ξτ(x)={1if1≤|x|≤τ,ξ(|x|τ)ifτ≤|x|≤2τ,0if|x|≥2τ. | (2.3) |
For k≫1, we introduce the function
ζτ(x)=H(x)ξkτ(x),x∈RN∖B1. | (2.4) |
We now introduce a second cut-off function G∈C∞(R) satisfying the following properties:
G≥0,supp(G)⊂⊂(0,1). | (2.5) |
For T>0 and k≫1, let
GT(t)=Gk(tT),t≥0. | (2.6) |
Let φ be the function defined by
φ(t,x)=GT(t)ζτ(x),(t,x)∈Q. | (2.7) |
By Lemma 2.1, (2.3)–(2.7), we obtain the following result.
Lemma 2.2. The function φ belongs to Φ.
For all λ>1, μ≥−2m, and φ∈Φ, we consider the integral terms
J(λ,μ,φ)=∫Q|x|−μλ−1φ−1λ−1|(−Δ)mφ|λλ−1dxdt | (2.8) |
and
K(λ,μ,φ)=∫Q|x|−μλ−1φ−1λ−1|φtt|λλ−1dxdt. | (2.9) |
Lemma 2.3. Let φ be the admissible test function defined by (2.7). Assume that
(i) J(p,a,φ),J(q,b,φ),K(p,a,φ),K(q,b,φ)<∞;
(ii) Ifm−1,Igm−1≥0.
If (u,v) is a weak solution to (1.1)-(1.2), then
Ifm−1≤CT−1([J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p)ppq−1([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q)pqpq−1 | (2.10) |
and
Igm−1≤CT−1([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q)qpq−1([J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p)pqpq−1. | (2.11) |
Proof. Let (u,v) be a weak solution to (1.1)-(1.2) and φ be the admissible test function defined by (2.7). By (1.8), we have
∫Q|x|a|v|pφdxdt−m−1∑i=0∫ΣQfi(x)∂((−Δ)m−1−iφ)∂νdσdt≤∫Qu(−Δ)mφdxdt+∫Quφttdxdt. |
On the other hand, by Lemma 2.1: (ⅴ), (ⅵ), (2.5)–(2.7), we have
m−1∑i=0∫ΣQfi(x)∂((−Δ)m−1−iφ)∂νdσdt=∫ΣQfm−1(x)∂φ∂νdσdt=−C∫ΣQfm−1(x)GT(t)dσdt=−C(∫∞0GT(t)dt)∫∂B1fm−1(x)dσ=−C(∫∞0Gk(tT)dt)Ifm−1=−CT(∫10Gk(s)ds)Ifm−1=−CTIfm−1. |
Consequently, we obtain
∫Q|x|a|v|pφdxdt+CTIfm−1≤∫Qu(−Δ)mφdxdt+∫Quφttdxdt. | (2.12) |
Similarly, by (1.9), we obtain
∫Q|x|b|u|qφdxdt+CTIgm−1≤∫Qv(−Δ)mφdxdt+∫Qvφttdxdt. | (2.13) |
Furthermore, by Hölder's inequality, we have
∫Qu(−Δ)mφdxdt≤∫Q|u||(−Δ)mφ|dxdt=∫Q(|x|bq|u|φ1q)(|x|−bq|(−Δ)mφ|φ−1q)dxdt≤(∫Q|x|b|u|qφdxdt)1q(∫Q|x|−bq−1|(−Δ)mφ|qq−1φ−1q−1dxdt)q−1q, |
that is,
∫Qu(−Δ)mφdxdt≤(∫Q|x|b|u|qφdxdt)1q[J(q,b,φ)]q−1q. | (2.14) |
Similarly, we obtain
∫Quφttdxdt≤(∫Q|x|b|u|qφdxdt)1q[K(q,b,φ)]q−1q. | (2.15) |
Thus, it follows from (2.12), (2.14), and (2.15) that
∫Q|x|a|v|pφdxdt+CTIfm−1≤(∫Q|x|b|u|qφdxdt)1q([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q). | (2.16) |
Using (2.13) and proceeding as above, we obtain
∫Q|x|b|u|qφdxdt+CTIgm−1≤(∫Q|x|a|v|pφdxdt)1p([J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p). | (2.17) |
Using (2.16)-(2.17), and taking into consideration that Igm−1≥0, we obtain
∫Q|x|a|v|pφdxdt+CTIfm−1≤(∫Q|x|a|v|pφdxdt)1pq([J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p)1q([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q). |
Then, by Young's inequality, it holds that
∫Q|x|a|v|pφdxdt+CTIfm−1≤1pq∫Q|x|a|v|pφdxdt+pq−1pq([J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p)pqq(pq−1)([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q)pqpq−1. |
Consequently, we have
(1−1pq)∫Q|x|a|v|pφdxdt+CTIfm−1≤pq−1pq([J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p)ppq−1([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q)pqpq−1, |
which yields (2.10). Similarly, using (2.16)-(2.17), and taking into consideration that Ifm−1≥0, we obtain
∫Q|x|b|u|qφdxdt+CTIgm−1≤(∫Q|x|b|u|qφdxdt)1pq([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q)1p([J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p), |
which implies by Young's inequality that
∫Q|x|b|u|qφdxdt+CTIgm−1≤1pq∫Q|x|b|u|qφdxdt+pq−1pq([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q)pqp(pq−1)([J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p)pqpq−1. |
Thus, it holds that
(1−1pq)∫Q|x|b|u|qφdxdt+CTIgm−1≤pq−1pq([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q)qpq−1([J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p)pqpq−1, |
which yields (2.11).
The aim of this subsection is to estimate the integral terms J(λ,μ,φ) and K(λ,μ,φ), where λ>1, μ≥−2m, and φ is the admissible test function defined by (2.7) with τ,k≫1.
The following result follows immediately from (2.5) and (2.6).
Lemma 2.4. We have
∫∞0GT(t)dt=CT. |
Lemma 2.5. We have
∫∞0G−1λ−1T|d2GTdt2|λλ−1dt≤CT1−2λλ−1. | (2.18) |
Proof. By (2.5) and (2.6), we have
∫∞0G−1λ−1T|d2GTdt2|λλ−1dt=∫T0G−1λ−1T|d2GTdt2|λλ−1dt | (2.19) |
and
d2GTdt2(t)=kT−2Gk−2(tT)((k−1)G′2(tT)+G(tT)G″(tT)) |
for all t∈(0,T). The above inequality yields
|d2GTdt2(t)|≤CT−2Gk−2(tT),t∈(0,T), |
which implies that
G−1λ−1T|d2GTdt2|λλ−1≤CT−2λλ−1Gk−2λλ−1(tT),t∈(0,T). |
Then, by (2.19), it holds that
∫∞0G−1λ−1T|d2GTdt2|λλ−1dt≤CT−2λλ−1∫T0Gk−2λλ−1(tT)dt=CT1−2λλ−1∫10Gk−2λλ−1(s)ds=CT1−2λλ−1, |
which proves (2.18).
To estimate J(λ,μ,φ) and K(λ,μ,φ), we consider separately the cases N≥3 and N=2.
Lemma 2.6. We have
∫RN∖B1|x|−μλ−1ζ−1λ−1τ|(−Δ)mζτ|λλ−1dx≤CτN−μ+2mλλ−1. | (2.20) |
Proof. Since H and ξτ are radial functions (see (2.1) and (2.3)), to simplify writing, we set
H(x)=H(r),ξτ(x)=ξτ(r), |
where r=|x|. By (2.4) and making use of Lemma 2.1 (ⅳ), one can show that for all x∈RN∖B1, we have
Δmζτ(x)=Δm(H(x)ξkτ(x))=2m−1∑i=0diHdri(r)2m−i∑j=1Ci,jdjξkτdrj(r)ri+j−2m, |
where Ci,j are some constants, which implies by (2.3) that
supp(Δmζτ)⊂{x∈RN:τ≤|x|≤2τ} | (2.21) |
and
|Δmζτ(x)|≤C2m−1∑i=0|diHdri(r)|2m−i∑j=1|djξkτdrj(r)|ri+j−2m,x∈supp(Δmζτ). | (2.22) |
On the other hand, for all x∈supp(Δmζτ), we have by (2.1) and (2.3) that
|diHdri(r)|={H(r)ifi=0,Cr2−N−iifi=1,⋯,2m−1 | (2.23) |
and (we recall that 0≤ξτ≤1)
|djξkτdrj(r)|≤Cτ−jξk−jτ(r)≤Cτ−jξk−2mτ(r),j=1,⋯,2m−i. | (2.24) |
Then, in view of (2.1), (2.21)–(2.24), we have
|Δmζτ(x)|≤Cξk−2mτ(r)(H(r)2m∑j=1τ−jrj−2m+r2−N2m−1∑i=12m−i∑j=1τ−jrj−2m)≤Cξk−2mτ(r)(τ−2m+τ2−N−2m)≤Cτ−2mξk−2mτ(x) |
for all x∈supp(Δmζτ). Taking into consideration that H≥C for all x∈supp(Δmζτ), the above estimate yields
|x|−μλ−1ζ−1λ−1τ|(−Δ)mζτ|λλ−1≤Cτ−2mλ−μλ−1ξk−2mλλ−1τ(x),x∈supp(Δmζτ). | (2.25) |
Finally, by (2.21) and (2.25), we obtain
∫RN∖B1|x|−μλ−1ζ−1λ−1τ|(−Δ)mζτ|λλ−1dx=∫τ<|x|<2τ|x|−μλ−1ζ−1λ−1τ|(−Δ)mζτ|λλ−1dx≤Cτ−2mλ−μλ−1∫τ<|x|<2τξk−2mλλ−1τ(x)dx≤Cτ−2mλ−μλ−1∫2τr=τrN−1dr=CτN−μ+2mλλ−1, |
which proves (2.20).
Lemma 2.7. We have
J(λ,μ,φ)≤CTτN−μ+2mλλ−1. |
Proof. By (2.7) and (2.8), we have
J(λ,μ,φ)=(∫∞0GT(t)dt)(∫RN∖B1|x|−μλ−1ζ−1λ−1τ|(−Δ)mζτ|λλ−1dx). |
Then, using Lemmas 2.4 and 2.6, we obtain the desired estimate.
Lemma 2.8. We have
∫RN∖B1|x|−μλ−1ζτ(x)dx≤C(τN−μλ−1+lnτ). | (2.26) |
Proof. By (2.1)–(2.4), we have
∫RN∖B1|x|−μλ−1ζτ(x)dx=∫1<|x|<2τ|x|−μλ−1(1−|x|2−N)ξκ(|x|τ)dx≤∫1<|x|<2τ|x|−μλ−1dx=C∫2τr=1rN−1−μλ−1dr≤{CτN−μλ−1ifN−μλ−1>0,[4pt]ClnτifN−μλ−1=0,[4pt]CifN−μλ−1<0≤C(τN−μλ−1+lnτ), |
which proves (2.26).
Lemma 2.9. We have
K(λ,μ,φ)≤CT1−2λλ−1(τN−μλ−1+lnτ). |
Proof. By (2.7) and (2.9), we have
K(λ,μ,φ)=(∫∞0G−1λ−1T|d2GTdt2|λλ−1dt)(∫RN∖B1|x|−μλ−1ζτ(x)dx). |
Then, using Lemmas 2.5 and 2.7, we obtain the desired estimate.
Lemma 2.10. We have
∫R2∖B1|x|−μλ−1ζ−1λ−1τ|(−Δ)mζτ|λλ−1dx≤Cτ2−2mλ+μλ−1lnτ. | (2.27) |
Proof. Proceeding as in the proof of Lemma 2.6, we obtain
supp(Δmζτ)⊂{x∈R2:τ≤|x|≤2τ} |
and
|Δmζτ(x)|≤Cτ−2mlnτξk−2mτ(x),x∈supp(Δmζτ). |
The above estimate yields
|x|−μλ−1ζ−1λ−1τ|(−Δ)mζτ|λλ−1≤Cτ−2mλ−μλ−1lnτξk−2mλλ−1τ(x),x∈supp(Δmζτ). |
Then, it holds that
∫R2∖B1|x|−μλ−1ζ−1λ−1τ|(−Δ)mζτ|λλ−1dx≤Cτ−2mλ−μλ−1lnτ∫τ<|x|<2τξk−2mλλ−1τ(x)dx≤Cτ−2mλ−μλ−1lnτ∫2τr=τrdr≤Cτ2−2mλ+μλ−1lnτ, |
which proves (2.27).
Using (2.7)-(2.8), Lemma 2.4, and Lemma 2.10, we obtain the following estimate of J(λ,μ,φ).
Lemma 2.11. We have
J(λ,μ,φ)≤CTτ2−2mλ+μλ−1lnτ. |
Lemma 2.12. We have
∫R2∖B1|x|−μλ−1ζτ(x)dx≤Clnτ(τ2−μλ−1+lnτ). | (2.28) |
Proof. By (2.1)–(2.4), we have
∫R2∖B1|x|−μλ−1ζτ(x)dx=∫1<|x|<2τ|x|−μλ−1ln|x|ξκ(|x|τ)dx≤∫1<|x|<2τ|x|−μλ−1ln|x|dx=C∫2τr=1r1−μλ−1lnrdr≤{Cτ2−μλ−1lnτif2−μλ−1>0,[4pt]C(lnτ)2if2−μλ−1=0,[4pt]Clnτif2−μλ−1<0≤Clnτ(τ2−μλ−1+lnτ), |
which proves (2.28).
Using (2.7), (2.9), Lemma 2.5, and Lemma 2.12, we obtain the following estimate of K(λ,μ,φ).
Lemma 2.13. We have
K(λ,μ,φ)≤CT1−2λλ−1lnτ(τ2−μλ−1+lnτ). |
This section is devoted to the proofs of Theorems 1.1 and 1.2.
By Remark 1.1, (1.10) is equivalent to (1.11). Without restriction of the generality, we assume that
N−2m<α,Ifm−1>0. | (3.1) |
Indeed, exchanging the roles of (Ifm−1,a,p) and (Igm−1,b,q), the case
N−2m<β,Igm−1>0 |
reduces to (3.1).
We use the contradiction argument. Namely, let us suppose that (u,v) is a weak solution to (1.1)-(1.2) (in the sense of Definition 1.2). For k,T,τ≫1, let φ be the admissible test function defined by (2.7). Then, by Lemma 2.3, we have
Ipq−1pfm−1≤CT−pq−1p([J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p)([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q)q. | (3.2) |
Making use of Lemmas 2.7 and 2.12, we obtain that for all N≥2,
J(λ,μ,φ)≤CTτN−μ+2mλλ−1lnτ,λ>1,μ≥−2m. | (3.3) |
Similarly, by Lemmas 2.9 and 2.13, we obtain that for all N≥2,
K(λ,μ,φ)≤CT1−2λλ−1(τN−μλ−1+lnτ)lnτ,λ>1,μ≥−2m. | (3.4) |
In particular, for (λ,μ)=(p,a), we obtain by (3.3) and (3.4) that
[J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p≤C[Tp−1pτ(N−a+2mpp−1)p−1p(lnτ)p−1p+T(1−2pp−1)p−1p(τN−ap−1+lnτ)p−1p(lnτ)p−1p]=CTp−1pτ(N−a+2mpp−1)p−1p(lnτ)p−1p[1+T−2(τ2mpp−1+τ−(N−a+2mpp−1)lnτ)p−1p]. | (3.5) |
Furthermore, taking T=τθ, where
θ>max{m,(a+2mpp−1−N)p−1p}, | (3.6) |
we obtain
1+T−2(τ2mpp−1+τ−(N−a+2mpp−1)lnτ)p−1p≤C. |
Then, from (3.5), we deduce that
[J(p,a,φ)]p−1p+[K(p,a,φ)]p−1p≤C[τθ+N−a+2mpp−1lnτ]p−1p. | (3.7) |
Similarly, for
θ>max{m,(b+2mqq−1−N)q−1q}, | (3.8) |
we obtain
([J(q,b,φ)]q−1q+[K(q,b,φ)]q−1q)q≤C[τθ+N−b+2mqq−1lnτ]q−1. | (3.9) |
Thus, for T=τθ, where θ satisfies (3.6) and (3.8), we obtain by (3.2), (3.7), and (3.9) that
Ipq−1pfm−1≤Cτ−θ(pq−1)p[τθ+N−a+2mpp−1lnτ]p−1p[τθ+N−b+2mqq−1lnτ]q−1, |
that is,
Ipq−1pfm−1≤Cτδ(lnτ)pq−1p, | (3.10) |
where
δ=pq−1p[N−(b+2mq)p+a+2mppq−1]=pq−1p(N−2m−α). |
Since N−2m<α, we have δ<0. Then, passing to the limit as τ→∞ in (3.10), we reach a contradiction with Ifm−1>0. This completes the proof of Theorem 1.1.
Let us introduce the family of polynomial functions {Pi}0≤i≤m, where
Pi(z)={1ifi=0,i−1∏j=0(z+2j)i∏j=1(N−2j−z)ifi=1,⋯,m. |
From (1.14), we deduce that
N−2j>max{α,β},j=1,⋯,m. |
Furthermore, because a,b≥−2m and (a,b)≠(−2m,−2m), we have α,β>0. Then,
Pi(z)>0,i=0,1,⋯,m,z∈{α,β}. | (3.11) |
For all
0<ε≤min{[Pm(α)]1p−1,[Pm(β)]1q−1}, | (3.12) |
we consider functions of the forms
uε(x)=ε|x|−α,x∈RN∖B1 | (3.13) |
and
vε(x)=ε|x|−β,x∈RN∖B1. | (3.14) |
Since uε and vε are radial functions, elementary calculations show that
(−Δ)iuε(x)=εPi(α)|x|−α−2i,i=0,1,⋯,m,x∈RN∖B1 | (3.15) |
and
(−Δ)ivε(x)=εPi(β)|x|−β−2i,i=0,1,⋯,m,x∈RN∖B1. | (3.16) |
Taking i=m in (3.15), using (3.11)–(3.14), we obtain
(−Δ)muε(x)=εPm(α)|x|−α−2m=|x|aεp|x|−βp(ε1−pPm(α)|x|−α−2m−a+βp)≥|x|avpε(x)|x|−α−2m−a+βp. |
On the other hand, by (1.12) and (1.13), one can show that
−α−2m−a+βp=0. |
Then, we obtain
(−Δ)muε(x)≥|x|avpε(x),x∈RN∖B1. | (3.17) |
Similarly, taking m=i in (3.16), using (3.11)–(3.14), we obtain
(−Δ)mvε(x)=εPm(β)|x|−β−2m=|x|bεq|x|−αq(ε1−qPm(β)|x|−β−2m−b+αq)≥|x|buqε(x)|x|−β−2m−b+αq. |
Using that
−β−2m−b+αq=0, |
we obtain
(−Δ)mvε(x)≥|x|buqε(x),x∈RN∖B1. | (3.18) |
Furthermore, by (3.11) and (3.15), for all i=0,⋯,m−1, we have
(−Δ)iuε(x)=εPi(α)>0,x∈∂B1. | (3.19) |
Similarly, by (3.11) and (3.16), for all i=0,⋯,m−1, we have
(−Δ)ivε(x)=εPi(β)>0,x∈∂B1. | (3.20) |
Finally, (3.17)–(3.20) show that for all ε satisfying (3.12), the pair of functions (uε,vε) given by (3.13) and (3.14) is a stationary solution to (1.1)-(1.2) with fi≡εPi(α) and gi≡εPi(β) for all i=0,⋯,m−1. The proof of Theorem 1.2 is then completed.
The system of polyharmonic wave inequalities (1.1) under the inhomogeneous Navier-type boundary conditions (1.2) was investigated. First, we established a nonexistence criterium for the nonexistence of weak solutions (see Theorem 1.1). Namely, under condition (1.10), we proved that (1.1)-(1.2) possesses no weak solution, provided Ifm−1,Igm−1≥0 and (Ifm−1,Igm−1)≠(0,0). Next, we proved the sharpness of the obtained criterium (1.10) by showing that under condition (1.14), (1.1)-(1.2) possesses weak solutions (stationary solutions) for some fi,gi∈L1(∂B1) (i=0,⋯,m−1) with Ifm−1,Igm−1>0 (see Theorem 1.2). From Theorem 1.1, we deduced an optimal criterium for the nonexistence of weak solutions to the corresponding stationary polyharmonic system (1.15) under the Navier-type boundary conditions (1.16) (see Corollary 1.1).
In this study, the critical case N≥2m+1,
N=max{sign(Ifm−1)×2mp(q+1)+pb+apq−1,sign(Igm−1)×2mq(p+1)+qa+bpq−1} |
is not investigated. It would be interesting to know whether there is existence or nonexistence of weak solutions in this case.
Manal Alfulaij: validation, investigation, writing review and editing; Mohamed Jleli: Conceptualization, methodology, investigation and formal analysis; Bessem Samet: Conceptualization, methodology, validation and investigation. All authors have read and approved the final version of the manuscript for publication.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Mohamed Jleli is supported by Researchers Supporting Project number (RSP2024R57), King Saud University, Riyadh, Saudi Arabia.
The authors declare no conflicts of interest.
[1] |
Q. H. Alqifiary, S. Jung, Laplace transform and generalized Hyers-Ulam stability of linear differential equations, Electron. J. Differ. Eq., 2014 (2014), 1-11. doi: 10.1186/1687-1847-2014-1
![]() |
[2] | Q. H. Alqifiary, J. K. Miljanovic, Note on the stability of system of differential equations ˙x(t)=f(t,x(t)), Gen. Math. Notes, 20 (2014), 27-33. |
[3] | C. Alsina, R. Ger, On some inequalities and stability results related to the exponential function, J. Inequal. Appl., 2 (1998), 373-380. |
[4] |
T. Aoki, On the stability of the linear transformation in Banach spaces, J. Math. Soc. Jpn., 2 (1950), 64-66. doi: 10.2969/jmsj/00210064
![]() |
[5] |
R. Fukutaka, M. Onitsuka, Best constant in Hyers-Ulam stability of first-order homogeneous linear differential equations with a periodic coefficient, J. Math. Anal. Appl., 473 (2019), 1432-1446. doi: 10.1016/j.jmaa.2019.01.030
![]() |
[6] |
D. H. Hyers, On the stability of a linear functional equation, P. Natl. Acad. Sci. USA., 27 (1941), 222-224. doi: 10.1073/pnas.27.4.222
![]() |
[7] |
S. M. Jung, Hyers-Ulam stability of linear differential equations of first order, Appl. Math. Lett., 17 (2004), 1135-1140. doi: 10.1016/j.aml.2003.11.004
![]() |
[8] |
S. M. Jung, Hyers-Ulam stability of linear differential equations of first order (III), J. Math. Anal. Appl., 311 (2005), 139-146. doi: 10.1016/j.jmaa.2005.02.025
![]() |
[9] | S. M. Jung, On the quadratic functional equation modulo a subgroup, Indian J. Pure Appl. Math., 36 (2005), 441-450. |
[10] |
S. M. Jung, Hyers-Ulam stability of linear differential equations of first order (II), Appl. Math. Lett., 19 (2006), 854-858. doi: 10.1016/j.aml.2005.11.004
![]() |
[11] |
S. M. Jung, Hyers-Ulam stability of a system of first order linear differential equations with constant coefficients, J. Math. Anal. Appl., 320 (2006), 549-561. doi: 10.1016/j.jmaa.2005.07.032
![]() |
[12] | S. M. Jung, Approximate solution of a linear differential equation of third order, B. Malays. Math. Sci. So., 35 (2012), 1063-1073. |
[13] | V. Kalvandi, N. Eghbali, J. M. Rassias, Mittag-Leffler-Hyers-Ulam stability of fractional differential equations of second order, J. Math. Ext., 13 (2019), 29-43. |
[14] |
T. Li, A. Zada, S. Faisal, Hyers-Ulam stability of nth order linear differential equations, J. Nonlinear Sci. Appl., 9 (2016), 2070-2075. doi: 10.22436/jnsa.009.05.12
![]() |
[15] |
Y. Li, Y. Shen, Hyers-Ulam stability of linear differential equations of second order, Appl. Math. Lett., 23 (2010), 306-309. doi: 10.1016/j.aml.2009.09.020
![]() |
[16] | K. Liu, M. Feckan, D. O'Regan, et al. Hyers-Ulam stability and existence of solutions for differential equations with Caputo-Fabrizio fractional derivative, Mathematics, 7 (2019), 333. |
[17] | T. Miura, S. Jung, S. E. Takahasi, Hyers-Ulam-Rassias stability of the Banach space valued linear differential equation y′=λy, J. Korean Math. Soc., 41 (2004), 995-1005. |
[18] | R. Murali, A. P. Selvan, On the generalized Hyers-Ulam stability of linear ordinary differential equations of higher order, Int. J. Pure Appl. Math., 117 (2017), 317-326. |
[19] | M. Ramdoss, P. S. Arumugan, Fourier transforms and Ulam stabilities of linear differential equations, In: G. Anastassiou, J. Rassias, editors, Frontiers in Functional Equations and Analytic Inequalities, Springer, Cham, 2019, 195-217. |
[20] | M. Obloza, Hyers stability of the linear differential equation, Rockznik Nauk-Dydakt. Prace Mat., 13 (1993), 259-270. |
[21] | M. Obloza, Connection between Hyers and Lyapunov stability of the ordinary differential equations, Rockznik Nauk-Dydakt. Prace Mat., 14 (1997), 141-146. |
[22] |
M. Onitsuka, Hyers-Ulam stability of first order linear differential equations of Carathéodory type and its application, Appl. Math. Lett., 90 (2019), 61-68. doi: 10.1016/j.aml.2018.10.013
![]() |
[23] |
M. Onitsuka, T. Shoji, Hyers-Ulam stability of first order homogeneous linear differential equations with a real valued coefficients, Appl. Math. Lett., 63 (2017), 102-108. doi: 10.1016/j.aml.2016.07.020
![]() |
[24] |
T. M. Rassias, On the stability of the linear mappings in Banach spaces, P. Am. Math. Soc., 72 (1978), 297-300. doi: 10.1090/S0002-9939-1978-0507327-1
![]() |
[25] | I. A. Rus, Ulam stabilities of ordinary differential equations in Banach space, Carpathian J. Math., 26 (2010), 103-107. |
[26] | S. E. Takahasi, T. Miura, S. Miyajima, On the Hyers-Ulam stability of the Banach space-valued differential equation y′=αy, Bull. Korean Math. Soc., 39 (2002), 309-315. |
[27] | S. M. Ulam, A Collection of Mathematical Problems, Interscience Publishers, New York, 1960. |
[28] |
G. Wang, M. Zhou, L. Sun, Hyers-Ulam stability of linear differential equations of first order, Appl. Math. Lett., 21 (2008), 1024-1028. doi: 10.1016/j.aml.2007.10.020
![]() |
[29] |
J. R. Wang, A. Zada, W. Ali, Ulam's type stability of first-order impulsive differential equations with variable delay in quasi-Banach spaces, Int. J. Nonlin. Sci. Num., 19 (2018), 553-560. doi: 10.1515/ijnsns-2017-0245
![]() |
[30] | X. Wang, M. Arif, A. Zada, β-Hyers-Ulam-Rassias stability of semilinear nonautonomous impulsive system, Symmetry, 11 (2019), 231. |
[31] | A. Zada, W. Ali, C. Park, Ulam's type stability of higher order nonlinear delay differential equations via integral inequality of Gronwall Bellman-Bihari's type, Appl. Math. Comput., 350 (2019), 60-65. |
[32] | A. Zada, S. O. Shah, Hyers-Ulam stability of first-order non-linear delay differential equations with fractional integrable impulses, Hacet. J. Math. Stat., 47 (2018), 1196-1205. |
[33] |
A. Zada, S. Shaleena, T. Li, Stability analysis of higher order nonlinear differential equations in β-normed spaces, Math. Method. Appl. Sci., 42 (2019), 1151-1166. doi: 10.1002/mma.5419
![]() |
[34] | A. Zada, M. Yar, T. Li, Existence and stability analysis of nonlinear sequential coupled system of Caputo fractional differential equations with integral boundary conditions, Ann. Univ. Paedagog. Crac. Stud. Math., 17 (2018), 103-125. |