Stabilization of a viscoelastic wave equation with boundary damping and variable exponents: Theoretical and numerical study

: In this work, we consider a viscoelastic wave equation with boundary damping and variable exponents source term. The damping terms and variable exponents are localized on a portion of the boundary. We ﬁrst, prove the existence of global solutions and then we establish optimal and general decay estimates depending on the relaxation function and the nature of the variable exponent nonlinearity. Finally, we run two numerical tests to demonstrate our theoretical decay results. This study generalizes and enhances existing literature results, and the acquired results are thus of signiﬁcant importance when compared to previous literature results with constant or variable exponents in the domain.


Introduction
The importance of partial differential equations in comprehending and explaining physical interpretation of problems that arise in numerous fields and engineering motivates many researchers to analyze and investigate the existence and stability of their solutions. Hyperbolic partial differential equations are the most interesting kind of partial differential equations, since they are utilized to simulate a wide and important collection of phenomena, such as aerodynamic flows, fluid and contaminant flows through porous media, atmospheric flows, and so on. Of the higher order hyperbolic equations, the wave equation is the most obvious. Klein-Gordon, Telegraph, sine-Gordon, Van der Pol, dissipative nonlinear wave and others are well-known hyperbolic equations that are important in the fields of wave propagation [1], random walk theory [2], signal analysis [3], relativistic quantum mechanics, dislocations in crystals and field theory [4], quantum field theory, solid-state physics, nonlinear optics [5], mathematical physics [6], solitons and condensed matter physics [7], interaction of solitons in collision-less plasma [8], fluxions propagation in Josephson junctions between two superconductors [9], motion of a rigid pendulum coupled to a stretched wire [10], material sciences [11] and non-uniform transmission lines [12] are some of the topics covered. For more related results, we refer to [13][14][15]. There is a vast range of publications for numerical solutions of hyperbolic partial differential equations, such as the one in [16][17][18][19][20][21]. In recent years, great efforts have been devoted to study problems with nonlinear dampings and source terms, and several existence, decay and blow up results have been established. Georgiev and Todorova [22] considered the following nonlinear problem on Ω × (0, T ) u = 0, on ∂Ω × (0, T ) u(x, 0) = u 0 (x), u t (x, 0) = u 1 (x), on Ω, where the damping term h(u t ) = |u t | m−2 u t and the source term F(u) = |u| q−2 u are localized on the domain and established global existence when q ≤ m and a blow up result when q > m. This work was improved by Levine and Serrin [23] to the case of negative energy and m > 1. For problems with boundary damping and source terms, we mention the work of Vitillaro [24] who considered the following problem on Ω × (0, T ) ∂u ∂n + h(u t ) = F(u), on Γ 1 × (0, T ) u(x, 0) = u 0 (x), u t (x, 0) = u 1 (x), on Ω.
(1. 2) where the damping term h(u t ) = |u t | m−2 u t and the source term F(u) = |u| q−2 u are localized on a part of the boundary. The author established local existence and global existence of the solutions under some suitable conditions on the initial data and the exponents. In the presence of the viscoelastic term, Cavalcanti et al. [25] discussed the following problem on Ω, In this work, a global existence result for strong and weak solutions was established and some uniform decay rates were proved under some assumptions on g and h. Al-Gharabli et al. [26] established general and optimal decay result for the same problem (1.3) considered in [25] where the relaxation g satisfies more general conditions than the one in [25]. For more results in this direction, we refer to [27][28][29][30][31]. In particular, Liu and Yu [31] investigated the following problem on Ω, (1.4) where the damping term h(u t ) = |u t | m−2 u t and the source term F(u) = |u| q−2 u are localized on a part of the boundary, and established several decay and blow up results under some suitable conditions on the initial data, the relaxation function and the exponents. Notice here that both the damping and source terms in [31] are localized on a part of the boundary, although, they are of constant nonlinearity. Moreover, the relaxation function g satisfies the condition where ξ is a positive differentiable function. In fact, Liu and Yu [31] used the Multiplier method for stability and the potential well technique to prove the existence of the global solution. Moreover, the authors established a general decay when m ≥ 2 and an exponential decay m = 2.
Many new real-world problems, such as electro-rheological fluid flows, fluids with temperature-dependent viscosity, filtration processes through porous media, image processing, hemorheological fluids, and others, came as a result of advances in science and technology, such as those problems which required modeling with non-standard mathematical functional spaces. The Lebesgue and Sobolev spaces with variable exponents [32][33][34][35] have shown to be very important and user-friendly tools to tackle such models. PDEs with variable exponents have recently attracted a lot of attention from researchers and academics. However, the majority of the findings for hyperbolic issues with variable exponents dealt with blow-up and non-global existence. On the stability of nonlinear damped wave equations with variable exponent nonlinearities, we only have a few results. It is worth mentioning the work of Messaoudi et al. [36], who explored the stability of the following equation where m(·) ≥ r(·) ≥ 2. The authors in their work showed that the solution energy decays exponentially if m ≡ 2 and when m 2 = esssup x∈Ω m(x) > 2, they obtained a polynomial decay at the rate of t 2/(m 2 −2) . Also, Ghegal et al. [37] established a stability result similar to that of [36] for the equation and proved under appropriate conditions on m(·), q(·), and the initial data, a global existence result. Messaoudi et al. [38] recently looked at the following problem and used the well-depth approach to verify global existence and provide explicit and general decay results under a very general relaxation function assumption.
In our present work, we are concerned with the following problem on Ω, on a bounded domain Ω ⊆ R n with a smooth boundary ∂Ω = Γ 0 ∪ Γ 1 where Γ 0 and Γ 1 are closed and disjoint and meas.(Γ 0 ) > 0. The vector n is the unit outer normal to ∂Ω. The function g is a relaxation function and u 0 and u 1 are given data. The functions m(·) and q(·) are the variable exponents. System 1.6 describes the spread of strain waves in a viscoelastic configuration. We first prove a global existence result for the solutions of problem (1.6) by using the potential well theory. Then we establish explicit and general decay results of problem (1.6) for a larger class of relaxation functions (see Assumption A1 below). To back up our theoretical decay results, we provide two numerical tests. Our decay results extend and improve some earlier results such as the one of Cavalcanti et al. [25], Al-Gharabli et al. [26], Liu and Yu [31] and the one of Messaoudi et al. [39]. In our work, we apply the energy approach (Multiplier Method), combined with various differential and integral inequalities equipped with the Lebesgue and Sobolev spaces with variable exponents. The multiplier method relies mostly on the construction of an appropriate Lyapunov functional L equivalent to the energy of the solution E. By equivalence L ∼ E, we mean for two positive constants α 1 and α 2 . To prove the exponential stability, we show that L satisfies for some c 1 > 0. A simple integration of (1.8) over (0, t) together with (1.7) gives the desired exponential stability result. In the case of a general decay result, we prove that L satisfies a differential inequality that combines the relaxation function and the other terms coming from the nonlinearites. Then we use some properties of the convex functions and other mathematical arguments to obtain general decay estimates depending on the relaxation function and the nature of the variable exponent nonlinearity. In fact, the Multiplier Method proved to be efficient in tackling such problems with dissipative terms either on the domain or in a part of the boundary. In the present paper, some properties of the convex functions are exploited. We also use the well-depth method to establish the global existence of the solutions. We show that the methods and tools used in this paper are sufficient to handle our problem and are less complicated than other methods which guide us to our target.

