Well-Posedness and Global Attractors for Viscous Fractional Cahn–Hilliard Equations with Memory

Abstract

We examine a viscous Cahn–Hilliard phase-separation model with memory and where the chemical potential possesses a nonlocal fractional Laplacian operator. The existence of global weak solutions is proven using a Galerkin approximation scheme. A continuous dependence estimate provides uniqueness of the weak solutions and also serves to define a precompact pseudometric. This, in addition to the existence of a bounded absorbing set, shows that the associated semigroup of solution operators admits a compact connected global attractor in the weak energy phase space. The minimal assumptions on the nonlinear potential allow for arbitrary polynomial growth.

1. Introduction

Let $\mathsf{\Omega }$ be a smooth (at least Lipschitz) bounded domain in ${\mathbb{R}}^N$, $N=3,2,1$, with boundary $\partial \mathsf{\Omega }$ and let $T>0$. We consider the following viscous fractional Cahn–Hilliard equation in the unknown (order parameter) u satisfying

$$\begin{array}{c}\hfill {\partial }_{t}u(t,x)={\int }_0^{\infty }k\left(s\right)\mathsf{\Delta }\mu (t-s,x)ds\mathrm{in}\mathsf{\Omega }\times (0,T),\end{array}$$

(1)

k is a so-called relaxation kernel, with a chemical potential $\mu$ given by

$$\begin{array}{c}\hfill \mu (t,x)=\alpha {\partial }_{t}u(t,x)+{(-\mathsf{\Delta })}^{\beta }u(t,x)+F^{\prime }\left(u(t,x)\right)\mathrm{in}\mathsf{\Omega }\times \mathbb{R},\end{array}$$

(2)

$\alpha >0,$ $\beta \in (0,1)$, and typically, F is a double-well potential (the precise assumptions on F are stated in (N1)–(N3) below), subject to the boundary conditions

$$\begin{array}{c}\hfill u=0\mathrm{on}{\mathbb{R}}^N\backslash \mathsf{\Omega }\times (0,T)\mathrm{and}{\partial }_{\mathbf{n}}\mu =0\mathrm{on}\partial \mathsf{\Omega }\times (0,T),\end{array}$$

(3)

with the given initial and past conditions

$$\begin{array}{c}\hfill u\left(0\right)=u_0\left(0\right)\mathrm{in}\mathsf{\Omega }\mathrm{and}u(-t)=u_0(-t)\mathrm{in}\mathsf{\Omega }\times [0,T),\end{array}$$

(4)

for

$$u_0:\mathsf{\Omega }\times (-\infty ,0)\to \mathbb{R}.$$

Here, we define ${(-\mathsf{\Delta })}^{\beta }$ with $0<\beta <1$ as the (nonlocal) fractional Laplace operator. In other words, let $\mathsf{\Omega }\subset {\mathbb{R}}^N$ be an arbitrary open set and fix

$${\mathcal{L}}^1\left(\mathsf{\Omega }\right):=\left\{u:\mathsf{\Omega }\to \mathbb{R}\mathrm{measurable},{\int }_{\mathsf{\Omega }}\frac{\left|u\right(x\left)\right|}{{(1+|x\left|\right)}^{N+2\beta }}dx<\infty \right\}.$$

For $u\in {\mathcal{L}}^1\left({\mathbb{R}}^N\right)$, $x\in {\mathbb{R}}^N$, and $\epsilon >0$, we write

$${(-\mathsf{\Delta })}_{\epsilon }^{\beta }u\left(x\right)=C_{N,\beta }{\int }_{\{y\in {\mathbb{R}}^N,|y-x|>\epsilon \}}\frac{u\left(x\right)-u\left(y\right)}{{|x-y|}^{N+2\beta }}dy$$

with the normalized constant $C_{N,\beta }$ given by

$$\begin{array}{c}\hfill C_{N,\beta }=\frac{\beta 2^{2\beta }\mathsf{\Gamma }\left(\frac{N+2\beta }{2}\right)}{{\pi }^{\frac{N}{2}}\mathsf{\Gamma }(1-\beta )},\end{array}$$

(5)

where $\mathsf{\Gamma }$ denotes the usual gamma function. The (restricted) fractional Laplacian ${(-\mathsf{\Delta })}^{\beta }u$ of the function u is defined by the formula

$$\begin{array}{c}\hfill {(-\mathsf{\Delta })}^{\beta }u\left(x\right)=C_{N,\beta }\mathrm{P}.\mathrm{V}.{\int }_{{\mathbb{R}}^N}\frac{u\left(x\right)-u\left(y\right)}{{|x-y|}^{N+2\beta }}dy=\underset{\epsilon \downarrow 0}{lim}{(-\mathsf{\Delta })}_{\epsilon }^{\beta }u\left(x\right),x\in {\mathbb{R}}^N,\end{array}$$

(6)

provided that the limit exists. We call $A^{\beta }$ the self-adjoint realization of the fractional Laplacian ${(-\mathsf{\Delta })}^{\beta }$ with Dirichlet boundary condition (3)${}_1$, see, e.g., [1] (Section 2.2) (see also [2]).

Some remarks: First, observe the chemical potential (2) involves the Neumann (no-flux) condition described by (3). Hence, when the memory function k is close to the Dirac delta function, we recover the usual parabolic equation associated with the Cahn–Hilliard equation with the flux-free chemical potential.

Naturally, we are also interested in the closely related problem to (1)–(4) whereby the fractional Laplace operator ${(-\mathsf{\Delta })}^{\beta }$ is replaced with the regional fractional Laplacian, $A_{\mathsf{\Omega }}^{\beta }$, defined by first setting

$$A_{\mathsf{\Omega },\epsilon }^{\beta }u\left(x\right)=C_{N,\beta }{\int }_{\{y\in \mathsf{\Omega },|y-x|>\epsilon \}}\frac{u\left(x\right)-u\left(y\right)}{{|x-y|}^{N+2\beta }}dy,$$

where $C_{N,\beta }$ is given by (5), then

$$\begin{array}{c}\hfill A_{\mathsf{\Omega }}^{\beta }u\left(x\right)=C_{N,\beta }\mathrm{P}.\mathrm{V}.{\int }_{\mathsf{\Omega }}\frac{u\left(x\right)-u\left(y\right)}{{|x-y|}^{N+2\beta }}dy=\underset{\epsilon \downarrow 0}{lim}{A}_{\mathsf{\Omega },\epsilon }^{\beta }u\left(x\right),x\in \mathsf{\Omega },\end{array}$$

(7)

provided that the limit exists. Assuming $u\in \mathcal{D}\left(\mathsf{\Omega }\right)$ (see [1] (page 1280)) then the two fractional Laplacian operators are related by

$$\begin{array}{c}\hfill {(-\mathsf{\Delta })}^{\beta }u\left(x\right)=A_{\mathsf{\Omega }}^{\beta }u\left(x\right)+V_{\mathsf{\Omega }}\left(x\right)u\left(x\right),\forall u\in \mathcal{D}\left(\mathsf{\Omega }\right)\end{array}$$

(8)

with the following potential

$$\begin{array}{c}\hfill V_{\mathsf{\Omega }}\left(x\right):=C_{N,\beta }{\int }_{{\mathbb{R}}^N\backslash \mathsf{\Omega }}\frac{dy}{{|x-y|}^{N+2\beta }},x\in \mathsf{\Omega }.\end{array}$$

(9)

The comparable Cahn–Hilliard problem with the regional fractional Laplacian is then (1) with the chemical potential

$$\begin{array}{c}\hfill \mu =\alpha {\partial }_{t}u+A_{\mathsf{\Omega }}^{\beta }u+F^{\prime }\left(u\right)\mathrm{in}\mathsf{\Omega }\times (0,T),\end{array}$$

(10)

now subject to the boundary conditions

$$\begin{array}{c}\hfill u=0\mathrm{on}\partial \mathsf{\Omega }\times (0,T)\mathrm{and}{\partial }_{\mathbf{n}}\mu =0\mathrm{on}\partial \mathsf{\Omega }\times (0,T),\end{array}$$

(11)

with the above initial and past conditions in (4). Our focus here is on obtaining results for the restricted fractional Laplacian, of which the regional counterpart can be view as a perturbation thanks to (8). The restricted fractional Laplacian appears in the context of nonlocal phase transitions with Dirichlet boundary conditions in [3,4]. On the other hand, the regional fractional Laplacian is generally better suited to treat problems with nonhomogeneous boundary data and even dynamic boundary conditions (see [1,5] and the references therein).

It should also be noted that we only consider the viscous case, where $\alpha >0$, since the nonviscous counterpart $\alpha =0$ inherits no added regularity for $\alpha {\partial }_t\varphi$.

Inside a bounded container $\mathsf{\Omega }\subset {\mathbb{R}}^3,$ the Cahn–Hilliard equation (see [6]) is a phase separation model for a binary solution (e.g., a cooling alloy, glass, or polymer),

$${\partial }_{t}u=\nabla ·[\kappa \left(u\right)\nabla \mu ],$$

where u is the order-parameter (the relative difference of the two phases), $\kappa$ is the mobility function (we set $\kappa \equiv 1$ throughout this article), and $\mu$ is the chemical potential (the first variation of the free-energy E with respect to u). In the classical model,

$$\mu =-\mathsf{\Delta }u+F^{\prime }\left(u\right)\mathrm{and}E\left(u\right)={\int }_{\mathsf{\Omega }}\left(\frac{1}{2}{|\nabla u|}^2+F\left(u\right)\right)dx,$$

where F describes the density of potential energy in $\mathsf{\Omega }$ (e.g., the double-well potential $F\left(s\right)={(1-s^2)}^2$).

Recently the nonlocal free-energy functional has appeared in the literature [7],

$$E\left(\varphi \right)={\int }_{\mathsf{\Omega }}{\int }_{\mathsf{\Omega }}\frac{1}{4}J(x-y){(\varphi \left(x\right)-\varphi \left(y\right))}^{2}dxdy+{\int }_{\mathsf{\Omega }}F\left(\varphi \right)dx,$$

hence, the chemical potential is,

$$\begin{array}{c}\hfill \mu =a\varphi -J\ast \varphi +F^{\prime }\left(\varphi \right),\end{array}$$

(12)

where

$$\begin{array}{c}\hfill a\left(x\right)={\int }_{\mathsf{\Omega }}J(x-y)dy\mathrm{and}(J\ast \varphi )\left(x\right)={\int }_{\mathsf{\Omega }}J(x-y)\varphi \left(y\right)dy.\end{array}$$

(13)

In view of [8,9], the nonlocality expressed in (12)–(13) (see also [10,11,12,13,14,15,16,17,18,19]) is termed weak while the type under consideration here in (2) and (6) is called strong. Under certain conditions the strong type reduces to the weak (see [8], and also see [7]). Recently there has been much interest in the nonlocal Cahn–Hilliard equation with strong interactions of the restricted fractional Laplacian type (6) and the regional fractional Laplacian type (7) (see [3,5,8,9,20]). The results in these references concern global well-posedness, and when available, the existence of finite dimensional global attractors and regularity.

Additionally, there has been exceptional growth concerning dissipative infinite-dimensional systems with memory including models arising in the theory of heat conduction in special materials (see, e.g., [21,22,23,24,25]) and the theory of phase-transitions (see, e.g., [26,27,28,29,30,31,32,33,34]). One feature of equations that undergo “memory relaxation” is the admissibility of a so-called inertia term. For example, (see, e.g., [35]) the first-order equation with memory

$$u_t\left(t\right)+{\int }_0^{\infty }{k}_{ϵ}\left(s\right)f\left(u(t-s)\right)ds=0$$

for

$$k_{ϵ}(s)=\frac{1}{ϵ}{e}^{-s/ϵ}$$

leads us (formally) to the “hyperbolic relaxation” equation

$$ϵu_{\mathrm{t}\mathrm{t}}(t)+u_t(t)+f(u(t))=0.$$

In this way, our model also includes the viscous Cahn–Hilliard equation with inertial term (see [36]). Hence, the novelty in the present work is a relaxation of a phase-field model with a strongly interacting nonlocal diffusion mechanism.

In this article, our aims were:

• To provide a framework to establish the global (in time) well-posedness of the model problems (1)–(4) and (1), (4), (10), and (11).
• To prove the semigroup of solution operators admits a compact global attractor.

In order to reach these aims, we require sufficient growth conditions on F (given below) in order to employ a Galerkin scheme with suitable a priori estimates. With a finite energy phase space identified, a one-parameter family of solution operators is defined, hence generating a semidynamical system. This semigroup is dissipative on the energy phase space and also defines an $\alpha$-contraction on the phase space. The existence of a compact global attractor follows.

2. Past History Formulation and Functional Setup

We now introduce the well-established past history approach from [37] (see also [27,29]) by defining the past history variable, for all $s>0$ and $t>0,$

$$\begin{array}{c}\hfill {\eta }^t(x,s)={\int }_0^s-\mathsf{\Delta }\mu (x,t-\sigma )d\sigma .\end{array}$$

(14)

Observe that $\eta$ satisfies the boundary condition

$$\begin{array}{c}\hfill {\eta }^t(x,0)=0\mathrm{on}\mathsf{\Omega }\times (0,\infty ).\end{array}$$

(15)

When k is sufficiently smooth and vanishes at $+\infty$ (these assumptions will be made more precise below), then integration by parts yields

$${\int }_0^{\infty }k\left(s\right)\mathsf{\Delta }\mu (x,t-s)ds=-{\int }_0^{\infty }\nu \left(s\right){\eta }^t(x,s)ds$$

where $\nu \left(s\right)=-k^{\prime }\left(s\right)$.

We may now formulate the model problem (1)–(4) as:

Problem P. Find $(u,\eta )=(u(x,t),{\eta }^t(x,s))$ on $(0,\infty )$ such that

$$\begin{array}{ccc}\hfill {\partial }_{t}u(x,t)+{\int }_0^{\infty }\nu \left(s\right){\eta }^t(x,s)ds=0& \mathrm{in}& \mathsf{\Omega }\times (0,\infty )\hfill \end{array}$$

(16)

$$\begin{array}{ccc}\hfill \mu (x,t)=\alpha {\partial }_{t}u(x,t)+{(-\mathsf{\Delta })}^{\beta }u(x,t)+F^{\prime }\left(u(x,t)\right)& \mathrm{in}& \mathsf{\Omega }\times (0,\infty )\hfill \end{array}$$

(17)

$$\begin{array}{ccc}\hfill {\partial }_t{\eta }^t(x,s)+{\partial }_s{\eta }^t(x,s)=-\mathsf{\Delta }\mu (x,t)& \mathrm{in}& \mathsf{\Omega }\times (0,\infty )\times (0,\infty )\hfill \end{array}$$

(18)

held subject to (3) and (15), and satisfying the initial conditions (4)${}_1$ and

$$\begin{array}{c}\hfill {\eta }^0(x,s)={\eta }_0(x,s)\mathrm{in}\mathsf{\Omega }\times (0,\infty ),\end{array}$$

(19)

whereby with (14),

$$\begin{array}{c}\hfill {\eta }_0(x,s)={\int }_0^s-\mathsf{\Delta }{\mu }_0(x,-y)dy\mathrm{in}\mathsf{\Omega }\times (0,\infty ),\end{array}$$

(20)

where in light of (4)${}_2$,

$$\begin{array}{c}\hfill {\mu }_0(x,t)=\alpha {\partial }_t{u}_0(x,t)+{(-\mathsf{\Delta })}^{\beta }{u}_0(x,t)+F^{\prime }\left(u_0(x,t)\right)\mathrm{for}t\le 0.\end{array}$$

(21)

Additionally, we are also interested in treating the related problem where the above fractional Laplace operator ${(-\mathsf{\Delta })}^{\beta }$ is replaced with the regional counterpart $A_{\mathsf{\Omega }}^{\beta }$. Hence, the formulation of the related regional Problem P is based on (1), (4), (10), and (11).

Here, we introduce some notation. From now on, we denote by ${\parallel ·\parallel }_X$ the norm in the specified (real) Banach space X, and ${(·,·)}_Y$ denotes the product on the specified (real) Hilbert space Y. The dual pairing between Y and the dual $Y^{*}$ is denoted by ${\langle u,v\rangle }_{Y^{*}\times Y}$. The set $\mathsf{\Omega }$ is omitted from the space when we indicate the norm. We denote the measure of the domain $\mathsf{\Omega }$ by $\left|\mathsf{\Omega }\right|$. In many calculations, functional notation indicating dependence on the variable t is dropped; for example, we write u in place of $u\left(t\right)$ or ${\eta }^t$ in place of ${\eta }^t\left(s\right)$. Throughout the paper, C denotes a generic positive constant, while $Q:{\mathbb{R}}_{+}\to {\mathbb{R}}_{+}$ denotes a generic increasing function. Such generic terms may or may not indicate dependencies on the (physical) parameters of the model problem, and may even change from line to line.

