Convergence of solutions of a rescaled evolution nonlocal cross-diffusion problem to its local diffusion counterpart

We prove that, under a suitable rescaling of the integrable kernel defining the nonlocal diffusion terms, the corresponding sequence of solutions of the Shigesada–Kawasaki–Teramoto nonlocal cross-diffusion problem converges to a solution of the usual problem with local diffusion. In particular, the result may be regarded as a new proof of existence of solutions for the local diffusion problem.

Terming the equation (1) as the evolution nonlocal p−Laplacian equation has its reasons.Firstly, because it arises as the gradient descent of the Euler-Lagrange equation of the energy functional E nl p (v) = Ω Ω J(x − y)|v(y)) − v(x)| p dxdy, in analogy to the usual evolution p−Laplacian equation, for which the energy is given by E p (v) = Ω |∇v| p .Secondly, and most importantly, because under the rescaling J n (x) = n N +p J(nx), the corresponding sequence of solutions, u n , of problem (1)-(2) converges to the usual weak solution of the local p−Laplacian problem with the same initial datum v(0, •) = u 0 and with homogeneous Neumann boundary conditions.
The main idea behind the choice of the rescaling is that, being n N J(nx) an approximation of the Dirac delta, the factor n p plays the role of the denominator of the continuous incremental spatial ratio v(t, x + 1 n z) − v(t, x) 1/n p , so that one expects that, for all w ∈ W 1,p (Ω), as n → ∞.
In fact, Andreu et al. [3,Theorem 6.12] show, among other properties, that the sequence {u n } of solutions of the nonlocal rescaled problems converge to v, the solution of the local problem, strongly in L ∞ (0, T ; L p (Ω).There are two main ingredients in their proof.The first is the precompactness result [3,Theorem 6.11] based on previous results by Bourgain et al. [4,Theorem 4], which shows that if w n ⇀ w weakly in L p (Q T ) and E nl p,n (w n ) is uniformly bounded then w n → w strongly in L p (Ω) and w ∈ W 1,p (Ω).The second ingredient is the monotonocity of |s| p−2 s, which plays an important role both in the theory of existence of solutions and in the identification of the limit of the solutions of the rescaled problems.
The objective of this article is to show that the same convergence property is true for a class of evolution nonlocal cross-diffusion systems.In [16] (see also [12] for related work), we introduced and proved the existence of solutions of the nonlocal version of the paradigmatic Shigesada-Kawasaki-Teramoto (SKT) population model [22].The nonlocal version of this model is the following: for i = 1, 2, find where u = (u 1 , u 2 ) and, for i, j = 1, 2, i = j, the diffusion and reaction functions are given by for some non-negative constant coefficients c i , a i , α i , β ij .The local diffusion problem, i.e., the original version of the SKT model, is to find The existence of solutions of evolution cross-diffusion problems like ( 5)- (7) has been addressed for a variety of problems [2,14,5,15,9,17], for which one can define an appropriate Lyapunov functional, also known as entropy functional, which in the case of the SKT model is given by Formally testing (5) with ln(v i ) yields the estimate which is the first step for a compactness argument for a sequence of solutions of regularized problems.The second step, with the horizon of applying Aubin-Lion's type compactness arguments, is to deduce suitable uniform estimates for the time derivatives ∂ t v i .For the local SKT problem, the fundamental tool to obtain these estimates is the Gagliardo-Niremberg inequality, although other approaches (for other cross-diffusion problems), like global regularity in Hölder spaces [2] or duality estimates [9], are also fruitful.
For the nonlocal SKT problem a nonlocal entropy estimate analogous to ( 9) is deduced by similar procedures than in the local case, obtaining an estimate which is independent of the kernel J, see (11).A fundamental difference between the local and the nonlocal models is that, in the latter, the L ∞ (Q T ) regularity of solutions is proven, although with a bound depending on J.As a consequence, the time derivative bound is simpler to achieve than in the local case, and moreover, other important properties, like the uniqueness of solutions, are deduced based on this regularity.
Turning to the rescaled problems, the entropy estimate is again deduced as in the proof of existence of solutions, since this estimate is independent of the kernel.However, the estimation of the time derivative is more problematic.In the proof of our main theorem, and in view of the lack of a suitable nonlocal Gagliardo-Nirenberg inequality, we have resorted to the use of duality estimates, see e.g.[20,9], which allow us to show an improved regularity of the sequence of rescaled solutions (with respect to that implied by the entropy estimate).
Finally, to close this introduction, let us notice that an interesting by-product of the main result of this article is to regard it as a new proof of the existence of solutions of the usual local diffusion SKT problem.Besides, since the solutions of the approximating problems do not need to have Sobolev regularity, our approach opens the possibility of approximating numerically the solutions of the local SKT problem by alternative methods [19,10,7,8,13].
We recall here the result on existence and uniqueness of solution of problem (3)-(4) proven in [16].Notice that Hypothesis (H) 2 , needed for the result on convergence of the rescaled problems, is more restrictive than the corresponding assumption in [16] used to prove the existence and uniqueness of solution.We merged both conditions for the sake of brevity.In particular, this assumption implies the existence of some r > 0 such that supp(J) = B r (0).(10) Regarding the restriction u 0i ∈ BV (Ω), unusual for the local diffusion problem, we must include it since it is needed for the compactness argument employed in the proof of the existence of solutions of the nonlocal diffusion model.
Theorem A (Existence and uniqueness of solution [16]).Assume (H) and Then, there exists a unique strong solution (u 1 , u 2 ) of problem ( 3)-( 4) with u i ≥ 0 a.e. in Q T and such that, for i = 1, 2 and t ∈ [0, T ], with E u (t) defined by (8), and for some constant C > 0 independent of J.
Although not mentioned in the above theorem, integrating (3) in (0, t)×Ω yields the uniform estimate where, here and in what follows, C denotes a generic positive constant independent of J.
For stating the main result contained in this article, we introduce the rescaled kernel The kernel J n satisfies the conditions of Theorem A and, therefore, there exists a unique solution, u n , of (3)-( 4), corresponding to J n .We also introduce here, for later reference, the approximation to the Dirac delta Theorem 1. Assume (H) and suppose that a i > 0 for i = 1, 2. Let q = 3 − δ, with δ > 0 a small number, and r = 2q/(q + 2).Let u n be the solution of problem ( 3)-( 4) corresponding to the kernel J n and consider the sequence {u n }.There exists a subsequence {u n j } such that u n j → u 5)- (7) in the following weak sense: In addition, u satisfies the entropy estimate (9).
Notice that the regularity implied by the definition of the exponents q and p is independent of the spatial dimension, N .Regarding Theorem 1 as a proof of existence of solutions of problem ( 1)-( 2), our result improves that of [5] for N ≥ 3.For N = 1, 2, the regularity proven in [14,5] may be also recovered in our solution by using the Sobolev imbbeding theorem (N = 1) or the Gagliardo-Niremberg inequality (N = 2).
Due to the lack of knowledge concerning the uniqueness of solution of the local diffusion problem (1)-( 2), the convergence of the full sequence of rescaled problems may not be ensured, but only that of a subsequence to some solution of ( 1)-( 2).

