Large Time Behavior for Inhomogeneous Damped Wave Equations with Nonlinear Memory

: We investigate the large time behavior for the inhomogeneous damped wave equation with nonlinear memory φ tt ( t , ω ) − ω ∈ R N imposing the condition ( φ ( 0, ω ) , φ t ( 0, ω )) = ( φ 0 ( ω ) , φ 1 ( ω )) in R N , where N ≥ 1, q > 1, 0 < ρ < 1, φ i ∈ L 1 loc ( R N ) , i = 0, 1, µ ∈ L 1 loc ( R N ) and µ (cid:54)≡ 0. Namely, it is shown that, if φ 0 , φ 1 ≥ 0, µ ∈ L 1 ( R N ) and (cid:90) R N µ ( ω ) d ω > 0, then for all q > 1, the considered problem has no global weak solution.


Introduction
We study the Cauchy problem of nonlinear damped wave equation in the form imposing the initial condition where N ≥ 1, 0 < ρ < 1, q > 1, φ i ∈ L 1 loc (R N ) with i = 0, 1, µ ∈ L 1 loc (R N ) and µ ≡ 0. The integral in the right-hand side of Equation (1) is known as "nonlinear memory".
We recall that the choices of both the domain and boundary conditions may influence significantly the properties and behavior of the physical system, which is mathematically represented by the above Equation (1). In general, the wave type equations are fundamental tools in recasting various propagation phenomena and in developing methods for numerical solving the physics problems. In detail, symmetries of wave type equations and their solutions have been pointed out and investigated by many contributors. We mention that the symmetry's properties were successfully used to obtain ortogonality's criteria for the existence of solutions in elastic and anisotropic media (see, for example, the pioneering papers of Love [1], Woodhouse [2], Chapman-Woodhouse [3], and the references therein). About the computational approach to the study of wave type equations, we recall that the nonlinear wave equations can be linked to linear wave equations, using symmetry transformations (in particular, non-local transformations). So, we can find a correspondence one-to-one between the solutions of nonlinear and linear wave equations. Resuming, each nonlinear wave equation can be linearized by a non-local symmetry analysis (for more details, see Bluman-Cheviakov [4]). For further discussion about the potential benefit of this procedure, we refer to the papers by Taylor-Kidder-Teukolsky [5] and Palacz [6] (spectral methods for propagation phenomena).
Here, we ask the question of whether the problem (1) and (2) admits global weak solutions. The interest in such a kind of results is motivated by a wide literature on what can be called "large time behavior of solutions" to wave type problems, which consists in providing sufficient criteria to the existence, nonexistence and blow-up of solutions to some classes of parabolic differential equations. For this purpose, we employ the test function method. Now, we recall some important results related to the blow-up of solutions to damped wave equations. First, we refer to the semilinear damped wave equation imposing (2). In the work of Todorova-Yordanov [7], can be found the following results: , then a unique global solution exists, under suitable initial values.
In the literature, the exponent "1 + 2 N " is known as the critical Fujita exponent. Indeed, it is critical for the problem (2) and (3), but it is also the critical exponent for the semilinear heat equation (see Fujita [8]) Further studies in Kirane-Qafsaoui [9] (semi-linear wave equation with linear damping) and Zhang [10] (nonlinear wave equation with damping), established that 1 + 2 N belongs to the nonexistence case.
It is worth pointing out that in the limit case ρ → 1 − , (1) reduces to The above equation was recently investigated by Jleli-Samet [11], who established the following results: , then there is no global weak solution to (4); • If N ≥ 3 and q ∈ (1 + 2 N−2 , ∞), then global solutions exist for suitable µ > 0.
This means that the critical exponent for Equation (4) is given by

Remark 1.
From the above result, one observes the considerable effect of the inhomogeneous term µ on the critical behavior of (2) and (3). Namely, for N ≥ 3, the critical exponent for (4) jumps from 1 + 2 N (the critical exponent for (2) and (3)) to the bigger exponent 1 + 2 N−2 . Notice that a similar phenomenon was observed for the heat equation [12] (Zhang, 1998) and the wave equation [13] (Zhang, 2000).
An interesting wave type problem is driven by the equation: assuming the initial condition (2). Problem (2) and (5) was first investigated by Fino [14], who proved the results as follows. Let Thus, we have: Moreover, a finite time blow-up occurs in the following cases (see again [14]): Finally, by considering compactly supported functions Studying the same problem (2) and (5), D'Abbicco [15] obtained a global existence result for q > q(N, ρ), where N ≤ 5 and ρ ∈ (0, 1), imposing suitable initial conditions.
The previous contributors give motivation to our work here. Indeed, we aim to study the effect of the inhomogeneous term µ on the large time behavior for problem (2) and (5).
Under sufficient conditions on the inhomogeneous term µ and the functions φ i ∈ L 1 loc (R N ), i = 0, 1, a nonexistence result is given in the following main result.
In the next Section 2 we collect the auxiliary mathematical tools which we will need in establishing the proof of Theorem 1 (see Section 3).