Let us define the linear operator $A_N:=-\mathsf{\Delta }$ on $D\left(A_N\right)=\{\psi \in H^2\left(\mathsf{\Omega }\right):{\partial }_{\mathbf{n}}\psi =0\mathrm{on}\partial \mathsf{\Omega }\}$, as the realization in $L^2\left(\mathsf{\Omega }\right)$ of the Laplace operator endowed with Neumann boundary conditions. Here, $-\mathsf{\Delta }$ denotes the usual (local) Laplace operator. It is well-known that $A_N$ is the generator of a bounded analytic semigroup $e^{-A_{N}t}$ on $L^2\left(\mathsf{\Omega }\right)$. Additionally, $A_N$ is nonnegative and self-adjoint on $L^2\left(\mathsf{\Omega }\right).$ With $H^{-r}\left(\mathsf{\Omega }\right):={\left(H^r\left(\mathsf{\Omega }\right)\right)}^{*}$, $r\in {\mathbb{N}}_{+}$, denote by $\langle ·\rangle$ the spatial average over $\mathsf{\Omega }$, i.e.,

$$\langle \psi \rangle :=\frac{1}{\left|\mathsf{\Omega }\right|}{\langle \psi ,1\rangle }_{H^{-r}\times H^r}.$$

We set $H_{\left(0\right)}^r\left(\mathsf{\Omega }\right)=\{\psi \in H^r\left(\mathsf{\Omega }\right):\langle \psi \rangle =0\}$, $H^0\left(\mathsf{\Omega }\right)=L^2\left(\mathsf{\Omega }\right),$ and we know that $A_N^{-1}:H_{\left(0\right)}^0\left(\mathsf{\Omega }\right)\to H_{\left(0\right)}^0\left(\mathsf{\Omega }\right)$ is a well-defined mapping. We refer to the following norms in $H^{-r}\left(\mathsf{\Omega }\right)$ (which are equivalent to the usual norms)

$${\parallel \psi \parallel }_{H^{-r}}^2=\parallel A_N^{-r/2}{(\psi -\langle \psi \rangle )\parallel }^2+{\left|\langle \psi \rangle \right|}^2.$$

(22)

The Sobolev space $H^1\left(\mathsf{\Omega }\right)$ is endowed with the norm,

$${\parallel \psi \parallel }_{H^1}^2:={\parallel \nabla \psi \parallel }^2+{\langle \psi \rangle }^2.$$

(23)

Denote by ${\lambda }_{\mathsf{\Omega }}>0$ the constant in the Poincaré–Wirtinger inequality,

$$\parallel \psi -\langle \psi \rangle \parallel \le \sqrt{{\lambda }_{\mathsf{\Omega }}}\parallel \nabla \psi \parallel .$$

(24)

Whence, for ${\lambda }_{\mathsf{\Omega }}^{*}:=max\{{\lambda }_{\mathsf{\Omega }},1\}$, there holds, for all $\psi \in H^1\left(\mathsf{\Omega }\right),$

$$\begin{array}{cc}\hfill {\parallel \psi \parallel }^2& \le {\lambda }_{\mathsf{\Omega }}{\parallel \nabla \psi \parallel }^2+{\langle \psi \rangle }^2\hfill \\ \hfill & \le {\lambda }_{\mathsf{\Omega }}^{*}{\parallel \psi \parallel }_{H^1}^2.\hfill \end{array}$$

(25)

We now more rigorously describe the fractional Laplacian with Dirichlet boundary conditions. For an arbitrary bounded domain $\mathsf{\Omega }\subset {\mathbb{R}}^N$ and for $\beta \in (0,1)$, denote the fractional-order Sobolev space by,

$$W^{\beta ,2}\left(\mathsf{\Omega }\right):=\left\{u\in L^2\left(\mathsf{\Omega }\right):{\int }_{\mathsf{\Omega }}{\int }_{\mathsf{\Omega }}\frac{{|u\left(x\right)-u\left(y\right)|}^2}{{|x-y|}^{N+2\beta }}dxdy<\infty \right\},$$

to be equipped with the norm

$${\parallel u\parallel }_{W^{\beta ,2}}:={\left({\int }_{\mathsf{\Omega }}{\left|u\left(x\right)\right|}^{2}dx+\frac{C_{N,\beta }}{2}{\int }_{\mathsf{\Omega }}{\int }_{\mathsf{\Omega }}\frac{{|u\left(x\right)-u\left(y\right)|}^2}{{|x-y|}^{N+2\beta }}dxdy\right)}^{1/2},$$

where $C_{N,\beta }$ is given by (5). Let

$$W_0^{\beta ,2}\left(\mathsf{\Omega }\right)={\overline{\mathcal{D}\left(\mathsf{\Omega }\right)}}^{W^{\beta ,2}\left(\mathsf{\Omega }\right)}.$$

Hence, $W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$ is a closed subspace of $W^{\beta ,2}\left(\mathsf{\Omega }\right)$ containing $\mathcal{D}\left(\mathsf{\Omega }\right)$. Moreover, thanks to [38] (Theorem 10.1.1),

$$W_0^{\beta ,2}\left(\mathsf{\Omega }\right)=\{u\in W^{\beta ,2}\left({\mathbb{R}}^N\right):\tilde{u}=0\mathrm{on}{\mathbb{R}}^N\backslash \mathsf{\Omega }\},$$

where $\tilde{u}$ is the quasi-continuous version (with respect to the capacity defined with the space $W^{\beta ,2}\left(\mathsf{\Omega }\right)$) of u. One may easily show that the following defines an equivalent norm on the space $W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$,

$$\begin{array}{cc}\hfill {|\parallel u\parallel |}_{W_0^{\beta ,2}}^2& =\frac{C_{N,\beta }}{2}{\int }_{\mathsf{\Omega }}{\int }_{\mathsf{\Omega }}\frac{{|u\left(x\right)-u\left(y\right)|}^2}{{|x-y|}^{N+2\beta }}dxdy+{\int }_{\mathsf{\Omega }}{V}_{\mathsf{\Omega }}\left(x\right){\left|u\left(x\right)\right|}^{2}dx\hfill \\ \hfill & =\frac{C_{N,\beta }}{2}{\int }_{{\mathbb{R}}^N}{\int }_{{\mathbb{R}}^N}\frac{{|u\left(x\right)-u\left(y\right)|}^2}{{|x-y|}^{N+2\beta }}dxdy.\hfill \end{array}$$

(26)

Here, $V_{\mathsf{\Omega }}$ is the potential (9).

Remark 1.

Either definition of the space $W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$ makes sense for any arbitrary open set $\mathsf{\Omega }\subset {\mathbb{R}}^3$ (not necessarily bounded). Furthermore, if Ω has a Lipschitz boundary, then by [39], $W_0^{\beta ,2}\left(\mathsf{\Omega }\right)=W^{\beta ,2}\left(\mathsf{\Omega }\right)$ for every $0<\beta \le \frac{1}{2}$.

From now on, we write $u\in W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$ to mean $u\in W^{\beta ,2}\left({\mathbb{R}}^N\right)$ and $u=0$ on ${\mathbb{R}}^N\backslash \mathsf{\Omega }.$ Let $a_{E,\beta }$ be the bilinear symmetric closed form with domain $D\left(a_{E,\beta }\right)=W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$ and defined for $u,v\in W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$ by

$$\begin{array}{cc}\hfill a_{E,\beta }(u,v)& =\frac{C_{N,\beta }}{2}{\int }_{\mathsf{\Omega }}{\int }_{\mathsf{\Omega }}\frac{\left(u\right(x)-u(y\left)\right)\left(v\right(x)-v(y\left)\right)}{{|x-y|}^{N+2\beta }}dxdy+{\int }_{\mathsf{\Omega }}{V}_{\mathsf{\Omega }}\left(x\right)u\left(x\right)v\left(x\right)dx\hfill \\ \hfill & =\frac{C_{N,\beta }}{2}{\int }_{{\mathbb{R}}^N}{\int }_{{\mathbb{R}}^N}\frac{\left(u\right(x)-u(y\left)\right)\left(v\right(x)-v(y\left)\right)}{{|x-y|}^{N+2\beta }}dxdy.\hfill \end{array}$$

(27)

Let $A_{E,\beta }$ be the closed linear self-adjoint operator on $L^2\left(\mathsf{\Omega }\right)$ associated with $a_{E,\beta }$ by

