Global existence and long-time behavior of solutions for fully nonlocal Boussinesq equations

: In this paper, we study initial boundary value problems for the following fully nonlocal Boussinesq equation


Introduction
Let Q T = Ω × (0, T ), B T = ∂Ω × (0, T ), where Ω ⊂ R N is a smooth bounded domain and T > 0. In this paper, we research initial boundary value problems for the following fully nonlocal Boussinesq in Ω, where 1 < β < 2, 0 < σ < 1 and N < 4σ. Moreover, f (u) is a given nonlinear function, and ϕ(x), φ(x) are initial value data. The Caputo fractional operator C 0 D β t is defined by where σ is the gamma function. It is worth noting that C 0 D β t u may turn into the usual derivative u tt when β → 2, see [1] for details. The fractional Laplacian operator (−∆) σ can be defined via spectral decomposition where µ k and w k , k ∈ N are eigenpairs of the following eigenvalue problem and v k = Ω v(x)w k dx, with ∥w k ∥ L 2 (Ω) = 1.
Therefore, it's called the spectral fraction Laplace operator, see [2,3] for details. The equation (1.1) is nonlocal both in space and time, so we call such a Boussinesq equation a fully nonlocal Boussinesq equation. Problem (1.1)'s widespread use as a model for anomalous diffusion in physical field serves as a significant incentive for study. Time fractional derivatives are generally exploited to model the omnipresent memory effects such as anomalous diffusion, wave propagations and neuronal transmission in Purkinje cells, etc. For example, in [4], the authors demonstrated how Caputo time fractional derivatives can be used to analyze turbulent eddies' trapping effects. In fact, β order time fractional derivatives have been used for "superdiffusion"-in which particles spread quickly against the laws of Brownian motion. Nevertheless, time fractional derivatives and "anomalous subdiffusion" are frequently linked when β ∈ (0, 1), see [5][6][7]. Furthermore, space fractional derivatives can be used to describe nonlocal effects, such as anomalous diffusion and Lévy processes. Recently, time or space fractional wave equations have drawn a lot of interest, see [8][9][10][11][12][13][14][15][16] for examples.
In 1872, J. Boussinesq [17] presented the Boussinesq equation which can illustrate how small amplitude long waves propagate on the surface of shallow water. The improved Boussinesq equation (IBq equation) may be written by which can describe the continuum limit of shallow water waves in a one-dimensional nonlinear lattice and other modes supporting linear waves with a negative dispersion. They can also explain the lowestorder nonlinear effects in the evolution of perturbations using a dispersion relation similar to that for sound waves (in terms of wave amplitudes). In [18], it was indicated that the IBq equation may be deduced from starting with the accurate hydrodynamical set of equations in plasma and modifying the IBq equation in a manner similar to modifying the Korteweg-de Vries equation to derive which is called IMBq (modified IBq) equation. During these years, the theory of Boussinesq equations has been developed significantly, see [19][20][21][22][23][24][25]. In [19,20], Wang and Chen studied Cauchy problems for the following generalized Boussinesq equation They discussed whether or not global solutions exist. Moreover, by Banach fixed point theorm, they obtained that small-amplitude solutions exist globally. In [21], the authors researched Cauchy problems for the following Boussinesq equation By using Banach fixed point theorem, they proved that the solution exist locally. Under different dimensions, they obtained global existence of smooth solutions using potential well method. In addition, they acquired the blow-up of solutions. In general, in [23] and [24], the authors studied Cauchy problems for the generalized Boussinesq equation with damping terms, respectively. By Banach fixed point theorem, they constructed a class of time-weighted Sobolev spaces, and obtained global existence and long-time behavior of small amplitude solutions. In [26], Li, Yan and Xie studied an extended (3 + 1)-dimensional B-type Kadomtsev-Petviashvili-Boussinesq equation, and obtained a family of rational solution through its bilinear form and symbolic computation. In addition, there are many results on nonlocal nonlinear problems, see [27,28] for examples. Recently, fully nonlocal evolution equations has received a lot of attention. In [29], Kemppainen, Siljander, Zacher studied classical solutions and large-time behavior for fully nonlocal diffusion equations. In [30], Li, Liu and Wang researched Cauchy problems for Keller-Segel type fully nonlocal diffusion equation. Therefore, the study of fully nonlocal Boussinesq equations has certain theoretical significance. Comparing the equation (1.1) and (1.2), we just replace u tt with C 0 D β t u and −∆ with (−∆) σ , as it comes to nonlocality and memory effect. In light of these works mentioned above, we aim to investigate Cauchy problems for the fully nonlocal Boussinesq equation in (1.1) and generalize their results in [19,20] by Wang and Chen. Nevertheless, the spectral fractional Laplacian operator (−∆) σ makes no sense in R N , because the Laplacian operator's spectrum in R N is continuous purely. As a result, we research the initial boundary value problem (1.1).
Nevertheless, as yet, there are few results on global existence and long-time behavior of solutions for problem (1.1). In reality, the corresponding linear problem has not received much attention. The major difficulty is the nonlocality and nonlinearity of (−∆) σ f (u). In addition, since the memory effect of ∂ β t u, the definition of weak solutions is difficult to introduce, and potential well method may be also ineffective for fully nonlocal Boussinesq equations. Inspired by [16], we first study the corresponding linear Boussinesq equation to obtain explicit expressions of Green operators. Further, we establish some rigorous estimates of the Green operators to acquire global existence and decay properties of weak solutions for linear problems. Second, we establish new chain and Leibniz rules concerning the spectral fractional derivatives. Based on these given results, under different dimensions N and small initial value condition, by Banach fixed point theorm, we obtain global existence and longtime behavior of weak solutions for problem (1.1) in the time-weighted fractional Sobolev spaces. Throughout this paper, we replace ∥ · ∥ H s (Ω) with ∥ · ∥ s , and the notation C ≲ D means that there is a constant M > 0 such that C ≤ MD.
The following are main results of this manuscript.
Suppose that all assumptions in Theorem 1.1 hold, and Let u be the global weak solution of problem (1.1), then there holds The paper is organized as follows. In Section 2, we introduce fractional Sobolev space briefly, and give several properties of Mittag-Leffler functions. In Section 3, we study global existence and decay estimates of weak solutions for the corresponding linear Boussinesq equation. In Section 4, for small initial values condition, we establish global existence and long-time behavior of weak solutions for problem (1.1) under different dimensions N.

