Nonlocal-to-Local Limiting Strategy

• The nonlocal-to-local limiting strategy is a framework where families of nonlocal models converge to a local PDE or operator as parameters like interaction kernels concentrate or horizons vanish.
• It employs structural identities, moment balance, and variational methods to ensure convergence while preserving key features such as diffusion, drift, and curvature.
• The approach finds applications in traffic models, viscoelasticity, minimal surfaces, and graph transport, with convergence verified via compactness and energy–dissipation estimates.

Nonlocal-to-local limiting strategy denotes a class of analytical and computational procedures in which a family of nonlocal models is shown to converge to a local partial differential equation, variational problem, or differential operator as the interaction kernel concentrates, the horizon vanishes, or a singular parameter approaches its local endpoint. In the cited literature, this passage is realized in several distinct regimes: convolution kernels are rescaled to a Dirac delta in nonlocal conservation laws and Cahn–Hilliard systems, the order parameter $a\uparrow 2$ in nonsymmetric jump processes, the fractional parameter $s\to 1^{-}$ in nonlocal minimal surfaces, the peridynamic horizon tends to $0$ in linearized viscoelasticity and elasticity coupling, and graph-localized transport weights converge to a tensor-mobility continuity equation in Euclidean space (Coclite et al., 2023, Davoli et al., 2019, Kassmann et al., 2022, Caffarelli et al., 2011, Friedrich et al., 2024, Esposito et al., 2023).

1. Representative limiting regimes

The underlying limit parameter varies by model class, but the objective is the same: identify a local object that retains the correct diffusion, drift, entropy, curvature, or mobility structure of the nonlocal system. In some settings the limit is a scalar conservation law, in others a Laplacian, a second-order divergence-form operator with drift, a local Kelvin–Voigt system, or a fourth-order aggregation–diffusion equation.

Nonlocal setting	Limiting parameter	Local limit
Nonlocal conservation laws with rescaled kernels	$\varepsilon\to 0^+$	$\partial_t \rho+\partial_x(V(\rho)\rho)=0$
Polygonally approximated nonlocal diffusion	$\delta\to 0$	$\mathcal L_0=\Delta$
Nonsymmetric nonlocal forms	$a\uparrow 2$	$\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$
Nonlocal minimal surfaces	$s\to 1^{-}$	classical perimeter / mean curvature
Nonlocal viscous Cahn–Hilliard	$s\to 1^{-}$0	local Cahn–Hilliard with $s\to 1^{-}$1
Peridynamic Kelvin–Voigt viscoelasticity	horizon $s\to 1^{-}$2	$s\to 1^{-}$3
Localizing infinite graphs	$s\to 1^{-}$4	$s\to 1^{-}$5

For one-dimensional traffic-type conservation laws, the local limit is typically produced by a kernel family $s\to 1^{-}$6 with $s\to 1^{-}$7; the limit equation is the scalar conservation law with flux $s\to 1^{-}$8 (Coclite et al., 2023). In nonlocal diffusion on $s\to 1^{-}$9, the normalization $0$0 identifies the Laplacian as the local limit and makes asymptotic compatibility a moment condition rather than a purely geometric one (Du et al., 2021). For nonsymmetric forms, the scaling $0$1 produces a second-order local energy, while the antisymmetric part survives as a first-order drift term (Kassmann et al., 2022).

A closely related pattern appears in short-range aggregation models and Gross–Pitaevskii equations. In the aggregation–diffusion setting, the scaling $0$2 and a Taylor expansion of the convolution produce a fourth-order thin-film/Cahn–Hilliard-type equation (Falcó et al., 13 May 2025). In the nonlocal Gross–Pitaevskii problem, a family $0$3 yields convergence of nonlocal dark solitons to the explicit local dark solitons of the one-dimensional Gross–Pitaevskii equation (Laire et al., 2024).

2. Structural devices that make the limit identifiable

A recurring feature of successful nonlocal-to-local arguments is the availability of a structural identity that converts the nonlocal quantity into an object with a local evolution equation. For the one-sided exponential kernel in traffic models, the “look-ahead” variable

$0$4

satisfies

$0$5

