Next Article in Journal
Heyting Locally Small Spaces and Esakia Duality
Next Article in Special Issue
Averaging Principle for ψ-Capuo Fractional Stochastic Delay Differential Equations with Poisson Jumps
Previous Article in Journal
Symmetry Perception and Psychedelic Experience
Previous Article in Special Issue
Numerical Simulation for a Hybrid Variable-Order Multi-Vaccination COVID-19 Mathematical Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Weakly Coupled System of Semi-Linear Fractional θ-Evolution Equations with Special Cauchy Conditions

by
Abdelhamid Mohammed Djaouti
Preparatory Year Deanship, King Faisal University, Hofuf 31982, Saudi Arabia
Symmetry 2023, 15(7), 1341; https://doi.org/10.3390/sym15071341
Submission received: 10 May 2023 / Revised: 15 June 2023 / Accepted: 25 June 2023 / Published: 30 June 2023
(This article belongs to the Special Issue Applied Mathematics and Fractional Calculus II)

Abstract

:
In this paper, we consider a weakly system of fractional θ -evolution equations. Using the fixed-point theorem, a global-in-time existence of small data solutions to the Cauchy problem is proved for one single equation. Using these results, we prove the global existence for the system under some mixed symmetrical conditions that describe the interaction between the equations of the system.

1. Introduction

In this paper, we show the existence of the global (in time) solutions with small data to the weakly coupled system of fractional wave equations
D 1 + λ 1 μ + ( Δ ) θ 1 2 u = | v | p , J 1 λ 1 u ( 0 , x ) = u λ 1 ( x ) , D λ 1 u ( 0 , x ) = 0 , D 1 + λ 2 v + ( Δ ) θ 2 2 v = | u | q , J 1 λ 2 u ( 0 , x ) = u λ 2 ( x ) , D λ 2 v ( 0 , x ) = 0 ,
where λ 1 , λ 2 ( 0 , 1 ) , θ 1 , θ 2 are real positive numbers and D 1 + λ is the Riemann–Liouville fractional derivative defined by
D 1 + λ f ( t ) : = t 2 ( J 1 λ f ) ( t )
with the Riemann–Liouville fractional integral operator
D a f ( t ) : = 1 Γ ( a ) 0 t ( t s ) a 1 f ( s ) d s , t > 0
for ( a ) > 0 , and Γ is the Euler Gamma function.
Such mathematical models have promising applications in engineering and in other physical sciences, as well as in numerical simulations of some fractional nonlinear viscoelastic flow problems, and they impact the bioconvection on the free stream flow of a pseudoplastic nanofluid past a rotating cone.
At the outset, since the fractional equation interpolates the heat equation for λ 0 and the wave equation for λ 1 we will provide briefly some previous results of the wave equations and heat equation.
On the one hand, we consider the Cauchy problem for the semi-linear heat equation
u t Δ u = | u | p , u ( 0 , x ) = u 0 ( x ) .
Fujita in [1] proved that the exponent p F u j : = 1 + 2 n is critical for the classical heat model, which means that we have the global (in time) existence of small data solutions for p > p c r i t , and the blow up if we have the inverse 1 < p < p F u j . In [2,3], the authors proved the blow-up for the critical case p = p F u j .
On the other hand, let us consider the Cauchy problem for the semi-linear wave equation
u t t Δ u = | u | p , u ( 0 , x ) = u 0 ( x ) , u t ( 0 , x ) = u 1 ( x ) ,
where the authors in [4] proved for n = 3 that the critical exponent is defined as a positive root of the quadratic equation
( n 1 ) p 2 ( n + 1 ) p 2 = 0 .
The defined exponent by the last equation is called the Strauss exponent and denoted by p S for further considerations, which means that we have the global (in time) existence of small data weak solutions for the above p S , whereas the local (in time) existence for p > 1 and large data can be only expected. In [5,6], the author proved in R 2 that the Strauss exponent p S is critical. After that, the global existence for n = 2 , 3 was treated in [7] and for n 4 in [8,9]. The nonexistence of solutions for data compactly supported was studied in [10] for 1 < p < n + 1 n 1 . For n = 3 , the authors proved some optimal results in [11] for p = 1 + 2 . For n > 3 , a nonexistence result with small data proved in [12] for 1 < p < p S .
In 2017, D’Abbicco et al. [13] considered the semi-linear fractional wave equation
t 1 + λ u Δ u = | u | p , u ( 0 , x ) = u 0 ( x ) , u t ( 0 , x ) = 0 ,
where λ ( 0 , 1 ) with the fractional Riemann–Liouville fractional derivative. They proved the critical exponent for the global existence of a small data solution in a low space dimension. The Caputo fractional order and the existence of non-null Cauchy data was studied in [14].
In [15], the authors proved the global (in time) existence of small data solutions to semi-linear fraction θ -evolution equations with mass or power nonlinearity. A similar problem was treated in [16] by considering a memory term instead of the power nonlinearity.
In the first part of our main results, we show the global existence of a small data solution to the fractional Riemann–Liouville order to the semi-linear θ -evolution problem (7).
For the systems, let us first consider the weakly coupled system of damped wave equations semi-linear heat equations
u t Δ u = | v | p , u ( 0 , x ) = u 0 ( x ) , u t ( 0 , x ) = u 1 ( x ) , v t Δ v = | u | q , v ( 0 , x ) = v 0 ( x ) , v t ( 0 , x ) = v 1 ( x ) ,
where t [ 0 , ) , x R n , p , q > 1 and p q > 1 . The authors of [17] showed that the exponents p and q satisfying
n 2 = max { p , q } + 1 p q 1
are critical, which means that the solutions exist globally for n 2 > max { p , q } + 1 p q 1 and blow-up for the inverse case. For more details about the system of damped wave equations semi-linear heat equations, the reader can also see [18,19,20,21].
Some papers are considered for the weakly coupled systems of semilinear classical damped wave equations with power non-linearities. The problem we have in mind is
u t t Δ u + u t = | v | p , u ( 0 , x ) = u 0 ( x ) , u t ( 0 , x ) = u 1 ( x ) , v t t Δ v + v t = | u | q , v ( 0 , x ) = v 0 ( x ) , v t ( 0 , x ) = v 1 ( x ) ,
where t [ 0 , ) , x R n . In 2007, Sun and Wang proved in [22] that if
λ : = max { p ; q } + 1 p q 1 < n 2 .
for n = 1 or 3 , then the solution exists globally in time for small initial data, while, if λ n 2 , then every solution having positive average value does not exist globally. In [23], the authors generalized the previous results to the case where n = 1 , 2 , 3 and improved the time decay estimates for n = 2 . In 2014, using the weighted energy method, Nishihara and Wakasugi proved, in [24], the critical exponent for any space dimensions. Considering the time-dependent dissipation terms, the authors of [25,26,27] proved the global (in time) existence of small data solutions under a plan condition, which presents the interplay between the exponents of power nonlinearities.
In our paper, we consider first the single equation from system (1) where we proved the global existence for some range of the exponent p under conditions related to the regularity of the data and the dimension. After that, we apply the results of the single equation to study the weakly coupled systems (1). We proved the global existence for the system with a loss of decay if one of the exponents of power nonlinearities did not satisfy the condition of the single equation.
The paper is organized as follows. In Section 2, we will show our main results of global (in time) existence with examples. Moreover, we mention some remarks of the interpolated cases of wave and heat equations. Next, in Section 3, we prove the existence of solution by applying Banach’s fixed point. Appendix A concludes the paper.

