On a time-space fractional di ff usion equation with a semilinear source of exponential type

: In the current paper, we are concerned with the existence and uniqueness of mild solutions to a Cauchy problem involving a time-space fractional di ff usion equation with an exponential semilin-ear source. By using the iteration method and some L p − L q -type estimates of fundamental solutions associated with the Mittag-Le ffl er function, we study the well-posedness of the problem in two di ff erent cases corresponding to two assumptions on the Cauchy data. On the one hand, when considering initial data in L p ( R N ) ∩ L ∞ ( R N ), the problem possesses a local-in-time solution. On the other hand, we obtain a global existence result for a mild solution with small data in an Orlicz space.


Introduction
In this paper, we consider the following initial value problem for a time-space fractional diffusion equation where N ≥ 1, 0 < σ ≤ 2, h is the initial data function, and the symbol ∂ α t stands for the Caputo derivative of fractional order α ∈ (0, 1) (Section 2).In Problem (1.1), we are mainly focus on the semilinear case in which the function F satisfies the following assumptions and where ν, p > 1 and L is a positive constant.The reason why we study this source function comes from the great interest of the PDEs community with a polynomial source of the form G p (u) = |u| p−1 u or G p (u) = u p and some similar forms.Many good papers about this topic have attracted our attention.Wang-Xu [1] and Xu-Su [2] used the potential well method to investigate the well-posedness of a pseudo-parabolic equation with nonlinear function G p .Lian et al. [3] studied a Schrödinger equation with polynomial nonlinearity.They used infinite Nehari manifolds with geometric features to provide infinite sharp conditions for global existence and blowup results of solutions.A modified form of G p was considered by Chen et al. [4] in the Gierer-Meinhardt system.The authors applied a functional method to obtain a bound of some ratios of the solution, and then, the existence of global and blowup solutions were proved.
and the value q c = N(p−1) 2 called the critical exponent plays an important role in investigating the existence and uniqueness results.Considering the case when q = p = q c = N N−2 in R 2 , the power exponent is approximated ∞.Therefore, it seems reasonable to replace the source function in (1.4) with the nonlinearity of the exponential type.As far as we know, the attention on nonlinear functions satisfying the assumptions (1.2) and (1.3) is derived and developed by several works in the literature [5][6][7][8][9].Two groups in [5,8] and [7,9] studied, respectively, parabolic equations and the Schrödinger equation (NLS) with exponential nonlinearities.More precisely, in [8], Ioku proved the global-in-time existence of a mild solution to a semilinear heat equation with exponential nonlinearity under some smallness assumptions on the initial data.Meanwhile, Furioli [5] showed that the notions of weak and mild solutions are equivalent and investigated decay estimates and the asymptotic behavior of small-data global solutions.Nakamura and Ozawa in [9] provided global-in-time results for solutions in homogeneous Sobolev spaces and homogeneous Besov space to a NLS with a source function of exponential type.A source function f (u) = (e 4π|u| 2 − 1 − 4π|u| 2 ) was considered for a two-dimensional NLS problem by Ibrahim et al. [7].The authors showed that the solution to their problem tends to be a free Schrödinger under certain conditions.Also, we refer to some other works by Nakamura and Ozawa [10] and Ibrahim et al. [6] for wave equations with the nonlinearity of exponential growth.
During the past decades, fractional calculus has received increased attention due to its wide applications in diverse fields of science and engineering such as stochastic processes [11], fluid mechanics [12,13], chemotaxis in biology [14], viscoelasticity [15], etc. Apart from that, many interesting mathematical models and results on this topic have been done [16][17][18][19][20][21][22][23][24][25].A strong inspiration for studying Problem (1.1) with the presence of the operator ∂ α t comes from the fact that many physical phenomena carry the substance of diffusion processes, many studies about diffusion equations have been done [26][27][28][29][30], and fractional calculus is very effective in modeling anomalous diffusion processes.In fact, while a diffusive particle in the usual diffusion process possesses the mean square displacement behaving like C 1 t for t → ∞, such behavior of a particle in an anomalous diffusion process is C 2 t α [31], C 1 , C 2 are positive constants.Starting from the above characteristics of an anomalous diffusion process, many good works about fractional diffusion equations have been done.Because it is very difficult and lengthy to present all the related works, we would like to present only the works, which motivated us.In [32], a general time-fractional diffusion equation subject to the Dirichlet boundary condition was studied by Vergara and Zacher.By using energy estimates and a powerful inequality for integrodifferential operators, the authors proved sharp estimates for the decay in time of solutions.Andrade et al. [33] considered the following non-local initial boundary value problem associated with a time-fractional heat equation where where p > 1 and c is a positive constant.The existence and regularity of mild solutions were established with some sufficient conditions.In [34], Tuan et al. were concerned in a terminal value problem for a time-space fractional diffusion equation.For the problem with linear source function, regularity properties of solutions were studied.The existence, uniqueness and regularity to solution were proved in the case of nonlinear source.
The main results of this paper are providing the existence and uniqueness of mild solutions to Problem (1.1) with function F satisfying (1.2) and (1.3).Corresponding to two different cases of initial data, we obtain a local-in-time solution and a global small-data solution.With usual initial data, by the Picard iteration method and some L p − L q or L p − L Ξ estimates of fundamental solutions involving the Mittag-Leffler function, we provide the existence and uniqueness of mild solutions to (1.1) on a reasonable time interval (0, T ].Apart from that, we also show that solutions are continuous from (0, T ] to L p (R N ).Global-in-time results are obtained by making use of the norm (3.14).From the technical point of view, we split the second term of the right-hand side of (2.1) into two parts, while the part with small-time is easy to handle, the controlling of the large time part requires small assumptions on the initial data in an Orlicz space to be achieved.
The structure of the paper is as follows.Firstly, we provided some preliminaries in section 2 including some function spaces, fractional settings and formula, linear estimates for a mild solution.The main results about the local and global well-posedness are stated in section 3.