Related results to our problem
• Cavalcanti et al. [25] and Al-Gharabli et al. [26] investigated the same problem. However, in [25] and [26], the nonlinear damping term is h(u t ) which satisfies some specific conditions. In our case h(u t ) = |u t | m(·)−2 u t where m(·) is a function of x where x is in a part of the boundary which makes our problem more complicated especially in the numerical computations . Additionally, in [25] the class of the relaxation function is a special case of the one in our paper. The decay results in both [25] and [26] were without numerical tests and without nonlinear source term.
• Liu and Yu [31] investigated the same problem where the exponents m and q are of constant nonlinearity and the relaxation function g satisfies the condition g (t) ≤ −ξ(t)g(t). In our paper, we extend the work of Liu and Yu [31] in which the exponents m(·) and q(·) are functions of x where x is in a part of the boundary. Moreover, we use a wider class of relaxation functions; that is g (t) ≤ −ξ(t)H (g(t)) so that the class of the relaxation function in [31] is a special case. In addition, we provide some numerical experiments to illustrate our decay theories. • Messaoudi et al. [38] investigated a similar problem. However, the nonlinear damping terms are in the domain. In our case the nonlinear damping and source terms are localized in the boundary, which makes the computations are more difficult. Also, [38] did not provide numerical computation.
The remainder of this work is arranged in the following manner: In Section 2, we write some of the assumptions and materials that are needed for our work. In Section 3, we establish and prove the global existence result. In Section 4, we present our main decay result as well as some examples. Section 5 presents and proves some technical lemmas. In Section 6, we prove the main decay results. Finally, in Section 7, we show numerical simulations to support our theoretical findings.