Preliminaries
We need some properties of fractional calculus to provide a definition of global weak solution to problem (1) and (2).
One of the tools of this study is the cut-off functions method. Here, we introduce a cut-off function In addition, for the sake of simplicity, fixing τ, L > 0, we introduce the functions: where δ, ς ≥ 2 are constants. Based on these functions, we construct a two-variable function ϑ as follows: Such a function ϑ satisfies the following result, whose proof is immediate.

Lemma 1.
For all τ, L > 0, the function ϑ defined by (10) belongs to V τ , where V τ is the test function space defined by (8).
We establish the following lemmata about the properties of function b L .

Lemma 2.
There exists a constant C > 0 such that Proof. Performing integration over the definition of b L , we have Acting with the definition of the cut-off function κ, we obtain Finally, taking z = L −1 ω, we deduce that which proves the desired result, where C is a constant equal to the integral in the right-hand side of the above equation.

Lemma 3.
Let ς > 2q q−1 . There exists a constant C > 0 such that Proof. Using the similar arguments as in the proof of Lemma 2 above, that is, performing integration over the definition of b L , and acting with the properties of the cut-off function κ, we obtain In order to manipulate the above equation, to get the aimed result, we perform some elementary calculations and have Invoking the properties of the cut-off function κ, we deduce that there exists a constant C κ > 0 such that Therefore, we deduce that Since ς > 2q q−1 , we derive easily the desired estimate.
And now, we have to consider the function a τ . The authors in a previous paper (see the proof of Theorem 1, p. 9, Equation (24), of [17]) obtained the following result: Here, we prove the following additional results: There exists a constant C > 0 such that Proof. By the definition of the function a τ , performing double derivation, we have Hence, using Lemma 4, we obtain Performing integration over [0, τ], the above equation gives us Taking s = t τ , it holds that which leads to the desired estimate, as δ > 2q+1−ρ q−1 .

Lemma 6.
Let δ > q+1−ρ q−1 . There exists a constant C > 0 such that Proof. Starting from the definition of a τ and performing derivation, we have Next, an application of Lemma 4 leads us to We integrate over [0, τ] to get Taking s = t τ , it holds that which yields the desired result (recall that δ > q+1−ρ q−1 ).

Lemma 7.
Let δ > 1−ρ q−1 . There exists a constant C > 0 such that Proof. Combining the definition of the function a τ with the estimate in Lemma 4, we obtain Integrating over the interval [0, τ], we deduce that Taking s = t τ , it follows easily that Since δ > 1−ρ q−1 , we have established our goal here.
The last lemma of this section provides an immediate property of the function a τ .

Proof of Theorem 1
For the sake of notational simplicity, in the sequel we denote by C any positive constant that is independent on τ. So, we give the complete proof of Theorem 1.

Proof.
We argue by contradiction. So, we suppose that φ ∈ L q loc (Ω) is a global weak solution to problem (1) and (2) (in the sense of Definition 1). For τ > 0, we define the set Then, by (9), for all τ > 0 and ϑ ∈ W τ , we have On the other hand, we use the -Young inequality with = 1 3 , to compute the following three estimates: and and Using (12)- (14) in (11), we get We recall that the inequality (15) holds for all τ > 0 and ϑ ∈ W τ . On the other hand, by Lemma 1, we know that for all τ, L > 0, the function ϑ defined by (10) belongs to V τ . Moreover, we point out that So, for all τ, L > 0, the function ϑ (defined by (10)) belongs to W τ . It follows that such a function ϑ is suitable as test function in (15), and (obviously) τ, L > 0. Moreover, we take δ and ς greater than two threshold values, that is and Now, we combine such "ingredients" to obtain suitable estimates of the three terms in the right-hand side of (15).
By the definition of the function b L , and using Lemma 8, we obtain Notice that since µ ∈ L 1 (R N ), by the properties of the cut-off function κ, and using the dominated convergence theorem, it holds that Since R N µ(ω) dω > 0, we deduce that there exists a constant L 0 > 0 such that Therefore, combining (22) and (23), we obtain On the other hand, by (10) we have Acting on the left-hand side of (21), it follows from (24) and (25) that Notice that Hence, fixing L > L 0 and passing to the limit as τ → ∞ in (26), we obtain R N µ(ω) dω ≤ 0, which contradicts the fact that R N µ(ω) dω > 0. We conclude that problem (1) and (2) admits no global weak solution.

Conclusions
The presented results confirm usefulness and simplicity of the test function method in analyzing different forms of wave equations. As a consequence of this approach in the previous sections, it is possible to see that the inhomogeneous damped wave equation with nonlinear memory (1) and (2) admits no global weak solution. In particular, we point out the absence of critical growth exponent q > 1 for the nonlinear memory, in proving such a result.
This topic may be significant for the study of the controllability of solutions to certain nonlinear models of physics systems, together with the symmetry analysis. We have already mentioned in Section 1 above, a possible relationship between the boundary conditions, and the properties and behavior of a physical system. Here we imposed a Cauchy condition, but it will be interesting to know how Neumann, Robin and mixed boundary conditions affect the analysis of the wave type problem (1) and (2).
Author Contributions: Investigation, all authors; Writing-original draft, all authors. All authors have read and agreed to the published version of the manuscript.
Funding: The first author is supported by Researchers Supporting Project number (RSP-2020/57), King Saud University, Riyadh, Saudi Arabia.