This study focused on a coupled nonlinear wave equation system featuring fractional damping and polynomial source terms within a bounded domain. We demonstrated that, under specific conditions, the solutions exhibited blow-up in a finite time-frame.
Citation: Abdelhadi Safsaf, Suleman Alfalqi, Ahmed Bchatnia, Abderrahmane Beniani. Blow-up dynamics in nonlinear coupled wave equations with fractional damping and external source[J]. Electronic Research Archive, 2024, 32(10): 5738-5751. doi: 10.3934/era.2024265
[1] | Yaning Li, Mengjun Wang . Well-posedness and blow-up results for a time-space fractional diffusion-wave equation. Electronic Research Archive, 2024, 32(5): 3522-3542. doi: 10.3934/era.2024162 |
[2] | Xiaoqiang Dai, Chao Yang, Shaobin Huang, Tao Yu, Yuanran Zhu . Finite time blow-up for a wave equation with dynamic boundary condition at critical and high energy levels in control systems. Electronic Research Archive, 2020, 28(1): 91-102. doi: 10.3934/era.2020006 |
[3] | Yuchen Zhu . Blow-up of solutions for a time fractional biharmonic equation with exponentional nonlinear memory. Electronic Research Archive, 2024, 32(11): 5988-6007. doi: 10.3934/era.2024278 |
[4] | Gongwei Liu . The existence, general decay and blow-up for a plate equation with nonlinear damping and a logarithmic source term. Electronic Research Archive, 2020, 28(1): 263-289. doi: 10.3934/era.2020016 |
[5] | Vo Van Au, Jagdev Singh, Anh Tuan Nguyen . Well-posedness results and blow-up for a semi-linear time fractional diffusion equation with variable coefficients. Electronic Research Archive, 2021, 29(6): 3581-3607. doi: 10.3934/era.2021052 |
[6] | Shuting Chang, Yaojun Ye . Upper and lower bounds for the blow-up time of a fourth-order parabolic equation with exponential nonlinearity. Electronic Research Archive, 2024, 32(11): 6225-6234. doi: 10.3934/era.2024289 |
[7] | Jorge A. Esquivel-Avila . Blow-up in damped abstract nonlinear equations. Electronic Research Archive, 2020, 28(1): 347-367. doi: 10.3934/era.2020020 |
[8] | Hui Jian, Min Gong, Meixia Cai . Global existence, blow-up and mass concentration for the inhomogeneous nonlinear Schrödinger equation with inverse-square potential. Electronic Research Archive, 2023, 31(12): 7427-7451. doi: 10.3934/era.2023375 |
[9] | Qianqian Zhu, Yaojun Ye, Shuting Chang . Blow-up upper and lower bounds for solutions of a class of higher order nonlinear pseudo-parabolic equations. Electronic Research Archive, 2024, 32(2): 945-961. doi: 10.3934/era.2024046 |
[10] | Xiao Su, Hongwei Zhang . On the global existence and blow-up for the double dispersion equation with exponential term. Electronic Research Archive, 2023, 31(1): 467-491. doi: 10.3934/era.2023023 |
This study focused on a coupled nonlinear wave equation system featuring fractional damping and polynomial source terms within a bounded domain. We demonstrated that, under specific conditions, the solutions exhibited blow-up in a finite time-frame.
In this study, we examine the following problem:
{φtt+∂α,ηtφ=Δφ+|φ+ψ|p−1(φ+ψ)+|φ|p−32φ|ψ|p+12,x∈Ω,t>0,ψtt+∂β,ηtψ=Δψ+|φ+ψ|p−1(φ+ψ)+|ψ|p+12ψ|φ|p−32,x∈Ω,t>0,φ(x,t)=ψ(x,t)=0,x∈∂Ω,t>0,φ(x,0)=φ0(x),φt(x,0)=φ1(x),x∈Ω,ψ(x,0)=ψ0(x),ψt(x,0)=ψ1(x),x∈Ω, | (1.1) |
where Ω is a bounded domain with a smooth boundary ∂Ω in Rn. The constants p and α satisfy
{2<p<∞,n=1,2,2<p≤2nn−2,n≥3, |
and 0<α,β<1. The functions φ0(x), φ1(x), ψ0(x), and ψ1(x) are given.
The symbol ∂α,ηt denotes the Caputo fractional derivative of order α with respect to the time variable t, where 0<α<1. It is defined as:
∂α,ηtW(t)=1Γ(1−α)∫t0(t−s)−αe−η(t−s)dWds(s)ds,η≥0, |
which was introduced by Michele Caputo in [1].
Problems of this type arise in material science and physics, particularly in the study of wave propagation in viscoelastic materials, fluid dynamics, and structural mechanics. These models are used to describe the behavior of materials with memory effects, where the damping terms are often nonlocal or fractional in nature. Such systems have been widely studied in the context of their stability, energy decay, and response to external forces (see, for example, [2,3,4]). The fractional damping terms in the model are used to capture more complex dissipative effects that arise in these physical systems, leading to more accurate descriptions of the material behavior under stress and deformation. Agre and Rammaha [2] established both existence and blow-up results for systems of the form (1.1), specifically when the initial energy is sufficiently small and positive. More recently, Han and Wang [4] extended these results to address systems involving a viscoelastic term, further generalizing the analysis conducted in [2].
In the case of a single wave equation given by
φtt−∂α,ηtφ=Δφ+a|φ|p−1φ, |
Kirane and Tatar [5] demonstrated that solutions exhibit exponential growth when the initial data is sufficiently large. Matignon et al. [6] explored the case where a=0, achieving results on well-posedness and asymptotic stability by reformulating the problem into a standard form. Tatar [7] employed Fourier transform techniques and the Hardy-Littlewood-Sobolev inequality to establish the finite-time blow-up of solutions. Additionally, [8] showed that finite-time blow-up can occur independently of the time variable T associated with the initial data.
Without internal fractional damping, a polynomial source can lead to finite-time blow-up for solutions with negative initial energy [9,10,11]. In [12], Georgiev and Todorova demonstrated that the solution remains global if p≤m, whereas solutions blow up in finite time if this condition is not met. Ball [9] established finite-time blow-up in certain scenarios using continuation theorems for ordinary differential equations within Banach spaces. Additional theorems in different contexts are discussed in [13].
Motivated by the aforementioned research, we consider a coupled system with fractional derivatives. This choice is driven by the growing interest in recent years within the scientific community to explore the complex dynamics and practical applications of wave equations. The extension to fractional derivatives is significant because they offer more accurate models for systems exhibiting memory effects and complex dissipation properties, which are commonly encountered in practical scenarios, such as material microstructure analysis. Introducing fractional derivatives into the coupled system is novel in the context of wave equations with external sources, allowing for a more realistic modeling of dissipative effects in various materials and structures. Our results provide new insights into the blow-up behavior of systems with fractional derivatives.
The rest of the paper is organized as follows: In Section 2, we present some preliminary results that will be useful throughout the paper. In Section 3, we provide the main result concerning the blow-up behavior of the solutions of the system and offer a detailed proof.
For a solution (φ,ψ) to the problem (1.1), the energy functional is given by
E(t)=12∫Ω(φ2t+ψ2t)dx+12∫Ω(|∇φ|2+|∇ψ|2)dx−1p+1∫Ω|φ+ψ|p+1dx−2p+1∫Ω|φ|p+12|ψ|p+12dx. | (2.1) |
In this definition, the classical energy functional E(t) is described in terms of kinetic and potential energy, along with additional interaction terms.
To analyze the behavior of this energy functional, we state and prove the following lemma.
Lemma 2.1. Assume (φ,ψ) is a regular solution to problem (1.1). The energy functional, as defined by (2.1), satisfies
dE(t)dt=−1Γ(1−α)∫Ωφt∫t0(t−s)−αe−η(t−s)φt(s)dsdx−1Γ(1−β)∫Ωψt∫t0(t−s)−βe−η(t−s)ψt(s)dsdx, | (2.2) |
and we deduce E(t)≤E(0),∀t≥0.
Proof. By multiplying the first equation in (1.1) by φt and the second equation by ψt, then integrating over Ω and applying integration by parts, we derive Eq (2.2). Determining the sign of this functional is not immediately apparent; however, Nohel and Shea demonstrated in [14] that a real function a(t,x)∈L1loc(R+;L∞(Ω)) is of positive type if it satisfies the following condition:
∫T0∫Ωh(t)∫t0a(t−s)h(s)dsdσdt≥0, | (2.3) |
for all h∈C(R+;H1(Ω)) and for every T>0.
Next, by substituting t with z and integrating Eq (2.2) from 0 to t, we obtain:
E(t)−E(0)=−1Γ(1−α)∫t0∫Ωφt(z)∫z0(z−s)−αe−η(z−s)φt(s)dsdxdz−1Γ(1−β)∫t0∫Ωψt(z)∫z0(z−s)−βe−η(z−s)ψt(s)dsdxdz. | (2.4) |
Using condition (2.3), we directly deduce the classical energy inequality:
∀t≥0,E(0)≥E(t). |
The following lemma will be beneficial for what follows.
Lemma 2.2. There exist positive constants c1 and c2 such that
c1(|x|p+1+|y|p+1)≤1p+1|x+y|p+1+2p+1|xy|p+12≤c2(|x|p+1+|y|p+1), |
for all x,y∈R.
Proof. To start, using the inequality 2ab≤a2+b2, we obtain:
|x|p+12|y|p+12≤12(|x|p+1+|y|p+1). |
Next, we use the fact that h(x):=|x|p+1 is convex over R+, and so
|x+y|p+1≤2p(|x|p+1+|y|p+1). |
By combining the results above, we obtain:
1p+1|x+y|p+1+2p+1|xy|p+12≤2p+1p+1(|x|p+1+|y|p+1). | (2.5) |
On the other hand, we can suppose that |y|≤|x|, and so:
|x|p+1=|x+y−y|p+1≤2p(|x+y|p+1+|y|p+1)≤2p(|x+y|p+1+|xy|p+12). |
Therefore,
|x|p+1+|y|p+1≤2p+1(|x+y|p+1+|xy|p+12). | (2.6) |
Finally, by combining (2.5) and (2.6), we have demonstrated that there exist positive constants c1 and c2 such that:
c1(|x|p+1+|y|p+1)≤1p+1|x+y|p+1+2p+1|xy|p+12≤c2(|x|p+1+|y|p+1). |
Denote by T∗ the maximal interval of time for which a solution is defined. When T∗≠+∞, the solution is said to blow up in finite time. Blow-up means that a singularity forms that prevents the solution from being extended beyond T∗. If we examine the ordinary differential equation related to the one-dimensional (1D) semi-linear wave equation: utt=|u|p−1u, it is known that its solutions exhibit blow-up in finite time. This observation supports our assertion that the solution of our system also experiences blow-up in finite time.
This phenomenon can have two interpretations. If it arises in the modeling of a physical phenomenon for which it has no meaning, it signifies that the model is oversimplified; this is linked with the instability of numerical schemes that one could employ. It can also be of physical relevance (finite time collapse of a star, shocks, etc.).
In this section, we are ready to state and prove our result. To begin, we present a local existence theorem, which can be established using [12] and easily adapted to the case of a system of equations:
Theorem 3.1. Let (φ0,φ1) and (ψ0,ψ1)∈H10(Ω)×L2(Ω) be given. Suppose p satisfies
{2<p<∞,forn=1,2,2<p≤2nn−2,forn≥3. |
Then, the problem (1.1) has a unique local solution
(φ,ψ)∈C([0,Tm);H10(Ω)),(φt,ψt)∈C([0,Tm);L2(Ω))∩L2(Ω), |
for some Tm>0.
We now present and prove the main result of our study.
Theorem 3.2. Let (φ,ψ) be the solution of the system (1.1), and assume that the initial data φ0 and ψ0 satisfies
∫Ω(φ20+ψ20)dx≠0. |
Then, for any T>0, there exist T∗≤T and sufficiently large initial data such that (φ,ψ) blows up at T∗.
Proof. To ensure clarity, we will break the proof into several steps.
Step 1. We introduce the functional H(t) defined by:
H(t)=−∫t0E(s)ds+(dt+l)∫Ω(φ20+ψ20)dx, |
where d and l are positive constants whose values will be specified later.
Since the proof relies on a contradiction argument, we will associate the function H with the function Ψ defined in Step 2. We assume that the solution exists up to time T and derive an inequality of the form Ψ′(t)≥CΨδ(t) with δ>1. Integrating this inequality demonstrates that, under certain conditions, the solution cannot exist for all time up to T.
Differentiating H(t) in relation to t, one gets:
H′(t)=−E(t)+d∫Ω(φ20+ψ20)dx≥d∫Ω(φ20+ψ20)dx−E(0). | (3.1) |
Selecting d ensures that H′(0)>0. It is sufficient to select d such that:
H′(0)=−E(0)+d∫Ω(φ20+ψ20)dx>0. |
By virtue of (2.4) and (3.1), it follows that:
H′(0)−H′(t)=E(t)−E(0)=−1Γ(1−α)∫t0∫Ωφt(z)∫z0(z−s)−αe−η(z−s)φt(s)dsdxdz−1Γ(1−β)∫t0∫Ωψt(z)∫z0(z−s)−βe−η(z−s)ψt(s)dsdxdz<0. | (3.2) |
Step 2. We define the functional
Ψ(t)=H1−γ(t)+ε2(∫Ω(φ2+ψ2)dx−∫Ω(φ20+ψ20)dx), | (3.3) |
where ε>0 and 0<γ<p−12(p+1).
Initially, we have
Ψ(0)=H1−γ(0)=(l∫Ω(φ20+ψ20)dx)1−γ. |
Differentiating (3.3) in relation to t gives
Ψ′(t)=(1−γ)H−γ(t)H′(t)+ε∫Ω(φφt+ψψt)dx. | (3.4) |
Further differentiating (3.4) and then integrating yields:
Ψ′(t)=(1−γ)H−γ(t)H′(t)+ε∫Ω(φ0φ1+ψ0ψ1)dx+ε∫t0∫Ω(φ2t+ψ2t)dxds+ε∫t0∫Ω(φφtt+ψψtt)dxds. | (3.5) |
To assess the last term on the righthand side of (3.5), we multiply the first and second equations of problem (1.1) by φ and ψ, respectively, and integrate over Ω×(0,t). This yields:
∫t0∫Ω(φφtt+ψψtt)dxds=−∫t0∫Ω(|∇φ|2+|∇ψ|2)dxds+1p+1∫t0∫Ω|φ+ψ|p+1dxds+1p+1∫t0∫Ω|φ|p+12|ψ|p+12dxds−1Γ(1−α)∫t0∫Ωφ∫s0(s−z)−αe−η(s−z)φt(z)dzdxds−1Γ(1−β)∫t0∫Ωψ∫s0(s−z)−βe−η(s−z)ψt(z)dzdxds. | (3.6) |
Inserting the relation (3.6) into (3.5), one gets
Ψ′(t)=(1−γ)H−γ(t)H′(t)+ε∫Ω(φ0φ1+ψ0ψ1)dx+ε∫t0∫Ω(φ2t+ψ2t)dxds−ε∫t0∫Ω(|∇φ|2+|∇ψ|2)dxds+εp+1∫t0∫Ω|φ+ψ|p+1dxds+2εp+1∫t0∫Ω|φ|p+12|ψ|p+12dxds−εΓ(1−α)∫t0∫Ωφ∫s0(s−z)−αe−η(s−z)φt(z)dzdxds−εΓ(1−β)∫t0∫Ωψ∫s0(s−z)−βe−η(s−z)ψt(z)dzdxds. | (3.7) |
We introduce the extension operators over the entire domain in the following manner:
Lω(s)={ωif s∈[0,t],0if s∈R∖[0,t], |
and
Lκβ(s)={κβ(s)if s>0,0if s<0, |
where the kernel is defined by
κγ(t)=tγ−1Γ(1−γ),0<γ<1. |
We start with the following equation:
1Γ(1−α)∫t0φ(s)∫s0(s−z)−αφt(z)dzds=∫RLφ(s)∫RLκα(s−z)(Lφt)(z)dzds. | (3.8) |
Using Parseval's theorem, we rewrite the left-hand side of (3.8) as:
∫RLφ(s)∫RLκα(s−z)(Lφt)(z)dzds=∫RF(Lφ)(σ)¯F(Lκα⋆Lφt)(σ)dσ, |
where F(f) represents the standard Fourier transform of f.
Using the convolution property, we obtain:
κγ+η(t)=κγ(t)⋆κη(t),0<γ,η<1. |
Applying the Cauchy-Schwarz and Young inequalities, we obtain:
∫RLφ(s)∫RLκα(s−z)(Lφt)(z)dzds≤(∫R|F(Lκα2)F(Lφ)|2dσ)12(∫R|F(Lκα2)F(Lφt)|2dσ)12≤δ∫R|F(Lκα2)F(Lφt)|2dσ+14δ∫R|F(Lκα2)F(Lφ)|2dσ, | (3.9) |
where δ>0. From (3.9), as stated in [15], we infer:
∫RLφ(s)∫RLκα(s−z)(Lφt)(z)dzds≤1cos(απ2)[δ∫RLφt(s)(Lκα⋆Lφt)(s)ds+14δ∫RLφ(s)(Lκα⋆Lφ)(s)ds]. | (3.10) |
Similarly, we have:
∫RLψ(s)∫RLκβ(s−z)(Lψt)(z)dzds≤1cos(βπ2)[δ∫RLψt(s)(Lκβ⋆Lψt)(s)ds+14δ∫RLψ(s)(Lκβ⋆Lψ)(s)ds]. | (3.11) |
Inserting the estimates (3.10) and (3.11) into (3.7), we deduce
Ψ′(t)≥(1−γ)H−γ(t)H′(t)+ε∫Ω(φ0φ1+ψ0ψ1)dx+ε∫t0∫Ω(φ2t+ψ2t)dxds+ε∫t0∫Ω(|∇φ|2+|∇ψ|2)dxds+εp+1∫t0∫Ω|φ+ψ|p+1dxds+2εp+1∫t0∫Ω|φ|p+12|ψ|p+12dxds−εcos(απ2)[δ∫Ω∫RLφt(s)(Lκα⋆Lφt)(s)dsdx+14δ∫Ω∫RLφ(s)(Lκα⋆Lφ)(s)dsdx]−εcos(βπ2)[δ∫Ω∫RLψt(s)(Lκβ⋆Lψt)(s)dsdx+14δ∫Ω∫RLψ(s)(Lκβ⋆Lψ)(s)dsdx]. | (3.12) |
If we define k=min(cos(απ2),cos(βπ2)), from (3.2) we see that
−εδcos(απ2)∫Ω∫RLφt(s)(Lκα⋆Lφt)(s)dsdx−εδcos(βπ2)∫Ω∫RLψt(s)(Lκβ⋆Lψt)(s)dsdx≥εδk[H′(0)−H′(t)]. | (3.13) |
Therefore, (3.12) and (3.13) imply:
Ψ′(t)≥(1−γ)H−γ(t)H′(t)+ε∫Ω(φ0φ1+ψ0ψ1)dx+ε∫t0∫Ω(φ2t+ψ2t)dxds+ε∫t0∫Ω(|∇φ|2+|∇ψ|2)dxds+εp+1∫t0∫Ω|φ+ψ|p+1dxds+2εp+1∫t0∫Ω|φ|p+12|ψ|p+12dxds+εδk[H′(0)−H′(t)]−ε4δk∫Ω∫RLφ(s)(Lκα⋆Lφ)(s)dsdx−ε4δk∫Ω∫RLψ(s)(Lκβ⋆Lψ)(s)dsdx. | (3.14) |
We now estimate the last two terms in (3.14). Using the Cauchy-Schwarz inequality, we have:
∫Ω∫RLφ(s)(Lκα⋆Lφ)(s)dsdx≤∫Ω(∫R|Lφ|2ds)12(∫R|Lκα⋆Lφ|2(s)ds)12dx. |
Then, by Young's inequality (see [8]), we obtain:
∫Ω∫RLφ(s)(Lκα⋆Lφ)(s)dsdx≤1Γ(α)(∫t0sα−1ds)(∫Ω∫R|Lφ|2dsdx)≤tαΓ(1+α)∫Ω∫R|Lφ|2dsdx. | (3.15) |
Similarly,
∫Ω∫t0Lφ(s)(Lκβ⋆Lψ)(s)dsdx≤1Γ(β)(∫t0sβ−1ds)(∫Ω∫R|Lψ|2dsdx)≤tβΓ(1+β)∫Ω∫R|Lψ|2dsdx. | (3.16) |
Substituting relations (3.15) and (3.16) into (3.14), one gets:
Ψ′(t)≥[(1−γ)H−γ(t)−εδk]H′(t)+εδkH′(0)+ε∫Ω(φ0φ1+ψ0ψ1)dx+ε∫t0∫Ω(φ2t+ψ2t)dxds+ε∫t0∫Ω(|∇φ|2+|∇ψ|2)dxds+εp+1∫t0∫Ω|φ+ψ|p+1dxds+2εp+1∫t0∫Ω|φ|p+12|ψ|p+12dxds−εtα4δkΓ(α+1)∫Ω∫t0|φ|2dsdx−εtβ4δkΓ(β+1)∫Ω∫t0|ψ|2dsdx. |
Choosing δ=MkH−γ(t) for some M>0, we obtain:
Ψ′(t)≥[(1−γ)−Mε]H−γ(t)H′(t)+εMH−γ(t)H′(0)+ε∫Ω(φ0φ1+ψ0ψ1)dx+ε∫t0∫Ω(φ2t+ψ2t)dxds+ε∫t0∫Ω(|∇φ|2+|∇ψ|2)dxds+εp+1∫t0∫Ω|φ+ψ|p+1dxds+2εp+1∫t0∫Ω|φ|p+12|ψ|p+12dxds−εtαHγ(t)4Mk2Γ(α+1)∫Ω∫t0|φ|2dsdx−εtβHγ(t)4Mk2Γ(β+1)∫Ω∫t0|ψ|2dsdx. | (3.17) |
To bound the last two terms in (3.17), we apply Hölder's inequality and Lemma 2.2:
∫t0∫Ω|φ|2dsdx≤|Ω|p−1p+1tp−1p+1(∫t0∫Ω|φ|p+1dxds)2p+1≤|Ω|p−1p+1Tp−1p+1(∫t0∫Ω(|φ|p+1+|ψ|p+1)dxds)2p+1≤C1(∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds)2p+1, | (3.18) |
where C1=|Ω|p−1p+1Tp−1p+1. Similarly, we have:
∫t0∫Ω|ψ|2dsdx≤C1(∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds)2p+1. | (3.19) |
By recalling the definition of Hγ(t), we obtain:
I=Hγ(t)(∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds)2p+1≤[1p+1∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds+(dt+l)∫Ω(φ20+ψ20)dx]γ×(∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds)2p+1≤[1(p+1)γ(∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds)γ+(dt+l)γ(∫Ω(φ20+ψ20)dx)γ]×(∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds)2p+1≤c1(∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds)γ+2p+1+(dt+l)γ(∫Ω(φ20+ψ20)dx)γ×(∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds)2p+1, | (3.20) |
where c1=1(p+1)γ. We choose γ<1−2p+1, and one gets:
I≤c1[1+∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds]+(dT+l)γ(∫Ω(φ20+ψ20)dx)γ×[1+∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds]. | (3.21) |
Substituting (3.18), (3.19), (3.20), and (3.21) into (3.17) yields:
Ψ′(t)≥[(1−γ)−Mε]H−γ(t)H′(t)+εMH−γ(t)H′(0)+ε∫Ω(φ0φ1+ψ0ψ1)dx+ε∫t0∫Ω(φ2t+ψ2t)dxds+ε∫t0∫Ω(|∇φ|2+|∇ψ|2)dxds+ε∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds−εBMk2[1+∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds], | (3.22) |
where B=14[c1TαΓ(α+1)+c1TβΓ(β+1)][c3+(dT+l)γ(∫Ω(φ20+ψ20)dx)γ].
We select ε>0 such that ε≤1−γM. Then, from (3.22), we deduce:
Ψ′(t)≥ε∫Ω(φ0φ1+ψ0ψ1)dx+ε∫t0∫Ω(φ2t+ψ2t)dxds+ε(1−BMK2)∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds−εBMk2. |
Moreover, by adding and subtracting λεH(t) with 0<λ<p+1, we obtain:
Ψ′(t)≥λεH(t)+ε∫Ω(φ0φ1+ψ0ψ1)dx+ε(1+λ2)∫t0∫Ω(φ2t+ψ2t)dxds+λε2∫t0∫Ω(|∇φ|2+|∇ψ|2)dxds+ε(1−BMK2−λ)∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds−λε(dt+l)∫Ω(φ20+ψ20)dx−εBMk2. | (3.23) |
We choose φ0,φ1,ψ0, and ψ1 such that
∫Ω(φ0φ1+ψ0ψ1)dx−λl∫Ω(φ20+ψ20)dx>0, |
and M is sufficiently large so that
∫Ω(φ0φ1+ψ0ψ1)dx−λ(dT+l)∫Ω(φ20+ψ20)dx>BMK2. |
Next, we select ˜b such that 0<˜b<1−BMK2−λ. Inequality (3.23) takes the form
Ψ′(t)≥λεH(t)+ε(λ+22)∫t0∫Ω(φ2t+ψ2t)dxds+ε˜b∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds. | (3.24) |
Step 3. Alternatively, we have:
Ψ11−γ(t)≤211−γ[H(t)+ε11−γ(∫t0∫Ω(φφt+ψt)dxds)11−γ], and (∫t0∫Ω(φφt+ψψt)dxds)11−γ≤(∫t0∫Ωφφtdxds)11−γ+(∫t0∫Ωψψtdxds)11−γ. | (3.25) |
Utilizing the Cauchy-Schwarz and Hölder's inequalities, one observes:
(∫t0∫Ωφφtdxds)11−γ≤(∫t0(∫Ω|φ|2dx)12(∫Ω|φt|2dx)12ds)11−γ≤c3(∫t0(∫Ω|φ|p+1dx)1p+1(∫Ω|φt|2dx)12ds)11−γ≤c3(∫t0(∫Ω|φ|p+1dx)2p+1ds)12(1−γ)(∫t0∫Ω|φt|2dxds)12(1−γ), |
for some c3>0. By applying Young's inequality and Hölder's inequality again, we obtain
(∫t0∫Ωφφtdxds)11−γ≤c4[∫t0∫Ω|φt|2dxds+(∫t0(∫Ω|φ|p+1dx)2p+1ds)11−2γ]≤c4[∫t0∫Ω|φt|2dxds+(∫t0ds)μ∫t0∫Ω|φ|p+1dxds], | (3.26) |
for some c4>0 and μ=p−1(p+1)(1−2γ). Similarly,
(∫t0∫Ωψψtdxds)11−γ≤c5[∫t0∫Ω|ψt|2dxds+Tμ∫t0∫Ω|ψ|p+1dxds]. | (3.27) |
From (3.26) and (3.27), we have
(∫t0∫Ω(φφt+ψψt)dxds)11−γ≤λ[∫t0∫Ω(φ2t+ψ2t)dxds+Tμ∫t0∫Ω(|φ|p+1+|ψ|p+1)dxds]≤λ[∫t0∫Ω(φ2t+ψ2t)dxds+Tμ∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds]. | (3.28) |
Substituting (3.28) in (3.25), we get:
Ψ11−γ(t)≤211−γH(t)+211−γε11−γλ[∫t0∫Ω(φ2t+ψ2t)dxds+Tμ∫t0∫Ω(1p+1|φ+ψ|p+1+2p+1|φ|p+12|ψ|p+12)dxds]. | (3.29) |
Obviously, from (3.24) and (3.29), the inequality
Ψ11−γ(t)≤RΨ′(t), | (3.30) |
holds for sufficiently large R. Integrating (3.30) implies
Ψ11−γ(t)≥1Ψ−11−γ(0)−γtR(1−γ), | (3.31) |
provided that the expression Ψ−11−γ(0)−γtR(1−γ) is positive. Therefore, ψ(t) blows up at a certain moment
T⋆≤R(1−γ)γΨ11−γ(0)<T. |
Noting that T∗ depends on the fractional orders α and β, we observe that R, which is essential in determining T∗, is related to ˜b as expressed in Eq (3.24). Moreover, ˜b depends on B, where B is explicitly a function of α and β, as defined after Eq (3.22). This establishes the connection between the blow-up time and the fractional parameters, further reinforcing the influence of these orders on the solution's behavior.
In this study, we have investigated the blow-up behavior of solutions to a system of nonlinear wave equations with fractional damping and source terms. Our analysis demonstrates that the introduction of fractional derivatives significantly impacts the stability and blow-up characteristics of the system. We provided a detailed examination of how fractional orders influence the blow-up time T∗, revealing that T∗ is intricately related to the fractional orders α and β. Specifically, we established that T∗ depends on the parameters involved in the system, including those introduced in the energy functional.
Our results offer new insights into the dynamics of systems with fractional damping, which are often encountered in real-world applications involving complex dissipation and memory effects. The approach of combining classical analysis with fractional calculus not only enriches the theoretical understanding but also opens new avenues for practical applications in material science and engineering.
Future research could explore further implications of fractional derivatives in other nonlinear systems and extend the analysis to higher dimensions or different types of sources and damping mechanisms. Overall, this work contributes to a deeper understanding of fractional wave equations and their behavior, providing a foundation for subsequent studies in this evolving field.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors extend their appreciation to the Deanship of Research and Graduate Studies at King Khaled University for funding this work through Large Research Project under grant number RGP2/37/45.
The authors declare there is no conflicts of interest.
[1] |
M. Caputo, Linear Models of Dissipation whose Q is almost Frequency Independent—Ⅱ, Geophys. J. Int., 13 (1967), 529–539. https://doi.org/10.1111/j.1365-246X.1967.tb02303.x doi: 10.1111/j.1365-246X.1967.tb02303.x
![]() |
[2] |
K. Agre, M.A. Rammaha, System of nonlinear wave equations with damping and source terms, Differ. Integr. Equations, 19 (2006), 1235–1270. https://doi.org/10.57262/die/1356050301 doi: 10.57262/die/1356050301
![]() |
[3] |
C. O. Alves, M. M. Cavalcanti, On existence, uniform decay rates and blow up for solutions of the 2-D wave equation with exponential source, Calc. Var. Partial Differ. Equations, 34 (2009), 377–411. https://doi.org/10.1007/s00526-008-0188-z doi: 10.1007/s00526-008-0188-z
![]() |
[4] |
X. Han, M. Wang, Global existence and blow-up of solution for a system of nonlinear viscoelastic wave equations with damping and source, Nonlinear Anal., 71 (2009), 5427–5450. https://doi.org/10.1016/j.na.2009.04.031 doi: 10.1016/j.na.2009.04.031
![]() |
[5] |
M. Kirane, N. Tatar, Exponential growth for fractionally damped wave equations, Z. Anal. Anwend., 22 (2003), 167–178. https://doi.org/10.4171/zaa/1137 doi: 10.4171/zaa/1137
![]() |
[6] | D. Matignon, J. Audounet, G. Montseny, Energy decay for wave equations with damping of fractional order, 1998. Available from: https://www.researchgate.net/publication/2757992_Energy_Decay_for_Wave_Equations_with_Damping_of_Fractional_order. |
[7] |
N. Tatar, A blow-up result for a fractionally damped wave equation, Nonlinear Differ. Equations Appl. NoDEA, 12 (2005), 215–226. https://doi.org/10.1007/s00030-005-0015-6 doi: 10.1007/s00030-005-0015-6
![]() |
[8] |
M. R. Alaimia, N. Tatar, Blow up for the wave equation with a fractional damping, J. Appl. Anal., 11 (2005), 133–144. https://doi.org/10.1515/JAA.2005.133 doi: 10.1515/JAA.2005.133
![]() |
[9] | J. M. Ball, Remarks on blow up and non-existence theorems for nonlinear evolution equations, Quart. J. Math. Oxford, 28 (1977), 473–486. |
[10] |
H. A. Levine, Instability and nonexistence of global solutions of nonlinear wave equations of the form Putt=Au+F(u), Trans. Am. Math. Soc., 192 (1974), 1–21. https://doi.org/10.2307/1996814 doi: 10.2307/1996814
![]() |
[11] |
H. A. Levine, Some additional remarks on the nonexistence of global solutions to nonlinear wave equations of the form Putt=Au+F(u), SIAM J. Math. Anal., 5 (1974), 138–146. https://doi.org/10.1137/0505015 doi: 10.1137/0505015
![]() |
[12] |
V. Georgiev, G. Todorova, Existence of solutions of the wave equation with nonlinear damping and source terms, J. Differ. Equations, 109 (1994), 295–308. https://doi.org/10.1006/jdeq.1994.1051 doi: 10.1006/jdeq.1994.1051
![]() |
[13] |
J. M. Ball, On the asymptotic behavior of generalized processes, with applications to nonlinear evolution equations, J. Differ. Equations, 27 (1978), 224–265. https://doi.org/10.1016/0022-0396(78)90032-3 doi: 10.1016/0022-0396(78)90032-3
![]() |
[14] |
J. A. Nohel, D. F. Shea, Frequency domain methods for Volterra equations, Adv. Math., 22 (1976), 278–304. https://doi.org/10.1016/0001-8708(76)90096-7 doi: 10.1016/0001-8708(76)90096-7
![]() |
[15] | G. Gripenberg, S. O. Londen, O. Staffans, Volterra Integral and Functional Equations, Cambridge University Press, Cambridge, 1990. https://doi.org/10.1017/CBO9780511662805 |