Linear estimates
In this section, we obtain explicit expressions and some rigorous estimates of the Green operators for problem (1.1). First, we study the corresponding linear problem in Ω, where the function h(t, x) is given. Inspired by [12,13,16], we try to find the solution of problem (3.1) as follows Therefore, it can be inferred that Then, we getũ Using the inverse Laplace transform, it is derived from Lemma 2.1 in [33] that In terms of (3.4), if (3.2) converges, then we may formally obtain the following weak solution where the Green operators are Then, we get from Lemma 2.2 Definition 3.1. We say that u is a weak solution of problem We can get the following estimations on the Green operators immediately from Lemma 2.1 and 2.2.
Therefore, u ∈ C([0, T ]; H s (Ω)) and u t ∈ C([0, T ]; H s (Ω)). □ Proposition 3.6. Suppose that h ∈ L ∞ (0, T ; H s (Ω)), ϕ ∈ H s (Ω), and φ ∈ H s (Ω). If u is a global weak solution of problem (3.1), then we have Remark 3.7. It is simple to observe that Proposition 3.6 is more rigorous than Proposition 3.5, which initially appears to be weak, when it comes to estimates of weak solutions. Proposition 3.6 cannot be ignored because Proposition 3.5 can be used to establish local existence theorems to problem (1.1), but Proposition 3.6 cannot.