Some basic setups
We first introduce some function spaces.Let (B, ∥•∥ B ) be a Banach space.For T > 0, we denote by C ([0, T ]; B) the space of all continuous functions u from [0, T ] to B and define the following space Recall that L ∞ (0, T ; B) is a Banach space with respect to the norm ∥u∥ L ∞ (0,T ;B) := ess sup t∈(0,T ) Let Ξ(z) := e z p − 1.The Orlicz space L Ξ (R N ) is defined as the space of all functions satisfying the following converging result for some κ > 0 The space L Ξ (R N ) is a Banach space with the Luxemburg norm given as follows In fact, we have the following embeddings For more details about Orlicz spaces, we refer the reader to [37,Chapter 8] and references given there.
Next, let us provide some fractional settings.For a, b > 0, the Beta function B and the Gamma function Γ are defined respectively as follows Let α ∈ (0, 1).By considering the following memory kernel for a smooth enough function u, we can define the Caputo derivative of order α by Also, for two real constants α 1 and α 2 , we define the Mittag-Leffler function E α 1 ,α 2 : C → C in the following way .

Fundamental solutions
The fractional Laplace operator can be defined via the Fourier multiplier [5, Section 2] (−∆) where the Fourier transform is recalled as follows and F −1 is the inverse Fourier transform.From the above definitions, for any α ∈ (0, 1) and σ ∈ (1, 2], we define two functions A 1 (α, σ)(t) and A 2 (α, σ)(t) by From [35, Section 1.3] or [14, Section 2], we see that the mild solution to Problem (1.1) satisfies the following Duhamel integral equality It turns out that handling the operators A 1 (α, σ)(t) and A 2 (α, σ)(t) plays an important role in controlling norms of the mild solution u.Therefore, we provide some linear estimates of A 1 (α, σ)(t) and A 2 (α, σ)(t) in the following lemma.
Then, there exists a positive constant C 0 such that the following statements hold Corollary 2.2.For any u ∈ L Ξ (R N ), we have Electronic Research Archive Volume 30, Issue 4, 1354-1373.
Proof of Corollary 2.2.First, for a positive real number a, we have the following observation inf κ > 0 Therefore, the following equality holds Next, from the expansion of the exponential function, for any κ > 0, we have Combining this result and (2.4) yields the desired estimate.Similarly, we can find the estimate of the norm for A 2 (α, σ)(t)u.The proof is completed.□