$$\begin{array}{c}\hfill \left\{\begin{array}{c}D\left(A_{E,\beta }\right):=\{u\in W_0^{\beta ,2}\left(\mathsf{\Omega }\right):\exists v\in L^2\left(\mathsf{\Omega }\right),a_{E,\beta }(u,\phi )=(v,\phi )\forall \phi \in W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\}\hfill \\ A_{E,\beta }u=v.\hfill \end{array}\right.\end{array}$$

(28)

According to [1] (Proposition 2.2), the operator $A_{E,\beta }$ on $L^2\left(\mathsf{\Omega }\right)$ associated with the bilinear form $a_{E,\beta }$ is given by

$$\begin{array}{c}\hfill D\left(A_{E,\beta }\right):=\{u\in W_0^{\beta ,2}\left(\mathsf{\Omega }\right):{(-\mathsf{\Delta })}_E^{\beta }u\in L^2\left(\mathsf{\Omega }\right)\}\mathrm{and}\forall u\in D\left(A_{E,\beta }\right),A_{E,\beta }u:={(-\mathsf{\Delta })}_E^{\beta }u.\end{array}$$

(29)

Observe that comparing (6) and (26)–(29) shows, for all $u\in D\left(A_{E,\beta }\right)$,

$$\begin{array}{cc}\hfill ({(-\mathsf{\Delta })}_E^{\beta }u,u)& =a_{E,\beta }(u,u)={|\parallel u\parallel |}_{W_0^{\beta ,2}}^2.\hfill \end{array}$$

(30)

Concerning the related regional problem discussed above, we let $a_{D,\beta }$ be the bilinear symmetric closed form with domain $D\left(a_{D,\beta }\right)=W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$ and defined for $u,v\in W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$ by

$$\begin{array}{cc}\hfill a_{D,\beta }(u,v)& =\frac{C_{N,\beta }}{2}{\int }_{\mathsf{\Omega }}{\int }_{\mathsf{\Omega }}\frac{\left(u\right(x)-u(y\left)\right)\left(v\right(x)-v(y\left)\right)}{{|x-y|}^{N+2\beta }}dxdy.\hfill \end{array}$$

(31)

Let $A_{D,\beta }$ be the closed linear self-adjoint operator on $L^2\left(\mathsf{\Omega }\right)$ associated with $a_{D,\beta }$ by

$$\begin{array}{c}\hfill \left\{\begin{array}{c}D\left(A_{D,\beta }\right):=\{u\in W_0^{\beta ,2}\left(\mathsf{\Omega }\right):\exists v\in L^2\left(\mathsf{\Omega }\right),a_{D,\beta }(u,\phi )=(v,\phi )\forall \phi \in W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\}\hfill \\ A_{D,\beta }u=v.\hfill \end{array}\right.\end{array}$$

(32)

Then, by [1] (Proposition 2.3), the operator $A_{D,\beta }$ on $L^2\left(\mathsf{\Omega }\right)$ associated with the bilinear form $a_{D,\beta }$ is given by

$$\begin{array}{c}\hfill D\left(A_{D,\beta }\right):=\{u\in W_0^{\beta ,2}\left(\mathsf{\Omega }\right):A_{\mathsf{\Omega }}^{\beta }u\in L^2\left(\mathsf{\Omega }\right)\}\mathrm{and}\forall u\in D\left(A_{D,\beta }\right),A_{D,\beta }u:=A_{\mathsf{\Omega }}^{\beta }u.\end{array}$$

(33)

We introduce the spaces for the memory variable $\eta .$ First, the product in $H^{\sigma }\left(\mathsf{\Omega }\right)$ for $\sigma \in \mathbb{R}$ and $u_1,u_2\in H^{\sigma }\left(\mathsf{\Omega }\right)$ is defined by

$$\begin{array}{c}\hfill {(u_1,u_2)}_{H^{\sigma }}=(A_N^{\sigma /2}{u}_1,A_N^{\sigma /2}{u}_2).\end{array}$$

(34)

For a nonnegative measurable function $\theta$ defined on ${\mathbb{R}}_{+}$ and for a Hilbert space W (with inner-product ${(·,·)}_W$), let $L_{\theta }^2({\mathbb{R}}_{+};W)$ be the Hilbert space of W-valued functions on ${\mathbb{R}}_{+}$ equipped with the following product,

$${({\varphi }_1,{\varphi }_2)}_{L_{\theta }^2({\mathbb{R}}_{+};W)}={\int }_0^{\infty }\theta \left(s\right){({\varphi }_1\left(s\right),{\varphi }_2\left(s\right))}_{W}ds.$$

Thus, we set

$${\mathcal{M}}_{\sigma }=L_{\nu }^2({\mathbb{R}}_{+};H^{\sigma }\left(\mathsf{\Omega }\right))\mathrm{and}{\mathcal{M}}_{\sigma }^{\left(0\right)}=L_{\nu }^2({\mathbb{R}}_{+};H_{\left(0\right)}^{\sigma }\left(\mathsf{\Omega }\right))\mathrm{for}\sigma \in \mathbb{R},$$

where $\nu =\nu \left(s\right)$ is the kernel from (16). Hence, for $\sigma \in \mathbb{R}$ and ${\varphi }_1,{\varphi }_2\in {\mathcal{M}}_{\sigma }$, using (34) the product in ${\mathcal{M}}_{\sigma }$ (and ${\mathcal{M}}_{\sigma }^{\left(0\right)}$) can be expressed as

$${({\varphi }_1,{\varphi }_2)}_{{\mathcal{M}}_{\sigma }}={\int }_0^{\infty }\nu \left(s\right)(A_N^{\sigma /2}{\varphi }_1\left(s\right),A_N^{\sigma /2}{\varphi }_2\left(s\right))ds.$$

Naturally, we may also consider spaces of the form $H_{\nu }^k({\mathbb{R}}_{+};H^{\sigma }\left(\mathsf{\Omega }\right))$ for $k\in \mathbb{N}$.

We mention that solutions of Problem P must also satisfy the mass conservation constraints,

$$\langle u\left(t\right)\rangle =\langle u_0\left(0\right)\rangle \mathrm{and}\langle {\eta }^t\left(s\right)\rangle =0\forall t>0,\forall s>0.$$

(35)

With this, it is important to realize that the norm of ${\eta }^t$ in the space ${\mathcal{M}}_{-1}^{\left(0\right)}$ may be expressed without writing the average value of ${\eta }_0$ in (22) by virtue of the second constraint of (35). Indeed, for ${\eta }^t\in {\mathcal{M}}_{-1}^{\left(0\right)}$,

$$\begin{array}{cc}\hfill \parallel {\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}& ={\left({\int }_0^{\infty }\nu \left(s\right){\parallel {\eta }^t\left(s\right)\parallel }_{H^{-1}}^{2}ds\right)}^{1/2}\hfill \\ \hfill & ={\left({\int }_0^{\infty }\nu \left(s\right){\parallel A_N^{-1/2}{\eta }^t\left(s\right)\parallel }^{2}ds\right)}^{1/2}.\hfill \end{array}$$

We now state the basic function spaces we intend to study Problem P in. For each $\beta \in (0,1)$ and $\sigma \in \mathbb{R}$, define the following (weak) energy Hilbertian phase-space ${\mathcal{H}}_{\beta ,\sigma }:=W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\times {\mathcal{M}}_{\sigma -1}^{\left(0\right)}$, equipped with the norm on $W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\times {\mathcal{M}}_{\sigma -1}^{\left(0\right)}$ whose square is given by, for all $\varphi ={(u,\eta )}^{tr}\in {\mathcal{H}}_{\beta ,\sigma }$,

$$\begin{array}{cc}\hfill {\parallel \varphi \parallel }_{{\mathcal{H}}_{\beta ,\sigma }}^2& :={\parallel u\parallel }_{W_0^{\beta ,2}}^2+{\parallel {\eta }^t\parallel }_{{\mathcal{M}}_{\sigma -1}}^2.\hfill \end{array}$$

Then, for each $M\ge 0$, define the closed subset

$$\begin{array}{c}\hfill {\mathcal{H}}_{\beta ,\sigma }^M{=\{\varphi =(u,\eta )}^{tr}\in {\mathcal{H}}_{\beta ,\sigma }:\left|\langle u\rangle \right|\le M\}.\end{array}$$

(36)

When we are concerned with the dynamical system associated with the model Problem P, we utilize the following metric space,

$$\begin{array}{c}\hfill {\mathcal{X}}_{\beta ,\sigma }^M:=\left\{\varphi ={(u,\eta )}^{tr}\in {\mathcal{H}}_{\beta ,\sigma }^M:F\left(u\right)\in L^1\left(\mathsf{\Omega }\right)\right\},\end{array}$$

endowed with the metric

$$\begin{array}{c}\hfill d_{{\mathcal{X}}_{\beta ,\sigma }^M}({\varphi }_1,{\varphi }_2):={\parallel {\varphi }_1-{\varphi }_2\parallel }_{{\mathcal{H}}_{\beta ,\sigma }^M}+{\left|{\int }_{\mathsf{\Omega }}F\left(u_1\right)dx-{\int }_{\mathsf{\Omega }}F\left(u_2\right)dx\right|}^{1/2}.\end{array}$$

Remark 2.

The embedding ${\mathcal{H}}_{\beta ,1}^M\hookrightarrow {\mathcal{H}}_{\beta ,0}^M$ is continuous but not compact, due to the presence of the second component ${\mathcal{M}}_{\sigma -1}^{\left(0\right)}$. Indeed, see [40] for a counterexample.

It is appropriate for us to state the various assumptions that may be used on the kernel $\nu$.

(K1) 

$\nu \in C^1\left({\mathbb{R}}_{+}\right)\cap L^1\left({\mathbb{R}}_{+}\right)$ and $\nu \left(s\right)\ge 0$ for all $s\in {\mathbb{R}}_{+}$.

(K2) 

${\nu }^{\prime }\left(s\right)\le 0$ for all $s\in {\mathbb{R}}_{+}$.

(K3) 

${\displaystyle k_0={\int }_0^{\infty }\nu \left(s\right)ds>0}$. (For the sake of simplicity, we now assume $k_0=1$ throughout the rest of the paper.)

(K4) 

${\displaystyle {\nu }_0=\underset{s\to 0^{+}}{lim}\nu \left(s\right)<\infty .}$

(K5) 

${\nu }^{\prime }\left(s\right)+\lambda \nu \left(s\right)\le 0$ for a.e. $s\in {\mathbb{R}}_{+}$, for some $\lambda >0.$

Some remarks for these assumptions: By assumption (K2), the inequality holds for all ${\eta }^t\in D\left(T_r\right)$

$${(T_r{\eta }^t,{\eta }^t)}_{{\mathcal{M}}_{-1}}\le 0.$$

(37)

We remind the reader that the assumption (K5) is only required when we examine the asymptotic behavior of the solutions (and in that case, (K2) is redundant).

In order to formulate a suitable (abstract) evolution equation for ${\eta }^t$, we define the linear operator $T_r=-{\partial }_s$ with the domain

$$D\left(T_r\right)=\{{\eta }^t\in {\mathcal{M}}_{-1}^{\left(0\right)}:{\partial }_s{\eta }^t\in {\mathcal{M}}_{-1}^{\left(0\right)},{\eta }^t\left(0\right)=0\}.$$

It is well-known that $T_r$ is the infinitesimal generator of the right-translation semigroup on ${\mathcal{M}}_{-1}$; indeed, the following result comes from [37] (Theorem 3.1).

Proposition 1.

The operator $T_r$ with domain $D\left(T_r\right)$ is an infinitesimal generator of a strongly continuous semigroup of contractions on ${\mathcal{M}}_{-1}$, denoted $e^{T_{r}t}$.

As a consequence, we also have (see, e.g., [41] (Corollary IV.2.2)).

Corollary 1.

Let $T>0$ and assume $g\in L^1(0,T;H^{-1}\left(\mathsf{\Omega }\right))$. Then, for every ${\eta }_0\in {\mathcal{M}}_{-1}$, the Cauchy problem for ${\eta }^t,$

$$\left\{\begin{array}{cc}{\partial }_t{\eta }^t=T_r{\eta }^t+g\left(t\right),\hfill & \mathrm{for}t>0,\hfill \\ {\eta }^0={\eta }_0,\hfill \end{array}\right.$$

(38)

has a unique (mild) solution $\eta \in C([0,T];{\mathcal{M}}_{-1})$ which can be explicitly given as

$${\eta }^t\left(s\right)=\left\{\begin{array}{cc}{\displaystyle {\int }_0^{s}g(t-y)dy,}\hfill & \mathrm{for}0<s\le t,\hfill \\ {\displaystyle {\eta }_0(s-t)+{\int }_0^{t}g(t-y)dy,}\hfill & \mathrm{for}s>t,\hfill \end{array}\right.$$

(39)

see also [21] (Section 3.2) and [37] (Section 3).

3. Variational Formulation and Well-Posedness

To begin this section, we state the assumptions on the nonlinear term F and report some important consequences of these assumptions. These assumptions on F are based on [13,15] and can be found in [5] (Section 3).

(N1) 

$F\in C_{loc}^2\left(\mathbb{R}\right)$ and there exists $c_F>0$ such that, for all $r\in \mathbb{R},$

$$\begin{array}{c}\hfill F^{″}\left(r\right)\ge -c_F.\end{array}$$

(N2) 

There exist $c_F>0$ and $p\in (1,2]$ such that, for all $r\in \mathbb{R},$

$$\begin{array}{c}\hfill |F^{\prime }{\left(r\right)|}^p\le c_F\left(\right|F\left(r\right)|+1).\end{array}$$

(N3) 

There exist $C_1,C_2>0$ such that, for all $r\in \mathbb{R},$

$$\begin{array}{c}\hfill F\left(r\right)\ge C_1{\left|r\right|}^{p/(p-1)}-C_2.\end{array}$$

The last assumption is not needed to obtain the existence of weak solutions, but it is relied upon later when we seek the existence of strong/regular solutions and uniqueness of these solutions.

(N4) 

There exist $\rho \ge 2$ and $C_3>0$ such that, for all $r\in \mathbb{R},$

$$\begin{array}{c}\hfill |F^{″}\left(r\right)|\le C_3{(1+|r|}^{\rho -2}).\end{array}$$

(40)

The following remarks are from [5]. Assumption (N1) implies that the potential F is a quadratic perturbation of some strictly convex function; i.e., there holds,

$$F\left(r\right)=G\left(r\right)-\frac{c_F}{2}{r}^2,$$

(41)

with $G\in C^2\left(\mathbb{R}\right)$ strictly convex as $G^{″}\ge 0$ in $\mathsf{\Omega }.$ Furthermore, with (N1), for each $M\ge 0$ there are constants $C_i>0$, $i=3,\dots ,6$, (with $C_4$ and $C_5$ depending on M and F) such that, for all $r\in \mathbb{R},$

$$\begin{array}{c}\hfill F\left(r\right)-C_3\le C_4{(r-M)}^2+F^{\prime }\left(r\right)(r-M),\end{array}$$

(42)

$$\begin{array}{c}\hfill \frac{1}{2}|F^{\prime }\left(s\right)\left|\right(1+\left|r\right|)\le F^{\prime }\left(r\right)(r-M)+C_5,\end{array}$$

(43)

(see [26] (Equations (4.7) and (4.8))) and

$$\begin{array}{c}\hfill \left|F\left(r\right)\right|-C_6\le |F^{\prime }\left(r\right)\left|\right(1+\left|r\right|).\end{array}$$

(44)

The last inequality appears in [42] (page 8). With the positivity condition (N3), it follows that, for all $r\in \mathbb{R},$

$$|F^{\prime }\left(r\right)|\le c_F\left(\right|F\left(r\right)|+1).$$

(45)

Assumption (N2) allows for arbitrary polynomial growth $\overline{p}=p/(p-1)$ in the potential F. Significantly, the double-well potential $F\left(r\right)={(r^2-1)}^2$ satisfies (N2) with $p=4/3$ and (N4) with $p=2$.

We are now ready to introduce the variational/weak formulation of Problem P.

Definition 1.

Let $T>0$ and ${\varphi }_0={(u_0,{\eta }_0)}^{tr}\in {\mathcal{H}}_{\beta ,0}^M=W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\times {\mathcal{M}}_{-1}^{\left(0\right)}$ be such that $F\left(u_0\right)\in L^1\left(\mathsf{\Omega }\right).$ A pair $\varphi =(u,\eta )$ satisfying

$$\begin{array}{cc}\hfill \varphi =(u,\eta )& \in L^{\infty }(0,T;{\mathcal{H}}_{\beta ,0}^M),\hfill \end{array}$$

(46)

$$\begin{array}{cc}\hfill {\partial }_{t}u& \in L^2(0,T;H^{-1}\left(\mathsf{\Omega }\right)),\hfill \end{array}$$

(47)

$$\begin{array}{cc}\hfill {\partial }_t\eta & \in L^2(0,T;H_{\nu }^{-1}({\mathbb{R}}_{+};H_{\left(0\right)}^{-1}\left(\mathsf{\Omega }\right))),\hfill \end{array}$$

(48)

$$\begin{array}{cc}\hfill \mu & \in L^2(0,T;W^{-\beta ,2}\left(\mathsf{\Omega }\right)),\hfill \end{array}$$

(49)

$$\begin{array}{cc}\hfill F^{\prime }\left(u\right)& \in L^{\infty }(0,T;L^p\left(\mathsf{\Omega }\right))\hfill \end{array}$$

(50)

is called aweak solutionto ProblemPon $[0,T]$ with initial data ${\varphi }_0=(u_0,{\eta }_0)\in {\mathcal{H}}_{\beta ,0}^M$ if the following identities hold almost everywhere in $(0,T)$, and for all $v\in H^1\left(\mathsf{\Omega }\right),$$\xi \in W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\cap L^p\left(\mathsf{\Omega }\right)$, and $\zeta \in {\mathcal{M}}_1$:

$$\begin{array}{cc}\hfill {\langle {\partial }_{t}u,v\rangle }_{H^{-1}\times H^1}+{\int }_0^{\infty }\nu \left(s\right){\langle {\eta }^t\left(s\right),v\rangle }_{H^{-1}\times H^1}ds& =0,\hfill \end{array}$$

(51)

$$\begin{array}{cc}\hfill a_{E,\beta }(u,\xi )+{\langle F^{\prime }\left(u\right),\xi \rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}}+\alpha {\langle {\partial }_{t}u,\xi \rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}}& ={\langle \mu ,\xi \rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}},\hfill \end{array}$$

(52)

$$\begin{array}{cc}\hfill {({\partial }_t{\eta }^t,\zeta )}_{{\mathcal{M}}_{-1}}-{(T_r{\eta }^t,\zeta )}_{{\mathcal{M}}_{-1}}& ={(\mu ,\zeta )}_{{\mathcal{M}}_0}.\hfill \end{array}$$

(53)

Furthermore, the initial conditions hold in the $L^2$-sense

$$\begin{array}{c}\hfill u\left(0\right)=u_0\mathrm{and}{\eta }^0={\eta }_0.\end{array}$$

(54)

Finally, we say that $\varphi ={(u,\eta )}^{tr}$ is aglobal weak solutionof ProblemPif it is a weak solution on $[0,T]$, for any $T>0.$

Remark 3.

It is important to note that although ${\eta }_0$ is defined by (14) and (21), ${\eta }_0$may be taken to be initial data independent of u. Henceforth we consider a more general problem with respect to the original one.

Remark 4.

Concerning Equation (53) and the representation Formula (39), we have

$$T_r{\eta }^t\left(s\right)=-{\partial }_s{\eta }^t\left(s\right)=\left\{\begin{array}{cc}\mathsf{\Delta }\mu (t-s)\hfill & for0<s\le t,\hfill \\ -{\partial }_s{\eta }_0(s-t)\hfill & fors>t.\hfill \end{array}\right.$$

Thus, when given ${\eta }_0\in {\mathcal{M}}_{-1}^{\left(0\right)},$ then $T_r{\eta }^t\in H_{\nu }^{-1}({\mathbb{R}}_{+};H^{-1}\left(\mathsf{\Omega }\right)),$ for each $t\in (0,T)$, by virtue of (49). Moreover, taking $\zeta =1$ in the variational equation

$$\begin{array}{cc}\hfill {({\partial }_t{\eta }^t,\zeta )}_{{\mathcal{M}}_{-1}}-{(T_r{\eta }^t,\zeta )}_{{\mathcal{M}}_{-1}}& =-{\int }_0^{\infty }\nu \left(s\right){(-\mathsf{\Delta }\mu ,\zeta )}_{H^{-1}\times H^1}ds,\hfill \end{array}$$

we find, for all $s>t$,

$$\begin{array}{c}\hfill \frac{\partial }{\partial t}\langle {\eta }^t\left(s\right)\rangle +\frac{\partial }{\partial s}\langle {\eta }_0(s-t)\rangle =k_0\langle \mathsf{\Delta }\mu (t-s)\rangle .\end{array}$$

We know that ${\eta }_0\in {\mathcal{M}}_{-1}^{\left(0\right)}$ and $k_0=1$, hence

$$\begin{array}{c}\hfill \frac{\partial }{\partial t}\langle {\eta }^t\left(s\right)\rangle =0,\end{array}$$

and it follows that

$$\langle {\eta }^t\left(s\right)\rangle =0\forall t\ge 0.$$

Remark 5.

In the Cahn–Hilliard model, it is well-known that the average value of u is conserved (see, e.g., [43] (Section III.4.2)). A similar property holds here for our problem. Indeed, we may choose the test function $v=1$ in (51) which yields

$$\begin{array}{c}\hfill \frac{\partial }{\partial t}\langle u\left(t\right)\rangle +{\int }_0^{\infty }\nu \left(s\right)\langle {\eta }^t\left(s\right)\rangle ds=0.\end{array}$$

By (4), there holds $\langle {\eta }^t\left(s\right)\rangle =0$ for all $t>0$ and for all $s>0$. Hence, we recover the conservation of mass

$$\begin{array}{c}\hfill \langle u\left(t\right)\rangle =\langle u_0\rangle \mathrm{and}\langle {\partial }_{t}u\left(t\right)\rangle =0,\forall t\ge 0.\end{array}$$

(55)

Remark 6.

Before we continue to the existence statement, it is worthwhile to recall Theorem A1 (d) in Appendix A for which the following embedding holds

$$D\left(A_{E,\beta }\right)\hookrightarrow L^{\infty }\left(\mathsf{\Omega }\right),\forall \beta \in (\frac{N}{4},1),\mathrm{for}N=1,2,3.$$

(56)

Theorem 1.

Let $T>0$ and ${\varphi }_0={(u_0,{\eta }_0)}^{tr}\in {\mathcal{H}}_{\beta ,0}^M=W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\times {\mathcal{M}}_{-1}^{\left(0\right)}$ for $\beta \in (\frac{N}{4},1)$, $N=1,2,3$, be such that $F\left(u_0\right)\in L^1\left(\mathsf{\Omega }\right).$ Assume $\alpha >0$ and that (K1)–(K4) and (N1)–(N3) hold. ProblemPadmits at least one weak solution $\varphi =(u,\eta )$ on $(0,T)$ according to Definition 1 with the additional regularity

$$\begin{array}{cc}\hfill u& \in L^{\infty }(0,T;L^{p/(p-1)}\left(\mathsf{\Omega }\right)),\hfill \end{array}$$

(57)

$$\begin{array}{cc}\hfill \sqrt{\alpha }{\partial }_{t}u& \in L^2(\mathsf{\Omega }\times (0,T)),\hfill \end{array}$$

(58)

$$\begin{array}{cc}\hfill \eta & \in L^2(0,T;L_{-{\nu }^{\prime }}^2({\mathbb{R}}_{+};H_{\left(0\right)}^{-1}\left(\mathsf{\Omega }\right))),\hfill \end{array}$$

(59)

$$\begin{array}{cc}\hfill F\left(u\right)& \in L^{\infty }(0,T;L^1\left(\mathsf{\Omega }\right)),F^{\prime }\left(u\right)\in L^{\infty }(0,T;L^1\left(\mathsf{\Omega }\right)).\hfill \end{array}$$

(60)

for any $T>0.$ Furthermore, setting

$$\begin{array}{c}\hfill \mathcal{E}\left(t\right):={|\parallel u\left(t\right)\parallel |}_{W_0^{\beta ,2}}^2+2(F\left(u\left(t\right)\right),1)+{\parallel {\eta }^t\parallel }_{{\mathcal{M}}_{-1}}^2+C\end{array}$$

(61)

for some $C>0$ sufficiently large, the following energy equality holds for every such weak solution,

$$\begin{array}{c}\hfill \mathcal{E}\left(t\right)+2{\int }_0^t\left(\alpha \parallel {\partial }_t{u\left(\tau \right)\parallel }^{2}d\tau -{\int }_0^{\infty }{\nu }^{\prime }\left(s\right){\parallel {\eta }^{\tau }\left(s\right)\parallel }_{H^{-1}}^{2}ds\right)d\tau =\mathcal{E}\left(0\right).\end{array}$$

(62)

Proof. 

The proof proceeds in several steps. The existence proof begins with a Faedo–Galerkin approximation procedure in which we later pass to the limit. We first assume that $u_0\in D\left(A_{E,\beta }\right)$. (This assumption will be used to show that there is a sequence ${\left\{u_{0n}\right\}}_{n=1}^{\infty }$ such that $u_{0n}\to u_0$ in $D\left(A_{E,\beta }\right)$ as well as $L^{\infty }\left(\mathsf{\Omega }\right)$ per (56), which will be important in light of the fact that $F\left(u_{0n}\right)$ is of arbitrary polynomial growth per assumptions (N1)–(N3).) The existence of a weak solution for $u_0\in W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$ with $F\left(u_0\right)\in L^1\left(\mathsf{\Omega }\right)$ follows from a density argument. To establish the equality in the energy identity, we exploit the fact that the potential F is a quadratic perturbation of some strictly convex function.

Step 1: The Galerkin approximation. To begin, we introduce the family ${\left\{v_j\right\}}_{j\ge 1}$ of eigenvectors of the fractional Laplacian $A_{E,\beta }$, which exist thanks to Theorem A1 in Appendix A. Moreover, there is a family ${\left\{w_j\right\}}_{j\ge 1}$ consisting of the eigenvectors of the Neumann–Laplacian $A_N$, and with this, we define the smooth sequence of ${\left\{z_j\right\}}_{j\ge 1}\subset D\left(T_r\right)\cap W_{\nu }^{1,2}({\mathbb{R}}_{+};H_{\left(0\right)}^1\left(\mathsf{\Omega }\right))$ by $z_j=b_j{w}_j$ such that ${\left\{b_j\right\}}_{j\ge 1}\subset C_c^{\infty }\left({\mathbb{R}}_{+}\right)$ is an orthonormal basis for $L_{\nu }^2\left({\mathbb{R}}_{+}\right)$. Using these, we define the following finite-dimensional spaces:

$$\begin{array}{c}\hfill V^n=\mathrm{span}\{v_1,v_2,\dots ,v_n\},W^n=\mathrm{span}\{w_1,w_2,\dots ,w_n\},{\mathcal{M}}^n=\mathrm{span}\{z_1,z_2,\dots ,z_n\},\end{array}$$

(63)

and set

$$V^{\infty }=\bigcup _{n=1}^{\infty }{V}^n,W^{\infty }=\bigcup _{n=1}^{\infty }{W}^n,{\mathcal{M}}^{\infty }=\bigcup _{n=1}^{\infty }{\mathcal{M}}^n.$$

Clearly, $V^{\infty }$ is a dense subspace of $W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$ and $W^{\infty }$ is a dense subspace of $H^1\left(\mathsf{\Omega }\right)$. In addition, ${\mathcal{M}}^{\infty }$ is a dense subspace of ${\mathcal{M}}_{-1}^{\left(0\right)}.$ For $T>0$ fixed, we look for two functions of the form on $(0,T)$,

$$u_n\left(t\right)=\sum _{k=1}^n{a}_k^{\left(n\right)}\left(t\right)v_k\mathrm{and}{\eta }_n^t\left(s\right)=\sum _{k=1}^n{c}_k^{\left(n\right)}\left(t\right)z_k,$$

(64)

where $a_j^{\left(n\right)}$ and $c_j^{\left(n\right)}$ are assumed to be (at least) $C^2\left([0,T]\right)$ for each $j=1,2,\dots$ an for each $n=1,2,\dots ,$ which solve the following approximating Problem ${\mathbf{P}}_n$:

$$\begin{array}{cc}\hfill ({\partial }_t{u}_n,v)+{\int }_0^{\infty }\nu \left(s\right)({\eta }_n^t\left(s\right),v)ds& =0\hfill \end{array}$$

(65)

$$\begin{array}{cc}\hfill a_{E,\beta }(u_n,\xi )+(F^{\prime }\left(u_n\right),\xi )+\alpha ({\partial }_t{u}_n,\xi )& =({\mu }_n,\xi )\hfill \end{array}$$

(66)

$$\begin{array}{cc}\hfill {({\partial }_t{\eta }_n^t,\zeta )}_{{\mathcal{M}}_{-1}}-{(T_r{\eta }_n^t,\zeta )}_{{\mathcal{M}}_{-1}}& ={({\mu }_n,\zeta )}_{{\mathcal{M}}_0}\hfill \end{array}$$

(67)

$$\begin{array}{c}\hfill u_n\left(0\right)=u_{0n},{\eta }_n^0={\eta }_{0n}\end{array}$$

(68)

for every $v\in V^n,$$\xi \in W^n$ and $\zeta \in {\mathcal{M}}^n$, and where $u_{0n}$, and ${\eta }_{0n}$ denote the finite-dimensional projections of $u_0$ and ${\eta }_0$ onto $V^n$ and ${\mathcal{M}}^n$, respectively. This approximating problem is equivalent to solving a Cauchy problem for a system of ordinary differential equations (indeed, see, e.g., [26] (page 131)). Hence, the Cauchy–Lipschitz theorem ensures that there exists a $T_n\in (0,\infty ]$ such that this approximating system has a unique maximal solution.

Step 2: A priori estimates. We now derive some a priori estimates in order to show that $T_n=\infty$ for every $n\ge 1$ and that the sequences of $u_n,{\eta }_n^t,{\mu }_n$ are bounded in suitable functional spaces. By using $v={\mu }_n$ as a test function in (65) and $\xi ={\partial }_t{u}_n$ as a test function in (66) we obtain

$$\begin{array}{cc}\hfill & ({\partial }_t{u}_n,{\mu }_n)+{\int }_0^{\infty }\nu \left(s\right)({\eta }_n^t\left(s\right),{\mu }_n)ds=0\hfill \end{array}$$

(69)

$$\begin{array}{cc}\hfill & ({\mu }_n,{\partial }_t{u}_n)=({(-\mathsf{\Delta })}_E^{\beta }{u}_n,{\partial }_t{u}_n)+(F^{\prime }\left(u_n\right),{\partial }_t{u}_n)+\alpha {\parallel {\partial }_t{u}_n\parallel }^2,\hfill \end{array}$$

(70)

and taking $\zeta ={\eta }_n^t$ as a test function in (67) yields (for the products in ${\mathcal{M}}_{-1}$, this is a multiplication by ${(-\mathsf{\Delta })}^{-1}{\eta }_n^t$ in ${\mathcal{M}}_0$)

$$\begin{array}{cc}\hfill {\int }_0^{\infty }\nu \left(s\right)\left({\int }_{\mathsf{\Omega }}{\partial }_t{\eta }_n^t(x,s){(-\mathsf{\Delta })}^{-1}{\eta }_n^t(x,s)dx\right)ds& +{\int }_0^{\infty }\nu \left(s\right)\left({\int }_{\mathsf{\Omega }}{\partial }_s{\eta }_n^t(x,s){(-\mathsf{\Delta })}^{-1}{\eta }_n^t(x,s)dx\right)ds\hfill \\ \hfill & ={\int }_0^{\infty }\nu \left(s\right)\left({\int }_{\mathsf{\Omega }}(-\mathsf{\Delta }){\mu }_n(x,t){(-\mathsf{\Delta })}^{-1}{\eta }_n^t(x,s)dx\right)ds,\hfill \end{array}$$

which is, after an integration by parts,

$$\begin{array}{c}\hfill {({\partial }_t{\eta }_n^t,{\eta }_n^t)}_{{\mathcal{M}}_{-1}}+{({\partial }_s{\eta }_n^t,{\eta }_n^t)}_{{\mathcal{M}}_{-1}}={({\mu }_n,{\eta }_n^t)}_{{\mathcal{M}}_0}.\end{array}$$

(71)

Then, combining the results produces the differential identity, which holds for almost all $t\in (0,T),$

$$\begin{array}{c}\hfill \frac{1}{2}\frac{d}{dt}\left\{|\parallel u_n{\parallel |}_{W_0^{\beta ,2}}^2+2(F\left(u_n\right),1)+{\parallel {\eta }^t\parallel }_{{\mathcal{M}}_{-1}}^2\right\}+\alpha {\parallel {\partial }_t{u}_n\parallel }^2-{(T_r{\eta }_n^t,{\eta }_n^t)}_{{\mathcal{M}}_{-1}}=0.\end{array}$$

(72)

For all $t\in (0,T_n),$ set

$${\mathcal{E}}_n\left(t\right):=|\parallel u_n{\left(t\right)\parallel |}_{W_0^{\beta ,2}}^2+2(F\left(u_n\left(t\right)\right),1)+{\parallel {\eta }_n^t\parallel }_{{\mathcal{M}}_{-1}}^2+C$$

(73)

where in light of (N3), the functional ${\mathcal{E}}_n\left(t\right)$ is nonnegative for all $t\in (0,T_n)$. We have

$$\begin{array}{c}\hfill \frac{d}{dt}{\mathcal{E}}_n+2\alpha \parallel {\partial }_t{u}_n{\parallel }^2-2{\int }_0^{\infty }{\nu }^{\prime }\left(s\right){\parallel {\eta }_n^t\left(s\right)\parallel }_{H^{-1}}^{2}ds=0\end{array}$$

(74)

for almost all $t\in (0,T_n)$. Hence, integrating the equation above with respect to time in $(0,t)$, we are led to the following integral equality (which does hold for the approximate solutions)

$$\begin{array}{c}\hfill {\mathcal{E}}_n\left(t\right)+2{\int }_0^t\left(\alpha \parallel {\partial }_t{u}_n{\left(\tau \right)\parallel }^2-{\int }_0^{\infty }{\nu }^{\prime }\left(s\right){\parallel {\eta }_n^{\tau }\left(s\right)\parallel }_{H^{-1}}^{2}ds\right)d\tau ={\mathcal{E}}_n\left(0\right).\end{array}$$

(75)

Furthermore, from (73) and assumption (N3), we find the lower bound

$$\begin{array}{c}\hfill |\parallel u_n{\left(t\right)\parallel |}_{W_0^{\beta ,2}}^2+2C_1\parallel u_n{\left(t\right)\parallel }_{L^{p/(p-1)}}^{p/(p-1)}+{\parallel {\eta }_n^t\parallel }_{{\mathcal{M}}_{-1}}^2\le {\mathcal{E}}_n\left(t\right).\end{array}$$

(76)

Using the fact that $F\left(u_0\right)\in L^1\left(\mathsf{\Omega }\right)$, we also obtain the upper bound

$$\begin{array}{cc}\hfill {\mathcal{E}}_n\left(t\right)\le {\mathcal{E}}_n\left(0\right)& \le |\parallel u_n{\left(0\right)\parallel |}_{W_0^{\beta ,2}}^2+(F\left(u_n\left(0\right)\right),1)+{\parallel {\eta }_n^0\parallel }_{{\mathcal{M}}_{-1}}^2\hfill \\ \hfill & \le Q(\parallel {\varphi }_n\left(0\right){\parallel }_{{\mathcal{H}}_{\beta ,0}^M})+C.\hfill \end{array}$$

(77)

In particular, the uniform bound derived from (75)–(77) implies that the local solution to Problem ${\mathbf{P}}_{\mathbf{n}}$ can be extended up to time T, that is $T_n=T$, for every n. Moreover, from (75) and (76) we deduce the following bounds for the approximate solution

$$\begin{array}{cc}\hfill \parallel u_n{\parallel }_{L^{\infty }(0,T;W_0^{\beta ,2})}& \le C\hfill \end{array}$$

(78)

$$\begin{array}{cc}\hfill \parallel {\eta }_n{\parallel }_{L^{\infty }(0,T;{\mathcal{M}}_{-1})}& \le C\hfill \end{array}$$

(79)

$$\begin{array}{cc}\hfill \parallel F\left(u_n\right){\parallel }_{L^{\infty }(0,T;L^1)}& \le C\hfill \end{array}$$

(80)

$$\begin{array}{cc}\hfill \sqrt{\alpha }{\parallel {\partial }_t{u}_n\parallel }_{L^2(\mathsf{\Omega }\times (0,T))}& \le C\hfill \end{array}$$

(81)

$$\begin{array}{cc}\hfill \parallel {\eta }_n{\parallel }_{L^2(0,T;L_{-{\nu }^{\prime }}^2({\mathbb{R}}_{+};H^{-1}))}& \le C\hfill \end{array}$$

(82)

$$\begin{array}{cc}\hfill \parallel u_n{\parallel }_{L^{\infty }(0,T;L^{p/(p-1)})}& \le C.\hfill \end{array}$$

(83)

Obviously, (45) and (80) immediately show us

$$\begin{array}{cc}\hfill \parallel F^{\prime }\left(u_n\right){\parallel }_{L^{\infty }(0,T;L^1)}& \le C.\hfill \end{array}$$

(84)

Next, since $\langle A_N^{-1}{\partial }_t{u}_n\rangle =0$ (recall (55)), we may (and do) take $v=A_N^{-1}{\partial }_t{u}_n$ in (65) which leads us to the estimate,

$$\begin{array}{c}\hfill \parallel A_N^{-\frac{1}{2}}{\partial }_t{u}_n{\parallel }^2\le {\int }_0^{\infty }\nu \left(s\right)\parallel A_N^{-\frac{1}{2}}{\eta }_n^t\left(s\right)\parallel \parallel A_N^{-\frac{1}{2}}{\partial }_t{u}_n\left(t\right)\parallel ds,\end{array}$$

(85)

that is,

$$\begin{array}{c}\hfill \parallel {\partial }_t{u}_n{\parallel }_{H^{-1}}^2\le {\int }_0^{\infty }\nu \left(s\right)\parallel {\eta }_n^t{\left(s\right)\parallel }_{H^{-1}}{\parallel {\partial }_t{u}_n\parallel }_{H^{-1}}ds.\end{array}$$

(86)

Using the Cauchy–Schwartz inequality and assumption (K3), we can write

$$\begin{array}{c}\hfill \parallel {\partial }_t{u}_n{\parallel }_{H^{-1}}\le {\parallel {\eta }_n^t\parallel }_{{\mathcal{M}}_{-1}}.\end{array}$$

(87)

Thus, (79) and (87) yield

$$\begin{array}{c}\hfill \parallel {\partial }_t{u}_n{\parallel }_{L^{\infty }(0,T;H^{-1})}\le C.\end{array}$$

(88)

We need to bound $F^{\prime }\left(u_n\right)$, then ${\mu }_n$. In light of (66), we apply (84), (88), and the fact that operator $A_{E,\beta }$ is bounded from $W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$ into $W^{-\beta ,2}\left(\mathsf{\Omega }\right)$ (in particular, $\parallel A_{E,\beta }{u}_n{\parallel }_{L^2(0,T;W^{-\beta ,2}\left(\mathsf{\Omega }\right))}$$\le C)$, to obtain the following uniform bounds for ${\mu }_n$

$$\begin{array}{c}\hfill |\langle {\mu }_n\rangle |\le C,\end{array}$$

(89)

and

$$\begin{array}{c}\hfill \parallel {\mu }_n{\parallel }_{L^2(0,T;W^{-\beta ,2}\left(\mathsf{\Omega }\right))}\le C.\end{array}$$

(90)

This completes Step 2.

Step 3: Passage to the limit. On account of the above uniform inequalities, we can argue that there are functions u, $\eta$, $\mu$, such that, up to subsequences,

$$\begin{array}{ccc}\hfill u_n\rightharpoonup u& \mathrm{weakly}-*\mathrm{in}& L^{\infty }(0,T;W_0^{\beta ,2}\left(\mathsf{\Omega }\right)),\hfill \end{array}$$

(91)

$$\begin{array}{ccc}\hfill u_n\rightharpoonup u& \mathrm{weakly}-*\mathrm{in}& L^{\infty }(0,T;L^{p/(p-1)}\left(\mathsf{\Omega }\right)),\hfill \end{array}$$

(92)

$$\begin{array}{ccc}\hfill {\partial }_t{u}_n\rightharpoonup {\partial }_{t}u& \mathrm{weakly}-*\mathrm{in}& L^{\infty }(0,T;H^{-1}\left(\mathsf{\Omega }\right)),\hfill \end{array}$$

(93)

$$\begin{array}{ccc}\hfill \sqrt{\alpha }{\partial }_t{u}_n\rightharpoonup \sqrt{\alpha }{\partial }_{t}u& \mathrm{weakly}\mathrm{in}& L^2(\mathsf{\Omega }\times (0,T)),\hfill \end{array}$$

(94)

$$\begin{array}{ccc}\hfill {\eta }_n\rightharpoonup \eta & \mathrm{weakly}-*\mathrm{in}& L^{\infty }(0,T;{\mathcal{M}}_{-1}),\hfill \end{array}$$

(95)

$$\begin{array}{ccc}\hfill {\eta }_n\rightharpoonup \eta & \mathrm{weakly}\mathrm{in}& L^2(0,T;L_{-{\nu }^{\prime }}^2({\mathbb{R}}_{+};H^{-1}(\mathsf{\Omega }))),\hfill \end{array}$$

(96)

$$\begin{array}{ccc}\hfill {\partial }_t{\eta }_n\rightharpoonup {\partial }_t\eta & \mathrm{weakly}\mathrm{in}& L^2(0,T;H_{\nu }^{-1}({\mathbb{R}}_{+};H^{-1}(\mathsf{\Omega }))),\hfill \end{array}$$

(97)

$$\begin{array}{ccc}\hfill {\mu }_n\rightharpoonup \mu & \mathrm{weakly}\mathrm{in}& L^2(0,T;W^{-\beta ,2}\left(\mathsf{\Omega }\right)).\hfill \end{array}$$

(98)

(Note that (97) is due to (67) and the definition of the operator $T_r$.) Using the above convergences (91) and (93), as well as the fact that the injection $W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\hookrightarrow L^2\left(\mathsf{\Omega }\right)$ is compact for any $\beta \in (0,1)$, we draw upon the conclusion of the Aubin–Lions Lemma (see Lemma A1 in Appendix A) to deduce the following embedding is compact

$$\begin{array}{c}\hfill W:=\{\chi \in L^2(0,T;W_0^{\beta ,2}(\mathsf{\Omega })):{\partial }_t\chi \in L^2(0,T;H^{-1}(\mathsf{\Omega }))\}\hookrightarrow L^2(\mathsf{\Omega }\times (0,T)).\end{array}$$

(99)

Hence,

$$\begin{array}{c}\hfill u_n\to u\mathrm{strongly}\mathrm{in}{L}^2(\mathsf{\Omega }\times (0,T)),\end{array}$$

(100)

and we deduce that $u_n$ converges to u, almost everywhere in $\mathsf{\Omega }\times (0,T)$. Using assumption (N1) with (100), we deduce

$$\begin{array}{c}\hfill F^{\prime }\left(u_n\right)\to F^{\prime }\left(u\right)\mathrm{strongly}\mathrm{in}{L}^2(0,T;L^1\left(\mathsf{\Omega }\right)).\end{array}$$

(101)

Thus, we now have all the sufficient convergence results to pass to the limit in Equations (65) and (66) in order to recover (16) and (17), respectively. It remains to recover Equation (67) after we pass to the limit. An integration by parts on the first term in (67) and then an application of (95) yields, for any $\zeta \in C_0^{\infty }((0,T);C_0^{\infty }((0,T);H^1\left(\mathsf{\Omega }\right)))$

$$\begin{array}{c}\hfill {\int }_0^T{({\partial }_t{\eta }_n^{\tau },\zeta )}_{{\mathcal{M}}_{-1}}d\tau =-{\int }_0^T{({\eta }_n^{\tau },{\partial }_t\zeta )}_{{\mathcal{M}}_{-1}}d\tau \to -{\int }_0^T{({\eta }^{\tau },{\partial }_t\zeta )}_{{\mathcal{M}}_{-1}}d\tau .\end{array}$$

(102)

With this, we have

$$\begin{array}{c}\hfill {\partial }_t{\eta }_n^t\rightharpoonup {\partial }_t{\eta }^t\mathrm{weakly}\mathrm{in}{L}^2(0,T;H_{\nu }^{-1}({\mathbb{R}}_{+};H^{-1}\left(\mathsf{\Omega }\right)))\end{array}$$

(103)

and that ${\eta }^t\in L^{\infty }(0,T;H_{\nu }^{-1}({\mathbb{R}}_{+};H^{-1}\left(\mathsf{\Omega }\right)))$. Furthermore, with the help of (96), we have

$$\begin{array}{c}\hfill -{\int }_0^T{(T_r{\eta }_n^{\tau },\zeta )}_{{\mathcal{M}}_{-1}}d\tau =-{\int }_0^T{\nu }^{\prime }\left(s\right){({\eta }_n^{\tau },\zeta )}_{H^{-1}}d\tau \to -{\int }_0^T{\nu }^{\prime }\left(s\right){({\eta }^{\tau },\zeta )}_{H^{-1}}d\tau .\end{array}$$

(104)

By using a density argument (see [37]) and the following distributional equality

$$\begin{array}{c}\hfill -{\int }_0^T{({\eta }_n^{\tau },{\partial }_t\zeta )}_{{\mathcal{M}}_{-1}}d\tau -{\int }_0^T{\nu }^{\prime }\left(s\right){({\eta }^{\tau },\zeta )}_{H^{-1}\left(\mathsf{\Omega }\right)}d\tau ={\int }_0^T{({\partial }_t{\eta }^{\tau }-T_r{\eta }^{\tau },\zeta )}_{{\mathcal{M}}_{-1}}d\tau ,\end{array}$$

(105)

we also get (67) on account of (95) and (98). This completes Step 3 of the proof.

Step 4: Energy equality. To begin, let $u_0\in D\left(A_{E,\beta }\right)$, ${\eta }_0\in {\mathcal{M}}_{-1}^{\left(0\right)}$ and let $\varphi ={(u,\eta )}^{tr}$ be the corresponding weak solution. Recall from (100), we have, for almost all $t\in (0,T),$

$$\begin{array}{c}\hfill u_n\left(t\right)\to u\left(t\right)\mathrm{strongly}\mathrm{in}{L}^2\left(\mathsf{\Omega }\right)\mathrm{and}\mathrm{a}.\mathrm{e}.\mathrm{in}\mathsf{\Omega }.\end{array}$$

(106)

Since F is measurable and positive (see (N1) and (N3), respectively), Fatou’s lemma implies

$${\int }_{\mathsf{\Omega }}F\left(u\left(t\right)\right)dx\le \underset{n\to +\infty }{lim\; inf}{\int }_{\mathsf{\Omega }}F\left(u_n\left(t\right)\right)dx.$$

(107)

Passing to the limit in (75), and while keeping in mind (91), (94), (95), (97), (98), and (101), as well as the weak lower-semicontinuity of the norm, we arrive at the integral inequality which holds for any weak solution

$$\begin{array}{c}\hfill \mathcal{E}\left(t\right)+2{\int }_0^t\left(\alpha \parallel {\partial }_t{u\left(\tau \right)\parallel }^{2}d\tau -{\int }_0^{\infty }{\nu }^{\prime }\left(s\right){\parallel {\eta }^{\tau }\left(s\right)\parallel }_{H^{-1}}^{2}ds\right)d\tau \le \mathcal{E}\left(0\right).\end{array}$$

We argue as in the proof of [12] (Corollary 2) to establish the energy equality. Indeed, take $\xi =\mu$ in (51). By (17), we need to treat the dual pairing ${\langle F^{\prime }\left(u\right),{\partial }_{t}u\rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}}$. It is here where we employ (41), where $F^{\prime }\left(u\right)=G^{\prime }\left(u\right)-c_{F}u$ and $G^{\prime }\in C^1\left(\mathbb{R}\right)$ is monotone increasing. Define the functional $\mathcal{G}:L^2\left(\mathsf{\Omega }\right)\to \mathbb{R}$ by

$$\mathcal{G}\left(\varphi \right):=\left\{\begin{array}{cc}{\displaystyle {\int }_{\mathsf{\Omega }}G\left(u\right)dx}\hfill & \mathrm{if}G\left(u\right)\in L^1\left(\mathsf{\Omega }\right),\hfill \\ +\infty \hfill & \mathrm{otherwise}.\hfill \end{array}\right.$$

Now, by [44] (Proposition 2.8, Chapter II), it follows that $\mathcal{G}$ is convex, lower-semi- continuous on $L^2\left(\mathsf{\Omega }\right)$, and $\chi \in \partial \mathcal{G}\left(u\right)$ if and only if $\chi =G^{\prime }\left(u\right)$ almost everywhere in $\mathsf{\Omega }$. Since we have (47), we apply [45] (Proposition 4.2) to find that there holds, for almost all $t\in (0,T),$

$$\begin{array}{cc}\hfill {\langle {\partial }_{t}u,F^{\prime }\left(u\right)\rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}}& ={\langle {\partial }_{t}u,G^{\prime }\left(u\right)\rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}}-c_F{\langle {\partial }_{t}u,u\rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}}\hfill \\ \hfill & =\frac{d}{dt}\left\{\mathcal{G}\left(u\right)-\frac{c_F}{2}{\parallel u\parallel }^2\right\}\hfill \\ \hfill & =\frac{d}{dt}{\int }_{\mathsf{\Omega }}F\left(u\right)dx.\hfill \end{array}$$

