Uniform Energy Decay Rates for the Fuzzy Viscoelastic Model with a Nonlinear Source

<jats:p>This paper considers the fuzzy viscoelastic model with a nonlinear source <jats:inline-formula>
                     <math xmlns="http://www.w3.org/1998/Math/MathML" id="M1">
                        <msub>
                           <mrow>
                              <mi>u</mi>
                           </mrow>
                           <mrow>
                              <mi>t</mi>
                              <mi>t</mi>
                           </mrow>
                        </msub>
                        <mo>+</mo>
                        <mi>L</mi>
                        <mi>u</mi>
                        <mo>+</mo>
                        <mstyle displaystyle="true">
                           <mrow>
                              <msubsup>
                                 <mo stretchy="false">∫</mo>
                                 <mn>0</mn>
                                 <mi>t</mi>
                              </msubsup>
                              <mrow>
                                 <mi>g</mi>
                                 <mfenced open="(" close=")" separators="|">
                                    <mrow>
                                       <mi>t</mi>
                                       <mo>−</mo>
                                       <mi>ζ</mi>
                                    </mrow>
                                 </mfenced>
                                 <mi mathvariant="normal">Δ</mi>
                                 <mi>u</mi>
                                 <mfenced open="(" close=")" separators="|">
                                    <mrow>
                                       <mi>ζ</mi>
                                    </mrow>
                                 </mfenced>
                                 <mtext>d</mtext>
                                 <mi>ζ</mi>
                              </mrow>
                           </mrow>
                        </mstyle>
                        <mo>−</mo>
                        <msup>
                           <mrow>
                              <mfenced open="|" close="|" separators="|">
                                 <mrow>
                                    <mi>u</mi>
                                 </mrow>
                              </mfenced>
                           </mrow>
                           <mrow>
                              <mi>γ</mi>
                           </mrow>
                        </msup>
                        <mi>u</mi>
                        <mo>−</mo>
                        <mi>η</mi>
                        <mi mathvariant="normal">Δ</mi>
                        <msub>
                           <mrow>
                              <mi>u</mi>
                           </mrow>
                           <mrow>
                              <mi>t</mi>
                           </mrow>
                        </msub>
                        <mo>=</mo>
                        <mn>0</mn>
                     </math>
                  </jats:inline-formula> in a bounded field Ω. Under weak assumptions of the function <jats:inline-formula>
                     <math xmlns="http://www.w3.org/1998/Math/MathML" id="M2">
                        <mi>g</mi>
                        <mfenced open="(" close=")" separators="|">
                           <mrow>
                              <mi>t</mi>
                           </mrow>
                        </mfenced>
                     </math>
                  </jats:inline-formula>, with the aid of Mathematica software, the computational technique is used to construct the auxiliary functionals and precise priori estimates. As time goes to infinity, we prove that the solution is global and energy decays to zero in two different ways: the exponential form and the polynomial form.</jats:p>


Introduction
In this paper, we take the following fuzzy viscoelastic model into account: where the fuzzy number η ∈ (1, 1000), c > 0: In R N (N ≥ 1), Ω is a domain which is well bounded. Besides, the boundary of Ω is smooth perfectly and expressed as Γ: � zΩ. Meanwhile, the memory kernel g (t) is positive and some assumptions will be given in detail.
is type of problems has been observed in many areas of scientific and engineering fields. For example, time analyticity for the viscoelastic equation is studied as follows [1]: S. Berrimi and S. A. Messaoudi applied weak conditions on the memory kernel g. Meanwhile, considering the condition that the energy is positive and relatively small, they obtained the existence of global solutions. Taking into account, they also obtained a decay rate exponentially in [2]. M.M. Cavalcanti and H.P. Oquendo improved this latter result in [3]. In their work, two situations, the internal dissipation and the viscoelastic dissipation, are considered to act on respective part separately. e authors in [3] expanded the internal dissipation to nonlinear cases as much as possible. Simultaneously, the system is well stabilized through the dissipation which is induced by the integral term.
As we know, the viscoelastic terms attract many mathematicians; for instance, the authors in [4] studied the energy decay rate for the global solution of a quasilinear viscoelastic model, and the authors in [5][6][7][8][9] considered time analysis of solutions for some viscoelastic models. Moreover, F. Y. Zhang et al. considered a nonlinear viscoelastic equation in [10]: e solution is perfectly stabilized through the dissipation, which is induced by the viscoelastic term. e modified energy functional in [10] has been used to prove the energy decay through two different ways: the exponential form and the polynomial form. Additionally, the construction of auxiliary functions is organized by the computational technique of undetermined coefficients. Besides, the authors in [11,12] discussed the adaptive fuzzy control of nonlinear systems, and the discussion of fuzzy coefficients involved in these papers is quite interesting.
Inspired by these works, we consider (1) in this paper, the two optimal decay rates, exponential decay and polynomial decay, are easily and directly established through the application of Mathematica software. e specific arrangement of this work is as follows: in Section 2, we present some notations and necessary materials; in Section 3, in view of the fuzzy number ƞ, we give the whole decay result, and our choice of the "Lyapunov" functional shows the extensive applicability and practical significance of the computational technique.

