Pith. sign in

REVIEW 3 major objections 5 minor 32 references

Numerical analysis of a semilinear fractional diffusion equation

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Finite element error bounds are proved for rough-data semilinear fractional diffusion.

desk verdict Worth refereeing, but the headline pointwise spatial estimate and the new Gronwall tool are not fully proven as written. read the letter →

arxiv 1909.00016 v1 pith:5TR3PKI7 submitted 2019-08-30 math.NA cs.NA

Binjie Li , Tao Wang , Xiaoping Xie This is my paper · ORCID
classification math.NAcs.NA MSC 65M6035B65
keywords semilinearfractionaldiffusionequationnonsmoothinitialdatafiniteelementmethodgradedtemporalgridGrönwallinequalitySobolevspaceserrorestimatesRiemann-Liouvillederivative
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper proves the first numerical error estimates for a semilinear fractional diffusion equation in which the initial data is only in $L^2(\Omega)$ rather than in a smoother Sobolev space. The central spatial result is a pointwise $L^2(\Omega)$ bound $\|(u-u_h)(t)\|_{L^2(\Omega)} \lesssim h^2(t^{-\alpha}+\ln(1/h))\|u_0\|_{L^2(\Omega)}$ for piecewise-linear finite elements, and the central temporal result is an $O(J^{-(1-\alpha)/2}\sqrt{\ln J})$ error bound for a discontinuous Galerkin time discretization on graded grids. Both estimates are new because earlier analyses required $u_0\in\dot H^2(\Omega)$. The paper attributes the gain to a new fractional Grönwall inequality and its discrete analogue, which convert fractional-derivative energy bounds into true error control. A sympathetic reader should care because these are the first rigorous convergence rates for a class of problems that are known to develop time singularities even from smooth data.

What carries the argument

The argument is carried by two interlocking devices. The first is a new fractional Grönwall inequality, Lemma 3.5: if $\|D^\beta_{0+} y\|^2_{L^2(0,t)} \le \epsilon + A\|y\|^2_{L^2(0,t)}$ for all $0<t<T$ with $0<\beta<1/2$, then $\|D^\beta_{0+}y\|_{L^2(0,t)} \le C_{\beta,A,T}\sqrt{\epsilon}$; Lemma 5.1 is its discrete, piecewise-constant analogue on graded grids. These inequalities let the authors turn the usual energy identity, which controls a fractional derivative of the error, into control of the error itself without assuming smoothness of $u_0$. The second device is the integral representation $Sg(t)=\int_0^t E(t-s)g(s)\,ds$ of the solution operator in terms of the kernel $E$ built from the resolvent of $-\Delta$, together with the spatial projection estimate Lemma 4.1, Eq. (48), $\|(S D^\alpha_{0+}v-S_h D^\alpha_{0+}P_h v)(t)\|_{L^2(\Omega)} \lesssim h^2 t^{-\alpha}\|v\|_{L^2(\Omega)}$, which transfers the nonsmooth linear estimate to the semidiscrete setting.

What would settle it

Evaluate the left side of Eq. (48) numerically for $v(x)=x^{-0.49}$ on $\Omega=(0,1)$ at a fixed $t>0$: if $\|(S D^\alpha_{0+}v - S_h D^\alpha_{0+}P_h v)(t)\|_{L^2(\Omega)}$ is not bounded by a constant times $h^2 t^{-\alpha}\|v\|_{L^2(\Omega)}$ on a sequence of uniformly refined meshes, the omitted lemma fails and Theorem 4.1 collapses.

Watch

Extended reading notes

Core claim

The paper's central claim is Theorem 4.1: for each $0<t\le T$, the semidiscrete piecewise-linear finite element solution $u_h$ satisfies $\|(u-u_h)(t)\|_{L^2(\Omega)} \lesssim h^2(t^{-\alpha}+\ln(1/h))\|u_0\|_{L^2(\Omega)}$ with $u_0\in L^2(\Omega)$. It also proves Sobolev-norm error estimates in ${}_0H^{\alpha/2}(0,T;L^2(\Omega))$, $L^2(0,T;\dot H^1(\Omega))$, and $L^2(0,T;L^2(\Omega))$, with rates that change at the thresholds $\alpha=1/3$ and $\alpha=1/2$ and are optimal with respect to the regularity of the solution. For the full discretization, Theorem 5.2 gives $O(J^{-(1-\alpha)/2}\sqrt{\ln J})$ convergence for $L^2$ data on graded grids and sharper rates, expressed through $\eta_1(\alpha,\sigma,J)$ and $\eta_2(\alpha,\sigma,J)$, for data in $\dot H^1(\Omega)$. The proof isolates the nonsmooth-data difficulty in the linear operator $S$ and its discrete counterpart $S_h$, controls the difference $u-\tilde u_h$ by projection estimates, and then uses the new Grönwall inequality to control the nonlinear feedback $u_h-\tilde u_h$.