Similar to Step 2 above, take $v=\mu$, $\xi ={\partial }_{t}u$, and $\zeta ={\eta }^t$ (now without the index n) in (51)–(53), respectively. Using the above result on the dual product with $F^{\prime }\left(u\right)$ and (47), we are led to the differential identity (74) with E, u, and $\eta$ in place of ${\mathcal{E}}_n$, $u_n$, and ${\eta }_n$, respectively. Integrating the resulting differential identity on $(0,t)$ produces (62) as claimed. This completes Step 4.

Step 5: $(u,\eta )$ weak solution to ProblemP. Now let us take ${\varphi }_0={(u_0,{\eta }_0)}^{tr}\in {\mathcal{H}}_{\beta ,0}^M$ where $F\left(u_0\right)\in L^1\left(\mathsf{\Omega }\right).$ Proceeding exactly as in [12] (page 440) the bounds (78)–(84) and (88)–(90) hold. Moreover, with the aid of the Aubin–Lions compact embedding (again see Lemma A1 in Appendix A below) we deduce the existence of functions u, $\eta$, and $\mu$ that satisfy (46), (49), (57), and (59). Thus, passing to the limit in the variational formulation for ${\varphi }_k={(u_k,{\eta }_k)}^{tr}$, we find $\varphi ={(u,\eta )}^{tr}$ is a solution corresponding to the initial data ${\varphi }_0={(u_0,{\eta }_0)}^{tr}\in {\mathcal{H}}_{\beta ,0}^M$ for which $F\left(u_0\right)\in L^1\left(\mathsf{\Omega }\right).$ This finishes the proof of the theorem. □

