Pith. sign in

REVIEW 3 major objections 3 minor 77 references

Exponential Time Differencing Runge-Kutta Discontinuous Galerkin (ETD-RKDG) Methods for Nonlinear Degenerate Parabolic Equations

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

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

Pith's one-line read This paper claims that ETD-RKDG methods resolve the stiffness of nonlinear degenerate parabolic equations, allowing time steps far larger than the explicit O(h^2) restriction while keeping high-order accuracy on unstructured simplex meshes.

desk verdict A solid engineering extension of ETD-RK to nodal DG for degenerate parabolic problems, with an honest but unproven stability heuristic; needs a revision with real comparisons and parameter reporting. read the letter →

arxiv 2506.04416 v1 pith:Q254JCQQ submitted 2025-06-04 math.NA cs.NA

Ziyao Xu , Yong-Tao Zhang This is my paper
classification math.NAcs.NA MSC 65M6065M2065L0635K65
keywords exponentialtimedifferencingRunge-KuttadiscontinuousGalerkinnodalformulationdegenerateparabolicequationsporousmediumequationRosenbrock-typetreatmentunstructuredsimplexmeshes
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

Nonlinear degenerate parabolic equations, such as the porous medium equation, develop sharp fronts where diffusion vanishes, while the parabolic part still forces explicit time steps of order h squared. The paper develops high-order exponential time differencing Runge-Kutta discontinuous Galerkin (ETD-RKDG) methods that absorb the Jacobian of the diffusion term into an exponential integrator and integrate the remaining nonlinearity explicitly. On triangular and tetrahedral meshes, the method captures fronts with limiters and allows much larger time steps while keeping high-order accuracy. A one-dimensional Fourier stability analysis identifies how much of the diffusion operator must be absorbed, and numerical tests in two and three dimensions, including a torus, support the claimed stability and efficiency.

What carries the argument

The central object is a Jacobian-based Rosenbrock-type ETD-RK splitting on a nodal DG discretization. The paper writes the semi-discrete ODE as $u_t = D(u) + A(u)$, sets $L = D'(u_n)$, and uses the $\phi$-functions of the exponential integrator to absorb $L$ exactly while integrating the residual $N(u)$ explicitly. The nodal formulation, with Lagrange basis functions at Legendre-Gauss-Lobatto nodes, makes the Jacobian computation vectorized and sparse via $L = G .* g'(u)^T + (\beta + \alpha/2) P$; the fully discrete stability analysis tracks the spectral radius of the Fourier growth factor.

What would settle it

Run the fully discrete growth-factor check $\rho(\hat{G}(\theta,h,\xi))$ for $\theta$ just above each $\theta_0$, for polynomial degrees $k=1,\ldots,5$ and mesh sizes spanning several orders of magnitude, with $\xi$ sampled much more densely than in the paper; finding any value above $1$ would disprove the stability claim. Alternatively, integrate the porous medium equation with ETD-RK3 and $P^2$ elements at time steps near the claimed limit and observe whether the solution remains bounded and the wave front is captured without growing oscillations.

Watch

Extended reading notes

Core claim

The paper's central claim is that the stiffness of a nonlinear degenerate parabolic equation can be resolved by splitting the semi-discrete discontinuous Galerkin system as $u_t = D(u) + A(u)$, taking $L = D'(u_n)$ as the linear operator absorbed by the exponential integrator, and treating $N(u) = D(u) + A(u) - L u$ explicitly. Because the Jacobian of the diffusion discretization is absorbed, the exponential decay of high-frequency modes is handled exactly, bypassing the $\tau \sim O(h^2)$ restriction. For the split model $u_t = a_0 u_{xx} + (a-a_0)u_{xx}$, the paper proves or argues that the ETD-RK schemes are stable whenever $a_0 \ge \theta_0 a$, with thresholds $1/2$, $1/2$, about $0.6034$, and $1/2$ for ETD-RK1 through ETD-RK4. A numerical search over mesh size and polynomial degree shows the same thresholds for the fully discrete scheme. The paper further claims that this stability is independent of the spatial discretization parameters, and that the method preserves high-order accuracy while allowing much larger time steps than explicit SSP-RK methods.

Load-bearing premise

The load-bearing premise is that the sufficiency of the threshold $\theta \ge \theta_0$, inferred from graphs of the growth factor and from an extensive numerical search over mesh sizes and polynomial degrees, covers every parameter regime. If some unexamined combination violates the stability bound, the claimed unconditional stability and mesh-independence are not established.

Editorial extensions

