Well-posedness and decay in a system of hyperbolic and biharmonic-wave equations with variable exponents and weak dampings

In this paper, we consider a coupled system of hyperbolic and biharmonic-wave equations with variable exponents in the damping and coupling terms. In each equation, the damping term is modulated by a time-dependent coefficient a(t) (or b(t)). First, we state and prove a well-posedness theorem of global weak solutions, by exploiting Galerkin’s method and some compactness arguments. Then, using the multiplier method, we establish the decay rates of the solution energy, under suitable assumptions on the time-dependent coefficients and the range of the variable exponents. We end our work with some illustrative examples.


Model
In this work, we consider the following coupled system of wave and plate equations, for the unknowns u and v, where T > 0, is a bounded domain of R n (n ∈ N * ) with a smooth boundary ∂ , ∂v ∂η denotes the external normal derivative of v on the boundary of ; a, b : [0, ∞) −→ (0, ∞) are two non-increasing C 1 -functions and m, r and p are given continuous functions on satisfying some conditions.See (H.1)-(H3).

Motivation
In recent years and with the advancement of sciences, equations and systems of nonlinear wave equations with variable exponents occur in the mathematical modeling of various physical phenomena such as flows of electro-rheological fluids or fluids with temperature dependent viscosity, filtration processes through a porous media and image processing, nonlinear viscoelasticity, or robotics, etc.For more details on the subject, the reader can see [1,8].Our coupled system of variable-exponent nonlinearities (P) can be regarded as a model of the interaction between two fields describing the motion of two nonlinear "smart" materials, such as the motion of a suspension bridge and the cables.This class of problems requires more sophisticated mathematical tools to be investigated and well understood.The Lebesgue and sobolev spaces with variable exponents proved to be the appropriate spaces for studying such problems.

Literature review
For a class of one wave equation, Antontsev [3] studied the equation u tt − div a|∇u| p(x,t)−2 ∇u − α u t − bu|u| σ (x,t)−2 = f, in × (0, T ) , where α > 0 is a constant and a, b, p, σ are given functions.Under specific conditions on the exponents, he proved the local and global existence of weak solutions and a blow-up result.In Guo and Gao [11] took σ (x, t) = r > 2 and established a finite-time blow-up result for certain solutions with positive initial energy.After that, Guo [12] applied an interpolation inequality and some energy inequalities to obtain an estimate of the lower bound for the blow-up time when the source is super-linear.In Sun et al. [22] looked into the following equation established a blow-up result and gave lower and upper bounds for the blow-up time, under some conditions on the initial data.In addition, they provided numerical illustrations for their results.In Messaoudi and Talahmeh [17] considered the equation where a, b > 0 are two constants and m, r, p are given functions.They proved a finite-time blow-up result.In the absence of source term (b = 0), Messaoudi et al. [18] obtained decay estimates of solutions and presented two numerical applications as illustration for their theoretical results.After that, they gave in [19] an overview of results concerning decay and blow up for nonlinear wave equations involving variable and constant exponents.Very recently, Xiaolei et al. [23] used some energy estimates and some komornik's inequality to establish an asymptotic stability of solutions to quasilinear hyperbolic equations with variable source and damping terms.
Concerning coupled systems of hyperbolic equations with variable exponents, we mention the work of Bouhoufani and Hamchi [5], where they proved the global existence of weak solution and established decay estimates of the energy depending on the variable exponents.See also thesis [6].Messaoudi and Talahmeh [20] considered a system of wave equations, with damping and source terms of variable-exponent nonlinearities, and proved a blow-up result for solutions with negative initial energy.Very recently, Messaoudi et al. [21], studied a coupled hyperbolic system with variable exponents.They, obtained an existence and uniqueness result of a weak solution, showed that certain solutions, with positive initial energy, blow up in finite time and gave some numerical applications.
For the case of systems with biharmonic operators and variable-exponents of nonlinearities, we cite the work of Bouhoufani et al. [4], in which the authors considered, in a bounded domain, two biharmonic-wave equations with nonlinear dampings and source terms.They established a local existence, uniqueness and blowup result for solutions with negative-intial energy, and illustrated their theoretical findings by presenting some numericals tests.In Bouhoufani et al. [7] proved a theorem of finite-time blow up, for certain solution with positive initial data, and obtained the global existence as well as the decay rates, for the same problem, under suitable assumptions on the exponents and the initial data.After that, Messaoudi et al. [16] studied a coupled system of Laplacian and bi-Laplacian equations with dampings and source terms.They established existence, uniqueness, a long-time asymptotic behavior and blow up for solutions with positive-initial energy.In their case, with the presence of source terms, the solutions are local and they could blow up or exist globally (in time), depending on the range of the variable exponents and the initial data.However, in our study, the local solution doesn't cease to exist in finite time, due to the nature of the coupling terms we have in our problem.