2. Main Results

2.1. Single Equation of Fractional Integral Equation

In this section, we will show our main results where we start with the global (in time) existence of solutions to the single equation of the Cauchy problem. Using the formal representation of the solution to our equation, we obtain the estimates of the solutions, and finally we prove the existence using fixed-point theorem explained in the Appendix A.
D λ + 1 u + ( Δ ) θ 2 u = | u | p , J 1 λ u ( 0 , x ) = u λ ( x ) , D λ u ( 0 , x ) = 0
where λ ( 0 , 1 ) , θ > 0 .
Theorem 1.
Let n 1 , and the data u λ are supposed to belong to L 1 L p . The following conditions are satisfied for the exponent p:
p > 1 + 1 + λ n θ ( 1 + λ ) λ ,
and
p < 1 + θ n θ i f n > θ .
Then, a small constant ϵ exists such that, if u λ L 1 L p ϵ , then there is a uniquely determined globally (in time) energy solution to (7) in C ( [ 0 , ) , L 1 L p ) .
Furthermore, the solution satisfies the estimates:
u L q ( 1 + t ) λ n θ ( 1 + λ ) 1 1 q u λ L 1 L p ,
where q [ 1 , p ] .
The new type of date has a strong influence in the representation of the solution of (1) after [28], which leads to a quite different admissible range of the exponent p compared with the classical equations presented in [14].
Remark 1.
If λ 0 , then the admissible range for the global (in time) existence corresponds with a Fujita like exponent 1 + θ n . On the contrary for λ 1 , we obtain a gap of continuity with respect to the Strauss exponent, which appeared in previous results as a critical exponent for the classical wave equation.
Remark 2.
One can obtain the optimal for the exponent p in (8) using the scaling argument similarity to prove of the critical exponent to (4) illustrated in [14].
Example 1.
We consider a concrete example by giving values to the parameters appearing in the theorem. Let us consider in R 3 the following model:
D 3 2 u + ( Δ ) 3 4 u = | u | p , J 1 2 u ( 0 , x ) = u λ ( x ) , D 1 2 u ( 0 , x ) = 0 .
Then, using Theorem 1, the admissible range for the global existence is 8 5 < p < 2 .

2.2. Weakly Coupled System of Fractional Integral Equations

In this section, we apply the results of the previous theorem to study systems of weakly coupled fractional θ -evolution equations.
Theorem 2.
Let n 1 , and the data u λ 1 , u λ 2 is supposed to belong to ( L 1 L p ) × ( L 1 L q ) . The following conditions are satisfied for the exponent p and q:
p < 1 + 1 + λ 2 n θ 2 ( 1 + λ 2 ) λ 2 , q > 1 + 1 + λ 1 n θ 1 ( 1 + λ 1 ) λ 1 ,
p < 1 + θ n θ 1 , q < 1 + θ n θ 2 i f n > min { θ 1 ; θ 2 }
and
Q ( λ 1 , λ 2 , θ 1 , θ 2 , q ) > o .
Then, a small constant ϵ exists such that, if u λ 1 L 1 L p + v λ 2 L 1 L q ϵ , then there is a uniquely determined globally (in time) energy solution to (1) in C ( [ 0 , ) , L 1 L p ) × C ( [ 0 , ) , L 1 L q ) .
Furthermore, the solution satisfies the estimates:
u L r 1 ( 1 + t ) λ + L ( p ) n θ 1 ( 1 + λ 1 ) 1 1 r 1 u λ L 1 L p ,
v L r 2 ( 1 + t ) λ 2 n θ 2 ( 1 + λ 2 ) 1 1 r 2 v λ L 1 L q ,
where L ( p ) = n θ 2 ( 1 + λ 2 ) ( p 1 ) + p λ 2 , Q ( λ 1 , λ 2 , θ 1 , θ 2 , q ) = n θ 2 ( 1 + λ 2 ) λ 2 q 2 n θ 1 ( 1 + λ 1 ) n θ 2 ( 1 + λ 2 ) λ 1 q n θ 1 ( 1 + λ 1 ) and r 2 [ 1 , p ] , r 2 [ 1 , q ] .
Remark 3.
If we take in Theorem 2 the condition p > 1 + 1 + λ 2 n θ 2 ( 1 + λ 2 ) λ 2 , then we cannot feel any interplay between the equations of the system since it will behave as a single equation.
Remark 4.
If we consider p = 1 + 1 + λ 2 n θ 2 ( 1 + λ 2 ) λ 2 then, after using Proposition A1, we obtain a new decay generated by the log term appearing in the estimate of u, exactly, ( 1 + t ) 1 log ( 1 + t ) ( 1 + t ) 1 + ε .
Example 2.
Let us consider θ 1 = θ 2 = 2 in R 2 and the parameter of the fractional derivative of the first equation λ 1 0 and the second λ 2 1 . Then, with the Cauchy condition the model, we obtain
t u + Δ u = | v | p , t t v + Δ v = | u | q .
Applying Theorem 2, we obtain the global (in time) existence of the solution for p < 3 and q > 2 .
Remark 5.
The reader can apply the last theorem for several examples. Giving values to some parameters such as the dimension or the order of the fractional derivative, we obtain the mixed condition that leads to the global existence.