Before we continue, we make some important remarks.

Remark 7.

The continuity property

$$\begin{array}{c}\hfill u\in C([0,T];W_0^{\beta -\iota ,2}\left(\mathsf{\Omega }\right)),\end{array}$$

for any $\iota >0$ sufficiently small follows from the conditions in Definition 1 after an application of the Aubin–Lions Lemma (see Lemma A1 in Appendix A). In addition, the property

$$\begin{array}{c}\hfill \eta \in C([0,T];{\mathcal{M}}_{-1}^{\left(0\right)})\end{array}$$

follows from the density argument in [37]. Thus, we deduce the continuity properties

$$\begin{array}{c}\hfill \varphi =(u,\eta )\in C([0,T];{\mathcal{H}}_{\beta ,0}^M).\end{array}$$

Remark 8.

From (62), we see that if there is a $t^{*}>0$ in which

$$\mathcal{E}\left(t^{*}\right)=\mathcal{E}\left(0\right),$$

then, for all $t\in (0,t^{*}),$

$${\int }_0^t\left(\alpha \parallel {\partial }_t{u\left(\tau \right)\parallel }^2+{\parallel {\eta }^{\tau }\parallel }_{L_{-{\nu }^{\prime }}^2({\mathbb{R}}_{+};H^{-1})}^2\right)d\tau =0.$$

(108)

We deduce ${\partial }_{t}u\left(t\right)=0$ for all $t\in (0,t^{*})$. Additionally, since $u\left(t\right)=u_0$ for all $t\in (0,t^{*}),$ Equation (17) shows

$$\mu \left(t\right)=A_{E,\beta }{u}_0+F^{\prime }\left(u_0\right)\forall t\in (0,t^{*}),$$

i.e., $\mu \left(t\right)={\mu }^{*}$ is also stationary. Thus, by the definition of ${\eta }^t$ given in (14), we find here that, for each $t\in (0,t^{*})$

$${\eta }^t\left(s\right)=s{A}_N{\mu }^{*}\forall s\ge 0.$$

Therefore, $\varphi ={(u,\eta )}^{tr}$ is a fixed point of the trajectory $\varphi \left(t\right)=\mathcal{S}\left(t\right){\varphi }_0$, where $\mathcal{S}$ is the solution operator defined below in Corollary 2.

The following result (see [26] (Theorem 3.4)) concerns the existence of strong/regular solutions which is utilized in the proof of the continuous dependence estimate. Note that we now employ the added assumption on the nonlinear term.

Theorem 2.

Let $T>0$, $\beta \in (0,1)$, and ${\varphi }_0={(u_0,{\eta }_0)}^{tr}\in {\mathcal{H}}_{\beta +1,\beta +1}^M:=W_0^{\beta +1,2}\left(\mathsf{\Omega }\right)\times L_{\nu }^2({\mathbb{R}}_{+};W_0^{\beta ,2}\left(\mathsf{\Omega }\right))$ be such that $F\left(u_0\right)\in L^1\left(\mathsf{\Omega }\right)$ and ${\eta }_0\in D\left(T_r\right)$. Assume $\alpha >0$ and that (K1)–(K4) and (N1)–(N3) hold. Additionally, assume that (N4) holds. ProblemPadmits at least one weak solution $\varphi =(u,\eta )$ on $(0,T)$ according to Definition (1) with the additional regularity, for any $T>0$,

$$\begin{array}{cc}\hfill \varphi =(u,\eta )& \in L^{\infty }(0,T;{\mathcal{H}}_{\beta +1,\beta +1}^M)\cap W^{1,\infty }(0,T;{\mathcal{H}}_{\beta ,0}^M),\hfill \end{array}$$

(109)

$$\begin{array}{cc}\hfill \sqrt{\alpha }{\partial }_{t}u& \in L^2(0,T;H^1\left(\mathsf{\Omega }\right))\hfill \end{array}$$

(110)

$$\begin{array}{cc}\hfill {\partial }_{tt}u& \in L^{\infty }(0,T;H^{-1}\left(\mathsf{\Omega }\right)),\hfill \end{array}$$

(111)

$$\begin{array}{cc}\hfill \sqrt{\alpha }{\partial }_{tt}u& \in L^2(\mathsf{\Omega }\times (0,T)),\hfill \end{array}$$

(112)

$$\begin{array}{cc}\hfill \mu & \in L^{\infty }(0,T;H^1\left(\mathsf{\Omega }\right)),\hfill \end{array}$$

(113)

$$\begin{array}{cc}\hfill \eta & \in L^{\infty }(0,T;D\left(T_r\right)).\hfill \end{array}$$

(114)

Proof. 

The proof relies on the Galerkin approximation scheme developed in the proof of Theorem 1. We seek ${\varphi }_n=(u_n,{\eta }_n)$ of the form (64) satisfying Problem ${\mathbf{P}}_n$:

$$\begin{array}{cc}\hfill ({\partial }_{tt}{u}_n,v)+{\int }_0^{\infty }\nu \left(s\right)({\partial }_t{\eta }_n^t\left(s\right),v)ds& =0\hfill \end{array}$$

