Next Article in Journal
Moderate Deviation Principle for Linear Processes Generated by Dependent Sequences under Sub-Linear Expectation
Next Article in Special Issue
Quadratic Phase Multiresolution Analysis and the Construction of Orthonormal Wavelets in L2(ℝ)
Previous Article in Journal
Sturm-Liouville Problem with Mixed Boundary Conditions for a Differential Equation with a Fractional Derivative and Its Application in Viscoelasticity Models
Previous Article in Special Issue
An Equivalent Form Related to a Hilbert-Type Integral Inequality
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence and General Energy Decay of Solutions to a Coupled System of Quasi-Linear Viscoelastic Variable Coefficient Wave Equations with Nonlinear Source Terms

1
School of Mathematics, Chengdu Normal University, Chengdu 611130, China
2
School of Mathematics, Suqian University, Suqian 223800, China
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(8), 780; https://doi.org/10.3390/axioms12080780
Submission received: 2 July 2023 / Revised: 9 August 2023 / Accepted: 10 August 2023 / Published: 11 August 2023

Abstract

:
Viscoelastic damping phenomena are ubiquitous in diverse kinds of wave motions of nonlinear media. This arouses extensive interest in studying the existence, the finite time blow-up phenomenon and various large time behaviors of solutions to viscoelastic wave equations. In this paper, we are concerned with a class of variable coefficient coupled quasi-linear wave equations damped by viscoelasticity with a long-term memory fading at very general rates and possibly damped by friction but provoked by nonlinear interactions. We prove a local existence result for solutions to our concerned coupled model equations by applying the celebrated Faedo-Galerkin scheme. Based on the newly obtained local existence result, we prove that solutions would exist globally in time whenever their initial data satisfy certain conditions. In the end, we provide a criterion to guarantee that some of the global-in-time-existing solutions achieve energy decay at general rates uniquely determined by the fading rates of the memory. Compared with the existing results in the literature, our concerned model coupled wave equations are more general, and therefore our theoretical results have wider applicability. Modified energy functionals (can also be viewed as certain Lyapunov functionals) play key roles in proving our claimed general energy decay result in this paper.

1. Introduction

We are concerned, in this paper, with the initial boundary value problem (IBVP) for a coupled system of two quasi-linear space-variable coefficient wave equations whose energy is inhibited by viscoelastic dampings with long-term memories and possibly inhibited by frictional dampings, but provoked by nonlinear interactions. More precisely, we consider
| t u | ρ 1 t 2 u μ 1 div ( A 1 u ) div ( A 1 t 2 u ) + t g 1 ( t s ) div ( A 1 u ) ( s ) d s + a 11 t u + a 12 t v = f 1 ( u , v ) in Ω × ( 0 , + ) , | t v | ρ 2 t 2 v μ 2 div ( A 2 v ) div ( A 2 t 2 v ) + t g 2 ( t s ) div ( A 2 v ) ( s ) d s + a 21 t u + a 22 t v = f 2 ( u , v ) in Ω × ( 0 , + ) , u = v = 0 on Γ × ( 0 , + ) , u = u 0 , v = v 0 in Ω × ( , 0 ) , t u ( 0 ) = u 1 , t v ( 0 ) = v 1 in Ω ,
in which: u = u ( x , t ) and v = v ( x , t ) , ( x , t ) Ω × ( 0 , + ) , are the unknowns of IBVP (1); Ω is a nonempty bounded open subset of the N-dimensional Euclidean space R N , of which, the boundary, denoted by Γ (i.e., Γ = Ω ), is smooth enough (say, is in the C 2 class); N is a given positive integer; ρ 1 and ρ 2 , as well as μ 1 and μ 2 , are given positive constants; A 1 and A 2 are given R N × N -valued functions depending merely on space variables; a 11 , a 12 , a 21 and a 22 are given functions which depend merely on space variables; g 1 and g 2 , the so-called relaxation functions, are given functions mapping R + (throughout this paper, R + denotes the closed interval [ 0 , + ) ; see our notational conventions at the rear of this section) into itself; f 1 and f 2 are given real-valued functions defined in R 2 ; the given functions u 0 , v 0 , u 1 and v 1 are initial data of the unknowns of IBVP (1); t denotes the partial differential operator t ; ∇ denotes, as usual, the gradient operator ( x 1 , , x N ) on the N-dimensional Euclidean space R N , with x k denoting the partial differential operator x k , k = 1 , , N ; div φ denotes formally the divergence of the vector field
φ = ( φ 1 , , φ N ) = ( φ 1 ( x 1 , , x N ) , , φ N ( x 1 , , x N ) ) ,
that is, div φ = k = 1 N x k φ k . We shall explain later the sense in which the coupled quasi-linear viscoelastic variable coefficient wave Equations (1) 1 and (1) 2 , as well as the homogeneous Dirichlet boundary condition ‘ u = v = 0 on Γ × ( 0 , + ) ’ (i.e., (1) 3 ) are satisfied.
As with coupled parabolic equations (see [1,2] and the references cited therein), the coupled wave equations in IBVP (1) have important implications in Physics. The assumption that the coefficients A 1 and A 2 depend on space variables indicates that the underlying media/material is inhomogeneous. The assumption that the constants ρ 1 and ρ 2 are positive indicates that some of the structural properties of the concerned media/material are influenced significantly by the vibrating velocity. To include the terms
t g 1 ( t s ) div ( A 1 u ) ( s ) d s , and t g 2 ( t s ) div ( A 2 v ) ( s ) d s
in the model Equation (1) 1 and (1) 2 , we stress that, in our concerned scenario in this paper, the wave motions of the concerned media/material are suppressed by its viscoelasticity (the kinetic energy is inhibited, the viscosity is influenced by the velocity, and the aftereffect or memory of div ( A 1 u ) sustains for infinitely long time in the media/material); see Reference [3] for the description of viscoelasticity phenomenon and the explanation of the inducing mechanism of this phenomenon. The terms a 11 t u + a 12 t v , a 21 t u + a 22 t v , f 1 ( u , v ) and f 2 ( u , v ) are incorporated to emphasize that the waves u and v are ‘strongly’ coupled to a certain extent; we shall impose some suitable conditions on the coefficients a 11 , a 12 , a 21 and a 22 (see Assumption 5) to ensure that the term a 11 t u + a 12 t v , together with the term a 21 t u + a 22 t v , plays a role as frictional damping.
Partial differential equations describing the dynamics of viscoelastic materials have enormous implications to applications of these materials in engineering and scientific communities; the governing equations incorporate hereditary terms to stress that the aftereffect in the materials can not be neglected (see References [3,4,5,6]). In theoretical study or engineering applications, the aftereffect of some materials could be neglected for sufficiently large time, while the aftereffect of the other materials could last in infinitely long time periods. Let us point out again that, as indicated by the structure of the model equations in IBVP (1), the aftereffect of the material concerned in this paper could last in infinitely long time periods.
Since viscoelastic materials play important roles in diverse application areas (as alluded before), many experts in mathematical communities have been, in the last two decades, attracted into studying the dynamics of viscoelastic materials from mathematical perspectives. Muñoz Rivera [4] studied the large time behaviour of a class of viscoelastic equation, defined in bounded open subset of Euclidean spaces, in which the aftereffect in large time was neglected, and proved that the associated energy decays exponentially as time approaches infinity. Muñoz Rivera, Lapa and Barreto [5] established later some similar energy decaying results for plate equations. Aassila, Cavalcanti and Soriano [7] established exponentially decaying and polynomially decaying estimates for the energy of a constant-coefficient wave equation governing the vibration of materials occupying a domain whose boundary is of viscoelasticity under different conditions, respectively; in the meanwhile, they justified that the assertion that the energy approaches zero as time goes to infinity holds for all linear viscoelastic wave equations on bounded domains subject to homogeneous Dirichlet boundary condition. The idea in References [5,7] is strikingly illuminating for later study of problems for viscoelastic wave equations; see References [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27] and the vast references cited therein. For example, Cavalcanti, Domingos Cavalcanti and Ferreira [8] considered the following initial boundary value problem
| t u | ρ t 2 u Δ u Δ t 2 u + 0 t g ( t s ) Δ u ( s ) d s γ t u = 0 in Ω × ( 0 , + ) , u = v = 0 on Γ × ( 0 , + ) , u = u 0 , t u ( 0 ) = u 1 in Ω ;
they proved, under certain conditions on the relaxation function g, that IBVP (2) admits global weak solutions in H 0 1 ( Ω ; R 2 ) whenever γ 0 , and that the energy E ( t ) associated to the corresponding u decays exponentially whenever γ > 0 , where
E ( t ) = 1 ρ + 2 Ω | t u ( t ) | ρ + 2 d x + 1 2 Ω | u ( t ) | 2 d x + 1 2 Ω | t u ( t ) | 2 d x , t R + .
The (quasi-)linear wave equation (when ρ = 0 , the model Equation (2) 1 is linear) in IBVP (2) includes two damping terms, namely, the viscoelastic damping 0 t g ( t s ) Δ u ( s ) d s and the frictional damping t u . If ρ = 0 , these two dampings seem to be equivalent in the sense that the energy of both IBVP (2) incorporating merely viscoelastic damping and IBVP (2) incorporating merely frictional damping decays exponentially. Hence it is interesting to compare the intensity of these two terms in inhibiting the energy of solutions IBVP (2) with ρ > 0 . In this direction, Cavalcanti and Portillo Oquendo [9] obtained some interesting results. Berrimi and Messaoudi [10] studied viscoelastic equations including nonlinear source terms, and proved that the associated energy decays to zero as time goes to infinity whenever the initial values is sufficiently small. Cavalcanti, Domingos Cavalcanti and Martinez [11] extended in a certain sense the results in Reference [8] and proved that the energy of viscoelastic equations with general relaxation function g (which has a slow decaying rate compared to the one in Reference [8]) could also approach zero as time goes to infinity. The system of coupled viscoelastic wave equations has also been studied by several mathematicians in recent years. Han and Wang [12] studied the initial boundary value problem for a coupled system of viscoelastic wave equations with two nonlinear frictional damping terms, that is
t 2 u Δ u + 0 t g 1 ( t s ) Δ u ( s ) d s + | t u | m 1 t u = f 1 ( u , v ) in Ω × ( 0 , + ) , t 2 v Δ v + 0 t g 2 ( t s ) Δ v ( s ) d s + | t v | r 1 t v = f 2 ( u , v ) in Ω × ( 0 , + ) , u = v = 0 on Γ × ( 0 , + ) , u ( 0 ) = u 0 , t u ( 0 ) = u 1 , v ( 0 ) = v 0 , t v ( 0 ) = v 1 in Ω ;
they proved, under some additional conditions, that IBVP (3) is globally well-posed and provided a blow up criterion for IBVP (3) under some other conditions; as alluded above, the nonlinear terms | t u | m 1 t u and | t v | m 1 t v , playing roles as frictional dampings, bring in dissipation mechanism in the energy of the system (3). Said-Houari, Messaoudi and Guesmia [13] and Mustafa [14] extended the results in References [8,9,10,11] to IBVP (3) with the nonlinear frictional dampings removed. As could be seen evidently from IBVP (3), the structural properties of the concerned materials do not depend on the velocity of the vibration of the materials. Liu [15] considered an initial boundary value problem which is more close to IBVP (1), our working model problem in this paper, namely
| t u | ρ t 2 u Δ u γ 1 Δ t 2 u + 0 t g 1 ( t s ) Δ u ( s ) d s + f 1 ( u , v ) = 0 in Ω × ( 0 , + ) , | t v | ρ t 2 v Δ v γ 2 Δ t 2 v + 0 t g 2 ( t s ) Δ v ( s ) d s + f 2 ( u , v ) = 0 in Ω × ( 0 , + ) , u = v = 0 on Γ × ( 0 , + ) , u ( 0 ) = u 0 , t u ( 0 ) = u 1 , v ( 0 ) = v 0 , t v ( 0 ) = v 1 in Ω ,
with ρ > 0 ; he established under certain additional conditions some uniform decaying estimates for the energy of IBVP (4). He [16] reported some uniform decaying results for the energy associated to IBVP (4) under some other conditions. The other more interesting existence and stability results concerning viscoelastic (quasi-)linear wave equations could be seen in References [17,18,19,20,21,22,23,24,25,26,27,28,29] and the references therein.
As can be infered from the above review: The model equations considered in the aforementioned references only reflect that the concerned materials are homogeneous in all directions and that the aftereffect is all neglected for large time. But it seems to be more realistic that the materials are inhomogeneous in directions and that the aftereffect could influence the materials all the time. This motivates us to study space-varying viscoelastic wave equations with infinitely long memory of which both improve the mathematical difficulty of the paper. As could be seen later, we shall not assume that ρ 1 = ρ 2 ; it is obvious that ρ 1 ρ 2 has significant physical implications. Aside from these innovations, the nonlinearity in IBVP (1) seems to be more general than those studied in the existing references. Our goal in this paper is to prove under some conditions that solutions to IBVP (1) exist globally in time whenever their initial values are sufficiently small, and prove under some additional conditions that the energy associated to some of the global in time solutions approaches zero as time escapes to infinity.
Assumption 1.
For i = 1 , 2 , A i C 2 ( Ω ¯ ; R N × N ) , the set of uniformly continuous functions of which partial derivatives whose orders not exceeding 2 are all uniformly continuous Ω. For every x Ω , the matrix A i ( x ) is symmetric, i = 1 , 2 . And ϑ i is assumed to be strictly positive, where the constant ϑ i is given by
ϑ i = inf x Ω , ξ R N { 0 } ξ A i ( x ) ξ ξ ξ , i = 1 , 2 .
Assumption 2.
The relaxation function g i is strictly monotonically decreasing, maps the closed interval R + into itself and satisfies
0 < 0 + g i ( s ) d s < μ i , i = 1 , 2 .
The derivative function g i , of the relaxation function g i , is locally Lebesgue integrable in R + (in other words, the function g i is absolutely continuous in the interval R + ), i = 1 , 2 . There exists a nonincreasing absolutely continuous function ξ i mapping R + into itself and a function K i C 2 ( [ 0 , r i ] ; R ) (with r i a given positive constant not less than g i ( 0 ) ), which is strictly increasing and strictly convex, and satisfies K i ( 0 ) = K i ( 0 ) = 0 , such that
K 1 ( t ) K 2 ( t ) or K 2 ( t ) K 1 ( t ) for all t [ 0 , min ( r 1 , r 2 ) ] ,
that
g i ( t ) ξ i ( t ) K i ( g i ( t ) ) for all t R + ,
and that there exists a positive constant k i satisfying
lim δ 0 + K i 1 ( δ ) max ( K 1 , K 2 ) 1 ( δ ) = k i , i = 1 , 2 ,
with max ( K 1 , K 2 ) 1 denoting the inverse of the derivative of the function max ( K 1 , K 2 ) .
Assumption 3.
The constants ρ 1 and ρ 2 satisfy min ( ρ 1 , ρ 2 ) > 0 and ( N 2 ) max ( ρ 1 , ρ 2 ) 2 .
Assumption 4.
f i is locally Lipschitz continuous in R 2 , and satisfies f i ( 0 , 0 ) = 0 , i = 1 , 2 . There exists a function F ( u , v ) , defined in the whole space R 2 , such that
F ( u , v ) = 0 u f 1 ( u ˜ , 0 ) d u ˜ + 0 v f 2 ( u , v ˜ ) d v ˜ , u , v R ,
or equivalently v f 1 ( u , v ) = u f 2 ( u , v ) , u, v R . f i satisfies the growth condition at infinity: There exist four absolute constants Ł 1 , Ł 2 , p 1 and p 2 satisfying min ( Ł 1 , Ł 2 , p 1 , p 2 ) > 0 and ( N 2 ) max ( p 1 , p 2 ) 2 , such that, for every pair ( u , v ) R 2 , it holds always that
| f i ( u , v ) | Ł i ( | u | p 1 + 1 + | u | p 2 + 1 + | u | ρ 1 + 1 + | v | p 1 + 1 + | v | p 2 + 1 + | v | ρ 2 + 1 ) , i = 1 , 2 .
Assumption 5.
For every x Ω : a 12 ( x ) = a 21 ( x ) ; and the matrix
a 11 ( x ) a 12 ( x ) a 21 ( x ) a 22 ( x )
is semi-positive definite. a i j belongs to the Banach space C ( Ω ¯ ) , the totality of uniformly continuous real-valued functions defined in Ω, i, j = 1 , 2 .
For the sake of convenience of our later presentation, we write
1 ζ i = inf φ H 0 1 ( Ω ) { 0 } 1 φ L 2 ( Ω ; R N ) 2 Ω φ A i φ d x , i = 1 , 2 ,
and write for every p 1 satisfying ( N 2 ) p 2 N :
κ p = sup φ H 0 1 ( Ω ) { 0 } φ L p ( Ω ) φ L 2 ( Ω ; R N ) .
For every given α R + , we denote henceforth
α = sup ( x , y ) R + 2 { ( 0 , 0 ) } x α + y α ( x + y ) α = 2 max ( 1 α , 0 ) .
Hereafter, we associate to a i j (i, j = 1 , 2 ) the following constant
ϰ a 11 , a 12 , a 21 , a 22 = sup ( φ 1 , ψ 1 ) , ( φ 2 , ψ 2 ) L 2 ( Ω ; R 2 ) , ( φ 1 , ψ 1 ) L 2 ( Ω ; R 2 ) = ( φ 2 , ψ 2 ) L 2 ( Ω ; R 2 ) = 1 | Ω φ 1 ψ 1 a 11 a 12 a 21 a 22 φ 2 ψ 2 d x | .
To improve the readability, we would like to give some remarks on our list of five standing assumptions (see Assumptions 1–5 for the details) of this paper.
Remark 1.
By applying the celebrated Sobolev–Poincaré inequality (see [1], THEOREM 1, p. 292), we can conclude that κ p given by (10) is a positive constant. In light of Assumption 1 on the coefficient matrix A i ( x ) , the term Ω φ A i φ d x is well-defined and ζ i given by (9) is also a positive constant obeying
1 tr ( A i ) 2 L ( Ω ) = 1 tr ( A i ) A i L ( Ω ) ζ i 1 ϑ i ,
where ϑ i is given by (5), and tr denotes the trace operator of square matrice, i = 1 , 2 .
Remark 2.
By imposing the restriction (6) in Assumption 2, our principal aim in this paper is to guarantee the twice continuous differentiability of the function max ( K 1 , K 2 ) in the interval [ 0 , min ( r 1 , r 2 ) ] . The restriction (6) can be probably removed via introducing the notion of subdifferential in convex analysis, or via utilizing Dini’s derivatives, together with some complicated calculations techniques.
Remark 3.
Let α R + . We can prove, based on the notation α given by (11), that
x α + y α α ( x + y ) α , ( x , y ) R + 2 .
Based on Assumption 4 ((7) and (8), in particular), we could obtain the following lemma by some routine but tedious calculations.
Lemma 1.
Let F ( u , v ) , f 1 ( u , v ) and f 2 ( u , v ) be three functions given as in Assumption 4. For every pair ( u , v ) R 2 , it holds always that
| F ( u , v ) | Ł 1 + Ł 2 ( p 1 + 1 ) p 1 + 2 | u | p 1 + 2 + Ł 1 + Ł 2 ( p 2 + 1 ) p 2 + 2 | u | p 2 + 2 + Ł 1 + Ł 2 ( ρ 1 + 1 ) ρ 1 + 2 | u | ρ 1 + 2 + 2 Ł 2 | v | p 1 + 2 p 1 + 2 + 2 Ł 2 | v | p 2 + 2 p 2 + 2 + Ł 2 | v | ρ 1 + 2 ρ 1 + 2 + Ł 2 | v | ρ 2 + 2 ρ 2 + 2 ,
where Ł i is exactly the one given in (8) in Assumption 4, i = 1 , 2 .
Remark 4.
By some routine calculations, we have immediately that
( ϰ a 11 , a 12 , a 21 , a 22 ) 2 ( a 11 ) 2 + ( a 12 ) 2 + ( a 21 ) 2 + ( a 22 ) 2 L ( Ω ) a 11 L ( Ω ) 2 + a 12 L ( Ω ) 2 + a 21 L ( Ω ) 2 + a 22 L ( Ω ) 2 .
This, together with Assumption 5, implies that ϰ a 11 , a 12 , a 21 , a 22 , given by (12), is indeed a non-negative constant.
Notational Conventions.
R is the field of real numbers; R + = [ 0 , + ) ; R = ( , 0 ] . For φ : Ω R , we write formally φ = ( φ ) . For 1 p + , L p ( Ω ) denotes, as usual, the classical Lebesgue space. The Hilbert space H 1 ( Ω ) , equipped with the inner product
( φ , ψ ) Ω ( φ ψ + φ ψ ) d x ,
denotes the totality of square-integrable functions defined in Ω whose first order partial derivatives, in the distributional sense, are all square-integrable functions in Ω . H 0 1 ( Ω ) denotes the totality of functions in H 1 ( Ω ) having zero as their boundary values in the trace sense, or equivalently, H 0 1 ( Ω ) is the completion of the totality C comp ( Ω ) of infinitely differentiable functions defined in Ω having compact support in the Hilbert space H 1 ( Ω ) ; inheriting the inner product (13) from H 1 ( Ω ) , H 0 1 ( Ω ) is also a Hilbert space. We write · H 1 ( Ω ) and · H 0 1 ( Ω ) for the norms induced by inner products of the Hilbert spaces H 1 ( Ω ) and H 0 1 ( Ω ) . Let H be an inner product space, we write , H for the inner product of the space H. Let X be a Banach space with X its topological dual, and J an interval; we write , X , X for the duality pairing (of the Banach space X and its dual X ) and denote by C w ( J ; X ) the following space
φ : J X ; J t ψ , φ ( t ) X , X R for every ψ X .
The rest of this paper is organized as follows. In Section 2, we prove that solutions to IBVP (1) exist globally in time whenever their initial data are sufficiently small. In Section 3, we provide a general decaying estimate on global-in-time solutions to IBVP (1); the estimate implies that global-in-time solutions to IBVP (1) decrease to zero as time goes to infinity, whenever their initial data satisfy some additional conditions. In Section 4, we provide several concluding remarks.