Preliminaries
In this section, L p (Ω) and H 1 0 (Ω) are understood and applied in their usual senses. We impose the following hypotheses and preliminaries on the memory kernel g (t). In addition, the definition of energy function plays a significant role to our main result: where both g (0) and ℓ are positive. H 2 : the existence of a positive constant ξ makes the following formula hold: Remark 1. From the assumption above, if p � 1, we have where d 1 � (p − 1)ξ > 0 and d 2 � g 1− p (0). Indeed, the condition p < (3/2) plays an important role to ensure that , then u (t) can be found as a unique solution to model (1) with Simultaneously, we get Taking the initial data, less regularity, and the priori estimate into consideration, the theorem above guarantees the existence of solutions for model (1) as in [7,13]. In addition, the Galerkin approximation method can be applied to accomplish the proof of the theorem above.
Our primary assignment is to find out the energy function ε(t). Combining the multiplier method, the integral subsection integration, and (H 1 , H 2 ), the calculation is provided: 2 Journal of Mathematics Now, (13) yields Now, it is obvious that the energy of model (1) can be given directly: where ereby, we can find which means ε(t) is decreasing.

Main Results
In this chapter, the main result of this work is put forward directly and clearly. We introduce the following auxiliary functionals firstly: where the coefficients α and β will be determined later.

Lemma 1.
e existence of a positive constant C makes the following conclusion: Proof. Recalling the inequalities of Cauchy and Poincaré, the following estimates can be arrived: where the symbol C p > 0 satisfies ‖u‖ 2 ≤ C p ‖∇u‖ 2 for all u ∈ H 1 0 (Ω). Here, it is necessary to indicate that the positive constants C and C p later in this article represent different meanings in different places.
From the representation of ε(t), it is naturally for us to get us, (20) follows. □ Lemma 2. For any t ≥ 0, the following estimate holds: where ε 1 and ε 2 are selected appropriately.
Proof. Deriving function φ(t) with respect to time, the following calculation will be given clearly: Considering the fourth item of (23), we adopt the following estimate from [10] without proof: Taking the sixth item of (23) into account, we adopt the following estimate from [1] as follows: For the last item of (23), we have us, (21) follows.
Proof. Deriving function χ(t) with respect to time, the expression will be given distinctly: (29)

Journal of Mathematics
For any α, β ≥ 0, we will estimate every item of the right hand of (29): Similar to the first term, considering the third item of (29), similar estimate will be given: Applying Young's inequality properly, the following estimates can be obtained from [10]: Taking (29)-(39) and H 2 into account, for all t ≥ 0, we get (28). □ Theorem 2. Let every pair (u 0 , u 1 ) ∈ H 1 0 (Ω) × H 1 0 (Ω), H 1 and H 2 hold, and t 0 ∈ (0, ∞) For every t ∈ [t 0 , ∞), there must be some positive constants K 1 , K 2 , and k, which would enable the solution of the model (1) to satisfy the following: Proof. We start the proof by selecting an appropriate auxiliary functional Journal of Mathematics where the positive constants M, a, and b will be chosen in the sequel.
Applying Lemma 1 and the hypothesis g ′ (t) ≤ − ξ · g p (t), the equivalence between ε(t) and L(t) is achieved. To this point, firstly, a simple calculation shows that By (25), we get 8 Journal of Mathematics us, the equivalence between L(t) and ε(t) is complete.
Secondly, differentiating (41) with respect to t, considering the assumption H 2 and Lemmas 1 and 2, the estimate of L ′ (t) is given naturally: where the coefficients c 1 , c 2 , c 3 , and c 4 are given as follows: Next, it is time to carry out a discussion on the different values of p.
In such a case, choosing we find (47) erefore, (45) yields Remark 2. Using the computational technique proposed in [10] and selecting some appropriate numerical values by means of the calculation software Mathematica, it is obvious that Considering the interval of the fuzzy number η, taking for example, with the aid of Mathematica, we easily verify that Taking another example into consideration, we find that Additionally, we may retrieve more effective data of c i , and the data provide an intuitive understanding in the discussion of the auxiliary functionals.
Considering the equivalence between L(t) and ε(t), it is easy for us to see that where the coefficients M 1 and M 2 in (53) have many possibilities to be chosen. e following expression can be presented through the combination of (48) and (53): After a simple integral on the interval (t 0 , t), (54) leads to As a consequence, taking k � − (h/M 1 ), (55) leads to where K 1 is some appropriate positive constant.  (9), we easily verify that Moreover, considering the inequality in [6], for some constant C, Taking θ � (1/2), c � 2p − 1, we find (cθ/p − 1 + θ) � 1. erefore, for some m > 1, we get that is, (60) Taking (48) into account, we have For each t ≥ t 0 , it is straight for us to get As we mentioned before, the symbol C denotes different constants in different places. Executing a simple integral on the interval (t 0 , t), (62) yields L(t) ≤ C a t + C b − (1/m− 1) � C a t + C b − (1/2(p− 1)) . (63) Considering the equivalence between L(t) and ε(t), it is apparently for us to get ε(t) ≤ K 2 (t + 1) (1/2(p− 1)) .
is completes the proof.

Conclusion
Under the assumptions on the relaxation function and the interval of the fuzzy number η, applying the computational technique, a lot of auxiliary functionals can be constructed numerically. Two decay results, the exponential one and the polynomial one, are derived for the model (1) eventually. e result shows a new way for the decay rates, which is quite different from other literatures.

Data Availability
All datasets generated for this study are included in the manuscript.

Conflicts of Interest
e author declares that there are no conflicts of interest.