Load-bearing premise

The load-bearing premise is the unproved spatial projection estimate $\|(S D^\alpha_{0+}v - S_h D^\alpha_{0+}P_h v)(t)\|_{L^2(\Omega)} \lesssim h^2 t^{-\alpha}\|v\|_{L^2(\Omega)}$ for all $v\in L^2(\Omega)$; if that transfer from a known linear result fails, the central pointwise error bound of Theorem 4.1 does not follow.

Editorial extensions

If this is right

  • Piecewise-linear finite element approximations of semilinear fractional diffusion converge in $L^2(\Omega)$ at the rate $h^2$ away from the initial singularity, even when $u_0$ is merely square-integrable.
  • Near $t=0$ the spatial error is instead $O(h^2 t^{-\alpha})$, so the estimate quantifies exactly how much the initial layer costs and where graded or adaptive meshes are needed.
  • The full discontinuous Galerkin time discretization on graded grids is guaranteed to converge at rate $J^{-(1-\alpha)/2}\sqrt{\ln J}$ in the fractional and $H^1$ norms for $L^2$ data, with the better graded-grid rates of Theorem 5.2 for $H^1$ data.
  • The new Grönwall inequality and its discrete version are reusable tools: any fractional evolution equation whose energy identity yields the same fractional-derivative bound now has a route from that bound to a convergence theorem.
  • The earlier smooth-data analyses become a special case of the same argument, and the new proof does not require an extra square of the logarithmic factor in the pointwise rough-data estimate.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The $t^{-\alpha}$ factor in the pointwise bound suggests that on a fixed spatial mesh the early-time error dominates; a time-adaptive finite element code that refines the first moments more finely should be able to recover uniform $O(h^2)$ accuracy in $t$, a design question the paper does not address.
  • The abstract Grönwall lemma is not tied to the specific nonlinearity $f(u)$; the same two-step semidiscrete argument should carry over to locally Lipschitz nonlinearities, time-dependent coefficients, or systems of fractional diffusion equations, giving testable analogues of Theorems 4.1 and 5.2.
  • If the omitted projection estimate, Eq. (48), is proved for higher-degree finite elements, the identical architecture should yield $O(h^{r+1}(t^{-\alpha}+\ln(1/h)))$ pointwise accuracy for $P_r$ elements with only cosmetic changes.
  • The numerical results for small $\alpha$ (for example, $\alpha=0.2$ in Table 5) show rates below the predicted ones, and the paper attributes this to instability at large $J$; a clean way to separate a proof gap from a numerical artifact is to rerun those cases with a stable solver on substantially larger graded grids and check whether $E_1$ approaches $J^{-1/2}$.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. This paper develops a priori error analysis for the semilinear fractional diffusion equation D^alpha_{0+}(u-u0) - Delta u = f(u) on a convex polytope, with alpha in (0,1) and u0 in L2(Omega). The authors prove existence, uniqueness, and temporal regularity of the weak solution, introduce a new continuous Gronwall inequality and a discrete analogue, and use them to derive spatial finite-element error estimates for a semidiscrete scheme (Theorem 4.1) and temporal error estimates for a piecewise-constant discontinuous Galerkin discretization on graded grids (Theorem 5.2). The headline results are the pointwise L2 spatial error O(h^2(t^{-alpha}+ln(1/h))) and temporal rates O(J^{-(1-alpha)/2} sqrt(ln J)) for nonsmooth data. Numerical experiments in one dimension verify the predicted orders.

Significance. If the proofs are completed, the paper would supply one of the first numerical analyses for a semilinear fractional diffusion problem with nonsmooth initial data, extending linear subdiffusion theory to the semilinear setting. The proposed Gronwall-type inequalities are a useful technical contribution, the spatial rates (42)-(44) are sharp with respect to the solution regularity, and the graded-grid temporal rates in Theorem 5.2 go beyond uniform-grid analyses. The numerical section is honest in that errors are measured against a fine reference solution rather than calibrated to the theory, and the experiments reproduce the predicted rates in most regimes, with the authors themselves noting the discrepancy for alpha=0.2 in Table 5. However, two load-bearing estimates are not fully proved: Lemma 4.1(48) and an inequality appearing in Lemma 3.5. Until those arguments are supplied, the central spatial and full-discretization claims are conditionally established rather than rigorously proven.