Main contribution
In this work, we intend to prove the local and global well-posedness for the problem (P) and establish explicit decay rates of the solution energy depending on the range of the variable exponents m, r and the time-dependent coefficients a and b.This paper consists of six sections.After the introduction, we recall the definitions of the variable-exponent Lebesgue and Sobolev spaces as well as some important lemmas related to these spaces.In section three, we state and establish the existence result of a weak solution of problem (P).Section four is devoted to the statement and the proof of our aim result of stabilty.In section five, we present some illustrative examples and end with a conclusion.

Definitions and essential tools
In this section, we presente some important facts from [2,10,13] related to the Lebesgue and Sobolev spaces with variable exponents.Let q : −→ [1, ∞) be a measurable function, where is a domain of R n .We define the Lebesgue space with a variable exponent by Endowed with the following Luxembourg-type norm L q(.) ( ) is a Banach space (see [2,13]).We, also, define the variable exponent Sobolev space W 1,q(.) ( ) as follows: This is a Banach space with respect to the norm Definition 2.1 We say that a function q : −→ R is log-Hölder continuous on , if there exists a constant θ > 0, such that for all 0 < δ < 1, we have In addition, for q satisfying the log-Hölder continuity, we denote by W 1,q(.) 0 ( ) the closure of C ∞ 0 ( ) in W 1,q(.) ( ) and by W −1,q (.) ( ) the dual space of W 1,q(.) 0 ( ), in the same way as the usual Sobolev spaces, where 1 q(.) + 1 q (.) = 1.Lemma 2.2 [2,13] (Young's Inequality) Let r, q, s ≥ 1 be measurable functions defined on , such that , for a.e y ∈ .
Lemma 2.4 (Embedding Property ) [9] Let q : −→ [1, ∞) be a measurable function and k ≥ 1 be an integer.Suppose that r is a log-Hölder continuous function on , such that, for all x ∈ , we have Then, the embedding W k,r (.) 0 ( ) → L q(.) ( ) is continuous and compact.
To establish our stability result, the following Lemma is necessary.

Notations and assumptions
Throughout this paper, we denote by V the following space So, V is a separable Hilbert space endowed with the inner product and norm, respectively, The assumptions on m, r and p, that will be used in the sequel, are as follows.For all x ∈ , we suppose that where Remark 2.6 Since m, r and p are C 1 ( ), then they satisfy the log-Hölder continuity condition.

Global well-posedness result
In this section, our goal is to prove a local and global existence theorem of weak solutions of (P).For this purpose, we introduce the definition of a weak solution of problem (P).
and (u, v) satisfies, We have the following well-posedness result.

Theorem 3.2 Assume that (H.1)-(H.3) are satisfied. Then, for any initial data
there exists a weak solution (u, v) of (P) (in the sense of Definition 3.1) defined in [0, T ), for all T > 0.
Proof We procced in several steps: Step 1.Consider T > 0 fixed but arbitrary.Let ω j ∞ j=1 be an orthonormal basis of V and V k = span {ω 1 , ω 2 , . . ., ω k } be the subspace generated by the first k vectors ω 1 , ω 2 , . . ., ω k .Consider for j = 1, 2, . . ., k, with the initial data such that (3.4) For any k ≥ 1, Eq. (3.2) generate a system of k nonlinear ordinary differential equations, which admits a unique local solution (u k , v k ) defined on [0, t k ), 0 < t k ≤ T , by the standard theory of ODE.In the following step, our purpose is to extend this solution to [0, T ), for any k ≥ 1.
Step 2. Multiplying both sides of (3.2) 1 and (3.2) 2 by a j (t) and b j (t), respectively, using Green's formula and the boundary conditions, and then summing each result over j, from 1 to k, we obtain, for all 0 Adding (3.5) and (3.6), we get We then integrate (3.7) over (0, t), with t ≤ t k , to arrive at p(x) dx thanks to the convergence (3.4).Note here that the last term in the right-hand side of (3.8) is finite by (H.3).
It can be established exactly by the same calculations of (3.16) below.Also, under the hypothesise (H.3), we conclude that where C T is a constant independante of t and k.Therefore, we can extend the (u k ) k and (v k ) k on [0, T ).Moreover, we have b ( × (0, T )) . (3.9) Step 3. From (3.9), there exist two subsequences of (u k ) k and (v k ) k , still denoted by (u k ) k and (v k ) k (for simplicity), respectively, and two functions u, (3.10) Next, we show that and |v| p(.)−2 v|u| p(.) in L 2 ( × (0, T )) . (3.12) By the convergences (3.10) 1 and (3.10) 2 , the fact that H i 0 ( ) → compact L 2 ( ) (i = 1, 2) and invoking Lions' Theorem [14], there exist two subsequences of (u k ) k and (v k ) k , denoted by (u k ) k and (v k ) k , respectively, such that for all T > 0. The continuity of the function: and the convergences (3.13) lead to and ) a.e. in × (0, T ).
On the other hand, applying Young's inequality, with we obtain, for a.e.x ∈ , where .
Since p is bounded on , C(x) is bounded too.Hence, it comes, for some C 1 > 0 and for a.e.x ∈ , where C > 0 is a generic positive constant.It follows, for some CT > 0, .) is bounded in L 2 ( × (0, T )).This result with (3.14) allow us to establish (3.11), by virtue of Lions' Lemma.In similar way, we obtain (3.12).For the damping terms, we claim that Indeed, using Hölder's inquality and the fact that and Therefore, there exist two subsequences of (|u By repeating the same steps of [21] for the sequences (S k ) k , ( Sk ) k defined, for all k ≥ 1, as and and establish that (u, v) satisfies the two differential equations of (P), on × (0, T ) (in the weak sense), for all T > 0.
Step 4. As in [21], we easily establish that (u, v) satisfies the initial conditions.Finally, we conclude that (u, v) is a global weak solution of (P).
Remark 3.3 Note that the uniqueness of the solution remains open.However, if a(.) = b(.),we can obtain uniqueness by repeating the same steps of [21].

The decay rates
In order to state and prove our stability result, we define the energy functional associated to problem (P), by Multiplying the first equation of (P) by u t , the second one by v t , integrating each result over , using Green's formula and the boundary conditions, and then summing up, we obtain for a.e t ∈ [0, T ).

123
Proof Let S ∈ (0, T ).Multiplying (P) 1 and (P) 2 by E q u and E q v, respectively, for q ≥ 0 (to be specified later), integrating each result over × (S, T ) and using Green's formula, we obtain and Adding and subtracting the following two terms to (4.2) and (4.3), respectively, and then adding the two resulting equations, we infer Recalling the expression of the energy, (4.4) leads to where ( E q + q E q−1 ) (uu t + vv t )dxdt, In what follows, we estimate I i , for i = 1, . . ., 4. First, using Young's and Poincaré's inequalities and the definition of E, we obtain where c e is the Poincaré constant.Therefore, recalling (4.1), we infer where C is a generic positive constant.Next, using E (t) ≤ 0, we get For the third term, we set with To estimate J 1 , we consider the following partition of Therefore, by Hölder's inequality and the definition of α + , we obtain This yields ) and the definition of .Similarly, we find Adding J 1 and J 2 , it results Two cases are possible: Case 1: if α + = 2 then, Case 2: if α + > 2, we exploit Young's inequality, with δ = q + 1 and δ = (q + 1) /q, to get, for all ε > 0, If we take ε = 1 2C and q = α + 2 − 1, then Therefore, for α + ≥ 2, Finally, we handle I 4 as follows.Since a and b are bounded functions on R + , then where Now, as in [18], applying Young's inequality with we obtain, for all ε > 0, where Similarly, where By addition, we find Under the hypothesis (H.1) and (H.2), and recalling that m − , r − ≥ 2, we have the following estimate 123 Moreover, C ε (.) and C ε (.), in (4.12), will be bounded since m(.) and r (.) are bounded.Consequently, inequality (4.12) leads to 1+t , then the solution energy of (P) decreases as for c 3 = ce ω > 0.

