pith. sign in

arxiv: 2306.10966 · v3 · submitted 2023-06-19 · 🧮 math.NA · cs.NA

Overcoming the order barrier two in splitting methods when applied to semilinear parabolic problems with non-periodic boundary conditions

Ramona H\"aberli This is my paper

Pith reviewed 2026-05-24 08:08 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords splitting methodsorder reductionsemilinear parabolic problemsnon-periodic boundary conditionstime integrationconvergence analysisnumerical experiments
0
0 comments X

The pith

A third-order splitting method for semilinear parabolic problems avoids order reduction with non-periodic boundary conditions.

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

The paper introduces a splitting method of order three for semilinear parabolic problems that prevents the order reduction typically seen in high-order splitting methods on PDEs with non-periodic boundary conditions. It proves third-order convergence in a simplified linear setting and confirms the result through numerical experiments. The work further shows numerically that an order-four variant and cases with nonlinear source terms also retain high-order convergence.

Core claim

We introduce a splitting method of order three for a class of semilinear parabolic problems that avoids order reduction in the context of non-periodic boundary conditions. We give a proof for the third order convergence of the method in a simplified linear setting and confirm the result by numerical experiments. Moreover, we show numerically that the high order convergence persists for an order four variant of a splitting method, and also for a nonlinear source term.

What carries the argument

The order-three splitting method constructed via corrector techniques to counteract boundary-induced order reduction.

If this is right

  • The method attains third-order accuracy on the target problems without the usual order drop.
  • Numerical results indicate the same high-order behavior carries over to fourth-order splitting variants.
  • High-order convergence holds when the source term is nonlinear.
  • The approach addresses the order barrier of two for splitting methods under non-periodic conditions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The corrector technique could be tested on other classes of evolution equations that exhibit order reduction.
  • Implementation in existing codes for parabolic problems would allow direct comparison of wall-clock time versus accuracy against lower-order methods.
  • Extension of the convergence proof from the linear case to the full nonlinear setting remains an open step suggested by the numerics.

Load-bearing premise

The convergence properties seen in the simplified linear setting and numerical experiments extend to the full class of semilinear parabolic problems with non-periodic boundary conditions without further order reduction.

What would settle it

A computation on the full nonlinear semilinear parabolic problem with non-periodic boundary conditions that yields observed order below three would falsify the central claim.

Figures

Figures reproduced from arXiv: 2306.10966 by Ramona H\"aberli.