More importantly, $0$6 obeys a transport equation with a structured nonlocal source term. This closure property is the backbone of the Oleĭnik-type approach and is specific to the exponential kernel (Coclite et al., 2023). The same exponential structure underlies the lane-changing balance-law system, where

$0$7

allows the system to be rewritten in terms of $0$8, which is the quantity for which uniform BV estimates are proved (Chiarello et al., 2023).

In nonlocal diffusion and elasticity, the corresponding structural ingredient is moment balance. The local limit depends on exact second-moment normalization, and the wrong geometry can destroy the limit even when the approximation looks visually close. In the polygonal-neighborhood problem, uniformly bounded side number yields

$0$9

so the approximated model can converge to the wrong multiple of the Laplacian; convergence is recovered when the polygonal refinement improves with $\varepsilon\to 0^+$0, for example when $\varepsilon\to 0^+$1 in the quasi-uniform inscribed case (Du et al., 2021). In peridynamic Kelvin–Voigt viscoelasticity, radial symmetry, kernel normalization, and concentration near zero identify the limit tensors $\varepsilon\to 0^+$2 and $\varepsilon\to 0^+$3 and make the nonlocal strain quadratic forms converge to the local expressions in $\varepsilon\to 0^+$4 and $\varepsilon\to 0^+$5 (Friedrich et al., 2024).

Operator-theoretic and graph limits use the same principle in a different language. For nonsymmetric jump kernels, the coefficients

$\varepsilon\to 0^+$6

identify the diffusion matrix and the drift, respectively (Kassmann et al., 2022). In the graph interaction equation, the limit tensor

$\varepsilon\to 0^+$7

encodes the anisotropy of the localizing graph and is obtained from the localized edge weights by second-moment asymptotics (Esposito et al., 2023). These examples show that nonlocal-to-local limits are usually governed by a precise moment mechanism rather than by weak convergence alone.

3. Compactness, entropy, and admissibility

The central analytical obstacle is that weak convergence is usually insufficient to pass to the nonlinear flux or to identify the correct local solution. The conservation-law literature makes this point explicitly: total variation of $\varepsilon\to 0^+$8 itself may blow up, so the compact quantity is often the averaged field $\varepsilon\to 0^+$9, not the original state (Colombo et al., 2023).

One major route is an Oleĭnik-type one-sided estimate. For the exponential-kernel traffic model, the nonlocal term $\partial_t \rho+\partial_x(V(\rho)\rho)=0$0 satisfies

$\partial_t \rho+\partial_x(V(\rho)\rho)=0$1

and under stronger structural assumptions the quantity $\partial_t \rho+\partial_x(V(\rho)\rho)=0$2 satisfies a one-sided bound of the form

$\partial_t \rho+\partial_x(V(\rho)\rho)=0$3

These estimates control oscillations, generate BV regularization for transformed quantities, yield precompactness in $\partial_t \rho+\partial_x(V(\rho)\rho)=0$4, and exclude non-entropic shocks in the limit. A notable feature is that the convergence to the local entropy solution is proved under the sole assumption $\partial_t \rho+\partial_x(V(\rho)\rho)=0$5, without assuming $\partial_t \rho+\partial_x(V(\rho)\rho)=0$6 (Coclite et al., 2023).

A second route is compensated compactness. For the class

$\partial_t \rho+\partial_x(V(\rho)\rho)=0$7

strong $\partial_t \rho+\partial_x(V(\rho)\rho)=0$8-convergence is obtained without total variation bounds or Oleĭnik estimates by proving compactness of the entropy productions

$\partial_t \rho+\partial_x(V(\rho)\rho)=0$9

in $\delta\to 0$0 and then applying Tartar–Murat compensated compactness. The two settings treated are the piecewise constant kernel $\delta\to 0$1 with Greenshields velocity $\delta\to 0$2, and strictly monotone kernels with decreasing velocities, including the exponential kernel (Coclite et al., 19 Nov 2025). In the anisotropic traffic overview, a different compactness estimate is proved: if $\delta\to 0$3 is supported in $\delta\to 0$4, nondecreasing, and convex on $\delta\to 0$5, then

$\delta\to 0$6