Existence and uniqueness
This section is used to present the main results of this paper including the existence and uniqueness of mild solutions satisfying Eq (2.1).We provide two different results about local well-posedness and global well-posedness according to two cases of initial data.
1) possesses a unique mild solution u on [0, T ] where T is small enough.In addition, this solution is also continuous on (0, T ]. • By making some small assumptions on the initial data in L Ξ (R N ), we can prove that the solution u to Problem (1.1) exists globally in time.
Then, we can find a reasonable number T such that Problem (1.1) possesses a unique mild solution Proof.To begin, we consider the sequence {u l } l∈N as follows Electronic Research Archive Volume 30, Issue 4, 1354-1373.
We aim to prove that {u l } l∈N is a Cauchy sequence in L ∞ (0, T ; L p (R N ) ∩ L ∞ (R N )).Then, the completeness of this space ensures the existence of a limit function u that can be shown to be the unique mild solution to Problem (1.1).To this end, the first task is to check whether {u l } l∈N is in or not.Indeed, we apply Lemma 2.1 to get and This result implies that . Before moving to the second step, we provide some nonlinear estimates of the source function.For functions w, v ∈ L p (R N ) ∩ L ∞ (R N ), we find that and We are now ready to consider the remaining elements of for any l ∈ N. From Lemma 2.1, the following estimate is satisfied for any t > 0 Apply (3.2) and the assumption (1.3), we find that By similar arguments, we also get a same result for the L p -norm as follows where we have used (3.3) with w = u l and (1.3).
Combining the above two estimates and choosing , we obtain the following result In view of (3.1) and (3.6), for any l > 1, if u l ∈ B(0, R 1 ), we obtain the estimate below ess sup From (3.1) and (3.7), the induction method can be applied to conclude that {u l } l∈N ⊂ B(0, R 1 ).
In addition, we can check that {u l } l∈N is a Cauchy sequence in B(0, R 1 ).In fact, presume for l ≥ 2 that u l and u l−1 are elements of B(0, R 1 ), the techniques as in (3.4) and (3.5) enable us to find for any t ∈ (0, T ) that Therefore, if we choose , the following estimate can be drawn from (3.8) and (3.9) , for any l ≥ 2. Based on this result, for any l 2 > l 1 ≥ 2, we have .
It means {u l } l∈N is a Cauchy sequence in B(0, R 1 ), provided that we have already shown that {u l } l∈N ⊂ B(0, R 1 ).By the completeness of the space L ∞ 0, T ; L p (R N ) ∩ L ∞ (R N ) and the dominated convergence theorem, there exists a unique limit function u satisfying In addition, we can also show that u ∈ C (0, T ]; L Ξ 0 (R N ) .For t, ε > 0, it is easy to check that where we define Electronic Research Archive Volume 30, Issue 4, 1354-1373.
By Theorem 3.2 and Remark 1.6 in [36], we deduce Furthermore, by using Lemma 2.1, we obtain .
From the fact that u ∈ B(0, R 1 ) In sum, we have We next consider the third term on the right hand side of (3.10) as follows where we apply Lemma 2.1.Therefore, there holds lim [Global small-data solution] Let ν > 4 3 and p ≥ 2. Suppose that one of the following assumptions is satisfied, , and there exists a constant q ≥ 2 satisfying and there exists a constant q ≥ 2 satisfying Remark 3.1.It's not too difficult to find a non-empty set of parameters meeting the assumptions of Theorem 3.2.Indeed, it can be pointed out some examples as follows (i) if α = 0.4, ν = 1.66, σ = 1.5 and p = 2 and N = 2, we can choose q = 2.1; (ii) if α = 0.7, ν = 2.5, σ = 2 and p = 3 and N = 3, we can choose q = 4.