Figure 1
Figure 1. Figure 1: Convergence error of the corrected method [PITH_FULL_IMAGE:figures/full_fig_p014_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Convergence error of the corrected method [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Convergence error of the splitting methods [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Convergence error of the splitting methods [PITH_FULL_IMAGE:figures/full_fig_p016_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: The pointwise error of the methods C3Naiv (left) and C3New (right) applied to problem (5.3) in the domain Ω = (0, 1)2 for time step τ = 10−2 . The error of the naiv splitting is concentrated on the boundary ∂Ω, where the error of the new method vanishes. biological feature are considered. The solution u represents the cell density. We compare in [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Convergence error of the corrected method [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
read the original abstract

In general, high order splitting methods suffer from an order reduction phenomena when applied to the time integration of partial differential equations with non-periodic boundary conditions. In the last decade, there were introduced several modifications to prevent the second order Strang splitting method from such a phenomena. In this article, inspired by these recent corrector techniques, we introduce a splitting method of order three for a class of semilinear parabolic problems that avoids order reduction in the context of non-periodic boundary conditions. We give a proof for the third order convergence of the method in a simplified linear setting and confirm the result by numerical experiments. Moreover, we show numerically that the high order convergence persists for an order four variant of a splitting method, and also for a nonlinear source term.

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 / 1 minor

Summary. The manuscript introduces a corrector-based splitting method of order three for semilinear parabolic problems with non-periodic boundary conditions to overcome the order-reduction barrier that typically affects high-order splitting methods. It supplies a rigorous proof of third-order convergence only in a simplified linear setting and relies on numerical experiments to support the claims for the full semilinear case, an order-four variant, and nonlinear source terms.

Significance. If the method attains the claimed order without reduction for the full problem class, the work would usefully extend existing corrector techniques (previously applied to Strang splitting) to higher orders. This could benefit applications requiring accurate time integration of parabolic PDEs on non-periodic domains. The provision of machine-checked elements is not mentioned, but the combination of a partial proof with targeted numerics is a standard approach in the field.

major comments (2)
  1. [Abstract] Abstract: the central claim is that the new splitting method 'avoids order reduction' for the class of semilinear parabolic problems with non-periodic boundary conditions. However, the convergence proof is supplied only for a simplified linear setting; no error expansion or stability argument is given that controls the additional commutator and boundary terms generated by the nonlinear source term. If those terms produce an uncancelled O(τ²) contribution, the headline result does not hold.
  2. [Numerical experiments] The numerical confirmation of order three (and persistence for the order-four variant) is presented without a supporting consistency analysis for the nonlinear case under the given boundary conditions. This leaves open whether the chosen test problems are representative of the general setting where order reduction can appear.
minor comments (1)
  1. The notation used for the splitting operators and the corrector could be introduced more explicitly at the beginning of the method description to aid readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback. We respond point-by-point to the major comments below, clarifying the scope of our results and proposing targeted revisions to avoid any overstatement.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim is that the new splitting method 'avoids order reduction' for the class of semilinear parabolic problems with non-periodic boundary conditions. However, the convergence proof is supplied only for a simplified linear setting; no error expansion or stability argument is given that controls the additional commutator and boundary terms generated by the nonlinear source term. If those terms produce an uncancelled O(τ²) contribution, the headline result does not hold.

    Authors: We agree that the rigorous third-order convergence proof is supplied only for the linear case. The abstract already states explicitly that the proof applies to a simplified linear setting and that the result for the semilinear problem is confirmed numerically. The corrector is constructed precisely to cancel the leading boundary-induced error terms that cause order reduction; the same cancellation mechanism is expected to control the additional commutators arising from the nonlinear source. Nevertheless, we acknowledge that an explicit error expansion for the nonlinear case is absent. We will revise the abstract and the concluding section to state more precisely that avoidance of order reduction for the full semilinear class is supported by numerical evidence rather than by a complete proof, and we will add a short remark outlining why the nonlinear commutator terms are anticipated to remain higher order under the given boundary conditions. revision: partial

  2. Referee: [Numerical experiments] The numerical confirmation of order three (and persistence for the order-four variant) is presented without a supporting consistency analysis for the nonlinear case under the given boundary conditions. This leaves open whether the chosen test problems are representative of the general setting where order reduction can appear.

    Authors: The test problems were deliberately chosen to reproduce the classic order-reduction scenario for splitting methods on non-periodic domains (non-homogeneous Dirichlet boundaries that generate the problematic boundary layers). Both linear and nonlinear source terms are included, and the observed convergence rates remain third order (and fourth order for the higher variant). While a full consistency analysis for the nonlinear case is not provided, the experiments are representative of the setting in which order reduction is known to occur. To strengthen the presentation we will expand the description of the test-problem selection criteria and add one further numerical example with a different nonlinearity to illustrate robustness. revision: partial

Circularity Check

0 steps flagged

No circularity; proof and numerics are independent of inputs.

full rationale

The manuscript introduces a new corrector-based splitting method and supplies an explicit convergence proof for order 3 in a simplified linear setting. Extension to the semilinear non-periodic case and to an order-4 variant rests on numerical experiments rather than on any fitted parameter, self-referential definition, or load-bearing self-citation. No equation or claim reduces by construction to its own inputs; the derivation chain therefore remains non-circular.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper introduces a modified splitting method but details on any new parameters or entities are not available from the abstract.

axioms (1)
  • domain assumption Convergence analysis assumptions for splitting methods in linear parabolic PDEs with non-periodic BCs
    The proof is given in a simplified linear setting as stated in the abstract.

pith-pipeline@v0.9.0 · 5657 in / 1142 out tokens · 32897 ms · 2026-05-24T08:08:43.532091+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

27 extracted references · 27 canonical work pages

  1. [1]

    Achouri, M

    T. Achouri, M. Ayadi, A. Habbal, and B. Yahyaoui. Numerical analysis for the two-dimensional fisher- kolmogorov-petrovski-piskunov equation with mixed boundary conditions. J. Appl. Math. and Comput., 68(6):3589–3614, 2022

    work page 2022

  2. [2]

    Bader, S

    P. Bader, S. Blanes, and F. Casas. Solving the schr ¨odinger eigenvalue problem by imaginary time propa- gation techniques using splitting methods with complex coefficients. J. Chem. Phys., 139(124117):1–11, 2013

    work page 2013

  3. [3]

    Bertoli and G

    G. Bertoli and G. Vilmart. Strang splitting method for semilinear parabolic problems with inhomoge- neous boundary conditions: a correction based on the flow of the nonlinearity. SIAM J. Sci. Comput. , 42(3):A1913–A1934, 2020

    work page 2020

  4. [4]

    Blanes and F

    S. Blanes and F. Casas. On the necessity of negative coefficients for operator splitting schemes of order higher than two. Appl. Num. Math., 54(1):23–37, 2005

    work page 2005

  5. [5]

    Blanes, F

    S. Blanes, F. Casas, P. Chartier, and A. Murua. High order splitting methods for some classes of parabolic equations. Math. of Comput., 82(283):1559–1576, 2013

    work page 2013

  6. [6]

    Blanes, F

    S. Blanes, F. Casas, and A. Murua. Splitting methods. J. Comput. Appl. Math., 142(2):313–330, 2013

    work page 2013

  7. [7]

    Castella, P

    F. Castella, P. Chartier, S. Decombes, and G. Vilmart. Splitting methods with complex times for parabolic equations. Bit Num. Math., 49:487–508, 2009. 18

    work page 2009

  8. [8]

    Descombes

    S. Descombes. Convergence of a splitting method of high order for reaction-diffusion systems. Math. of Comput., 70(236):1481–1501, 2001

    work page 2001

  9. [9]

    Einkemmer, M

    L. Einkemmer, M. Moccaldi, and A. Ostermann. Efficient boundary corrected Strang splitting. Appl. Math. Comput., 332:76–89, 2018

    work page 2018

  10. [10]

    Einkemmer and A

    L. Einkemmer and A. Ostermann. Overcoming order reduction in diffusion-reaction splitting. part 1: Dirichlet boundary conditions. J. on Sc. Comput., 37(3):A1577–A1592, 2015

    work page 2015

  11. [11]

    Einkemmer and A

    L. Einkemmer and A. Ostermann. Overcoming order reduction in diffusion-reaction splitting. part 2: Oblique boundary conditions. J. on Sc. Comput., 38(6):A3741–A3757, 2016

    work page 2016

  12. [12]

    L. C. Evans. Partial Differential Equations, volume 19 of Graduate Studies in Mathematics. American Mathematical Society, 1997

    work page 1997

  13. [13]

    Gerisch and J

    A. Gerisch and J. G. Verwer. Operator splitting and approximate factorization for taxis-diffusion-reaction models. Arch. Rational Mech. Anal., 42(1-3):159–176, 2002

    work page 2002

  14. [14]

    Goldman and T

    D. Goldman and T. J. Kaper. N-th order operator splitting schemes and nonreversible systems. SIAM J. Numer. Anal., 33(1):349–367, 1996

    work page 1996

  15. [15]

    Hairer, C

    E. Hairer, C. Lubich, and G. Wanner. Geometric numerical integration, volume 31 of Springer Series in Computational Mathematics. Springer, Heidelberg, 2010. Structure-preserving algorithms for ordinary differential equations, Reprint of the second (2006) edition

    work page 2010

  16. [16]

    Hansen, F

    E. Hansen, F. Kramer, and A. Ostermann. A second-order positivity preserving scheme for semilinear parabolic problems. Appl. Num. Math., 62(10):1428–1435, 2012

    work page 2012

  17. [17]

    Hansen and A

    E. Hansen and A. Ostermann. High order splitting methods for analytic semigroups exist. Bit Num. Math., 49:527–542, 2009

    work page 2009

  18. [18]

    Hundsdorfer and J

    W. Hundsdorfer and J. Verwer. Numerical solution of time-dependent advection-diffusion-reaction equa- tions, volume 33 of Springer Series in Computational Mathematics. Springer-Verlag, Berlin, 2003

    work page 2003

  19. [19]

    Lorenzi, A

    L. Lorenzi, A. Lunardi, G. Metafune, and D. Pallara. Analytic semigroups and reaction-diffusion prob- lems. Internet Seminar, 2005

    work page 2005

  20. [20]

    A. Lunardi. Analytic semigroups and optimal regularity in parabolic problems . Modern Birkh ¨auser Classics. Birkh¨auser/Springer Basel AG, Basel, 1995

    work page 1995

  21. [21]

    R. I. McLachlan and G. R. W. Quispel. Splitting methods. Acta Numerica, pages 341–434, 2002

    work page 2002

  22. [22]

    J. Ne ˘cas. Direct Methods in the Theory of Elliptic Equations . Springer Monographs in Mathematics. Springer Berlin, Heidelberg, 1967

    work page 1967

  23. [23]

    Niesen and W

    J. Niesen and W. M. Wright. Algorithm 919: a Krylov subspace algorithm for evaluating the ϕ-functions appearing in exponential integrators. ACM Trans. Math. Software, 38(3):Art. 22, 19, 2012. 19

    work page 2012

  24. [24]

    M. J. Simpson and K. A. Landmann. Characterizing and minimizing the operator split error for fisher’s equation. Appl. Math. Letters 19, pages 604–612, 2006

    work page 2006

  25. [25]

    M. Suzuki. Fractial decomposition of exponential operators with applications to many-body theories and monte carlo simulations. Physics Letters A, 146(6):319–323, 1990

    work page 1990

  26. [26]

    V . Thom´ee. Galerkin finite element methods for parabolic problems , volume 25 of Springer Series in Computational Mathematics. Springer-Verlag, Berlin, second edition, 2006

    work page 2006

  27. [27]

    H. Yoshida. Construction of higher order symplectic operators. Physics Letters A , 150(5-7):262–268, 1990. 20

    work page 1990