If this is right

  • For linear diffusion problems, the scheme is unconditionally stable in the weak sense $\tau \le C$ with $C$ independent of the spatial mesh, removing the parabolic time-step restriction.
  • For nonlinear degenerate problems, absorbing at least a fraction $\theta_0$ of the diffusion Jacobian gives the same qualitative stability, allowing time steps far larger than the explicit $\tau \sim O(h^2)$ limit.
  • High-order accuracy is preserved with large time steps, as shown by the convergence tables for $P^2$ and $P^3$ spaces with $\tau$ proportional to $h$ or $0.2h$.
  • The method works on unstructured triangular and tetrahedral meshes, enabling simulations on complex domains such as a three-dimensional torus.
  • Adding a reaction term to the equation is straightforward in the nodal formulation, so the method extends to convection-diffusion-reaction problems without new algorithmic machinery.

Reading between the lines

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

  • A natural extension would use the threshold $\theta_0$ as a practical guide for choosing the background diffusion coefficient $a_0$ adaptively at each time step to minimize numerical cost while preserving stability.
  • The graph-based sufficiency argument for the thresholds could probably be made rigorous by bounding the derivative of the growth factor, or it could fail in an unexpected regime; a dense random search over $\theta$, $h$, and $\xi$ would settle the question.
  • The method's structure should carry over to other degenerate equations such as thin-film or Stefan-type problems, but positivity preservation near the front, which the paper handles only through a limiter, would need separate attention.
  • On very large meshes, the Krylov subspace evaluation of $\phi$-functions may become the dominant cost, so the practical efficiency gain over implicit methods will depend on the linear algebra solver as much as on the relaxed time-step restriction.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 3 minor

Summary. The paper develops a class of exponential time differencing Runge-Kutta discontinuous Galerkin (ETD-RKDG) methods for nonlinear degenerate parabolic equations. The spatial discretization uses a nodal DG formulation on simplex meshes, with an interior-penalty treatment of the diffusion term and a Lax-Friedrichs flux for convection; the time integration uses ETD-RK1 through ETD-RK4 with a Rosenbrock-type splitting in which the Jacobian of the DG diffusion operator is absorbed into the exponential integrator and the residual is treated explicitly. A vectorized matrix formulation and the associated Jacobian computation are described in detail for triangular and tetrahedral meshes. The paper includes a one-dimensional Fourier stability analysis of the semi-discrete split equation, a numerical study of fully discrete matrix growth factors, and numerical tests including smooth accuracy problems, Barenblatt solutions in two and three dimensions, and the porous medium equation on a torus. The central claim is that the Jacobian-based Rosenbrock-type treatment resolves the stiffness of the degenerate parabolic term and permits time steps much larger than the explicit restriction tau ~ O(h^2).

Significance. If the stability claims hold, this is a practically useful contribution: it extends exponential integrator ideas to discontinuous Galerkin methods for degenerate parabolic problems on unstructured meshes, with an explicit and implementable Jacobian assembly procedure. The accuracy tables show design-order convergence with tau = h and tau = 0.2h, which is credible evidence that the O(h^2) step restriction is avoided at least for smooth test problems, and the Barenblatt and torus simulations demonstrate qualitative behavior of the method. The stability thresholds in Section 3 are concrete and falsifiable. However, the higher-order stability sufficiency is asserted from graphs rather than proved, and the fully discrete mesh-independence conclusion rests on an underspecified numerical search, so the theoretical component needs substantial strengthening before the main stability claim can be accepted as stated.

major comments (3)
  1. [Section 3.2, Eq. (3.14)] For ETD-RK2, ETD-RK3, and ETD-RK4, the paper derives only a necessary condition from the limits at xi = infinity and then states that the condition theta >= theta0 is also sufficient 'from the graphs' of the growth factor, because the maximum of |G(theta, xi)| is always attained at xi = 0 or xi = infinity. This is load-bearing for the claimed unconditional stability: the growth factor is a transcendental function of xi, and a graphical observation cannot exclude an intermediate interval where |G(theta, xi)| > 1. Please replace this assertion with an analytic proof, or with a rigorous verified-interval argument covering all xi in [0, infinity) for each theta0, or explicitly reclassify the higher-order stability condition as numerical evidence and adjust the abstract and conclusions accordingly.
  2. [Section 3.3, Eq. (3.19)] The fully discrete stability conclusion — that the spectral radius of the matrix growth factor is at most 1 for all theta >= theta0, all h > 0, and all polynomial degrees k — is justified only by 'an extensive numerical search', with no reported ranges of xi, h, or k, no sampling counts, no penalty parameter beta values, and no reproducibility data. Since this is the step that yields mesh-independence of the stability condition, the claim cannot be stated as a conclusive result. Please provide either a proof of the mesh-independence or a precisely documented numerical study, and phrase the conclusion as numerical evidence rather than as an established property.
  3. [Section 4, Examples 2 and 3] The central practical claim of the paper is 'significant improvements in stability and large time-step sizes' for nonlinear degenerate parabolic equations, but the nonlinear tests in Examples 2 and 3 report only one successful time step per setup and do not measure the stability limit or compare with the explicit tau ~ O(h^2) restriction. The smooth accuracy test in Table 2 uses tau = 0.2h and is useful evidence, but the Barenblatt and torus tests do not quantify how much larger tau can actually be taken. Please add stability-limit measurements or a comparison table for the nonlinear degenerate cases, or restrict the claim accordingly.