major comments (3)
  1. [Section 4, Lemma 4.1, Eq. (48)] The estimate (48) is the only ingredient that produces the O(h^2 t^{-alpha}) term in the pointwise error (42), since Lemma 4.2 applies it with v=u0 and Theorem 4.1 uses (51). Its proof is omitted ('similar to that of [17, Theorem 2.1]'), and the only supporting citation given in Remark 4.2, [5, Theorem 3.7], provides the strictly weaker bound O(h^2 ln(1/h) t^{-alpha}). If (48) holds only with the logarithmic factor, then the headline bound (42) and the spatial rates (43)-(44) degrade unless the theorems are restated. Please supply a complete proof of (48), or weaken the statement and propagate the change through Lemma 4.2 and Theorem 4.1.
  2. [Section 3, Lemma 3.5] In the proof of Lemma 3.5, the displayed chain of inequalities contains the step ||D^{-beta}_{t-}y||_{L2(0,t)} <= C_beta ||D^{-beta}_{0+}y||_{L2(0,t)}, justified by 'Lemma 2.2'. Lemma 2.2, however, concerns positive-order fractional derivatives and does not, as stated, imply a comparison of negative-order right- and left-fractional integrals. This inequality is load-bearing: Lemma 3.5 is invoked in the uniqueness proof in Theorem 3.1 and in the derivation of (54)-(55) in Theorem 4.1, and its discrete version underlies Lemma 5.1 and therefore (64)-(65). Please add a proof of this inequality or replace Lemma 3.5 with a version whose every step is justified.
  3. [Sections 3-5, delegated proofs] Several results that are used as theorems are not proved in the text: the cases 0<alpha<=1/2 of Theorem 3.1 ('the rest... can be proved analogously'), the semidiscrete regularity bound (41) ('Following the proof of Theorems 3.1 and 3.2, we can easily prove'), the proof of (83) in Lemma B.1 ('Since the proof of (83) is similar, it is omitted'), and estimate (64) in Theorem 5.2 ('Since (64) can be proved analogously'). The omitted cases carry the spatial rates for alpha<=1/2 in (43)-(44) and the temporal rate (64), so these are not purely presentational omissions. Please expand these proofs or state explicitly which existing argument is being modified, where, and why the modifications do not change the estimates.
minor comments (5)
  1. [Section 2, conventions] The text 'l = 1, 2.3, 4' should read 'l = 1, 2, 3, 4'.
  2. [Section 6, Tables 1-3] The headers 'J = 2^{-16}' and the reference-solution description should involve J = 2^{16} (or 65536), since J is the number of temporal intervals and must be a positive integer.
  3. [Section 5, Lemma 5.1] The statement of Lemma 5.1 has the exponent (tau^* - tau)^{-1/2} in the exponential, while the proof yields (tau^* - tau)^{-1}; please reconcile the statement with the proof.
  4. [Section 3, proof of Lemma 3.5] The notation '2kbeta' is ambiguous; it should be written as '2^k beta' to distinguish powers from products.
  5. [Abstract and Introduction] The phrase 'optimal with respect to the regularity of the solution' is used without a definition of optimality or any lower-bound discussion; consider rephrasing to 'sharp up to logarithmic factors' where that is what is proved.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity; the error estimates are derived from stated assumptions and external operator estimates, not from fitted data or self-referential definitions.

full rationale

The derivation chain is not circular. The spatial semidiscrete bound (42) follows from the auxiliary comparison function ~u_h defined in (50), Lemma 4.2, and a fractional Grönwall inequality, none of which encodes the target error as an input. The temporal estimates (64) and (65) follow from approximation properties of the projection P_tau and the semidiscrete regularity estimate (41), again without calibrating any unknown constant against the numerical results. The numerical experiments are verifications, not fitting sources. The paper does cite prior works by its own authors (e.g., [9], [10], [18]) for technical building blocks, but those are independent published estimates used as lemmas, not assumptions equivalent to the conclusions. The one load-bearing gap is Lemma 4.1, Eq. (48): the proof is omitted with the note "The proof of (48) is similar to that of [17, Theorem 2.1] and so is omitted," and Remark 4.2 concedes that the cited [5, Theorem 3.7] yields only the weaker h^2 ln(1/h) t^{-alpha} bound. This is a genuine correctness risk because (48) is stronger than what Remark 4.2 attributes to the cited result, but it is a missing proof or unverified external transfer, not a circular reduction: (48) is not shown to be equivalent to its own input, and nothing in the paper defines or fits the target error in terms of itself. Accordingly, no circular step can be exhibited under the required standard.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No free parameters are fitted to data; the analysis is a priori. The paper borrows several operator regularity and spatial approximation results from prior work, and those are treated as explicit assumptions. No new physical or mathematical entities are postulated.