(115)

$$\begin{array}{cc}\hfill a_{E,\beta }({\partial }_t{u}_n,\xi )+(F^{″}\left(u_n\right){\partial }_{t}u,\xi )+\alpha ({\partial }_{tt}{u}_n,\xi )& =({\partial }_t{\mu }_n,\xi )\hfill \end{array}$$

(116)

$$\begin{array}{cc}\hfill {({\partial }_{tt}{\eta }_n^t,\zeta )}_{{\mathcal{M}}_{-1}}-{(T_r{\partial }_t{\eta }_n^t,\zeta )}_{{\mathcal{M}}_{-1}}& ={({\partial }_t{\mu }_n,\zeta )}_{{\mathcal{M}}_0}\hfill \end{array}$$

(117)

for every $t\in (0,T),$$v\in V^n$, $\xi \in W^n$, and $\zeta \in {\mathcal{M}}^n,$ and which satisfy the initial conditions

$$\begin{array}{c}\hfill u_n\left(0\right)={\tilde{u}}_{0n}\mathrm{and}{\eta }_n^0={\tilde{\eta }}_{0n},\end{array}$$

(118)

where we set

$${\tilde{u}}_{0n}:=-{\int }_0^{\infty }\nu \left(s\right){\eta }_{0n}\left(s\right)ds,$$

(119)

and

$${\tilde{\eta }}_{0n}:=T_r{\eta }_{0n}+A_N{\mu }_{0n},$$

(120)

and also

$${\mu }_{0n}=-\alpha {\int }_0^{\infty }\nu \left(s\right){\eta }_{0n}\left(s\right)ds+A_{E,\beta }{u}_{0n}+F^{\prime }\left(u_{0n}\right).$$

(121)

It is important to note that when ${\varphi }_0=(u_0,{\eta }_0)$ satisfies the assumptions of Theorem 2, then it is guaranteed that $({\tilde{u}}_0,{\tilde{\eta }}_0)\in {\mathcal{H}}_{1,0}^M$. Indeed, relying on the fact that $\parallel (u_{0n},{\eta }_{0n}){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}\le {\parallel (u_0,{\eta }_0)\parallel }_{{\mathcal{H}}_{\beta ,0}^M}$, we easily obtain the estimate $\parallel ({\partial }_t{u}_n\left(0\right),{\partial }_t{\eta }_n^0){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}\le Q(\parallel (u_0,{\eta }_0){\parallel }_{{\mathcal{H}}_{\beta +1,\beta +1}^M})$. Now, for any fixed $n\in \mathbb{N}$, we find a unique local maximal solution ${\varphi }_n=(u_n,{\eta }_n)\in C^2([0,T_n];{\mathcal{H}}_{\beta +1,2}^M).$ Next, we integrate (115) and (116) with respect to time on $(0,t)$ and argue as in the proof of Theorem 1 to find the uniform bounds (78)–(84), (88), and (90). In order to obtain the required higher-order estimates, let us begin by labeling

$$\tilde{u}\left(t\right)={\partial }_{t}u\left(t\right),{\tilde{\eta }}^t={\partial }_t{\eta }^t,\tilde{\mu }\left(t\right)={\partial }_t\mu \left(t\right),$$

where we are also dropping the index n for the sake of simplicity. Then, $(\tilde{u},\tilde{\eta })$ solves the system

$$\begin{array}{cc}\hfill {\langle {\partial }_t\tilde{u},v\rangle }_{H^{-1}\times H^1}+{\int }_0^{\infty }\nu \left(s\right){\langle {\tilde{\eta }}^t\left(s\right),v\rangle }_{H^{-1}\times H^1}ds& =0,\hfill \end{array}$$

(122)

$$\begin{array}{cc}\hfill a_{E,\beta }(\tilde{u},\xi )+(F^{″}\left(u\right)\tilde{u},\xi )+\alpha ({\partial }_t\tilde{u},\xi )& ={\langle \mu ,\xi \rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}},\hfill \end{array}$$

(123)

$$\begin{array}{cc}\hfill {({\partial }_t{\tilde{\eta }}^t,\zeta )}_{{\mathcal{M}}_{-1}}-{(T_r{\tilde{\eta }}^t,\zeta )}_{{\mathcal{M}}_{-1}}& ={(\tilde{\mu },\zeta )}_{{\mathcal{M}}_0},\hfill \end{array}$$

(124)

for all $v\in H^1\left(\mathsf{\Omega }\right)$, $\xi \in W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$, and $\zeta \in {\mathcal{M}}_1$, with the initial conditions

$$\begin{array}{c}\hfill \tilde{u}\left(0\right)={\tilde{u}}_0\mathrm{and}{\tilde{\eta }}^0={\tilde{\eta }}_0.\end{array}$$

Let us now take $v=\tilde{\mu }$, $\xi ={\partial }_t\tilde{u}$, and $\zeta ={\tilde{\eta }}^t$ in (122)–(124), respectively. Summing the resulting identities together, we obtain, for all $t\in (0,T)$,

$$\begin{array}{c}\hfill \frac{1}{2}\frac{d}{dt}\left\{\parallel \tilde{u}{\parallel }_{W_0^{\beta ,2}}^2+{\parallel {\tilde{\eta }}^t\parallel }_{{\mathcal{M}}_{-1}}^2\right\}-{\int }_0^{\infty }{\nu }^{\prime }\left(s\right)\parallel {\tilde{\eta }}^t{\left(s\right)\parallel }_{H^{-1}}^{2}ds+\alpha {\parallel {\partial }_t\tilde{u}\parallel }^2=-(F^{″}\left(u\right)\tilde{u},{\partial }_t\tilde{u}).\end{array}$$

Here, we apply (K5) as well as (N4) with (83) and the embedding $W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\hookrightarrow L^2\left(\mathsf{\Omega }\right)$ to find

$$\begin{array}{cc}\hfill \frac{1}{2}\frac{d}{dt}\left\{\parallel \tilde{u}{\parallel }_{W_0^{\beta ,2}}^2+{\parallel {\tilde{\eta }}^t\parallel }_{{\mathcal{M}}_{-1}}^2\right\}+\lambda \parallel {\tilde{\eta }}^t{\parallel }_{{\mathcal{M}}_{-1}}^2+\alpha {\parallel {\partial }_t\tilde{u}\parallel }^2& \le C_{\alpha }\parallel \tilde{u}{\parallel }^2+\frac{\alpha }{2}\parallel {\partial }_t\tilde{u}\parallel \hfill \\ \hfill & \le C_{\alpha }\parallel \tilde{u}{\parallel }_{W_0^{\beta ,2}}^2+\frac{\alpha }{2}\parallel {\partial }_t\tilde{u}\parallel ,\hfill \end{array}$$

(125)

where $C_{\alpha }\sim {\alpha }^{-1}$ is a positive constant. Integrating (125) over $(0,t)$ produces

$$\begin{array}{cc}\hfill & \parallel \tilde{u}{\left(t\right)\parallel }_{W_0^{\beta ,2}}^2+{\parallel {\tilde{\eta }}^t\parallel }_{{\mathcal{M}}_{-1}}^2+{\int }_0^t\left(2\lambda \parallel {\tilde{\eta }}^{\tau }{\parallel }_{{\mathcal{M}}_{-1}}^2+\alpha {\parallel {\partial }_t\tilde{u}\left(\tau \right)\parallel }^2\right)d\tau \hfill \\ \hfill & \le \parallel \tilde{u}{\left(0\right)\parallel }_{W_0^{\beta ,2}}^2+\parallel {\tilde{\eta }}^0{\parallel }_{{\mathcal{M}}_{-1}}^2+C_{\alpha }{\int }_0^t{\parallel \tilde{u}\left(\tau \right)\parallel }_{W_0^{\beta ,2}}^{2}d\tau ,\hfill \end{array}$$

(126)

and an application of Grönwall’s (integral) inequality shows, for all $t\ge 0$,

$$\begin{array}{c}\hfill \parallel (\tilde{u}\left(t\right),{\tilde{\eta }}^t){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}\le Q(\parallel ({\tilde{u}}_0,{\tilde{\eta }}_0){\parallel }_{{\mathcal{H}}_{\beta ,0}^M})\end{array}$$

(127)

and

$$\begin{array}{c}\hfill \sqrt{\alpha }\parallel {\partial }_t\tilde{u}{\left(t\right)\parallel }_{L^2(\mathsf{\Omega }\times (0,T))}\le Q(\parallel ({\tilde{u}}_0,{\tilde{\eta }}_0){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}).\end{array}$$

(128)

Through (119)–(121), we find $\parallel ({\tilde{u}}_0,{\tilde{\eta }}_0){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}$ depends on

$${\int }_0^{\infty }\nu \left(s\right)\parallel {\eta }_0{\left(s\right)\parallel }_{W_0^{\beta ,2}}^{2}ds,\parallel A_N{\mu }_0{\parallel }_{{\mathcal{M}}_{-1}}\mathrm{and}{\parallel T_r{\eta }_0\parallel }_{{\mathcal{M}}_{-1}},$$

hence the assumption on the initial data is justified.

Furthermore, we now consider (67) and take $\zeta =A_N\overline{\mu }\left(t\right)$ where $\overline{\mu }=\mu -\langle \mu \rangle ,$ so that, with (79), (82), and (127), we obtain, for all $t\ge 0$ and for every $ϵ>0,$

$$\begin{array}{}\mathrm{(129)}& {\parallel \nabla \mu \parallel }^2& ={({\partial }_t{\eta }^t,\mu )}_{{\mathcal{M}}_0}-{(T_r{\eta }^t,\mu )}_{{\mathcal{M}}_0}\mathrm{(130)}& & ={\int }_0^{\infty }\nu \left(s\right)({\partial }_t{\eta }^t\left(s\right),\mu \left(t\right))ds-{\int }_0^{\infty }{\nu }^{\prime }\left(s\right)({\eta }^t\left(s\right),\mu \left(t\right))ds\mathrm{(131)}& & \le Cϵ\left.(\parallel {\partial }_t{\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}^2-{\int }_0^{\infty }{\nu }^{\prime }\left(s\right)\parallel {\eta }^t\left(s\right){\parallel }_{H^{-1}}^{2}ds\right)+ϵ{\parallel \nabla \mu \parallel }^2\mathrm{(132)}& & \le Cϵ\left.(1-{\int }_0^{\infty }{\nu }^{\prime }\left(s\right)\parallel {\eta }^t\left(s\right){\parallel }_{H^{-1}}^{2}ds\right)+ϵ{\parallel \nabla \mu \parallel }^2\mathrm{(133)}& & \le C_{ϵ}+ϵ{\parallel \nabla \mu \parallel }^2\end{array}$$

where $C_{ϵ}\sim {ϵ}^{-1}.$ Together (89) and (133) show us, for all $t\ge 0$

$${\parallel \mu \left(t\right)\parallel }_{H^1}\le C.$$

(134)

At this point we can reason as is in the proof of Theorem 1 to find that there is a solution $\varphi =(u,\eta )\in W^{1,\infty }(0,T;{\mathcal{H}}_{\beta ,0}^M)$ to Problem P satisfying (111) and (112). Additionally, thanks to (134), the condition (113) holds. It remains to show that

$$\varphi =(u,\eta )\in L^{\infty }\left(0,T;\left[W_0^{\beta +1,2}\left(\mathsf{\Omega }\right)\times L_{\nu }^2({\mathbb{R}}_{+};W_0^{\beta ,2}(\mathsf{\Omega }))\right]\right).$$

First, in light of (127), we multiply (16) by $A_{E,\beta }{\eta }^t$ in $L^2\left(\mathsf{\Omega }\right)$ which yields

$$\begin{array}{c}\hfill \parallel {\eta }^t{\parallel }_{L_{\nu }^2({\mathbb{R}}_{+};W_0^{\beta ,2}(\mathsf{\Omega }))}^2=-{\int }_0^{\infty }\nu \left(s\right)(A_{E,\beta }^{\frac{1}{2}}{\partial }_{t}u\left(t\right),A_{E,\beta }^{\frac{1}{2}}{\eta }^t\left(s\right))ds.\end{array}$$

Hence, $\eta \in L^{\infty }(0,T;L_{\nu }^2({\mathbb{R}}_{+};W_0^{\beta ,2}\left(\mathsf{\Omega }\right)))$. Next, we consider the identity (52) whereby we may now rely on the regularity properties of ${\partial }_{t}u$ and $\mu$. We take $\xi =A_N{\partial }_{t}u$ to produce

$$\begin{array}{c}\hfill \frac{1}{2}\frac{d}{dt}{|\parallel u\parallel |}_{W_0^{\beta +1,2}}^2+\langle F^{″}\left(u\right)\nabla u,\nabla u\rangle +\alpha {\parallel {\partial }_{t}u\parallel }_{H^1}^2=\langle \nabla \mu ,\nabla u\rangle .\end{array}$$

After applying (N1) and integrating the resulting differential inequality with respect to t over $(0,t)$, we obtain for all $t\ge 0,$

$$\begin{array}{c}\hfill {|\parallel u\left(t\right)\parallel |}_{W_0^{\beta +1,2}}^2+2{\int }_0^{\infty }\alpha \parallel {\partial }_t{u\left(\tau \right)\parallel }_{H^1}^{2}d\tau \le {|\parallel u\left(0\right)\parallel |}_{W_0^{\beta +1,2}}^2+Q(\parallel (u_0,{\eta }_0){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}).\end{array}$$

We now deduce

$$u\in L^{\infty }(0,T;W_0^{\beta +1,2}\left(\mathsf{\Omega }\right))\mathrm{and}\sqrt{\alpha }{\partial }_{t}u\in L^2(0,T;H^1\left(\mathsf{\Omega }\right)).$$

This completes the proof. □

The following proposition provides continuous dependence and uniqueness for the solutions constructed above.

Proposition 2.

Let the assumptions of Theorem 1 hold. Additionally, assume (N4) holds. Let $T>0$ and let ${\varphi }_i={(u_i,{\eta }_i)}^{tr}$, $i=1,2$, be two solutions to ProblemPon $(0,T)$ corresponding to the initial data ${\varphi }_{0i}={(u_{0i},{\eta }_{0i})}^{tr}\in {\mathcal{H}}_{\beta ,0}^M=W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\times {\mathcal{M}}_{-1}^{\left(0\right)}$, such that $F\left(u_{0i}\right)\in L^1\left(\mathsf{\Omega }\right)$, $i=1,2$. Then, for each $\alpha >0$, there is a positive constant $C_{\alpha }\sim {\alpha }^{-1}$ such that the following estimate holds, for any $t\in (0,T)$,

$$\begin{array}{cc}\hfill \parallel {\varphi }_1\left(t\right)-{\varphi }_2{\left(t\right)\parallel }_{{\mathcal{H}}_{\beta ,0}^M}^2& +{\int }_0^t\left(\alpha \parallel {\partial }_t{u}_1\left(\tau \right)-{\partial }_t{u}_2{\left(\tau \right)\parallel }^2+{\parallel {\eta }_1^{\tau }-{\eta }_2^{\tau }\parallel }_{L_{-{\nu }^{\prime }}^2({\mathbb{R}}_{+};H^{-1})}^2\right)d\tau \hfill \\ \hfill & \le e^{C_{\alpha }t}{\parallel {\varphi }_{01}-{\varphi }_{02}\parallel }_{{\mathcal{H}}_{\beta ,0}^M}^2.\hfill \end{array}$$

(135)

Proof. 

To begin, we assume $(u_{0i},{\eta }_{0i})$, $i=1,2$, satisfy the assumptions of Theorem 2 (recall, above we are assuming (N4) holds), and we work with the more regular solutions to obtain (135). For all $t\in [0,T]$, we then set

$$\begin{array}{c}\hfill \varphi \left(t\right):={\varphi }_1\left(t\right)-{\varphi }_2\left(t\right),u\left(t\right):=u_1\left(t\right)-u_2\left(t\right),{\eta }^t:={\eta }_1^t-{\eta }_2^t\mathrm{and}\mu :={\mu }_1-{\mu }_2\end{array}$$

where ${\varphi }_i\left(t\right)=(u_i\left(t\right),{\eta }_i^t)$ is a solution corresponding to $(u_{0i},{\eta }_{0i})$, $i=1,2.$ Then, formally, $\varphi =(u,\eta )$ solves the equations for all $v\in H^1\left(\mathsf{\Omega }\right),$$\xi \in W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\cap L^p\left(\mathsf{\Omega }\right)$, and $\zeta \in {\mathcal{M}}_1$:

$$\begin{array}{cc}\hfill {\langle {\partial }_{t}u,v\rangle }_{H^{-1}\times H^1}+{\int }_0^{\infty }\nu \left(s\right){\langle {\eta }^t\left(s\right),v\rangle }_{H^{-1}\times H^1}ds& =0,\hfill \end{array}$$

(136)

$$\begin{array}{cc}\hfill a_{E,\beta }(u,\xi )+{\langle F^{\prime }\left(u_1\right)-F^{\prime }\left(u_2\right),\xi \rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}}+\alpha {\langle {\partial }_{t}u,\xi \rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}}& ={\langle \mu ,\xi \rangle }_{W^{-\beta ,2}\times W_0^{\beta ,2}},\hfill \end{array}$$

(137)

$$\begin{array}{cc}\hfill {({\partial }_t{\eta }^t,\zeta )}_{{\mathcal{M}}_{-1}}-{(T_r{\eta }^t,\zeta )}_{{\mathcal{M}}_{-1}}& ={(\mu ,\zeta )}_{{\mathcal{M}}_0}\hfill \end{array}$$

(138)

with the initial data

$$\begin{array}{c}\hfill u\left(0\right)=u_{01}-u_{02},{\eta }^0={\eta }_{01}-{\eta }_{02}.\end{array}$$

In (136), we choose $v=\mu$ and in (137), we choose $\xi ={\partial }_{t}u$. Owing to Theorem 2, for each $t\in [0,T]$, these elements are in $H^1\left(\mathsf{\Omega }\right)$ and $W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$, respectively, then we sum the results to obtain

$$\begin{array}{c}\hfill (A_{E,\beta }u,{\partial }_{t}u)+(F^{\prime }\left(u_1\right)-F^{\prime }\left(u_2\right),{\partial }_{t}u)+\alpha {\parallel {\partial }_{t}u\parallel }^2+{\int }_0^{\infty }\nu \left(s\right)(\mu ,{\eta }^t\left(s\right))ds=0.\end{array}$$

(139)

Further, multiplying (138) by $A_N^{-1}{\eta }^t$ in ${\mathcal{M}}_0$, then adding the obtained relation to (139), we have