2. Global Existence Results Concerning Solutions to IBVP (1)

In this section, our main aim is to prove the global existence of solutions, whose initial data belonging to a certain function space, to IBVP (1). We shall first demonstrate the local existence of solutions to IBVP (1) via utilizing the Faedo-Galerkin method, and shall then prove the desired global existence by establishing a priori estimates and a standard continuation argument. For the sake of convenience of our later presentation, we write, in the sequel, for every T ( 0 , + ) :
S [ 0 , T ] = { ( u , v ) ; u , v : Ω × ( , T ] R , u | Ω × [ 0 , T ] , v | Ω × [ 0 , T ] C w ( [ 0 , T ] ; H 0 1 ( Ω ) ) , t u | Ω × [ 0 , T ] , t v | Ω × [ 0 , T ] C w ( [ 0 , T ] ; H 0 1 ( Ω ) ) , u | Ω × R , v | Ω × R L ( R ; H 0 1 ( Ω ) ) } ,
and for every 0 < T + , we write similarly
S [ 0 , T ) = { ( u , v ) ; u , v : Ω × ( , T ) R , ( u , v ) | Ω × ( , T ˜ ] S [ 0 , T ˜ ] , T ˜ ( 0 , T ) } .
Definition 1.
Let T ( 0 , + ) . The pair ( u , v ) S [ 0 , T ] is said to be a local weak solution, in the interval [ 0 , T ] , to IBVP (1) provided that the following two equalities hold for every pair ( φ , ψ ) H 0 1 ( Ω ; R 2 ) of test functions:
1 ρ 1 + 1 Ω | t u ( t ) | ρ 1 t u ( t ) φ d x + Ω φ A 1 t u ( t ) d x 1 ρ 1 + 1 Ω | u 1 | ρ 1 u 1 φ d x Ω φ A 1 u 1 d x + 0 t Ω φ A 1 ( u ( s ) s g 1 ( s τ ) u ( τ ) d τ ) d x d s = 0 t Ω φ ( f 1 ( u ( s ) , v ( s ) ) a 11 t u ( s ) a 12 t v ( s ) ) d x d s , t ( 0 , T ) ,
and
1 ρ 2 + 1 Ω | t v ( t ) | ρ 2 t v ( t ) ψ d x + Ω ψ A 2 t v ( t ) d x 1 ρ 2 + 1 Ω | v 1 | ρ 2 v 1 ψ d x Ω ψ A 2 v 1 d x + 0 t Ω ψ A 2 ( v ( s ) s g 2 ( s τ ) v ( τ ) d τ ) d x d s = 0 t Ω ψ ( f 2 ( u ( s ) , v ( s ) ) a 21 t u ( s ) a 22 t v ( s ) ) d x d s , t ( 0 , T ) .
Definition 2.
Let 0 < T + . ( u , v ) : Ω × ( , T ) R 2 is said to be a weak solution to IBVP (1), in the interval [ 0 , T ) , if for every 0 < T < T , ( u , v ) | Ω × ( , T ] is a local weak solution, in the interval [ 0 , T ] , to IBVP (1). In the case that T = + (or equivalently, ( u , v ) S R + ; see (15) for the definition of S R + ), ( u , v ) is called a global weak solution to IBVP (1); otherwise, ( u , v ) is still called a local weak solution to IBVP (1).
To every solution pair ( u , v ) S [ 0 , T ) (see (15) for the definition of S [ 0 , T ) ), with 0 < T + , to IBVP (1), we associate the following functional (a certain Lyapunov functional candidate)
E u , v ( t ) = 1 ρ 1 + 2 t u ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + 1 ρ 2 + 2 t v ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 + 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + 1 2 Ω t u ( t ) A 1 t u ( t ) d x + 1 2 ( g 1 A 1 u ) ( t ) + 1 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + 1 2 Ω t v ( t ) A 2 t v ( t ) d x + 1 2 ( g 2 A 2 v ) ( t ) Ω F ( u ( t ) , v ( t ) ) d x , t [ 0 , T ) ,
where F is given by (7), and the operation “⋄”, associated to two given functions Φ L ( Ω ; R N × N ) and g L ( 0 , + ) , is defined, in a formal way, as follows: For every ψ L loc 2 ( [ , T ) ; L 2 ( Ω ; R N ) ) ,
( g Φ ψ ) ( t ) = t g ( t s ) Ω ( ψ ( t ) ψ ( s ) ) Φ ( ψ ( t ) ψ ( s ) ) d x d s .
Lemma 2.
Let Φ W 1 , ( Ω ; R N × N ) , g W 1 , ( 0 , + ) . For every
φ L loc 2 ( [ , T ) ; H 0 1 ( Ω ) )
satisfying φ | Ω × [ 0 , T ) H loc 1 ( [ 0 , T ) ; H 0 1 ( Ω ) ) , it holds that
d d t ( g Φ φ ) ( t ) 0 + g ( s ) d s Ω φ ( t ) Φ φ ( t ) d x = ( g Φ φ ) ( t ) 2 t g ( t s ) Ω t φ ( t ) Φ φ ( s ) d x d s , t ( 0 , T ) .
Proof. 
Conduct some routine calculations, to obtain
d d t ( g Φ φ ) ( t ) 0 + g ( s ) d s Ω φ ( t ) Φ φ ( t ) d x = ( g Φ φ ) ( t ) 2 0 + g ( s ) d s Ω t φ ( t ) Φ φ ( t ) d x + 2 t g ( t s ) Ω t φ ( t ) Φ ( φ ( t ) φ ( s ) ) d x d s = ( g Φ φ ) ( t ) 2 t g ( t s ) Ω t φ ( t ) Φ φ ( s ) d x d s , t ( 0 , T ) ,
in which the “=” in the second line follows immediately from the very definition of the expression ( g Φ φ ) ( t ) (see (19) for the detail). □
The differential identity in Lemma 2 is of great importance in our later calculations and will be used frequently in proving our main results in this paper. For example, we shall employ Lemma 2 as one of the main ingredients in the procedure of deducing an a priori inequality which plays a key role in proving the following local existence theorem.
Theorem 1.
Suppose that Assumptions 1–4 hold true, and that the space-varying coefficient a i j belongs to the Banach space C ( Ω ¯ ) , i, j = 1 , 2 . Then for every initial datum pair
( u 0 , v 0 ) L ( R ; H 0 1 ( Ω ; R 2 ) )
and every initial datum pair ( u 1 , v 1 ) H 0 1 ( Ω ; R 2 ) , IBVP (1) admits a local weak solution ( u , v ) S [ 0 , T ˜ ] (see (14) for the detailed definition of the notation S [ 0 , T ˜ ] ), in the interval [ 0 , T ˜ ] , in which T ˜ is a certain postive time instant depending merely on Ω, A 1 , A 2 , f 1 , f 2 , g 1 , g 2 , ρ 1 , ρ 2 , a 11 , a 12 , a 21 , a 22 , u 0 , v 0 , u 1 and v 1 .
Please notice that we do not use Assumption 5 in Theorem 1 temporarily, instead, we used a weaker condition that ‘ a i j belongs to the Banach space C ( Ω ¯ ) , i, j = 1 , 2 ’. We shall prove Theorem 1 via the very standard Faedo-Galerkin procedure.
Proof. 
Thanks to Assumption 1, by recalling theory on elliptic partial differential equations, one can find: One orthonormal basis, designated by { e i n } n = 1 , of the Hilbert space L 2 ( Ω ) is composed of the solutions of the following eigenvalue problems
div ( A i e i ) = λ i e i in Ω , e i = 0 on Γ .
By using mainly the divergence theorem, we have
Ω e i k A i e i d x = Ω e i k div ( A i e i ) d x = λ i Ω e i k A i e i d x = λ i if k = , 0 if k .
Since A i C 2 ( Ω ¯ ; R N × N ) (see Assumption 1), for any positive integer n, we have e i n C 2 ( Ω ¯ ) , i = 1 , 2. Let us introduce the following two sequences of approximate solutions
u k ( x , t ) = j = 1 k e 1 j u ˜ k j ( x , t ) = j = 1 k u ˜ k j ( t ) e 1 j ( x ) , v k ( x , t ) = j = 1 k e 2 j v ˜ k j ( x , t ) = j = 1 k v ˜ k j ( t ) e 2 j ( x ) , k N .
As with (20), to construct sequences of approximate solutions is one of the main steps in applying the Faedo-Galerkin scheme to prove local existence result of evolution partial differential equations. After some calculations, we can find that ( u k , v k ) given by (20) is approximate solution pair to IBVP (1) if and only if
U k ( t ) = ( u ˜ k 1 ( t ) , , u ˜ k k ( t ) , u ˜ k 1 ( t ) , , u ˜ k k ( t ) ,   v ˜ k 1 ( t ) , , v ˜ k k ( t ) , v ˜ k 1 ( t ) , , v ˜ k k ( t ) )
is the solution to the following Cauchy problem
Ω | j = 1 k u ˜ k j ( t ) e 1 j | ρ 1 j = 1 k u ˜ k j ( t ) e 1 j e 1 d x + λ 1 u ˜ k ( t ) + λ 1 u ˜ k ( t ) λ 1 t g 1 ( t s ) u ˜ k ( s ) d s + j = 1 k u ˜ k j ( t ) Ω a 11 e 1 j e 1 d x + j = 1 k v ˜ k j ( t ) Ω a 12 e 2 j e 1 d x = Ω f 1 ( j = 1 k u ˜ k j ( t ) e 1 j , j = 1 k v ˜ k j ( t ) e 2 j ) e 1 d x , t R + , = 1 , , k , Ω | j = 1 k v ˜ k j ( t ) e 2 j | ρ 2 j = 1 k v ˜ k j ( t ) e 2 j e 2 m d x + λ 2 m v ˜ k m ( t ) + λ 2 m v ˜ k m ( t ) λ 2 m t g 2 ( t s ) v ˜ k m ( s ) d s + j = 1 k v ˜ k j ( t ) Ω a 21 e 2 j e 2 m d x + j = 1 k v ˜ k j ( t ) Ω a 22 e 2 j e 2 m d x = Ω f 2 ( j = 1 k u ˜ k j ( t ) e 1 j , j = 1 k v ˜ k j ( t ) e 2 j ) e 2 m d x , t R + , m = 1 , , k , u ˜ k ( t ) = Ω u 0 ( t ) e 1 d x , v ˜ k ( t ) = Ω v 0 ( t ) e 2 d x , t R , = 1 , , k , u ˜ k ( 0 ) = Ω u 1 e 1 d x , v ˜ k ( 0 ) = Ω v 1 e 2 d x , = 1 , , k .
In accordance with (21), we write
U k 0 ( t ) = ( Ω u 0 ( t ) e 11 d x , , Ω u 0 ( t ) e 1 k d x , Ω u 1 e 11 d x , , Ω u 1 e 1 k d x , Ω v 0 ( t ) e 21 d x , , Ω v 0 ( t ) e 2 k d x , Ω v 1 e 21 d x , , Ω v 1 e 2 k d x ) , t R .
Then the Cauchy problem (22) can be recast into the following Cauchy problem for a functional differential equation
U k ( t ) = F k ( U k ( t ) , U k t ) , t R + , U k = U k 0 in R ,
where U k t ( τ ) = U k ( t + τ ) , τ R , and F k : R 4 k × L ( R ; R 4 k ) R 4 k is locally Lipschitz continuous. Then by a variant of the classical Cauchy–Lipschitz existence theory, the Cauchy problem (23), or equivalently, the Cauchy problem (22) is locally well-posed in Hadamard’s sense. As a consequence, by applying a standard continuation argument, we can prove, based on the aforementioned local well-posedness result, that the Cauchy problem (22) (for system of ordinary differential equations) admits a unique solution
( u ˜ k 1 ( t ) , , u ˜ k k ( t ) , v ˜ k 1 ( t ) , , v ˜ k k ( t ) )
in the classical sense in [ 0 , T * ) , the maximal interval of existence. Thus the pair ( u k , v k ) in the form (20) is well-defined. Let us now introduce the following auxilliary functional
Ξ k ( t ) = E u k , v k ( t ) + Ω F ( u k ( t ) , v k ( t ) ) d x = 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + 1 ρ 2 + 2 t v k ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 + 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u k ( t ) A 1 u k ( t ) d x + 1 2 Ω t u k ( t ) A 1 t u k ( t ) d x + 1 2 ( g 1 A 1 u k ) ( t ) + 1 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v k ( t ) A 2 v k ( t ) d x + 1 2 Ω t v k ( t ) A 2 t v k ( t ) d x + 1 2 ( g 2 A 2 v k ) ( t ) , t [ 0 , T * ) .
Differentiate Ξ k ( t ) and simplify further the obtained result, to yield
Ξ k ( t ) = 1 2 ( g 1 A 1 u k ) ( t ) + 1 2 ( g 2 A 2 v k ) ( t ) t g 1 ( t s ) Ω t u k ( t ) A 1 u k ( s ) d x d s t g 2 ( t s ) Ω t v k ( t ) A 2 v k ( s ) d x d s t g 1 ( t s ) Ω t u k ( t ) div ( A 1 u k ( s ) ) d x d s t g 2 ( t s ) Ω t v k ( t ) div ( A 2 v k ( s ) ) d x d s + Ω t u k ( t ) f 1 ( u k ( t ) , v k ( t ) ) d x + Ω t v k ( t ) f 2 ( u k ( t ) , v k ( t ) ) d x Ω t u k ( t ) t v k ( t ) a 11 a 12 a 21 a 22 t u k ( t ) t v k ( t ) d x Ω t u k ( t ) f 1 ( u k ( t ) , v k ( t ) ) d x + Ω t v k ( t ) f 2 ( u k ( t ) , v k ( t ) ) d x + ϰ a 11 , a 12 , a 21 , a 22 ( t u k ( t ) L 2 ( Ω ) 2 + t v k ( t ) L 2 ( Ω ) 2 ) , t [ 0 , T * ) .
Before continuing our proof, it is worth noticing that the unique existence of solutions to the Cauchy problem (22) and solutions to the Cauchy problem (22) depending continuously on their initial data are both attributed to the local Lipschitz continuity of F k in the right hand side of the partial differential equation in the Cauchy problem (22), and that the maximal existence time instant T * , independent of k, depends on Ω , A 1 , A 2 , f 1 , f 2 , g 1 , g 2 , ρ 1 , ρ 2 , a 11 , a 12 , a 21 , a 22 , u 0 , v 0 , u 1 and v 1 .
Thanks to
ρ 1 ρ 1 + 2 + 2 ρ 1 + 2 = 1 ,
by the Fenchel–Young inequality, it holds that
| a b | ρ 1 ρ 1 + 2 | a | ρ 1 + 2 ρ 1 + 2 ρ 1 + 2 | b | ρ 1 + 2 2 , a , b R ,
from which it follows further that
Ω a 11 | t u k ( t ) | 2 d x ρ 1 ρ 1 + 2 a 11 L ρ 1 + 2 ρ 1 ( Ω ) ρ 1 + 2 ρ 1 + 2 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 , t [ 0 , T * ) .
And similarly, we have
Ω a 22 | t v k ( t ) | 2 d x ρ 2 ρ 2 + 2 a 22 L ρ 2 + 2 ρ 2 ( Ω ) ρ 2 + 2 ρ 2 + 2 ρ 2 + 2 t v k ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 , t [ 0 , T * ) .
In view of the Fenchel–Young inequality
| a b c | 1 ρ 1 + 2 | a | ρ 1 + 2 + 1 ρ 2 + 2 | b | ρ 2 + 2 + ρ 1 ρ 2 + ρ 1 + ρ 2 ( ρ 1 + 2 ) ( ρ 2 + 2 ) | c | ( ρ 1 + 2 ) ( ρ 2 + 2 ) ρ 1 ρ 2 + ρ 1 + ρ 2 , a , b , c R ,
which follows directly from the identity
ρ 1 ρ 2 + ρ 1 + ρ 2 ( ρ 1 + 2 ) ( ρ 2 + 2 ) + 1 ρ 1 + 2 + 1 ρ 2 + 2 = 1 ,
with the aid of the experience gained in the procedure of deriving (26) and (27), we have, after some routine but careful calculations, that
| Ω a 12 t u k ( t ) t v k ( t ) d x | 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + 1 ρ 2 + 2 t v k ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 + ρ 1 ρ 2 + ρ 1 + ρ 2 ( ρ 1 + 2 ) ( ρ 2 + 2 ) a 12 L ( ρ 1 + 2 ) ( ρ 2 + 2 ) ρ 1 ρ 2 + ρ 1 + ρ 2 ( Ω ) ( ρ 1 + 2 ) ( ρ 2 + 2 ) ρ 1 ρ 2 + ρ 1 + ρ 2 , t [ 0 , T * ) .
And analogously, we have
| Ω a 21 t u k ( t ) t v k ( t ) d x | 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + 1 ρ 2 + 2 t v k ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 + ρ 1 ρ 2 + ρ 1 + ρ 2 ( ρ 1 + 2 ) ( ρ 2 + 2 ) a 21 L ( ρ 1 + 2 ) ( ρ 2 + 2 ) ρ 1 ρ 2 + ρ 1 + ρ 2 ( Ω ) ( ρ 1 + 2 ) ( ρ 2 + 2 ) ρ 1 ρ 2 + ρ 1 + ρ 2 , t [ 0 , T * ) .
By utilizing the growth condition (8) upon the nonlinearity f i ( i = 1 , 2 ) in Assumption 4, we obtain immediately
Ω t u k ( t ) f 1 ( u k ( t ) , v k ( t ) ) d x Ł 1 Ω | t u k ( t ) | | u k ( t ) | p 1 + 1 d x + Ł 1 Ω | t u k ( t ) | | u k ( t ) | p 2 + 1 d x + Ł 1 Ω | t u k ( t ) | | u k ( t ) | ρ 1 + 1 d x + Ł 1 Ω | t u k ( t ) | | v k ( t ) | p 1 + 1 d x + Ł 1 Ω | t u k ( t ) | | v k ( t ) | p 2 + 1 d x + Ł 1 Ω | t u k ( t ) | | v k ( t ) | ρ 2 + 1 d x , t [ 0 , T * ) .
By virtue of observing
1 ρ 1 + 2 + ρ 1 + 1 ρ 1 + 2 = 1 ,
inspired by the experience gained in deducing (26)–(29), we have, by applying also the Fenchel–Young inequality, immediately that
Ω | t u k ( t ) | | u k ( t ) | p 1 + 1 d x 1 ρ 1 + 2 Ω | t u k ( t ) | p 1 + 2 d x + ρ 1 + 1 ρ 1 + 2 Ω | u k ( t ) | ( p 1 + 2 ) ( ρ 1 + 2 ) ρ 1 + 1 d x = 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + ρ 1 + 1 ρ 1 + 2 u k ( t ) L ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( Ω ) ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 , t [ 0 , T * ) .
In view of the coercivity condition on A 1 (see Assumption 1), we have by applying the Sobolev–Poincaré inequality (see Remark 1) that
u k ( t ) L ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( Ω ) ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( κ ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ) ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 u k ( t ) L 2 ( Ω ; R N ) ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( κ ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ) ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( ζ 1 ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 ;
see (9) and (10) for the detailed explanations on the notations ζ 1 and κ ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 , respectively. Combine (31) and (32), to arrive at
Ω | t u k ( t ) | | u k ( t ) | p 1 + 1 d x 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + M 11 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , t [ 0 , T * ) ,
in which the positive constant M 11 1 is given by
M 11 1 = ρ 1 + 1 ρ 1 + 2 ( κ ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ) ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( ζ 1 ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 .
With the help of the experience of deriving (33), we conlude similarly that
Ω | t u k ( t ) | | v k ( t ) | p 1 + 1 d x 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + ρ 1 + 1 ρ 1 + 2 v k ( t ) L ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( Ω ) ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + M 12 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , t [ 0 , T * ) ,
in which the positive constant M 12 1 is given by
M 12 1 = ρ 1 + 1 ρ 1 + 2 ( κ ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ) ( p 1 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( ζ 2 ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 .
With the help of the experience of deriving (33) and (34), we arrive at
Ω | t u k ( t ) | | u k ( t ) | p 2 + 1 d x 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + ρ 1 + 1 ρ 1 + 2 u k ( t ) L ( p 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( Ω ) ( p 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + M 13 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , t [ 0 , T * ) ,
in which the positive constant M 13 1 is given by
M 13 1 = ρ 1 + 1 ρ 1 + 2 ( κ ( p 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ) ( p 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( ζ 1 ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 .
With the help of the experience of deriving (33)–(35), we conlude that
Ω | t u k ( t ) | | v k ( t ) | p 2 + 1 d x 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + ρ 1 + 1 ρ 1 + 2 v k ( t ) L ( p 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( Ω ) ( p 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + M 14 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , t [ 0 , T * ) ,
in which the positive constant M 14 1 is given by
M 14 1 = ρ 1 + 1 ρ 1 + 2 ( κ ( p 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ) ( p 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( ζ 2 ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 .
Based on the algebraic identity (30), we apply the Fenchel–Young inequality, to obtain
Ω | t u k ( t ) | | u k ( t ) | ρ 1 + 1 d x 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + ρ 1 + 1 ρ 1 + 2 u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 , t [ 0 , T * ) .
Mimicking steps in deriving (32), having the notations in (9) and (10) at our disposal, and based on the coercivity condition on A 1 (see Assumption 1), we apply the Sobolev–Poincaré inequality (see Remark 1), to arrive at
u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 ( κ ρ 1 + 2 ) ρ 1 + 2 u k ( t ) L 2 ( Ω ; R N ) ρ 1 + 2 ( κ ρ 1 + 2 ) ρ 1 + 2 ( ζ 1 ) ρ 1 + 2 2 ( Ω u k ( t ) A 1 u k ( t ) d x ) ρ 1 + 2 2 ,
which, together with (37), implies directly
Ω | t u k ( t ) | | u k ( t ) | ρ 1 + 1 d x 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + M 15 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ρ 1 + 2 2 , t [ 0 , T * ) ,
in which the positive constant M 15 1 is given by
M 15 1 = ρ 1 + 1 ρ 1 + 2 ( κ ρ 1 + 2 ) ρ 1 + 2 ( ζ 1 ) ρ 1 + 2 2 .
With the aid of the experience of deriving (33)–(36), we conlude that
Ω | t u k ( t ) | | v k ( t ) | ρ 2 + 1 d x 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + ρ 1 + 1 ρ 1 + 2 v k ( t ) L ( ρ 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( Ω ) ( ρ 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + M 16 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , t [ 0 , T * ) ,
in which the positive constant M 16 1 is given by
M 16 1 = ρ 1 + 1 ρ 1 + 2 ( κ ( ρ 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ) ( ρ 2 + 1 ) ( ρ 1 + 2 ) ρ 1 + 1 ( ζ 2 ) ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 .
Combine (33)–(39), to arrive, after some simple calculations, at
Ω t u k ( t ) f 1 ( u k ( t ) , v k ( t ) ) d x 6 Ł 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + M 11 1 Ł 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 12 1 Ł 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 13 1 Ł 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 14 1 Ł 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 15 1 Ł 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ρ 1 + 2 2 + M 16 1 Ł 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , t [ 0 , T * ) .
Mimicking the steps as in deducing (40) from (33)–(39), we could prove similarly that
Ω t v k ( t ) f 2 ( u k ( t ) , v k ( t ) ) d x 6 Ł 2 ρ 2 + 2 t u k ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 + M 21 1 Ł 2 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 22 1 Ł 2 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 23 1 Ł 2 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 24 1 Ł 2 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 25 1 Ł 2 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 26 1 Ł 2 ( Ω v k ( t ) A 2 v k ( t ) d x ) ρ 2 + 2 2 , t [ 0 , T * ) .
in which the positive constants M 21 1 , M 22 1 , M 23 1 , M 24 1 , M 25 1 and M 26 1 are given by
M 21 1 = ρ 2 + 1 ρ 2 + 2 ( κ ( p 1 + 1 ) ( ρ 2 + 2 ) ρ 2 + 1 ) ( p 1 + 1 ) ( ρ 2 + 2 ) ρ 2 + 1 ( ζ 1 ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 ,
M 22 1 = ρ 2 + 1 ρ 2 + 2 ( κ ( p 1 + 1 ) ( ρ 2 + 2 ) ρ 2 + 1 ) ( p 1 + 1 ) ( ρ 2 + 2 ) ρ 2 + 1 ( ζ 2 ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 ,
M 23 1 = ρ 2 + 1 ρ 2 + 2 ( κ ( p 2 + 1 ) ( ρ 2 + 2 ) ρ 2 + 1 ) ( p 2 + 1 ) ( ρ 2 + 2 ) ρ 2 + 1 ( ζ 1 ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 ,
M 24 1 = ρ 2 + 1 ρ 2 + 2 ( κ ( p 2 + 1 ) ( ρ 2 + 2 ) ρ 2 + 1 ) ( p 2 + 1 ) ( ρ 2 + 2 ) ρ 2 + 1 ( ζ 2 ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 ,
M 25 1 = ρ 2 + 1 ρ 2 + 2 ( κ ( ρ 1 + 1 ) ( ρ 2 + 2 ) ρ 2 + 1 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) ρ 2 + 1 ( ζ 1 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2
and
M 26 1 = ρ 2 + 1 ρ 2 + 2 ( κ ρ 2 + 2 ) ρ 2 + 2 ( ζ 2 ) ρ 2 + 2 2 ,
respectively. Plug (26), (27), (28), (29), (40) and (41) into (25) and simplify the obtained result further, to arrive at finally
Ξ k ( t ) M 1 2 + 4 + 6 Ł 1 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + 4 + 6 Ł 2 ρ 2 + 2 t v k ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 + M 11 1 Ł 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 12 1 Ł 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 13 1 Ł 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 14 1 Ł 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 15 1 Ł 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ρ 1 + 2 2 + M 16 1 Ł 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 21 1 Ł 2 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 22 1 Ł 2 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 23 1 Ł 2 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 24 1 Ł 2 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 25 1 Ł 2 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 26 1 Ł 2 ( Ω v k ( t ) A 2 v k ( t ) d x ) ρ 2 + 2 2 + ϰ a 11 , a 12 , a 21 , a 22 ( t u k ( t ) L 2 ( Ω ) 2 + t v k ( t ) L 2 ( Ω ) 2 ) , t [ 0 , T * ) ,
in which the positive constant M 1 2 is given by
M 1 2 = ρ 1 ρ 1 + 2 a 11 L ρ 1 + 2 ρ 1 ( Ω ) ρ 1 + 2 ρ 1 + ρ 1 ρ 2 + ρ 1 + ρ 2 ( ρ 1 + 2 ) ( ρ 2 + 2 ) a 12 L ( ρ 1 + 2 ) ( ρ 2 + 2 ) ρ 1 ρ 2 + ρ 1 + ρ 2 ( Ω ) ( ρ 1 + 2 ) ( ρ 2 + 2 ) ρ 1 ρ 2 + ρ 1 + ρ 2 + ρ 2 ρ 2 + 2 a 22 L ρ 2 + 2 ρ 2 ( Ω ) ρ 2 + 2 ρ 2 + ρ 1 ρ 2 + ρ 1 + ρ 2 ( ρ 1 + 2 ) ( ρ 2 + 2 ) a 21 L ( ρ 1 + 2 ) ( ρ 2 + 2 ) ρ 1 ρ 2 + ρ 1 + ρ 2 ( Ω ) ( ρ 1 + 2 ) ( ρ 2 + 2 ) ρ 1 ρ 2 + ρ 1 + ρ 2 .
With the aid of the definition Ξ k ( t ) (see (24) for the details), by some routine but careful calculations, we arrive at
M 11 1 Ł 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 12 1 Ł 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 M 11 3 ( 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u k ( t ) A 1 u k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 11 3 ( 1 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v k ( t ) A 2 v k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 M 11 3 ( 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u k ( t ) A 1 u k ( t ) d x + 1 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v k ( t ) A 2 v k ( t ) d x ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 M 11 3 Ξ k ( t ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 ,
in which the positive constant M 11 3 is given by
M 11 3 = Ł 1 ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 max ( M 11 1 ( 1 2 ( μ 1 0 + g 1 ( s ) d s ) ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , M 12 1 ( 1 2 ( μ 2 0 + g 2 ( s ) d s ) ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 ) ;
see (11) for the detailed explanation on the notation ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 and see Remark 3 for its applications. Take similar steps as in deriving (43), to get
M 13 1 Ł 1 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 14 1 Ł 1 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 M 12 3 Ξ k ( t ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 ,
in which the positive constant M 12 3 is given by
M 12 3 = Ł 1 ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 max ( M 13 1 ( 1 2 ( μ 1 0 + g 1 ( s ) d s ) ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , M 14 1 ( 1 2 ( μ 2 0 + g 2 ( s ) d s ) ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 ) ;
see (11) for the detailed explanation on the notation ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 and see Remark 3 for its applications. Take similar steps as in deriving (43) and (44), to get
M 21 1 Ł 2 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 22 1 Ł 2 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 M 21 3 Ξ k ( t ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 ,
in which the positive constant M 21 3 is given by
M 21 3 = Ł 2 ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 max ( M 21 1 ( 1 2 ( μ 1 0 + g 1 ( s ) d s ) ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 , M 22 1 ( 1 2 ( μ 2 0 + g 2 ( s ) d s ) ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 ) ;
see (11) for the detailed explanation on the notation ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 and see Remark 3 for its applications. Take similar steps as in deriving (43), (44) and (45), to get
M 23 1 Ł 2 ( Ω u k ( t ) A 1 u k ( t ) d x ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 24 1 Ł 2 ( Ω v k ( t ) A 2 v k ( t ) d x ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 M 22 3 Ξ k ( t ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 ,
in which the positive constant M 22 3 is given by
M 22 3 = Ł 2 ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 max ( M 23 1 ( 1 2 ( μ 1 0 + g 1 ( s ) d s ) ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 , M 24 1 ( 1 2 ( μ 2 0 + g 2 ( s ) d s ) ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 ) ;
see (11) for the detailed explanation on the notation ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 and see Remark 3 for its applications. By the Fenchel–Young inequality, we have
t u k ( t ) L 2 ( Ω ) 2 + t v k ( t ) L 2 ( Ω ) 2 ( ρ 1 ρ 1 + 2 + ρ 2 ρ 2 + 2 ) meas Ω + 2 ρ 1 + 2 t u k ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + 2 ρ 2 + 2 t v k ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 .
Plug (43), (44), (45), (46) and (47) into (42) and perform some simple computations, to arrive at finally the semi-linear differential inequality
Ξ k ( t ) M 1 4 + M 2 4 Ξ k ( t ) + M 3 4 ( Ξ k ( t ) ) ρ 1 + 2 2 + M 4 4 ( Ξ k ( t ) ) ρ 2 + 2 2 + M 5 4 ( Ξ k ( t ) ) ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 6 4 ( Ξ k ( t ) ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 11 3 Ξ k ( t ) ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 12 3 Ξ k ( t ) ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 21 3 Ξ k ( t ) ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 22 3 Ξ k ( t ) ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 , t [ 0 , T * ) ,
in which M 1 4 , M 2 4 , M 3 4 , M 4 4 , M 5 4 and M 6 4 are given by
M 1 4 = M 1 2 + ϰ a 11 , a 12 , a 21 , a 22 ( ρ 1 ρ 1 + 2 + ρ 2 ρ 2 + 2 ) meas Ω ,
M 2 4 = 4 + 6 max ( Ł 1 , Ł 2 ) + 2 ϰ a 11 , a 12 , a 21 , a 22 ,
M 3 4 = M 15 1 Ł 1 ( 1 2 ( μ 1 0 + g 1 ( s ) d s ) ) ρ 1 + 2 2 ,
M 4 4 = M 26 1 Ł 2 ( 1 2 ( μ 2 0 + g 2 ( s ) d s ) ) ρ 2 + 2 2 ,
M 5 4 = M 16 1 Ł 1 ( 1 2 ( μ 2 0 + g 2 ( s ) d s ) ) ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2
and
M 6 4 = M 25 1 Ł 2 ( 1 2 ( μ 1 0 + g 1 ( s ) d s ) ) ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 2 + 2 ,
respectively. By careful calculations, we have
Ξ k ( 0 ) = 1 ρ 1 + 2 u k 0 L ρ 1 + 2 ( Ω ) ρ 1 + 2 + 1 ρ 2 + 2 v k 0 L ρ 2 + 2 ( Ω ) ρ 2 + 2 + 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u k 0 A 1 u k 0 d x + 1 2 Ω u k 1 A 1 u k 1 d x + 1 2 ( g 1 A 1 u k 0 ) ( 0 ) + 1 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v k 0 A 2 v k 0 d x + 1 2 Ω v k 1 A 2 v k 1 d x + 1 2 ( g 2 A 2 v k 0 ) ( 0 ) 1 ρ 1 + 2 ( κ ρ 1 + 2 ) ρ 1 + 2 ( ζ 1 ) ρ 1 + 2 2 ( Ω u k 0 A 1 u k 0 d x ) ρ 1 + 2 2 + 1 ρ 2 + 2 ( κ ρ 2 + 2 ) ρ 2 + 2 ( ζ 2 ) ρ 2 + 2 2 ( Ω v k 0 A 2 v k 0 d x ) ρ 2 + 2 2 + 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u k 0 A 1 u k 0 d x + 1 2 Ω u k 1 A 1 u k 1 d x + 1 2 ( g 1 A 1 u k 0 ) ( 0 ) + 1 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v k 0 A 2 v k 0 d x + 1 2 Ω v k 1 A 2 v k 1 d x + 1 2 ( g 2 A 2 v k 0 ) ( 0 ) 1 ρ 1 + 2 ( κ ρ 1 + 2 ) ρ 1 + 2 ( ζ 1 ) ρ 1 + 2 2 ( Ω u 0 A 1 u 0 d x ) ρ 1 + 2 2 + 1 ρ 2 + 2 ( κ ρ 2 + 2 ) ρ 2 + 2 ( ζ 2 ) ρ 2 + 2 2 ( Ω v 0 A 2 v 0 d x ) ρ 2 + 2 2 + 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u 0 A 1 u 0 d x + 1 2 Ω u 1 A 1 u 1 d x + 1 2 ( g 1 A 1 u 0 ) ( 0 ) + 1 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v 0 A 2 v 0 d x + 1 2 Ω v 1 A 2 v 1 d x + 1 2 ( g 2 A 2 v 0 ) ( 0 ) G 0 ,
where the nonnegative constant G 0 is given by
G 0 = 1 ρ 1 + 2 ( κ ρ 1 + 2 ) ρ 1 + 2 ( ζ 1 tr ( ( A 1 ) 2 ) L ( Ω ) ) ρ 1 + 2 2 u 0 L 2 ( Ω ; R N ) ρ 1 + 2 + 1 ρ 2 + 2 ( κ ρ 2 + 2 ) ρ 2 + 2 ( ζ 2 tr ( ( A 2 ) 2 ) L ( Ω ) ) ρ 2 + 2 2 v 0 L 2 ( Ω ; R N ) ρ 2 + 2 + 1 2 ( μ 1 0 + g 1 ( s ) d s ) tr ( ( A 1 ) 2 ) L ( Ω ) u 0 L 2 ( Ω ; R N ) 2 + 1 2 tr ( ( A 1 ) 2 ) L ( Ω ) u 1 L 2 ( Ω ; R N ) 2 + 0 + g 1 ( s ) d s u 0 L ( R ; L 2 ( Ω ; R N ) ) 2 + 1 2 ( μ 2 0 + g 2 ( s ) d s ) tr ( ( A 2 ) 2 ) L ( Ω ) v 0 L 2 ( Ω ; R N ) 2 + 1 2 tr ( ( A 2 ) 2 ) L ( Ω ) v 1 L 2 ( Ω ; R N ) 2 + 0 + g 2 ( s ) d s v 0 L ( R ; L 2 ( Ω ; R N ) ) 2 .
Let us introduce an auxiliary Cauchy problem
G = M 1 4 + M 2 4 G + M 3 4 G ρ 1 + 2 2 + M 4 4 G ρ 2 + 2 2 + M 5 4 G ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 6 4 G ( ρ 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 11 3 G ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 12 3 G ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 + M 21 3 G ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 + M 22 3 G ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 for t [ 0 , T * ) , G ( 0 ) = G 0 .
Since the right hand side of the differential equation in the problem (49) is smooth (and hence, is locally Lipschitz continuous), by the classical Cauchy–Lipschitz existence theory of ordinary differential equations, the Cauchy problem (49) admits a unique solution, denoted by G ( t ) , in an interval [ 0 , T * * ) , where T * * is a certain constant fulfilling T * * ( 0 , T * ] . It is not difficult to conclude: Both G ( t ) and T * * are uniquely determined by Ω , A 1 , A 2 , f 1 , f 2 , g 1 , g 2 , ρ 1 , ρ 2 , a 11 , a 12 , a 21 , a 22 , u 0 , v 0 , u 1 and v 1 ; G ( t ) is strictly increasing in the interval [ 0 , T * * ) ; and T * * can be chosen arbitrarily in the interval ( 0 , T ˘ ] with T ˘ given by
T ˘ = min ( T * , 1 ⨿ ( G 0 ) α 1 ( α 1 ) )
in which ⨿ = max ( ς , ϶ ) with ς given by
ς = M 1 4 + M 2 4 ( 1 1 α ) + M 3 4 ( 1 ρ 1 + 2 2 α ) + M 4 4 ( 1 ρ 2 + 2 2 α ) + M 5 4 ( 1 ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 α ( ρ 1 + 1 ) ) + M 6 4 ( 1 ( ρ 1 + 1 ) ( ρ 2 + 2 ) 2 α ( ρ 2 + 1 ) ) + M 11 3 ( 1 ( p 1 + 1 ) ( ρ 1 + 2 ) 2 α ( ρ 1 + 1 ) ) + M 12 3 ( 1 ( p 2 + 1 ) ( ρ 1 + 2 ) 2 α ( ρ 1 + 1 ) ) + M 21 3 ( 1 ( p 1 + 1 ) ( ρ 2 + 2 ) 2 α ( ρ 2 + 1 ) ) + M 22 3 ( 1 ( p 2 + 1 ) ( ρ 2 + 2 ) 2 α ( ρ 2 + 1 ) ) ,
϶ given by
϶ = M 2 4 α + M 3 4 ( ρ 1 + 2 ) 2 α + M 4 4 ( ρ 2 + 2 ) 2 α + M 5 4 ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 α ( ρ 1 + 1 ) + M 6 4 ( ρ 1 + 1 ) ( ρ 2 + 2 ) 2 α ( ρ 2 + 1 ) + M 11 3 ( p 1 + 1 ) ( ρ 1 + 2 ) 2 α ( ρ 1 + 1 ) + M 12 3 ( p 2 + 1 ) ( ρ 1 + 2 ) 2 α ( ρ 1 + 1 ) + M 21 3 ( p 1 + 1 ) ( ρ 2 + 2 ) 2 α ( ρ 2 + 1 ) + M 22 3 ( p 2 + 1 ) ( ρ 2 + 2 ) 2 α ( ρ 2 + 1 ) ,
and α given by
α = max ( ρ 1 + 2 2 , ρ 2 + 2 2 , ( ρ 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , ( ρ 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 , ( p 1 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , ( p 2 + 1 ) ( ρ 1 + 2 ) 2 ρ 1 + 2 , ( p 1 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 , ( p 2 + 1 ) ( ρ 2 + 2 ) 2 ρ 2 + 2 ) .
Thanks to 0 Ξ k ( 0 ) G 0 , by standard comparison theory of ordinary differential equations, we deduce from (48) and (49) that
Ξ k ( t ) G ( t ) for t [ 0 , T ˘ ) .
Now let us put T ˜ = 1 2 min ( 2 , T ˘ ) . In light of (50), we have by recalling that the function G ( t ) is strictly increasing in the interval [ 0 , T ˘ ) [ 0 , T ˜ ]
Ξ k ( t ) G ( t ) G ( T ˜ ) for t [ 0 , T ˜ ] ,
and have further by recalling (24) and performing some more elementary calculations
max ( u k ( t ) H 0 1 ( Ω ) 2 , t u k ( t ) H 0 1 ( Ω ) 2 , v k ( t ) H 0 1 ( Ω ) 2 , t v k ( t ) H 0 1 ( Ω ) 2 ) 2 ( 1 + ( κ 2 ) 2 ) M 1 5 Ξ k ( t ) G ( T ˜ ) for t [ 0 , T ˜ ] ,
in which: the positive constant M 1 5 is given by
M 1 5 = min ( 1 ζ 1 ( μ 1 0 + g 1 ( s ) d s ) , 1 ζ 1 , 1 ζ 2 ( μ 2 0 + g 2 ( s ) d s ) , 1 ζ 2 ) ,
and G ( T ˜ ) , a constant in the sense that it is independent of k and t, depends on Ω , A 1 , A 2 , f 1 , f 2 , g 1 , g 2 , ρ 1 , ρ 2 , a 11 , a 12 , a 21 , a 22 , u 0 , v 0 , u 1 and v 1 .
Since ρ i > 0 and p i > 0 (see Assumptions 3 and 4), we have
ρ i + 2 < 2 ( ρ i + 1 ) , p i + 2 < 2 ( p i + 1 ) , i = 1 , 2 .
This, together with Assumptions 3 and 4, and the Rellich–Kondrachov theorem (see [1], Theorem 1, p. 286), implies that the Sobolev embeddings
H 0 1 ( Ω ) L ρ i + 2 ( Ω ) , and H 0 1 ( Ω ) L p i + 2 ( Ω )
are both compact, i = 1 , 2. By applying the Banach–Alaoglu theorem (see [30], Theorem 3.16, p. 66), we can prove, via utilizing (51) and by applying the aforementioned compact embeddings H 0 1 ( Ω ) L ρ i + 2 ( Ω ) and H 0 1 ( Ω ) L p i + 2 ( Ω ) , i = 1 , 2, that there exists a pair ( u , v ) (whose restriction to Ω × [ 0 , T ˜ ] could be proved to be weakly continuous with respect to time t in the Hilbert space H 0 1 ( Ω ; R 2 ) ), and a subsequence { i k } of { k } (that is, a strictly increasing sequence in N ) such that
u i k u weakly * in L ( 0 , T ˜ ; H 0 1 ( Ω ) ) as k , v i k v weakly * in L ( 0 , T ˜ ; H 0 1 ( Ω ) ) as k , t u i k t u weakly * in L ( 0 , T ˜ ; H 0 1 ( Ω ) ) as k , t v i k t v weakly * in L ( 0 , T ˜ ; H 0 1 ( Ω ) ) as k ,
and such that
| t u i k ( t ) | ρ 1 t u i k ( t ) | t u ( t ) | ρ 1 t u ( t ) in L ρ 1 + 2 ρ 1 + 1 ( Ω ) as k , | t v i k ( t ) | ρ 2 t v i k ( t ) | t v ( t ) | ρ 2 t v ( t ) as k , f 1 ( u i k ( t ) , v i k ( t ) ) f 1 ( u ( t ) , v ( t ) ) in L min ( ρ 1 + 2 ρ 1 + 1 , ρ 2 + 2 ρ 2 + 1 , p 1 + 2 p 1 + 1 , p 2 + 2 p 2 + 1 ) ( Ω ) as k , f 2 ( u i k ( t ) , v i k ( t ) ) f 2 ( u ( t ) , v ( t ) ) in L min ( ρ 1 + 2 ρ 1 + 1 , ρ 2 + 2 ρ 2 + 1 , p 1 + 2 p 1 + 1 , p 2 + 2 p 2 + 1 ) ( Ω ) as k .
In addition, for every pair ( φ , ψ ) H 0 1 ( Ω ; R 2 ) , it holds that
1 ρ 1 + 1 Ω | t u i k ( t ) | ρ 1 t u i k ( t ) φ d x + Ω φ A 1 t u i k ( t ) d x 1 ρ 1 + 1 Ω | j = 1 i k Ω u 1 e 1 j d x e 1 j | ρ 1 j = 1 i k Ω u 1 e 1 j d x e 1 j φ d x Ω φ A 1 ( j = 1 i k Ω u 1 e 1 j d x e 1 j ) d x + 0 t Ω φ A 1 ( u i k ( s ) s g 1 ( s τ ) u i k ( τ ) d τ ) d x d s = 0 t Ω φ ( f 1 ( u i k ( s ) , v i k ( s ) ) a 11 t u i k ( s ) a 12 t v i k ( s ) ) d x d s , t [ 0 , T ˜ ] ,
and
1 ρ 2 + 1 Ω | t v i k ( t ) | ρ 2 t v i k ( t ) ψ d x + Ω ψ A 2 t v i k ( t ) d x 1 ρ 2 + 1 Ω | j = 1 i k Ω v 1 e 2 j d x e 2 j | ρ 2 j = 1 i k Ω v 1 e 2 j d x e 2 j ψ d x Ω ψ A 2 ( j = 1 i k Ω v 1 e 2 j d x e 2 j ) d x + 0 t Ω ψ A 2 ( v i k ( s ) s g 2 ( s τ ) v i k ( τ ) d τ ) d x d s = 0 t Ω ψ ( f 2 ( u i k ( s ) , v i k ( s ) ) a 21 t u i k ( s ) a 22 t v i k ( s ) ) d x d s , t [ 0 , T ˜ ] .
By recalling that
j = 1 i k Ω u 0 ( t ) e 1 j d x e 1 j u 0 in L ( R ; H 0 1 ( Ω ) ) as k , j = 1 i k Ω v 0 ( t ) e 2 j d x e 2 j v 0 in L ( R ; H 0 1 ( Ω ) ) as k , j = 1 i k Ω u 1 e 1 j d x e 1 j u 1 in H 0 1 ( Ω ) as k , j = 1 i k Ω v 1 e 2 j d x e 2 j v 1 in H 0 1 ( Ω ) as k ,
we could conclude, based on the idea of passing to the limit of (53) and (54), that ( u , v ) , the limit of ( u i k , v i k ) , satisfies (16), (17), (1) 3 and (1) 4 .
Lastly, we can mimick steps in [9,23], to show that ( u , v ) S [ 0 , T ˜ ] . To summarize, ( u , v ) , the limit of ( u i k , v i k ) , is indeed a local solution, in the interval [ 0 , T ˜ ] , to IBVP (1) in the sense of Definition 1. The proof is complete. □
Remark 5.
Illuminated by the integral identity
0 t Ω φ A 1 t 2 u i k ( s ) d x d s = Ω φ A 1 t u i k ( t ) d x Ω φ A 1 ( j = 1 i k Ω u 1 e 1 j d x e 1 j ) d x , t [ 0 , T ˜ ] , φ H 0 1 ( Ω ) ,
and by the integral identity
0 t Ω ψ A 2 t 2 v i k ( s ) d x d s = Ω ψ A 2 t v i k ( t ) d x Ω ψ A 2 ( j = 1 i k Ω v 1 e 2 j d x e 2 j ) d x , t [ 0 , T ˜ ] , ψ H 0 1 ( Ω ) ,
we conclude that
t 2 u i k t 2 u weakly * in D ( 0 , T ˜ ; H 0 1 ( Ω ) ) , as k , and t 2 v i k t 2 v weakly * in D ( 0 , T ˜ ; H 0 1 ( Ω ) ) , as k .
Besides, enlightened by the integral identity
0 t Ω | t u i k ( s ) | ρ 1 t 2 u i k ( s ) φ d x d s = 1 ρ 1 + 1 Ω | t u i k ( t ) | ρ 1 t u i k ( t ) φ d x 1 ρ 1 + 1 Ω | j = 1 i k Ω u 1 e 1 j d x e 1 j | ρ 1 j = 1 i k Ω u 1 e 1 j d x e 1 j φ d x , t [ 0 , T ˜ ] , φ L 2 ( Ω ) ,
and by the integral identity
0 t Ω | t v i k ( s ) | ρ 2 t 2 v i k ( s ) ψ d x d s = 1 ρ 2 + 1 Ω | t v i k ( t ) | ρ 2 t v i k ( t ) ψ d x 1 ρ 2 + 1 Ω | j = 1 i k Ω v 1 e 2 j d x e 2 j | ρ 2 j = 1 i k Ω v 1 e 2 j d x e 2 j ψ d x , t [ 0 , T ˜ ] , ψ L 2 ( Ω ) ,
we conclude that
| t u i k | ρ 1 t 2 u i k | t u | ρ 1 t 2 u weakly * in D ( 0 , T ˜ ; L 2 ( Ω ) ) , as k , and | t v i k | ρ 2 t 2 v i k | t v | ρ 2 t 2 v weakly * in D ( 0 , T ˜ ; L 2 ( Ω ) ) , as k .
To sum up, the solution ( u , v ) S [ 0 , T ˜ ] (see (14) for the definition of S [ 0 , T ˜ ] ) in the sense of Definition 1 whose existence justified by Theorem 1 satisfies automatically
t 2 u , t 2 v D ( 0 , T ˜ ; H 0 1 ( Ω ) ) , and | t u | ρ 1 t 2 u , | t v | ρ 2 t 2 v D ( 0 , T ˜ ; L 2 ( Ω ) ) .
Remark 6.
By re-checking the proof of Theorem 1, we may find that the restriction on the symmetric matrices-valued A i function, defined in the domain Ω, (see Assumption 1 for the details) could be weakened to: A i W 1 , ( Ω ) satisfies
ess inf x Ω , ξ R N { 0 } ξ A i ( x ) ξ ξ ξ > 0 , i = 1 , 2 ,
and may find that the restriction that a i j C ( Ω ¯ ) on the coefficient a i j in Assumption 5 could be weakened to a i j L ( Ω ) , i, j = 1 , 2 . For the sake of convenience of our calculations, unless stated otherwise, we abide by Assumptions 1 and 5 in the rest of the paper.
Theorem 2.
Suppose that Assumptions 1–5 hold true. For every pair
( u 0 , v 0 ) L ( R ; H 0 1 ( Ω ; R 2 ) )
and every pair ( u 1 , v 1 ) H 0 1 ( Ω ; R 2 ) , IBVP (1) admits a solution ( u , v ) S [ 0 , T ) in the sense of Definition 2, in which [ 0 , T ) , is the maximal existence time interval of the solution pair ( u , v ) with the maximal existence time instant 0 < T + independent of ( u , v ) but depending on Ω, A 1 , A 2 , f 1 , f 2 , g 1 , g 2 , ρ 1 , ρ 2 , a 11 , a 12 , a 21 , a 22 , u 0 , v 0 , u 1 and v 1 .
Proof. 
Theorem 2 can be proved by a standard continuation procedure. And therefore we choose to leave out the detailed steps in this paper. □
Lemma 3.
Suppose that Assumptions 1–5 hold true. E u , v ( t ) defined as in (18), associated to every pair ( u , v ) of solution, in the interval [ 0 , T ) , to IBVP (1) in the sense of Definition 2 (see Theorem 2 for the existence of ( u , v ) ), is non-increasing in [ 0 , T ) .
Proof. 
Mimicking steps conducted in (25), we have
d d t E u , v ( t ) = 1 2 ( g 1 A 1 u ) ( t ) + 1 2 ( g 2 A 2 v ) ( t ) Ω t u ( t ) t v ( t ) a 11 a 12 a 21 a 22 t u ( t ) t v ( t ) d x , t [ 0 , T ) ,
which, together with Assumptions 2 and 5, implies immediately that E u , v ( t ) is indeed non-increasing in the interval [ 0 , T ) . The proof is complete. □
Theorem 2, having Theorem 1 as its basis, states that to every initial-datum pair
( u 0 , v 0 ) L ( R ; H 0 1 ( Ω ; R 2 ) )
and every initial-datum pair ( u 1 , v 1 ) H 0 1 ( Ω ; R 2 ) , there corresponds a weak solution pair ( u , v ) S [ 0 , T ) in the sense of Definition 2 to IBVP (1), in which [ 0 , T ) is the maximal existence time interval of the solution pair ( u , v ) . Now it is natural to start to investigate the global existence of solutions to IBVP (1).
To make it convenient to present our other results in the rest of this paper, let us introduce two auxilliary functionals
I ( t ) = 1 4 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + 1 4 Ω t u ( t ) A 1 t u ( t ) d x + 1 4 ( g 1 A 1 u ) ( t ) + 1 4 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + 1 4 Ω t v ( t ) A 2 t v ( t ) d x + 1 4 ( g 2 A 2 v ) ( t ) Ω F ( u ( t ) , v ( t ) ) d x , t [ 0 , T )
and
J ( t ) = 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + 1 2 Ω t u ( t ) A 1 t u ( t ) d x + 1 2 ( g 1 A 1 u ) ( t ) + 1 2 ( μ 1 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + 1 2 Ω t v ( t ) A 2 t v ( t ) d x + 1 2 ( g 2 A 2 v ) ( t ) Ω F ( u ( t ) , v ( t ) ) d x , t [ 0 , T ) .
We shall see below that our global existence and general energy decay results hold true only for solutions to IBVP (1) having small initial data. To measure such smallness, we need introduce the following two constants
β 1 = 4 M 11 6 Ł 1 + Ł 2 ( p 1 + 1 ) p 1 + 2 + 4 M 12 6 Ł 1 + Ł 2 ( p 2 + 1 ) p 2 + 2 + 4 M 13 6 Ł 1 + Ł 2 ( ρ 1 + 1 ) ρ 1 + 2 ,
and
β 2 = 8 Ł 2 M 21 6 p 1 + 2 + 8 Ł 2 M 22 6 p 2 + 2 + 4 Ł 2 M 23 6 ρ 1 + 2 + 4 Ł 2 M 24 6 ρ 2 + 2 ,
in which the constants M 11 6 , M 12 6 , M 13 6 , M 14 6 , M 21 6 , M 22 6 , M 23 6 and M 24 6 are given by
M 11 6 = ( ζ 1 ( κ p 1 + 2 ) 2 μ 1 0 + g 1 ( s ) d s ) p 1 + 2 2 ( 4 E u , v ( 0 ) ) p 1 2 ,
M 12 6 = ( ζ 1 ( κ p 2 + 2 ) 2 μ 1 0 + g 1 ( s ) d s ) p 2 + 2 2 ( 4 E u , v ( 0 ) ) p 2 2 ,
M 13 6 = ( ζ 1 ( κ ρ 1 + 2 ) 2 μ 1 0 + g 1 ( s ) d s ) ρ 1 + 2 2 ( 4 E u , v ( 0 ) ) ρ 1 2 ,
M 14 6 = ( ζ 1 ( κ ρ 2 + 2 ) 2 μ 1 0 + g 1 ( s ) d s ) ρ 2 + 2 2 ( 4 E u , v ( 0 ) ) ρ 2 2 ,
M 21 6 = ( ζ 2 ( κ p 1 + 2 ) 2 μ 2 0 + g 2 ( s ) d s ) p 1 + 2 2 ( 4 E u , v ( 0 ) ) p 1 2 ,
M 22 6 = ( ζ 2 ( κ p 2 + 2 ) 2 μ 2 0 + g 2 ( s ) d s ) p 2 + 2 2 ( 4 E u , v ( 0 ) ) p 2 2 ,
M 23 6 = ( ζ 2 ( κ ρ 1 + 2 ) 2 μ 2 0 + g 2 ( s ) d s ) ρ 1 + 2 2 ( 4 E u , v ( 0 ) ) ρ 1 2 ,
and
M 24 6 = ( ζ 2 ( κ ρ 2 + 2 ) 2 μ 2 0 + g 2 ( s ) d s ) ρ 2 + 2 2 ( 4 E u , v ( 0 ) ) ρ 2 2 ,
respectively, where the constant ζ i is given as in (9), i = 1 , 2, and E u , v ( 0 ) = E u , v ( t ) | t = 0 can be formulated explicitly as
E u , v ( 0 ) = 1 ρ 1 + 2 u 1 L ρ 1 + 2 ( Ω ) ρ 1 + 2 + 1 ρ 2 + 2 v 1 L ρ 2 + 2 ( Ω ) ρ 2 + 2 + 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u 0 ( 0 ) A 1 u 0 ( 0 ) d x + 1 2 Ω u 1 A 1 u 1 d x + 1 2 ( g 1 A 1 u 0 ) ( 0 ) + 1 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v 0 ( 0 ) A 2 v 0 ( 0 ) d x + 1 2 Ω v 1 A 2 v 1 d x + 1 2 ( g 2 A 2 v 0 ) ( 0 ) Ω F ( u 0 , v 0 ) d x .
Please consult (18) for the detailed expression of E u , v ( t ) .
Lemma 4.
Suppose that Assumptions 1–5 hold true. For every weak solution pair ( u , v ) S [ 0 , T ) to IBVP (1) with [ 0 , T ) the maximal existence interval of ( u , v ) , if the associated functional I ( t ) given by (56) satisfies I ( 0 ) > 0 and the associated constants β 1 and β 2 , given by (58) and (59), respectively, satisfy max ( β 1 , β 2 ) < 1 , then I ( t ) > 0 holds for all t [ 0 , T ) .
Proof. 
Thanks to the assumption that I ( 0 ) > 0 and to the continuity of the function I ( t ) , there exists a time T 1 ( 0 , T ) , such that
I ( t ) 0 , t [ 0 , T 1 ] .
By the very definition of I ( t ) (see (56) for the details), we have
J ( t ) = I ( t ) + 1 4 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + 1 4 Ω t u ( t ) A 1 t u ( t ) d x + 1 4 ( g 1 A 1 u ) ( t ) + 1 4 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + 1 4 Ω t v ( t ) A 2 t v ( t ) d x + 1 4 ( g 2 A 2 v ) ( t ) , t [ 0 , T 1 ] ,
where the functional J ( t ) , associated to the solution ( u , v ) to IBVP (1), is defined as in (57). Combine (68) and (69), to arrive at directly
1 4 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + 1 4 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x J ( t ) , t [ 0 , T 1 ] .
In light of the definition (57) of the functional J ( t ) , we have
J ( t ) E u , v ( t ) E u , v ( 0 ) , t [ 0 , T 1 ] ,
where the second ‘⩽’ follows from the non-increasing monotonicity (see Lemma 3 for the details). Substitute (71) into (70), to arrive at immediately
( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x 4 E u , v ( 0 ) , t [ 0 , T 1 ] .
By Lemma 1, we have
| Ω F ( u ( t ) , v ( t ) ) d x | Ł 1 + Ł 2 ( p 1 + 1 ) p 1 + 2 u ( t ) L p 1 + 2 ( Ω ) p 1 + 2 + 2 Ł 2 p 1 + 2 v ( t ) L p 1 + 2 ( Ω ) p 1 + 2 + Ł 1 + Ł 2 ( p 2 + 1 ) p 2 + 2 u ( t ) L p 2 + 2 ( Ω ) p 2 + 2 + 2 Ł 2 p 2 + 2 v ( t ) L p 2 + 2 ( Ω ) p 2 + 2 + Ł 1 + Ł 2 ( ρ 1 + 1 ) ρ 1 + 2 u ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + Ł 2 ρ 1 + 2 v ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + Ł 2 ρ 2 + 2 v ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 , t [ 0 , T 1 ] .
Having the notations in (9) and (10) at our disposal and based on the coercivity condition on A 1 (see Assumption 1), mimicking steps in deducing (32), we apply the Sobolev–Poincaré inequality (see Remark 1), to conclude
u ( t ) L p 1 + 2 ( Ω ) p 1 + 2 M 11 6 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x .
And similarly, we have
u ( t ) L p 2 + 2 ( Ω ) p 2 + 2 M 12 6 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x ,
u ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 M 13 6 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x ,
and
u ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 M 14 6 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x .
Having the notations in (9) and (10) at our disposal and based on the coercivity condition on A 2 (see Assumption 1), mimicking steps in deducing (32) and (74)–(77), we apply the Sobolev–Poincaré inequality (see Remark 1), to obtain
v ( t ) L p 1 + 2 ( Ω ) p 1 + 2 M 21 6 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x .
And analogously, we have
v ( t ) L p 2 + 2 ( Ω ) p 2 + 2 M 22 6 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x ,
v ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 M 23 6 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x ,
and
v ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 M 24 6 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x .
Plug (74)–(81) into (73), to obtain
| Ω F ( u ( t ) , v ( t ) ) d x | β 1 4 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + β 2 4 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x , t [ 0 , T 1 ] ,
where the constants β 1 and β 2 are defined as in (58) and (59), respectively. Owing to the assumption that max ( β 1 , β 2 ) < 1 , it follows that
| Ω F ( u ( t ) , v ( t ) ) d x | < 1 4 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + 1 4 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x , t [ 0 , T 1 ] .
This, together with the definition (56) of I ( t ) , implies that
I ( t ) > 0 for t [ 0 , T 1 ] .
Lastly, to finish the proof, we introduce
T 2 = sup { 0 < T T ; I ( t ) > 0 , t [ 0 , T ) } .
Obviously, T 2 does not exceed T. By a contradiction argument, we shall show that T 2 coincides actually with T. We assume to the contrary that T 2 < T , then by the definition of T 2 as well as the continuity of I ( t ) , we have immediately that I ( T 2 ) 0 . Enlightened by the procedure of deriving (83) from (68), we have therefore I ( T 2 ) > 0 . Thanks to the continuity of I ( t ) , there exists a T 2 < T 3 < T such that I ( t ) > 0 holds for all t [ 0 , T 3 ] . This contradicts the definition of T 2 . This implies indeed that T 2 coincides actually with T. In other words, I ( t ) > 0 holds for all t [ 0 , T ) . □
Theorem 3.
Suppose that Assumptions 1–5 hold true. For every weak solution pair ( u , v ) S [ 0 , T ) to IBVP (1) with [ 0 , T ) the maximal existence interval of ( u , v ) , if the associated functional I ( t ) given by (56) satisfies I ( 0 ) > 0 and the associated constants β 1 and β 2 , given by (58) and (59), respectively, satisfy max ( β 1 , β 2 ) < 1 , then T = + . In other words, weak solutions to IBVP (1) exist globally in time whenever their initial data satisfy I ( 0 ) > 0 and max ( β 1 , β 2 ) < 1 .
Proof. 
We shall prove, by a contradiction argument, that T is exactly + . Let us assume that T < + . Following the idea used to obtain (72) in the proof of Lemma 4, we have
( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + Ω t u ( t ) A 1 t u ( t ) d x + ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + Ω t v ( t ) A 2 t v ( t ) d x 4 I ( t ) + ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + Ω t u ( t ) A 1 t u ( t ) d x + ( g 1 A 1 u ) ( t ) + ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + Ω t v ( t ) A 2 t v ( t ) d x + ( g 2 A 2 v ) ( t ) = 4 J ( t ) 4 E u , v ( t ) 4 E u , v ( 0 ) , t [ 0 , T ) .
This, together with (9) (see Remark 1 for the details), implies that
max ( u ( t ) L 2 ( Ω ; R N ) 2 , t u ( t ) L 2 ( Ω ; R N ) 2 , v ( t ) L 2 ( Ω ; R N ) 2 , t v ( t ) L 2 ( Ω ; R N ) 2 ) 4 M 1 5 E u , v ( 0 ) , t [ 0 , T ) ,
which, together with (10) (see also Remark 1 for the details), implies further
max ( u ( t ) H 0 1 ( Ω ) 2 , t u ( t ) H 0 1 ( Ω ) 2 , v ( t ) H 0 1 ( Ω ) 2 , t v ( t ) H 0 1 ( Ω ) 2 ) 4 ( 1 + ( κ 2 ) 2 ) M 1 5 E u , v ( 0 ) , t [ 0 , T ) ,
where the positive constant M 1 5 , given by (52), is independent of the time variable t. Now, by a standard continution procedure, we could obtain a time T 4 > T such that IBVP (1) admits a solution in [ 0 , T 4 ) . This contradicts the assumption that [ 0 , T ) is the maximal existence interval of IBVP (1). This implies immediately that [ 0 , T ) , the maximal existence time interval of the solution ( u , v ) , coincides actually with R + . In other words, the solution ( u , v ) exists globally in time. The proof is complete. □

3. General Energy Decay Results Concerning Solutions to IBVP (1)

Let us now associate to every solution ( u , v ) to IBVP (1) the functional
L 1 ( t ) = 1 ρ 1 + 1 Ω u ( t ) | t u ( t ) | ρ 1 t u ( t ) d x + Ω u ( t ) A 1 t u ( t ) d x + 1 ρ 2 + 1 Ω v ( t ) | t v ( t ) | ρ 2 t v ( t ) d x + Ω v ( t ) A 2 t v ( t ) d x + 1 2 Ω u ( t ) v ( t ) a 11 a 12 a 21 a 22 u ( t ) v ( t ) d x .
As we shall see, L 1 ( t ) plays a role of the energy perturbation, and one differential inequality concerning the functional L 1 ( t ) plays an important role in proving our main results in this paper. This inequality can only be established for those solutions whose initial data are small. To give a precise sense by which we mean the smallness, we need to introduce two useful constants 1 and 2 :
1 = ( 4 E u , v ( 0 ) ) p 1 2 ( Ł 1 + Ł 1 p 1 + 2 + Ł 2 ( p 1 + 1 ) p 1 + 2 ) ( ζ 1 ( κ p 1 + 2 ) 2 μ 1 0 + g 1 ( s ) d s ) p 1 + 2 2 + ( 4 E u , v ( 0 ) ) p 2 2 ( Ł 1 + Ł 1 p 2 + 2 + Ł 2 ( p 2 + 1 ) p 2 + 2 ) ( ζ 1 ( κ p 2 + 2 ) 2 μ 1 0 + g 1 ( s ) d s ) p 2 + 2 2 + ( 4 E u , v ( 0 ) ) ρ 1 2 ( Ł 1 + Ł 2 ( ρ 1 + 1 ) ρ 1 + 2 ) ( ζ 1 ( κ ρ 1 + 2 ) 2 μ 1 0 + g 1 ( s ) d s ) ρ 1 + 2 2 + Ł 1 ( ρ 2 + 1 ) ( 4 E u , v ( 0 ) ) ρ 2 2 ρ 2 + 2 ( ζ 1 ( κ ρ 2 + 2 ) 2 μ 1 0 + g 1 ( s ) d s ) ρ 2 + 2 2 ,
2 = ( 4 E u , v ( 0 ) ) p 1 2 ( Ł 2 + Ł 2 p 1 + 2 + Ł 1 ( p 1 + 1 ) p 1 + 2 ) ( ζ 2 ( κ p 1 + 2 ) 2 μ 2 0 + g 2 ( s ) d s ) p 1 + 2 2 + ( 4 E u , v ( 0 ) ) p 2 2 ( Ł 2 + Ł 2 p 2 + 2 + Ł 1 ( p 2 + 1 ) p 2 + 2 ) ( ζ 2 ( κ p 2 + 2 ) 2 μ 2 0 + g 2 ( s ) d s ) p 2 + 2 2 + Ł 2 ( 4 E u , v ( 0 ) ) ρ 1 2 ρ 1 + 2 ( ζ 2 ( κ ρ 1 + 2 ) 2 μ 2 0 + g 2 ( s ) d s ) ρ 1 + 2 2 + ( 4 E u , v ( 0 ) ) ρ 2 2 ( Ł 2 + Ł 1 ( ρ 2 + 1 ) ρ 2 + 2 ) ( ζ 2 ( κ ρ 2 + 2 ) 2 μ 2 0 + g 2 ( s ) d s ) ρ 2 + 2 2 .
Lemma 5.
Suppose that Assumptions 1–5 hold true. If the associated functional I ( t ) given by (56) satisfies I ( 0 ) > 0 , the associated constants β 1 and β 2 , given by (58) and (59), respectively, satisfy max ( β 1 , β 2 ) < 1 , and the associated constants 1 and 2 , given by (86) and (87), respectively, satisfy max ( 1 , 2 ) < 1 4 , then weak solutions to IBVP (1) exist globally in time, and render the associated functional L 1 ( t ) given by (85) to satisfy
L 1 ( t ) 1 ρ 1 + 1 Ω | t u ( t ) | ρ 1 + 2 d x + 1 ρ 2 + 1 Ω | t v ( t ) | ρ 2 + 2 d x + Ω t u ( t ) A 1 t u ( t ) d x + Ω t v ( t ) A 2 t v ( t ) d x 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x 1 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + 1 μ 1 0 + g 1 ( s ) d s 0 + ( g 1 ( s ) ) 2 g 1 ( s ) G g 1 ( s ) d s ( ( g 1 G g 1 ) A 1 u ) ( t ) + 1 μ 2 0 + g 2 ( s ) d s 0 + ( g 2 ( s ) ) 2 g 2 ( s ) G g 2 ( s ) d s ( ( g 2 G g 2 ) A 2 v ) ( t ) , t R + ,
where G is a positive constant yet to be determined later.
Proof. 
As mentioned in Lemma 5, thanks to the assumptions I ( 0 ) > 0 and max ( β 1 , β 2 ) < 1 , it follows from Theorem 3 that L 1 ( t ) exists globally in time. Recalling (85), we differentiate L 1 ( t ) , to arrive at
L 1 ( t ) = 1 ρ 1 + 1 Ω | t u ( t ) | ρ 1 + 2 d x + 1 ρ 2 + 1 Ω | t v ( t ) | ρ 2 + 2 d x + Ω t u ( t ) A 1 t u ( t ) d x + Ω t v ( t ) A 2 t v ( t ) d x μ 1 Ω u ( t ) A 1 u ( t ) d x μ 2 Ω v ( t ) A 2 v ( t ) d x + t g 1 ( t s ) Ω u ( t ) A 1 u ( s ) d x d s + t g 2 ( t s ) Ω v ( t ) A 2 v ( s ) d x d s + Ω ( u ( t ) f 1 ( u ( t ) , v ( t ) ) + v ( t ) f 2 ( u ( t ) , v ( t ) ) ) d x , t R + .
By some routine calculations, we have
t g 1 ( t s ) Ω u ( t ) A 1 u ( s ) d x d s μ 1 Ω u ( t ) A 1 u ( t ) d x = ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x t g 1 ( t s ) Ω u ( t ) A 1 ( u ( t ) u ( s ) ) d x d s , t R + .
Since A 1 satisfies the coercivity condition that ϑ 1 > 0 with ϑ 1 given by (5) (see Assumption 1 for the details), from the Cauchy–Schwarz inequality it follows that
t g 1 ( t s ) Ω u ( t ) A 1 ( u ( t ) u ( s ) ) d x d s t Ω ( g 1 ( t s ) ) 2 g 1 ( t s ) G g 1 ( t s ) u ( t ) A 1 u ( t ) d x d s 1 2 · t Ω g 1 ( t s ) G g 1 ( t s ) ( u ( t ) u ( s ) ) A 1 ( u ( t ) u ( s ) ) d x d s 1 2 = Ω u ( t ) A 1 u ( t ) d x 1 2 0 + ( g 1 ( s ) ) 2 g 1 ( s ) G g 1 ( s ) d s ( ( g 1 G g 1 ) A 1 u ) ( t ) 1 2 μ 1 0 + g 1 ( s ) d s 4 Ω u ( t ) A 1 u ( t ) d x + 1 μ 1 0 + g 1 ( s ) d s 0 + ( g 1 ( s ) ) 2 g 1 ( s ) G g 1 ( s ) d s ( ( g 1 G g 1 ) A 1 u ) ( t ) , t R + ,
which, together with (90), implies
t g 1 ( t s ) Ω u ( t ) A 1 u ( s ) d x d s μ 1 Ω u ( t ) A 1 u ( t ) d x 3 4 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + 1 μ 1 0 + g 1 ( s ) d s 0 + ( g 1 ( s ) ) 2 g 1 ( s ) G g 1 ( s ) d s ( ( g 1 G g 1 ) A 1 u ) ( t ) , t R + ,
where G is a sufficiently large positive constant. Since A 2 satisfies the coercivity condition that ϑ 2 > 0 with ϑ 2 given by (5) (see Assumption 1 for the details), it follows by mimicking steps in deducing (91) that
t g 2 ( t s ) Ω v ( t ) A 2 v ( s ) d x d s μ 2 Ω v ( t ) A 2 v ( t ) d x 3 4 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + 1 μ 2 0 + g 2 ( s ) d s 0 + ( g 2 ( s ) ) 2 g 2 ( s ) G g 2 ( s ) d s ( ( g 2 G g 2 ) A 2 v ) ( t ) , t R + .
By recalling Assumption 4 (especially (8)), by recalling the experience of deriving (31), (33)–(36), and by utilizing mainly the Fenchel–Young inequality, we have
Ω ( u ( t ) f 1 ( u ( t ) , v ( t ) ) + v ( t ) f 2 ( u ( t ) , v ( t ) ) ) d x Ł 1 Ω ( | u ( t ) | p 1 + 2 + | u ( t ) | p 2 + 2 + | u ( t ) | ρ 1 + 2 ) d x + Ł 1 Ω | u ( t ) | | v ( t ) | p 1 + 1 d x + Ł 1 Ω | u ( t ) | | v ( t ) | p 2 + 1 d x + Ł 1 Ω | u ( t ) | | v ( t ) | ρ 2 + 1 d x + Ł 2 Ω ( | v ( t ) | p 1 + 2 + | v ( t ) | p 2 + 2 + | v ( t ) | ρ 2 + 2 ) d x + Ł 2 Ω | v ( t ) | | u ( t ) | p 1 + 1 d x + Ł 2 Ω | v ( t ) | | u ( t ) | p 2 + 1 d x + Ł 2 Ω | v ( t ) | | u ( t ) | ρ 1 + 1 d x ( Ł 1 + Ł 1 p 1 + 2 + Ł 2 ( p 1 + 1 ) p 1 + 2 ) u ( t ) L p 1 + 2 ( Ω ) p 1 + 2 + ( Ł 1 + Ł 1 p 2 + 2 + Ł 2 ( p 2 + 1 ) p 2 + 2 ) u ( t ) L p 2 + 2 ( Ω ) p 2 + 2 + ( Ł 1 + Ł 2 ( ρ 1 + 1 ) ρ 1 + 2 ) u ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + Ł 1 ( ρ 2 + 1 ) ρ 2 + 2 u ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 + ( Ł 2 + Ł 2 p 1 + 2 + Ł 1 ( p 1 + 1 ) p 1 + 2 ) v ( t ) L p 1 + 2 ( Ω ) p 1 + 2 + ( Ł 2 + Ł 2 p 2 + 2 + Ł 1 ( p 2 + 1 ) p 2 + 2 ) v ( t ) L p 2 + 2 ( Ω ) p 2 + 2 + Ł 2 ρ 1 + 2 v ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 + ( Ł 2 + Ł 1 ( ρ 2 + 1 ) ρ 2 + 2 ) v ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 , t R + .
With the aid of (74), (75), (76), (77), (78), (79), (80) and (81), we combine (89), (91), (92) and (93), to obtain
L 1 ( t ) 1 ρ 1 + 1 Ω | t u ( t ) | ρ 1 + 2 d x + 1 ρ 2 + 1 Ω | t v ( t ) | ρ 2 + 2 d x + Ω t u ( t ) A 1 t u ( t ) d x + Ω t v ( t ) A 2 t v ( t ) d x + ( 1 3 4 ) ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + ( 2 3 4 ) ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + 1 μ 1 0 + g 1 ( s ) d s 0 + ( g 1 ( s ) ) 2 g 1 ( s ) G g 1 ( s ) d s ( ( g 1 G g 1 ) A 1 u ) ( t ) + 1 μ 2 0 + g 2 ( s ) d s 0 + ( g 2 ( s ) ) 2 g 2 ( s ) G g 2 ( s ) d s ( ( g 2 G g 2 ) A 2 v ) ( t ) , t R + .
This, together with the assumption
max ( 1 , 2 ) < 1 4 ,
implies that the proof of Lemma 5 is complete. □
As with L 1 ( t ) defined by (85), the following energy perturbation functional associated to each solution pair ( u , v ) to IBVP (1) is of great importance:
L 2 ( t ) = Ω div ( A 1 t u ( t ) ) | t u ( t ) | ρ 1 t u ( t ) ρ 1 + 1 t g 1 ( t s ) u ( t ) u ( s ) d s d x + Ω div ( A 2 t v ( t ) ) | t v ( t ) | ρ 2 t v ( t ) ρ 2 + 1 t g 2 ( t s ) v ( t ) v ( s ) d s d x .
As with what we did in Lemma 5 for L 1 ( t ) , we shall establish a useful differential inequality for L 2 ( t ) in the following lemma.
Lemma 6.
Suppose that Assumptions 1–5 hold true. If the associated functional I ( t ) given by (56) satisfies I ( 0 ) > 0 , the associated constants β 1 and β 2 , given by (58) and (59), respectively, satisfy max ( β 1 , β 2 ) < 1 , then weak solutions to IBVP (1) exist globally in time, and render the associated functional L 2 ( t ) given by (94) to satisfy
L 2 ( t ) δ 0 + g 1 ( s ) d s Ω u ( t ) A 1 u ( t ) d x + δ 0 + g 2 ( s ) d s Ω v ( t ) A 2 v ( t ) d x ( 0 + g 1 ( s ) d s δ g 1 ( 0 ) ) Ω t u ( t ) A 1 t u ( t ) d x ( 0 + g 2 ( s ) d s δ g 2 ( 0 ) ) Ω t v ( t ) A 2 t v ( t ) d x 1 ρ 1 + 1 0 + g 1 ( s ) d s Ω | t u ( t ) | ρ 1 + 2 d x 1 ρ 2 + 1 0 + g 2 ( s ) d s Ω | t v ( t ) | ρ 2 + 2 d x + Ω t u ( t ) t v ( t ) a 11 a 12 a 21 a 22 t u ( t ) t v ( t ) d x + W 11 ( g 1 A 1 u ) ( t ) + W 21 ( g 2 A 2 v ) ( t ) W 12 ( g 1 A 1 u ) ( t ) W 22 ( g 2 A 2 v ) ( t ) , t R + ,
where the positive constant δ is given arbitrarily, the positive constants W 11 , W 12 , W 21 and W 22 are given respectively by
W 11 = 1 2 δ ( μ 1 0 + g 1 ( s ) d s ) 2 + 0 + g 1 ( s ) d s + ( κ 2 ) 2 ϰ a 11 , a 12 , a 21 , a 22 ζ 1 4 0 + g 1 ( s ) d s + 3 ( κ 2 Ł 1 ) 2 ζ 1 δ ( κ 2 p 1 + 2 ) 2 p 1 + 2 ( ζ 1 ) p 1 + 1 ( 4 E u , v ( 0 ) μ 1 0 + g 1 ( s ) d s ) p 1 + ( κ 2 p 2 + 2 ) 2 p 2 + 2 ( ζ 1 ) p 2 + 1 ( 4 E u , v ( 0 ) μ 1 0 + g 1 ( s ) d s ) p 2 + ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 1 ( 4 E u , v ( 0 ) μ 1 0 + g 1 ( s ) d s ) ρ 1 ,
W 12 = 1 2 δ + ( κ 2 ) 2 ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 2 ( 4 E u , v ( 0 ) ) ρ 1 2 δ ( ρ 1 + 1 ) 2 ,
W 21 = 1 2 δ ( μ 2 0 + g 2 ( s ) d s ) 2 + 0 + g 2 ( s ) d s + ( κ 2 ) 2 ϰ a 11 , a 12 , a 21 , a 22 ζ 2 4 0 + g 2 ( s ) d s + 3 ( κ 2 Ł 2 ) 2 ζ 2 δ ( κ 2 p 1 + 2 ) 2 p 1 + 2 ( ζ 2 ) p 1 + 1 ( 4 E u , v ( 0 ) μ 2 0 + g 2 ( s ) d s ) p 1 + ( κ 2 p 2 + 2 ) 2 p 2 + 2 ( ζ 2 ) p 2 + 1 ( 4 E u , v ( 0 ) μ 2 0 + g 2 ( s ) d s ) p 2 + ( κ 2 ρ 2 + 2 ) 2 ρ 2 + 2 ( ζ 2 ) ρ 2 + 1 ( 4 E u , v ( 0 ) μ 2 0 + g 2 ( s ) d s ) ρ 2
and
W 22 = 1 2 δ + ( κ 2 ) 2 ( κ 2 ρ 2 + 2 ) 2 ρ 2 + 2 ( ζ 2 ) ρ 2 + 2 ( 4 E u , v ( 0 ) ) ρ 2 2 δ ( ρ 2 + 1 ) 2 .
Proof. 
As pointed in Lemma 6, in view of the assumptions I ( 0 ) > 0 and max ( β 1 , β 2 ) < 1 , we conclude by Theorem 3 that L 2 ( t ) exists globally in time. Differentiate both sides of the Equation (94), to yield
L 2 ( t ) = μ 1 t g 1 ( t s ) Ω u ( t ) A 1 ( u ( t ) u ( s ) ) d x d s Ω t g 1 ( t s ) u ( s ) d s A 1 t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x t g 1 ( t s ) Ω t u ( t ) A 1 ( u ( t ) u ( s ) ) d x d s t g 1 ( t s ) Ω | t u ( t ) | ρ 1 t u ( t ) ρ 1 + 1 ( u ( t ) u ( s ) ) d x d s t g 1 ( t s ) Ω t u ( t ) A 1 t u ( t ) d x d s 1 ρ 1 + 1 0 + g 1 ( s ) d s Ω | t u ( t ) | ρ 1 + 2 d x + μ 2 t g 2 ( t s ) Ω v ( t ) A 2 ( v ( t ) v ( s ) ) d x d s Ω t g 2 ( t s ) v ( s ) d s A 2 t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x t g 2 ( t s ) Ω t v ( t ) A 2 ( v ( t ) v ( s ) ) d x d s t g 2 ( t s ) Ω | t v ( t ) | ρ 2 t v ( t ) ρ 2 + 1 ( v ( t ) v ( s ) ) d x d s t g 2 ( t s ) Ω t v ( t ) A 2 t v ( t ) d x d s 1 ρ 2 + 1 0 + g 2 ( s ) d s Ω | t v ( t ) | ρ 2 + 2 d x Ω f 1 ( u ( t ) , v ( t ) ) t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x Ω f 2 ( u ( t ) , v ( t ) ) t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x + Ω ( a 11 t u ( t ) + a 12 t v ( t ) ) t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x + Ω ( a 21 t u ( t ) + a 22 t v ( t ) ) t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x , t R + .
We shall split (100) into several parts, and we shall treat each part separately. By some routine calculations, it is not difficult to find that
μ 1 t g 1 ( t s ) Ω u ( t ) A 1 ( u ( t ) u ( s ) ) d x d s Ω t g 1 ( t s ) u ( s ) d s A 1 t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x = Ω t g 1 ( t s ) ( u ( t ) u ( s ) ) d s A 1 t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x + ( μ 1 0 + g 1 ( s ) d s ) t g 1 ( t s ) Ω u ( t ) A 1 ( u ( t ) u ( s ) ) d x d s , t R + .
But by exploiting the positive definiteness of A 1 (see Assumption 1, especially (5), for the details) and by mainly utilizing the Cauchy–Schwarz inequality, we have
φ A 1 ( x ) ψ δ 2 ( μ 1 0 + g 1 ( s ) d s ) φ A 1 ( x ) φ + 1 2 δ ( μ 1 0 + g 1 ( s ) d s ) ψ A 1 ( x ) ψ , φ , ψ R N , x Ω .
This implies immediately
t g 1 ( t s ) Ω u ( t ) A 1 ( u ( t ) u ( s ) ) d x d s δ 0 + g 1 ( s ) d s 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + 1 2 δ ( μ 1 0 + g 1 ( s ) d s ) ( g 1 A 1 u ) ( t ) , t R + .
Besides, by using again the positive definiteness of A 1 , as in (102), and by employing Jensen’s inequality, we have
t g 1 ( t s ) ( U ( s ) ) d s A 1 t g 1 ( t s ) U ( s ) d s t g 1 ( t s ) d s t g 1 ( t s ) ( U ( s ) ) A 1 U ( s ) d s , U L 2 ( , t ; R N ) , t R + ,
which implies directly
Ω t g 1 ( t s ) ( u ( t ) u ( s ) ) d s A 1 t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x Ω t g 1 ( t s ) d s t g 1 ( t s ) ( u ( t ) u ( s ) ) A 1 ( u ( t ) u ( s ) ) d s d x = 0 + g 1 ( s ) d s ( g 1 A 1 u ) ( t ) , t R + .
Plug (102) and (103) into (101), to obtain
μ 1 t g 1 ( t s ) Ω u ( t ) A 1 ( u ( t ) u ( s ) ) d x d s Ω t g 1 ( t s ) u ( s ) d s A 1 t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x δ 2 0 + g 1 ( s ) d s Ω u ( t ) A 1 u ( t ) d x + M 11 7 ( g 1 A 1 u ) ( t ) , t R + ,
with the positive constant M 11 7 given by
M 11 7 = 1 2 δ ( μ 1 0 + g 1 ( s ) d s ) 2 + 0 + g 1 ( s ) d s .
And similarly, we have also
μ 2 t g 2 ( t s ) Ω v ( t ) A 2 ( v ( t ) v ( s ) ) d x d s Ω t g 2 ( t s ) v ( s ) d s A 2 t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x δ 2 0 + g 2 ( s ) d s Ω v ( t ) A 2 v ( t ) d x + M 21 7 ( g 2 A 2 v ) ( t ) , t R + ,
where the positive constant M 21 7 is given by
M 21 7 = 1 2 δ ( μ 2 0 + g 2 ( s ) d s ) 2 + 0 + g 2 ( s ) d s .
Thanks to the positive definiteness of A 1 (see Assumption 1, especially (5)), taking steps similar to those in the derivation of (102), we can prove
t g 1 ( t s ) Ω t u ( t ) A 1 ( u ( t ) u ( s ) ) d x d s δ g 1 ( 0 ) 2 Ω t u ( t ) A 1 t u ( t ) d x 1 2 δ ( g 1 A 1 u ) ( t ) , t R + ,
and
t g 2 ( t s ) Ω t v ( t ) A 2 ( v ( t ) v ( s ) ) d x d s δ g 2 ( 0 ) 2 Ω t v ( t ) A 2 t v ( t ) d x 1 2 δ ( g 2 A 2 v ) ( t ) , t R + .
Apply the Fenchel–Young inequality, use the notations introduced in (9) and (10), and conduct some routine calculations, to obtain
t g 1 ( t s ) Ω | t u ( t ) | ρ 1 t u ( t ) ρ 1 + 1 ( u ( t ) u ( s ) ) d x d s t g 1 ( t s ) Ω δ 2 ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 1 ( 4 E u , v ( 0 ) ) ρ 1 | t u ( t ) | 2 ρ 1 + 2 + ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 1 ( 4 E u , v ( 0 ) ) ρ 1 2 δ ( ρ 1 + 1 ) 2 | u ( t ) u ( s ) | 2 d x d s = δ g 1 ( 0 ) 2 ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 1 ( 4 E u , v ( 0 ) ) ρ 1 t u ( t ) L 2 ρ 1 + 2 ( Ω ) 2 ρ 1 + 2 ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 1 ( 4 E u , v ( 0 ) ) ρ 1 2 δ ( ρ 1 + 1 ) 2 t g 1 ( t s ) Ω | u ( t ) u ( s ) | 2 d x d s δ g 1 ( 0 ) 2 ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 1 ( 4 E u , v ( 0 ) ) ρ 1 t u ( t ) L 2 ρ 1 + 2 ( Ω ) 2 ρ 1 + 2 M 12 7 ( g 1 A 1 u ) ( t ) , t R + ,
with the positive constant M 12 7 given by
M 12 7 = ( κ 2 ) 2 ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 2 ( 4 E u , v ( 0 ) ) ρ 1 2 δ ( ρ 1 + 1 ) 2 .
To proceed further, we need the help of the following inequality
Ω t u ( t ) A 1 t u ( t ) d x + Ω t v ( t ) A 2 t v ( t ) d x 4 I ( t ) + ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + Ω t u ( t ) A 1 t u ( t ) d x + ( g 1 A 1 u ) ( t ) + ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + Ω t v ( t ) A 2 t v ( t ) d x + ( g 2 A 2 v ) ( t ) = 4 J ( t ) 4 E u , v ( 0 ) , t R + ,
which can be proved via applying Lemma 4 and using (56) as well as (57). With the help of (112), we can prove easily
t u ( t ) L 2 ρ 1 + 2 ( Ω ) 2 ρ 1 + 2 ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 1 ( Ω t u ( t ) A 1 t u ( t ) d x ) ρ 1 + 1 ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 1 ( 4 E u , v ( 0 ) ) ρ 1 Ω t u ( t ) A 1 t u ( t ) d x .
This, together with (110), implies
t g 1 ( t s ) Ω | t u ( t ) | ρ 1 t u ( t ) ρ 1 + 1 ( u ( t ) u ( s ) ) d x d s δ g 1 ( 0 ) 2 Ω t u ( t ) A 1 t u ( t ) d x M 12 7 ( g 1 A 1 u ) ( t ) , t R + .
By applying (112), we take steps similar to those used to obtain (113), to arrive at
t v ( t ) L 2 ρ 2 + 2 ( Ω ) 2 ρ 2 + 2 ( κ 2 ρ 2 + 2 ) 2 ρ 2 + 2 ( ζ 2 ) ρ 2 + 1 ( 4 E u , v ( 0 ) ) ρ 2 Ω t v ( t ) A 2 t v ( t ) d x .
With this at our hand, we can use the idea similar to the one utilized to establish the inequality (114), to prove successfully that
t g 2 ( t s ) Ω | t v ( t ) | ρ 2 t v ( t ) ρ 2 + 1 ( v ( t ) v ( s ) ) d x d s δ g 2 ( 0 ) 2 Ω t v ( t ) A 2 t v ( t ) d x M 22 7 ( g 2 A 2 v ) ( t ) , t R + ,
where the positive constant M 22 7 is given by
M 22 7 = ( κ 2 ) 2 ( κ 2 ρ 2 + 2 ) 2 ρ 2 + 2 ( ζ 2 ) ρ 2 + 2 ( 4 E u , v ( 0 ) ) ρ 2 2 δ ( ρ 2 + 1 ) 2 .
With the help of Assumption 5 and the notations in (9), (10) and (12), by mainly exploiting the Cauchy–Schwarz inequality and Jensen’s inequality, we have
Ω ( a 11 t u ( t ) + a 12 t v ( t ) ) t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x + Ω ( a 21 t u ( t ) + a 22 t v ( t ) ) t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x = Ω t u ( t ) t v ( t ) a 11 a 12 a 21 a 22 t g 1 ( t s ) ( u ( t ) u ( s ) ) d s t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x Ω t u ( t ) t v ( t ) a 11 a 12 a 21 a 22 t u ( t ) t v ( t ) d x + 1 4 Ω t g 1 ( t s ) ( u ( t ) u ( s ) ) d s t g 2 ( t s ) ( v ( t ) v ( s ) ) d s a 11 a 12 a 21 a 22 t g 1 ( t s ) ( u ( t ) u ( s ) ) d s t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x Ω t u ( t ) t v ( t ) a 11 a 12 a 21 a 22 t u ( t ) t v ( t ) d x + ϰ a 11 , a 12 , a 21 , a 22 4 Ω t g 1 ( t s ) ( u ( t ) u ( s ) ) d s 2 + t g 2 ( t s ) ( v ( t ) v ( s ) ) d s 2 d x Ω t u ( t ) t v ( t ) a 11 a 12 a 21 a 22 t u ( t ) t v ( t ) d x + ϰ a 11 , a 12 , a 21 , a 22 4 Ω t g 1 ( t s ) d s t g 1 ( t s ) | u ( t ) u ( s ) | 2 d s + t g 2 ( t s ) d s t g 2 ( t s ) | v ( t ) v ( s ) | 2 d s d x Ω t u ( t ) t v ( t ) a 11 a 12 a 21 a 22 t u ( t ) t v ( t ) d x + M 13 7 ( g 1 A 1 u ) ( t ) + M 23 7 ( g 2 A 2 v ) ( t ) , t R + ,
where the positive constants M 13 7 and M 23 7 are given respectively by
M 13 7 = ( κ 2 ) 2 ϰ a 11 , a 12 , a 21 , a 22 ζ 1 4 0 + g 1 ( s ) d s ,
and
M 23 7 = ( κ 2 ) 2 ϰ a 11 , a 12 , a 21 , a 22 ζ 2 4 0 + g 2 ( s ) d s .
With the growth condition (8) in Assumption 4 at our hand, we can prove via applying the Cauchy–Schwarz inequality and using the notations in (9) and (10) that
Ω f 1 ( u ( t ) , v ( t ) ) t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x Ω f 2 ( u ( t ) , v ( t ) ) t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x δ 1 2 t g 1 ( t s ) Ω | f 1 ( u ( t ) , v ( t ) ) | 2 d x d s + 1 2 δ 1 t g 1 ( t s ) Ω | u ( t ) u ( s ) | 2 d x d s + δ 2 2 t g 2 ( t s ) Ω | f 2 ( u ( t ) , v ( t ) ) | 2 d x d s + 1 2 δ 2 t g 2 ( t s ) Ω | v ( t ) v ( s ) | 2 d x d s 3 δ 1 ( Ł 1 ) 2 t g 1 ( t s ) u ( t ) L 2 p 1 + 2 2 p 1 + 2 + u ( t ) L 2 p 2 + 2 2 p 2 + 2 + u ( t ) L 2 ρ 1 + 2 2 ρ 1 + 2 + v ( t ) L 2 p 1 + 2 2 p 1 + 2 + v ( t ) L 2 p 2 + 2 2 p 2 + 2 + v ( t ) L 2 ρ 2 + 2 2 ρ 2 + 2 d s + ( κ 2 ) 2 ζ 1 2 δ 1 ( g 1 A 1 u ) ( t ) + 3 δ 2 ( Ł 2 ) 2 t g 2 ( t s ) u ( t ) L 2 p 1 + 2 2 p 1 + 2 + u ( t ) L 2 p 2 + 2 2 p 2 + 2 + u ( t ) L 2 ρ 1 + 2 2 ρ 1 + 2 + v ( t ) L 2 p 1 + 2 2 p 1 + 2 + v ( t ) L 2 p 2 + 2 2 p 2 + 2 + v ( t ) L 2 ρ 2 + 2 2 ρ 2 + 2 d s + ( κ 2 ) 2 ζ 2 2 δ 2 ( g 2 A 2 u ) ( t ) = 3 δ 1 ( Ł 1 ) 2 0 + g 1 ( s ) d s u ( t ) L 2 p 1 + 2 2 p 1 + 2 + u ( t ) L 2 p 2 + 2 2 p 2 + 2 + u ( t ) L 2 ρ 1 + 2 2 ρ 1 + 2 + v ( t ) L 2 p 1 + 2 2 p 1 + 2 + v ( t ) L 2 p 2 + 2 2 p 2 + 2 + v ( t ) L 2 ρ 2 + 2 2 ρ 2 + 2 + ( κ 2 ) 2 ζ 1 2 δ 1 ( g 1 A 1 u ) ( t ) + 3 δ 2 ( Ł 2 ) 2 0 + g 2 ( s ) d s u ( t ) L 2 p 1 + 2 2 p 1 + 2 + u ( t ) L 2 p 2 + 2 2 p 2 + 2 + u ( t ) L 2 ρ 1 + 2 2 ρ 1 + 2 + v ( t ) L 2 p 1 + 2 2 p 1 + 2 + v ( t ) L 2 p 2 + 2 2 p 2 + 2 + v ( t ) L 2 ρ 2 + 2 2 ρ 2 + 2 + ( κ 2 ) 2 ζ 2 2 δ 2 ( g 2 A 2 u ) ( t ) , t R + .
In view of the assumptions I ( 0 ) > 0 and max ( β 1 , β 2 ) < 1 , we conclude by Theorem 3 the inequality (72) holds true in R + . This, together with the notations in (9) and (10), implies
u ( t ) L 2 p 1 + 2 2 p 1 + 2 ( κ 2 p 1 + 2 ) 2 p 1 + 2 ( ζ 1 ) p 1 + 1 ( Ω u ( t ) A 1 u ( t ) d x ) p 1 + 1 ( κ 2 p 1 + 2 ) 2 p 1 + 2 ( ζ 1 ) p 1 + 1 ( 4 E u , v ( 0 ) μ 1 0 + g 1 ( s ) d s ) p 1 Ω u ( t ) A 1 u ( t ) d x ,
u ( t ) L 2 p 2 + 2 2 p 2 + 2 ( κ 2 p 2 + 2 ) 2 p 2 + 2 ( ζ 1 ) p 2 + 1 ( 4 E u , v ( 0 ) μ 1 0 + g 1 ( s ) d s ) p 2 Ω u ( t ) A 1 u ( t ) d x ,
u ( t ) L 2 ρ 1 + 2 2 ρ 1 + 2 ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 1 ( 4 E u , v ( 0 ) μ 1 0 + g 1 ( s ) d s ) ρ 1 Ω u ( t ) A 1 u ( t ) d x ,
v ( t ) L 2 p 1 + 2 2 p 1 + 2 ( κ 2 p 1 + 2 ) 2 p 1 + 2 ( ζ 2 ) p 1 + 1 ( 4 E u , v ( 0 ) μ 2 0 + g 2 ( s ) d s ) p 1 Ω v ( t ) A 2 v ( t ) d x ,
v ( t ) L 2 p 2 + 2 2 p 2 + 2 ( κ 2 p 2 + 2 ) 2 p 2 + 2 ( ζ 2 ) p 2 + 1 ( 4 E u , v ( 0 ) μ 2 0 + g 2 ( s ) d s ) p 2 Ω v ( t ) A 2 v ( t ) d x ,
and
v ( t ) L 2 ρ 2 + 2 2 ρ 2 + 2 ( κ 2 ρ 2 + 2 ) 2 ρ 2 + 2 ( ζ 2 ) ρ 2 + 1 ( 4 E u , v ( 0 ) μ 2 0 + g 2 ( s ) d s ) ρ 2 Ω v ( t ) A 2 v ( t ) d x .
Plug (121), (122), (123), (124), (125) and (126) into (120), to get
Ω f 1 ( u ( t ) , v ( t ) ) t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x Ω f 2 ( u ( t ) , v ( t ) ) t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x = 3 δ 1 ( Ł 1 ) 2 M 14 7 0 + g 1 ( s ) d s Ω u ( t ) A 1 u ( t ) d x + ( κ 2 ) 2 ζ 1 2 δ 1 ( g 1 A 1 u ) ( t ) + 3 δ 2 ( Ł 2 ) 2 M 24 7 0 + g 2 ( s ) d s Ω v ( t ) A 2 v ( t ) d x + ( κ 2 ) 2 ζ 2 2 δ 2 ( g 2 A 2 u ) ( t ) , t R + ,
where the positive constant δ i is given in an arbitrary way, the positive constants M 14 7 and M 24 7 are given respectively by
M 14 7 = ( κ 2 p 1 + 2 ) 2 p 1 + 2 ( ζ 1 ) p 1 + 1 ( 4 E u , v ( 0 ) μ 1 0 + g 1 ( s ) d s ) p 1 + ( κ 2 p 2 + 2 ) 2 p 2 + 2 ( ζ 1 ) p 2 + 1 ( 4 E u , v ( 0 ) μ 1 0 + g 1 ( s ) d s ) p 2 + ( κ 2 ρ 1 + 2 ) 2 ρ 1 + 2 ( ζ 1 ) ρ 1 + 1 ( 4 E u , v ( 0 ) μ 1 0 + g 1 ( s ) d s ) ρ 1 ,
and
M 24 7 = ( κ 2 p 1 + 2 ) 2 p 1 + 2 ( ζ 2 ) p 1 + 1 ( 4 E u , v ( 0 ) μ 2 0 + g 2 ( s ) d s ) p 1 + ( κ 2 p 2 + 2 ) 2 p 2 + 2 ( ζ 2 ) p 2 + 1 ( 4 E u , v ( 0 ) μ 2 0 + g 2 ( s ) d s ) p 2 + ( κ 2 ρ 2 + 2 ) 2 ρ 2 + 2 ( ζ 2 ) ρ 2 + 1 ( 4 E u , v ( 0 ) μ 2 0 + g 2 ( s ) d s ) ρ 2 .
Pick δ 1 and δ 2 so that
3 δ 1 ( Ł 1 ) 2 M 14 7 = 3 δ 2 ( Ł 2 ) 2 M 24 7 = δ 2 ,
to yield
Ω f 1 ( u ( t ) , v ( t ) ) t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x Ω f 2 ( u ( t ) , v ( t ) ) t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x δ 2 0 + g 1 ( s ) d s Ω u ( t ) A 1 u ( t ) d x + δ 2 0 + g 2 ( s ) d s Ω v ( t ) A 2 v ( t ) d x + 3 ( κ 2 Ł 1 ) 2 ζ 1 M 14 7 δ ( g 1 A 1 u ) ( t ) + 3 ( κ 2 Ł 2 ) 2 ζ 2 M 24 7 δ ( g 2 A 2 v ) ( t ) , t R + .
Plug (104), (106), (108), (109), (114), (115), (117) and (129) into (100) and perform some routine but tedious calculations, to finish the proof of Lemma 6. □
Remark 7.
In some occasions, we are inclined to establish an inequality similar to (117) without the term
Ω t u ( t ) t v ( t ) a 11 a 12 a 21 a 22 t u ( t ) t v ( t ) d x .
To this end, we actually have several different approaches. For example, by using the Cauchy–Schwarz inequalities (in several different forms), we have
Ω ( a 11 t u ( t ) + a 12 t v ( t ) ) t g 1 ( t s ) ( u ( t ) u ( s ) ) d s d x + Ω ( a 21 t u ( t ) + a 22 t v ( t ) ) t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x = Ω t u ( t ) t v ( t ) a 11 a 12 a 21 a 22 t g 1 ( t s ) ( u ( t ) u ( s ) ) d s t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x 1 2 ( κ 2 ) 2 Ω t u ( t ) t v ( t ) 0 + g 1 ( s ) d s δ g 1 ( 0 ) ζ 1 0 0 0 + g 2 ( s ) d s δ g 2 ( 0 ) ζ 2 t u ( t ) t v ( t ) d x + ( κ 2 ) 2 2 Ω t g 1 ( t s ) ( u ( t ) u ( s ) ) d s t g 2 ( t s ) ( v ( t ) v ( s ) ) d s a 11 a 12 a 21 a 22 · ζ 1 0 + g 1 ( s ) d s δ g 1 ( 0 ) 0 0 ζ 2 0 + g 2 ( s ) d s δ g 2 ( 0 ) a 11 a 12 a 21 a 22 t g 1 ( t s ) ( u ( t ) u ( s ) ) d s t g 2 ( t s ) ( v ( t ) v ( s ) ) d s d x 1 2 ( 0 + g 1 ( s ) d s δ g 1 ( 0 ) ) Ω t u ( t ) A 1 t u ( t ) d x + 1 2 ( 0 + g 2 ( s ) d s δ g 2 ( 0 ) ) Ω t v ( t ) A 2 t v ( t ) d x + ( κ 2 ) 2 ϰ a 11 , a 12 , a 21 , a 22 8 max ( ζ 1 0 + g 1 ( s ) d s δ g 1 ( 0 ) , ζ 2 0 + g 2 ( s ) d s δ g 2 ( 0 ) ) · Ω t g 1 ( t s ) ( u ( t ) u ( s ) ) d s 2 + t g 2 ( t s ) ( v ( t ) v ( s ) ) d s 2 d x 1 2 ( 0 + g 1 ( s ) d s δ g 1 ( 0 ) ) Ω t u ( t ) A 1 t u ( t ) d x + 1 2 ( 0 + g 2 ( s ) d s δ g 2 ( 0 ) ) Ω t v ( t ) A 2 t v ( t ) d x + ( κ 2 ) 2 ϰ a 11 , a 12 , a 21 , a 22 8 max ( ζ 1 0 + g 1 ( s ) d s δ g 1 ( 0 ) , ζ 2 0 + g 2 ( s ) d s δ g 2 ( 0 ) ) Ω t g 1 ( t s ) d s t g 1 ( t s ) | u ( t ) u ( s ) | 2 d s + t g 2 ( t s ) d s t g 2 ( t s ) | v ( t ) v ( s ) | 2 d s d x 1 2 ( 0 + g 1 ( s ) d s δ g 1 ( 0 ) ) Ω t u ( t ) A 1 t u ( t ) d x + 1 2 ( 0 + g 2 ( s ) d s δ g 2 ( 0 ) ) Ω t v ( t ) A 2 t v ( t ) d x + M ˜ 13 7 ( g 1 A 1 u ) ( t ) + M ˜ 23 7 ( g 2 A 2 v ) ( t ) , t R + ,
where the positive constants M ˜ 13 7 and M ˜ 23 7 , slightly different from M 13 7 and M 23 7 in (117) (see (118) and (119) for the details), are given by
M ˜ 13 7 = ( κ 2 ) 4 ϰ a 11 , a 12 , a 21 , a 22 ζ 1 8 max ( ζ 1 0 + g 1 ( s ) d s δ g 1 ( 0 ) , ζ 2 0 + g 2 ( s ) d s δ g 2 ( 0 ) ) 0 + g 1 ( s ) d s ,
and
M ˜ 23 7 = ( κ 2 ) 4 ϰ a 11 , a 12 , a 21 , a 22 ζ 2 8 max ( ζ 1 0 + g 1 ( s ) d s δ g 1 ( 0 ) , ζ 2 0 + g 2 ( s ) d s δ g 2 ( 0 ) ) 0 + g 2 ( s ) d s .
As mentioned above, compared to the estimate (117), the last three lines of the sequence (131) of inequalities do not include the term (130). As will be seen, this could be beneficial to relatively wide choice of portions of the energy functional E u , v ( t ) (see (18)), in the procedure of constructing modified energy functionals.
As with L 1 ( t ) and L 2 ( t ) (see (85) and (94), respectively), the following energy perturbation functional will be useful in our later presentation: To every pair ( u , v ) S R + (see (15) for the definition of S R + ), we associate the functional
L 3 ( t ) = μ 1 0 + g 1 ( s ) d s 0 + g 1 ( s ) d s 0 + g 1 ( s ) t s t Ω u ( r ) A 1 u ( r ) d x d r d s + μ 2 0 + g 2 ( s ) d s 0 + g 2 ( s ) d s 0 + g 2 ( s ) t s t Ω v ( r ) A 2 v ( r ) d x d r d s , t R + .
Lemma 7.
Suppose that Assumptions 1–5 hold true. If the associated functional I ( t ) given by (56) satisfies I ( 0 ) > 0 , the associated constants β 1 and β 2 , given by (58) and (59), respectively, satisfy max ( β 1 , β 2 ) < 1 , then weak solutions to IBVP (1) exist globally in time, and render the associated functional L 3 ( t ) given by (134) to satisfy
L 3 ( t ) 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x μ 1 0 + g 1 ( s ) d s 2 0 + g 1 ( s ) d s ( g 1 A 1 u ) ( t ) μ 2 0 + g 2 ( s ) d s 2 0 + g 2 ( s ) d s ( g 2 A 2 v ) ( t ) , t R + .
Proof. 
Differentiate directly L 3 ( t ) (see (134)) with respect to t, to arrive at
L 3 ( t ) = μ 1 0 + g 1 ( s ) d s 0 + g 1 ( s ) d s 0 + g 1 ( s ) d s Ω u ( t ) A 1 u ( t ) d x μ 1 0 + g 1 ( s ) d s 0 + g 1 ( s ) d s 0 + g 1 ( s ) Ω u ( t s ) A 1 u ( t s ) d x d s + μ 2 0 + g 2 ( s ) d s 0 + g 2 ( s ) d s 0 + g 2 ( s ) d s Ω v ( t ) A 2 v ( t ) d x μ 2 0 + g 2 ( s ) d s 0 + g 2 ( s ) d s 0 + g 2 ( s ) Ω v ( t s ) A 2 v ( t s ) d x d s = ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x μ 1 0 + g 1 ( s ) d s 0 + g 1 ( s ) d s t g 1 ( t s ) Ω u ( s ) A 1 u ( s ) d x d s + ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x μ 2 0 + g 2 ( s ) d s 0 + g 2 ( s ) d s t g 2 ( t s ) Ω v ( s ) A 2 v ( s ) d x d s , t R + .
By direct calculations, we have
t g 1 ( t s ) Ω u ( s ) A 1 u ( s ) d x d s = 1 2 t g 1 ( t s ) Ω ( u ( t ) u ( s ) ) A 1 ( u ( t ) u ( s ) ) d x d s + 1 2 t g 1 ( t s ) Ω ( u ( t ) + u ( s ) ) A 1 ( u ( t ) + u ( s ) ) d x d s t g 1 ( t s ) Ω u ( t ) A 1 u ( t ) d x d s = 1 2 t g 1 ( t s ) Ω ( u ( t ) + u ( s ) ) A 1 ( u ( t ) + u ( s ) ) d x d s + 1 2 ( g 1 A 1 u ) ( t ) 0 + g 1 ( s ) d s Ω u ( t ) A 1 u ( t ) d x , t R + .
And we have analogously
t g 2 ( t s ) Ω v ( s ) A 2 v ( s ) d x d s = 1 2 t g 2 ( t s ) Ω ( v ( t ) + v ( s ) ) A 2 ( v ( t ) + v ( s ) ) d x d s + 1 2 ( g 2 A 2 v ) ( t ) 0 + g 2 ( s ) d s Ω v ( t ) A 2 v ( t ) d x , t R + .
This, together with (136) and (137), implies immediately
L 3 ( t ) = 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x μ 1 0 + g 1 ( s ) d s 2 0 + g 1 ( s ) d s t g 1 ( t s ) Ω ( u ( t ) + u ( s ) ) A 1 ( u ( t ) + u ( s ) ) d x d s μ 2 0 + g 2 ( s ) d s 2 0 + g 2 ( s ) d s t g 2 ( t s ) Ω ( v ( t ) + v ( s ) ) A 2 ( v ( t ) + v ( s ) ) d x d s μ 1 0 + g 1 ( s ) d s 2 0 + g 1 ( s ) d s ( g 1 A 1 u ) ( t ) μ 2 0 + g 2 ( s ) d s 2 0 + g 2 ( s ) d s ( g 2 A 2 v ) ( t ) , t R + .
Thanks to Assumptions 1 and 2, it holds that
t g 1 ( t s ) Ω ( u ( t ) + u ( s ) ) A 1 ( u ( t ) + u ( s ) ) d x d s 0 , t g 2 ( t s ) Ω ( v ( t ) + v ( s ) ) A 2 ( v ( t ) + v ( s ) ) d x d s 0 , t R + ,
which, together with (138), implies that the proof of Lemma 7 is complete. □
To obtain our claimed general energy decay result, we need to design various modified (perturbed) energy functionals by adding energy perturbation functionals to the conventional energy functional E u , v ( t ) (see (18) for the precise definition). Here, we are in a position to introduce the modified energy functional
( t ) = m 0 E u , v ( t ) + m 1 L 1 ( t ) + m 2 L 2 ( t ) , t R + ,
where m k is a positive constant yet to be determined later, k = 0 , 1 , 2 .
Lemma 8.
Suppose that Assumptions 1–5 hold true. Then there exists a triple
( m 0 , m 1 , m 2 ) ( 0 , + ) 3 ,
such that each solution ( u , v ) to IBVP (1) makes the following differential inequality hold:
( t ) 1 ρ 1 + 1 Ω | t u ( t ) | ρ 1 + 2 d x 1 ρ 2 + 1 Ω | t v ( t ) | ρ 2 + 2 d x 3 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x 3 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x Ω t u ( t ) A 1 t u ( t ) d x Ω t v ( t ) A 2 t v ( t ) d x + μ 1 0 + g 1 ( s ) d s 3 0 + g 1 ( s ) d s ( g 1 A 1 u ) ( t ) + μ 2 0 + g 2 ( s ) d s 3 0 + g 2 ( s ) d s ( g 2 A 2 v ) ( t ) for a . e . t R + ,
whenever its initial datum ( u 0 , v 0 , u 1 , v 1 ) L ( R ; H 0 1 ( Ω ; R 2 ) ) × H 0 1 ( Ω ; R 2 ) render the associated functional I ( t ) given by (56) to satisfy I ( 0 ) > 0 , the associated constants β 1 and β 2 , given by (58) and (59), respectively, to satisfy max ( β 1 , β 2 ) < 1 , the associated constants 1 and 2 , given by (86) and (87), respectively, satisfy max ( 1 , 2 ) < 1 4 , and the constants W 11 and W 21 , given by (96) and (98), respectively, to satisfy
W 11 μ 1 0 + g 1 ( s ) d s 6 0 + g 1 ( s ) d s max ( 10 0 + g 1 ( s ) d s , 10 0 + g 2 ( s ) d s ) , W 21 μ 2 0 + g 2 ( s ) d s 6 0 + g 2 ( s ) d s max ( 10 0 + g 1 ( s ) d s , 10 0 + g 2 ( s ) d s ) .
Proof. 
Thanks to (139), it follows that
( t ) = m 0 d d t E u , v ( t ) + m 1 L 1 ( t ) + m 2 L 2 ( t ) for a . e . t R + ,
where ( m 0 , m 1 , m 2 ) ( 0 , + ) 3 , as in (139), is yet to be determined later. By Lemmas 3, 5 and 6, this, together with (55), (88) and (95), implies
( t ) 1 ρ 1 + 1 m 1 m 2 0 + g 1 ( s ) d s Ω | t u ( t ) | ρ 1 + 2 d x + 1 ρ 2 + 1 m 1 m 2 0 + g 2 ( s ) d s Ω | t v ( t ) | ρ 2 + 2 d x + δ m 2 0 + g 1 ( s ) d s m 1 2 ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x + δ m 2 0 + g 2 ( s ) d s m 1 2 ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + m 1 m 2 ( 0 + g 1 ( s ) d s δ g 1 ( 0 ) ) Ω t u ( t ) A 1 t u ( t ) d x + m 1 m 2 ( 0 + g 2 ( s ) d s δ g 2 ( 0 ) ) Ω t v ( t ) A 2 t v ( t ) d x + m 2 W 11 ( g 1 A 1 u ) ( t ) + m 1 μ 1 0 + g 1 ( s ) d s 0 + ( g 1 ( s ) ) 2 g 1 ( s ) G g 1 ( s ) d s ( ( g 1 G g 1 ) A 1 u ) ( t ) + m 2 W 21 ( g 2 A 2 v ) ( t ) + m 1 μ 2 0 + g 2 ( s ) d s 0 + ( g 2 ( s ) ) 2 g 2 ( s ) G g 2 ( s ) d s ( ( g 2 G g 2 ) A 2 v ) ( t ) + m 0 2 m 2 W 12 ( g 1 A 1 u ) ( t ) + m 0 2 m 2 W 22 ( g 2 A 2 v ) ( t ) + m 2 m 0 Ω t u ( t ) t v ( t ) a 11 a 12 a 21 a 22 t u ( t ) t v ( t ) d x for a . e . t R + .
Since g i : R + ( 0 , + ) is strictly decreasing, g i is of bounded variation and is therefore Lebesgue measurable, i = 1 , 2 . For any G > 0 , it is easy to show that
0 < ( g i ( s ) ) 2 g i ( s ) G g i ( s ) < g i ( s ) g i ( 0 ) , s R + , i = 1 , 2 .
It is also obvious to see that
lim G + ( g i ( s ) ) 2 g i ( s ) G g i ( s ) = 0 , s R + , i = 1 , 2 .
By Lebesgue’s dominated convergence theorem, we combine (143) and (144), to obtain
lim G + 0 + ( g i ( s ) ) 2 g i ( s ) G g i ( s ) d s = 0 , i = 1 , 2 .
By recalling the limit theory of one-variable functions and in view of (145), we conclude that there exists a G 0 > 0 such that
0 + ( g i ( s ) ) 2 g i ( s ) G g i ( s ) d s < ( μ i 0 + g i ( s ) d s ) 2 24 0 + g i ( s ) d s , G G 0 , i = 1 , 2 .
With (146) as one of the main tools, we have here by some routine computations
4 μ 1 0 + g 1 ( s ) d s 0 + ( g 1 ( s ) ) 2 g 1 ( s ) G 0 g 1 ( s ) d s ( ( g 1 G 0 g 1 ) A 1 u ) ( t ) = 4 μ 1 0 + g 1 ( s ) d s 0 + ( g 1 ( s ) ) 2 g 1 ( s ) G 0 g 1 ( s ) d s ( μ 1 0 + g 1 ( s ) d s ) 2 24 0 + g 1 ( s ) d s ( ( g 1 G 0 g 1 ) A 1 u ) ( t ) + μ 1 0 + g 1 ( s ) d s 6 0 + g 1 ( s ) d s ( g 1 A 1 u ) ( t ) G 0 ( μ 1 0 + g 1 ( s ) d s ) 6 0 + g 1 ( s ) d s ( g 1 A 1 u ) ( t ) μ 1 0 + g 1 ( s ) d s 6 0 + g 1 ( s ) d s ( g 1 A 1 u ) ( t ) G 0 ( μ 1 0 + g 1 ( s ) d s ) 6 0 + g 1 ( s ) d s ( g 1 A 1 u ) ( t ) for a . e . t R + .
With the aid of (146), we take several steps similar to (147), to arrive at
4 μ 2 0 + g 2 ( s ) d s 0 + ( g 2 ( s ) ) 2 g 2 ( s ) G 0 g 2 ( s ) d s ( ( g 2 G 0 g 2 ) A 2 v ) ( t ) μ 2 0 + g 2 ( s ) d s 6 0 + g 2 ( s ) d s ( g 2 A 2 v ) ( t ) G 0 ( μ 2 0 + g 2 ( s ) d s ) 6 0 + g 2 ( s ) d s ( g 2 A 2 v ) ( t ) for a . e . t R + .
Now we are ready to choose appropriate values of the parameters δ and m i , i = 0 , 1 , 2 . Actually, there are many ways to accomplish this goal. For instance, we can pick
m 1 = 4 , m 2 = max ( 10 0 + g 1 ( s ) d s , 10 0 + g 2 ( s ) d s ) ,
δ = min ( 0 + g 1 ( s ) d s 2 g 1 ( 0 ) , μ 1 0 + g 1 ( s ) d s 2 0 + g 1 ( s ) d s max ( 10 0 + g 1 ( s ) d s , 10 0 + g 2 ( s ) d s ) , 0 + g 2 ( s ) d s 2 g 2 ( 0 ) , μ 2 0 + g 2 ( s ) d s 2 0 + g 2 ( s ) d s max ( 10 0 + g 1 ( s ) d s , 10 0 + g 2 ( s ) d s ) ) ,
and
m 0 = max ( max ( 10 0 + g 1 ( s ) d s , 10 0 + g 2 ( s ) d s ) , 2 sup { max ( W 12 max ( 10 0 + g 1 ( s ) d s , 10 0 + g 2 ( s ) d s ) + G 0 ( μ 2 0 + g 2 ( s ) d s ) 6 0 + g 2 ( s ) d s , W 22 max ( 10 0 + g 1 ( s ) d s , 10 0 + g 2 ( s ) d s ) + G 0 ( μ 2 0 + g 2 ( s ) d s ) 6 0 + g 2 ( s ) d s ) ; max ( β 1 , β 2 ) < 1 , max ( 1 , 2 ) < 1 4 and ( 141 ) is satisfied } ) .
It is not difficult to find that when the parameters δ and m i , i = 0 , 1 , 2 , take the values shown above, it holds that
m 1 m 2 0 + g 1 ( s ) d s 1 , m 1 m 2 0 + g 2 ( s ) d s 1 , δ m 2 0 + g 1 ( s ) d s m 1 2 ( μ 1 0 + g 1 ( s ) d s ) 3 2 ( μ 1 0 + g 1 ( s ) d s ) , δ m 2 0 + g 2 ( s ) d s m 1 2 ( μ 2 0 + g 2 ( s ) d s ) 3 2 ( μ 2 0 + g 2 ( s ) d s ) , m 1 m 2 ( 0 + g 1 ( s ) d s δ g 1 ( 0 ) ) 4 10 0 + g 1 ( s ) d s ( 0 + g 1 ( s ) d s 0 + g 1 ( s ) d s 2 g 1 ( 0 ) · g 1 ( 0 ) ) 1 , m 1 m 2 ( 0 + g 2 ( s ) d s δ g 2 ( 0 ) ) 4 10 0 + g 2 ( s ) d s ( 0 + g 2 ( s ) d s 0 + g 2 ( s ) d s 2 g 2 ( 0 ) · g 2 ( 0 ) ) 1 , m 0 2 m 2 W 12 G 0 ( μ 1 0 + g 1 ( s ) d s ) 6 0 + g 1 ( s ) d s 0 , m 0 2 m 2 W 22 G 0 ( μ 2 0 + g 2 ( s ) d s ) 6 0 + g 2 ( s ) d s 0 , m 2 m 0 0 , m 2 W 11 μ 1 0 + g 1 ( s ) d s 6 0 + g 1 ( s ) d s , m 2 W 21 μ 2 0 + g 2 ( s ) d s 6 0 + g 2 ( s ) d s .
In view of the fact that the term (130) in (142) is non-negative, by some routine calculations, we conclude immediately: This, together with (142), (147) and (148), implies that the proof of Lemma 8 is complete. □
Theorem 4.
Suppose that Assumptions 1–5 hold true. For every quadruple
( u 0 , v 0 , u 1 , v 1 ) L ( R ; H 0 1 ( Ω ; R 2 ) ) × H 0 1 ( Ω ; R 2 )
of initial datum, if it renders the associated functional I ( t ) given by (56) to satisfy I ( 0 ) > 0 , the associated constants β 1 and β 2 , given respectively by (58) and (59), to satisfy max ( β 1 , β 2 ) < 1 , the associated constants 1 and 2 , given respectively by (86) and (87), to satisfy max ( 1 , 2 ) < 1 4 , and the associated constants W 11 and W 21 , given respectively by (96) and (98), to satisfy (141), then the corresponding global in time solution ( u , v ) S R + (see (15) for the definition of S R + ) to IBVP (1) makes the associated energy E u , v ( t ) , given by (18), satisfy
E u , v ( t ) ϖ 1 K 1 ( ϖ 2 0 t min ( ξ 1 ( s ) , ξ 2 ( s ) ) d s ) ) for all t R + ,
where the positive constants ϖ 1 and ϖ 2 , depending on ( u 0 , v 0 , u 1 , v 1 ) , is independent of the time variable t, and K 1 is the inverse function of K ( t ) which is given by
K ( t ) = t min ( r 1 , r 2 ) d s s K ( s ) for all t ( 0 , min ( r 1 , r 2 ) ] .
Remark 8.
By Assumption 2 (in particular, the restriction (6)), we conclude that the K ( t ) , given by (153), is well-defined as a continuous function, that it takes non-negative real numbers as its values, and that it obeys
K ( t ) = 1 s K ( s ) for all t ( 0 , min ( r 1 , r 2 ) ] .
This implies, in particular, that the function K ( t ) is strictly monotonically decreasing in the interval ( 0 , min ( r 1 , r 2 ) ] . We have directly by applying the definition (153) of K ( t )
lim t min ( r 1 , r 2 ) K ( t ) = K ( min ( r 1 , r 2 ) ) = 0 .
We apply Assumption 2 again and the definition (153) of K ( t ) , to get lim t 0 + K ( t ) = + . In conclusion, K ( t ) , given by (153), is a strictly decreasing continuous function mapping the interval ( 0 , min ( r 1 , r 2 ) ] onto the closed interval R + . This implies, among other things, the inverse function K 1 of K , appearing in (152), is well-defined, and K 1 actually maps the closed interval R + onto the interval ( 0 , min ( r 1 , r 2 ) ] .
Proof of Theorem 4. 
In view of (18), (82) and (140), we have
( t ) ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) E u , v ( t ) ( 1 ρ 1 + 1 ρ 2 + 2 ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) ) Ω | t u ( t ) | ρ 1 + 2 d x ( 1 ρ 2 + 1 ρ 1 + 2 ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) ) Ω | t v ( t ) | ρ 2 + 2 d x 1 2 ( 3 ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) ) ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x 1 2 ( 3 ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) ) ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x 1 2 ( 2 ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) ) Ω t u ( t ) A 1 t u ( t ) d x 1 2 ( 2 ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) ) Ω t v ( t ) A 2 t v ( t ) d x ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) Ω F ( u ( t ) , v ( t ) ) d x + M 11 8 ( g 1 A 1 u ) ( t ) + M 12 8 ( g 2 A 2 v ) ( t ) ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) E u , v ( t ) 1 4 ( 6 3 ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) ) ( μ 1 0 + g 1 ( s ) d s ) Ω u ( t ) A 1 u ( t ) d x 1 4 ( 6 3 ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) ) ( μ 2 0 + g 2 ( s ) d s ) Ω v ( t ) A 2 v ( t ) d x + M 11 8 ( g 1 A 1 u ) ( t ) + M 12 8 ( g 2 A 2 v ) ( t ) M 2 8 E u , v ( t ) + M 11 8 ( g 1 A 1 u ) ( t ) + M 12 8 ( g 2 A 2 v ) ( t ) for a . e . t R + ,
in which the constants M 11 8 , M 12 8 and M 2 8 are given respectively by
M 11 8 = μ 1 0 + g 1 ( s ) d s 3 0 + g 1 ( s ) d s + ( ρ 1 + 2 ) ( ρ 2 + 2 ) 2 ( ρ 1 + 1 ) ( ρ 2 + 2 ) + 2 ( ρ 1 + 2 ) ( ρ 2 + 1 ) ,
M 12 8 = μ 2 0 + g 2 ( s ) d s 3 0 + g 2 ( s ) d s + ( ρ 1 + 2 ) ( ρ 2 + 2 ) 2 ( ρ 1 + 1 ) ( ρ 2 + 2 ) + 2 ( ρ 1 + 2 ) ( ρ 2 + 1 ) ,
and
M 2 8 = ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) .
Case I. K 1 and K 2 are both linear.
By some routine but tedious calculations, we have
( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) M 2 8 E u , v ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) + M 11 8 min ( ξ 1 ( t ) , ξ 2 ( t ) ) ( g 1 A 1 u ) ( t ) + M 12 8 min ( ξ 1 ( t ) , ξ 2 ( t ) ) ( g 2 A 2 v ) ( t ) , M 2 8 E u , v ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) + M 11 8 ξ 1 ( t ) ( g 1 A 1 u ) ( t ) + M 12 8 ξ 2 ( t ) ( g 2 A 2 v ) ( t ) , M 2 8 E u , v ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) M 31 8 ( g 1 A 1 u ) ( t ) M 32 8 ( g 2 A 2 v ) ( t ) , M 2 8 E u , v ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) M 4 8 d d t E u , v ( t ) for a . e . t R + ,
or equivalently, we have
d d t ( M 4 8 E u , v ( t ) + ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) ) = M 4 8 d d t E u , v ( t ) + ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) + ( t ) d d t min ( ξ 1 ( t ) , ξ 2 ( t ) ) M 2 8 E u , v ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) for a . e . t R + .
M 4 8 E u , v ( t ) + ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) E u , v ( t ) , more precisely, there exist two positive constants M 5 8 and M 6 8 , such that
1 M 6 8 E u , v ( t ) M 4 8 E u , v ( t ) + ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) 1 M 5 8 E u , v ( t ) , t R + .
Therefore, we deduce from (158) immediately that
d d t ( M 4 8 E u , v ( t ) + ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) ) M 2 8 M 5 8 min ( ξ 1 ( t ) , ξ 2 ( t ) ) ( M 4 8 E u , v ( t ) + ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) ) .
This, together with Gronwall’s Lemma, implies
( M 4 8 E u , v ( t ) + ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) ) ( M 4 8 E u , v ( 0 ) + ( 0 ) min ( ξ 1 ( 0 ) , ξ 2 ( 0 ) ) ) e M 2 8 M 5 8 0 t min ( ξ 1 ( s ) , ξ 2 ( s ) ) d s , t R + .
By recalling the afore-mentioned equivalence (159), we realize that this implies
E u , v ( t ) M 7 8 e M 2 8 M 5 8 0 t min ( ξ 1 ( s ) , ξ 2 ( s ) ) d s , t R + ,
where the positive constant M 7 8 is given by
M 7 8 = M 6 8 ( M 4 8 E u , v ( 0 ) + ( 0 ) min ( ξ 1 ( 0 ) , ξ 2 ( 0 ) ) ) .
Case II. K 1 or K 2 is non-linear.
We combine ( t ) (with δ and m i , i = 0 , 1 , 2 given as in (149), (150) and (151), respectively) and L 3 ( t ) , to associate with each solution ( u , v ) to IBVP (1) a new functional
𝟊 ( t ) = ( t ) + 17 24 L 3 ( t ) , t R + .
Obviously, 𝟊 ( t ) 0 for all t R + . Besides, by Lemmas 7 and 8, there exists a positive constant M 1 9 such that
𝟊 ( t ) E u , v ( t ) M 1 9 for a . e . t R + .
For example, we could put M 1 9 = 1 M 2 9 with
M 2 9 = min ( 1 9 , ( ρ 1 + 2 ) ( ρ 2 + 2 ) ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) , ( μ 1 0 + g 1 ( s ) d s ) ( μ 2 0 + g 2 ( s ) d s ) 24 0 + g 2 ( s ) d s ( μ 1 0 + g 1 ( s ) d s ) + 24 0 + g 1 ( s ) d s ( μ 2 0 + g 2 ( s ) d s ) ) .
In this situation, we have
𝟊 ( t ) = ( t ) + 17 24 L 3 ( t ) E u , v ( t ) M 1 9 ( 1 ρ 1 + 1 ρ 2 + 2 ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) ) t u ( t ) L ρ 1 + 2 ( Ω ) ρ 1 + 2 ( 1 ρ 2 + 1 ρ 1 + 2 ( ρ 1 + 1 ) ( ρ 2 + 2 ) + ( ρ 1 + 2 ) ( ρ 2 + 1 ) ) t v ( t ) L ρ 2 + 2 ( Ω ) ρ 2 + 2 17 18 Ω t u ( t ) A 1 t u ( t ) d x 17 18 Ω t v ( t ) A 2 t v ( t ) d x 1 48 ( g 1 A 1 u ) ( t ) ( μ 1 0 + g 1 ( s ) d s 0 + g 1 ( s ) d s ( μ 1 0 + g 1 ( s ) d s ) ( μ 2 0 + g 2 ( s ) d s ) 0 + g 2 ( s ) d s ( μ 1 0 + g 1 ( s ) d s ) + 0 + g 1 ( s ) d s ( μ 2 0 + g 2 ( s ) d s ) ) 1 48 ( g 2 A 2 v ) ( t ) ( μ 2 0 + g 2 ( s ) d s 0 + g 2 ( s ) d s ( μ 1 0 + g 1 ( s ) d s ) ( μ 2 0 + g 2 ( s ) d s ) 0 + g 2 ( s ) d s ( μ 1 0 + g 1 ( s ) d s ) + 0 + g 1 ( s ) d s ( μ 2 0 + g 2 ( s ) d s ) ) E u , v ( t ) M 1 9 for a . e . t R + ,
which implies directly
0 t E u , v ( s ) d s M 1 9 0 t 𝟊 ( s ) d s = M 1 9 ( 𝟊 ( 0 ) 𝟊 ( t ) ) M 1 9 𝟊 ( 0 ) for all t R + .
This, in turn, implies
0 + E u , v ( t ) d t = lim t + 0 t E u , v ( s ) d s lim   sup t + 0 t E u , v ( s ) d s M 1 9 𝟊 ( 0 ) .
Fix provisionally t (sufficiently large if necessary) in the interval ( 0 , + ) , write
U ( t ) = t Ω ( u ( t ) u ( s ) ) A 1 ( u ( t ) u ( s ) ) d x d s , V ( t ) = t Ω ( v ( t ) v ( s ) ) A 2 ( v ( t ) v ( s ) ) d x d s ,
and pick a q ( 0 , + ) such that
U ( t ) < min ( 1 , r 1 , r 2 ) q g 1 ( 0 ) , V ( t ) < min ( 1 , r 1 , r 2 ) q g 2 ( 0 ) , t ( 0 , + ) .
We assume first that U ( t ) 0 and V ( t ) 0 . In this case, we deduce from (154) that
( t ) M 2 8 E u , v ( t ) + M 11 8 ( g 1 A 1 u ) ( t ) + M 12 8 ( g 2 A 2 v ) ( t ) M 2 8 E u , v ( t ) + g 1 ( 0 ) M 11 8 U ( t ) + g 2 ( 0 ) M 12 8 V ( t ) = M 2 8 E u , v ( t ) for a . e . t R + .
This, together with the equivalence ( t ) E u , v ( t ) , implies that the energy E u , v ( t ) associated to the system (1) the decays exponentially as time t escapes to infinity. Let us now assume that U ( t ) > 0 and V ( t ) > 0 when t is sufficiently large. Without loss of generality, we assume that U ( t ) > 0 and V ( t ) > 0 for all t R + . With Jensen’s inequality as one of our main tools, by some routine but tedious computations, we arrive at
K 1 q ( g 1 A 1 u ) ( t ) = K 1 t ( q A ( t ) g 1 ( t s ) ) Ω ( u ( t ) u ( s ) ) A 1 ( u ( t ) u ( s ) ) d x d s t Ω ( u ( t ) u ( s ) ) A 1 ( u ( t ) u ( s ) ) d x d s t K 1 ( q A ( t ) g 1 ( t s ) ) Ω ( u ( t ) u ( s ) ) A 1 ( u ( t ) u ( s ) ) d x d s t Ω ( u ( t ) u ( s ) ) A 1 ( u ( t ) u ( s ) ) d x d s q t ξ 1 ( t s ) K 1 g 1 ( t s ) Ω ( u ( t ) u ( s ) ) A 1 ( u ( t ) u ( s ) ) d x d s ξ 1 ( t ) q t g 1 ( t s ) Ω ( u ( t ) u ( s ) ) A 1 ( u ( t ) u ( s ) ) d x d s ξ 1 ( t ) = q ( g 1 A 1 u ) ( t ) ξ 1 ( t ) for a . e . t R + .
Since K 1 is strictly increasing (see Assumption 2 for the details), it follows immediately
( g 1 A 1 u ) ( t ) 1 q K 1 1 q ( g 1 A 1 u ) ( t ) ξ 1 ( t ) for a . e . t R + .
And analogously, we can prove
( g 2 A 2 v ) ( t ) 1 q K 2 1 q ( g 2 A 2 v ) ( t ) ξ 2 ( t ) for a . e . t R + .
Combine (154), (166) and (167), to obtain
( t ) M 2 8 E u , v ( t ) + M 11 8 ( g 1 A 1 u ) ( t ) + M 12 8 ( g 2 A 2 v ) ( t ) M 2 8 E u , v ( t ) + M 11 8 q K 1 1 q ( g 1 A 1 u ) ( t ) ξ 1 ( t ) + M 12 8 q K 2 1 q ( g 2 A 2 v ) ( t ) ξ 2 ( t ) for a . e . t R + .
Let us define a new function K ( t ) by giving
K ( t ) = max ( K 1 ( t ) , K 2 ( t ) ) for all t ( 0 , min ( r 1 , r 2 ) ] .
It is not difficult to check that K ( t ) is strictly increasing and strictly convex. Put
( t ) = K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) ( t ) + E u , v ( t ) for a . e . t R + ,
where the positive constant ε 0 is suitably chosen so that
ε 0 E u , v ( t ) E u , v ( 0 ) min ( r 1 , r 2 ) for all t R + .
By direct computations, we deduce from (168) and (169) that
( t ) = ε 0 ( t ) E u , v ( 0 ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) d d t E u , v ( t ) + K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) ( t ) + d d t E u , v ( t ) d d t E u , v ( t ) M 2 8 E u , v ( t ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) + M 11 8 q K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) K 1 1 q ( g 1 A 1 u ) ( t ) ξ 1 ( t ) + M 12 8 q K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) K 2 1 q ( g 2 A 2 v ) ( t ) ξ 2 ( t ) for a . e . t R + .
By applying Assumption 2, we have
lim δ 0 + K i 1 ( δ ) K 1 ( δ ) = k i , i = 1 , 2 .
By applying the method of change variable, we have furthermore
lim δ 0 + K i 1 ( K ( δ ) ) δ = k i , i = 1 , 2 .
By recalling the limit theory of one-variable functions, we conclude that there exists a positive constant min ( r 1 , r 2 ) , such that
K i 1 ( K ( δ ) ) < 2 k i δ , i = 1 , 2 .
whenever the positive variable δ does not exceed . With (171) as one of our main tools, by applying mainly the Fenchel–Young inequality, we have
K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) K 1 1 q ( g 1 A 1 u ) ( t ) ξ 1 ( t ) K 1 * ( K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) ) + K 1 ( K 1 1 q ( g 1 A 1 u ) ( t ) ξ 1 ( t ) ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) K 1 1 ( K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) ) K 1 ( K 1 1 ( K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) ) ) q ( g 1 A 1 u ) ( t ) ξ 1 ( t ) 2 k 1 ε 0 E u , v ( t ) E u , v ( 0 ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) q ( g 1 A 1 u ) ( t ) ξ 1 ( t ) for a . e . t R + .
By analogy to (172), we have
K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) K 2 1 q ( g 2 A 2 v ) ( t ) ξ 2 ( t ) 2 k 2 ε 0 E u , v ( t ) E u , v ( 0 ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) q ( g 2 A 2 v ) ( t ) ξ 2 ( t ) for a . e . t R + .
We substitute (172) and (173) into (170), to conclude
( t ) d d t E u , v ( t ) M 2 8 E u , v ( t ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) + 2 k 1 ε 0 M 11 8 q E u , v ( t ) E u , v ( 0 ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) M 11 8 ( g 1 A 1 u ) ( t ) ξ 1 ( t ) + 2 k 2 ε 0 M 12 8 q E u , v ( t ) E u , v ( 0 ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) M 12 8 ( g 2 A 2 v ) ( t ) ξ 2 ( t ) 2 ε 0 k 1 M 11 8 q E u , v ( t ) E u , v ( 0 ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) M 11 8 ( g 1 A 1 u ) ( t ) ξ 1 ( t ) + 2 ε 0 k 2 M 12 8 q E u , v ( t ) E u , v ( 0 ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) M 12 8 ( g 2 A 2 v ) ( t ) ξ 2 ( t ) M 2 8 E u , v ( t ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) for a . e . t R + .
With Lemma 3 (in particular, (55)) and (174) as our main tools, we perform some routine computations, to arrive at
( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) 2 ε 0 k 1 M 11 8 q E u , v ( t ) E u , v ( 0 ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) M 11 8 ( g 1 A 1 u ) ( t ) ξ 1 ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) + 2 ε 0 k 2 M 12 8 q E u , v ( t ) E u , v ( 0 ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) M 12 8 ( g 2 A 2 v ) ( t ) ξ 2 ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) M 2 8 E u , v ( t ) K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) M 1 10 K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) M 2 10 d d t E u , v ( t ) for a . e . t R + ,
in which, M 1 10 , M 2 10 and K ( t ) are given respectively by
M 1 10 = M 2 8 E u , v ( 0 ) ε 0 2 q ( k 1 M 11 8 + k 2 M 12 8 ) ,
M 2 10 = 2 max ( M 11 8 , M 12 8 ) ,
and
K ( t ) = t K ( t ) for all t ( 0 , min ( r 1 , r 2 ) ] .
The strict increasing monotonicity of K ( t ) , the strict convexity of K ( t ) and the continuity of K ( t ) imply that the function K ( t ) is strictly increasing. For the sake our later presentation, we introduce the auxiliary functional
˜ ( t ) = ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) + M 2 10 E u , v ( t ) for all t R + .
Since M 2 10 > 0 , ˜ ( t ) E u , v ( t ) . More exactly, there exists, as with (159), positive constants M 1 11 and M 2 11 , such that
M 1 11 ˜ ( t ) E u , v ( t ) M 2 11 ˜ ( t ) for all t R + .
With the help of (175) and the definition (177) of ˜ ( t ) , we have
˜ ( t ) = ( t ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) + ( t ) d d t min ( ξ 1 ( t ) , ξ 2 ( t ) ) + M 2 10 d d t E u , v ( t ) M 1 10 K ( ε 0 E u , v ( t ) E u , v ( 0 ) ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) for a . e . t R + ,
which, together with the equivalence (178), implies
˜ ( t ) M 1 10 K ( ε 0 M 1 11 ˜ ( t ) E u , v ( 0 ) ) min ( ξ 1 ( t ) , ξ 2 ( t ) ) for a . e . t R + ,
where the positive constant ε 0 is appropriately picked so that
ε 0 M 1 11 ˜ ( t ) E u , v ( 0 ) min ( r 1 , r 2 ) for all t R + .
Based on the differential inequality (179), we perform some routine calculations, to yield
M 1 10 0 t min ( ξ 1 ( s ) , ξ 2 ( s ) ) d s 0 t ˜ ( s ) K ( ε 0 M 1 11 ˜ ( s ) E u , v ( 0 ) ) d s = E u , v ( 0 ) ε 0 M 1 11 ε 0 M 1 11 ˜ ( t ) E u , v ( 0 ) ε 0 M 1 11 ˜ ( 0 ) E u , v ( 0 ) d s K ( s ) E u , v ( 0 ) ε 0 M 1 11 ε 0 M 1 11 ˜ ( t ) E u , v ( 0 ) min ( r 1 , r 2 ) d s K ( s ) = E u , v ( 0 ) ε 0 M 1 11 K ( ε 0 M 1 11 ˜ ( t ) E u , v ( 0 ) ) for all t R + ,
where the function K ( t ) is defined as in (153). It is not difficult to observe that
K ( t ) = t min ( r 1 , r 2 ) d s K ( s ) for all t ( 0 , min ( r 1 , r 2 ) ] ,
where the function K ( t ) is given as in (176). By recalling that the function K ( t ) is strictly decreasing (see Remark 8 for the detailed explanation), we deduce from (180) that
ε 0 M 1 11 ˜ ( t ) E u , v ( 0 ) K 1 ( ε 0 M 1 10 M 1 11 E u , v ( 0 ) 0 t min ( ξ 1 ( s ) , ξ 2 ( s ) ) d s ) ) for all t R + ,
which, together with the equivalence (178), implies directly
E u , v ( t ) M 2 11 E u , v ( 0 ) ε 0 M 1 11 K 1 ( ε 0 M 1 10 M 1 11 E u , v ( 0 ) 0 t min ( ξ 1 ( s ) , ξ 2 ( s ) ) d s ) ) for all t R + .
This, together with (160), implies that the proof of Theorem 4 is complete. □

4. Conclusions

In this paper, we studied the initial boundary value problem (that is, IBVP (1)) for a coupled system of two quasi-linear viscoelastic space-variable coefficient wave equations. We proved (see Theorems 1 and 2 for the details), under some seemingly natural conditions on A 1 , A 2 , μ 1 , μ 2 , ρ 1 , ρ 2 , f 1 , f 2 , g 1 , g 2 , a 11 , a 12 , a 21 and a 22 , via the celebrated Faedo–Galerkin method, that IBVP (1) admits a local solution in the sense of Definitions 1 and 2 for every initial datum in the space
L ( R ; H 0 1 ( Ω ; R 2 ) ) × H 0 1 ( Ω ; R 2 ) .
Based on our new obtained local existence results, we proved, via establishing a priori inequalities, a global existence result for solutions, having small initial data, to IBVP (1) (see Theorem 3 for the details). Based on our new established global existence result, we proved via constructing various modified energy functionals (functionals, equivalent to the energy functional E u , v ( t ) , defined by (18), of IBVP (1) and can be seen as Lyapunov functional from other perspectives), that if the initial data satisfy some additional conditions, then global in time solutions would decrease to zero, at the optimal decaying rate in a sense given by Remark 2.3 in Reference [22], as time escapes to infinity; see Theorem 4.

Author Contributions

Conceptualization, C.W. (Chengqiang Wang), X.Z. and Z.L.; methodology, C.W. (Chengqiang Wang), X.Z. and Z.L.; validation, C.W. (Chengqiang Wang) and X.Z.; formal analysis, C.W. (Chengqiang Wang); investigation, C.W. (Chengqiang Wang), C.W. (Can Wang), X.Z. and Z.L.; resources, C.W. (Chengqiang Wang), C.W. (Can Wang), X.Z. and Z.L.; writing—original draft preparation, C.W. (Chengqiang Wang), C.W. (Can Wang); writing—review and editing, C.W. (Chengqiang Wang), X.Z. and Z.L.; supervision, C.W. (Can Wang), X.Z. and Z.L.; funding acquisition, C.W. (Chengqiang Wang). All authors have read and agreed to the published version of the manuscript.

Funding

Chengqiang Wang is supported partially by KJZD Programme(#CS19ZA10) of Chengdu Normal University, by Startup Foundation for Newly Recruited Employees and Xichu Talents Foundation of Suqian University (#2022XRC033), NSFC (#11701050), and Jiangsu Qin–Lan Project of Fostering Excellent Teaching Team ‘University Mathematics Teaching Team’.

Data Availability Statement

Our results are theoretical, therefore we have no data to offer.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Evans, L. Partial Differential Equations, 2nd ed.; American Mathematical Society: Providence, RI, USA, 2010. [Google Scholar]
  2. Ali, I.; Saleem, M.T. Spatiotemporal dynamics of reaction–diffusion system and its application to turing pattern formation in a gray–scott model. Mathematics 2023, 11, 1459. [Google Scholar] [CrossRef]
  3. Creus, G.J. Viscoelasticity–Basic Theory and Applications to Concrete Structures; Springer: New York, NY, USA, 1986. [Google Scholar]
  4. Muñoz Rivera, J.E. Asymptotic behaviour in linear viscoelasticity. Quart. Appl. Math. 1994, 52, 628–648. [Google Scholar] [CrossRef]
  5. Muñoz Rivera, J.E.; Lapa, E.C.; Barreto, R. Decay rates for viscoelastic plates with memory. J. Elast. 1996, 44, 61–87. [Google Scholar] [CrossRef]
  6. Durdiev, D.K.; Totieva, Z.D. Kernel Determination Problems in Hyperbolic Integro-Differential Equations; Springer: Singapore, 2023. [Google Scholar]
  7. Aassila, M.; Cavalcanti, M.M.; Soriano, J.A. Asymptotic stability and energy decay rates for solutions of the wave equation with memory in a starsd shaped domain. SIAM J. Control Optim. 2000, 38, 1581–1602. [Google Scholar] [CrossRef]
  8. Cavalcanti, M.M.; Domingos Cavalcanti, V.N.; Ferreira, J. Existence and uniform decay for nonlinear viscoelastic equation with strong damping. Math. Methods Appl. Sci. 2001, 24, 1043–1053. [Google Scholar] [CrossRef]
  9. Cavalcanti, M.M.; Portillo Oquendo, H. Frictional versus viscoelastic damping in a semilinear wave equation. SIAM J. Control Optim. 2003, 42, 1310–1324. [Google Scholar] [CrossRef]
  10. Berrimi, S.; Messaoudi, S.A. Existence and decay of solutions of a viscoelastic equation with a nonlinear source. Nonlinear Anal. Theory Methods Appl. 2006, 64, 2314–2331. [Google Scholar] [CrossRef]
  11. Cavalcanti, M.M.; Domingos Cavalcanti, V.N.; Martinez, P. General decay rate estimates for viscoelastic dissipative systems. Nonlinear Anal. Theory Methods Appl. 2008, 68, 177–193. [Google Scholar] [CrossRef]
  12. Han, X.; Wang, M. Global existence and blow up of solutions for a system of nonlinear viscoelastic equation with damping and source. Nonlinear Anal. Theory Methods Appl. 2009, 71, 5427–5450. [Google Scholar] [CrossRef]
  13. Said-Houari, B.; Messaoudi, S.A.; Guesmia, A. General decay of solutions of a nonlinear system of viscoelastic wave equations. Nonlinear Differ. Equ. Appl. NoDEA 2011, 18, 659–684. [Google Scholar] [CrossRef]
  14. Mustafa, M. Well posedness and asymptotic behavior of a coupled system of nonlinear viscoelastic equations. Nonlinear Anal. Real World Appl. 2012, 13, 452–463. [Google Scholar] [CrossRef]
  15. Liu, W. Uniform decay of solutions for a quasilinear system of viscoelastic equations. Nonlinear Anal. Theory Methods Appl. 2009, 71, 2257–2267. [Google Scholar] [CrossRef]
  16. He, L. On decay of solutions for a system of coupled viscoelastic equations. Acta Appl. Math. 2020, 167, 171–198. [Google Scholar] [CrossRef]
  17. Guesmia, A.; Messaoudi, S.A. A general decay result for a viscoelastic equation in the presence of past and finite history memories. Nonlinear Anal. Real World Appl. 2012, 13, 476–485. [Google Scholar] [CrossRef]
  18. Park, J.Y.; Park, S.H. General decay for a quasilinear system of viscoelastic equations with nonlinear damping. Acta Math. Sci. 2012, 32, 1321–1332. [Google Scholar] [CrossRef]
  19. Feng, B.; Qin, Y.; Zhang, M. General decay for a system of nonlinear viscoelastic wave equations with weak damping. Bound. Value Probl. 2012, 2012, 146. [Google Scholar] [CrossRef]
  20. Araujo, R.D.; Ma, T.F.; Qin, Y. Long-time behavior of a quasilinear viscoelastic equation with past history. J. Differ. Equ. 2013, 254, 4066–4087. [Google Scholar] [CrossRef]
  21. Messaoudi, S.A.; Al-Khulaifi, W. General and optimal decay for a quasilinear viscoelastic equation. Appl. Math. Lett. 2017, 66, 16–22. [Google Scholar] [CrossRef]
  22. Mustafa, M.I. General decay result for nonlinear viscoelastic equations. J. Math. Anal. Appl. 2018, 457, 134–152. [Google Scholar] [CrossRef]
  23. Li, F.; Jia, Z. Global existence and stability of a class of nonlinear evolution equations with hereditary memory and variable density. Bound. Value Probl. 2019, 2019, 37. [Google Scholar] [CrossRef]
  24. Gheraibia, B.; Boumaza, N. General decay result of solutions for viscoelastic wave equation with Balakrishnan-Taylor damping and a delay term. Z. Angew. Math. Phys. 2020, 71, 198. [Google Scholar] [CrossRef]
  25. Kelleche, A.; Feng, B. On general decay for a nonlinear viscoelastic equation. Appl. Anal. 2021, 102, 1582–1600. [Google Scholar] [CrossRef]
  26. Li, C.; Jin, K. General decay results for viscoelastic systems with memory and time-varying delay. Math. Methods Appl. Sci. 2022, 45, 4397–4407. [Google Scholar] [CrossRef]
  27. Youkana, A.; Al-Mahdi, A.M.; Messaoudi, S.A. General energy decay result for a viscoelastic swelling porous-elastic system. Z. Angew. Math. Phys. 2022, 73, 88. [Google Scholar] [CrossRef]
  28. Liang, F.; Gao, H. Exponential energy decay and blow-up of solutions for a system of nonlinear viscoelastic wave equations with strong damping. Bound. Value Probl. 2011, 2011, 22. [Google Scholar] [CrossRef]
  29. Lv, M.; Hao, J. General decay and blow-up for coupled Kirchhoff wave equations with dynamic boundary conditions. Math. Control. Relat. Fields 2021, 13, 303–329. [Google Scholar] [CrossRef]
  30. Brezis, H. Functional Analysis, Sobolev Spaces and Partial Differential Equations; Springer: New York, NY, USA, 2011. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, C.; Wang, C.; Zhao, X.; Lv, Z. Existence and General Energy Decay of Solutions to a Coupled System of Quasi-Linear Viscoelastic Variable Coefficient Wave Equations with Nonlinear Source Terms. Axioms 2023, 12, 780. https://doi.org/10.3390/axioms12080780

AMA Style

Wang C, Wang C, Zhao X, Lv Z. Existence and General Energy Decay of Solutions to a Coupled System of Quasi-Linear Viscoelastic Variable Coefficient Wave Equations with Nonlinear Source Terms. Axioms. 2023; 12(8):780. https://doi.org/10.3390/axioms12080780

Chicago/Turabian Style

Wang, Chengqiang, Can Wang, Xiangqing Zhao, and Zhiwei Lv. 2023. "Existence and General Energy Decay of Solutions to a Coupled System of Quasi-Linear Viscoelastic Variable Coefficient Wave Equations with Nonlinear Source Terms" Axioms 12, no. 8: 780. https://doi.org/10.3390/axioms12080780

APA Style

Wang, C., Wang, C., Zhao, X., & Lv, Z. (2023). Existence and General Energy Decay of Solutions to a Coupled System of Quasi-Linear Viscoelastic Variable Coefficient Wave Equations with Nonlinear Source Terms. Axioms, 12(8), 780. https://doi.org/10.3390/axioms12080780

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop