Stability of stochastic 2D Navier-Stokes equations with memory and Poisson jumps

1. Introduction

In this paper, we will investigate the stability of the weak solutions of stochastic 2D Navier-Stokes equations with memory and Poisson jumps of the form:

where \(X\) is the velocity field of the fluid, \(\nu>0\) is the kinematic viscosity, \(p\) is the associated hydrostatic pressure, \(f\) is a non-delayed external force field, \(g(X(t-\rho(t)))\) is another external force field with delay, \(\sigma(t,X(t-\delta(t)))dW(t)+\int_{Z}k(t,X(t-\gamma(t)),z)\widetilde{\eta}(dt,dz)\) is a random external force field with delays, where \(W\) is an infinite dimensional Wiener process and \(\widetilde{\eta}\) is a compensated time homogeneous Poisson random measure, \(\rho,\delta,\gamma: \mathbb{R}^+\rightarrow[0,r]\) are continuous functions, \(D\) is a regular open bounded domain of \(\mathbb{R}^2\) with boundary \(\partial D\), \(X_0\) is the initial velocity, \(\phi\) is the initial datum in the interval \([-r,0]\).

Navier-Stokes equations are the fundamental model of the fluid mechanics and turbulence. These equations have been the object of numerous works (see, for instance [1,2] and the references therein), even with unbounded domains (see [2,3]) since the first paper of Leray was published in 1933 [4]. Beside, noise or stochastic perturbation is unavoidable and omnipresent in nature as well as in man-made systems. Therefore, it is of great significance to import the stochastic effects into the investigation of Navier-Stokes equations. To our best knowledge, the theory of stochastic Navier-Stokes equations apparently has its roots in the 1959 edition of the Landau and Lifshitz [5] and the first work on the stochastic Navier-Stokes equations written from the mathematical point of view is the paper [6]. We mention here the works ([7,8,9,10] and the references cited therein) concerning the two-dimensional (2D in short) stochastic Navier-Stokes equation.

Partial differential equations with memory (e.g. delay) have attracted great interest due to their applications in describing many sophisticated dynamical systems in physics, chemistry, biology, economics and social sciences. On this matter, we refer the reader to [11,12,13,14,15,16] and the references therein. In the past few years, many papers have studied the Navier-Stokes equations with a forcing term which contains some hereditary characteristics (see, for instance [17,18,19] and the references therein). In addition, the long-time behavior and exponential stability of the Navier-Stokes equations is an interesting and challenging problem, since it can provide useful information on the future evolution of the system (see, [20,21]). On the exponential behavior for stochastic 2D Navier-Stokes equations with variable delay, we refer the reader to [22,23,24], and recently, in [25] and [26], the authors investigated the asymptotic behavior of solutions of stochastic evolution equations for second grade fluids and non-Newtonian fluids, respectively.

On the other hand, the world is more complicated and models which are allowed to have jumps – both big and small – are desirable. Hence, it is necessary and important to study stochastic systems with Poisson jumps or Lévy processes (see, for instance [27,28,29,30,31] and the references therein). In recent years, the stochastic 2D(3D) Navier-Stokes equations with Lévy noise has attracted much attention of researchers. To be more precise, in [32], Motyl considered the existence of solutions, in a probabilistic sense, to the stochastically driven, viscous, incompressible Navier-Stokes equations driven by the Lévy noise consisting of the compensated time homogeneous Poisson random measure and the Wiener process in two and three spatial dimensions, Dong et al., [33] discussed the existence of stationary weak solutions of stochastic 3D Navier-Stokes equations involving jumps and compared the Galerkin stationary probability measures for the case driven by Lévy noise and the one driven by Wiener process, and in [34], by using an abstract setting, Brzezniak et al. studied the existence and uniqueness of the solution of an abstract nonlinear equation driven by a multiplicative noise of Lévy type. More specifically, just recently, Taniguchi [35] obtained the existence and exponential stability of energy solutions to 2D stochastic functional Navier-Stokes equation perturbed by the Lévy process. However, to prove the exponential stability results in [35], the author have imposed that the delay function \((\rho(t))\) appears only in the non-random external force and furthermore this delay function have required a strong assumption, i.e. \(\rho(t)\) is differentiable and satisfies \(0\leq \rho'(t)< M_{\ast}0\). Thus, we will make the first attempt to establish some results for more general forcing terms and relax this restriction. Under some suitable assumptions, by using the Itô formula for jumps and Burkholder-Davis-Gundy inequalities in stochastic analysis, we establish the weak solution to (1) converges to the stationary solution of its stationary versions exponentially stable in the mean square. Further, by establishing a lemma for compensated Poisson random measures (Lemma 3), we give the result that the weak solution to (1) converges to the stationary solution of its stationary version almost surely exponentially. Assumptions given in our article do not require the monotone decreasing behaviour of the delays and satisfying \(0\leq \rho'(t),\delta'(t),\gamma'(t)< 1\). Therefore, the current paper can be regarded as the extension of the work of Caraballo and Real [20] to the stochastic settings, simultaneously extend and improve the one of Taniguchi [35] (where the random external force field does not include delays and the external force field contains delay but the memory function \(\rho(t)\) is a differentiable function) as well as the asymptotic behavior results published in [7] (where the external force fields does not contain delays and the stochastic Navier-Stokes equation without non-Gaussian Lévy noise perturbation) and the papers announced by Chen [22], Wan and Zhou [23] (where the random external force field does not contain discontinuous multiplicative noise).

The remainder of this paper is organized as follows: In Section 2, we briefly present some basic notations, preliminaries. The main results in Section 3 are devoted to studying the asymptotic behavior for the weak solutions of the system (1) with their proofs.

2. Preliminaries

In this section, we introduce notations and preliminary results need to establish our results. For more details on this section, we refer the reader to [2,8,36,37,38].

We first introduce the following function spaces, which are usual in the study of Navier-Stokes equations:

\[\mathcal{V}:=\{u\in C_{0}^{\infty}(D,\mathbb{R}^2):\text{div \(u =0\)}\}.\] \(\mathbb H:=\) the closure of \(\mathcal{V}\) in \(L^2(D,\mathbb{R}^2)\) with the norm \(|u|=(u,u)^{\frac{1}{2}}\), where for \(u,v\in L^2(D,\mathbb{R}^2),\) \[(u,v)=\sum_{j=1}^{2}\int_{D}u_j(x)v_j(x)dx.\] \(\mathbb V:=\) the closure of \(\mathcal{V}\) in \(\mathbb H_{0}^{1}(D,\mathbb{R}^2)\) with the norm \(\|u\|=((u,u))^{\frac{1}{2}}\), where for \(u,v\in \mathbb H_{0}^{1}(D,\mathbb{R}^2),\) \[((u,v))=\sum_{i,j=1}^{2}\int_{D}\frac{\partial u_i}{\partial x_j}\frac{\partial v_i}{\partial x_j}dx.\] It follows that \(\mathbb H\) and \(\mathbb V\) are separable Hilbert spaces with associated inner products \((\cdot,\cdot)\) and \(((\cdot,\cdot))\) and the following is satisfied \[\mathbb V\subset \mathbb H\equiv \mathbb H^\ast\subset \mathbb V^\ast,\] where injections are dense, continuous, and compact; \(\mathbb H^\ast\) and \(\mathbb V^\ast\) stand for the topological dual of \(\mathbb H\) and \(\mathbb V\), respectively.

Now, let \(P_{ \mathbb H}\) be an orthogonal projector from \(L^2(D,\mathbb{R}^2)\) onto \(\mathbb H\). We define the operator \(A:L^2(D,\mathbb{R}^2)\rightarrow \mathbb H\) by \(Au=-P_{ \mathbb H}\triangle u\) and \(B:\mathbb V\times \mathbb V\rightarrow\mathbb V^\ast\) by \(\langle B(u,v),w\rangle=b(u,v,w)\), \(\forall u,v,w\in \mathbb V,\) where \(\langle\cdot,\cdot\rangle\) denotes the duality \(\langle\mathbb V^\ast,\mathbb V\rangle\) and \(b\) is the trilinear form defined by

\[b(u,v,w)=\sum_{i,j=1}^{2}\int_{D}u^j(x)\frac{\partial v^j}{\partial x_i}w^j(x)dx.\] We also set \(B(u)=B(u,u)\), \(\forall u\in \mathbb V.\)

Furthermore, we shall need some properties of the trilinear form \(b\), and we list below the ones that we will used later on (see, [38]),

\begin{align*} &|b(u,v,w)|\leq C_1|u|^{\frac{1}{2}}\|u\|^{\frac{1}{2}}\|v\||w|^{\frac{1}{2}}\|w\|^{\frac{1}{2}}, \quad \forall u,v,w\in \mathbb V,\&b(u,v,v)=0,\quad \forall u,v\in \mathbb V,\\ &b(u,u,v-u)-b(v,v,v-u)=-b(v-u,u,v-u),\quad \forall u,v\in \mathbb V, \end{align*} where \(C_1>0\) is an appropriate constant which depends on the regular open domain \(D\) (see [39]).

Let \((\Omega,\mathcal{F},\mathbf{P})\) be a complete probability space equipped with some filtration \((\mathcal F_{t})_{t\geqslant 0}\) satisfying the usual conditions (i.e., it is right continuous and \(\mathcal F_{0}\) contains all \(\mathbf{P}\)-null sets).

With the symbol \(\{W(t)\}_{t\geq 0}\), we denote a \(\mathbb{K}\)-valued \((\mathcal F_{t})_{t\geqslant 0}\)-Wiener process defined on the probability space \((\Omega,\mathcal{F},\mathbf{P})\) with covariance operator \(Q\), i.e.

\[\mathbf{E}\langle W(t),a\rangle_{\mathbb{K}}\langle W(s),b\rangle_{\mathbb{K}}=\min\{t,s\} \langle Qa,b\rangle_{\mathbb{K}},\quad \forall a, b\in\mathbb{K}, \] where \(Q\) is a positive, self-adjoint, trace class operator on \(\mathbb{K}\). In particular, we call such \(\{W(t)\}_{t\geq 0}\) a \(\mathbb{K}\)-valued \(Q\)-Wiener process relative to \((\mathcal F_{t})_{t\geqslant 0}\). We assume that there exist a complete orthonormal system \(\{e_n\}_{n\in \mathbb{N}}\) in \(\mathbb{K}\), a bounded sequence of nonnegative real number \(\{\lambda_n\}_{n\in \mathbb{N}}\) such that \(Qe_{n}=\lambda_{n}e_n,\quad n=1,2,3, …\), and a sequence \(\{\beta_{n}\}_{n\geq 1}\) of independent standard Brownian motions such that \[\langle W(t),e\rangle_{\mathbb{K}}=\Big\langle\sum_{n=1}^{\infty}\sqrt{\lambda_{n}} e_n\beta_{n},e\Big\rangle_ {\mathbb{K}},\quad t\geq 0,\quad e\in \mathbb{K},\] and let \(\mathcal F_{t}=\sigma\{W(s): 0\leq s\leq t\}\) be the \(\sigma\)-algebra generated by \(W\).

In order to define stochastic integrals with respect to the \(Q\)-Wiener process \(W(t)\), we introduce the subspace \(\mathbb K_{0}=Q^{\frac{1}{2}}\mathbb{K}\) of \(\mathbb{K}\), which endowed with the inner product,

\[\langle a,b\rangle_{\mathbb K_{0}}=\langle Q^{-\frac{1}{2}}a,Q^{-\frac{1}{2}}b\rangle_{\mathbb{K}},\quad \forall a, b\in\mathbb{K}_{0}\] is a Hilbert space. Let \({\mathcal L_{2}^{0}}=\mathcal L_{2}(\mathbb K_{0};\mathbb{H})\) denote the space of all Hilbert-Schmidt operators from \(\mathbb K_0\) into \(\mathbb{H}\). It turns out to be a separable Hilbert space, equipped with the norm \[\|\psi\|_{{\mathcal L_{2}^{0}}}^{2}=tr((\psi Q^{\frac{1}{2}})(\psi Q^{\frac{1}{2}})^{\ast}),\] for any \(\psi \in {\mathcal L_{2}^{0}}\). Obviously, for any bounded operators \(\psi \in\mathcal L(\mathbb K;\mathbb{H})\) – the set of all linear bounded operators from \(\mathbb{K}\) into \(\mathbb{H}\), this norm reduces to \[\|\psi\|_{{\mathcal L_{2}^{0}}}^{2}=tr(\psi Q\psi^{\ast})=\sum_{n=1}^{\infty}|\sqrt{\lambda_{n}}\psi e_{n}|^{2}.\] Let \(\Phi : (0,\infty)\rightarrow{\mathcal L_{2}^{0}}\) be a predictable, \(\mathcal F_t\)-adapted process such that \[\int_{0}^{t}\mathbf{E}\|\Phi (s)\|_{\mathcal L_{2}^{0}}^{2}ds0.\] Then, we can define the \(\mathbb{H}\)-valued stochastic integral \(\int_{0}^{t}\Phi(s)dW(s)\) (which is a continuous square-integrable martingale) of \(\Phi\) with respect to the \(\mathcal F_{t}\)-valued \(Q\)-Wiener process \(W(t)\) by \[\Big\langle\int_{0}^{t}\Phi(s)dW(s),e\Big\rangle:=\sum_{n=1}^{\infty}\int_{0}^{t}\Big\langle\Phi(s)\sqrt{\lambda_{n}} e_n,e\Big\rangle d\beta_{n}(s)\] for any \(e\in \mathbb H\) using the Itô integral with respect to \(\beta_{n}(s)\). For the construction, we can refer to Da Prato and Zabczyk [8].

Let \(L=(L_t)_{t\geq 0}\) is a \(\mathbb K\)-valued Lévy processes such that \(L\) has stationary and independent increments, is stochastically continuous and satisfies \(L_0=0\) almost surely. Let \(p(t)\), \(t\geq 0\) be the law of \(L_t\), then \((p(t))_{t\geq 0}\) is a weakly continuous convolution semigroup of probability measures on \(\mathbb K\). We have the Lévy-Khinchin formula [36] which yields for all \(t\geq 0\), \(x\in \mathbb K\),

\[\mathbf E\Big(e^{i\langle x,L_t\rangle_{\mathbb K}}\Big)=e^{t\zeta(x)},\] where \begin{align*} \zeta(x):=\exp\Big\{i\langle a,x\rangle_{\mathbb K}-\frac{1}{2}\langle Qx,x\rangle_{\mathbb K} +\int_{\mathbb K-\{0\}}\big[e^{i\langle x,y\rangle_{\mathbb K}}-1-\chi_{\{|y|_{\mathbb K}< 1\}}(y)i\langle x,y\rangle_{\mathbb K}\big]\lambda(dy)\Big\}, \end{align*} where \(a\in \mathbb K\) and \(\lambda\) is a Lévy measure or a jump intensity measure of \(L\) on \(\mathbb K-\{0\}\), i.e., \(\int_{\mathbb K-\{0\}}\min(|y|_{\mathbb K}^{2},1)\lambda(dy)< \infty\); \(\chi_Z\) denotes the characteristic function on set \(Z\subset \mathbb K\); the triple \((a,Q,\lambda)\) is the characteristics of \(L\) and the mapping \(\zeta\) is the characteristic exponent of \(L\). We can also define the Lévy measure on the whole of \(\mathbb K\) via the assignment \(\lambda(\{0\})=0\).

Now, we shall write \(\triangle L_t:=L_t-L_{t-}\), \(\forall t\geq 0,\) where \(L_{t-}:=\lim_{s\uparrow t}L_s.\) Then, almost surely for any \(Z\in\mathcal{B}(\mathbb K-\{0\})\), which denotes the Borel \(\sigma\)-field of \((\mathbb{K}-\{0\})\) and with \(0\notin\) the closure of \(Z\), we get a counting Poisson random measure \(\eta\) on \((\mathbb{K}-\{0\})\):

\[\eta(t,Z)=\# \{0\leq s\leq t,\triangle L_s\in Z\}< \infty,\quad t\geq 0.\] Let \[\widetilde{\eta}(t,dy):=\eta(t,dy)-\lambda(dy)t\] be the compensated Poisson measure that is independent of \(W(t).\)

Let \(\lambda_Z\) denotes the restriction of the measure \(\lambda\) to \(Z\), still denoted by \(\lambda\), such that \(\lambda\) is finite on \(Z\). Denote by \(\mathcal{ P}^{2}([0,T]\times Z;\mathbb{H})\) the space of all predictable mappings \(k:[0,T]\times Z\rightarrow\mathbb{H}\) for which

\[\int_{0}^{T}\int_{Z}\mathbf{E}|k(t,y)|_{\mathbb{H}}^{2}\lambda(dy)dt< \infty.\] We may then define the \(\mathbb{H}\)-valued stochastic integral \[\int_{0}^{T}\int_{Z}k(t,y){\eta}(dt,dy):=\sum_{0\leq t\leq T}k(t,\triangle Y_t)\chi_{Z}(\triangle Y_t),\] where \[Y_t:=\int_{Z}y\eta(t,dy)=\sum_{0\leq t\leq t}\triangle L_s\chi_{Z}(\triangle L_t)\] as a random finite sum which enables us to define \[\int_{0}^{T}\int_{Z}k(t,y)\widetilde{\eta}(dt,dy):=\int_{0}^{T}\int_{Z}k(t,y){\eta}(dt,dy)-\int_{0}^{T}\int_{Z}k(t,y))\lambda(dy)dt.\] Furthermore, we can see that \(\int_{0}^{t}\int_{Z}k(s,y)\widetilde{\eta}(ds,dy)\) is an \(\mathbb H\)-valued centered square-integrable martingale such that \[\mathbf E \Big(\Big|\int_{0}^{T}\int_{Z}k(t,y)\widetilde{\eta}(dt,dy)\Big|_{\mathbb H}^{2}\Big) =\int_{0}^{T}\int_{Z}\mathbf{E}|k(t,y)|_{\mathbb{H}}^{2}\lambda(dy)dt.\] We can refer to Protter [37] for a systematic theory about stochastic integrals of this kind.

Thus, the stochastic 2D Navier-Stokes equations with with memory and Poisson jumps (1) can be rewritten as follows in the abstract mathematical setting:

where \(f\in \mathbb V^\ast\), \(g:\mathbb V\rightarrow\mathbb V^\ast\), \(\sigma:[0,+\infty)\times \mathbb V\rightarrow{\mathcal L_{2}^{0}}(\mathbb K,\mathbb H)\), and \(k:[0,+\infty)\times\mathbb H\times(\mathbb{K}-\{0\})\rightarrow\mathbb H\) are progressively measurable; \(L^{2}(\Omega,\mathcal{C}([-r,0],\mathbb H))\) denotes the family of all almost surely bounded \((\mathcal F_{t})_{t\geqslant 0}\)-measurable and \(\mathcal{C}([-r,0],\mathbb H)\)-valued stochastic process, here \(\mathcal{C}([-r,0],\mathbb H)\) denote the family of all right-continuous functions with left-hand limits \(\phi\) from \([-r,0]\) to \(\mathbb H\) which is equipped with the norm \(\|\phi\|_{\mathcal{C}}:=\sup_{-r\leq \theta\leq 0}|\phi(\theta)|\).

Now, we give the definition of the weak solution of system (2).

Definition 1. An \(\mathcal F_{t}\)-adapted process \(X(t)\) is called the weak solution to (2) if the following conditions are satisfied:

• (i) \(X(t)\in \mathcal{C}(-r,T;\mathbb H)\cap L^{2}(-r,T;\mathbb V)\), a.s., \(\forall T>0\);
• (ii) the following integral equation holds as an identity in \(\mathbb V^\ast\) a.s., \(\forall t\in [0,T]\),
  \begin{align} \label{eq2.2} X(t)=&X(0)+\int_{0}^{t}\big[-\nu AX(s)-B(X(s))+f(s)+g(X(s-\rho(s)))\big]ds\nonumber\\ &\;\;\;+\int_{0}^{t}\sigma(s,X(s-\delta(s)))dW(s)+\int_{0}^{t}\int_{Z}k(s,X(s-\gamma(s)),z)\widetilde{\eta}(ds,dz). \end{align}

  (3)

For our purpose, we recall the Itô formula, which will play a key role in what follows. Let \(C^{2}(\mathbb H;\mathbb{R}^+)\) denote the space of all real-valued nonnegative functions \(\Upsilon\) on \(\mathbb H\) with the following properties:

• (a) \(\Upsilon(x)\) is twice (Fréchet) differentiable in \(x\);
• (b) Both \(\Upsilon_{x}(x)\) and \(\Upsilon_{xx}(x)\) are continuous in \(\mathbb H\) and \(\mathcal{L}(\mathbb H)\).

Lemma 1.[40] Suppose \(\Upsilon\in C^{2}(\mathbb H;\mathbb{R}^+)\) and \(X(t)\), \(t\geq 0\), is a weak solution to (2). Then \begin{align*} \Upsilon(X(t))=&\Upsilon(X(0))+\int_{0}^{t}\mathcal{L}\Upsilon(X(s))ds+\int_{0}^{t}\langle\Upsilon_{x}(X(s)),\sigma(s,X(s-\delta(s)))dW(s)\rangle \nonumber\\ &\;\;\;+\int_{0}^{t}\int_{Z}[\Upsilon(X(s)+k(s,X(s-\gamma(s)),z))-\Upsilon(X(s))]\widetilde{\eta}(ds,dz), \end{align*} where \(\mathcal{L}\) is the associated diffusion operator defined, for any \(x\in \mathbb V\), by \begin{align*} \mathcal{L}\Upsilon(x(t))= &\langle-\nu Ax(t)-B(x(t))+f(t)+g(x(t-\rho(t))),\Upsilon_{x}(x(t))\rangle\\ &\;\;\;+\frac{1}{2}trace\Big(\Upsilon_{xx}(x(t))\sigma(t,x(t-\delta(t)))Q\sigma(t,x(t-\delta(t)))^{\ast}\Big)\\ &\;\;\;+\int_{Z}[\Upsilon(x(t)+k(s,x(s-\gamma(s)),z))-\Upsilon(x(t))-\langle\Upsilon_{x}(x(t)),k(s,x(s-\gamma(s)),z)\rangle]\lambda(dz). \end{align*}

Definition 2. We say that a weak solution \(X(t)\) of system (2) converges to \(u_\infty\in\mathbb H\) exponentially in the mean square if there exist \(c>0\) and \(M_0>0\) such that for all \(t\geq 0\) \[\mathbf E|X(t)-u_\infty|^2\leq M_0e^{-ct}.\] In particular, if \(u_\infty\) is a solution to (2), then it is said that \(u_\infty\) is exponentially stable in the mean square provided that every weak solution to (2) converges to \(u_\infty\) exponentially in the mean square with the same exponential order \(c>0\).

Definition 3. We say that a weak solution \(X(t)\) of system (2) converges to \(u_\infty\in\mathbb H\) almost surely exponentially if there exist \(\gamma>0\) such that \[\lim_{t\to\infty}\frac{\log |X(t)-u_\infty|}{t}\leq -\gamma,\quad a.s..\] In particular, if \(u_\infty\) is a solution to (2), then it is said that \(u_\infty\) is almost surely exponentially stable provided that every weak solution to (2) converges to \(u_\infty\) almost surely exponentially with the same constant \(\gamma\).

3. Main results

In this section, we will discuss the asymptotic behavior for the weak solutions of stochastic 2D Navier-Stokes equations with finite memory and Poisson jumps.

Let \(\lambda_1>0\) be the first eigenvalue of \(A\), then

\[\|v\|^2\geq \lambda_1|v|^2,\quad\forall v\in V.\] To investigate the asymptotic behavior for the weak solutions of (2), we assume the following hypotheses:

• (\(\mathbf H1\)) \(g(0)=0\) and there exists a positive number \(C_g\) such that \[\|g(u)-g(v)\|_{\mathbb V^\ast}\leq C_g|u-v|,\quad \forall u,v\in \mathbb H.\]
• (\(\mathbf H2\)) There exist integrable functions \(\alpha_1,\gamma_1:[0,\infty)\rightarrow \mathbb{R}^+\) such that, for certain constant \(\beta_1\geq 0\) and \(u\in \mathbb H\), \[\|\sigma(t,u)\|_{\mathcal{L}_{2}^{0}}^{2}\leq \alpha_1(t)+(\beta_1+\gamma_1(t))|u-u_\infty|^2.\]
• (\(\mathbf H3\)) There exist integrable functions \(\alpha_2,\gamma_2:[0,\infty)\rightarrow \mathbb{R}^+\) such that, for certain constant \(\beta_2\geq 0\) and \(u\in \mathbb H\), \[\int_{Z}|k(t,u,z)|^2\lambda(dz)\leq \alpha_2(t)+(\beta_2+\gamma_2(t))|u-u_\infty|^2.\]
• (\(\mathbf H4\)) There exists \(\theta>0\) such that, for \(t\geq 0\), \[\int_{0}^{\infty}e^{\theta s}\alpha_i(t)dt< \infty,\quad \int_{0}^{\infty}e^{\theta s}\gamma_i(t)dt< \infty,\quad i=1,2.\]

We first consider the existence of the stationary solution to the equation We have the following lemma:

Lemma 2.[20] Suppose that \(g\) satisfies the condition \((\mathbf {H1})\) and \(\nu>\lambda_{1}^{-1}C_g\). Then we have the following:

• (i) For all \(f\in \mathbb V^\ast\), there exists a stationary solution \(u_\infty\) to (4).
• (ii) There exists a constant \(C(D)>0\) such that, if \((\nu-\lambda_{1}^{-1}C_g)^2>C(D)\|f\|_{\mathbb V^\ast}\), then the stationary solution to (4) is unique.

Now, using above lemma, we will discuss the asymptotic behavior for the weak solutions of (2). Hence, throughout this paper we assume that there exists a unique stationary solution \(u_\infty\in \mathbb V\) to (4).

Set \(y(t):=X(t)-u_\infty\) and \(\Upsilon(y(t))=y^{2}(t)\). Then the function \(y(t)\) satisfies the following equation: \begin{align*} d(X(t)-u_\infty)&=[-\nu A(X(t)-u_\infty))-(B(X(t)-B(u_\infty))+(g(X(t-\rho(t)))-g(u_\infty))]dt\\ &\quad+\sigma(t,X(t-\delta(t)))dW(t) +\int_{Z}k(t,X(t-\gamma(t)),v)\widetilde{\eta}(dt,dv). \end{align*} Similar to the articles of Caraballo et al., [7], Liu [41] and Wan [42], in which, they studied Gaussian white noise by Itô formula, utilizing Itô formula for L\(\acute{e}\)vy noise in Lemma 1 to the function \(\Upsilon(y(t))=y^{2}(t)\) and taking expectation, we easily obtain the following result.