for $\delta\to 0$7, which again shifts compactness from $\delta\to 0$8 to its nonlocal average (Colombo et al., 2023).

Variational limits can also encode admissibility and boundary information. In the nonlocal viscous Cahn–Hilliard equation, the nonlocal operator $\delta\to 0$9 converges to $\mathcal L_0=\Delta$0 in the variational sense, and the boundary condition $\mathcal L_0=\Delta$1 is not imposed at the nonlocal level but emerges from the $\mathcal L_0=\Delta$2-limit of the energies (Davoli et al., 2019). In the control problem for conservation laws, the improved nonlocal-to-local convergence theorem allows simultaneous convergence of the kernels and the initial data, which is the key step in the subsequent $\mathcal L_0=\Delta$3-convergence of the objective functionals (Friedrich et al., 1 Oct 2025).

4. Variational, evolutionary, and operator-theoretic formulations

A large part of the modern theory formulates the nonlocal-to-local passage at the level of energies, forms, semigroups, or gradient structures rather than only at the level of weak solutions. In nonsymmetric nonlocal forms, the convergence is proved in the sense of Mosco–Hino: for $\mathcal L_0=\Delta$4, the forms $\mathcal L_0=\Delta$5 converge to $\mathcal L_0=\Delta$6, which implies convergence of resolvents, semigroups, and dual semigroups. The limiting local form is

$\mathcal L_0=\Delta$7

making precise the statement that antisymmetry behaves like a lower-order drift (Kassmann et al., 2022).

A broader abstraction is provided by nonlocal $\mathcal L_0=\Delta$8-convergence on a closed exact complex $\mathcal L_0=\Delta$9. There the coefficients are bounded linear operators $a\uparrow 2$0, not necessarily multiplication operators, and convergence is defined through the solution operators of two dual elliptic problems. The theory proves uniqueness of the nonlocal $a\uparrow 2$1-limit, relative compactness in the topology $a\uparrow 2$2, a block-matrix characterization through $a\uparrow 2$3, $a\uparrow 2$4, $a\uparrow 2$5, and the Schur complement, and the coincidence of local and nonlocal $a\uparrow 2$6-convergence for multiplication operators on the standard $a\uparrow 2$7 complex (Waurick, 2018).

Gradient-flow formulations provide the evolutionary analogue of these static convergence principles. In linearized viscoelasticity, the quasistatic nonlocal Kelvin–Voigt evolution is written as

$a\uparrow 2$8

and studied through the exact energy–dissipation balance

$a\uparrow 2$9

Following an evolutionary $\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$0-convergence argument, the nonlocal solutions converge strongly in $\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$1 to the unique weak solution of the local Kelvin–Voigt problem (Friedrich et al., 2024). The graph interaction equation uses the same philosophy in a different geometry: graph dynamics are Finslerian gradient flows, the continuum limit is a weighted Wasserstein gradient flow with mobility tensor $\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$2, and the passage is formulated as EDP convergence of gradient structures (Esposito et al., 2023).

Optimization problems also fit naturally into this framework. In the basis-pursuit limit for optimal design, the antisymmetric nonlocal mixed-norm problem provides only a one-sided estimate toward the local measure-valued limit, whereas dropping antisymmetry yields genuine $\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$3-convergence to the relaxed local dual problem (Evgrafov et al., 2021). In the control of conservation laws, the objective functionals

$\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$4

$\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$5-converge to the corresponding local functional, and minimizers of the nonlocal problems converge, up to subsequences, to minimizers of the local problem (Friedrich et al., 1 Oct 2025).

5. Boundary conditions, coupling, and asymptotic compatibility

A major computational manifestation of the nonlocal-to-local strategy is the construction of boundary or interface conditions that remain correct as the horizon shrinks. One approach converts local surface data into nonlocal volume constraints. In the nonlocal diffusion problem based on the nonlocal gradient/divergence calculus, the procedure is: solve a local Poisson problem with the available surface data, define the nonlocal Neumann data by

$\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$6

and then solve the nonlocal problem with those generated volume constraints. The resulting nonlocal solution $\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$7 satisfies the quadratic estimates

$\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$8