assumptions (4)
  • standard math Standard fractional calculus, Sobolev space interpolation, and spectral properties of the Dirichlet Laplacian are used without proof.
    Background for the operator S, fractional Sobolev spaces, and the energy estimates in Sections 2 and 3.
  • domain assumption The solution operator S is well defined and satisfies the smoothing estimate (4), taken from references [12, 9, 18].
    Invoked in Section 3 to define weak solutions and to prove regularity estimates.
  • domain assumption The spatial error estimate ||(S D^alpha_{0+} v - S_h D^alpha_{0+} P_h v)(t)|| is bounded by C h^2 t^{-alpha} ||v|| for v in L2(Omega), as stated in Lemma 4.1 (48) without proof.
    This result, attributed to [17, Theorem 2.1], is the basis for the pointwise spatial error estimate (42).
  • domain assumption Omega is a convex d-polytope with d = 1, 2, 3, u0 is in L2(Omega), and f is Lipschitz continuous with f(0) = 0.
    These are the stated problem conditions in Section 1 and define the class of problems analyzed.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Numerical analysis of a semilinear fractional diffusion equation." pith.science (2026). https://pith.science/paper/5TR3PKI7

@misc{pith2026190900016,
  author       = {Pith},
  title        = {Pith review of: Numerical analysis of a semilinear fractional diffusion equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5TR3PKI7}},
  note         = {Machine review of arXiv:1909.00016}
}
abstract

This paper considers the numerical analysis of a semilinear fractional diffusion equation with nonsmooth initial data. A new Gr\"onwall's inequality and its discrete version are proposed. By the two inequalities, error estimates in three Sobolev norms are derived for a spatial semi-discretization and a full discretization, which are optimal with respect to the regularity of the solution. A sharp temporal error estimate on graded temporal grids is also rigorously established. In addition, the spatial accuracy $\scriptstyle O(h^2(t^{-\alpha} + \ln(1/h)\!)\!) $ in the pointwise $ \scriptstyle L^2(\Omega) $-norm is obtained for a spatial semi-discretization. Finally, several numerical results are provided to verify the theoretical results.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 31 canonical work pages

  1. [17]

    Lubich, I

    C. Lubich, I. Sloan, and V. Thom´ ee. Nonsmooth data error es timates for approximations of an evolution equation with a positive-type memory term. Math. Comput. , 65(213):1–17, 1996

    work page 1996

  2. [1]

    Cuesta, C

    E. Cuesta, C. Lubich, and C. Palencia. Convolution quadrature t ime discretization of fractional diffusion-wave equations. Math. Comput. , 75(254):673–696, 2006

    work page 2006

  3. [2]

    Ervin and J

    V. Ervin and J. Roop. Variational formulation for the stationary fractional advection dispersion equation. Numer. Meth. Part. D. E. , 22(3):558–576, 2006

    work page 2006

  4. [3]

    L. C. Evans. Partial differential equations . American Mathematical Society, 2 edition, 2010

    work page 2010

  5. [4]

    B. Jin, R. Lazarov, J. Pasciak, and Z. Zhou. Error analysis of se midiscrete finite element methods for inhomogeneous time-fractional diffusion . IMA. J. Numer. Anal. , 35(2):561–582, 2015

    work page 2015

  6. [5]

    B. Jin, R. Lazarov, and Z. Zhou. Error estimates for a semidiscr ete finite element method for fractional order parabolic equations. SIAM J. Numer. Anal., 51(1):445–466, 2013

    work page 2013

  7. [6]

    B. Jin, R. Lazarov, and Z. Zhou. Two fully discrete schemes for f ractional diffusion and diffusion-wave equations with nonsmooth data. SIAM J. Sci Comput., 38(1):A146–A170, 2016

    work page 2016

  8. [7]

    B. Jin, B. Li, and Z. Zhou. Numerical analysis of nonlinear subdiffus ion equations. SIAM J. Numer. Anal. , 56(1):1–23, 2018

    work page 2018

