Citation: Swarnali Sharma, Vitaly Volpert, Malay Banerjee. Extended SEIQR type model for COVID-19 epidemic and data analysis[J]. Mathematical Biosciences and Engineering, 2020, 17(6): 7562-7604. doi: 10.3934/mbe.2020386
[1] | Salim A. Messaoudi, Mohammad M. Al-Gharabli, Adel M. Al-Mahdi, Mohammed A. Al-Osta . A coupled system of Laplacian and bi-Laplacian equations with nonlinear dampings and source terms of variable-exponents nonlinearities: Existence, uniqueness, blow-up and a large-time asymptotic behavior. AIMS Mathematics, 2023, 8(4): 7933-7966. doi: 10.3934/math.2023400 |
[2] | Adel M. Al-Mahdi, Mohammad M. Al-Gharabli, Maher Nour, Mostafa Zahri . Stabilization of a viscoelastic wave equation with boundary damping and variable exponents: Theoretical and numerical study. AIMS Mathematics, 2022, 7(8): 15370-15401. doi: 10.3934/math.2022842 |
[3] | Xiongmei Fan, Sen Ming, Wei Han, Zikun Liang . Lifespan estimate of solution to the semilinear wave equation with damping term and mass term. AIMS Mathematics, 2023, 8(8): 17860-17889. doi: 10.3934/math.2023910 |
[4] | Mohammad Kafini, Jamilu Hashim Hassan, Mohammad M. Al-Gharabli . Decay result in a problem of a nonlinearly damped wave equation with variable exponent. AIMS Mathematics, 2022, 7(2): 3067-3082. doi: 10.3934/math.2022170 |
[5] | Adel M. Al-Mahdi . The coupling system of Kirchhoff and Euler-Bernoulli plates with logarithmic source terms: Strong damping versus weak damping of variable-exponent type. AIMS Mathematics, 2023, 8(11): 27439-27459. doi: 10.3934/math.20231404 |
[6] | Noufe H. Aljahdaly . Study tsunamis through approximate solution of damped geophysical Korteweg-de Vries equation. AIMS Mathematics, 2024, 9(5): 10926-10934. doi: 10.3934/math.2024534 |
[7] | Mohammad Kafini, Mohammad M. Al-Gharabli, Adel M. Al-Mahdi . Existence and stability results of nonlinear swelling equations with logarithmic source terms. AIMS Mathematics, 2024, 9(5): 12825-12851. doi: 10.3934/math.2024627 |
[8] | Zhiqiang Li . The finite time blow-up for Caputo-Hadamard fractional diffusion equation involving nonlinear memory. AIMS Mathematics, 2022, 7(7): 12913-12934. doi: 10.3934/math.2022715 |
[9] | Muhammad Shakeel, Amna Mumtaz, Abdul Manan, Marouan Kouki, Nehad Ali Shah . Soliton solutions of the nonlinear dynamics in the Boussinesq equation with bifurcation analysis and chaos. AIMS Mathematics, 2025, 10(5): 10626-10649. doi: 10.3934/math.2025484 |
[10] | Khaled Kefi, Nasser S. Albalawi . Three weak solutions for degenerate weighted quasilinear elliptic equations with indefinite weights and variable exponents. AIMS Mathematics, 2025, 10(2): 4492-4503. doi: 10.3934/math.2025207 |
We consider the Cauchy problem for the nonlinear damped wave equation with linear damping
utt−Δu+ut=|u|p+|∇u|q+w(x),t>0,x∈RN, | (1.1) |
where N≥1, p,q>1, w∈L1loc(RN), w≥0 and w≢0. Namely, we are concerned with the existence and nonexistence of global weak solutions to (1.1). We mention below some motivations for studying problems of type (1.1).
Consider the semilinear damped wave equation
utt−Δu+ut=|u|p,t>0,x∈RN | (1.2) |
under the initial conditions
(u(0,x),ut(0,x))=(u0(x),u1(x)),x∈RN. | (1.3) |
In [12], a Fujita-type result was obtained for the problem (1.2) and (1.3). Namely, it was shown that,
(ⅰ) if 1<p<1+2N and ∫RNuj(x)dx>0, j=0,1, then the problem (1.2) and (1.3) admits no global solution;
(ⅱ) if 1+2N<p<N for N≥3, and 1+2N<p<∞ for N∈{1,2}, then the problem (1.2) and (1.3) admits a unique global solution for suitable initial values.
The proof of part (ⅰ) makes use of the fundamental solution of the operator (∂tt−Δ+∂t)k. In [15], it was shown that the exponent 1+2N belongs to the blow-up case (ⅰ). Notice that 1+2N is also the Fujita critical exponent for the semilinear heat equation ut−Δu=|u|p (see [3]). In [6], the authors considered the problem
utt+(−1)m|x|αΔmu+ut=f(t,x)|u|p+w(t,x),t>0,x∈RN | (1.4) |
under the initial conditions (1.3), where m is a positive natural number, α≥0, and f(t,x)≥0 is a given function satisfying a certain condition. Using the test function method (see e.g. [11]), sufficient conditions for the nonexistence of a global weak solution to the problem (1.4) and (1.3) are obtained. Notice that in [6], the influence of the inhomogeneous term w(t,x) on the critical behavior of the problem (1.4) and (1.3) was not investigated. Recently, in [5], the authors investigated the inhomogeneous problem
utt−Δu+ut=|u|p+w(x),t>0,x∈RN, | (1.5) |
where w∈L1loc(RN) and w≢0. It was shown that in the case N≥3, the critical exponent for the problem (1.5) jumps from 1+2N (the critical exponent for the problem (1.2)) to the bigger exponent 1+2N−2. Notice that a similar phenomenon was observed for the inhomogeneous semilnear heat equation ut−Δu=|u|p+w(x) (see [14]). For other results related to the blow-up of solutions to damped wave equations, see, for example [1,2,4,7,8,9,13] and the references therein.
In this paper, motivated by the above mentioned works, our aim is to study the effect of the gradient term |∇u|q on the critical behavior of problem (1.5). We first obtain the Fujita critical exponent for the problem (1.1). Next, we determine its second critical exponent in the sense of Lee and Ni [10].
Before stating the main results, let us provide the notion of solutions to problem (1.1). Let Q=(0,∞)×RN. We denote by C2c(Q) the space of C2 real valued functions compactly supported in Q.
Definition 1.1. We say that u=u(t,x) is a global weak solution to (1.1), if
(u,∇u)∈Lploc(Q)×Lqloc(Q) |
and
∫Q(|u|p+|∇u|q)φdxdt+∫Qw(x)φdxdt=∫Quφttdxdt−∫QuΔφdxdt−∫Quφtdxdt, | (1.6) |
for all φ∈C2c(Q).
The first main result is the following Fujita-type theorem.
Theorem 1.1. Let N≥1, p,q>1, w∈L1loc(RN), w≥0 and w≢0.
(i) If N∈{1,2}, for all p,q>1, problem (1.1) admits no global weak solution.
(ii) Let N≥3. If
1<p<1+2N−2or1<q<1+1N−1, |
then problem (1.1) admits no global weak solution.
(iii) Let N≥3. If
p>1+2N−2andq>1+1N−1, |
then problem (1.1) admits global solutions for some w>0.
We mention below some remarks related to Theorem 1.1. For N≥3, let
pc(N,q)={1+2N−2ifq>1+1N−1,∞ifq<1+1N−1. |
From Theorem 1.1, if 1<p<pc(N,q), then problem (1.1) admits no global weak solution, while if p>pc(N,q), global weak solutions exist for some w>0. This shows that pc(N,q) is the Fujita critical exponent for the problem (1.1). Notice that in the case N≥3 and q>1+1N−1, pc(N,q) is also the critical exponent for utt−Δu=|u|p+w(x) (see [5]). This shows that in the range q>1+1N−1, the gradient term |∇u|q has no influence on the critical behavior of utt−Δu=|u|p+w(x).
It is interesting to note that the gradient term |∇u|q induces an interesting phenomenon of discontinuity of the Fujita critical exponent pc(N,q) jumping from 1+2N−2 to ∞ as q reaches the value 1+1N−1 from above.
Next, for N≥3 and σ<N, we introduce the sets
I+σ={w∈C(RN):w≥0,|x|−σ=O(w(x)) as |x|→∞} |
and
I−σ={w∈C(RN):w>0,w(x)=O(|x|−σ) as |x|→∞}. |
The next result provides the second critical exponent for the problem (1.1) in the sense of Lee and Ni [10].
Theorem 1.2. Let N≥3, p>1+2N−2 and q>1+1N−1.
(i) If
σ<max{2pp−1,qq−1} |
and w∈I+σ, then problem (1.1) admits no global weak solution.
(ii) If
σ≥max{2pp−1,qq−1}, |
then problem (1.1) admits global solutions for some w∈I−σ.
From Theorem 1.2, if N≥3, p>1+2N−2 and q>1+1N−1, then problem (1.1) admits a second critical exponent, namely
σ∗=max{2pp−1,qq−1}. |
The rest of the paper is organized as follows. In Section 2, some preliminary estimates are provided. Section 3 is devoted to the study of the Fujita critical exponent for the problem (1.1). Namely, Theorem 1.1 is proved. In Section 4, we study the second critical exponent for the problem (1.1) in the sense of Lee and Ni. Namely, we prove Theorem 1.2.
Consider two cut-off functions f,g∈C∞([0,∞)) satisfying
f≥0,f≢0,supp(f)⊂(0,1) |
and
0≤g≤1,g(s)={1if0≤s≤1,0ifs≥2. |
For T>0, we introduce the function
ξT(t,x)=f(tT)ℓg(|x|2T2ρ)ℓ=F(t)G(x),(t,x)∈Q, | (2.1) |
where ℓ≥2 and ρ>0 are constants to be chosen. It can be easily seen that
ξT≥0andξT∈C2c(Q). |
Throughout this paper, the letter C denotes various positive constants depending only on known quantities.
Lemma 2.1. Let κ>1 and ℓ>2κκ−1. Then
∫Q|ξT|−1κ−1|ΔξT|κκ−1dxdt≤CT−2ρκκ−1+Nρ+1. | (2.2) |
Proof. By (2.1), one has
∫Q|ξT|−1κ−1|ΔξT|κκ−1dxdt=(∫∞0F(t)dt)(∫RN|G|−1κ−1|ΔG|κκ−1dx). | (2.3) |
Using the properties of the cut-off function f, one obtains
∫∞0F(t)dt=∫∞0f(tT)ℓdt=∫T0f(tT)ℓdt=T∫10f(s)ℓds, |
which shows that
∫∞0F(t)dt=CT. | (2.4) |
Next, using the properties of the cut-off function g, one obtains
∫RN|G|−1κ−1|ΔG|κκ−1dx=∫Tρ<|x|<√2Tρ|G|−1κ−1|ΔG|κκ−1dx | (2.5) |
On the other hand, for Tρ<|x|<√2Tρ, an elementary calculation yields
ΔG(x)=Δ[g(r2T2ρ)ℓ],r=|x|=(d2dr2+N−1rddr)g(r2T2ρ)ℓ=2ℓT−2ρg(r2T2ρ)ℓ−2θ(x), |
where
θ(x)=Ng(r2T2ρ)g′(r2T2ρ)+2(ℓ−1)T−2ρr2g(r2T2ρ)g′(r2T2ρ)2+2T−2ρr2g(r2T2ρ)g″(r2T2ρ). |
Notice that since Tρ<r<√2Tρ, one deduces that
|θ(x)|≤C. |
Hence, it holds that
|ΔG(x)|≤CT−2ρg(|x|2T2ρ)ℓ−2,Tρ<|x|<√2Tρ |
and
|G(x)|−1κ−1|ΔG(x)|κκ−1≤CT−2ρκκ−1g(|x|2T2ρ)ℓ−2κκ−1,Tρ<|x|<√2Tρ. |
Therefore, by (2.5) and using the change of variable x=Tρy, one obtains
∫RN|G|−1κ−1|ΔG|κκ−1dx≤CT−2ρκκ−1∫Tρ<|x|<√2Tρg(|x|2T2ρ)ℓ−2κκ−1dx=CT−2ρκκ−1+Nρ∫1<|y|<√2g(|y|2)ℓ−2κκ−1dy, |
which yields (notice that ℓ>2κκ−1)
∫RN|G|−1κ−1|ΔG|κκ−1dx≤CT−2ρκκ−1+Nρ. | (2.6) |
Finally, using (2.3), (2.4) and (2.6), (2.2) follows.
Lemma 2.2. Let κ>1 and ℓ>κκ−1. Then
∫Q|ξT|−1κ−1|∇ξT|κκ−1dxdt≤CT−ρκκ−1+Nρ+1. | (2.7) |
Proof. By (2.1) and the properties of the cut-off function g, one has
∫Q|ξT|−1κ−1|∇ξT|κκ−1dxdt=(∫∞0F(t)dt)(∫Tρ<|x|<√2Tρ|G|−1κ−1|∇G|κκ−1dx). | (2.8) |
On the other hand, for Tρ<|x|<√2Tρ, an elementary calculation shows that
|∇G(x)|=2ℓT−2ρ|x||g′(|x|2T2ρ)|g(|x|2T2ρ)ℓ−1≤2√2ℓT−ρ|g′(|x|2T2ρ)|g(|x|2T2ρ)ℓ−1≤CT−ρg(|x|2T2ρ)ℓ−1, |
which yields
|G|−1κ−1|∇G|κκ−1≤CT−ρκκ−1g(|x|2T2ρ)ℓ−κκ−1,Tρ<|x|<√2Tρ. |
Therefore, using the change of variable x=Tρy and the fact that ℓ>κκ−1, one obtains
∫Tρ<|x|<√2Tρ|G|−1κ−1|∇G|κκ−1dx≤CT−ρκκ−1∫Tρ<|x|<√2Tρg(|x|2T2ρ)ℓ−κκ−1dx=CT−ρκκ−1+Nρ∫1<|y|<√2g(|y|2)ℓ−κκ−1dy, |
i.e.
∫Tρ<|x|<√2Tρ|G|−1κ−1|∇G|κκ−1dx≤CT−ρκκ−1+Nρ. | (2.9) |
Using (2.4), (2.8) and (2.9), (2.7) follows.
Lemma 2.3. Let κ>1 and ℓ>2κκ−1. Then
∫Q|ξT|−1κ−1|ξTtt|κκ−1dxdt≤CT−2κκ−1+1+Nρ. | (2.10) |
Proof. By (2.1) and the properties of the cut-off functions f and g, one has
∫Q|ξT|−1κ−1|ξTtt|κκ−1dxdt=(∫T0F(t)−1κ−1|F″(t)|κκ−1dt)(∫0<|x|<√2TρG(x)dx). | (2.11) |
An elementary calculation shows that
F(t)−1κ−1|F″(t)|κκ−1≤CT−2κκ−1f(tT)ℓ−2κκ−1,0<t<T, |
which yields (since ℓ>2κκ−1)
∫T0F(t)−1κ−1|F″(t)|κκ−1dt≤CT−2κκ−1+1. | (2.12) |
On the other hand, using the change of variable x=Tρy, one obtains
∫0<|x|<√2TρG(x)dx=∫0<|x|<√2Tρg(|x|2T2ρ)ℓdx=TNρ∫0<|y|<√2g(|y|2)ℓdy, |
which yields
∫0<|x|<√2TρG(x)dx≤CTNρ. | (2.13) |
Therefore, using (2.11), (2.12) and (2.13), (2.10) follows.
Lemma 2.4. Let κ>1 and ℓ>κκ−1. Then
∫Q|ξT|−1κ−1|ξTt|κκ−1dxdt≤CT−κκ−1+1+Nρ. | (2.14) |
Proof. By (2.1) and the properties of the cut-off functions f and g, one has
∫Q|ξT|−1κ−1|ξTt|κκ−1dxdt=(∫T0F(t)−1κ−1|F′(t)|κκ−1dt)(∫0<|x|<√2TρG(x)dx). | (2.15) |
An elementary calculation shows that
F(t)−1κ−1|F′(t)|κκ−1≤CT−κκ−1f(tT)ℓ−κκ−1,0<t<T, |
which yields (since ℓ>κκ−1)
∫T0F(t)−1κ−1|F′(t)|κκ−1dt≤CT1−κκ−1. | (2.16) |
Therefore, using (2.13), (2.15) and (2.16), (2.14) follows.
Proposition 3.1. Suppose that problem (1.1) admits a global weak solution u with (u,∇u)∈Lploc(Q)×Lqloc(Q). Then there exists a constant C>0 such that
∫Qw(x)φdxdt≤Cmin{A(φ),B(φ)}, | (3.1) |
for all φ∈C2c(Q), φ≥0, where
A(φ)=∫Qφ−1p−1|φtt|pp−1dxdt+∫Qφ−1p−1|φt|pp−1dxdt+∫Qφ−1p−1|Δφ|pp−1dxdt |
and
B(φ)=∫Qφ−1p−1|φtt|pp−1dxdt+∫Qφ−1p−1|φt|pp−1dxdt+∫Qφ−1q−1|∇φ|qq−1dxdt. |
Proof. Let u be a global weak solution to problem (1.1) with
(u,∇u)∈Lploc(Q)×Lqloc(Q). |
Let φ∈C2c(Q), φ≥0. Using (1.6), one obtains
∫Q|u|pφdxdt+∫Qw(x)φdxdt≤∫Q|u||φtt|dxdt+∫Q|u||Δφ|dxdt+∫Q|u||φt|dxdt. | (3.2) |
On the other hand, by ε-Young inequality, 0<ε<13, one has
∫Q|u||φtt|dxdt=∫Q(|u|φ1p)(φ−1p|φtt|)dxdt≤ε∫Q|u|pφdxdt+C∫Qφ−1p−1|φtt|pp−1dxdt. | (3.3) |
Here and below, C denotes a positive generic constant, whose value may change from line to line. Similarly, one has
∫Q|u||Δφ|dxdt≤ε∫Q|u|pφdxdt+C∫Qφ−1p−1|Δφ|pp−1dxdt | (3.4) |
and
∫Q|u||φt|dxdt≤ε∫Q|u|pφdxdt+C∫Qφ−1p−1|φt|pp−1dxdt | (3.5) |
Therefore, using (3.2), (3.3), (3.4) and (3.5), it holds that
(1−3ε)∫Q|u|pφdxdt+∫Qw(x)φdxdt≤CA(φ). |
Since 1−3ε>0, one deduces that
∫Qw(x)φdxdt≤CA(φ). | (3.6) |
Again, by (1.6), one has
∫Q(|u|p+|∇u|q)φdxdt+∫Qw(x)φdxdt=∫Quφttdxdt−∫QuΔφdxdt−∫Quφtdxdt=∫Quφttdxdt+∫Q∇u⋅∇φdxdt−∫Quφtdxdt, |
where ⋅ denotes the inner product in RN. Hence, it holds that
∫Q(|u|p+|∇u|q)φdxdt+∫Qw(x)φdxdt≤∫Q|u||φtt|dxdt+∫Q|∇u||∇φ|dxdt+∫Q|u||φt|dxdt. | (3.7) |
By ε-Young inequality (ε=12), one obtains
∫Q|u||φtt|dxdt≤12∫Q|u|pφdxdt+C∫Qφ−1p−1|φtt|pp−1dxdt, | (3.8) |
∫Q|u||φt|dxdt≤12∫Q|u|pφdxdt+C∫Qφ−1p−1|φt|pp−1dxdt, | (3.9) |
∫Q|∇u||∇φ|dxdt≤12∫Q|∇u|qφdxdt+C∫Qφ−1q−1|∇φ|qq−1dxdt. | (3.10) |
Next, using (3.7), (3.8), (3.9) and (3.10), one deduces that
12∫Q|∇u|qφdxdt+∫Qw(x)φdxdt≤CB(φ), |
which yields
∫Qw(x)φdxdt≤CB(φ). | (3.11) |
Finally, combining (3.6) with (3.11), (3.1) follows.
Proposition 3.2. Suppose that problem (1.1) admits a global weak solution u with (u,∇u)∈Lploc(Q)×Lqloc(Q). Then
∫Qw(x)ξTdxdt≤Cmin{A(ξT),B(ξT)}, | (3.12) |
for all T>0, where ξT is defined by (2.1).
Proof. Since for all T>0, ξT∈C2c(Q), ξT≥0, taking φ=ξT in (3.1), (3.12) follows.
Proposition 3.3. Let w∈L1loc(RN), w≥0 and w≢0. Then, for sufficiently large T,
∫Qw(x)ξTdxdt≥CT, | (3.13) |
where ξT is defined by (2.1).
Proof. By (2.1) and the properties of the cut-off functions f and g, one has
∫Qw(x)ξTdxdt=(∫T0f(tT)ℓdt)(∫0<|x|<√2Tρw(x)g(|x|2T2ρ)ℓdx). | (3.14) |
On the other hand,
∫T0f(tT)ℓdt=T∫10f(s)ℓds. | (3.15) |
Moreover, for sufficiently large T (since w,g≥0),
∫0<|x|<√2Tρw(x)g(|x|2T2ρ)ℓdx≥∫0<|x|<1w(x)g(|x|2T2ρ)ℓdx. | (3.16) |
Notice that since w∈L1loc(RN), w≥0, w≢0, and using the properties of the cut-off function g, by the domianted convergence theorem, one has
limT→∞∫0<|x|<1w(x)g(|x|2T2ρ)ℓdx=∫0<|x|<1w(x)dx>0. |
Hence, for sufficiently large T,
∫0<|x|<1w(x)g(|x|2T2ρ)ℓdx≥C. | (3.17) |
Using (3.14), (3.15), (3.16) and (3.17), (3.13) follows.
Now, we are ready to prove the Fujita-type result given by Theorem 1.1. The proof is by contradiction and makes use of the nonlinear capacity method, which was developed by Mitidieri and Pohozaev (see e.g. [11]).
Proof of Theorem 1.1. (ⅰ) Suppose that problem (1.1) admits a global weak solution u with (u,∇u)∈Lploc(Q)×Lqloc(Q). By Propositions 3.2 and 3.3, for sufficiently large T, one has
CT≤min{A(ξT),B(ξT)}, | (3.18) |
where ξT is defined by (2.1),
A(ξT)=∫Qξ−1p−1T|ξTtt|pp−1dxdt+∫Qξ−1p−1T|ξTt|pp−1dxdt+∫Qξ−1p−1T|ΔξT|pp−1dxdt | (3.19) |
and
B(ξT)=∫Qφ−1p−1|ξTtt|pp−1dxdt+∫Qξ−1p−1T|ξTt|pp−1dxdt+∫Qξ−1q−1T|∇ξT|qq−1dxdt. | (3.20) |
Taking ℓ>max{2pp−1,qq−1}, using Lemmas 2.1, 2.3, 2.4 with κ=p, and Lemma 2.2 with κ=q, one obtains the following estimates
∫Qξ−1p−1T|ξTtt|pp−1dxdt≤CT−2pp−1+1+Nρ, | (3.21) |
∫Qξ−1p−1T|ξTt|pp−1dxdt≤CT−pp−1+1+Nρ, | (3.22) |
∫Qξ−1p−1T|ΔξT|pp−1dxdt≤CT−2ρpp−1+Nρ+1, | (3.23) |
∫Qξ−1q−1T|∇ξT|qq−1dxdt≤CT−ρqq−1+Nρ+1. | (3.24) |
Hence, by (3.19), (3.21), (3.22) and (3.23), one deduces that
A(ξT)≤C(T−2pp−1+1+Nρ+T−pp−1+1+Nρ+T−2ρpp−1+Nρ+1). |
Observe that
−2pp−1+1+Nρ<−pp−1+1+Nρ. | (3.25) |
So, for sufficiently large T, one deduces that
A(ξT)≤C(T−pp−1+1+Nρ+T−2ρpp−1+Nρ+1). |
Notice that the above estimate holds for all ρ>0. In particular, when ρ=12, it holds that
A(ξT)≤CT−pp−1+1+N2 | (3.26) |
Next, by (3.20), (3.21), (3.22), (3.24) and (3.25), one deduces that
B(ξT)≤C(T−pp−1+1+Nρ+T−ρqq−1+Nρ+1). |
Similarly, the above inequality holds for all ρ>0. In particular, when ρ=p(q−1)q(p−1), it holds that
B(ξT)≤CT−pp−1+1+Np(q−1)q(p−1). | (3.27) |
Therefore, it follows from (3.18), (3.26) and (3.27) that
0<C≤min{T−pp−1+N2,T−pp−1+Np(q−1)q(p−1)}, |
which yields
\begin{equation} 0 \lt C\leq T^{\frac{-p}{p-1}+\frac{N}{2}}: = T^{\alpha(N)} \end{equation} | (3.28) |
and
\begin{equation} 0 \lt C\leq T^{\frac{-p}{p-1}+\frac{Np(q-1)}{q(p-1)}}: = T^{\beta(N)}. \end{equation} | (3.29) |
Notice that for N\in\{1, 2\} , one has
\alpha(N) = \left\{\begin{array}{lll} \frac{-p-1}{2(p-1)} &\mbox{if}& N = 1,\\ \frac{-1}{p-1} &\mbox{if}& N = 2, \end{array} \right. |
which shows that
\alpha(N) \lt 0,\quad N\in \{1,2\}. |
Hence, for N\in\{1, 2\} , passing to the limit as T\to \infty in (3.28), a contradiction follows ( 0 < C\leq 0 ). This proves part (ⅰ) of Theorem 1.1.
(ⅱ) Let N\geq 3 . In this case, one has
\alpha(N) \lt 0 \Longleftrightarrow p \lt 1+\frac{2}{N-2}. |
Hence, if p < 1+\frac{2}{N-2} , passing to the limit as T\to \infty in (3.28), a contradiction follows. Furthermore, one has
\beta(N) \lt 0 \Longleftrightarrow q \lt 1+\frac{1}{N-1}. |
Hence, if q < 1+\frac{1}{N-1} , passing to the limit as T\to \infty in (3.29), a contradiction follows. Therefore, we proved part (ⅱ) of Theorem 1.1.
(ⅲ) Let
\begin{equation} p \gt 1+\frac{2}{N-2}\quad\mbox{and}\quad q \gt 1+\frac{1}{N-1}. \end{equation} | (3.30) |
Let
\begin{equation} u(x) = \varepsilon (1+r^2)^{-\delta},\quad r = |x|,\quad x\in \mathbb{R}^N, \end{equation} | (3.31) |
where
\begin{equation} \max\left\{\frac{1}{p-1},\frac{2-q}{2(q-1)}\right\} \lt \delta \lt \frac{N-2}{2} \end{equation} | (3.32) |
and
\begin{equation} 0 \lt \varepsilon \lt \min\left\{1,\left(\frac{2\delta(N-2\delta-2)}{1+2^q\delta^q}\right)^{\frac{1}{\min\{p,q\}-1}}\right\}. \end{equation} | (3.33) |
Note that due to (3.30), the set of \delta satisfying (3.32) is nonempty. Let
w(x) = -\Delta u-|u|^p-|\nabla u|^q,\quad x\in \mathbb{R}^N. |
Elementary calculations yield
\begin{eqnarray} w(x)& = &2\delta \varepsilon \left[(1+r^2)^{-\delta-1}-2r^2(\delta+1)(1+r^2)^{-\delta-2}+(N-1)(1+r^2)^{-\delta-1}\right]\\ && -\varepsilon^p(1+r^2)^{-\delta p}-2^q\delta^q\varepsilon^qr^q(1+r^2)^{-(\delta+1)q}. \end{eqnarray} | (3.34) |
Hence, one obtains
\begin{eqnarray*} w(x) &\geq & 2\delta \varepsilon \left[(1+r^2)^{-\delta-1}-2(\delta+1)(1+r^2)^{-\delta-1}+(N-1)(1+r^2)^{-\delta-1}\right]\\ &&-\varepsilon^p(1+r^2)^{-\delta p}-2^q\delta^q\varepsilon^q(1+r^2)^{-(\delta+1)q+\frac{q}{2}}\\ & = & 2\delta \varepsilon (N-2\delta-2)(1+r^2)^{-\delta-1}-\varepsilon^p(1+r^2)^{-\delta p}-2^q\varepsilon^q\delta^q(1+r^2)^{-\left(\delta+\frac{1}{2}\right)q}. \end{eqnarray*} |
Next, using (3.32) and (3.33), one deduces that
w(x)\geq \varepsilon \left[2\delta (N-2\delta-2)-\varepsilon^{p-1}-2^q\varepsilon^{q-1}\delta^q\right](1+r^2)^{-\delta-1} \gt 0. |
Therefore, for any \delta and \varepsilon satisfying respectively (3.32) and (3.33), the function u defined by (3.31) is a stationary solution (then global solution) to (1.1) for some w > 0 . This proves part (ⅲ) of Theorem 1.1.
Proposition 4.1. Let N\geq 3 , \sigma < N and w\in I_\sigma^+ . Then, for sufficiently large T ,
\begin{equation} \int_Q w(x)\xi_T\,dx\,dt \geq C T^{\rho(N-\sigma)+1}. \end{equation} | (4.1) |
where \xi_T is defined by (2.1).
Proof. By (3.14) and (3.15), one has
\begin{equation} \int_Q w(x)\xi_T\,dx\,dt = CT\int_{0 \lt |x| \lt \sqrt{2}T^\rho} w(x) g\left(\frac{|x|^2}{T^{2\rho}}\right)^\ell\,dx. \end{equation} | (4.2) |
On the other hand, by definition of I_\sigma^+ , and using that g(s) = 1 , 0\leq s\leq 1 , for sufficiently large T , one obtains (since w, g\geq 0 )
\begin{eqnarray*} \int_{0 \lt |x| \lt \sqrt{2}T^\rho} w(x) g\left(\frac{|x|^2}{T^{2\rho}}\right)^\ell\,dx &\geq & \int_{0 \lt |x| \lt T^\rho} w(x) g\left(\frac{|x|^2}{T^{2\rho}}\right)^\ell\,dx\\ & = & \int_{0 \lt |x| \lt T^\rho} w(x)\,dx\\ &\geq & \int_{\frac{T^\rho}{2} \lt |x| \lt T^\rho} w(x)\,dx\\ &\geq & C \int_{\frac{T^\rho}{2} \lt |x| \lt T^\rho} |x|^{-\sigma}\,dx \\ & = & C T^{\rho(N-\sigma)}. \end{eqnarray*} |
Hence, using (4.2), (4.1) follows.
Now, we are ready to prove the new critical behavior for the problem (1.1) stated by Theorem 1.2.
Proof of Theorem 1.2. (ⅰ) Suppose that problem (1.1) admits a global weak solution u with (u, \nabla u)\in L^p_{loc}(Q)\times L^q_{loc}(Q) . By Propositions 3.2 and 4.1, for sufficiently large T , one obtains
\begin{equation} CT^{\rho(N-\sigma)+1}\leq \min\{A(\xi_T),B(\xi_T)\}, \end{equation} | (4.3) |
where \xi_T , A(\xi_T) and B(\xi_T) are defined respectively by (2.1), (3.19) and (3.20). Next, using (3.26) and (4.3) with \rho = \frac{1}{2} , one deduces that
\begin{equation} 0 \lt C\leq T^{\frac{-p}{p-1}+\frac{\sigma}{2}}. \end{equation} | (4.4) |
Observe that
\frac{-p}{p-1}+\frac{\sigma}{2} \lt 0 \Longleftrightarrow \sigma \lt \frac{2p}{p-1}. |
Hence, if \sigma < \frac{2p}{p-1} , passing to the limit as T\to \infty in (4.4), a contradiction follows. Furthermore, using (3.27) and (4.3) with \rho = \frac{p(q-1)}{q(p-1)} , one deduces that
\begin{equation} 0 \lt C\leq T^{\frac{p}{p-1}\left(\frac{\sigma(q-1)}{q}-1\right)}. \end{equation} | (4.5) |
Observe that
\frac{p}{p-1}\left(\frac{\sigma(q-1)}{q}-1\right) \lt 0\Longleftrightarrow \sigma \lt \frac{q}{q-1}. |
Hence, if \sigma < \frac{q}{q-1} , passing to the limit as T\to \infty in (4.5), a contradiction follows. Therefore, part (ⅰ) of Theorem 1.2 is proved.
(ⅱ) Let
\begin{equation} \sigma \geq \max\left\{\frac{2p}{p-1},\frac{q}{q-1}\right\}. \end{equation} | (4.6) |
Let u be the function defined by (3.31), where
\begin{equation} \frac{\sigma-2}{2} \lt \delta \lt \frac{N-2}{2} \end{equation} | (4.7) |
and \varepsilon satisfies (3.33). Notice that since \sigma < N , the set of \delta satisfying (4.7) is nonempty. Moreover, due to (4.6) and (4.7), \delta satisfies also (3.32). Hence, from the proof of part (iii) of Theorem 1.2, one deduces that
w(x) = -\Delta u-|u|^p-|\nabla u|^q \gt 0,\quad x\in \mathbb{R}^N. |
On the other hand, using (3.34) and (4.7), for |x| large enough, one obtains
w(x)\leq C (1+|x|^2)^{-\delta-1} \leq C |x|^{-2\delta-2}\leq C|x|^{-\delta}, |
which proves that w\in I_\sigma^- . This proves part (ⅱ) of Theorem 1.2.
We investigated the large-time behavior of solutions to the nonlinear damped wave equation (1.1). In the case when N\in\{1, 2\} , we proved that for all p > 1 , problem (1.1) admits no global weak solution (in the sense of Definition 1). Notice that from [5], the same result holds for the problem without gradient term, namely problem (1.5). This shows that in the case N\in\{1, 2\} , the nonlinearity |\nabla u|^q has no influence on the critical behavior of problem (1.5). In the case when N\geq 3 , we proved that, if 1 < p < 1+\frac{2}{N-2} or 1 < q < 1+\frac{1}{N-1} , then problem (1.1) admits no global weak solution, while if p > 1+\frac{2}{N-2} and q > 1+\frac{1}{N-1} , global solutions exist for some w > 0 . This shows that in this case, the Fujita critical exponent for the problem (1.1) is given by
p_c(N,q) = \left\{\begin{array}{lll} 1+\frac{2}{N-2} &\mbox{if}& q \gt 1+\frac{1}{N-1},\\ \infty &\mbox{if}& q \lt 1+\frac{1}{N-1}. \end{array} \right. |
From this result, one observes two facts. First, in the range q > 1+\frac{1}{N-1} , from [5], the critical exponent p_c(N, q) is also equal to the critical exponent for the problem without gradient term, which means that in this range of q , the nonlinearity |\nabla u|^q has no influence on the critical behavior of problem (1.5). Secondly, one observes that the gradient term induces an interesting phenomenon of discontinuity of the Fujita critical exponent p_c(N, q) jumping from 1+\frac{2}{N-2} to \infty as q reaches the value 1+\frac {1}{N-1} from above. In the same case N\geq 3 , we determined also the second critical exponent for the problem (1.1) in the sense of Lee and Ni [10], when p > 1+\frac{2}{N-2} and q > 1+\frac{1}{N-1} . Namely, we proved that in this case, if \sigma < \max\left\{\frac{2p}{p-1}, \frac{q}{q-1}\right\} and w\in I_\sigma^+ , then there is no global weak solution, while if \max\left\{\frac{2p}{p-1}, \frac{q}{q-1}\right\}\leq \sigma < N , global solutions exist for some w\in I_\sigma^- . This shows that the second critical exponent for the problem (1.1) in the sense of Lee and Ni is given by
\sigma^* = \max\left\{\frac{2p}{p-1},\frac{q}{q-1}\right\}. |
We end this section with the following open questions:
(Q1). Find the first and second critical exponents for the system of damped wave equations with mixed nonlinearities
\begin{eqnarray*} \left\{\begin{array}{lllll} u_{tt}-\Delta u +u_t & = & |v|^{p_1}+|\nabla v|^{q_1} +w_1(x) &\mbox{in}& (0,\infty)\times \mathbb{R}^N,\\ v_{tt}-\Delta v +v_t& = & |u|^{p_2}+|\nabla u|^{q_2} +w_2(x) &\mbox{in}& (0,\infty)\times \mathbb{R}^N, \end{array} \right. \end{eqnarray*} |
where p_i, q_i > 1 , w_i\in L^1_{loc}(\mathbb{R}^N) , w_i\geq 0 and w_i\not\equiv 0 , i = 1, 2 .
(Q2). Find the Fujita critical exponent for the problem (1.1) with w\equiv 0 .
(Q3) Find the Fujita critical exponent for the problem
u_{tt}-\Delta u +u_t = \frac{1}{\Gamma(1-\alpha)}\int_0^t (t-s)^{-\alpha}|u(s,x)|^{p}\,ds +|\nabla u|^{q} +w(x),\quad t \gt 0,\,x\in \mathbb{R}^N, |
where N\geq 1 , p, q > 1 , 0 < \alpha < 1 and w\in L^1_{loc}(\mathbb{R}^N) , w\geq 0 . Notice that in the limit case \alpha \to 1^- , the above equation reduces to (1.1).
(Q4) Same question as above for the problems
u_{tt}-\Delta u +u_t = |u|^p+\frac{1}{\Gamma(1-\alpha)}\int_0^t (t-s)^{-\alpha}|\nabla u(s,x)|^{q}\,ds+w(x) |
and
u_{tt}-\Delta u +u_t = \frac{1}{\Gamma(1-\alpha)}\int_0^t (t-s)^{-\alpha}|u(s,x)|^{p}\,ds +\frac{1}{\Gamma(1-\beta)}\int_0^t(t-s)^{-\beta}|\nabla u(s,x)|^{q} +w(x), |
where 0 < \alpha, \beta < 1 .
The author is supported by Researchers Supporting Project RSP-2020/4, King Saud University, Saudi Arabia, Riyadh.
The author declares no conflict of interest.
[1] | World Health Organization, Coronavirus disease 2019. cited March 15, 2020. Available from: https://www.who.int/health-topics/coronavirus. |
[2] | Editorial, The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health- The latest 2019 novel coronavirus outbreak in Wuhan, China, Int. J. Infect. Dis., 91 (2020), 264-266. |
[3] | Worldometer. Available from: https://www.worldometers.info/coronavirus. |
[4] | World Health Organization, Population-based age-stratified seroepidemiological investigation protocol for covid-19 virus infection, 2020. |
[5] | N. M. Ferguson, D. Laydon, G. Nedjati-Gilani, N. Imai, K. Ainslie, M. Baguelin, et al., Impact of non-pharmaceutical interventions (npis) to reduce covid-19 mortality and healthcare demand, London: Imperial College COVID-19 Response Team, 10 (2020), 10.25561/77482. |
[6] | B. Tang, N. L. Bragazzi, Q. Li, S. Tang, Y. Xiao, J. Wu, An updated estimation of the risk of transmission of the novel coronavirus (2019-ncov), Infect. Dis. Model., 5 (2020), 248-255. |
[7] | B. J. Quilty, S. Clifford, S. Flasche, R. M. Eggo, Effectiveness of airport screening at detecting travellers infected with novel coronavirus (2019-nCoV). Euro. Surveil., 25 (2020), 200080. |
[8] | M. Shen, Z. Peng, Y. Xiao, L Zhang, Modelling the epidemic trend of the 2019 novel coronavirus outbreak in china, bioRxiv, 2020. |
[9] |
A. J. Kucharski, T. W. Russell, C. Diamond, Y. Liu, J. Edmunds, S. Funk, et al., Early dynamics of transmission and control of COVID-19: a mathematical modelling study, Lancet Infect. Dis., 20 (2020), 553-558. doi: 10.1016/S1473-3099(20)30144-4
![]() |
[10] |
J. Yuan, M. Li, G. Lv, Z. K. Lu, Monitoring transmissibility and mortality of COVID-19 in Europe, Int. J. Infec. Dis., 95 (2020), 311-325. doi: 10.1016/j.ijid.2020.03.050
![]() |
[11] | G. Giordano, F. Blanchini, R. Bruno, P. Colaneri, A. D. Filippo, A. D. Matteo, et al., Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy, Nat. Med., (2020), 1-6. |
[12] |
Q. Lin, S. Zhao, D. Gao, Y. Lou, S. Yang, S. S. Musa, et al., A conceptual model for the coronavirus disease 2019 (COVID-19) outbreak in Wuhan, China with individual reaction and governmental action, Int. J. Infect. Dis., 93 (2020), 211-216. doi: 10.1016/j.ijid.2020.02.058
![]() |
[13] | T. Chen, J. Rui, Q. Wang, Z. Zhao, J. Cui, L. Yin, A mathematical model for simulating the phase-based transmissibility of a novel coronavirus, Infect. Dis. Poverty, 9 (2020) 1-8. |
[14] | J. T. Wu, K. Leung, G. M. Leung, Nowcasting and forecasting the potential domestic and international spread of the 2019-ncov outbreak originating in wuhan, china: a modelling study, Lancet, 395 (2020), 689-697. |
[15] | B. Tang, X. Wang, Q. Li, N. L. Bragazzi, S. Tang, Y. Xiao, et al., Estimation of the transmission risk of the 2019-nCov and its implication for public health interventions, J. Clin. Med., 9 (2020), 462. |
[16] | J. M. Read, J. R. Bridgen, D. A. Cummings, A. Ho, C. P. Jewell, Novel coronavirus 2019-nCov: early estimation of epidemiological parameters and epidemic predictions, medRxiv, 2020. |
[17] |
S. Zhao, Q. Lin, J. Ran, S. S. Musa, G. Yang, W. Wang, et al., Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak, Int. J. Infect. Dis., 92 (2020), 214-217. doi: 10.1016/j.ijid.2020.01.050
![]() |
[18] |
J. Chen, Pathogenicity and transmissibility of 2019-nCoV - a quick overview and comparison with other emerging viruses, Microb. infect., 22 (2020), 69-71. doi: 10.1016/j.micinf.2020.01.004
![]() |
[19] | R. Singh, R. Adhikari, Age-structured impact of social distancing on the covid-19 epidemic in India, arXiv preprint, (2020), arXiv: 2003.12055. |
[20] | U, Avila-Ponce de Leon, A, G. C. Perez, E, Avila-Vales, An SEIARD epidemic model for COVID-19 in Mexico: mathematical analysis and state-level forecast, Chaos Solitons Fractals, 140 (2020), 110165. |
[21] | J. H. Rojas, M. Paredes, M. Banerjee, O. Akman, A. Mubayi, Mathematical Modeling & the Transmission Dynamics of SARS-CoV-2 in Cali, Colombia: Implications to a 2020 Outbreak & public health preparedness, medRxiv, (2020), https://doi.org/10.1101/2020.05.06.20093526. |
[22] | E. Shim, G. Chowell, Regional variability in time-varying transmission potential of COVID-19 in South Korea, medRxiv, (2020), https://doi.org/10.1101/2020.07.21.20158923. |
[23] | A. Srivastava, G. Chowell, Understanding Spatial Heterogeneity of COVID-19 Pandemic Using Shape Analysis of Growth Rate Curves, medRxiv, (2020), https://doi.org/10.1101/2020.05.25.20112433. |
[24] | Y. Belgaid, M. Helal, E. Venturino, Analysis of a model for Coronavirus spread, Mathematics, 8(5), (2020), 820. |
[25] | J. Dolbeault, G. Turinici, Heterogeneous social interactions and the Covid-19 lockdown outcome in a multi-group SEIR model, Math. Model. Nat. Phenom., 15 (2020), 36. |
[26] | M. Kochanczyk, F. Grabowski, T. Lipniacki, Dinamics of Covid-19 pandemic at constant and time-dependent contact rates, Math. Model. Nat. Phenom., 15 (2020), 28. |
[27] | S. G. Krantz, P. Polyakov, A. S. R. S. Rao, True epidemic growth construction through harmonic analysis, J. Theor. Biol., 494, (2020), 110243. |
[28] | S. Sinha, Epidemiological dynamics of the COVID-19 pandemic in India: an interim assessment, Stat. Appl., 18 (2020), 333-350. |
[29] | H. R. Thieme, Mathematics in Population Biology, Princeton University Press, Princeton, 2003. |
[30] | A. V. Emmanuelle, Lifting the COVID-19 lockdown: different scenarios for France, Math. Model. Nat. Phenom., (2020), In press. |
[31] | L. D. Domenico, G. Pullano, C. E. Sabbatini, P. Y. Boelle, V. Colizza, Expected impact of reopening schools after lockdown on COVID-19 epidemic in Ile-de-France, medRxiv, (2020), https://doi.org/10.1101/2020.05.08.20095521. |
[32] |
P. van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29-48. doi: 10.1016/S0025-5564(02)00108-6
![]() |
[33] | O. Diekmann, J. A. P. Heesterbeek, J. A. J. Metz, On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations, J. Math. Biol., 28 (1990), 365-382. |
[34] | F. Brauer, C. Castillo-Chavez, Z. Feng, Mathematical Models in Epidemiology, Springer, New York, 2019. |
[35] | M. Martcheva, An Introduction to Mathematical Epidemiology, Springer, New York, 2015. |
[36] |
V. Andreasen, The final size of an epidemic and its relation to the basic reproduction number, Bull. Math.Biol., 73 (2011), 2305-2321. doi: 10.1007/s11538-010-9623-3
![]() |
[37] | R. J. Freund, W. J. Wilson, D. L. Mohr, Statistical Methods, Elsevier, Canada, 2010. |
[38] | C. T. Kelley, Iterative Methods for Optimization, SIAM, Philadelphia, USA, 1999. |
[39] | D. Caccavo, Chinese and Italian COVID-19 outbreaks can be correctly described by a modified SIRD model, medRxiv, (2020), https://doi.org/10.1101/2020.03.19.20039388. |
[40] | N. M. Duggan, S. M. Ludy, B. C. Shannon, A. T. Reisner, S. R. Wilcox, Is novel coronavirus 2019 reinfection possible? Interpreting dynamic SARS-CoV-2 test results through a case report, Am. J. Emerg. Med., (2020), in press. |
[41] |
V. Capasso, S. Gabriella, A generalization of the Kermack-McKendrick deterministic epidemic model, Math. Biosci., 42 (1978), 43-61. doi: 10.1016/0025-5564(78)90006-8
![]() |
1. | Sen Ming, Jiayi Du, Bo Du, Existence of Blow-Up Solution to the Cauchy Problem of Inhomogeneous Damped Wave Equation, 2025, 17, 2073-8994, 1009, 10.3390/sym17071009 |