Theorem 1. Suppose that the conditions \((\mathbf {H1})-(\mathbf {H4})\) hold. Then there exists a unique weak solution \(X(t)-u_\infty\in \mathcal{C}(-r,T;\mathbb H)\cap L^{2}(-r,T;\mathbb V)\) a.s.. Furthermore, the following identity holds: \begin{align*} \frac{d}{dt}\mathbf E|X(t)-u_\infty|^2 &=-2\mathbf E\langle\nu A(X(t)-u_\infty),X(t)-u_\infty\rangle -2\mathbf E\langle B(X(t)-B(u_\infty),X(t)-u_\infty\rangle\\ &\;\;\;+2\mathbf E\langle g(X(t-\rho(t)))-g(u_\infty),X(t)-u_\infty\rangle+\mathbf E\|\sigma(t,X(t-\delta(t)))\|_{\mathcal{L}_{2}^{0}}^{2}\\ &\;\;\;+\mathbf E\int_{Z}|k(t,X(t-\gamma(t)),z)|^{2}\lambda(dz). \end{align*} The first main result of this section is the following theorem.

Theorem 2. Suppose that the conditions \((\mathbf {H1})-(\mathbf {H4})\) hold. Then the weak solution \(X(t)\) to (2) converges to the stationary solution \(u_{\infty}\) to (4) exponentially stable in the mean square provided that the following inequality

holds.

Proof. From (5), there exists a positive constant \(\upsilon\) such that \begin{align*} C(D)\sqrt{\lambda_1}\|u_\infty\|+C_g+\frac{\beta_1+\beta_2}{2} < C(D)\sqrt{\lambda_1}\|u_\infty\|+\frac{\beta_1+\beta_2}{2}+\frac{1}{2\upsilon}+\frac{\upsilon C_{g}^{2}}{2}< \nu\lambda_1. \end{align*} Furthermore, there exists a constant \(c\in(0,\theta)\) sufficiently small such that

For convenience, we shall denote \[\alpha(t)=\alpha_1(t)+\alpha_2(t),\quad \eta(t)=(\gamma_1(t)+\gamma_2(t))e^{cr}.\] Since \(\alpha_i,\gamma_i\), \(i=1,2\), is integrable, together with the assumption \((\mathbf {H4})\), we deduce that Define the function \begin{equation*} F(t):= \begin{cases} \mathbf E|X(t)-u_\infty|^2e^{ct}\exp\big(-\int_{0}^{t}[\eta(s)+\alpha(s)e^{cs}]ds\big),\quad t\geq 0,\\ \mathbf E|X(t)-u_\infty|^2e^{ct},\quad t\in[-r,0). \end{cases} \end{equation*} Clearly, \(F(t)\) is right-continuous functions with left-hand limits on \([-r,+\infty)\) and Now, from definition on the operator \(B\), we get and from the properties on trilinear form \(b\), we obtain Furthermore, by Young’s inequality and the assumption \((\mathbf {H1})\), we have Hence, \begin{align*} \frac{dF(t)}{dt} &\leq e^{ct}\exp\Big(-\int_{0}^{t}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\Bigg\{\Big[ c-\eta(t)-\alpha(t)e^{ct}-2\nu\lambda_1 \\ &\quad+2C(D)\sqrt{\lambda_1}\|u_\infty\|+\frac{1}{\upsilon}\Big]\mathbf E|X(t)-u_\infty|^2+\upsilon C_{g}^{2}\mathbf E|X(t-\rho(t))-u_\infty|^2\\ &\quad+\mathbf E\|\sigma(t,X(t-\delta(t)))\|_{\mathcal{L}_{2}^{0}}^{2}+\mathbf E\int_{Z}|k(t,X(t-\gamma(t)),z)|^{2}\lambda(dz)\Bigg\}\\ &\leq\Big(c+2C(D)\sqrt{\lambda_1}\|u_\infty\|+\frac{1}{\upsilon}-\eta(t)-2\nu\lambda_1\Big)F(t)+\alpha(t)e^{ct}-\alpha(t)e^{ct}F(t)\\ &\quad+e^{ct}\exp\Big(-\int_{0}^{t}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\upsilon C_{g}^{2}\mathbf E|X(t-\rho(t))-u_\infty|^2\\ &\quad+e^{ct}\exp\Big(-\int_{0}^{t}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\big[\beta_1+\gamma_1(t)\big]\mathbf E|X(t-\delta(t))-u_\infty|^2 \\&\quad+e^{ct}\exp\Big(-\int_{0}^{t}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\big[\beta_2+\gamma_2(t)\big]\mathbf E|X(t-\gamma(t))-u_\infty|^2. \end{align*} In what follows, we claim that for any \(t\geq 0\) If the inequality (12) does not hold, then there exists \(t^{\ast}>0\) such that, for any \(\varepsilon>0\) This, in addition to (8), it can be shown that Furthermore, \begin{align*} \frac{dF(t^{\ast})}{dt} &\leq\Big(c+2C(D)\sqrt{\lambda_1}\|u_\infty\|+\frac{1}{\upsilon}-\eta(t^{\ast})-2\nu\lambda_1\Big)F(t^{\ast}) +\alpha(t^{\ast})e^{ct^{\ast}}-\alpha(t^{\ast})e^{ct^{\ast}}F(t^{\ast})\\ &\quad+e^{ct^{\ast}}\exp\Big(-\int_{0}^{t^{\ast}}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\upsilon C_{g}^{2}\mathbf E|X(t^{\ast}-\rho(t^{\ast}))-u_\infty|^2\\ &\quad+e^{ct^{\ast}}\exp\Big(-\int_{0}^{t^{\ast}}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\big[\beta_1+\gamma_1(t^{\ast})\big]\mathbf E|X(t^{\ast}-\delta(t^{\ast}))-u_\infty|^2 \\&\quad+e^{ct^{\ast}}\exp\Big(-\int_{0}^{t^{\ast}}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\big[\beta_2+\gamma_2(t^{\ast})\big]\mathbf E|X(t^{\ast}-\gamma(t^{\ast}))-u_\infty|^2. \end{align*} Next, we split the following cases to derive the desired assertion.

• Case 1: If \(t^{\ast}-\rho(t^{\ast})\geq 0, t^{\ast}-\delta(t^{\ast})\geq 0, t^{\ast}-\gamma(t^{\ast})\geq 0\), we then have from (6) and (13) that \begin{align*} \frac{dF(t^{\ast})}{dt} &\leq\Big(c+2C(D)\sqrt{\lambda_1}\|u_\infty\|+\frac{1}{\upsilon}-\eta(t^{\ast})-2\nu\lambda_1\Big)F(t^{\ast}) +\alpha(t^{\ast})e^{ct^{\ast}}-\alpha(t^{\ast})e^{ct^{\ast}}F(t^{\ast})\\ &\quad+e^{c\rho(t^{\ast})}\exp\Big(-\int_{t^\ast-\rho(t^{\ast})}^{t^{\ast}}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\upsilon C_{g}^{2} F\big(X(t^{\ast}-\rho(t^{\ast}))\big)\\ &\quad+e^{c\delta(t^{\ast})}\exp\Big(-\int_{t^\ast-\delta(t^{\ast})}^{t^{\ast}}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\big[\beta_1+\gamma_1(t^{\ast})\big] F\big(X(t^{\ast}-\delta(t^{\ast}))\big) \\&\quad+e^{c\gamma(t^{\ast})}\exp\Big(-\int_{t^\ast-\gamma(t^{\ast})}^{t^{\ast}}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\big[\beta_2+\gamma_2(t^{\ast})\big] F\big(X(t^{\ast}-\gamma(t^{\ast}))\big) \end{align*} \begin{align*} &\leq\Big(c+2C(D)\sqrt{\lambda_1}\|u_\infty\|+\beta_1+\beta_2+\frac{1}{\upsilon}+\upsilon C_{g}^{2}e^{cr}-2\nu\lambda_1\Big)\widetilde{M}< 0,\end{align*} which contradicts with (14). That is, the desired assertion (12) must hold.
• Case 2: If \(t^{\ast}-\rho(t^{\ast})\leq 0, t^{\ast}-\delta(t^{\ast})\leq 0, t^{\ast}-\gamma(t^{\ast})\leq 0\), we then have from (6) and (13) that \begin{align*} \frac{dF(t^{\ast})}{dt} &\leq\Big(c+2C(D)\sqrt{\lambda_1}\|u_\infty\|+\frac{1}{\upsilon}-\eta(t^{\ast})-2\nu\lambda_1\Big)F(t^{\ast}) +\alpha(t^{\ast})e^{ct^{\ast}}-\alpha(t^{\ast})e^{ct^{\ast}}F(t^{\ast})\\ &\quad+e^{c\rho(t^{\ast})}\exp\Big(-\int_{0}^{t^{\ast}}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\upsilon C_{g}^{2} F\big(X(t^{\ast}-\rho(t^{\ast}))\big)\\ &\quad+e^{c\delta(t^{\ast})}\exp\Big(-\int_{0}^{t^{\ast}}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\big[\beta_1+\gamma_1(t^{\ast})\big] F\big(X(t^{\ast}-\delta(t^{\ast}))\big) \\&\quad+e^{c\gamma(t^{\ast})}\exp\Big(-\int_{0}^{t^{\ast}}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\big[\beta_2+\gamma_2(t^{\ast})\big] F\big(X(t^{\ast}-\gamma(t^{\ast}))\big)\\ &\leq\Big(c+2C(D)\sqrt{\lambda_1}\|u_\infty\|+\beta_1+\beta_2+\frac{1}{\upsilon}+\upsilon C_{g}^{2}e^{cr}-2\nu\lambda_1\Big)\widetilde{M}< 0. \end{align*}

This is a contradiction. Hence, (12) holds true for any \(t\geq 0\). For other cases, for example: \(t^{\ast}-\rho(t^{\ast})\geq 0, t^{\ast}-\delta(t^{\ast})\leq 0, t^{\ast}-\gamma(t^{\ast})\leq 0\),…, in the same way as Cases 1 and 2 were done, we can show (12).

Therefore, from (12), we infer that

\[\mathbf E|X(t)-u_\infty|^2\leq \widetilde{M}e^{-ct}\exp\Big(\int_{0}^{t}[\eta(s)+\alpha(s)e^{cs}]ds\Big)\stackrel{(3.4)}{\leq}\widetilde{M}e^{\Lambda_1+\Lambda_3}e^{-ct},\quad \forall t\geq 0.\] This completes the proof of the theorem.

In the following, \(C\) will denote a generic constant whose values might change from line to line.

In order to prove the almost surely exponential stability of the weak solution of (2) we shall establish the following lemma.

Lemma 3. For any \(t\geq 0\), there exists a constant \(C>0\) such that \begin{align*} &\mathbf E\sup_{\tau\in[0,t]}\Big|\int_{0}^{\tau}\int_{Z}\Big(\big|k\big(s,X(s-\gamma(s)),z\big)\big|_{\mathbb H}^{2}+2\big\langle X(s-),k\big(s,X(s-\gamma(s)),z\big)\big\rangle_{\mathbb H}\Big)\widetilde{\eta}(ds,dz)\Big|\\ &\leq C\mathbf E\int_{0}^{t}\int_{Z}\big|k\big(\tau,X(\tau-\gamma(\tau)),z\big)\big|_{\mathbb H}^{2}\lambda(dz)d\tau+\frac{1}{4}\mathbf E\sup_{\tau\in[0,t]}|X(\tau)|_{\mathbb H}^{2}. \end{align*}

Proof. Set \[J_t:=\int_{0}^{\tau}\int_{Z}\Big(\|\gamma\big(s,X(s-\rho(s)),z\big)\|_{\mathbb H}^{2}+2\big\langle X(s-),\gamma\big(s,X(s-\rho(s)),z\big)\big\rangle_{\mathbb H}\Big)\widetilde{\eta}(ds,dz).\] Then, \begin{align*} [J,J]_{t}^{\frac{1}{2}} &=\Big\{\sum_{s\in D_p,s\leq t}\Big(|\gamma\big(s,X(s-\rho(s)),p(s)\big)|_{\mathbb H}^{2}+2\big\langle X(s-),\gamma\big(s,X(s-\rho(s)),p(s)\big)\big\rangle_{\mathbb H}\Big)^2\Big\}^{\frac{1}{2}}\\ &\leq C\Big(\sum_{s\in D_p,s\leq t}\big|\gamma\big(s,X(s-\rho(s)),p(s)\big)\big|_{\mathbb H}^{4}\Big)^{\frac{1}{2}} +C\Big(\sum_{s\in D_p,s\leq t}\big|X(s-)\big|_{\mathbb H}^{2}\big|\gamma\big(s,X(s-\rho(s)),p(s)\big)\big|_{\mathbb H}^{2}\Big)^{\frac{1}{2}}\\ &\leq C\sum_{s\in D_p,s\leq t}\big|\gamma\big(s,X(s-\rho(s)),p(s)\big)\big|_{\mathbb H}^{2} +C\sup_{0\leq s\leq t}(|X(s-)|_{\mathbb H})\Big(\sum_{s\in D_p,s\leq t}\big|\gamma\big(s,X(s-\rho(s)),p(s)\big)\big|_{\mathbb H}^{2}\Big)^{\frac{1}{2}}\\ &\leq C\sum_{s\in D_p,s\leq t}\big|\gamma\big(s,X(s-\rho(s)),p(s)\big)\big|_{\mathbb H}^{2}+\frac{1}{4}\sup_{0\leq s\leq t}(|X(s-)|_{\mathbb H}^{2}). \end{align*} By Burkholder-Davis-Gundy inequality [43], we obtain that \begin{align*} \mathbf E[\sup_{\tau\in[0,t]}|J_\tau|]&\leq C\mathbf E([J,J]_{t}^{\frac{1}{2}})\\ &\leq C\mathbf E\Big(\sum_{s\in D_p,s\leq t}\big|\gamma\big(s,X(s-\rho(s)),p(s)\big)\big|_{\mathbb H}^{2}\Big)+ \frac{1}{4}\mathbf E\Big(\sup_{0\leq s\leq t}(|X(s-)|_{\mathbb H}^{2})\Big)\\ &=C\mathbf E\int_{0}^{t}\int_{Z}\big|k\big(\tau,X(\tau-\gamma(\tau)),z\big)\big|_{\mathbb H}^{2}\lambda(dz)d\tau+\frac{1}{4}\mathbf E\sup_{\tau\in[0,t]}|X(\tau)|_{\mathbb H}^{2}. \end{align*} The proof is therefore complete.

We have the following theorem:

Theorem 3. Assume that all the assumptions of Theorem 2 are satisfied. Then, any weak solution \(X(t)\) to (2) converges to the stationary solution \(u_\infty\) of (4) almost surely exponentially.

Proof. Let \(n_1,n_2\) and \(n_3\) be positive integers such that \[n_1-\rho(n_1)\geq n_1-r\geq 1,\quad n_2-\delta(n_2)\geq n_2-r\geq 1,\quad n_3-\gamma(n_3)\geq n_3-r\geq 1.\] Set \(n=\max\{n_1,n_2,n_3\}\). By the Itô formula, it follows for any \(t\geq n,\) \begin{align*} |X(t)-u_\infty|^2 &=|X(n)-u_\infty|^2-2\int_{n}^{t}\langle\nu A(X(s)-u_\infty),X(s)-u_\infty\rangle ds-2\int_{n}^{t}\langle B(X(s))-B(u_\infty),X(s)-u_\infty\rangle ds\\ &\;\;\;+2\int_{n}^{t}\langle g(X(s-\rho(s)))-g(u_\infty), X(s)-u_\infty\rangle ds+\int_{n}^{t}\|\sigma(s,X(s-\delta(s)))\|_{\mathcal{L}_{2}^{0}}^{2}ds\\ &\;\;\;+2\int_{n}^{t}\langle X(s)-u_\infty,\sigma(s,X(s-\delta(s))) dW(s)\rangle +\int_{n}^{t}\int_{Z}|k(s,X(s-\gamma(s)),z)|^2\lambda(dz)ds\\ &\;\;\;+\int_{n}^{t}\int_{Z}\big[|k(s,X(s-\gamma(s)),z)|^2+2\int_{n}^{t}\langle X(s)-u_\infty,k(s,X(s-\gamma(s)),z)\rangle \big] \widetilde{\eta}(ds,dz). \end{align*} In view of Burkholder-Davis-Gundy inequality and the Young inequality, we have

Applying Lemma 2, for any \(t\geq 0\) we can get \begin{align*} &\mathbf E\sup_{t\in[n,n+1]}\int_{n}^{t}\int_{Z}\Big(\big|k\big(s,X(s-\gamma(s)),z\big)\big|^{2}+2\big\langle X(s-u_\infty),k\big(s,X(s-\gamma(s)),z\big)\big\rangle\Big)\widetilde{\eta}(ds,dz)\\ &\leq C\mathbf E\int_{n}^{n+1}\int_{Z}\big|k\big(s,X(s-\gamma(s)),z\big)\big|^{2}\lambda(dz)ds+\frac{1}{4}\mathbf E\Big(\sup_{t\in[n,n+1]}|X(t)-u_\infty|^{2}\Big). \end{align*} Hence, \begin{align*} \mathbf E\Big(\sup_{t\in[n,n+1]}|X(t)-u_\infty|^{2}\Big) &\leq \mathbf E|X(n)-u_\infty|^2+\Big(2C(D)\sqrt{\lambda_1}\|u_\infty\|+\frac{1}{\upsilon}-2\nu\lambda_1\Big)\int_{n}^{n+1}\mathbf E|X(s)-u_\infty|^2ds\\ &\quad+\upsilon C_{g}^{2}\int_{n}^{n+1}\mathbf E|X(s-\rho(s))-u_\infty|^2ds+33\int_{n}^{n+1}\mathbf E\|\sigma(s,X(s-\delta(s)))\|_{\mathcal{L}_{2}^{0}}^{2}ds\\ &\quad+(C+1)\mathbf E\int_{n}^{n+1}\int_{Z}\big|k\big(s,X(s-\gamma(s)),z\big)\big|^{2}\lambda(dz)ds+\frac{3}{4}\mathbf E\Big(\sup_{t\in[n,n+1]}|X(t)-u_\infty|^{2}\Big). \end{align*} This implies that \begin{align*} \mathbf E\Big(\sup_{t\in[n,n+1]}|X(t)-u_\infty|^{2}\Big) &\leq 4\mathbf E|X(n)-u_\infty|^2+4\Big(2C(D)\sqrt{\lambda_1}\|u_\infty\|+\frac{1}{\upsilon}-2\nu\lambda_1\Big)\int_{n}^{n+1}\mathbf E|X(s)-u_\infty|^2ds\\ &\quad+4\upsilon C_{g}^{2}\int_{n}^{n+1}\mathbf E|X(s-\rho(s))-u_\infty|^2ds\\ &\quad+132\int_{n}^{n+1}\big[\alpha_1(s)+(\beta_1+\gamma_1(s))\mathbf E|X(s-\delta(s))-u_\infty|^2\big]ds\\ &\quad+4(C+1)\int_{n}^{n+1}\big[\alpha_2(s)+(\beta_2+\gamma_2(s))\mathbf E|X(s-\gamma(s))-u_\infty|^2\big]ds \end{align*} On the other hand, from Theorem 2, it is easy to show that \begin{align*} \mathbf E\Big(\sup_{t\in[n,n+1]}|X(t)-u_\infty|^{2}\Big) \leq 4M_0e^{-cn}+4\int_{n}^{n+1}M_0e^{-cs}[\gamma^\ast(s)+\alpha^\ast(s)e^{cs}]ds, \end{align*} where \begin{eqnarray*} \alpha^\ast(t):=33\alpha_1(t)+4(C+1)\alpha_2(t), \end{eqnarray*} \begin{eqnarray*} \gamma^\ast(t):=2C(D)\sqrt{\lambda_1}\|u_\infty\|+\frac{1}{\upsilon}-2\nu\lambda_1+\big[\upsilon C_{g}^{2}+33(\beta_1+\gamma_1(t))+(C+1)(\beta_2+\gamma_2(t))\big]e^{cr}. \end{eqnarray*} In view of assumption \((\mathbf {H4})\), there exists a positive constant \(\Lambda\) such that we obtain that \[\mathbf E\Big(\sup_{t\in[n,n+1]}|X(t)-u_\infty|^{2}\Big)\leq 4M_0e^{-cn}\Big(1+\frac{\Lambda}{c}\Big).\] Let \(\epsilon_n>0\) be any fixed positive real number. Then by Chebychev’s inequality, we deduce that \[\mathbf P\Big\{\sup_{t\in[n,n+1]}|X(t)-u_\infty|>\epsilon_n\Big\}\leq \frac{4M_0e^{-cn}\Big(1+\frac{\Lambda}{c}\Big)}{\epsilon_{n}^{2}}.\] Therefore, since \(\epsilon_n\) is any fixed real number, let \(\epsilon_n=e^{-\frac{(c-\varepsilon) n}{4}}\), where \(\varepsilon\in (0,c)\). Then by the Borel-Cantelli lemma [44], we can yield that \[\overline{\lim}_{t\rightarrow\infty}\frac{\log|X(t)-u_\infty|}{t}\leq -\frac{c-\varepsilon}{4},\quad a.s..\] Letting \(\varepsilon\rightarrow 0^+\), this completes the proof of the theorem.

Remark 1. We consider a special case of the system (1) with variable delays when \(\rho\equiv\delta\) and \(\gamma\equiv 0\). That is, our system (1) reduces to the following system

The system (16) has been recently studied by Chen [22]. Furthermore, Wan and Zhou [23] discussed the following system \begin{equation*} \begin{cases} dX(t)=\big[-\nu AX(t)-B(X(t))+f(X(t))+g(X(\rho(t)))\big]dt +\sigma(t,X(t))dW(t), \\ \text{div \(X=0\)} \quad in \quad (0,+\infty)\times D,\quad X(t,x)=0 \quad on \quad(0,+\infty)\times \partial D,\quad \\ X(t,x)=\phi(t,x), \text{ for \(x\in D\), and \(t\in [-r,0]\), with \(r>0\)}, \end{cases} \end{equation*} and Caraballo et al., [7] also studied the stability of the stationary solutions of the following stochastic 2D Navier-Stokes without memory: \begin{equation*} \begin{cases} dX(t)=\big[-\nu AX(t)-B(X(t))+g(X(t))\big]dt +\sigma(t,X(t))dW(t), \\ \text{div \(X=0\)} \quad in \quad (0,+\infty)\times D,\quad X(t,x)=0 \quad on \quad(0,+\infty)\times \partial D,\\ X(0,x)=X_{0}(x),\quad x\in D. \end{cases} \end{equation*} By using method in our paper, the conclusions of some theorems in the works [7,22,23] also easily obtained. Obviously, our work have extended the asymptotic behavior results of above works to cover a class of much more general stochastic 2D Navier-Stokes equations with memory and discontinuous multiplicative noise.

Remark 2. If \(\sigma\equiv 0\) in the system (16), then by utilizing the direct method, in [20], Caraballo and Real have considered the asymptotic behavior for the weak solutions and Taniguchi [35] have investigated exponential stability of energy solutions to 2D stochastic functional Navier-Stokes equation perturbed by the Lévy process. However, unlike the works [20,35], we need not require the function \(\rho(t)\) to be differentiable and satisfies \(0\leq \rho'(t)< 1\). Therefore, our results extend and improve the one of Caraballo and Real [20] and Taniguchi [35].

Conflict of Interests

The author declares no conflict of interest.