Remark 3.2. By the embedding L
Proof.We first introduce a function space for the existence of solutions as follows Next, we consider the estimate of the source term.Suppose that w, v ∈ L Ξ (R N ), we can deduce from Taylor's expansion of the exponential function and Hölder's inequality with 1 q = 1 3q + 1 3q + 1 3q that Thanks to the definition of the Luxemburg norm and the monotone convergence theorem, for any u ∈ L Ξ (R N ) there holds provided that z q Γ(q + 1) < e z − 1 for any q > 1 and z > 0. Therefore, we obtain the following estimate Applying (3.16) to (3.15), we get immediately that Using [5, Lemma 3.3], we derive where we use the fact that for j ∈ N, there holds Γ( j + 1) = j! and denote , where C 1 is a positive constant that is independent of w, v. Let R 2 be a sufficiently small constant.
for any l ∈ N, we can show that u l+1 ∈ B(0, R 2 ).In fact, on the one hand, by applying Lemma 2.1, we derive dm.
Then, we use (3.17) with w = u l and (1.3) to find where we presume that R 2 < (3q) − 1 p .On the other hand, repeat application of Lemma (2.1) with s = p and r = pq p+q yields where we note that if N < σq < 2σq, we deduce p < Nr N−2rσ .As a consequence, if R 2 < (3r) − 1 p , it follows from the above estimate and (3.17) that provided that max{3r(ν−1), 3r} ≥ p. Combining (3.18), (3.19) and the embedding where C 2 is a positive constant coming from the embedding.According to the definition of the Beta function, we have In view of (3.20) and (3.21), for any t > 0, the following estimate is satisfied provided that we get immediately for all t > 0.
Next, we set Then, for any t ≤ T , (3.20) implies At the same time, if t > T , we deduce u l+1 (t) − u 1 (t) where we apply (3.22).If R 2 satisfies there holds immediately u l+1 (t) − u 1 (t) , for all t > T .
From the above results, whether t is greater than T or t is less than T , we always get the following result as long as u l ∈ B(0, R 2 ) and R 2 is small enough.For the purpose of proving {u l } l∈N is a subset of B(0, R 2 ), we also need to consider the initial data h.On the one hand, by using Corollary 2.2, we get easily that On the other hand, Lemma 2.1 shows that where η = αN σβ > 1.Then, we get Presume that the initial data is small enough, precisely, Then, we can conclude that u 1 ∈ B(0, R 2 ).Hence, if u l ∈ B(0, R 2 ) for any l ≥ 2, we have u l+1 ∈ B(0, R 2 ), provided that R 2 and the initial data are sufficiently small.Summarily, we have {u l } l∈N ∈ B(0, R 2 ).
To complete the Banach principle argument, we need also to show that {u l } l∈N is a Cauchy sequence in B(0, R 2 ).Since the techniques are not too different from those in the results above, we only briefly present the main estimates.For u l and u l−1 in B(0, R 2 ), l ≥ 2, Lemma 2.1 yields It follows immediately that We then find that ess sup t∈(0,T ) ) .
On the one hand, for any t ≤ T and l ≥ 2, we derive (3.23) that u l+1 (t) − u l (t) .
On the other hand, if t > T , the following estimate holds u l+1 (t) − u l (t) .
Then, by similar arguments of Theorem 3.1, we can show that {u l } l∈N is a Cauchy sequence.In addition, since L ∞ β 0, T ; L Ξ (R N ) is a complete space with the metric , there exists a unique limit function of {u l } l∈N in B(0, R 2 ), which is the unique mild solution to Problem (1.1).The proof is completed.□

Conclusions
This paper considers a Cauchy Problem for a time-space fractional diffusion equation with exponential source term.By iteration method, a unique local mild solution is derived for initial data in L p (R N ) ∩ L ∞ (R N ).The existence and uniqueness are extended to be global in time when we suppose additionally that initial data in an Orlicz space are small enough.However, since the space C ∞ 0 (R N ) is not dense in L Ξ (R N ), the continuity of solutions in the term of time-variable is not considered for the global case.This will be a potential approach to improve the results of this work in the future.