Preliminaries
In this section, we present some background information on the Lebesgue and Sobolev spaces with variable exponents (see [40,41]) as well as some assumptions for the main result proofs . We will use the letter c to denote a generic positive constant.

Lemma 2.7. [42] [Embedding Property]
Let Ω be a bounded domain in R n with a smooth boundary ∂Ω. Assume that p, k ∈ C(Ω) such that then we have continuous and compact embedding W 1,p(.

Assumptions
The following assumptions are essential in the proofs of the main results in this work.
(A1) The relaxation function g : R + → R + is a C 1 nonincreasing function satisfying and there exists a C 1 function Ψ : (0, ∞) → (0, ∞) which is linear or it is strictly increasing and strictly convex C 2 function on (0, r] for some 0 and there exists a a C 1 nonincreasing function ϑ such that is a continuous function such that . (A4) The variable exponents m and q are given continuous functions on Γ 1 satisfying the log-Hölder continuity condition: where c > 0 and 0 < δ < 1.
We introduce the "modified energy" associated to our problem A direct differentiation, using (2.6), leads to Lemma 2.9. [43] Under the assumptions in (A1), we have, for any t ≥ t 0 ,

Existence
The local existence theorem is stated in this section, and its proof can be demonstrated by combining the arguments of [44][45][46]. We also state and show a global existence result on the initial data under smallness conditions on (u 0 , u 1 ).
(Ω) and assume that (A1) − (A4) hold . Then, there exists T > 0, such that problem (1.6) has a weak solution We will now go over the following functionals: and Clearly, we have Proof. Sine I is continuous and I(u 0 ) > 0, then there exists T m < T such that which gives Now, Using Youngs and Poincaré inequalities and the trace theorem, we get ∀t ∈ [0, T m ], where Notice that (3.7) shows that u ∈ L q(·) (Γ 1 × (0, T )).
. Then the solution of (1.6) is global and bounded.
Proof. It suffices to show that ∇u 2 2 + u t 2 2 is bounded independently of t. To achieve this, we use (2.7), (3.2) and (3.5) to get (3.8) Since I(t) and (g • ∇u)(t) are positive, Therefore where C is a positive constant, which depends only on q 1 and and the proof is completed.
• Consider g(t) = αe −σ(1+t) ν , t ≥ 0, 0 < ν < 1, α, σ > 0, and α is chosen so that (A1) holds, then g (t) = −σΨ(g(t)) with ϑ(t) = ν(1 + t) ν−1 and Ψ(s) = s. We next infer that the solution of (1.6) satisfies the following energy estimate under the conditions of Theorem 4.2 • Let and α is chosen such that hypothesis (A1) remains valid. Then where σ is a fixed constant. Then, we conclude for t large enough and some constant C > 0 that the solution of (1.6) satisfies the following energy estimate under the conditions of Theorem 4.2 The proofs of Theorem 4.1 and Theorem 4.2 will be done through several Lemmas.
To prove (5.3), we apply Young's and Poincaré's inequalities and the trace theorem to obtain where Next, for any δ we have, by the case m(x) ≥ 2, As a result of combining the estimates above, we arrive at By choosing δ = 8c 0 , then c δ (x) is bounded and hence (5.3) is obtained.
and for 1 < m 1 < 2, we have where the constant c m > 0 depends on m 1 , m 2 and , and h ε is defined earlier in Lemma (5.1).
Proof. Direct differentiation of F 2 and using (1.6) leads to Using Young's inequality and Lemma 5.1, we get (5.16) From Lemma (5.1) and Young's inequality, we get Now, for almost every x ∈ Ω, we get Next, for almost every x ∈ Ω, we obtain Using Young's, Hölder's, Poincaré's inequalities and Lemma 5.1, we have Further, we have .
To estimate the last term in (5.15), we use Young's inequality and Lemma 5.1, to obtain (5.23) The first term in (5.23) can be estimated as follows: Collecting all the above estimates with (5.15), we see that (5.13) is archived.
To prove (5.14), we start by re-estimating the fifth term in (5.15) as follows: Hence, (5.14) is established.
and for any t ≥ t 0 , Proof. The equivalence L ∼ E can be proved straightforward. For the proof of (5.28), we start combining (2.6), (2.7), (5.2) and (5.13) and recalling g (t) := εg(t) − h ε (t), to get: (5.29) Now, set g 0 = t 0 0 g(s)ds and select δ small enough so that Once δ is fixed, then c δ (x) is bounded and the choice of ε 1 = 3 8 g 0 ε 2 yields Then (5.29) becomes < g(s) and using the Lebesgue Dominated Convergence Theorem, we conclude So, there exists 0 < ε 0 < 1 such that if ε < ε 0 , then Now, choosing N large enough so that L ∼ E and

Proofs of the decay theorems
In this section, we prove Theorem 4.1 and Theorem 4.2.

Numerical tests
We give numerical simulations in this section to support our theoretical results in Theorems 4.1 and 4.2. We use the conservative Lax-Wendroff strategy presented in [49] to demonstrate the decay of two tests. To discretize the system (1.1), we use a second-order finite difference method (FDM) in time and space for the space-time domain Ω × (0, T ) = [0, 1] × (0, 25). The mixed boundary conditions in the system (1.1) could be viewed as a Dirichlet boundary condition on one hand and a Neuman boundary condition on the other. Let for instance, u 0 (x) = x and u 1 (x) = 1 − x. Then, the condition ∂u ∂n − t 0 g(t − s) ∂u ∂n ds + |u t | m(x)−2 u t = |u| q(x)−2 u, in Γ 1 × (0, T ) will apply to the following two tests: • TEST 1: In the first test, we set m(x) = q(x) = 2. We use the boundary condition at x = 1 (the term u 1 will be vanish at x = 1, while the right-hand side condition will have a nonzero starting value).
• TEST 2: In the second numerical test, we examine the case m(x) 2 and q(x) 2 for all x ∈ [0, 1]. We use the boundary at x = 0 (the term u 1 term will not vanish at x = 0 and the right hand side condition will be canceled). For this, we use the functions m(x) = q(x) = 2 + 1 1 + x for all x ∈ [0, 1].
To check that the implemented method and the run code are numerically stable, we use ∆t < 0.5∆x, satisfying the stability condition according to the Courant-Friedrichs-Lewy (CFL) inequality, where ∆t = 0.0025 represents the time step and ∆x = 0.01 the spatial step. The spatial interval [0, 1] is subdivided into 100 subintervals, whereas the temporal interval [0, T ] = [0, 25] is deduced from the stability condition above. We run our code for 10, 000 time steps using the following initial conditions: u(x, 0) = x(1 − x) and u t (x, 0) = 0, in [0, 1].
In Tests 1 and 2, we demonstrated the decay under the initial and boundary conditions. The plots in Figure 1 show the temporal wave evolution in cross sections. The three cross sections are taken at x = 0.75, 1.5, 2.25 (see Figure 1. The corresponding energies given by the "modified" equation (3.1) are presented in Figure 2. The damping behavior is well seen in both tests. The result shown in Figure  2 is equally important. As a result, the similarity decrease for the energy decay rates obtained in Test 1 and Test 2 can be clearly observed. We normalized the output by dividing the maximums value in order to compare the asymptotic convergence of the energy.
Finally, it should be stressed that our intention focuses is to show the energy decay represented in Figure 2. However, we remarked that there are some similarities in the energy decay behavior. Both functions have at least a polynomial decay. This is due to the initial conditions used for the problem. We believe that, for other choices of the initial solution, we could obtain a clear difference between the outputs of the energy function.

Conclusions
In this work, we considered a viscoelastic wave equation with boundary damping and variable exponents. We first proved the existence of global solutions and then we established optimal and general decay estimates depending on the behavior of the relaxation function and the nature of the variable exponent nonlinearity. We finally end our paper with some numerical illustrations. Working with variable exponents in the boundary is totally different from the earlier results and of much challenging. We compared our results with other related results and showed that our results improved and extended some earlier results in the literature.