Proof of Theorem 1
We start recalling a fundamental result by Bourgain et al. [4,Theorem 4].We use the variant introduced by Andreu et al. [3,Theorem 6.11] and state a straightforward extension to deal with time-dependent functions.First, we introduce the notation for the extension by zero of a function ψ : Ω → R: Let ρ : R N → R be a non-negative, continuous, radial, radially non-increasing function with compact support, and set ρ n (x) = n N ρ(nx).Let {f n } and {g n } be sequences in L p (Ω) and L p (Q T ), respectively, such that , for δ > 0.Then, there exists a constant C ≡ C(p, Ω, ρ) independent of n, and a number n δ ∈ N such that, for n ≥ n δ , (a) In consequence, there exists a subsequence {f n k } and a function f ∈ W 1,p (Ω) (BV (Ω), if p = 1) such that In addition, if Proof.Estimates ( 17) and (18), as well as their consequence (19), are proven in [4,3].Their extension to the time-dependent functions of (b) is straightforward once we notice that the results of (a) are valid pointwise for a.e.t ∈ (0, T ).Finally, ( 21) is also a direct consequence of [3,Theorem 6.11 (1)].✷ Remark 1.The inequality ( 17) is a Poincaré's type inequality, since the constant C 0 may be replaced by the nonlocal energy appearing in (16), see the proof of [4,Theorem 4] for details.Thus, this provides an alternative (constructive) proof of the result stated in [3,Proposition 6.19].
Step 1. Uniform bound in L 3 (Q T ).The entropy inequality (11) and the L 1 (Q T ) estimate (12) together with (20) applied with p = 2, ρ n = Jn , and g n = u in imply that the sequences {u in }, for i = 1, 2, are uniformly bounded in L 2 (Q T ).However, this bound is not enough to define the weak limit in an appropriate reflexive space of test functions since the nonlinear parts of the limit diffusive term are expected to be of the form with the regularity ∇u i ∈ L 2 (Q T ).Thus, we start improving the uniform bounds of {u in } to the space L 3 (Q T ).This bound is obtained using an estimate of the dual problem corresponding to the nonlocal rescaled problem.The existence of solutions of this dual problem is ensured by the following lemma.
The proof of this and the other lemmas used for proving Theorem 1 are given in Section 4. In the following, we shall use the notation , where λ > 0 is a constant and φ is the non-negative solution of ( 22)-( 23 with the same regularity than φ, see (24).
Proof of Corollary 1.By Theorem A, we have that pi (u n ) ∈ L ∞ (Q T ) is non-negative, so that ρ, a, b, c ∈ L ∞ are non-negative.By Lemma 1, there exists a unique non-negative solution φ of problem ( 22)-( 23) corresponding to this data.A simple calculation shows that ϕ in is then the non-negative solution of (25)-(26) inheriting the same regularity than φ.✷ Now we proceed to obtain the L 3 (Q T ) uniform bound of the sequences {u in }, for i = 1, 2. We multiply the equation (3) of u in by the solution ϕ in of problem (25)-( 26) and integrate to get, for i = 1, 2, Using the equation (25), the explicit expression of f i and the non-negativity of u in and ϕ in , we obtain Noticing that ψ ≤ 0 and using Hölder's inequality and the uniform estimate of u in in L 2 (Q T ), see Step 1, we obtain Our objective is to estimate ϕ in (0, Multiplying the equation (25) of ϕ in by ∆ n ϕ in and using the nonlocal integration by parts formula and Young's inequality yields We have so that from (28), we obtain Therefore, using Theorem B (a) with f n = ϕ in (0, •) and (b) with g n = ϕ in we deduce from (29) Integrating the equation (25

Convergence in rescaled nonlocal cross-diffusion problems
Observing that ϕ in and λ are non-negative and using Hölder's inequality, we get and, therefore, .
Returning to (30)-( 31) and taking into account the uniform L 1 (Q T ) estimate ( 12), we deduce Finally, (27) yields, by duality, an uniform estimate for u in pi (u n ) L 2 (Q T ) , or, in other words, the estimate In particular, if a i > 0, we obtain uniform estimates of u in in L 3 (Q T ).
Step 2. Strong convergence.The following lemma is a consequence of two results concerning compactness: the precompactness result of Bourgain et al. [4,Theorem 4] for sequences defined on the spatial domain Ω, and the compensated compactness result of P.L. Lions [18,Lemma 5.1] for the product of sequences defined in Q T .
Lemma 2. Let ρ : R N → R be like in Theorem B, and let {f n } be a sequence in L 3 (Q T ) such that (ii) for some m ≥ 0 independent of n.Then there exists a subsequence {f n k } and a function Taking ρ = J and f n = u in for i = 1, 2, the uniform estimate of u in in L 3 (Q T ) obtained in Step 1 and the entropy inequality satisfied by these functions, see (11), imply (i) and (ii) of Lemma 2. We now check that the uniform time estimate (iii) is also satisfied.For any smooth function ξ : Q T → R, we have Then there exist a constant C independent of n and a constant n J ∈ N such that if n > n J then . Using this lemma with p = 3 and noting that u in is uniformly bounded in L 3 (Q T ) we obtain from (36), by duality, Therefore, (iii) of Lemma 2 is satisfied and (35) follows, this is, there exist subsequences (not relabeled) such that u in → u i strongly in L q (Q T ), for i = 1, 2 and for any 1 ≤ q < 3.
Step 3. Time derivative estimate.Once we have proven the strong convergence of {u in } in L q (Q T ), we may improve the uniform time estimate obtained in (37).For ξ smooth, we have Clearly, For the linear term, we have, for i = 1, 2, where we used the entropy estimate (11) and an straightforward modification of [4, Theorem 1] for including the time variable.For the quadratic terms, we have, for i, j = 1, 2, where, for k, ℓ = 1, 2, Jn (z), being Jn an approximation of the Dirac delta, we have .
The first factor is bounded due to the uniform L 3 (Q T ) estimate found in Step 1.The second factor is bounded due to the entropy estimate (11).Finally, a new use of [4, Theorem 1] yields that the third factor is bounded by ξ L 6 (0,T ;W 1,6 (Ω)) .This is, we obtain Thus, by duality, we deduce from (38) Step 4. Identification of the limit.Since u in is a strong solution of ( 3)-( 4), we have, for ξ smooth, Using the strong convergence (35), we easily justify the passing to the limit n → ∞ for the reaction terms if ξ ∈ L q/(q−2) (Q T ).In view of (40), the terms involving the time derivative are well defined and some subsequences (not relabeled) pass to their corresponding limits (weakly) if ξ ∈ L 6 (0, T ; W 1,6 (Ω)).Regarding the diffusion term, we rewrite it as where ε = 1/n.We, again, examine the terms of p i (u) = c i u i + a i u 2 i + u i u j separately.For the linear term, we have, using (21), where we defined, for v ∈ R N , according to [3,Lemma 6.16].Observe that the convergence (41) may be extended, by density, to ξ ∈ L 2 (0, T ; H 1 (Ω)).Similarly, for u i u j we use, in addition to the weak convergence (21), the strong convergence deduced in Step 2. Splitting this term as in (39), we only have to examine, by symmetry, the following expression Since r ′ = 2q/(q − 2), we have where we used (42) and (21).Observe that, in this case, we may extend the functional space of test functions, by density, to ξ ∈ L r ′ (0, T ; W 1,r ′ (Ω)), and that, in relation to the weak convergence of the time derivatives, see (40), we have r ′ > 6.
Finally, a new and straightforward duality calculation shows that the initial data may be interpreted in the sense of (15).This finishes the proof of Theorem 1.
Regarding the regularity of the solution, since φ , and then, from the equation ( 22), we deduce that Finally, assume that a, b, c are non-negative and suppose that φ is negative somewhere.Let ξ(t, x) = φ(t, x) + ǫt, with ǫ > 0 small enough so that ξ is negative somewhere.Let (t 0 , x 0 ) be a point where ξ attains its negative minimum.Then, t 0 > 0 and Estimate (32) implies the existence of a subsequence {f n k } ⊂ L 3 (Q T ) and a function f ∈ L 3 (Q T ) such that f n k ⇀ f weakly in L 3 (Q T ).Moreover, estimates (34) and (43) ensure, see [18, Lemma 5.1], the convergence f 2 n k → f 2 in the sense of distributions in Q T (we take g n = h n = f n in [18, Lemma 5.1]).This is, Since f n k ⇀ f weakly in L 3 (Q T ), we have so that, passing to a new subsequence if required, we have f 2 n k ⇀ f 2 weakly in L 3/2 (Q T ), by density.As a consequence, we deduce that since 1 ∈ L 3 (Q T ) and f ∈ L 3/2 (Q T ).Then, the uniform bound (32) allows to obtain (35).✷ Proof of Lemma 3. We prove the result for a general power p, with 1 < p < ∞.By density, it is enough to show that ∆ n ψ L p (Ω) ≤ C ψ W 2,p 0 (Ω) for ψ ∈ C ∞ c (Ω). .