$$\begin{array}{c}\hfill \frac{1}{2}\frac{d}{dt}{\left\{\right|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+\parallel {\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}^2\}+\alpha \parallel {\partial }_t{u\parallel }^2-{\int }_0^{\infty }{\nu }^{\prime }\left(s\right){\parallel {\eta }^t\left(s\right)\parallel }_{H^{-1}}^{2}ds+(F^{\prime }\left(u_1\right)-F^{\prime }\left(u_2\right),{\partial }_{t}u)=0.\end{array}$$

(140)

Using Hölder’s inequality, (N4), Young’s inequality, and the embedding $L^{\infty }\left(\mathsf{\Omega }\right)\hookrightarrow W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$, we estimate the remaining product as

$$\begin{array}{cc}\hfill \left|\right(F^{\prime }\left(u_1\right)-F^{\prime }\left(u_2\right),{\partial }_{t}u\left)\right|& \le \parallel F^{\prime }\left(u_1\right)-F^{\prime }\left(u_2\right)\parallel \parallel {\partial }_{t}u\parallel \hfill \\ \hfill & \le C\parallel (1+|u_1{|}^{\rho -2}+|u_2{|}^{\rho -2})u\parallel \parallel {\partial }_{t}u\parallel \hfill \\ \hfill & \le C(1+\parallel u_1{\parallel }_{L^{2(\rho -2)}}^{\rho -2}+\parallel u_2{\parallel }_{L^{2(\rho -2)}}^{\rho -2}{)\parallel u\parallel }_{L^{\infty }}\parallel {\partial }_{t}u\parallel \hfill \\ \hfill & \le Q_{\alpha }(\parallel (u_{0i},{\eta }_{0i}){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}{\left)\right|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+\frac{\alpha }{2}{\parallel {\partial }_{t}u\parallel }^2,\hfill \end{array}$$

(141)

where the positive monotone increasing function $Q_{\alpha }(·)\sim {\alpha }^{-1}$ (we remind the reader $\parallel (u_{0i},{\eta }_{0i}){\parallel }_{{\mathcal{H}}_{\beta +1,\beta +1}^M}\le Q{\parallel (u_{0i},{\eta }_{0i})\parallel }_{{\mathcal{H}}_{\beta ,0}^M}$, for $i=1,2$ and the bounds on $u_1$ and $u_2$ follow from (61) and (62)). With (140) and (141), we obtain the following differential inequality which holds for almost all $t\in [0,T]$

$$\begin{array}{cc}\hfill \frac{d}{dt}& {\left\{\right|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+\parallel {\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}^2\}+\alpha \parallel {\partial }_t{u\parallel }^2+\parallel {\eta }^t{\parallel }_{L_{-{\nu }^{\prime }}^2({\mathbb{R}}_{+};H^{-1})}^2\le Q_{\alpha }(\parallel (u_{0i},{\eta }_{0i}){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}{\left)\right|\parallel u\parallel |}_{W_0^{\beta ,2}}^2\hfill \\ \hfill & \le Q_{\alpha }(\parallel (u_{0i},{\eta }_{0i}){\parallel }_{{\mathcal{H}}_{\beta ,0}^M})\left({|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+{\parallel {\eta }^t\parallel }_{{\mathcal{M}}_{-1}}^2\right).\hfill \end{array}$$

(142)

Applying a Grönwall inequality to (142), we obtain, for all $t\in [0,T],$

$$\begin{array}{cc}\hfill {|\parallel u\left(t\right)\parallel |}_{W_0^{\beta ,2}}^2+{\parallel {\eta }^t\parallel }_{{\mathcal{M}}_{-1}}^2& +{\int }_0^t\left(\alpha \parallel {\partial }_t{u\left(\tau \right)\parallel }^2+{\parallel {\eta }^{\tau }\parallel }_{L_{-{\nu }^{\prime }}^2({\mathbb{R}}_{+};H^{-1})}^2\right)d\tau \hfill \\ \hfill & \le e^{C_{\alpha }}\left({|\parallel u\left(0\right)\parallel |}_{W_0^{\beta ,2}}^2+{\parallel {\eta }^0\parallel }_{{\mathcal{M}}_{-1}}^2\right).\hfill \end{array}$$

(143)

This shows the claim (135) holds for the regular solutions. Since none of the above constants due to the above estimate actually depend on the assumptions of Theorem 2, then standard approximation arguments can be employed to obtain (135) for the weak solutions as well. □

Remark 9.

It is quite important to remark that when $N=3$, the uniqueness for the nonviscous problem (where $\alpha =0$) remains an open problem (indeed, see [36,46,47]).

We now formalize the semidynamical system generated by Problem P.

Corollary 2.

Let the assumptions of Theorem 1 be satisfied. Additionally, assume (N4) holds. We can define a strongly continuous semigroup of solution operators $\mathcal{S}={\left(\mathcal{S}\left(t\right)\right)}_{t\ge 0}$, for each $\alpha >0$ and $\beta \in (0,1)$,

$$\mathcal{S}\left(t\right):{\mathcal{X}}_{\beta ,0}^M\to {\mathcal{X}}_{\beta ,0}^M$$

by setting, for all $t\ge 0,$

$$\mathcal{S}\left(t\right){\varphi }_0:=\varphi \left(t\right)$$

where $\varphi \left(t\right)=(u\left(t\right),{\eta }^t)$ is the unique global weak solution to ProblemP. Furthermore, as a consequence of (135), the semigroup $\mathcal{S}\left(t\right):{\mathcal{X}}_{\beta ,0}^M\to {\mathcal{X}}_{\beta ,0}^M$ is Lipschitz continuous on ${\mathcal{X}}_{\beta ,0}^M$, uniformly in t on compact intervals.

4. Absorbing Sets and Global Attractors

We now give a dissipation estimate for Problem P from which we deduce the existence of a bounded absorbing set and an important uniform bound on the solutions of Problem P. The existence of an absorbing set is also used later to show that the semigroup of solution operators $\mathcal{S}$ admits a compact global attractor in the metric space ${\mathcal{X}}_{\beta ,0}^M$.

Lemma 1.

Let ${\varphi }_0={(u_0,{\eta }_0)}^{tr}\in {\mathcal{H}}_{\beta ,0}^M=W_0^{\beta ,2}\left(\mathsf{\Omega }\right)\times {\mathcal{M}}_{-1}^{\left(0\right)}$ for $\beta \in (\frac{N}{4},1),$$N=1,2,3,$ be such that $F\left(u_0\right)\in L^1\left(\mathsf{\Omega }\right).$ Assume (K1), (K3)–(K5), and (N1)–(N3) hold. Assume $\varphi ={(u,\eta )}^{tr}$ is a weak solution to ProblemP.There are positive constants ${\kappa }_1$ and C, each depending on Ω but independent of t, α, and ${\varphi }_0$, such that, for all $t\ge 0$, the following holds

$$\begin{array}{cc}\hfill & {\parallel \varphi \left(t\right)\parallel }_{{\mathcal{H}}_{\beta ,0}^M}^2+{\int }_t^{t+1}\alpha \parallel {\partial }_t{u\left(\tau \right)\parallel }^{2}d\tau \le Q(\parallel {\varphi }_0{\parallel }_{{\mathcal{H}}_{\beta ,0}^M})e^{-{\kappa }_{1}t}+C,\hfill \end{array}$$

(144)

for some monotonically increasing function Q independent of t and α.

Proof. 

The idea of the proof is from [26]. We give a formal calculation that can be justified by a suitable Faedo–Galerkin approximation based on the proof of Theorem 1 above. To begin, define the functional, for all $t\ge 0,$

$$\begin{array}{c}\hfill \mathcal{Y}\left(t\right):=\mathcal{E}\left(t\right)+{ϵ}^{\alpha }{\parallel u\left(t\right)\parallel }^2-2ϵ{\int }_0^{\infty }\nu \left(s\right)\left(u\left(t\right),A_N^{-1}{\eta }^t\left(s\right)\right)ds,\end{array}$$

(145)

where ${ϵ}^{\in }(0,\lambda )$ will be chosen sufficiently small later. From (16)–(18), we find

$$\begin{array}{cc}\hfill -\frac{d}{dt}& {\int }_0^{\infty }\nu \left(s\right)(u,A_N^{-1}{\eta }^t\left(s\right))ds\hfill \\ \hfill & =\parallel {\partial }_t{u\parallel }_{H^{-1}}^2-{\int }_0^{\infty }\nu \left(s\right)(u,A_N^{-1}{\partial }_t{\eta }^t\left(s\right))ds\hfill \\ \hfill & =\parallel {\partial }_t{u\parallel }_{H^{-1}}^2-{\int }_0^{\infty }{\nu }^{\prime }\left(s\right)(u,A_N^{-1}{\eta }^t\left(s\right))ds-{\int }_0^{\infty }\nu \left(s\right)(u,\mu )ds\hfill \\ \hfill & =\parallel {\partial }_t{u\parallel }_{H^{-1}}^2-{\int }_0^{\infty }{\nu }^{\prime }\left(s\right)(u,A_N^{-1}{\eta }^t\left(s\right))ds-\frac{\alpha }{2}\frac{d}{dt}{\parallel u\parallel }^2-{|\parallel u\parallel |}_{W_0^{\beta ,2}}^2-(F^{\prime }\left(u\right),u).\hfill \end{array}$$

(146)

Differentiating $\mathcal{Y}$ with respect to t while keeping in mind (73), (74) (without the index n), and (146), we find

$$\begin{array}{c}\hfill \frac{d}{dt}\mathcal{Y}+{ϵ}_0\mathcal{Y}-2{\int }_0^{\infty }{\nu }^{\prime }\left(s\right){\parallel {\eta }^t\left(s\right)\parallel }_{H^{-1}}^{2}ds=h\left(t\right),\end{array}$$

(147)

for ${ϵ}_0\in (0,ϵ)$ where

$$\begin{array}{cc}\hfill h\left(t\right)=& -2\alpha \parallel {\partial }_{t}u\left(t\right){{\parallel }^2+2ϵ\parallel {\partial }_{t}u\left(t\right)\parallel }_{H^{-1}}^2-2ϵ{\int }_0^{\infty }{\nu }^{\prime }\left(s\right)(u\left(t\right),A_N^{-1}{\eta }^t\left(s\right))ds\hfill \\ \hfill & -2{ϵ}_0(F^{\prime }\left(u\left(t\right)\right)u\left(t\right)-F\left(u\left(t\right)\right),1)-2(ϵ-{ϵ}_0)(F^{\prime }\left(u\left(t\right)\right),u\left(t\right))+{ϵ}_0{\parallel {\eta }^t\parallel }_{{\mathcal{M}}_{-1}}^2\hfill \\ \hfill & -(2ϵ-{ϵ}_0)|\parallel u\left(t\right){\parallel |}_{W_0^{\beta ,2}}^2+{ϵ}_0{ϵ}^{\alpha }{\parallel u\left(t\right)\parallel }^2-2{ϵ}_0ϵ{\int }_0^{\infty }\nu \left(s\right)(u\left(t\right),A_N^{-1}{\eta }^t\left(s\right))ds+{ϵ}_{0}C.\hfill \end{array}$$

(148)

From (42) and (43) (with $M=0$), it follows that

$$\begin{array}{cc}\hfill & -2{ϵ}_0(F^{\prime }\left(u\left(t\right)\right)u\left(t\right)-F\left(u\left(t\right)\right),1)-2(ϵ-{ϵ}_0)(F^{\prime }\left(u\left(t\right)\right),u\left(t\right))\hfill \\ \hfill & \le -(ϵ-{ϵ}_0)\left(\right|F\left(u\right)|,1)+{ϵ}_0{C|\parallel u\parallel |}_{W_0^{\beta ,2}}^2.\hfill \end{array}$$

(149)

Next, using assumption (K4) and the embeddings $H^{-1}\left(\mathsf{\Omega }\right)\hookleftarrow L^2\left(\mathsf{\Omega }\right)\hookleftarrow W_0^{\beta ,2}\left(\mathsf{\Omega }\right)$, we find

$$\begin{array}{cc}\hfill -2ϵ{\int }_0^{\infty }{\nu }^{\prime }\left(s\right)(u,A_N^{-1}{\eta }^t\left(s\right))ds& =-2ϵ{\int }_0^{\infty }{\nu }^{\prime }\left(s\right)(A_N^{-1/2}u,A_N^{-1/2}{\eta }^t\left(s\right))ds\hfill \\ \hfill & \le -ϵ{\int }_0^{\infty }{\nu }^{\prime }\left(s\right)\left(\frac{1}{{\nu }_0}{|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+C{\nu }_0{\parallel {\eta }^t\left(s\right)\parallel }_{H^{-1}}^2\right)ds\hfill \\ \hfill & \le ϵ|{\parallel u\parallel |}_{W_0^{\beta ,2}}^2-ϵC{\int }_0^{\infty }{\nu }^{\prime }\left(s\right){\parallel {\eta }^t\left(s\right)\parallel }_{H^{-1}}^{2}ds,\hfill \end{array}$$

(150)

and, with (K3) and (87) (without the index n),

$$\begin{array}{cc}\hfill -2{ϵ}_0ϵ{\int }_0^{\infty }\nu \left(s\right)(u,A_N^{-1}{\eta }^t\left(s\right))ds& \le {ϵ}_0ϵC{|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+{ϵ}_0ϵ\parallel {\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}^2.\hfill \end{array}$$

(151)

Together, (148)–(151) make the following estimate

$$\begin{array}{cc}\hfill h\le & -2\alpha \parallel {\partial }_t{u\parallel }^2+2{ϵ}^{\parallel }{\partial }_t{u\parallel }_{H^{-1}}^2-(ϵ-{ϵ}_0(1+C+{ϵ}^{\alpha }C)){|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+2{ϵ}_0{\parallel {\eta }^t\parallel }_{{\mathcal{M}}_{-1}}^2\hfill \\ \hfill & -ϵC{\int }_0^{\infty }{\nu }^{\prime }\left(s\right){\parallel {\eta }^t\left(s\right)\parallel }_{H^{-1}}^{2}ds+C.\hfill \end{array}$$

(152)

Here, we employ assumption (K5) so that from (147) and (152), we are able to fix ${ϵ}^{\in }(0,\lambda )$ and ${ϵ}_0\in (0,ϵ)$ sufficiently small to, in turn, find positive constants ${ϵ}_1,{ϵ}_2,{ϵ}_3$ so that there holds

$$\begin{array}{c}\hfill \frac{d}{dt}\mathcal{Y}+{ϵ}_1\mathcal{Y}+2\parallel {\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}^2+{ϵ}_2\alpha \parallel {\partial }_t{u\parallel }^2+{ϵ}_3{|\parallel u\parallel |}_{W_0^{\beta ,2}}^2\le C.\end{array}$$

(153)

It is important to note that C on the right-hand side of (150) is independent of t and ${\varphi }_0.$ One can readily show (see (73), (76)–(77)) that there holds, for all $t\ge 0,$

$$\begin{array}{c}\hfill C_1{\parallel \varphi \left(t\right)\parallel }_{{\mathcal{H}}_{\beta ,0}^M}^2-C_2\le \mathcal{Y}\left(t\right)\le Q(\parallel {\varphi }_0{\parallel }_{{\mathcal{H}}_{\beta ,0}^M}),\end{array}$$

(154)

for some positive constants $C_1,C_2$, and for some monotone nondecreasing function Q independent of t. Finally, by applying a Grönwall type inequality to (153) (see, e.g., [34] (Lemma 2.5)), then integrating the result and applying (154) yield the claim (144). This finishes the proof. □

We immediately deduce the existence of a bounded absorbing set from Lemma 1.

Proposition 3.

Let the assumptions of Lemma 1 hold. Additionally, assume (N4) holds. Then, there exists $R_0>0$, independent of t and ${\varphi }_0$, such that $\mathcal{S}\left(t\right)$ possesses an absorbing ball ${\mathcal{B}}_{\beta ,0}^M\left(R_0\right)\subset {\mathcal{H}}_{\beta ,0}^M,$ bounded in ${\mathcal{H}}_{\beta ,0}^M$. Precisely, for any bounded subset $B\subset {\mathcal{H}}_{\beta ,0}^M,$ there exists $t_0=t_0\left(B\right)>0$ such that $\mathcal{S}\left(t\right)B\subset {\mathcal{B}}_{\beta ,0}^M\left(R_0\right)$, for all $t\ge t_0.$ Moreover, for every $R>0$, there exists $C_{*}=C_{*}\left(R\right)\ge 0,$ such that, for any ${\varphi }_0\in {\mathcal{B}}_{\beta ,0}^M\left(R\right),$

$$\begin{array}{c}\hfill \underset{t\ge 0}{sup}\parallel \mathcal{S}\left(t\right){\varphi }_0{\parallel }_{{\mathcal{H}}_{\beta ,0}^M}+{\int }_0^{\infty }{\parallel {\partial }_{t}u\left(\tau \right)\parallel }^{2}d\tau \le C_{*},\end{array}$$

(155)

where ${\mathcal{B}}_{\beta ,0}^M\left(R\right)$ denotes the ball in ${\mathcal{H}}_{\beta ,0}^M$ of radius R, centered at $\mathbf{0}$.

Throughout the remainder of the article, we simply write ${\mathcal{B}}_{\beta ,0}^M$ in place of ${\mathcal{B}}_{\beta ,0}^M\left(R_0\right)$ to denote the bounded absorbing set admitted by the semigroup of solution operators $\mathcal{S}\left(t\right)$.

For the rest of this section, our aim is to prove the following.

Theorem 3.

Let the assumptions of Lemma 1 hold. Additionally, assume (N4) holds. The dynamical system$({\mathcal{X}}_{\beta ,0}^M,\mathcal{S}\left(t\right))$(see Corollary 2) possesses a connected global attractor${\mathcal{A}}_{\beta ,0}^M$in${\mathcal{H}}_{\beta ,0}^M.$Precisely:

1 

For each$t\ge 0$, $\mathcal{S}\left(t\right){\mathcal{A}}_{\beta ,0}^M={\mathcal{A}}_{\beta ,0}^M$;

2 

For every nonempty bounded subset B of${\mathcal{H}}_{\beta ,0}^M$,

$$\underset{t\to \infty }{lim}{\mathrm{dist}}_{{\mathcal{H}}_{\beta ,0}^M}(\mathcal{S}\left(t\right)B,{\mathcal{A}}_{\beta ,0}^M):=\underset{t\to \infty }{lim}\underset{\zeta \in B}{sup}\underset{\xi \in {\mathcal{A}}_{\beta ,0}^M}{inf}{\parallel \mathcal{S}\left(t\right)\zeta -\xi \parallel }_{{\mathcal{H}}_{\beta ,0}^M}=0.$$

Additionally:

3 

The global attractor is the unique maximal compact invariant subset in ${\mathcal{H}}_{\beta ,0}^M$ given by

$${\mathcal{A}}_{\beta ,0}^M:=\omega \left({\mathcal{B}}_{\beta ,0}^M\right):=\bigcap _{s\ge 0}{\overline{\bigcup _{t\ge s}\mathcal{S}\left(t\right){\mathcal{B}}_{\beta ,0}^M}}^{{\mathcal{H}}_{\beta ,0}^M}.$$

Furthermore:

4 

The global attractor ${\mathcal{A}}_{\beta ,0}^M$ is connected and given by the union of the unstable manifolds connecting the equilibria of $\mathcal{S}\left(t\right)$.

5 

For each ${\zeta }_0={({\varphi }_0,{\theta }_0)}^{tr}\in {\mathcal{H}}_{\beta ,0}^M$, the set $\omega \left({\zeta }_0\right)$ is a connected compact invariant set, consisting of the fixed points of $\mathcal{S}\left(t\right).$

With the existence of a bounded absorbing set ${\mathcal{B}}_{\beta ,0}^M$ (in Lemma 1), the existence of a global attractor now depends on the precompactness of the semigroup of solution operators $\mathcal{S}$. To this end we show there is a $t_{*}>0$ such that the map $\mathcal{S}\left(t_{*}\right)$ is a so-called $\alpha$-contraction on ${\mathcal{B}}_{\beta ,0}^M$; that is, there exist a time $t_{*}>0$, a constant $0<\kappa <1$, and a precompact pseudometric $M_{*}$ on ${\mathcal{B}}_{\beta ,0}^M$ such that, for all ${\varphi }_{01},{\varphi }_{02}\in {\mathcal{B}}_{\beta ,0}^M$,

$$\parallel \mathcal{S}\left(t_{*}\right){\varphi }_{01}-\mathcal{S}\left(t_{*}\right){\varphi }_{02}{\parallel }_{{\mathcal{H}}_{\beta ,0}^M}\le \kappa {\parallel {\varphi }_{01}-{\varphi }_{02}\parallel }_{{\mathcal{H}}_{\beta ,0}^M}+M_{*}({\varphi }_{01},{\varphi }_{02}).$$

(156)

Such a contraction is commonly used in connection with phase-field-type equations as an alternative to establish the precompactness of a semigroup; for some particular recent results see [16,48,49].

Lemma 2.

Under the assumptions of Proposition 2 where ${\varphi }_{01},{\varphi }_{02}\in {\mathcal{B}}_{\beta ,0}^M$, there are positive constants ${\kappa }_2,C_1$, and $C_{2\alpha }\sim {\alpha }^{-1}$, each depending on $\mathsf{\Omega }$ but independent of t and ${\varphi }_{01},{\varphi }_{02}$, such that, for all $t\ge 0,$

$$\begin{array}{cc}\hfill \parallel {\varphi }_1\left(t\right)-{\varphi }_2{\left(t\right)\parallel }_{{\mathcal{H}}_{\beta ,0}^M}^2& \le C_1{e}^{-{\kappa }_{2}t}{\parallel {\varphi }_1\left(0\right)-{\varphi }_2\left(0\right)\parallel }_{{\mathcal{H}}_{\beta ,0}^M}^2\hfill \\ \hfill & +C_{2\alpha }{\int }_0^t\left(\parallel \nabla {\mu }_1\left(\tau \right)-\nabla {\mu }_2{\left(\tau \right)\parallel }^2+{\parallel u_1\left(\tau \right)-u_2\left(\tau \right)\parallel }^2\right)d\tau .\hfill \end{array}$$

(157)

Proof. 

The proof is based on the proof of Proposition 2. We begin by recovering (140) by multiplying (136) and (137) by $\mu$ and ${\partial }_{t}u$, respectively, in $L^2\left(\mathsf{\Omega }\right)$, and multiplying (138) by $A_N^{-1}{\eta }^t$ in ${\mathcal{M}}_0$, then adding the obtained relations together to find

$$\begin{array}{c}\hfill \frac{1}{2}\frac{d}{dt}{\left\{\right|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+\parallel {\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}^2\}+\alpha \parallel {\partial }_t{u\parallel }^2-{\int }_0^{\infty }{v}^{\prime }\left(s\right){\parallel {\eta }^t\left(s\right)\parallel }_{H^{-1}}^{2}ds+(F^{\prime }\left(u_1\right)-F^{\prime }\left(u_2\right),{\partial }_{t}u)=0.\end{array}$$

(158)

Recall ${\varphi }_1=(u_1,{\eta }_1)$, ${\varphi }_2=(u_2,{\eta }_2)$ are the unique weak solutions corresponding to the initial data ${\varphi }_{01}$ and ${\varphi }_{02}$, respectively; also, $u=u_1-u_2$ and ${\eta }^t={\eta }_1^t-{\eta }_2^t$ formally satisfy (136) and (137). Applying Assumption (K5) and the estimate based on (N4),

$$\begin{array}{}& \left|\right(F^{\prime }\left(u_1\right)-F^{\prime }\left(u_2\right),{\partial }_{t}u\left)\right|& \le \parallel F^{\prime }\left(u_1\right)-F^{\prime }\left(u_2\right)\parallel \parallel {\partial }_{t}u\parallel & & \le C\parallel (1+|u_1{|}^{\rho -2}+|u_2{|}^{\rho -2})u\parallel \parallel {\partial }_{t}u\parallel & & \le C(1+\parallel u_1{\parallel }_{L^{2(\rho -2)}}^{\rho -2}+\parallel u_2{\parallel }_{L^{2(\rho -2)}}^{\rho -2}{)\parallel u\parallel }_{L^{\infty }}\parallel {\partial }_{t}u\parallel \mathrm{(159)}& & \le Q_{\alpha }(\parallel (u_{0i},{\eta }_{0i}){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}{\left)\right|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+\frac{\alpha }{2}{\parallel {\partial }_{t}u\parallel }^2\mathrm{(160)}& & \le Q_{\alpha }(\parallel (u_{0i},{\eta }_{0i}){\parallel }_{{\mathcal{H}}_{\beta ,0}^M})+\frac{\alpha }{2}{\parallel {\partial }_{t}u\parallel }^2,\end{array}$$

where the positive monotone increasing function $Q_{\alpha }(·)\sim {\alpha }^{-1},$ and we find the differential inequality

$$\begin{array}{c}\hfill \frac{1}{2}\frac{d}{dt}{\left\{\right|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+\parallel {\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}^2\}+\frac{\alpha }{2}\parallel {\partial }_t{u\parallel }^2+\lambda \parallel {\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}^2\le Q_{\alpha }(\parallel (u_{0i},{\eta }_{0i}){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}).\end{array}$$

(161)

In addition, we now multiply (137) by u in $L^2\left(\mathsf{\Omega }\right)$ to obtain

$$\begin{array}{c}\hfill {|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+(F^{\prime }\left(u_1\right)-F^{\prime }\left(u_2\right),u)+\frac{\alpha }{2}\frac{d}{dt}{\parallel u\parallel }^2=(\mu ,u).\end{array}$$

(162)

Estimating the first product above using (N1) yields

$$\begin{array}{cc}\hfill (F^{\prime }\left(u_1\right)-F^{\prime }\left(u_2\right),u)& \ge -c_F{\parallel u\parallel }^2.\hfill \end{array}$$

(163)

We also estimate with Young’s inequality

$$\begin{array}{cc}\hfill (\mu ,u)& \le \frac{1}{2}{\parallel \mu \parallel }^2+\frac{1}{2}{\parallel u\parallel }^2.\hfill \end{array}$$

(164)

Combining (161)–(164) yields

$$\begin{array}{cc}\hfill \frac{1}{2}\frac{d}{dt}& \left\{{|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+\parallel {\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}^2+\frac{\alpha }{2}{\parallel u\parallel }^2\right\}+\frac{\alpha }{2}\parallel {\partial }_t{u\parallel }^2+{|\parallel u\parallel |}_{W_0^{\beta ,2}}^2+\lambda {\parallel {\eta }^t\parallel }_{{\mathcal{M}}_{-1}}^2\hfill \\ \hfill & \le \frac{1}{2}{\parallel \mu \parallel }^2+Q_{\alpha }(\parallel (u_{0i},{\eta }_{0i}){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}{\left)\right|\parallel u\parallel |}_{W_0^{\beta ,2}}^2.\hfill \end{array}$$

(165)

Then, adding $\frac{\alpha }{2}{\parallel u\parallel }^2$ to each side of (165), we find

$$\begin{array}{c}\hfill \frac{d}{dt}\mathcal{N}+c\mathcal{N}+\alpha \parallel {\partial }_t{u\parallel }^2\le {\parallel \mu \parallel }^2+Q_{\alpha }(\parallel (u_{0i},{\eta }_{0i}){\parallel }_{{\mathcal{H}}_{\beta ,0}^M}),\end{array}$$

(166)

where $c=min\{2,2\lambda ,\alpha \}$ and

$$\begin{array}{c}\hfill \mathcal{N}\left(t\right):={|\parallel u\left(t\right)\parallel |}_{W_0^{\beta ,2}}^2+\parallel {\eta }^t{\parallel }_{{\mathcal{M}}_{-1}}^2+\frac{\alpha }{2}{\parallel u\left(t\right)\parallel }^2.\end{array}$$

(167)

Applying Grönwall’s inequality to (166) after omitting the term $\alpha \parallel {\partial }_t{u\parallel }^2$, we obtain the claim (157). □

Consequently, we deduce the following precompactness result for the semigroup $\mathcal{S}$.

Proposition 4.

Let the assumptions of Lemma 2 hold. There is $t_{*}>0$ such that the operator $\mathcal{S}\left(t_{*}\right)$ is a strict contraction up to the precompact pseudometric on ${\mathcal{B}}_{\beta ,0}^M$, in the sense of (156), where

$$\begin{array}{cc}\hfill M_{*}({\varphi }_{01},{\varphi }_{02})& :=C_{2\alpha }{\left({\int }_0^{t_{*}}\left(\parallel \nabla {\mu }_1\left(\tau \right)-\nabla {\mu }_2{\left(\tau \right)\parallel }^2+{\parallel u_1\left(\tau \right)-u_2\left(\tau \right)\parallel }^2\right)d\tau \right)}^{1/2},\hfill \end{array}$$

(168)

with $C_{\alpha }\sim {\alpha }^{-1}$. Furthermore, $\mathcal{S}$ is precompact on ${\mathcal{B}}_{\beta ,0}^M$.

Proof. 

Naturally, we follow from the conclusion of Lemma 2. Clearly, there is a $t_{*}>0$ so that $C_1{e}^{-{\kappa }_2{t}_{*}/2}<1.$ Thus, the operator $\mathcal{S}\left(t_{*}\right)$ is a strict contraction up to the pseudometric $M_{*}$ defined by (168). The pseudometric $M_{*}$ is precompact thanks to the Aubin–Lions compact embedding (99). This completes the proof. □

Proof of Theorem 3.

The precompactness of the solution operators $\mathcal{S}$ follows via the method of precompact pseudometrics (see Proposition 4). With the existence of a bounded absorbing set ${\mathcal{B}}_{\beta ,0}^M$ in ${\mathcal{H}}_{\beta ,0}^M$ (Lemma 1), the existence of a global attractor in ${\mathcal{H}}_{\beta ,0}^M$ is well-known and can be found in [50,51] for example. Additional characteristics of the attractor follow thanks to the gradient structure of Problem P (Remark 8). In particular, the first three claims in the statement of Theorem 3 are a direct result of the existence of an absorbing set, a Lyapunov functional $\mathcal{E}$, and the fact that the system $({\mathcal{X}}_{\beta ,0}^M,\mathcal{S}\left(t\right),\mathcal{E})$ is a gradient. The fourth property is a direct result of [51] (Theorem VII.4.1), and the fifth follows from [52] (Theorem 6.3.2). This concludes the proof. □