3. Philosophy of Our Approach

In this section, we will prove results for the Cauchy problems (1) and (7). Our main interest is to prove the global (in time) existence of small data solutions, which means the global existence after the perturbation of the null Cauchy condition u λ L 1 L p ϵ . Such results imply immediate stability results for the zero solution.

3.1. Proof of Theorem 1

In this section, we deal with the following single equation:
t λ + 1 u + ( Δ ) θ 2 u = | u | p , J 1 λ u ( 0 , x ) = u λ ( x ) , D λ u ( 0 , x ) = 0 .
We define the norm of the solution space X ( t ) , which we will propose in all of the proofs of the above theorems by
u X ( t ) = sup τ [ 0 , t ] ( 1 + t ) λ u ( τ , · ) L 1 + ( 1 + t ) n θ ( 1 + λ ) 1 1 p u ( τ , · ) L p ,
We introduce the operator N by
N : u X ( t ) N u = N u ( t , x ) : = u l n ( t , x ) + u n l ( t , x ) ,
where u l n is a Sobolev solution to the Cauchy problem
t λ + 1 u + ( Δ ) θ 2 u = 0 , J 1 λ u ( 0 , x ) = u λ ( x ) , D λ u ( 0 , x ) = 0 ,
and u n l is a Sobolev solution to the Cauchy problem
t λ + 1 u + ( Δ ) θ 2 u = | u | p , J 1 λ u ( 0 , x ) = u λ ( x ) , D λ u ( 0 , x ) = 0 .
Using Fourier analysis together with Theorem A1 from Appendix A, we can show that the solutions of the previous problems can be presented by u ( t , x ) = u l n ( t , x ) + u n l ( t , x ) as follows:
u l n ( t , x ) = t λ 1 F 1 E 1 + λ , λ t 1 + λ | ξ | θ ( t , x ) ( x ) u λ ( x ) ,
and
u n l ( t , x ) = 0 t ( t s ) λ F 1 E 1 + λ , 1 + λ t 1 + λ | ξ | θ ( t s , x ) ( x ) | u ( s , x ) | p d s .
Following Proposition A2, our aim is to prove the following inequalities:
N u X ( t ) u λ L 1 L p + u X ( t ) p ,
N u N v X ( t ) u v X ( t ) u X ( t ) p 1 + v X ( t ) p 1 .
After proving these both inequalities, we apply Banach’s fixed-point theorem. In this way, we obtain the local (in time) existence of large data Sobolev solutions and the global (in time) existence of small data Sobolev solutions as well.
We split the prove of the first inequality (17) into the following inequalities:
u l n X ( t ) u λ L 1 L p ,
and
u n l X ( t ) u X ( t ) p .
To prove inequality (19) we have to derive the estimate of F 1 E 1 + λ , λ t 1 + λ | ξ | θ L p in order to use Young’s inequality. Using the scaling property, we obtain
F 1 E 1 + λ , λ t 1 + λ | ξ | θ L p = t n θ ( 1 + λ ) 1 1 p F 1 E 1 + λ , λ | ξ | θ L p .
Indeed, after change of variable ξ 1 = t 1 + λ | ξ | we obtain
F 1 G t 1 + λ | ξ | θ = t n θ ( 1 + λ ) R n e i t 1 + λ θ x ξ 1 G ( | ξ 1 | θ ) d ξ 1 = t n θ ( 1 + λ ) F 1 G ( | ξ | θ ) ( t 1 + λ θ x ) .
Using the last equality, we obtain
F 1 G t 1 + λ | · | θ L p p = t n θ ( 1 + λ ) p F 1 G ( | · | θ ) ( t 1 + λ θ x ) L p = t n θ ( 1 + λ ) p R n F 1 G ( | · | θ ) ( t 1 + λ θ x ) p d x .
The change of variable y = t 1 + λ θ x leads to
F 1 G t 1 + λ | · | θ L p p = t n θ ( 1 + λ ) p + n θ ( 1 + λ ) F 1 G | · | θ L p p ,
which completes the proof of 21.
Then, we restrict ourselves to the estimates of F 1 E 1 + λ , λ | ξ | θ L p . After applying Theorem A2 from the Appendix A, we obtain
E 1 + λ , λ | ξ | θ = 2 1 + λ | ξ | θ 1 2 1 + λ e θ 1 + λ cos π 1 + λ c o s | ξ | θ 1 + λ s i n π 1 + λ + π 1 | ξ | θ 1 2 1 + λ 0 s 2 + λ s 2 ( 1 + λ ) + 2 c o s ( π ( 1 + λ ) ) + 1 e s | ξ | θ 1 + λ d s sin ( λ π ) ,
which leads to
F 1 E 1 + λ , λ | ξ | θ = 2 1 + λ A ( s , x ) + π 1 sin ( λ π ) 0 s 2 + λ s 2 ( 1 + λ ) + 2 c o s ( π ( 1 + λ ) ) + 1 B ( s , x ) d s ,
where
A ( s , x ) = F 1 | ξ | θ 1 2 1 + λ e θ 1 + λ cos π 1 + λ c o s | ξ | θ 1 + λ s i n π 1 + λ , B ( s , x ) = F 1 | ξ | θ 1 2 1 + λ e s | ξ | θ 1 + λ .
First, we consider B ( s , x ) . Similarly to (21), we have
B ( s , · ) L p = F 1 | ξ | θ 1 2 1 + λ e s | ξ | θ 1 + λ L p = F 1 s ( 1 + λ ) 1 2 1 + λ ( s 1 + λ | ξ | θ ) 1 2 1 + λ e ( s 1 + λ | ξ | θ ) 1 1 + λ L p = s λ 1 n θ ( 1 + λ ) 1 1 p F 1 | ξ | θ 1 2 1 + λ e | ξ | θ 1 + λ L p = s λ 1 n θ ( 1 + λ ) 1 1 p B ( 1 , · ) L p .
Then,
B ( s , · ) L p = s λ 1 n θ ( 1 + λ ) 1 1 p B ( 1 , · ) L p .
Then,
F 1 E 1 + λ , λ | · | θ L p = 2 1 + λ A ( s , · ) L p + π 1 sin ( λ π ) 0 s 1 + 2 λ n θ ( 1 + λ ) 1 1 p s 2 ( 1 + λ ) + 2 c o s ( π ( 1 + λ ) ) + 1 B ( 1 , · ) L p d s ,
Using the last estimate together with (A5) from Remark A1, one can obtain the following estimate from Lemma 2.1 in [14] for d = θ 1 2 1 + λ :
F 1 E 1 + λ , λ | ξ | θ L p i f n θ 1 1 p < 2 ,
which satisfied (9).
From (15) with (21), and after using Young’s inequality, we obtain
u l n ( t , x ) L 1 ( 1 + t ) λ 1 u λ L 1 ,
u l n ( t , x ) L p ( 1 + t ) λ 1 n θ ( 1 + λ ) 1 1 p u λ L 1 + u λ L p .
Replacing last estimates in the definition of the norm of solution space (14) leads to the desired estimate (19).
For the second estimate (20), under the same conditions requested for (23) we have
F 1 E 1 + λ , 1 + λ t 1 + λ | ξ | θ L p ( 1 + t ) n θ ( 1 + λ ) 1 1 p .
From (16), we obtain
u n l ( t , x ) L 1 0 t ( t s ) λ | u ( s , x ) | p L 1 d s ,
and
u n l ( t , x ) L p 0 t ( t s ) λ n θ ( 1 + λ ) 1 1 p | u ( s , x ) | p L 1 d s .
Using the definition of solution space from (14), we obtain
u n l ( t , x ) L 1 u X ( t ) p 0 t ( t s ) λ ( 1 + s ) n θ ( 1 + λ ) ( p 1 ) + p λ d s , u n l ( t , x ) L p u X ( t ) p 0 t ( t s ) λ n θ ( 1 + λ ) 1 1 p ( 1 + s ) n θ ( 1 + λ ) ( p 1 ) + p λ d s .
Using Proposition A1, we obtain
u n l ( t , x ) L 1 ( 1 + t ) λ u X ( t ) p ,
u n l ( t , x ) L p ( 1 + t ) λ n θ ( 1 + λ ) 1 1 p u X ( t ) p ,
provided that n θ ( 1 + λ ) 1 1 p λ < 1 and n θ ( 1 + λ ) ( p 1 ) p λ > 1 , which are equivalent to (8) and (9), respectively.
Replacing the last estimates in the norm of solution space, we obtain (20), which complete, together with (19), the proof of the first inequality (17).
For the second condition (18), we assume that u and v belong to X ( t ) . Then,
N u N v = 0 t ( t s ) λ F 1 E 1 + λ , 1 + λ t 1 + λ | ξ | θ ( t s , x ) ( x ) | u ( s , x ) | p | v ( s , x ) | p d s .
We control all norms appearing in N u N v X ( t ) . These are the norms N u N v L 1 and | D | s ( N u N v ) L p .
Similarly to (26), we have
N u N v L 1 0 t ( t s ) λ | u ( s , x ) | p | v ( s , x ) | p L 1 d s .
Hölder’s inequality implies
| u ( s , · ) | p | v ( s , · ) | p L 1 u ( s , · ) v ( s , · ) L p u ( s , · ) L p p 1 + v ( s , · ) L p p 1 ,
Using the norm of the solution space X ( t ) , we obtain
u ( s , · ) v ( s , · ) L p ( 1 + s ) n θ ( 1 + λ ) ( 1 1 p ) + λ u ( s , · ) v ( s , · ) X ( t ) , u ( s , · ) L p p 1 ( 1 + s ) n θ ( 1 + λ ) ( 1 1 p ) + λ ( p 1 ) v ( s , · ) X ( t ) ( p 1 ) , u ( s , · ) L p p 1 ( 1 + s ) n θ ( 1 + λ ) ( 1 1 p ) + λ ( p 1 ) v ( s , · ) X ( t ) ( p 1 ) .
Using the last estimates, we can obtain similarly to (28) and (29)
N u N v L 1 ( 1 + t ) λ u v X ( t ) u X ( t ) p 1 + v X ( t ) p 1 ,
and
N u N v L p ( 1 + t ) λ n θ ( 1 + λ ) 1 1 p u v X ( t ) u X ( t ) p 1 + v X ( t ) p 1 .
Then, the proof of the second condition and the theorem is completed.