Proofs of main results
In this section, by constructing time-weighted fractional Sobolev spaces and Banach fixed point theorem, we get global existence and long-time behavior of weak solutions for problem (1.1). Now, we provide the definition of weak solutions for problem (1.1).
Moreover, if T > 0 can be chosen as any positive number, u is called a global weak solution for problem (1.1).
First, we need the following lemmas to establish chain and Leibnitz rules concerning (−∆) σ .
Since Ω is regular enough, it has the so-called extension property, namely: For any s ∈ (0, 1), there exists an extensionȗ of u ∈ H s (Ω) such thatȗ ∈ H s (R N ) andȗ| Ω = u where Moreover, ∥ȗ∥ H s (R N ) ≤ C∥u∥ H s (Ω) . In particular, taking such extension is the trivial one, namely the extension by zero outside Ω, there holds and A ≍ B represents that there are two constants c 1 , c 2 > 0 satisfying c 1 A ≤ B ≤ c 2 A, see [3,35] for details. Moreover, from [3], the space H s (Ω) is redefined by and these two norms ∥ · ∥ H s (Ω) and ∥ · ∥ H s (Ω) on H s (Ω) are equivalent. Therefore, we conclude that Further, let s ∈ [1, 2) and u ∈ H s (Ω). Taking s = 1 + δ, we derive Therefore, we obtain ∥ȗ∥ H s (R N ) ≍ ∥u∥ H s (Ω) .
Therefore, we obtain the following chain and Leibniz rules concerning the spectral fractional derivatives, which plays an fundamental role on existence of weak solutions.
Proof. By Lemma 4.2, we have Suppose that s ∈ (0, 2) and u, v ∈ H s (Ω), then there holds Proof. By Lemma 4.3, we have □ Now, we prove main results of this manuscript.
Proof of Theorem 1.1. Define and ρ(v, w) = ∥v − w∥ V for any v, w ∈ V. Consequently, it is evident that the metric space (V, ρ) is complete. Moreover, The operator P on V is defined as By Sobolev embedding theorem, it is easy to get P(u) ∈ L ∞ (0, T ; H s (Ω)) from Lemma 3.2. For any u ∈ V, using Proposition 3.6, Lemma 4.3 and Sobolev embedding theorem, we obtain Hence, when ε and ζ are small enough, we get Thus, we acquire P(u) ∈ V. Next, we prove that P : V → V is contractive. Taking any u, v ∈ V, by Sobolev embedding theorem, we obtain from Lemma 4.4 for some ϑ ∈ (0, 1). Then, ρ(P(u), P(v)) ≲ ζ q−1 ρ(u, v).
Taking ζ small enough, we conclude that P is contractive. Using Banach fixed point theorem, we derive that P has a unique fixed point u ∈ V.
Remarkably, we have get the solution u ∈ C((0, T ]; H s (Ω)) and sup 0<t≤T t αβ ∥u(t)∥ s ≤ ζ. In reality, we just have to verify that there exist T 0 > 0 small enough such that problem (1.1) admits a weak solution in Using Lemma 3.2 and Sobolev embedding theorem, we get P(u) ∈ C([0, T ]; H s (Ω)). For any u ∈ X, using Lemma 3.2, Lemma 4.3 and Sobolev embedding theorem, we acquire Then, when ε and ζ are small enough, we derive ∥P(u)∥ X ≤ 2ζ.

Conclusions
In this paper, we study initial boundary value problems for fully nonlocal Boussinesq equations. We overcome full nonlocal effects generated by C 0 D β t and (−∆) σ , and obtain some new results as follows: (a) We obtain explicit expressions and some rigorous estimates of the Green operators for the corresponding linear equation;(b) We establish new chain and Leibnitz rules concerning (−∆) σ ;(c) We establish time-wighted fractional Sobolev spaces and obtain global existence and long-time behavior of weak solutions. What's more, our work adds some novelty results to the subject of Boussinesq equations, which may provide a certain theoretical support for the study of fully nonlocal wave equations and have certain theoretical significance.