minor comments (3)
  1. [Section 3.1] The text 'Rosenborg-type treatment' should read 'Rosenbrock-type treatment'.
  2. [Section 3.3] The numerical search for fully discrete stability does not state the value or scaling of the penalty parameter beta used; since the local matrices in Appendix C depend on beta, the claimed independence of the stability condition from h and k also needs to address beta.
  3. [Remark 2.1] Remark 2.1 refers to the prior work [72] for the weak unconditional stability of the linear convection-diffusion case; because [72] is an arXiv preprint, the relevant stability statement should either be proved in this paper or stated with sufficient detail to be checked independently.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular derivation found: the stability thresholds come from an independent growth-factor analysis, and the paper's self-citations are contextual rather than load-bearing. The main caveats are proof gaps (graph-based sufficiency, numerical search), not circularity.

full rationale

The claimed derivation chain is not circular. The ETD-RK schemes are standard exponential integrators obtained from (2.19)-(2.20), and the Rosenbrock-type splitting L=D'(u_n), N(u)=D(u)+A(u)-D'(u_n)u defines the algorithm without fitting any parameter to the quantities being predicted. Section 3 derives the growth factors (3.9)-(3.12) for the linearized split equation u_t=a0 u_xx+(a-a0)u_xx and obtains the thresholds theta0 from the xi->infinity limits in (3.13)-(3.14), with ETD-RK1's threshold proven in Theorem 3.1. No equation is defined in terms of the target stability conclusion, and no fitted value is renamed as a prediction. The numerical tests are checked against exact manufactured and Barenblatt solutions, providing independent benchmarks. The two caveats are: (i) the sufficiency of theta>=theta0 for ETD-RK2-4 is asserted 'from the graphs' in Section 3.2 rather than proved, and (ii) the fully discrete mesh-independence claim in Section 3.3 rests on 'an extensive numerical search' with no reported parameter ranges for h and k. These are gaps in proof or correctness risk, not circular reductions. Self-citations to [72] and [73] motivate the semilinearization and the semi-discrete/full-discrete alignment, but the present paper's own linear analysis and numerical experiments carry the central claim; the citations are not load-bearing in an equivalence-by-construction sense.

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

The method introduces no new physical entities. Its central claims rest on several unproved or only numerically supported ingredients: the imported DG diffusion discretization, the limiter, the transfer of semi-discrete stability to the fully discrete scheme, and the constant-coefficient linearization. Free parameters are limited to stabilization and time-step choices, but their reported values are incomplete.

free parameters (3)
  • Penalty parameter beta = O(1/h), exact values per test not reported
    Stabilizes the diffusion flux in equations (2.9) and (2.10); chosen sufficiently large by hand, and the numerical tests do not report the specific values used.
  • Time-step size tau in accuracy tests = tau = h for linear, tau = 0.2h for nonlinear
    Chosen by hand in Example 1; the claimed large-time-step advantage is demonstrated only at these fixed ratios, not at maximum stable steps.
  • Time-step size in Examples 2 and 3 = not reported
    The Barenblatt and torus simulations do not state tau, mesh size, or error magnitudes, so the stability and efficiency claims cannot be independently quantified.
assumptions (4)
  • domain assumption The ultra-weak DG diffusion discretization (2.8) from [17] is stable and convergent for the degenerate problems considered.
    The paper adopts the scheme without analysis and then builds stability claims on it; no proof is given for degenerate solutions with low regularity.
  • domain assumption The TVB limiter from [24] controls oscillations without affecting the claimed order or the stability analysis.
    The limiter is invoked to suppress Gibbs phenomena, but its interaction with the exponential integrator is not analyzed.
  • domain assumption Semi-discrete stability of ETD-RK aligns with fully discrete ETD-RKDG stability.
    Assumed following [72] and stated in Section 3.1; used to transfer conclusions from the Fourier analysis to the full scheme.
  • domain assumption The constant-coefficient linear analysis transfers to nonlinear degenerate problems through pointwise Jacobian freezing.
    The split model (3.1) assumes a0 and a are constant, while in the method a0 is derived from the Jacobian at each step; the extension to the nonlinear scheme is heuristic.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Exponential Time Differencing Runge-Kutta Discontinuous Galerkin (ETD-RKDG) Methods for Nonlinear Degenerate Parabolic Equations." pith.science (2026). https://pith.science/paper/Q254JCQQ

@misc{pith2026250604416,
  author       = {Pith},
  title        = {Pith review of: Exponential Time Differencing Runge-Kutta Discontinuous Galerkin (ETD-RKDG) Methods for Nonlinear Degenerate Parabolic Equations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Q254JCQQ}},
  note         = {Machine review of arXiv:2506.04416}
}
abstract