3.2. Proof of Theorem 2

We define the norm of the solution space X ( t ) by
( u , v ) X ( t ) = sup τ [ 0 , t ] M ( τ , u ) + M ( τ , v )
where
M ( τ , u ) = ( 1 + t ) λ 1 L ( p ) u ( τ , · ) L 1 + ( 1 + t ) n θ 1 ( 1 + λ 1 ) 1 1 p u ( τ , · ) L p ,
M ( τ , v ) = ( 1 + t ) λ 2 v ( τ , · ) L 1 + ( 1 + t ) n θ 2 ( 1 + λ 2 ) 1 1 q v ( τ , · ) L q .
Then, we introduce the operator N by
N : ( u , v ) X ( t ) N ( u , v ) = ( u l n + u n l , v l n + v n l ) ,
where
u l n ( t , x ) : = t λ 1 1 F 1 E 1 + λ 1 , λ 1 t 1 + λ 1 | ξ | θ 1 ( t , x ) ( x ) u λ 1 ( x ) , u n l ( t , x ) : = 0 t ( t s ) λ 1 F 1 E 1 + λ 1 , 1 + λ 1 t 1 + λ 1 | ξ | θ 1 ( t s , x ) ( x ) | v ( s , x ) | p d s , v l n ( t , x ) : = t λ 2 1 F 1 E 1 + λ 2 , λ 1 t 1 + λ 2 | ξ | θ 2 ( t , x ) ( x ) v λ 2 ( x ) , v n l ( t , x ) : = 0 t ( t s ) λ 2 F 1 E 1 + λ 2 , 1 + λ 2 t 1 + λ 2 | ξ | θ 2 ( t s , x ) ( x ) | u ( s , x ) | q d s .
If we consider the results Proposition A3, then our aim is to prove the following inequalities, which imply, among other things, the global existence of small data solutions:
N ( u , v ) X ( t ) u λ 1 L 1 L p + v λ 2 L 1 L q + ( u , v ) X ( t ) p + ( u , v ) X ( t ) q ,
N ( u , v ) N ( u ˜ , v ˜ ) X ( t ) ( u , v ) ( u ˜ , v ˜ ) X ( t ) ( ( u , v ) X ( t ) p 1 + ( u ˜ , v ˜ ) X ( t ) p 1 + ( u , v ) X ( t ) q 1 + ( u ˜ , v ˜ ) X ( t ) q 1 ) .
Let us start by the first condition. Similarly to (24) and (25), we obtain
u l n ( t , x ) L 1 ( 1 + t ) λ 1 1 u λ 1 L 1 , u l n ( t , x ) L p ( 1 + t ) λ 1 1 n θ 1 ( 1 + λ 1 ) 1 1 p u λ 1 L 1 + u λ 1 L p , v l n ( t , x ) L 1 ( 1 + t ) λ 2 1 v λ 2 L 1 , v l n ( t , x ) L q ( 1 + t ) λ 2 1 n θ 2 ( 1 + λ 2 ) 1 1 q v λ 2 L 1 + u λ 2 L q .
The last estimates, together with the definition of the norm in (31), lead to
( u l n , v l n ) X ( t ) u λ 1 L 1 L p + v λ 2 L 1 L q .
Then, we complete the proof by showing the inequality
( u n l , v n l ) X ( t ) ( u , v ) X ( t ) p + ( u , v ) X ( t ) q .
For u n l , we have
u n l ( t , x ) L 1 0 t ( t s ) λ 1 | v ( s , x ) | p L 1 d s ,
and
u n l ( t , x ) L p 0 t ( t s ) λ 1 n θ 1 ( 1 + λ 1 ) 1 1 p | v ( s , x ) | p L 1 d s .
Using the definition of solution space from (31), we obtain
u n l ( t , x ) L 1 ( u , v ) X ( t ) p 0 t ( t s ) λ 1 ( 1 + s ) n θ 2 ( 1 + λ 2 ) ( p 1 ) + p λ 2 d s , u n l ( t , x ) L p ( u , v ) X ( t ) p 0 t ( t s ) λ 1 n θ 1 ( 1 + λ 1 ) 1 1 p ( 1 + s ) n θ 2 ( 1 + λ 2 ) ( p 1 ) + p λ 2 d s .
From Proposition A1, one can obtain
u n l ( t , x ) L 1 ( u , v ) X ( t ) p ( 1 + t ) λ 1 n θ 2 ( 1 + λ 2 ) ( p 1 ) + p λ 2 = ( u , v ) X ( t ) p ( 1 + t ) λ 1 + L ( p ) ,
u n l ( t , x ) L p ( u , v ) X ( t ) p ( 1 + t ) λ 1 n θ 1 ( 1 + λ 1 ) 1 1 p n θ 2 ( 1 + λ 2 ) ( p 1 ) + p λ 2 = ( u , v ) X ( t ) p ( 1 + t ) λ 1 n θ 1 ( 1 + λ 1 ) 1 1 p + L ( p ) ,
provided that n θ ( 1 + λ 1 ) 1 1 p λ 1 < 1 , which is equivalent to (11).
For u n l , we have
v n l ( t , x ) L 1 0 t ( t s ) λ 2 | u ( s , x ) | q L 1 d s ,
and
v n l ( t , x ) L q 0 t ( t s ) λ 2 n θ 2 ( 1 + λ 2 ) 1 1 q | u ( s , x ) | q L 1 d s .
Using the norm of the solution space, we obtain
v n l ( t , x ) L 1 ( u , v ) X ( t ) p 0 t ( t s ) λ 2 ( 1 + t ) Q ( λ 1 , λ 2 , θ 1 , θ 2 , q ) d s ,
and
v n l ( t , x ) L q ( u , v ) X ( t ) p 0 t ( t s ) λ 2 n θ 2 ( 1 + λ 2 ) 1 1 q ( 1 + t ) Q ( λ 1 , λ 2 , θ 1 , θ 2 , q ) d s .
Proposition A1, together with (12), leads to
v n l ( t , x ) L 1 ( u , v ) X ( t ) p ( 1 + t ) λ 2 ,
and
v n l ( t , x ) L q ( u , v ) X ( t ) p ( 1 + t ) λ 2 n θ 2 ( 1 + λ 2 ) 1 1 q d s ,
provided that n θ ( 1 + λ 2 ) 1 1 q λ 2 < 1 , which is equivalent to (11).
From (38) to (41), we obtain (37), which implies, together with (36), the first condition (34).
To prove (35), we assume that ( u , v ) and ( u ˜ , v ˜ ) are two elements from the function space X ( t ) . Then, we have
N ( u , v ) N ( u ˜ , v ˜ ) = u n l ( t , x ) u ˜ n l ( t , x ) , v n l ( t , x ) v ˜ n l ( t , x )
= ( 0 t F 1 E 1 + λ 1 , 1 + λ 1 t 1 + λ 1 | ξ | θ 1 ( t s , x ) ( x ) | v ( s , x ) | p | v ˜ ( s , x ) | p d s ,
0 t F 1 E 1 + λ 2 , 1 + λ 2 t 1 + λ 2 | ξ | θ 2 ( t s , x ) ( x ) | u ( s , x ) | q | u ˜ ( s , x ) | q d s ) .
Similarly to the proof of the estimates (30), we can derive the following estimates for 0 τ t :
| v ( τ , · ) | p | v ˜ ( τ , · ) | p L 1 ( 1 + t ) n θ 2 ( 1 + λ 2 ) ( p 1 ) + p λ 2 v v ˜ X ( t ) v X ( t ) p 1 + v ˜ X ( t ) p 1 ,
| u ( τ , · ) | q | u ˜ ( τ , · ) | q L 1 ( 1 + t ) Q ( λ 1 , λ 2 , θ 1 , θ 2 , q ) u u ˜ X ( t ) u X ( t ) q 1 + u ˜ X ( t ) q 1 .
Using the last estimates, one may finally conclude, similarly to (38) to (41), the following estimates:
0 t F 1 E 1 + λ 1 , 1 + λ 1 t 1 + λ 1 | ξ | θ 1 ( t s , x ) ( x ) | v ( s , x ) | p | v ˜ ( s , x ) | p d s L 1 ( 1 + t ) λ 1 + L ( p ) v v ˜ X ( t ) v X ( t ) p 1 + v ˜ X ( t ) p 1 , 0 t F 1 E 1 + λ 1 , 1 + λ 1 t 1 + λ 1 | ξ | θ 1 ( t s , x ) ( x ) | v ( s , x ) | p | v ˜ ( s , x ) | p d s L p ( 1 + t ) λ 1 n θ 1 ( 1 + λ 1 ) 1 1 p + L ( p ) v v ˜ X ( t ) v X ( t ) p 1 + v ˜ X ( t ) p 1 , 0 t F 1 E 1 + λ 2 , 1 + λ 2 t 1 + λ 2 | ξ | θ 2 ( t s , x ) ( x ) | u ( s , x ) | q | u ˜ ( s , x ) | q d s ) L 1 ( 1 + t ) λ 2 v v ˜ X ( t ) v X ( t ) p 1 + v ˜ X ( t ) p 1 , 0 t F 1 E 1 + λ 2 , 1 + λ 2 t 1 + λ 2 | ξ | θ 2 ( t s , x ) ( x ) | u ( s , x ) | q | u ˜ ( s , x ) | q d s ) L q ( 1 + t ) λ 2 n θ 2 ( 1 + λ 2 ) 1 1 q v v ˜ X ( t ) v X ( t ) p 1 + v ˜ X ( t ) p 1 ,
In this way, we can conclude the proof of the last condition (35) and the theorem.