Conclusion
In this work, we considered a coupled system of two weakly damped wave and plate equations, with Laplacian and bi-Laplacian operators and with variable exponent in the damping and coupling terms.We proved a theorem of global well-posedness and established different decay rates of the energy, depending on the variable exponents m(x) and r (x), and on the function coefficients a(t) and b(t).To illustrate our theoretical findings, we also gave some examples.This work generalizes many other works in the literature, in particular, those concerning the special case, where a and b are constants.

2 2 − 5 Examples
α + , i f α + > 2, ce −ωt , i f α + = 2.Based on Theorem 4.1, we present the following examples to illustrate different types of the energy decay, depending on the value of α + and on the damping coefficients a(t) and b(t).Example 5.1 For a(t) = 1 c 1 +t , b(t) = 1 c 2 +t , with 0 < c 1 ≤ c 2 , the estimate in Theorem 4.1 leads to
.5) On the other hand, we have for a.e.t ∈ [S, T ] + vv t ) dx = E q(uu t + vv t ) dx + E q (uu t + vv t ) t dx; I 1 + I 2 + I 3 + I 4 , ( q (uu t + vv t ) t dx = d dt E q (uu t + vv t ) dx − E q (uu t + vv t ) dx.
2 |u t | m(x) + |v t | r (x) dxdt Remark 4.2 As a special case, when a and b are constants, we have the following corollary.