Nonlocal-to-Local Convergence of Cahn–Hilliard Equations: Neumann Boundary Conditions and Viscosity Terms

We consider a class of nonlocal viscous Cahn–Hilliard equations with Neumann boundary conditions for the chemical potential. The double-well potential is allowed to be singular (e.g. of logarithmic type), while the singularity of the convolution kernel does not fall in any available existence theory under Neumann boundary conditions. We prove well-posedness for the nonlocal equation in a suitable variational sense. Secondly, we show that the solutions to the nonlocal equation converge to the corresponding solutions to the local equation, as the convolution kernels approximate a Dirac delta. The asymptotic behaviour is analyzed by means of monotone analysis and Gamma convergence results, both when the limiting local Cahn–Hilliard equation is of viscous type and of pure type.


Introduction
The aim of the present paper is to study the well-posedness and the asymptotic behaviour as ε 0 of a family of nonlocal viscous Cahn-Hilliard equations with Neumann boundary conditions in the following form: where is a smooth bounded domain in R d (d = 2, 3), T > 0 is a fixed final time, and represents the derivative a double-well potential. Moreover, ε > 0 is a fixed parameter, τ ε > 0 is a positive viscosity coefficient, K ε : × → R is a suitable symmetric convolution kernel, and g ε represents a distributed forcing term. The variables u ε and μ ε are referred to as "order parameter" and "chemical potential", respectively.
The evolution problem (1.1)-(1.4) is related to the gradient flow (in the H −1metric) associated to a nonlocal free energy functional of the form Indeed, the contributions (K ε * 1)u ε − K ε * u ε + (u ε ) in the definition of the chemical potential are obtained exactly from the (sub)differentiation of the functional (1.5). The extra term τ ε ∂ t u ε represents on the other side a viscosity regularization, acting on the dissipation of the system. The analysis of nonlocal models dates back to the early 90's, when Giacomin and Lebowitz investigated, in their seminal paper [36], a hydrodynamic limit of a microscopic model for a d-dimensional lattice gas evolving via a Poisson nearestneighbor process. In that work, the authors derived a free energy functional in nonlocal form (1.5), and proposed the corresponding gradient flow to model phasechange in binary alloys. The viscous regularization in the definition of the chemical potential was originally introduced in the context of the local Cahn-Hilliard equation by Novick-Cohen in [50]. The mathematical literature on the nonlocal Cahn-Hilliard equation is widely developed: we can mention, among many others, the contributions [2,5,27,28,39] and the references therein.
The rapidly growing attention to the nonlocal Cahn-Hilliard equation is due on the one hand to its microscopic justification, and on the other hand to its connection with the corresponding local model. Indeed, at least in a formal way, the nonlocal dynamics approach the local ones when the family of interaction kernels (K ε ) ε concentrates around the origin. The main issue we assess in this paper is the asymptotic convergence of solutions to the nonlocal system (1.1)- (1.4) to the corresponding local one, as the data (g ε ) ε approximate a new source g and the coefficients τ ε converge to a certain new viscosity parameter τ . The local form of the limiting Cahn-Hilliard equation reads is (1.7) ∂ n u = 0 and ∂ n μ = 0 on(0, T ) × ∂ , (1.8) u(0) = u 0 in , (1.9) where τ 0 is the limiting viscosity parameter, which is allowed to vanish. The choices τ > 0 and τ = 0 correspond to the viscous case and pure case, respectively.
As its nonlocal counterpart, the local Cahn-Hilliard equation is related to the gradient flow in the H −1 metric of the Ginzburg-Landau free energy functional (1.10) in the sense that the contribution − u + (u) results from the subdifferentiation of E. Again, the viscosity term τ ∂ t u acts on the dissipation of the system: if τ = 0, one recovers the so-called pure Cahn-Hilliard equation, while if τ > 0 one obtains the viscous Cahn-Hilliard equation. In our analysis, the nonlocal viscosity coefficients (τ ε ) ε are assumed to be strictly positive, while the local coefficient τ is allowed to vanish. The local Cahn-Hilliard equation was first proposed in [9] in relation to phasechange in metallic alloys and to spinodal decomposition (see [44]). Nowadays, the model is a widely used in various contexts such as diffuse interface modelling in physics and biology, with several applications to tumor growth dynamics, image processing, and population dynamics. From the mathematical point of view, the local Cahn-Hilliard equation has been studied thoroughly in the last decades, also in much more complex settings. We mention, among many others, the works [11][12][13]15,16,18,37,43] on well-posedness also under more general dynamic boundary conditions. Some studies on nonlinear viscosity contributions have been proposed in [6,49,56]. We also recall the contributions [14,19,20,23,40] dealing with optimal control problems, as well as [17,22,38] on asymptotics. The local Cahn-Hilliard equation has also been widely studied recently in connection to diffuse-interface models for fluid-dynamics: we refer to [1,3,10,29,30] and the references therein.
As has already been mentioned, the behaviour of the nonlocal Cahn-Hilliard equation "approaches" the one of the local equation when the family of convolution kernels is sufficiently peaked around 0. The study of nonlocal-to-local convergence of energy functionals in relation to Sobolev spaces theory had been carried out originally by by Bourgain et al. [7,8], and by Mazy'a and Shaposhnikova [45,46]. This asymptotic analysis was also extended by Ponce [51,52], with studies on Gamma convergence and nonlocal Poincaré-type inequalities. A first criterion for the convergence of gradient flows from the Gamma-convergence of the respective energies was given by Sandier and Serfaty [55] in a abstract setting and for smooth energies, with applications to Ginzburg-Landau functionals (see also [42,53,57] for further details in this direction).
In particular, the above-mentioned results [51,52] provide the pointwise convergence as soon as the convolution kernels (K ε ) ε are chosen as where (ρ ε ) ε is a suitable family of mollifiers converging to a Dirac delta. Building upon these variational convergences, in a previous contribution of ours [24] we rigorously derived some nonlocal-to-local asymptotics of solutions to Cahn-Hilliard equations in the setting of periodic boundary conditions and with no viscosity effects. The periodic setting adopted in [24] was fundamental to overcome the singular behaviour of the convolution kernel (1.11). Indeed, kernels in the form (1.11) do not possess any W 1,1 regularity (see for example [21,Remark 1]), which is the usual minimum requirement in the whole literature on nonlocal Cahn-Hilliard systems. This resulted in the impossibility of framing the nonlocal problem in any available existence theory, and required an ad-hoc analysis. In this direction, the arguments strongly relied on the assumption of periodic boundary conditions. The results in [24] (see also [47] for a simpler case) are very satisfactory since they provide a novel contribution in the direction of local asymptotics of Cahn-Hilliard equations. Nevertheless, the most natural choice of boundary conditions in phase-field modelling if of no-flux type. Consequently, it is crucial in this direction to generalize the periodic framework to other settings more suited for applications. The nonlocal-to-local convergence of pure Cahn-Hilliard equations with Neumann boundary conditions was, to the authors' knowledge, still an open problem. The main novelty of the present paper is to finally extend some rigorous nonlocal-tolocal convergence results for Cahn-Hilliard equations to the case of homogeneous Neumann boundary conditions. Let us briefly describe now the main difficulties arising in the case of Neumann boundary conditions. The first hurdle has been already anticipated and concerns the regularity of the convolution kernel. Indeed, in the form (1.11) the kernel K ε is not W 1,1 , and not even L 1 in dimension d = 2. This results in the necessity of rigorously formulate the nonlocal problem without relying on any available existence theory. The main idea here is that even if the convolution operator ϕ → K ε * ϕ may be ill-defined under (1.11), the nonlocal operator B ε : ϕ → (K ε * 1)ϕ − (K ε * ϕ) appearing in the equation (1.2) can be rigorously defined instead.
The second main problem consists in the (im)possibility of proving space regularity for the solutions to the nonlocal equation (i.e. when ε > 0 is fixed). If the convolution kernel is W 1,1 this follows directly from the properties of the convolution, i.e. formally shifting the gradient operator on the kernel as ∇(K ε * u ε ) = (∇ K ε ) * u ε . However, for singular kernels as in (1.11) this procedure fails. Under periodic boundary conditions (i.e. working on the d-dimensional flat torus) the main idea to overcome this problem was to use a certain integration-by-parts formula, which hinges in turn on some compatibility conditions between the convolution operator and the Laplace operator. More specifically, in [24] the periodic setting allowed to prove a (formal) relation in the form ∇(K ε * u ε ) = K ε * ∇u ε , from which one could deduce H 1 -regularity of the nonlocal solutions. Nevertheless, under Neumann boundary conditions (i.e. working on a bounded domain ⊂ R d ), in order to prove an analogous compatibility relation one is forced to extend the nonlocal solution u ε to 0 outside . Clearly, H 1 -regularity in does not imply H 1 -regularity on the whole R d for such extension. This gives rise to several extra boundary contribution terms which blow up as the approximating parameter vanishes.
The main consequence is that in the case of Neumann boundary conditions one loses any H 1 -estimate on the nonlocal solutions. It follows that the natural variational setting to frame the nonlocal problem (1.1)-(1.4) is not the usual one given by the triple (H 1 ( ), L 2 ( ), H 1 ( ) * ), but instead an abstract one (V ε , L 2 ( ), V * ε ), depending on ε, where V ε represents, roughly speaking, the domain of the nonlocal energy contribution in (1.5). As the inclusion V ε → L 2 ( ) is not compact, one loses any reasonable compactness property on the approximated solutions in order to pass to the limit in the nonlinearity. This issue is overcome by the introduction of the viscosity term τ ε ∂ t u ε . Indeed, if τ ε is strictly positive one can show "by hand" a strong convergence in L 2 ( ) for some regularized solutions, even without relying on any H 1 estimates.
The third main problem concerns the boundary conditions of Neumann type for u in the limiting local problem. Indeed, while the nonlocal system is of order 2 in space, hence it only needs one boundary condition (for the chemical potential), the limiting local equation is of order 4 in space and requires two boundary conditions instead: one for μ and one for u. One of the major point is to understand which is the natural extra boundary condition for u, and how this one emerges when ε 0. It is clear that the Neumann boundary condition for the chemical potential is preserved by the local asymptotics. On the other hand, the scenario for u is more subtle: the answer is implicitly given by studying the Gamma convergence of the nonlocal energies. Indeed, in [52] Ponce proved a Gamma convergence result in the form whenever ϕ ε → ϕ in L 2 ( ). Note that the limiting energy contribution on the righthand side is the potential associated to the negative Laplacian with homogeneous Neumann boundary conditions. Hence, this implicitly reveals that the "correct" choice of boundary condition arising for u in the local limit is of Neumann type. Such idea is indeed proved rigorously performing the local asymptotics on the variational formulation for the nonlocal problem (1.1)-(1.4). The advantage of working using a variational approach is that the boundary conditions are implicitly contained in the variational formulation itself, and they have not to be tracked explicitly performing a pointwise analysis on the boundary. We are now in a position to present the two main theorems that we prove in this paper.
The first main result is the well-posedness for the nonlocal system (1.1)-(1.4) with Neumann boundary conditions when ε > 0 is fixed. Here, the viscosity coefficient τ ε is assumed to be strictly positive, the convolution kernel is of the form (1.11), and the double-well potential may be singular. In particular, we include in our analysis all the typical examples of polynomial, logarithmic, and double-obstacle potentials: In view of this, the derivative of is interpreted as a subdifferential in the sense of convex analysis, and equation (1.2) becomes a differential inclusion. The proof of well-posedness is based on a suitable approximation of the problem, given by a Yosida-type regularization on the nonlinearity and an additional elliptic local regularization in the chemical potential. A novel abstract variational setting (V ε , L 2 ( ), V * ε ) is introduced and uniform estimates on the approximated solutions are obtained. Using the viscous contribution in the chemical potential, strong compactness in L 2 is recovered even with no H 1 -estimates on the solutions. Strong convergences are then proved and a passage to the limit provides solutions to the original nonlocal problem.
The second main result of this paper is the asymptotic analysis of the nonlocal system as ε 0. Here, we assume that the forcing terms (g ε ) ε converge to a certain source g, and that the viscosity coefficients satisfy Here, the coefficient τ is allowed to be nonnegative: when τ > 0 we obtain then nonlocal-to-local convergence of viscous Cahn-Hilliard equations, while if τ = 0 we obtain the local asymptotics of nonlocal viscous Cahn-Hilliard equations with vanishing viscosities. The proof is based on uniform estimates in ε on the nonlocal solutions. Here, the strong compactness in L 2 is obtained by proving an ad-hoc compactness inequality involving the family on functional spaces (V ε ) ε>0 . The identification of the local limit − u is obtained through the combination of monotone analysis techniques and Gamma-convergence results for the nonlocal energy functional (1.5).
We conclude by highlighting some possible applications of our results to phasefield modelling. The relevance of nonlocal-to-local convergence of Cahn-Hilliard equations with Neumann boundary conditions is significant: among many others, we can mention here possible connections with optimal control of tumor growth models. In the recent years, phase-field models have been widely used in tumor growth dynamics, both in the local case (see [25,[31][32][33][34][35] and the references therein) and in the nonlocal case (see [26] and [48,54] for nonlocal Cahn-Hilliard equations with reaction terms). One of the main advantages of the nonlocal setting is that regularity results on the solutions are usually easier to obtain, not needing to rely on elliptic regularity properties. As a consequence, the availability of rigorous nonlocal-to-local convergence results would give the opportunity to approximate solutions to local phase-field systems with the solutions to the corresponding nonlocal ones, which are indeed simpler to handle on the mathematical side. For example, refined regularity on the solutions are fundamental when dealing with optimal control problems, in order to write first-order conditions for optimality. Hence, possible outcomes of nonlocal-to-local asymptotics concern refined analysis of optimal control of phase-field systems, in terms of passing to the (local) limit within first-order conditions for optimality for the nonlocal system.
The paper is structured in the following way: in Sect. 2 we state the assumptions, and we introduce the abstract variational settings. Section 3 is devoted to present the two main results. Section 4 contains the proof of well-posedness of the nonlocal system (1.1)-(1.4), while Sect. 5 focuses on the proof of nonlocal-to-local asymptotics.

Assumptions
Throughout the paper, is a smooth bounded domain in R d , with d = 2, 3, and T > 0 is a fixed final time. We will use the notation Q t := (0, t) × for every t ∈ (0, T ], and set Q := Q T , and := (0, T ) × ∂ . Moreover, (ρ ε ) ε>0 is a family of mollifiers with the following properties (see [51,52]): We define the family of convolution kernels as Throughout the paper, γ : R → 2 R is a maximal monotone graph with 0 ∈ γ (0) and : R → R is C -Lipschitz-continuous with (0) = 0. It follows in particular that there exists a proper, convex, lower semicontinuous functionγ : R → [0, +∞] withγ (0) = 0 and ∂γ = γ in the sense of convex analysis. Similarly, we set (s) := s 0 (r ) dr for every s ∈ R. With these notations, the double-well potential entering the system is represented by the sumγ +ˆ .

Variational Setting and Preliminaries
We introduce the functional spaces endowed with their natural norms, and we identify H with its dual space in the usual way, so that where all the inclusions are continuous, dense, and compact. The Laplace operator with homogeneous Neumann conditions will be intended both as a bounded linear operator and as unbounded linear operator on H with domain W . For every ϕ ∈ V * , we use the notation ϕ := 1 | | ϕ, 1 V for the mean value on . As a direct consequence of the Poincaré-Wirtinger inequality it holds that is a linear isomorphism. We will denote its inverse by For every ε > 0, we set Denoting by a ε : we also define Let us note that with such definitions the symmetry of K ε yields We point out that E ε : V ε → [0, +∞) is convex and B ε : H → H is a linear unbounded operator with domain W ε . Additionally, we define the maps and the bilinear forms We collect some properties in the next lemma.
(4) The following inclusions are continuous and dense: Moreover, such extension coincides with the linear operator A ε : V ε → V * ε associated to the bilinear form a ε , defined as Proof.
Step 1: properties (1)- (2). It is clear that · V ε and · W ε are norms on V ε and W ε , respectively. Let now (y n ) n be a Cauchy sequence in V ε : then in particular it is a Cauchy sequence in H , so there exists y ∈ H such that y n → y in H . By lower semicontinuity it follows that y ∈ V ε as well, and that y n → y in V ε . A similar argument shows that W ε is complete as well. A direct computation shows that (·, ·) V ε and (·, ·) W ε are scalar products inducing the norms above.
Step 2: property (3). For every ϕ ∈ C 0,σ ( ), we have The first term on the right-hand side is finite since 2 − σ < d, while the second term can be written as by the assumptions on (ρ ε ) ε . The thesis follows by the arbitrariness of x ∈ .
Step 3: property (4). First of all the fact that the inclusion V ε → H is continuous is trivial by the definition of · V ε . Second, for ϕ ∈ W ε , a direct computation shows that The monotonicity of B ε is a direct consequence of its definition. We proceed by showing that it is maximal monotone. Let ϕ ∈ H . For every λ, δ > 0 the elliptic problem admits a unique weak solution ϕ δλ ∈ W → C 0,1/4 ( ) → W ε , in the sense that Taking arbitrary ψ ∈ C ∞ c ( ) in this variational formulation, we infer that 2 ϕ δλ ∈ H in the sense of distributions. Consequently, the classical elliptic regularity theory yields also that ϕ δλ ∈ H 4 ( ), with ∂ n ϕ δλ = 0 almost everywhere on ∂ . Fix now δ > 0. Testing (2.2) by ϕ δλ and using the monotonicity of B ε and the Young inequality, it follows that Thus, there exists a positive constant M such that we have Noting that for all ζ ∈ V ε , by the symmetry of the kernel K ε and the Hölder inequality we have the estimates just proved ensure also that We infer that there exist ϕ δ ∈ V ε and η δ ∈ V * ε such that, as λ 0, λϕ δλ → 0 in W , ϕ δλ ϕ δ in V ε , and B ε (ϕ δλ ) η δ in V * ε . It follows that Now, for all ζ ∈ V ε , by the symmetry of B ε and the bilinearity of a ε if holds that This shows that so we conclude that ϕ δ ∈ W ε and η δ = B ε (ϕ δ ). Hence, Let us note that since we have just proved that ϕ δ ∈ W ε , in particular we have that B ε (ϕ δ ) ∈ H . Hence, if additionally ϕ ∈ V ε , testing (2.3) by B ε (ϕ δ ) and using Hölder and Young inequalities yields We deduce that, as δ 0, δ B ε (ϕ δ ) → 0 in H . Hence, by (2.3), ϕ δ → ϕ in H . By combining (2.4) and (2.5), we obtain that ϕ δ V ε ϕ V ε for every δ > 0. As V ε is uniformly convex, this implies that ϕ δ → ϕ in V ε , so that W ε → V ε densely.
Step 4: property (5). For every ϕ ∈ W ε and ζ ∈ V ε , by the Hölder inequality we have This implies that for every ϕ ∈ W ε , the operator is linear and continuous on V ε , and such that Since W ε → V ε is dense, we deduce that B ε extends to a bounded linear operator from V ε to V * ε , and the thesis follows.
Step 5: property (6). We observe that E ε : V ε → [0, +∞) is convex and lower semicontinuous. A direct computation also shows that DE ε = B ε in the sense of Gâteaux: since B ε : V ε → V * ε is linear and continuous, the thesis follows.
The next lemma shows some boundedness properties of the family (B ε ) ε , uniformly in ε.

Lemma 2. The following inclusion is continuous:
and there exists a constant C, independent of ε, such that For every ϕ, ζ ∈ V , if holds that (2.6) Finally, for every ϕ ∈ H and for every sequence (ϕ ε ) ε>0 ⊂ H with ϕ ε → ϕ in H , we have In other words, (E ε ) ε>0 -converges to E with respect to the norm-topology of H .
Proof. By [52], there is a constant C > 0 independent of ε such that from which the first part of the thesis follows directly. The first limit in (2.6) is also a direct consequence of [52], the second limit in (2.6) can be proved by choosing ϕ = ϕ 1 ± ϕ 2 in the first limit. Finally, by the -convergence result in [51, Thm. 8], we know that where sc −Ẽ is the lower semicontinuous envelope of It is a standard matter to check that sc −Ẽ = E, so that the thesis follows.
The last result of this section is a compactness criterion involving the family of operators (E ε ) ε . The following lemma is fundamental as we do not have any compactness properties for the inclusions of the spaces V ε and W ε . For the proof we refer to [24,Lemma. 4].

Main Results
Before stating our main results, we recall that the local Cahn-Hilliard equation is well-posed in the following sense: Then, there exists a triple (u, μ, ξ) such that

Moreover, the solution component u is unique, and the solution components μ and ξ are unique if γ is single-valued.
Proof. We refer to [18] for a proof in a more general setting.
The first result of this paper is the well-posedness of the nonlocal viscous Cahn-Hilliard equation complemented by Neumann boundary conditions for the chemical potential.
Our second contribution concerns the nonlocal-to-local convergence. In particular, we show that, under suitable assumptions on the initial data (u 0,ε ) ε and on the forcing terms (g ε ) ε , if the viscosities (τ ε ) ε converge to a coefficient τ 0, then the solutions to the respective viscous nonlocal Cahn-Hilliard equations converge, in suitable topologies, to the solutions to the limiting local Cahn-Hilliard equation with viscosity parameter τ 0. Note that the viscosities (τ ε ) ε are required to be strictly positive for all ε > 0, whereas the limiting viscosity parameter τ is also allowed to vanish. Hence, such result has a duplex formulation. Indeed, if τ > 0 this shows the asymptotic convergence of the nonlocal viscous equation to the corresponding local viscous equation, while if τ = 0 this proves the approximability of solutions to the local pure equation by solutions to nonlocal equations with vanishing viscosities.

Proof of Theorem 3.2
This section is devoted to the proof of well-posedness of the nonlocal viscous Cahn-Hilliard equation. Throughout the section, ε > 0 and τ ε > 0 are fixed.

Approximation
For every λ > 0, let γ λ : R → R be the Yosida approximation of γ , having Lipschitz constant 1/λ, and setγ λ (s) := s 0 γ λ (r ) dr for every s ∈ R. We consider the approximated problem where the initial datum u λ 0,ε satisfies for a certain λ 0 > 0 (possibly depending on ε). The existence of an approximating sequence (u λ 0,ε ) λ satisfying (4.5)-(4.6) is guaranteed by assumption (3.1): for example, one can check that the classical elliptic regularization given by the unique solution to the problem is a possible choice. The existence of a unique approximated solution (u λ ε , μ λ ε ) for every λ > 0 relies on a fixed-point argument, as in [24,Section 3.1]. For every v ∈ L 2 (0, T ; W ), since W → C 0, 1 4 ( ) by the Sobolev embeddings, thanks to the properties of B ε proved in Lemma 1 we have that B ε (v) ∈ L 2 (0, T ; H ). Hence, by the classical literature on the local viscous Cahn-Hilliard equation (see again [18]), the map is the unique solution to the local viscous Cahn-Hilliard equation Now, arguing as in [24, Section 3.1], exploiting the Lipschitz-continuity of γ λ , the Sobolev embeddings, and the properties of B ε contained in Lemma 1, we deduce that there exist constants L λ ε > 0 and σ > 0 such that, for every It follows that one can choose T 0 ∈ (0, T ] sufficiently small so that λ ε is a contraction on the respective functional spaces defined in (0, T 0 ). Performing then a classical patching argument (we refer again to [24, Section 3.1] for details), we infer that λ ε has a unique fixed point on the whole interval [0, T ]. This proves that the approximated system (4.1)-(4.4) has a unique solution

Uniform Estimates
We prove here some uniform estimates independently of λ and ε. In what follows we will always assume that λ ∈ [0, 1]. Moreover, ε > 0 and τ ε > 0 are still fixed.

Passage to the Limit as λ 0
In this section we analyze the passage to the limit as λ 0, with ε > 0 and τ ε > 0 still fixed. In view of the uniform bounds (4.8)-(4.12) and the Aubin-Lions lemma, up to the extraction of (not relabeled) subsequences we have the following convergences: for some From (4.14) and the fact that B ε ∈ L (V ε , V * ε ), it is readily seen that Moreover, from (4.15) and (4.19), it follows by comparison that We deduce in particular that B ε (u ε ) = η ε ∈ L 2 (0, T ; H ), so that also u ε ∈ L 2 (0, T ; W ε ). The strong convergence (4.13) implies also that u ε (0) = u 0,ε . Passing to the limit in (4.1)-(4.4) in the weak topology of L 2 (0, T ; H ), we obtain We proceed now providing an identification of the nonlinear terms ξ ε and ε : we adapt an argument performed in [22,Subsection 3.6]. To this end, since is Lipschitz-continuous, there exists α > 0 such that the operator γ + + ατ ε Id : R → 2 R is maximal monotone. For example, one can choose α := 2 τ ε L ∞ (R) (recall that τ ε > 0 is fixed). Multiplying (4.2) by e −αt , we obtain Thus, testing the previous equation by e −αt u λ ε and integrating in time yields On the one hand, owing to (4.13) and (4.16), On the other hand, by the weak lower semicontinuity of the norms, the convergence (4.14), and the assumption (4.5), we have ε (s, ·)) ds .
Hence, we deduce that Testing (4.21) by e −2αt u ε and integrating in time, the right-hand side of (4.24) rewrites as Since the bilinear form is an equivalent scalar product on L 2 (Q), by the maximal monotonicity of γ + + ατ ε Id we conclude that ξ ε + ε + ατ ε u ε ∈ (γ + + ατ ε Id)(u ε ) a.e. in Q . (4.25) This allows us to show the further strong convergences Indeed, taking the difference between (4.2) and (4.21), multiplying again by e −αt , and testing by e −αt (u λ ε − u ε ), we get We use now the notation J γ λ := (Id + λγ ) −1 : R → R for the resolvent of γ . Summing and subtracting J γ λ (u λ ε ) in the last term on the left-hand side, rearranging the terms, and recalling that u λ ε − J γ λ (u λ ε ) = λγ λ (u λ ε ), we infer that, for every t ∈ [0, T ], Recalling that γ λ (r ) ∈ γ (J γ λ (r )) for every r ∈ R, by (4.25) and the monotonicity of the operator γ + + ατ ε Id, the third term on the left-hand side is nonnegative. Let us show that the right-hand side converges to 0, analyzing each term separately. The first two terms on the right-hand side converge to 0 thanks to (4.5), (4.13) and (4.16). Moreover, thanks to (4.14), (4.19), and the fact that u ε ∈ L 2 (0, T ; W ε ), we have Finally, since (γ λ (u λ ε )) λ is bounded in L 2 (0, T ; H ) by (4.12), using the Lipschitzcontinuity of Pi, the last two terms on the right-hand side can be handled by Since t ∈ [0, T ] is arbitrary, the strong convergences (4.26) follows. In particular, this readily implies that ε = (u ε ) and ξ ε ∈ γ (u ε ) almost everywhere in Q by the Lipschitz-continuity of and by the maximal monotonicity of γ , respectively.
It is then clear that (u ε , μ ε , ξ ε ) is a solution to the nonlocal viscous Cahn-Hilliard equation in the sense of (3.11)- (3.16). This completes the proof of the first assertion of Theorem 3.2.
We observe that their difference solves By the assumption on the initial data, we have that (u 1 ε − u 2 ε ) = 0. Therefore, we can test the first equation by N (u 1 ε − u 2 ε ), the second by u 1 ε − u 2 ε , and take the difference: by performing classical computations we get The last term on the left-hand side is nonnegative by the monotonicity of γ . Hence, the continuous-dependence property stated in Theorem 3.2 follows from the Lipschitz-continuity of and the Gronwall lemma.
We test (3.14) by μ ε , (3.15) by ∂ t u ε , take the difference, and integrate on Q t : recalling (3.18) and using the Young inequality, we deduce that Note that 1 τ ε 1 τ * by (5.1). Hence, rearranging the terms and using (3.17) we infer that there exists a constant C > 0, independent of ε, such that hence also, by comparison in (3.14), Now, we can proceed as in the previous Sect. 4.2. Since (u ε ) = (u 0,ε ) , we can test (3.14) by N (u ε − (u 0,ε ) ), (3.15) by u ε − (u 0,ε ) , and sum: we obtain, for almost every t ∈ (0, T ), Again, by the estimates already performed, all the terms are bounded in L 2 (0, T ) except Thanks to assumption (3.19), there are two constants c, c > 0, independent of ε, such that Hence, we deduce that ξ ε L 2 (0,T ;L 1 ( )) C , which implies, by comparison in (3.15), that We deduce that Thus, by comparison in (3.15) and by monotonicity of γ , we obtain that By the Aubin-Lions compactness theorem we infer that, up to the extraction of (not relabeled) subsequences, as ε 0, We proceed by showing the strong convergence To this end, we show that the sequence (u ε ) ε is Cauchy in C 0 ([0, T ]; H ). For any arbitrary σ > 0, we apply Lemma 3 with the choice δ := σ 4C , where C > 0 is the constant obtained in the estimates above. We deduce that there existsε =ε σ and C σ > 0 such that for every ε 1 , ε 2 ∈ (0,ε σ ), for every t ∈ [0, T ]. Thanks to (5.2), there exists ε σ ∈ (0,ε σ ) such that Hence, taking the supremum in time and using the estimates above we infer that for every ε 1 , ε 2 ∈ (0,ε σ ). Since σ > 0 is arbitrary, we obtain the strong convergence (5.7). Now, from (5.7) and the Lipschitz continuity of , it follows that while the strong-weak closure of γ readily ensures that ξ ε ∈ γ (u ε ) almost everywhere in Q.
To conclude the proof of the theorem, it remains to prove additional spatial regularity for u and to provide an identification of η. First of all, note that since (u ε ) ε is bounded in L ∞ (0, T ; V ε ), by the Ponce criterion [51, Theorem 1.2] we have that u ∈ L ∞ (0, T ; V ).
The Aubin-Lions theorems ensure then that, up to not relabeled subsequence, as ε 0, Arguing as in the previous Sect. 5.1 thanks to the Lemma 3, the convergence (5.10) and the boundedness of (E ε (u ε )) ε in L ∞ (0, T ) imply the strong convergence Hence, by the Lipschitz continuity of we have while the strong-weak closure of γ yields ξ ε ∈ γ (u ε ) almost everywhere in Q. Moreover, still arguing as in the previous section we obtain that u ∈ L ∞ (0, T ; V ), η = − u, and u ∈ L 2 (0, T ; W ) by elliptic regularity.
This concludes the proof of Theorem 3.3 also in the case τ = 0.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/ licenses/by/4.0/.
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.