4. Concluding Remarks

  • We need to prove the blow-up for the system an interaction between the exponents of both equations. However, the method of scaling is not suitable to prove the blow-up result for the system since we have no interactions between the exponents. Moreover, the influence of each equation to the other one generated a condition presented by several parameters, fractional derivatives, dimensions, and others. For this reason, we will devote the blow-up problem in a forthcoming project using another approach.
  • The applications of our results in real world problems and phenomena can be investigated after mathematical modeling by choosing the suitable parameters involved in our problem, such as dimension, and by taking the experimental values into consideration.

Funding

This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Project No. GRANT3371).

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Project No. GRANT3371).

Conflicts of Interest

The author declares that there is no competing interest.

Appendix A

Theorem A1.
Let λ ( 0 , 1 ) , a λ R . Then, the unique solution solution to
t λ + 1 f + | ξ | θ f = g ( t ) , J 1 λ f ( 0 ) = a λ , D λ g ( 0 ) = 0 .
is given by
f ( t ) = t λ 1 E 1 + λ , λ t 1 + λ | ξ | θ a λ + 0 t ( t s ) λ E 1 + λ , 1 + λ t 1 + λ | ξ | θ ( t s , · ) g ( t ) d s ,
where E 1 + λ , μ are the Mittag–Leffler functions defined by
E 1 + λ , μ ( z ) = k = 0 z k Γ ( k + λ k + μ ) .
For the proof, see [28].
Theorem A2.
Let 0 < λ < 2 , μ R , and m N , with m μ 1 + λ 1 . Then, for the real number z > 0 , the following holds:
E 1 + λ , μ ( z 1 + λ ) = 2 1 + λ z 1 μ e z cos ( π 1 + λ ) cos ( z sin ( π 1 + λ ) π 1 + λ ( μ 1 ) )
+ k = 1 m ( 1 ) k 1 Γ ( μ k ( 1 + λ ) ) z k ( 1 + λ ) + Ω m ( z ) ,
where
Ω m ( z ) = ( 1 ) m z 1 μ π I 1 , m ( z ) sin ( π ( μ ( m + 1 ) ( 1 + λ ) ) + I 2 , m ( z ) sin ( π ( μ m ( 1 + λ ) ) ) ,
and
I j , m ( z ) = 0 s ( m + j ) ( 1 + λ ) μ s 2 ( 1 + λ ) 2 cos ( π ( 1 + λ ) ) s 1 + λ + 1 e s z d s .
Remark A1.
The integral I j , m ( z ) is uniformly bounded if
1 < m + j 1 + 1 μ 1 + λ < 1 .
For the proof, see [29].
Proposition A1.
Let a R < 1 and b R . Then,
0 t ( t s ) a ( 1 + s ) b d s ( 1 + t ) a i f a < 1 < b ; ( 1 + t ) 1 log ( 1 + t ) i f a < 1 = b ; ( 1 + t ) 1 a b i f a , b < 1 .
The reader can find the proof of Proposition A1 in [14].
Proposition A2.
The operator N maps X ( t ) into itself and has one and only one fixed point u X ( t ) if the following inequalities hold:
N u X ( t ) C 0 ( t ) ( u 0 , u 1 ) A m , s + C 1 ( t ) u X ( t ) p ,
N u N v X ( t ) C 2 ( t ) u v X ( t ) u X ( t ) p 1 + v X ( t ) p 1 ,
where C 1 ( t ) , C 2 ( t ) 0 for t + 0 and C 0 ( t ) , C 1 ( t ) , C 2 ( t ) C for all t [ 0 , ) .
For the proof, see [30].
Proposition A3.
Let us suppose that for any ( u 0 , u 1 ) , ( v 0 , v 1 ) A m , s 1 × A m , s 2 , the mapping N satisfies the following estimates:
N ( u , v ) X ( t ) C 0 ( t ) ( u 0 , u 1 ) A m 1 , s 1 + ( v 0 , v 1 ) A m 2 , s 2 + C 1 ( t ) ( u , v ) X ( t ) p + ( u , v ) X ( t ) q ,
N ( u , v ) N ( u ˜ , v ˜ ) X ( t ) C 2 ( t ) ( u , v ) ( u ˜ , v ˜ ) X ( t ) × ( u , v ) X ( t ) p 1 + ( u ˜ , v ˜ ) X ( t ) p 1 + ( u , v ) X ( t ) q 1 + ( u ˜ , v ˜ ) X ( t ) q 1 ,
where C 1 ( t ) , C 2 ( t ) 0 for t + 0 and C 0 ( t ) , C 1 ( t ) , C 2 ( t ) C for all t [ 0 , ) .
Then, N maps X ( t ) into itself and has one and only one fixed point ( u , v ) X ( t ) .
For the proof, see [26].

References

  1. Fujita, H.T. On the blowing-up of solutions of the Cauchy problem for tu=Δu+u1+λ. J. Fac. Sci. Univ. Tokyo Sect. 1966, 13, 109–124. [Google Scholar]
  2. Hayakawa, K. On the growing up problem for semi-linear heat equations. Proc. Jpn. Acad. 1973, 49, 503–505. [Google Scholar]
  3. Kobayashi, K.; Sirao, T.; Tanaka, H. The Critical Exponent(s) for the Semilinear Fractional Diffusive Equation. J. Math. Soc. Jpn. 1977, 29, 407–424. [Google Scholar]
  4. Strauss, W.A. Nonlinear scattering theory at low energy. J. Funct. Anal. 1981, 41, 209–409. [Google Scholar] [CrossRef] [Green Version]
  5. Glassey, R.T. Existence in the large for □u=F(u) in two space dimensions. Math. Z. 1981, 178, 233–261. [Google Scholar] [CrossRef]
  6. Glassey, R.T. Finite-time blow-up for solutions of nonlinear wave equations. Math. Z. 1981, 177, 323–340. [Google Scholar] [CrossRef]
  7. Schaeffer, J. The equation ttu+-Δu=|u|p for the critical value of p. Proc. R. Soc. Edinb. Sect. A 1985, 101, 31–44. [Google Scholar] [CrossRef]
  8. Yordanov, B.T.; Zhang, Q.S. Finite time blow up for critical wave equations in high dimensions. J. Funct. Anal. 2006, 231, 361–374. [Google Scholar] [CrossRef] [Green Version]
  9. Zhou, Y. Blow up of solutions to semilinear wave equations with critical exponent in high dimensions. Chin. Ann. Math. Ser. B 2007, 28, 205–212. [Google Scholar] [CrossRef] [Green Version]
  10. Kato, T. Blow-up of solutions of some nonlinear hyperbolic equations. Commun. Pure Appl. Math. 1980, 33, 501–505. [Google Scholar] [CrossRef]
  11. John, F. Blow-up of solutions of nonlinear wave equations in three space dimensions. Manuscr. Math. 1979, 28, 235–268. [Google Scholar] [CrossRef]
  12. Sideris, T.C. Nonexistence of global solutions to semilinear wave equations in high dimensions. J. Differ. Equ. 1984, 52, 378–406. [Google Scholar] [CrossRef] [Green Version]
  13. D’Abbicco, M.; Ebert, M.R.; Picon, T. Global Existence of Small Data Solutions to the Semilinear Fractional Wave Equation. In Trends in Mathematics, 1st ed.; Springer: Boston, MA, USA, 2017; pp. 465–471. [Google Scholar]
  14. D’Abbicco, M.; Ebert, M.R.; Picon, T.H. The Critical Exponent(s) for the Semilinear Fractional Diffusive Equation. J. Fourier Anal. Appl. 2019, 25, 696–731. [Google Scholar] [CrossRef] [Green Version]
  15. Kainane, M.A.; Reissig, M. Semi-linear fractional σ- evolution equations with mass or power non-linearity. Nonlinear Differ. Equ. Appl. 2018, 25, 42. [Google Scholar] [CrossRef]
  16. Kainane, M.A. Global Existence of Small Data Solutions to Semi-linear Fractional σ- Evolution Equations with Mass and Nonlinear Memory. Mediterr. J. Math. 2020, 17, 159. [Google Scholar] [CrossRef]
  17. Escobedo, M.; Herrero, A. Boundedness and blow up for a semilinear reaction-diffusion system. J. Differ. Equ. 1991, 89, 176–202. [Google Scholar] [CrossRef] [Green Version]
  18. Andreucci, D.; Herrero, M.A.; Velázquez, J.L. Liouville theorems and blow up behaviour in semilinear reaction diffusion systems. Ann. Inst. Poincaré Anal. Non Linéaire 1997, 14, 1–53. [Google Scholar] [CrossRef] [Green Version]
  19. Rencławowicz, J. Blow up, global existence and growth rate estimates in nonlinear parabolic systems. Colloq. Math. 2000, 86, 43–66. [Google Scholar] [CrossRef]
  20. Escobedo, M.; Levine, H.A. Critical blowup and global existence numbers for a weakly coupled system of reaction-diffusion equations. Arch. Rational Mech. Anal. 1995, 129, 47–100. [Google Scholar] [CrossRef]
  21. Snoussi, S.; Tayachi, S. Global existence, asymptotic behavior and self-similar solutions for a class of semilinear parabolic systems. Nonlinear Anal. 2002, 48, 13–35. [Google Scholar] [CrossRef]
  22. Sun, F.; Wang, M. Existence and nonexistence of global solutions for a non-linear hyperbolic system with damping. Nonlinear Anal. 2007, 66, 2889–2910. [Google Scholar] [CrossRef]
  23. Narazaki, T. Global solutions to the Cauchy problem for the weakly coupled of damped wave equations. Discrete Contin. Dyn. Syst. 2009, 2009, 592–601. [Google Scholar]
  24. Nishihara, K.; Wakasugi, Y. Critical exponant for the Cauchy problem to the weakly coupled wave system. Nonlinear Anal. 2014, 108, 249–259. [Google Scholar] [CrossRef]
  25. Mohammed Djaouti, A.; Reissig, M. Weakly Coupled Systems of Semilinear Effectively Damped Waves with Different Time-Dependent Coefficients in the Dissipation Terms and Different Power Nonlinearities; D’Abbicco, M., Ebert, M., Georgiev, V., Ozawa, T., Eds.; Birkhäuser: Basel, Switzerland, 2019; pp. 209–409, New Tools for Nonlinear PDEs and Application; Trends in Mathematics. [Google Scholar]
  26. Mohammed Djaouti, A. Semilinear Systems of Weakly Coupled Damped Waves. Ph.D. Thesis, TU Bergakademie Freiberg, Freiberg, Germany, 2018. [Google Scholar]
  27. Mohammed Djaouti, A. Modified different nonlinearities for weakly coupled systems of semilinear effectively damped waves with different time-dependent coefficients in the dissipation terms. Adv. Differ. Equ. 2021, 66. [Google Scholar] [CrossRef]
  28. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006; p. 204. [Google Scholar]
  29. Popov, A.Y.; Sedletskii, A.M. Distribution of roots of Mittag-Leffler functions. J. Math. Sci. 2013, 190, 209–409. [Google Scholar] [CrossRef]
  30. Mohammed Djaouti, A.; Reissig, M. Weakly coupled systems of semilinear effectively damped waves with time-dependent coefficient, different power nonlinearities and different regularity of the data. Nonlinear Anal. 2018, 175, 28–55. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mohammed Djaouti, A. Weakly Coupled System of Semi-Linear Fractional θ-Evolution Equations with Special Cauchy Conditions. Symmetry 2023, 15, 1341. https://doi.org/10.3390/sym15071341

AMA Style

Mohammed Djaouti A. Weakly Coupled System of Semi-Linear Fractional θ-Evolution Equations with Special Cauchy Conditions. Symmetry. 2023; 15(7):1341. https://doi.org/10.3390/sym15071341

Chicago/Turabian Style

Mohammed Djaouti, Abdelhamid. 2023. "Weakly Coupled System of Semi-Linear Fractional θ-Evolution Equations with Special Cauchy Conditions" Symmetry 15, no. 7: 1341. https://doi.org/10.3390/sym15071341

APA Style

Mohammed Djaouti, A. (2023). Weakly Coupled System of Semi-Linear Fractional θ-Evolution Equations with Special Cauchy Conditions. Symmetry, 15(7), 1341. https://doi.org/10.3390/sym15071341

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