Show all 32 references
  1. [8]

    B. Li, H. Luo, and X. Xie. A time-spectral algorithm for fractiona l wave problems. J. Sci. Comput. , 77(2):1164–1184, 2018. 28

    2018

  2. [9]

    B. Li, H. Luo, and X. Xie. Analysis of a time-stepping scheme for tim e fractional diffusion problems with nonsmooth data. SIAM J. Numer. Anal. , 57(2):779–798, 2019

    2019

  3. [10]

    B. Li, T. Wang, and X. Xie. Analysis of the L1 scheme for fraction al wave equations with nonsmooth data. submitted, axXiv:1908.09145

    1908 arXiv

  4. [11]

    D. Li, C. Wu, and Z. Zhang. Linearized galerkin fems for nonlinear time fractional parabolic problems with non-smooth solutions in time direc tion. J. Sci. Comput. , 80(1):403–419, 2019

    2019

  5. [12]

    Li and C

    X. Li and C. Xu. A space-time spectral method for the time fra ctional diffusion equation. SIAM J. Numer. Anal. , 47(3):2108–2131, 2009

    2009

  6. [13]

    H. Liao, W. McLean and J. Zhang. A discrete Gr¨ onwall inequality with ap- plication to numerical schemes for subdiffusion problems SIAM J. Numer. Anal., 57(1):218–237, 2019

    2019

  7. [14]

    Lin and C

    Y. Lin and C. Xu. Finite difference/spectral approximations for the time- fractional diffusion equation. J. Comput. Phys. , 225(2):1533–1552, 2007

    2007

  8. [15]

    C. Lubich. Discretized fractional calculus. SIAM J. Math. Anal. , 17(3):704– 719, 1986

    1986

  9. [16]

    C. Lubich. Convolution quadrature and discretized operationa l calculus. Numer. Math. , 52(4):129–145, 1988

    1988

  10. [18]

    H. Luo, B. Li., and X. Xie. Convergence analysis of a petrov–gale rkin method for fractional wave problems with nonsmooth data. J. Sci. Com- put., 80(2):957–992, 2019

    2019

  11. [19]

    McLean and K

    W. McLean and K. Mustapha. Convergence analysis of a discont inuous galerkin method for a sub-diffusion equation. Numer. Algor. , 52(1):69–88, 2009

    2009

  12. [20]

    McLean and K

    W. McLean and K. Mustapha. Time-stepping error bounds for f rac- tional diffusion problems with non-smooth initial data. J. Comput. Phys. , 293:201–217, 2015

    2015

  13. [21]

    Mustapha and W

    K. Mustapha and W. McLean. Discontinuous Galerkin method for an evolution equation with a memory term of positive type. Math. Comput. , 78(268):1975–1995, 2009

    1975

  14. [22]

    Mustapha and W

    K. Mustapha and W. McLean. Uniform convergence for a discon tinuous Galerkin, time-stepping method applied to a fractional diffusion equa tion. IMA J. Numer. Anal. , 32(3):906–925(20), 2012

    2012

  15. [23]

    Mustapha and H

    K. Mustapha and H. Mustapha. A second-order accurate num erical method for a semilinear integro-differential equations with a weakly singular k ernel. IMA J. Numer. Anal. , 30:555–578, 2010. 29

    2010

  16. [24]

    Podlubny

    I. Podlubny. Fractional differential equations . Academic Press, 1998

    1998

  17. [25]

    Sun and X

    Z. Sun and X. Wu. A fully discrete difference scheme for a diffusion -wave system. Appl. Numer. Math. , 56(2):193–209, 2006

    2006

  18. [26]

    L. Tartar. An Introduction to Sobolev Spaces and Interpolation Spaces . Springer, Berlin, 2007

    2007

  19. [27]

    R. Wang, D. Chen, and T. Xiao. Abstract fractional cauchy pr oblems with almost sectorial operators. J. Differ. Equations , 252(1):202–235, 2012

    2012

  20. [28]

    A. Yagi. Abstract parabolic evolution equations and their applica tions. Springer, Berlin, 2010

    2010

  21. [29]

    Zayernouri and G

    M. Zayernouri and G. E. Karniadakis. Fractional spectral co llocation method. SIAM J. Sci. Comput. , 36(1):A40–A62, 2014

    2014

  22. [30]

    Zayernouri and G

    M. Zayernouri and G. E. Karniadakis. Exponentially accurate s pectral and spectral element methods for fractional odes. J. Comput. Phys. , 257(2):460– 480, 2014

    2014

  23. [31]

    Zayernouri and G

    M. Zayernouri and G. E. Karniadakis. Discontinuous spectral element methods for time- and space-fractional advection equations. SIAM J. Sci. Comput., 36(4):B684–B707, 2014

    2014

  24. [32]

    Zheng, F

    M. Zheng, F. Liu, I. Turner, and V. Anh. A novel high order spa ce-time method for the time fractional fokker-planck equation. SIAM J. Sci. Com- put., 37(2):A701–A724, 2015. 30

    2015

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.