pith. sign in

arxiv: 2606.04049 · v1 · pith:J7B4ZBXXnew · submitted 2026-06-02 · 🧮 math.NA · cs.NA

Boundary-compatible interacting approximations of quasilinear PDEs on bounded domains

Thierry Paul (LYSM) , Emmanuel Tr\'elat (LJLL (UMR\_7598) , CaGE) This is my paper

Pith reviewed 2026-06-28 09:19 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords quasilinear PDEsregularizationdiscretizationerror estimatesboundary conditionsmollifiersinteracting approximationsBanach scales
0
0 comments X

The pith

Regularization and discretization yield quantitative error bounds separating scales for quasilinear PDE approximations on bounded domains.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops an operator-theoretic construction that converts Kato-type quasilinear evolution systems on a Banach scale into finite-dimensional interacting ODEs. It proceeds by first regularizing the drift via a family indexed by ε so that it lands in a space suitable for discretization, then applying a sampling-reconstruction pair to produce the discrete system. The central theorem bounds the discrepancy between the lifted discrete solution and the exact solution by isolating a regularization error term from a discretization error controlled by the regularized drift size times N to a negative power. In the concrete setting of quasilinear PDEs on bounded Lipschitz domains, boundary-compatible mollifiers keep the drift size uniformly bounded and deliver algebraic convergence rates for quasi-uniform grids.

Core claim

The main abstract theorem provides a quantitative estimate of the discrepancy y_ε^N - y between the lifted discrete solution and the exact one, separating the regularization error χ(ε) from the discretization error (1 + L_ε) N^{-γ}, where L_ε measures the size of the regularized drift in the output norm. This makes explicit the trade-off between the regularization scale ε, the discretization scale N, and the possible deterioration of L_ε as ε → 0. For the running example of quasilinear PDEs on bounded Lipschitz domains, Burenkov's variable-step mollifiers provide a boundary-compatible kernelization that regularizes differential operators into explicit integral-interaction operators supported

What carries the argument

The regularized family (A_ε, f_ε) indexed by ε > 0 that makes the drift A_ε[t,z]z + f_ε[t,z] take values in an output space Y suitable for discretization, together with the sampling-reconstruction pair (P_N, R_N) that produces an interacting ODE on the finite-dimensional space V_N.

Load-bearing premise

A regularized family (A_ε, f_ε) exists such that the drift lands in a discretization-friendly output space while the size of the regularized drift stays controlled as the regularization parameter approaches zero.

What would settle it

A numerical test in which the total observed error y_ε^N - y does not approach zero for sequences where χ(ε) is driven to zero but L_ε grows, or where the mollifiers fail to preserve boundary traces of the fields.

Figures

Figures reproduced from arXiv: 2606.04049 by CaGE), Emmanuel Tr\'elat (LJLL (UMR\_7598), Thierry Paul (LYSM).

