pth Moment Stability of a Stationary Solution for a Reaction Diffusion System with Distributed Delays

In this paper, the Sobolev embedding theorem, Holder inequality, the Lebesgue contrl convergence theorem, the operator norm estimation technique, and critical point theory are employed to prove the existence of nontrivial stationary solution for p-Laplacian diffusion system with distributed delays. Furthermore, by giving the definition of pth moment stability, the authors use the Lyapunovfunctional method and Kamke function to derive the stability of nontrivialstationary solution. Moreover, a numerical example illuminates the effectiveness of the proposed methods. Finally, an interesting further thought is put forward, which is conducive to the in-depth study of the problem.

In 2018, the asymptotic behavior of a p-Laplacian reaction-diffusion dynamic system is studied in [15]: However, in practical engineering, the time delay is unavoidable, which may lead to chaos and instability of the system [16][17][18][19][20][21][22][23][24]. Thus, in this paper, we are to investigate the delayed p-Laplacian reaction-diffusion dynamic system. In addition, similar research has already begun. In [24], the p-Laplacian diffusion was firstly introduced in a time-delay dynamical system. Inspired by some methods of [24], the authors in [25] employed an impulsive differential inequality lemma to further study the time-delay neural networks with pulse perturbation. In [16], Ruofeng Rao and Shouming Zhong employed the Ekeland variational principle and Lyapunov stability theory to derive a globally exponential pth moment stability criterion for a Markovian jumping T-S fuzzy diffusion system with time-delays: (2) On the other hand, fuzzy logic theory has shown to be an appealing and efficient approach to deal with the analysis and synthesis problems for complex nonlinear systems. Among various kinds of fuzzy methods, Takagi-Sugeno (T-S) fuzzy models provide a successful method to describe certain complex nonlinear systems using some local linear subsystems [26][27][28][29][30][31].

Preliminaries
Consider the following p-Laplacian system with time delays: where N {1, 2, · · · , n} is a finite index set, and Ω ∈ R n is a bounded set with a smooth boundary ∂Ω. {µ jk (j = 1, 2, · · · , N; k = 1, 2, · · · , m)} is a fuzzy set, ω k (t) represents a premise variable, m is the number of premise variables, and N is the number of IF-THEN rules. a i (t, x, u 1 , · · · , u n ) represents a diffusion operator. Time delay τ(t) satisfies τ(t) ∈ (0, τ], where τ > 0 is a constant. b i represents the behavior function, dependent on t and x. c ij and d ij denote the strength of state links between neurons, f i denotes the activation function of the neurons at i, and J i denotes the external input of the neurons at i. Assume in this paper that J i is a bounded quantity. In view of a standard fuzzy inference method, the system (3) can be inferred as follows: where ω(t) = (ω 1 (t), · · · , ω m (t)), and Υ j (ω(t)) is the membership function corresponding to rule j.
In this paper, we assume (H1) for any given q ∈ (1, p), (H2) for any given i ∈ N , there is the corresponding positive number a i such that where u = (u 1 , u 2 , · · · , u n ); (H3) for any given i ∈ N , there is the corresponding positive number F i such that (H4) for any given i ∈ N , there is the corresponding positive numberb i such that (4) is said to be pth moment stable if, for any given ε > 0 and any given initial value Here, the above definition mainly imitates the Definition 7.1 of [41] (Chapter 1).