In this paper, we study high-order exponential time differencing Runge-Kutta (ETD-RK) discontinuous Galerkin (DG) methods for nonlinear degenerate parabolic equations. This class of equations exhibits hyperbolic behavior in degenerate regions and parabolic behavior in non-degenerate regions, resulting in sharp wave fronts in the solution profiles and a parabolic-type time-step restriction, $\tau \sim O(h^2)$, for explicit time integration. To address these challenges and solve such equations in complex domains, we employ DG methods with appropriate stabilizing limiters on unstructured meshes to capture the wave fronts and use ETD-RK methods for time integration to resolve the stiffness of parabolic terms. We extract the system's stiffness using the Jacobian matrix of the DG discretization for diffusion terms and adopt a nodal formulation to facilitate its computation. The algorithm is described in detail for two-dimensional triangular meshes. We also conduct a linear stability analysis in one spatial dimension and present computational results on three-dimensional simplex meshes, demonstrating significant improvements in stability and large time-step sizes.

Figures

Figures reproduced from arXiv: 2506.04416 by the authors.

Figure 1.1
Figure 1.1. The Barenblatt solutions (1.3) with m = 3 at t = 1 in different dimensions d. exhibits parabolic behaviors in non-degenerate regions (u > 0), with the diffusion coefficient g ′ (u) = mum−1 > 0, leading to smooth solution profiles and infinite propagation speed for perturbations. Both the hyperbolic and parabolic features of the PME pose challenges for computational methods. The sharp wave front at the boundary of th… view at source ↗
Figure 2.1
Figure 2.1. Triangulation of a local region in the computational domain Ω. The DG finite element space is defined as V k h = {v ∈ L 2 (Ω) : v|K ∈ Pk (K), ∀K ∈ T }, (2.1) where P k (K) denotes the space of polynomials of degree at most k on the element K. There are two primary choices for the basis of the finite element space V k h : V k h = M K∈T span{φ K i (x), i = 1, 2, . . . , Nk} = M K∈T span{ℓ K i (x), i = 1, 2, . . . , Nk… view at source ↗
Figure 2.2
Figure 2.2. Distribution of the two-dimensional Lagrange nodes (blue dots) on Kb. The Lagrange nodes coincide with the (k + 1)-point Legende-Gauss-Lobatto (LGL) nodes on each edge. Finally, we introduce the inner products for f, g ∈ C∞(K) and ¯f ∈ C∞(K′ ), g¯ ∈ C∞(K′′), which will be used in the DG formulation [PITH_FULL_IMAGE:figures/full_fig_p006_2_2.png] view at source ↗
Figures from the paper (7 more)
Figure 3.1
Figure 3.1. Figure 3.1: The growth factor Gb(θ0, ξ), with θ0 = 1 2 , 1 2 , 2 3 + 3 q − 1 27 + √ 78 36 + 3 q − 1 27 − √ 78 36 ≈ 0.6034, and 1 2 for ETD-RK1, ETD-RK2, ETD-RK3, and ETD-RK4 respectively. 17 [PITH_FULL_IMAGE:figures/full_fig_p017_3_1.png]
Figure 4.1
Figure 4.1. Figure 4.1: Example 1. Accuracy test. Computational mesh of level 0 (coarsest). We test different combinations of ETD-RK and DG methods on various meshes with different levels of refinement. The coarsest mesh (level 0) is shown in [PITH_FULL_IMAGE:figures/full_fig_p021_4_1.png]
Figure 4.2
Figure 4.2. Figure 4.2: Example 2. Barenblatt solution. Computational meshes. 21 [PITH_FULL_IMAGE:figures/full_fig_p021_4_2.png]
Figure 4.3
Figure 4.3. Figure 4.3: Example 2: Barenblatt solution in 2D. Numerical approximations and errors at T = 2 to the Barenblatt solution of the two-dimensional PME. 22 [PITH_FULL_IMAGE:figures/full_fig_p022_4_3.png]
Figure 4.4
Figure 4.4. Figure 4.4: Example 2: Barenblatt solution in 3D. Numerical approximation and error at T = 3 to the Barenblatt solution of the three-dimensional PME with m = 3. Example 3. PME on a torus In this example, we solve the porous medium equation (1.1) on a three-dimensional, doughnut￾…
Figure 4.5
Figure 4.5. Figure 4.5: Example 3: PME on a torus. Computational mesh. The initial condition on the torus is given by u0(x, y, z) = 2.5 − r(x, y, z)  B [PITH_FULL_IMAGE:figures/full_fig_p023_4_5.png]
Figure 4.6
Figure 4.6. Figure 4.6: Example 3. Numerical solutions of the PME on a torus at different times. 5 Conclusions In this paper, we develop a class of efficient exponential time differencing Runge-Kutta discontin￾uous Galerkin methods for solving nonlinear degenerate parabolic equations. The p…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

77 extracted references · 76 canonical work pages

  1. [72]

    Z. Xu, Z. Sun and Y.-T. Zhang, Stability and Time-Step Constraints of Exponential Time Differencing Runge–Kutta Discontinuous Galerkin Methods f or Advection-Diffusion Equations , (2025), arXiv preprint arXiv:2503.03019

    arXiv 2025

  2. [73]

    Xu and Y.-T

    Z. Xu and Y.-T. Zhang, High-order exponential time differencing multi-resolutio n alterna- tive finite difference WENO methods for nonlinear degenerate parabolic equations, J. Comput. Phys., v529 (2025), p.113838

    work page 2025

  3. [1]

    Abedian, H

    R. Abedian, H. Adibi and M. Dehghan, A high-order weighted essentially non-oscillatory (WENO) finite difference scheme for nonlinear degenerate par abolic equations, Comput. Phys. 29 Commun., 184 (2013), pp.1874-1888

    work page 2013

  4. [2]

    Abedian and M

    R. Abedian and M. Dehghan, A RBF-WENO finite difference scheme for non-linear degenerat e parabolic equations, J. Sci. Comput., 93 (2022), 60

    work page 2022

  5. [3]

    Ahmat and J

    M. Ahmat and J. Qiu, Hybrid HWENO method for nonlinear degenerate parabolic equ ations, J. Sci. Comput., 96 (2023), 83

    work page 2023

  6. [4]

    Arbogast, C.S

    T. Arbogast, C.S. Huang and X. Zhao, Finite volume WENO schemes for nonlinear parabolic problems with degenerate diffusion on non-uniform meshes , J. Comput. Phys., 399 (2019), p.108921

    work page 2019

  7. [5]

    D. N. Arnold, An interior penalty finite element method with discontinuou s elements , SIAM J. Numer. Anal., 19 (1982), pp. 742–760

    work page 1982

  8. [6]

    Aregba-Driollet, R

    D. Aregba-Driollet, R. Natalini and S. Tang, Explicit diffusive kinetic schemes for nonlinear degenerate parabolic systems, Math. Comp., 73 (2004), pp. 63-94

    work page 2004

Show all 77 references
  1. [7]

    D. G. Aronson, The porous medium equation , in Nonlinear Diffusion Problems, Lecture Notes in Math. 1224, Springer Berlin, 1986, pp. 1–46

    1986

  2. [8]

    Ascher, S

    U. Ascher, S. Ruuth and B. Wetton, Implicit-explicit methods for time-dependent PDE’s, SIAM J. Numer. Anal., 32 (1995), pp. 797-823

    1995

  3. [9]

    Baines, M.E

    M.J. Baines, M.E. Hubbard and P.K. Jimack, A moving mesh finite element algorithm for the adaptive solution of time-dependent partial differenti al equations with moving boundaries , Appl. Numer. Math., 54 (2005), pp.450-469

    2005

  4. [10]

    Barenblatt, On some unsteady motions of a liquid or a gas on a porous medium , Prikl

    G.I. Barenblatt, On some unsteady motions of a liquid or a gas on a porous medium , Prikl. Mat. Meh., 16 (1952), pp. 67-78

    1952

  5. [11]

    Bessemoulin-Chatard and F

    M. Bessemoulin-Chatard and F. Filbet, A finite volume scheme for nonlinear degenerate parabolic equations, SIAM J. Sci. Comput. 34 (2012), pp. B559-B583

    2012

  6. [12]

    Beylkin, J.M

    G. Beylkin, J.M. Keiser and L. Vozovoi, A new class of time discretization schemes for the solution of nonlinear PDEs , J. Comput. Phys., 147 (1998), pp. 362-387

    1998

  7. [13]

    Boscarino, R

    S. Boscarino, R. B¨ urger, P. Mulet, G. Russo and L.M. Vil lada, Linearly implicit IMEX Runge– Kutta methods for a class of degenerate convection-diffusio n problems, SIAM J. Sci. Comput., 37 (2015), pp.B305-B331

    2015

  8. [14]

    B¨ urger, P

    R. B¨ urger, P. Mulet and L.M. Villada, Regularized nonlinear solvers for IMEX methods applied to diffusively corrected multispecies kinematic flow models , SIAM J. Sci. Comput., 35 (2013), pp.B751-B777

    2013

  9. [15]

    Cavalli, G

    F. Cavalli, G. Naldi, G. Puppo and M. Semplice, High-order relaxation schemes for nonlinear degenerate diffusion problems , SIAM J. Numer. Anal., 45 (2007), pp. 2098-2119

    2007

  10. [16]

    Chen and C.-W

    T. Chen and C.-W. Shu, Entropy stable high order discontinuous Galerkin methods w ith suitable quadrature rules for hyperbolic conservation laws , J. Comput. Phys., 345 (2017), pp.427-461

    2017

  11. [17]

    Chen and Y.-T

    S. Chen and Y.-T. Zhang, Krylov implicit integration factor methods for spatial dis cretization on high dimensional unstructured meshes: application to di scontinuous Galerkin methods , J. Comput. Phys., 230 (2011), pp. 4336-4352. 30

    2011

  12. [18]

    Cheng and C.-W

    Y. Cheng and C.-W. Shu, A discontinuous Galerkin finite element method for time depe ndent partial differential equations with higher order derivativ es, Math. Comput., 77 (2008), pp.699- 730

    2008

  13. [19]

    explicit

    A. Christlieb, W. Guo, Y. Jiang and H. Yang, Kernel based high order “explicit” uncondition- ally stable scheme for nonlinear degenerate advection-diff usion equations , J. Sci. Comput., 82 (2020), 52

    2020

  14. [20]

    Cockburn, S

    B. Cockburn, S. Hou, and C.-W. Shu, TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws IV: The mu ltidimensional case , Math. Comput., 54 (1990), pp. 545–581

    1990

  15. [21]

    Cockburn, S.-Y

    B. Cockburn, S.-Y. Lin, and C.-W. Shu, TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: One di mensional systems, J. Comput. Phys., 84 (1989), pp. 90–113

    1989

  16. [22]

    Cockburn and C.-W

    B. Cockburn and C.-W. Shu, TVB Runge-Kutta local projection discontinuous Galerkin fin ite element method for scalar conservation laws II: General frame work, Math. Comput., 52 (1989), pp. 411–435

    1989

  17. [23]

    Cockburn and C.-W

    B. Cockburn and C.-W. Shu, The Runge-Kutta local projection P1-discontinuous-Galerk in finite element method for scalar conservation laws , ESAIM Math. Model. Numer. Anal., 25 (1991), pp. 337–361

    1991

  18. [24]

    Cockburn and C.-W

    B. Cockburn and C.-W. Shu, The Runge-Kutta discontinuous Galerkin finite element metho d for conservation laws V: Multidimensional systems , J. Comput. Phys., 141 (1998), pp. 199–224

    1998

  19. [25]

    Cockburn and C.-W

    B. Cockburn and C.-W. Shu, The local discontinuous Galerkin method for time-dependent convection-diffusion systems , SIAM J. Numer. Anal., 35 (1998), pp.2440-2463

    1998

  20. [26]

    Cox and P.C

    S.M. Cox and P.C. Matthews, Exponential time differencing for stiff systems , J. Comput. Phys., 176 (2002), pp.430-455

    2002

  21. [27]

    Dobrev, R.D

    V.A. Dobrev, R.D. Lazarov, P.S. Vassilevski and L.T. Zi katanov, Two-level preconditioning of discontinuous Galerkin approximations of second-order elliptic equations , Numer. Linear Algebra Appl., 13 (2006), pp.753-770

    2006

  22. [28]

    Q. Du, L. Ju, X. Li and Z. Qiao, Maximum bound principles for a class of semilinear paraboli c equations and exponential time-differencing schemes , SIAM Rev., 63 (2021), pp. 317-359

    2021

  23. [29]

    Gaudreault and J.A

    S. Gaudreault and J.A. Pudykiewicz, An efficient exponential time integration method for the numerical solution of the shallow water equations on the sph ere, J. Comput. Phys., 322 (2016), pp.827-848

    2016

  24. [30]

    Gottlieb, D.I

    S. Gottlieb, D.I. Ketcheson and C.-W. Shu, High order strong stability preserving time dis- cretizations, J. Sci. Comput., 38 (2009), pp.251-289

    2009

  25. [31]

    Gottlieb, D.I

    S. Gottlieb, D.I. Ketcheson and C.-W. Shu, Strong stability preserving Runge-Kutta and mul- tistep time discretizations , World Scientific, 2011

    2011

  26. [32]

    Gottlieb, C.-W

    S. Gottlieb, C.-W. Shu and E. Tadmor, Strong stability-preserving high-order time discretiza- tion methods, SIAM Rev., 43 (2001), pp.89–112. 31

    2001

  27. [33]

    Hesthaven and T

    J.S. Hesthaven and T. Warburton, Nodal discontinuous Galerkin methods: algorithms, analy- sis, and applications , Springer Science & Business Media, (2007)

    2007

  28. [34]

    Hochbruck and A

    M. Hochbruck and A. Ostermann, Exponential integrators, Acta Numer., 19 (2010), pp. 209- 286

    2010

  29. [35]

    Hochbruck, A

    M. Hochbruck, A. Ostermann and J. Schweitzer, Exponential Rosenbrock-type methods, SIAM J. Numer. Anal., 47 (2009), pp.786-803

    2009

  30. [36]

    Jiang, High order finite difference multi-resolution WENO method fo r nonlinear degenerate parabolic equations, J

    Y. Jiang, High order finite difference multi-resolution WENO method fo r nonlinear degenerate parabolic equations, J. Sci. Comput., 86 (2021), pp.1-20

    2021

  31. [37]

    Jiang and Y.-T

    T. Jiang and Y.-T. Zhang, Krylov single-step implicit integration factor WENO metho ds for advection–diffusion–reaction equations, J. Comput. Phys., 311 (2016), pp.22-44

    2016

  32. [38]

    Kassam and L.N

    A.K. Kassam and L.N. Trefethen, Fourth-order time-stepping for stiff PDEs , SIAM J. Sci. Comput., 26 (2005), pp.1214-1233

    2005

  33. [39]

    C. A. Kennedy and M. H. Carpenter, Additive Runge–Kutta schemes for convection-diffusion- reaction equations, Appl. Numer. Math., 44 (2003), pp. 139-181

    2003

  34. [40]

    Ketcheson, C.B

    D.I. Ketcheson, C.B. Macdonald and S. Gottlieb, Optimal implicit strong stability preserving Runge-Kutta methods, Appl. Numer. Math., 59 (2009), pp. 373–392

    2009

  35. [41]

    Kraus and S.K

    J.K. Kraus and S.K. Tomar, Multilevel preconditioning of two-dimensional elliptic p roblems discretized by a class of discontinuous Galerkin methods , SIAM J. Sci. Comput., 30 (2008), pp.684-706

    2008

  36. [42]

    Kurganov and E

    A. Kurganov and E. Tadmor, New high-resolution central schemes for nonlinear conserv ation laws and convection-diffusion equations , J. Comput. Phys., 160 (2000), pp. 241-282

    2000

  37. [43]

    Li and Z

    J. Li and Z. Du, A two-stage fourth order time-accurate discretization for Lax–Wendroff type flow solvers I. Hyperbolic conservation laws , SIAM J. Sci. Comput., 38 (2016), pp.A3046- A3069

    2016

  38. [44]

    Liu and X

    S. Liu and X. Liu, Exponential Time Differencing Method for a Reaction-Diffusio n System with Free Boundary , Commun. Appl. Math. Comput., 6 (2024), pp.354-371

    2024

  39. [45]

    Y. Liu, Y. Cheng, S. Chen and Y.-T. Zhang, Krylov implicit integration factor discontinuous Galerkin methods on sparse grids for high dimensional react ion-diffusion equations , J. Comput. Phys., 388 (2019), pp.90-102

    2019

  40. [46]

    Liu, C.-W

    Y. Liu, C.-W. Shu and M. Zhang, High order finite difference WENO schemes for nonlinear degenerate parabolic equations, SIAM J. Sci. Comput., 33 (2011), pp.939-965

    2011

  41. [47]

    Liu and J

    H. Liu and J. Yan, The direct discontinuous Galerkin (DDG) methods for diffusio n problems, SIAM J. Numer. Anal., 47 (2009), pp. 675–698

    2009

  42. [48]

    Liu and J

    H. Liu and J. Yan, The direct discontinuous Galerkin (DDG) method for diffusion with interface corrections, Commun. Comput. Phys., 8 (2010), p.541

    2010

  43. [49]

    J. Lu, Y. Liu and C.-W. Shu, An oscillation-free discontinuous Galerkin method for sca lar hyperbolic conservation laws , SIAM J. Numer. Anal., 59 (2021), pp.1299-1324. 32

    2021

  44. [50]

    Lu and Y.-T

    D. Lu and Y.-T. Zhang, Krylov integration factor method on sparse grids for high sp atial dimension convection-diffusion equations , J. Sci. Comput., 69 (2016), pp. 736-763

    2016

  45. [51]

    Luan, J.A

    V.T. Luan, J.A. Pudykiewicz and D.R. Reynolds, Further development of efficient and accurate time integration schemes for meteorological models , J. Comput. Phys., 376 (2019), pp.817-837

    2019

  46. [52]

    Magenes, R

    E. Magenes, R. H. Nochetto, and C. Verdi, Energy error estimates for a linear scheme to approximate nonlinear parabolic problems, ESAIM Math. Model. Numer. Anal., 21 (1987), pp. 655-678

    1987

  47. [53]

    Moler and C

    C. Moler and C. Van Loan, Nineteen dubious ways to compute the exponential of a matrix , twenty-five years later , SIAM Rev., 45 (2003), pp.3-49

    2003

  48. [54]

    Nie, Y.-T

    Q. Nie, Y.-T. Zhang and R. Zhao, Efficient semi-implicit schemes for stiff systems , J. Comput. Phys., 214 (2006), pp.521-537

    2006

  49. [55]

    Ngo and W

    C. Ngo and W. Huang, A study on moving mesh finite element solution of the porous me dium equation, J. Comput. Phys. 331 (2017), 357-380

    2017

  50. [56]

    Niesen and W.M

    J. Niesen and W.M. Wright, Algorithm 919: A Krylov subspace algorithm for evaluating t he ϕ-functions appearing in exponential integrators , ACM Trans. Math. Software (TOMS), 38 (2012), pp.1-19

    2012

  51. [57]

    R. H. Nochetto, A. Schmidt, and C. Verdi, A posteriori error estimation and adaptivity for degenerate parabolic problems, Math. Comp., 69 (2000), pp. 1-24

    2000

  52. [58]

    Peraire and P.O

    J. Peraire and P.O. Persson, The compact discontinuous Galerkin (CDG) method for elliptic problems, SIAM J. Sci. Comput., 30 (2008), pp.1806-1824

    2008

  53. [59]

    M. Peng, Z. Sun and K. Wu, OEDG: Oscillation-eliminating discontinuous Galerkin me thod for hyperbolic conservation laws , Math. Comput., 94 (2025), pp.1147-1198

    2025

  54. [60]

    Qin and C.-W

    T. Qin and C.-W. Shu, Implicit positivity-preserving high-order discontinuous Galerkin methods for conservation laws , SIAM J. Sci. Comput., 40 (2018), pp.A81-A107

    2018

  55. [61]

    J. Qiu, M. Dumbser and C.-W. Shu The discontinuous Galerkin method with Lax–Wendroff type time discretizations , Comput. Methods Appl. Mech. Eng., 194 (2005), pp.4528-454 3

    2005

  56. [62]

    W. H. Reed and T. R. Hill, Triangular mesh methods for the neutron transport equation (No. LA-UR-73-479; CONF-730414-2), Los Alamos Scientific Lab., New Mexico, (1973)

    1973

  57. [63]

    Shu, Total-Variation-Diminishing time discretizations, SIAM J

    C.-W. Shu, Total-Variation-Diminishing time discretizations, SIAM J. Sci. and Stat. Comput., 9 (1988), pp.1073-1084

    1988

  58. [64]

    Q. Tao, Y. Liu, Y. Jiang, and J. Lu, An oscillation free local discontinuous Galerkin method for nonlinear degenerate parabolic equations , Numer. Methods Partial Differ. Equ., 39 (2023), pp. 3145-3169

    2023

  59. [65]

    Tokman, Efficient integration of large stiff systems of ODEs with expon ential propagation iterative (EPI) methods , J

    M. Tokman, Efficient integration of large stiff systems of ODEs with expon ential propagation iterative (EPI) methods , J. Comput. Phys., 213 (2006), pp.748-776

    2006

  60. [66]

    Trefethen and D

    L.N. Trefethen and D. Bau III, Numerical Linear Algebra , SIAM, 1997. 33

    1997

  61. [67]

    Vijaywargiya and G

    A. Vijaywargiya and G. Fu, Two finite element approaches for the porous medium equation that are positivity preserving and energy stable , J. Sci. Comput., 100 (2024), p.86

    2024

  62. [68]

    M. F. Wheeler, An elliptic collocation-finite element method with interio r penalties , SIAM J. Numer. Anal., 15 (1978), pp. 152–161

    1978

  63. [69]

    Xu and C.-W

    Z. Xu and C.-W. Shu, Third order maximum-principle-satisfying and positivity- preserving Lax- Wendroff discontinuous Galerkin methods for hyperbolic con servation laws , J. Comput. Phys., 470 (2022), p.111591

    2022

  64. [70]

    Xu and C.-W

    Z. Xu and C.-W. Shu, On the conservation property of positivity-preserving dis continuous Galerkin methods for stationary hyperbolic equations , J. Comput. Phys., 490 (2023), 112304

    2023

  65. [71]

    Xu and C.-W

    Z. Xu and C.-W. Shu, A high-order well-balanced discontinuous Galerkin method for hyperbolic balance laws based on the Gauss-Lobatto quadrature rules , J. Sci. Comput., 101 (2024), 39

    2024

  66. [74]

    Zhang and T

    P. Zhang and T. Xiong, High order implicit finite difference schemes with a semi-imp licit WENO reconstruction for nonlinear degenerate parabolic eq uations, J. Comput. Phys., 467 (2022), p.111442

    2022

  67. [75]

    Zhang and Z.-L

    Q. Zhang and Z.-L. Wu, Numerical simulation for porous medium equation by local di scontin- uous Galerkin finite element method , J. Sci. Comput., 38 (2009), 127-148

    2009

  68. [76]

    Zhang, X

    Y. Zhang, X. Zhang, and C.-W. Shu, Maximum-principle-satisfying second order discontinuous Galerkin schemes for convection-diffusion equations on tri angular meshes , J. Comput. Phys., 234 (2013), pp. 295–316

    2013

  69. [77]

    Zhang and C.-W

    M. Zhang and C.-W. Shu, An analysis of three different formulations of the discontin uous Galerkin method for diffusion equations , Math. Models Methods Appl. Sci., 13 (2003), pp. 395–413. 34

    2003

Pith tools

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