Figure 1
Figure 1. Figure 1: Discretization diagram The discrete system on VN is u˙ N ε (t) = A N ε [t, uN ε (t)]u N ε (t) + f N ε [t, uN ε (t)] (27) with initial condition u N ε (0) = PN y 0 . Its lift in Z is y N ε (t) = RN u N ε (t) ∈ Z (28) and satisfies y˙ N ε (t) = QN Aε[t, yN ε (t)]y N ε (t) + QN fε[t, yN ε (t)] (29) 17 [PITH_FULL_IMAGE:figures/full_fig_p017_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Periodic transport interacting system (translation-invariant counterpart of Example [PITH_FULL_IMAGE:figures/full_fig_p029_2.png] view at source ↗
read the original abstract

We develop a general operator-theoretic route that turns Kato-type quasilinear evolution systems on a Banach scale $(Z,X)$ into finite-dimensional interacting approximations. The construction proceeds in two steps. First, one introduces a regularized family $(A_\varepsilon,f_\varepsilon)$ indexed by a scale parameter $\varepsilon>0$, for which the drift $A_\varepsilon[t,z]z+f_\varepsilon[t,z]$ takes values in an output space $Y$ suitable for discretization. Second, one discretizes this regularized dynamics by a sampling-reconstruction pair $(P_N,R_N)$ and obtains an interacting ODE on a finite-dimensional state space $V_N\simeq\R^{dN}$. Our main abstract theorem provides a quantitative estimate of the discrepancy $y_\varepsilon^N-y$ between the lifted discrete solution and the exact one, separating the regularization error $\chi(\varepsilon)$ from the discretization error $(1+L_\varepsilon)N^{-\gamma}$, where $L_\varepsilon$ measures the size of the regularized drift in the output norm. This makes explicit the trade-off between the regularization scale $\varepsilon$, the discretization scale $N$, and the possible deterioration of $L_\varepsilon$ as $\varepsilon\to 0$. As a running example, we focus on quasilinear PDEs on bounded Lipschitz domains with boundary conditions. We show that Burenkov's variable-step mollifiers provide a boundary-compatible kernelization: they regularize differential operators into explicit integral-interaction operators supported inside the domain and preserve boundary traces of sufficiently regular fields. In this setting one can choose an output space $Y$ for which $L_\varepsilon$ remains uniformly bounded, leading to algebraic convergence rates in $N$ for quasi-uniform discretizations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops an operator-theoretic framework for turning Kato-type quasilinear evolution systems on a Banach scale (Z,X) into finite-dimensional interacting approximations. It first regularizes the drift via a family (A_ε, f_ε) so that A_ε[t,z]z + f_ε[t,z] lands in a discretizable output space Y, then applies a sampling-reconstruction pair (P_N, R_N) to obtain an ODE on V_N ≃ ℝ^{dN}. The central abstract theorem bounds the lifted discrepancy y_ε^N - y by separating the regularization error χ(ε) from the discretization error (1 + L_ε)N^{-γ}, where L_ε controls the size of the regularized drift in Y. As a running example on bounded Lipschitz domains, Burenkov variable-step mollifiers are used to produce boundary-compatible integral operators while preserving traces, with the claim that a suitable Y exists making L_ε uniformly bounded and yielding algebraic rates in N.

Significance. If the uniform boundedness of L_ε is rigorously established, the result supplies an explicit, quantitative trade-off between regularization scale ε and discretization scale N that is independent of ε in the leading term. This is potentially useful for the numerical analysis of quasilinear PDEs on domains with boundary conditions, as it converts abstract evolution equations into concrete interacting particle or finite-element-type systems with controllable error.

major comments (2)
  1. [Main abstract theorem] Main abstract theorem: the claimed separation χ(ε) + (1 + L_ε)N^{-γ} is load-bearing only if L_ε remains bounded independently of ε. The manuscript asserts that Burenkov mollifiers achieve this on Lipschitz domains for a suitable Y, but the verification that the variable-step kernel support and trace-preservation properties produce ||A_ε[t,z]z + f_ε[t,z]||_Y ≤ L_ε with no ε-dependent blow-up (or only controllable factors) is not supplied with explicit estimates; without it the algebraic rate independent of ε does not follow.
  2. [Running example] Running example (quasilinear PDEs on bounded domains): the choice of output space Y and the assertion that mollification keeps the regularized drift in Y with uniform norm bound must be tied directly to the specific quasilinear coefficients and the Lipschitz character of the domain. If the mollification estimates only yield L_ε ≲ ε^{-α} for α > 0, the discretization term forces N to grow with ε, undermining the stated convergence claim for quasi-uniform discretizations.
minor comments (2)
  1. The notation for the Banach scale (Z,X), the output space Y, and the precise meaning of the lifted discrete solution y_ε^N should be introduced with a short diagram or table of spaces and operators before the abstract theorem.
  2. Clarify whether the sampling-reconstruction pair (P_N, R_N) is required to commute with the boundary-trace operator or only to be stable in Y; this affects how the boundary compatibility of the mollifiers is used.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on the manuscript. The points raised correctly identify that the main result's utility hinges on uniform control of L_ε, and we address each concern below with plans to strengthen the presentation.

read point-by-point responses

  1. Referee: [Main abstract theorem] Main abstract theorem: the claimed separation χ(ε) + (1 + L_ε)N^{-γ} is load-bearing only if L_ε remains bounded independently of ε. The manuscript asserts that Burenkov mollifiers achieve this on Lipschitz domains for a suitable Y, but the verification that the variable-step kernel support and trace-preservation properties produce ||A_ε[t,z]z + f_ε[t,z]||_Y ≤ L_ε with no ε-dependent blow-up (or only controllable factors) is not supplied with explicit estimates; without it the algebraic rate independent of ε does not follow.

    Authors: We agree that the separation in the main theorem yields the stated algebraic rates independent of ε only when L_ε is bounded uniformly. The manuscript constructs Y precisely so that the boundary-compatible properties of the Burenkov kernels (support strictly inside the domain and trace preservation) map the regularized drift into Y with controlled norm. However, the referee is right that the explicit estimates establishing ||A_ε[t,z]z + f_ε[t,z]||_Y ≤ L with L independent of ε are only outlined rather than derived in full detail with all constants. We will insert a new lemma in the running-example section that supplies these estimates, using the Lipschitz character of the domain, the normalization of the variable-step kernels, and the local Lipschitz assumption on the coefficients to obtain the uniform bound. revision: yes

  2. Referee: [Running example] Running example (quasilinear PDEs on bounded domains): the choice of output space Y and the assertion that mollification keeps the regularized drift in Y with uniform norm bound must be tied directly to the specific quasilinear coefficients and the Lipschitz character of the domain. If the mollification estimates only yield L_ε ≲ ε^{-α} for α > 0, the discretization term forces N to grow with ε, undermining the stated convergence claim for quasi-uniform discretizations.

    Authors: We concur that the uniform bound must be verified in terms of the given quasilinear coefficients and the Lipschitz geometry. The output space Y is chosen (as a suitable fractional Sobolev or weighted space adapted to the boundary) exactly to absorb the mollified terms without ε-growth for the assumed class of coefficients. The variable-step construction prevents the negative powers of ε that would otherwise appear near the boundary. We will add an explicit remark and a short calculation in the example section that traces the dependence on the Lipschitz constant and the coefficient bounds, confirming that L_ε remains O(1) and that the discretization rate therefore holds for quasi-uniform N independent of ε. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper presents a general abstract theorem that quantifies the discrepancy between lifted discrete solutions and exact solutions by separating regularization error χ(ε) from discretization error (1 + L_ε) N^{-γ}, under the assumption that a regularized family (A_ε, f_ε) exists with drift in a suitable output space Y. The running example asserts that Burenkov variable-step mollifiers achieve boundary-compatible regularization with L_ε uniformly bounded on Lipschitz domains, but this is derived from the properties of the mollifiers rather than by fitting parameters to data or reducing to self-citations. No steps reduce predictions to inputs by construction, invoke load-bearing self-citations, or smuggle ansatzes; the framework is self-contained against the stated assumptions and external mollifier estimates.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The abstract relies on the existence of suitable Banach scales, output spaces Y, and the boundary-trace preservation property of the chosen mollifiers; no numerical free parameters or new postulated entities are introduced.

axioms (2)
  • domain assumption Existence of a Banach scale (Z, X) and output space Y such that the regularized drift lands in Y
    Invoked to enable the discretization step after regularization
  • domain assumption Burenkov's variable-step mollifiers preserve boundary traces of sufficiently regular fields while keeping the regularized drift bounded in the output norm
    Central to obtaining uniform L_ε and algebraic rates on Lipschitz domains

pith-pipeline@v0.9.1-grok · 5858 in / 1525 out tokens · 21595 ms · 2026-06-28T09:19:53.122954+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

46 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    Babuˇ ska, J

    I. Babuˇ ska, J. M. Melenk,The partition of unity method, Internat. J. Numer. Methods Eng. 40 (1997), no. 4, 727–758

    1997

  2. [2]

    B´ enilan, M

    P. B´ enilan, M. Crandall, A. Pazy,Nonlinear evolution equations governed by accretive operators, book manuscript, Besan¸ con, 2001

    2001

  3. [3]

    S. C. Brenner, L. R. Scott, The Mathematical Theory of Finite Element Methods , 3rd ed., Texts in Applied Mathematics, vol. 15, Springer, New York, 2008

    2008

  4. [4]

    Brezis, Op´ erateurs maximaux monotones et semigroupes de contractions dans les espaces de Hilbert, North-Holland Math

    H. Brezis, Op´ erateurs maximaux monotones et semigroupes de contractions dans les espaces de Hilbert, North-Holland Math. Stud., No. 5, North-Holland Publishing Co., Amsterdam-London; American Elsevier Publishing Co., Inc., New York, 1973, vi+183 pp

    1973

  5. [5]

    V. I. Burenkov, Mollifying operators with variable step and their application to approximation by infinitely differentiable functions, Nonlinear analysis, function spaces and applications, Vol. 2 (Pisek, 1982), pp. 5–37. Teubner-Texte Math., 49, Leipzig, 1982

    1982

  6. [6]

    Burenkov, Sobolev spaces on domains , Teubner-Texte Math., 137, B

    V.I. Burenkov, Sobolev spaces on domains , Teubner-Texte Math., 137, B. G. Teubner Verlagsge- sellschaft mbH, Stuttgart, 1998. 312 pp

    1998

  7. [7]

    Carrillo, K

    J.A. Carrillo, K. Craig, F. Patacchini, A blob method for diffusion , Calc. Var. Partial Differential Equations 58 (2019), no. 2, Paper No. 53, 53 pp

    2019

  8. [8]

    Carrillo, A

    J.A. Carrillo, A. Esposito, J. S.-H. Wu, Nonlocal approximation of nonlinear diffusion equations, Calc. Var. Partial Differential Equations 63 (2024), no. 4, Paper No. 100, 44 pp

    2024

  9. [9]

    A. Chertock, A practical guide to deterministic particle methods, in Handbook of Numerical Methods for Hyperbolic Problems: Applied and Modern Issues , Handbook of Numerical Analysis, Vol. 18, Elsevier, 2017, pp. 177–202

    2017

  10. [10]

    Cottet and P

    G.-H. Cottet and P. D. Koumoutsakos, Vortex Methods: Theory and Practice , Cambridge University Press, 2000

    2000

  11. [11]

    Crandall, T.M

    M.G. Crandall, T.M. Liggett, Generation of semigroups of nonlinear transformations on general Banach spaces, Amer. J. Math. 93 (1971), 265–298

    1971

  12. [12]

    Degond, S

    P. Degond, S. Mas-Gallic, The weighted particle method for convection-diffusion equations. Part 1 and Part 2 , Math. Comp. 53 (1989), no. 188, 485–525. 46

    1989

  13. [13]

    Di Francesco, M.D

    M. Di Francesco, M.D. Rosini, Rigorous derivation of nonlinear scalar conservation laws from follow- the-leader type models via many particle limit , Arch. Ration. Mech. Anal. 217 (2015), no. 3, 831–871

    2015

  14. [14]

    Di Francesco, S

    M. Di Francesco, S. Fagioli, M.D. Rosini, Deterministic particle approximation of scalar conservation laws, Boll. Unione Mat. Ital. 10 (2017), no. 3, 487–501

    2017

  15. [15]

    Di Francesco, S

    M. Di Francesco, S. Fagioli, E. Radici, Deterministic particle approximation for nonlocal transport equations with nonlinear mobility , J. Differential Equations 266 (2019), no. 5, 2830–2868

    2019

  16. [16]

    D¨ orfler, H

    W. D¨ orfler, H. Gerner, R. Schnaubelt, Local well-posedness of a quasilinear wave equation , Appl. Anal. 95 (2016), no. 9, 2110–2123

    2016

  17. [17]

    D¨ orich, M

    B. D¨ orich, M. Hochbruck,Exponential integrators for quasilinear wave-type equations , SIAM J. Nu- mer. Anal. 60 (2022), no. 3, 1472–1493

    2022

  18. [18]

    D¨ orich,Strong norm error bounds for quasilinear wave equations under weak CFL-type conditions , Found

    B. D¨ orich,Strong norm error bounds for quasilinear wave equations under weak CFL-type conditions , Found. Comput. Math. 25 (2025), 303–350

    2025

  19. [19]

    Engel, R

    K.-J. Engel, R. Nagel, One-parameter semigroups for linear evolution equations , Graduate Texts Math. 194, Springer-Verlag, 2000

    2000

  20. [20]

    Eymard, T

    R. Eymard, T. Gallou¨ et, R. Herbin, Finite volume methods, in Handbook of Numerical Analysis , vol. VII, P. G. Ciarlet and J. L. Lions, eds., North-Holland, Amsterdam, 2000, pp. 713-1020

    2000

  21. [21]

    Godlewski, P.-A

    E. Godlewski, P.-A. Raviart, Numerical Approximation of Hyperbolic Systems of Conservation Laws , Applied Mathematical Sciences, vol. 118, Springer, New York, 1996

    1996

  22. [22]

    Hinterm¨ uller, K

    M. Hinterm¨ uller, K. Papafitsoros, C. N. Rautenberg,Variable step mollifiers and applications, Integral Equations Operator Theory 92 (2020), no. 6, Paper No. 53, 34 pp

    2020

  23. [23]

    Hochbruck, T

    M. Hochbruck, T. Paˇ zur,Error analysis of implicit Euler methods for quasilinear hyperbolic evolution equations, Numer. Math. 135 (2017), no. 2, 547–569

    2017

  24. [24]

    Hochbruck, T

    M. Hochbruck, T. Paˇ zur, R. Schnaubelt, Error analysis of implicit Runge-Kutta methods for quasi- linear hyperbolic evolution equations, Numer. Math. 138 (2018), no. 3, 557–579

    2018

  25. [25]

    R. W. Hockney and J. W. Eastwood, Computer Simulation Using Particles , Hilger, Bristol, 1988

    1988

  26. [26]

    Hughes, T

    T. Hughes, T. Kato, J.E. Marsden, Well-posed quasi-linear second-order hyperbolic systems with applications to nonlinear elastodynamics and general relativity, Arch. Rational Mech. Anal. 63 (1976), no. 3, 273–294

    1976

  27. [27]

    K. Ito, F. Kappel, The Trotter-Kato theorem and approximation of PDEs , Math. Comp. 67 (1998), no. 221, 21–44

    1998

  28. [28]

    K. Ito, F. Kappel, Evolution equations and approximations , Ser. Adv. Math. Appl. Sci., 61, World Scientific Publishing Co., Inc., River Edge, NJ, 2002, xiv+498 pp

    2002

  29. [29]

    Kato, Nonlinear semigroups and evolution equations , J

    T. Kato, Nonlinear semigroups and evolution equations , J. Math. Soc. Japan 19 (1967), 508–520

    1967

  30. [30]

    Kato, Quasi-linear equations of evolution, with applications to partial differential equations , Spec- tral theory and differential equations (Proc

    T. Kato, Quasi-linear equations of evolution, with applications to partial differential equations , Spec- tral theory and differential equations (Proc. Sympos., Dundee, 1974; dedicated to Konrad J¨ orgens), pp. 25–70, Lecture Notes in Math., Vol. 448, Springer-Verlag, Berlin-New York, 1975

    1974

  31. [31]

    T. Kato, Abstract evolution equations, linear and quasilinear, revisited, Functional analysis and related topics, 1991 (Kyoto), 103–125, Lecture Notes in Math., 1540, Springer-Verlag, Berlin, 1993

    1991

  32. [32]

    Kov` acs, C

    B. Kov` acs, C. Lubich,Stability and convergence of time discretizations of quasi-linear evolution equa- tions of Kato type , Numer. Math. 138 (2018), no. 2, 365–388

    2018

  33. [33]

    R. J. LeVeque, Finite Volume Methods for Hyperbolic Problems , Cambridge Texts in Applied Math- ematics, Cambridge University Press, Cambridge, 2002

    2002

  34. [34]

    Lions, S

    P.-L. Lions, S. Mas-Gallic, Une m´ ethode particulaire d´ eterministe pour des ´ equations diffusives non lin´ eaires, C. R. Acad. Sci. Paris S´ er. I Math.332 (2001), no. 4, 369–376

    2001

  35. [35]

    Medvedev, The nonlinear heat equation on dense graphs and graph limits , SIAM J

    G.S. Medvedev, The nonlinear heat equation on dense graphs and graph limits , SIAM J. Math. Anal. 46 (2014), no. 4, 2743–2766. 47

    2014

  36. [36]

    J. J. Monaghan, Smoothed particle hydrodynamics, Annual Review of Astronomy and Astrophysics 30 (1992), 543–574

    1992

  37. [37]

    M¨ uller,Well-posedness for a general class of quasilinear evolution equations with applications to Maxwell’s equations, PhD thesis, Karlsruhe Institute of Technology, 2014

    D. M¨ uller,Well-posedness for a general class of quasilinear evolution equations with applications to Maxwell’s equations, PhD thesis, Karlsruhe Institute of Technology, 2014

    2014

  38. [38]

    T. Paul, E. Tr´ elat,Mean field, hydrodynamic and graph limits for deterministic interacting particle systems: a survey with quantitative estimates , Preprint arXiv:2209.08832v4

    work page internal anchor Pith review Pith/arXiv arXiv

  39. [39]

    Pazy, Semigroups of linear operators and applications to partial differential equations , Applied Mathematical Sciences, 44, Springer-Verlag, New York, 1983

    A. Pazy, Semigroups of linear operators and applications to partial differential equations , Applied Mathematical Sciences, 44, Springer-Verlag, New York, 1983. viii+279 pp

    1983

  40. [40]

    Russo, Deterministic diffusion of particles , Comm

    G. Russo, Deterministic diffusion of particles , Comm. Pure Appl. Math. 43 (1990), no. 6, 697–733

    1990

  41. [41]

    Sanekata, Abstract quasi-linear equations of evolution in nonreflexive Banach spaces , Hiroshima Math

    N. Sanekata, Abstract quasi-linear equations of evolution in nonreflexive Banach spaces , Hiroshima Math. J. 19 (1989), no. 1, 109–139

    1989

  42. [42]

    Schnaubelt, M

    R. Schnaubelt, M. Spitz, Local well-posedness of quasilinear Maxwell equations with conservative interface conditions, Nonlinearity 36 (2023), no. 12, 6712–6746

    2023

  43. [43]

    L. R. Scott, S. Zhang, Finite element interpolation of nonsmooth functions satisfying boundary con- ditions, Math. Comp. 54 (1990), no. 190, 483–493

    1990

  44. [44]

    Sogge, Fourier integrals in classical analysis , Second edition, Cambridge Tracts in Math., 210, Cambridge University Press, Cambridge, 2017

    C. Sogge, Fourier integrals in classical analysis , Second edition, Cambridge Tracts in Math., 210, Cambridge University Press, Cambridge, 2017. xiv+334 pp

    2017

  45. [45]

    Thom´ ee, Galerkin Finite Element Methods for Parabolic Problems , 2nd ed., Springer Series in Computational Mathematics, vol

    V. Thom´ ee, Galerkin Finite Element Methods for Parabolic Problems , 2nd ed., Springer Series in Computational Mathematics, vol. 25, Springer, Berlin, 2006

    2006

  46. [46]

    Wendland, Scattered Data Approximation, Cambridge University Press, 2005

    H. Wendland, Scattered Data Approximation, Cambridge University Press, 2005. 48

    2005