Main Result
Theorem 1. If the conditions (H1) and (H2) hold, there is a nontrivialstationary solution for fuzzy system (4).
If, in addition, the conditions (H3) and (H4) and the following condition hold: then, for any given p > 1, the nontrivialstationary solution is pth moment stable.
Proof. Denote a i (t) = a i (t, u 1 , · · · , u n ). We shall complete the proof after two steps.
Step 1. We claim that there is a nontrivialstationary solution for fuzzy system (4). Indeed, for any given i ∈ N , we consider the following functional: For convenience, we denote Combining the continuity hypothesis and (H1) results in )dx, we can see it from the differential mean value theorem that the Gateaux differential of where θ ∈ [0, 1]. Now, we try to employ the Lebesgue control convergence theorem to deal with the limit. For |s| 1, the Young inequality derives |g(x, ξ i + θsv)v| in which the function 2 q |v| q is a Lebesgue Integrable function due to Sobolev embedding theorem [42]. In addition, then the Lebesgue control convergence theorem yields On the other hand, it follows from |g(x, ξ i )| C 0 + M 0 |ξ i | q−1 that the operator S : is the bounded continuous operator of L q (Ω) → L q q−1 (Ω). In addition, the Sobolev embedding theorem yields Below, we shall show that the Gateaux derivative DΓ i (ξ i ) is continuous on ξ i . In fact, which implies the norm of the operator This means the operator T : In addition, the Sobolev embedding theorem yields that the operator I : In addition, it follows from [43,44] and the condition of a i (t) that Λ i ∈ C 1 (W 1,p 0 (Ω), R). Therefore, we have proved I i ∈ C 1 (W 1,p 0 (Ω), R). If, for any given i ∈ N , the critical point of I i exists, say ξ * = (ξ * 1 , ξ * 2 , · · · , ξ * n ), a nontrivial stationary solution for the fuzzy system, where ξ * i = ξ * i (x). Below, we shall prove the existence of the critical point. In fact, the condition (H1) yields where C 1 > 0. By the arbitrariness of q, we select a suitable constant q ∈ (1, p) such that Letting C 1 be big enough, we can prove In fact, combining the continuity of b i and the hypothesis (H1) yields that there existsC 1 which deduces the inequality (8) due to the big C 1 . From the Sobolev embedding theorem, we know that there are C 2 , C 3 > 0 such that where ξ i = Ω |∇ξ i | p dx 1 p . In addition, (3.5) and 1 < q < p yield that the lower bound of I i exists. We shall prove that I i is coercive on W 1,p 0 (Ω). Due to η k → ∞ and I i (η k ) C, (9) leads to a contradiction. Hence, I i must be coercive on W (Ω) such that I i (η k ) → c i whenever k → ∞. In addition, we know from the coercivity and the lower bound of the functional I i that c i is the global minimum of I i on the Sobolev space W 1,p 0 (Ω). Moreover, if η k → ξ * i ∈ W 1,p 0 (Ω), then I i (ξ * i ) = c i with I i (ξ * i ) = 0, and ξ * i must be in the bounded subset of W 1,p 0 (Ω). Thus, the minimization sequence {η k } ⊂ W 1,p 0 (Ω) satisfies I i (η k ) → c i and I i (η k ) → 0. In the inequality (9), let ξ i = η k . Then, (9) p   0 (Ω)) * as follows: and It follows from (H1) that there are positive numbers C 4 , C 5 such that From [4,5], the operators A i and A −1 i are continuous, and B i is compact operator. On one hand, and On the other hand, {η k } ∞ k=1 is bounded on the space W 1,p 0 (Ω). Thus, it follows by the reflexivity of W 1,p 0 (Ω) that there is a subsequence of {η k } ∞ k=1 that is weak convergent, say, {η k } ∞ k=1 . Since B is compact operator, there is a subsequence of {η k } ∞ k=1 such that {B i (η k )} ∞ k=1 is convergent. In addition, {A i (η k )} ∞ k=1 owns a convergent subsequence. Moreover, the continuity of A −1 yields that {η k } ∞ k=1 owns a convergent subsequence, say, η k → ξ * i ∈ W 1,p 0 (Ω). Hence, we have proved I i (ξ * i ) = c i and I i (ξ * i ) = 0. By the arbitrariness of i, we have also proved that there is a nontrivial stationary solution ξ * = (ξ * 1 , ξ * 2 , · · · , ξ * n ) for the system (4).
Step 2. To prove that ξ * is pth moment stable.

Remark 1.
In this numerical example, τ = 3 is the bigger upper bound of time delays, which shows the effectiveness of the proposed methods.

Conclusions and Further Considerations
Mainly inspired by some methods and ideas of literature [16,27,40] related to p-Laplace, we employed the critical point theory, variational technique, and Lyapunov functional method to derive the existence theorem of pth moment stable non-trivial stationary solutions for a class of p-Laplacian reaction-diffusion delay systems. The theorem holds for all p > 1, and the methods used in this paper are different from those in previous related literature to some extent. For example, we proved I i ∈ C 1 (W 1,p 0 (Ω), R) while it was not proved in related literature. A numerical example illustrates the effectiveness of the proposed methods. In [45] (Theorem 4.3), Ruofeng Rao and Shouming Zhong proposed a stability criterion for the delayed feedback financial system, in which the pulse effect occurs at a long time (see [45] (Remark 7) for details). How can the impulse control method be used to derive the stability criterion for the p-Laplacian diffusion system (4)? This is an interesting problem.