so the vanishing-horizon limit recovers the local solution in both the energy seminorm and $\partial_i(a_{ij}\partial_j u)+b_i\partial_i u$9 (D'Elia et al., 2019).

A second approach is local-to-nonlocal coupling by optimization. In the peridynamic–FEM setting, a static Linear Peridynamic Solid model is coupled with classical Navier–Cauchy elasticity on overlapping subdomains by minimizing

$s\to 1^{-}$0

subject to the two state equations and the virtual control data on the interface. The method relies on the facts that $s\to 1^{-}$1 as $s\to 1^{-}$2 and that the LPS and NC equations are exactly equivalent for quadratic displacement fields (D'Elia et al., 2021).

The term asymptotically compatible is used explicitly in the local-to-nonlocal heat coupling literature. In the nonoverlapping Robin–Dirichlet coupling problem, a nonlocal Neumann-type boundary treatment is constructed by converting local flux into a collar-volume correction, then extended to a nonlocal Robin-type condition. The continuous formulation recovers the local Neumann problem with $s\to 1^{-}$3 accuracy in $s\to 1^{-}$4, while the fully discrete coupled scheme shows asymptotic convergence to the local limit with observed $s\to 1^{-}$5 rate when $s\to 1^{-}$6 is fixed (You et al., 2019). In one-dimensional diffusion, the quasinonlocal coupling method based on geometric reconstruction gives a variationally consistent mixed operator, passes the patch test, preserves flux balance and the maximum principle, and satisfies

$s\to 1^{-}$7

as the interaction horizon shrinks (Du et al., 2017).

These results also delimit what asymptotic compatibility does not permit. A polygonal approximation of a ball-shaped interaction neighborhood is not automatically safe in the vanishing-horizon regime: if the number of sides stays uniformly bounded, the operator can converge to the wrong local limit; geometric refinement must therefore accompany horizon refinement (Du et al., 2021).

6. Scope, limitations, and open directions

The nonlocal-to-local limiting strategy is used across conservation laws, diffusion, viscoelasticity, minimal-surface theory, optimal design, graph transport, aggregation–diffusion, and dispersive wave problems. In traffic applications, it has been used for single-lane conservation laws, nonlocal lane-changing balance laws, and optimal control with the initial datum as control (Coclite et al., 2023, Chiarello et al., 2023, Friedrich et al., 1 Oct 2025). In geometry, it provides the bridge from $s\to 1^{-}$8-minimal surfaces to classical minimal surfaces: as $s\to 1^{-}$9, $s\to 1^{-}$00 converges to the classical perimeter and the normalized nonlocal curvature converges to mean curvature, enabling an improvement-of-flatness theorem uniform in $s\to 1^{-}$01 near $s\to 1^{-}$02 (Caffarelli et al., 2011). In dispersive problems, it yields convergence of nonlocal gray and black solitons to the explicit local Gross–Pitaevskii dark solitons (Laire et al., 2024).

Several limitations are structural rather than technical. The local limit can be wrong if geometric fidelity is too coarse, as in polygonal interaction neighborhoods with bounded side number (Du et al., 2021). In basis-pursuit optimal design, antisymmetry is physically natural and mathematically useful for compactness, but it weakens the nonlocal-to-local passage to a one-sided estimate; genuine $s\to 1^{-}$03-convergence requires dropping the antisymmetry constraint (Evgrafov et al., 2021). In anisotropic conservation laws, the compactness mechanism depends critically on one-sided kernels and the sign structure of $s\to 1^{-}$04; symmetric or nonconvex kernels need different arguments (Colombo et al., 2023). These examples show that the local limit is not a universal consequence of “small horizon” alone.

Open problems remain model-dependent. In the local fourth-order aggregation–diffusion model obtained from short-range adhesion, the cited work emphasizes unresolved questions on uniqueness of weak solutions, existence of minimizers on the whole space, long-time asymptotics, and the extent of validity of the local approximation for real co-culture experiments (Falcó et al., 13 May 2025). A plausible implication is that future progress will continue to depend on choosing the transformed quantity, variational structure, or interface formulation that makes the small-horizon limit rigid enough to identify the correct local object. Across the cited literature, that choice is the defining content of the nonlocal-to-local limiting strategy.