pith. sign in

arxiv: 2604.26336 · v1 · submitted 2026-04-29 · 🧮 math.NA · cs.NA

A Second-Order Maximum-Principle-Preserving Crouzeix-Raviart Finite Element Method for Time-dependent Transport Equation

Shipeng Mao , Mingyang Zhang This is my paper

Pith reviewed 2026-05-07 12:56 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Crouzeix-Raviart elementmaximum principletransport equationfinite element methodflux-corrected transportexplicit time steppingviscosity method
0
0 comments X

The pith

A Crouzeix-Raviart finite element method for the time-dependent transport equation is made explicit, second-order accurate, and maximum-principle-preserving through its diagonal mass matrix and controlled viscosity additions.

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

The paper constructs a finite element discretization for the two-dimensional time-dependent transport equation that remains explicit in time, reaches second-order accuracy, and never produces values outside the initial range of the data. The central device is the Crouzeix-Raviart element whose mass matrix is diagonal on triangular meshes; this property removes the need to solve linear systems at each step and lets the authors add simple viscosity terms that enforce the maximum principle. Low-order schemes are built first with minimum or bilinear viscosities, then lifted to second order by greedy or flux-corrected-transport corrections, with a special limiter added at inflow boundaries. A Wachspress-coordinate reconstruction produces a continuous piecewise-linear field on a refined mesh that obeys the bound everywhere. When the velocity is divergence-free the schemes are also conservative, and tests on both smooth and discontinuous data confirm the claimed properties.

Core claim

We construct an explicit, second-order, and maximum-principle-preserving Crouzeix-Raviart finite element method for the two-dimensional time-dependent transport equation. The key observation is that the mass matrix of the CR element has diagonal structure, which allows us to avoid solving a large linear system at each time step and helps to construct a low-order scheme that preserves the maximum principle in a simple way. We first introduce low-order schemes based on minimum and bilinear viscosities, then recover second-order accuracy by means of greedy and flux-corrected transport viscosities. For inflow boundary conditions we further design a modified FCT limiter. In addition we propose a

What carries the argument

The diagonal mass matrix of the Crouzeix-Raviart element on triangular meshes, which permits fully explicit time stepping and direct construction of minimum or bilinear viscosity terms that enforce the maximum principle.

If this is right

  • The method steps forward in time without solving linear systems at each step.
  • Solutions remain bounded by the initial data range even when discontinuities are present.
  • The scheme is conservative whenever the velocity field is divergence-free.
  • A continuous piecewise-linear reconstruction that also obeys the maximum principle is obtained on the refined mesh.
  • Both solenoidal and non-solenoidal velocities are handled without loss of the stated properties.

Where Pith is reading between the lines

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

  • The same diagonal-mass-matrix property may allow similar explicit maximum-principle schemes for other linear hyperbolic problems on triangular meshes.
  • The reconstruction step could be replaced by other averaging techniques if the goal is only nodal bounds rather than a globally continuous field.
  • Extension to three space dimensions would require checking whether an analogous diagonal mass matrix exists for the three-dimensional Crouzeix-Raviart element.

Load-bearing premise

The mass matrix of the Crouzeix-Raviart element on the triangular mesh used is exactly diagonal.

What would settle it

A run on a smooth test problem with known exact solution in which the computed field at any node exceeds the maximum or minimum value attained by the exact initial data, or in which the observed convergence rate falls below two.

Figures

Figures reproduced from arXiv: 2604.26336 by Mingyang Zhang, Shipeng Mao.

Figure 1
Figure 1. Figure 1: The CR interpolation attains a value of −3 at the lower-left vertex of the blue triangle. 𝐲6 𝐲1 𝐲2 𝐲3 𝐲4 𝐲5 𝐲0 view at source ↗
Figure 3
Figure 3. Figure 3: Solid body rotation at 𝑡 = 1: CR method with global FCT viscosity. Mesh sizes ℎ = 0.05, 0.025, 0.0125 (left to right). Reconstruction via Wachspress coordinates. S. Mao, M. Zhang: Preprint submitted to Elsevier Page 20 of 19 view at source ↗
Figure 4
Figure 4. Figure 4: Solid body rotation at 𝑡 = 1: continuous 𝑃1 (left) vs. CR method (right), both with global FCT viscosity, ℎ = 0.0125 view at source ↗
Figure 5
Figure 5. Figure 5: Numerical solutions for the non-solenoidal compressive rotational test at 𝑡 = 2. The left and right panels show the CR and continuous 𝑃1 solutions, respectively, both with global FCT viscosity on the mesh with 𝑁 = 80. The CR method produces a more robust solution with less dispersion compared to the continuous 𝑃1 method. S. Mao, M. Zhang: Preprint submitted to Elsevier Page 21 of 19 view at source ↗
read the original abstract

In this paper, we construct an explicit, second-order, and maximum-principle-preserving Crouzeix-Raviart (CR) finite element method for two-dimensional time-dependent transport equation. The key observation is that the mass matrix of the CR element is with diagonal structure, which allows us to avoid the need to solve a large linear system for each time step and help to construct a low-order scheme that preserves the maximum principle in a simple way. We first introduce low-order schemes based on minimum and bilinear viscosities, and then recover second-order accuracy by means of greedy and flux-corrected transport viscosities. For inflow boundary conditions, we further design a modified FCT limiter. In addition, we propose a simple reconstruction based on Wachspress coordinates to obtain a continuous piecewise linear approximation on the $\frac{h}2$-mesh that satisfies the maximum principle on the whole domain. Under divergence-free velocity fields, the proposed schemes are conservative. Numerical experiments on both smooth and discontinuous test cases, with both solenoidal and non-solenoidal velocity fields, confirm the accuracy and robustness of the proposed schemes.

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

3 major / 2 minor

Summary. The paper constructs an explicit, second-order, maximum-principle-preserving Crouzeix-Raviart finite element method for the two-dimensional time-dependent transport equation. It relies on the claimed diagonal structure of the CR mass matrix to enable explicit forward-Euler stepping without solving linear systems at each step. Low-order schemes are built using minimum and bilinear artificial viscosities that preserve the maximum principle; second-order accuracy is recovered via greedy and flux-corrected transport (FCT) viscosities, with a modified FCT limiter for inflow boundaries. A Wachspress-coordinate reconstruction produces a continuous piecewise-linear approximation satisfying the maximum principle on the refined mesh. The schemes are conservative for divergence-free velocities. Numerical experiments on smooth and discontinuous solutions with solenoidal and non-solenoidal fields are reported to confirm accuracy and robustness.

Significance. If the diagonal-mass-matrix property and the accompanying proofs hold, the work supplies a practical explicit high-order MP-preserving discretization for hyperbolic transport that avoids implicit solves while retaining second-order accuracy. The combination of CR elements with FCT and the boundary reconstruction is a reasonable extension of existing viscosity and limiter techniques, and the conservation property under divergence-free advection is a useful feature for long-time simulations.

major comments (3)
  1. [Introduction] Introduction (key observation paragraph): The assertion that 'the mass matrix of the CR element is with diagonal structure' is load-bearing for the explicit stepping and the simple viscosity construction, yet no proof, mesh restriction, or quadrature rule is supplied. On a general triangulation the local mass matrix for standard CR basis functions (edge-midpoint degrees of freedom) has nonzero off-diagonal entries; the global matrix is therefore not diagonal. This directly affects whether the low-order scheme remains explicit and whether the discrete maximum principle proof carries over without additional lumping arguments.
  2. [Low-order schemes] Section on low-order schemes (viscosity definitions): The minimum and bilinear viscosity constructions are presented as preserving the maximum principle, but the proofs appear to invoke the diagonal mass matrix at the algebraic level. If the mass matrix is not diagonal, the algebraic maximum-principle argument must be re-established after any lumping or quadrature is introduced; the current derivation does not address this case.
  3. [Numerical experiments] Numerical experiments section: While the abstract states that experiments confirm accuracy and robustness, the manuscript provides no tabulated L^∞ or L^2 errors versus mesh size, no explicit values for the viscosity parameters or limiter thresholds, and no mesh-quality metrics. Without these data it is impossible to verify the claimed second-order convergence or the strict satisfaction of the discrete maximum principle on the reported test cases.
minor comments (2)
  1. [Abstract] The phrase 'is with diagonal structure' is grammatically awkward; 'has a diagonal structure' or 'is diagonal' would be clearer.
  2. [Reconstruction] The reconstruction step using Wachspress coordinates is described only briefly; a short algorithmic outline or pseudocode would improve reproducibility.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments provided. We address each of the major comments in detail below and outline the revisions we intend to make.

read point-by-point responses

  1. Referee: [Introduction] Introduction (key observation paragraph): The assertion that 'the mass matrix of the CR element is with diagonal structure' is load-bearing for the explicit stepping and the simple viscosity construction, yet no proof, mesh restriction, or quadrature rule is supplied. On a general triangulation the local mass matrix for standard CR basis functions (edge-midpoint degrees of freedom) has nonzero off-diagonal entries; the global matrix is therefore not diagonal. This directly affects whether the low-order scheme remains explicit and whether the discrete maximum principle proof carries over without additional lumping arguments.

    Authors: We appreciate the referee's identification of this critical point. The manuscript's claim refers to the use of a lumped mass matrix via a suitable quadrature rule (specifically, the midpoint quadrature on edges for CR elements), which renders the mass matrix diagonal. This is a standard technique for obtaining explicit schemes with nonconforming elements. However, we acknowledge that the current version lacks a detailed explanation and proof of this property. In the revised manuscript, we will include a new subsection detailing the quadrature rule, proving the diagonal structure, and showing how it preserves the algebraic conditions for the discrete maximum principle without additional arguments beyond lumping. revision: yes

  2. Referee: [Low-order schemes] Section on low-order schemes (viscosity definitions): The minimum and bilinear viscosity constructions are presented as preserving the maximum principle, but the proofs appear to invoke the diagonal mass matrix at the algebraic level. If the mass matrix is not diagonal, the algebraic maximum-principle argument must be re-established after any lumping or quadrature is introduced; the current derivation does not address this case.

    Authors: We agree that the proofs should explicitly account for the lumping procedure. The minimum and bilinear viscosities are designed such that, when combined with the lumped diagonal mass matrix, the scheme satisfies the discrete maximum principle through positive coefficients and row-sum unity. We will revise the relevant section to first introduce the lumping, then derive the viscosity terms, and provide a complete proof that incorporates these steps, ensuring the maximum principle holds for the explicit scheme. revision: yes

  3. Referee: [Numerical experiments] Numerical experiments section: While the abstract states that experiments confirm accuracy and robustness, the manuscript provides no tabulated L^∞ or L^2 errors versus mesh size, no explicit values for the viscosity parameters or limiter thresholds, and no mesh-quality metrics. Without these data it is impossible to verify the claimed second-order convergence or the strict satisfaction of the discrete maximum principle on the reported test cases.

    Authors: We concur that additional quantitative information would strengthen the numerical section. In the revised version, we will add tables with L^∞ and L^2 error norms for successive mesh refinements to demonstrate the second-order accuracy, list the specific viscosity coefficients and limiter parameters employed in each test case, and include mesh quality measures such as the minimum angle in the triangulations. This will facilitate verification of the convergence rates and the maximum principle preservation. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation uses standard element property and FCT techniques without self-referential reduction

full rationale

The paper's central construction begins from the asserted diagonal mass-matrix property of the CR element (used to enable explicit stepping and simple viscosity), then applies standard minimum/bilinear viscosity for a low-order scheme, followed by greedy/FCT correction for second-order accuracy and a Wachspress-based reconstruction for continuity. None of these steps define the target properties (second-order accuracy or DMP) in terms of themselves, nor rename fitted quantities as predictions, nor rely on self-citation chains for uniqueness. The mass-matrix claim is presented as an observation rather than derived from the method's outputs, so no self-definitional loop exists. The overall chain remains independent of the final claims.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Reviewed from abstract only; full derivation, parameter choices, and mesh assumptions are not visible.

axioms (1)
  • domain assumption The mass matrix of the Crouzeix-Raviart element is diagonal
    Stated as the key observation that enables explicit time stepping and simple viscosity construction.

pith-pipeline@v0.9.0 · 5499 in / 1304 out tokens · 64565 ms · 2026-05-07T12:56:00.437784+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

41 extracted references · 41 canonical work pages

  1. [1]

    The maximum angle condition for mixed and nonconforming elements: application to the Stokes equations

    Acosta, G., Durán, R.G., 1999. The maximum angle condition for mixed and nonconforming elements: application to the Stokes equations. SIAM J. Numer. Anal. 37, 18–36

    work page 1999

  2. [2]

    A new upwind weak galerkin finite element method for linear hyperbolic equations

    AL-Taweel, A., Mu, L., 2021. A new upwind weak galerkin finite element method for linear hyperbolic equations. J. Comput. Appl. Math. 390, Article ID 113376, 12 p

    work page 2021

  3. [3]

    On nonconforming linear-constant elements for some variants of the Stokes equations

    Arnold, D.N., 1993. On nonconforming linear-constant elements for some variants of the Stokes equations. Istit. Lombardo Accad. Sci. Lett. Rend. A 127, 83–93

    work page 1993

  4. [4]

    On locking and robustness in the finite element method

    Babuška, I., Suri, M., 1992. On locking and robustness in the finite element method. SIAM J. Numer. Anal. 29, 1261–1293

    work page 1992

  5. [5]

    Analysis of a combined CG1-DG2 method for the transport equation

    Becker, R., Bittl, M., Kuzmin, D., 2015. Analysis of a combined CG1-DG2 method for the transport equation. SIAM J. Numer. Anal. 53, 445–463

    work page 2015

  6. [6]

    Quasi-optimality of adaptive nonconforming finite element methods for the Stokes equations

    Becker, R., Mao, S., 2011. Quasi-optimality of adaptive nonconforming finite element methods for the Stokes equations. SIAM J. Numer. Anal. 49, 970–991

    work page 2011

  7. [7]

    Flux-corrected transport

    Boris, J.P., Book, D.L., 1997. Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works [J. Comput. Phys.11(1973), no. 1, 38–69], volume 135, pp. 170–186. With an introduction by Steven T. Zalesak, Commemoration of the 30th anniversary {of J. Comput. Phys.}

    work page 1997

  8. [8]

    Anonconformingmixedmultigridmethodforthepuredisplacementprobleminplanarlinearelasticity

    Brenner,S.C.,1993. Anonconformingmixedmultigridmethodforthepuredisplacementprobleminplanarlinearelasticity. SIAMJ.Numer. Anal. 30, 116–135

    work page 1993

  9. [9]

    Forty years of the Crouzeix-Raviart element

    Brenner, S.C., 2015. Forty years of the Crouzeix-Raviart element. Numer. Methods Partial Differential Equations 31, 367–396. S. Mao, M. Zhang:Preprint submitted to ElsevierPage 18 of 19 CR FEM for Time-dependent transport equations

    work page 2015

  10. [10]

    Linear finite element methods for planar linear elasticity

    Brenner, S.C., Sung, L.Y., 1992. Linear finite element methods for planar linear elasticity. Math. Comp. 59, 321–338

    work page 1992

  11. [11]

    Brooks,A.N.,Hughes,T.J.R.,1982.Streamlineupwind/Petrov-Galerkinformulationsforconvectiondominatedflowswithparticularemphasis ontheincompressibleNavier-Stokesequations.Comput.MethodsAppl.Mech.Engrg.32,199–259.FENOMECH”81,PartI(Stuttgart,1981)

    work page 1982

  12. [12]

    Optimal adaptive nonconforming FEM for the Stokes problem

    Carstensen, C., Peterseim, D., Rabus, H., 2013. Optimal adaptive nonconforming FEM for the Stokes problem. Numer. Math. 123, 291–308

    work page 2013

  13. [13]

    The Runge-Kutta discontinuous Galerkin method for conservation laws

    Cockburn, B., Shu, C.W., 1998. The Runge-Kutta discontinuous Galerkin method for conservation laws. V. Multidimensional systems. J. Comput. Phys. 141, 199–224

    work page 1998

  14. [14]

    Methods Appl

    Codina,R.,1993.Adiscontinuity-capturingcrosswind-dissipationforthefiniteelementsolutionoftheconvection-diffusionequation.Comput. Methods Appl. Mech. Engrg. 110, 325–342

    work page 1993

  15. [15]

    ConformingandnonconformingfiniteelementmethodsforsolvingthestationaryStokesequations.I

    Crouzeix,M.,Raviart,P.A.,1973. ConformingandnonconformingfiniteelementmethodsforsolvingthestationaryStokesequations.I. Rev. Française Automat. Informat. Recherche Opérationnelle Sér. Rouge 7, 33–75

    work page 1973

  16. [16]

    Finite elements I—Approximation and interpolation

    Ern, A., Guermond, J.L., [2021]©2021. Finite elements I—Approximation and interpolation. volume 72 ofTexts in Applied Mathematics. Springer, Cham

    work page 2021

  17. [17]

    Convergence of a numerical method for the compressible Navier-Stokes system on general domains

    Feireisl, E., Karper, T., Michálek, M., 2016. Convergence of a numerical method for the compressible Navier-Stokes system on general domains. Numer. Math. 134, 667–704

    work page 2016

  18. [18]

    Wachspressandmeanvaluecoordinates,in:ApproximationtheoryXIV:SanAntonio2013.Springer,Cham.volume83 ofSpringer Proc

    Floater,M.S.,2014. Wachspressandmeanvaluecoordinates,in:ApproximationtheoryXIV:SanAntonio2013.Springer,Cham.volume83 ofSpringer Proc. Math. Stat., pp. 81–102

    work page 2014

  19. [19]

    A convergent finite element-finite volume scheme for the compressible Stokes problem

    Gallouët, T., Herbin, R., Latché, J.C., 2009. A convergent finite element-finite volume scheme for the compressible Stokes problem. I. The isothermal case. Math. Comp. 78, 1333–1352

    work page 2009

  20. [20]

    A difference method for numerical calculation of discontinuous solutions of the equations of hydrodynamics

    Godunov, S.K., 1959. A difference method for numerical calculation of discontinuous solutions of the equations of hydrodynamics. Mat. Sb. (N.S.) 47(89), 271–306

    work page 1959

  21. [21]

    Strong stability-preserving high-order time discretization methods

    Gottlieb, S., Shu, C.W., Tadmor, E., 2001. Strong stability-preserving high-order time discretization methods. SIAM Rev. 43, 89–112

    work page 2001

  22. [22]

    A maximum-principle preserving𝐶0 finite element method for scalar conservation equations

    Guermond, J.L., Nazarov, M., 2014. A maximum-principle preserving𝐶0 finite element method for scalar conservation equations. Comput. Methods Appl. Mech. Engrg. 272, 198–213

    work page 2014

  23. [23]

    A second-order maximum principle preserving Lagrange finite element technique for nonlinear scalar conservation equations

    Guermond, J.L., Nazarov, M., Popov, B., Yang, Y., 2014. A second-order maximum principle preserving Lagrange finite element technique for nonlinear scalar conservation equations. SIAM J. Numer. Anal. 52, 2163–2182

    work page 2014

  24. [24]

    Invariant domains and second-order continuous finite element approximation for scalar conservation equations

    Guermond, J.L., Popov, B., 2017. Invariant domains and second-order continuous finite element approximation for scalar conservation equations. SIAM J. Numer. Anal. 55, 3120–3146

    work page 2017

  25. [25]

    Matrix-free subcell residual distribution for bernstein finite element discretizations of linear advection equations

    Hajduk, H., Kuzmin, D., Kolev, T., Abgrall, R., 2020. Matrix-free subcell residual distribution for bernstein finite element discretizations of linear advection equations. Comput. Methods Appl. Mech. Eng. 359, Article ID 112658, 26 p

    work page 2020

  26. [26]

    Amultidimensionalupwindschemewithnocrosswinddiffusion,in:Finiteelementmethodsforconvection dominated flows (Papers, Winter Ann

    Hughes,T.J.R.,Brooks,A.,1979. Amultidimensionalupwindschemewithnocrosswinddiffusion,in:Finiteelementmethodsforconvection dominated flows (Papers, Winter Ann. Meeting Amer. Soc. Mech. Engrs., New York, 1979). Amer. Soc. Mech. Engrs. (ASME), New York, pp. 19–35

    work page 1979

  27. [27]

    On spurious oscillations at layers diminishing (sold) methods for convection-diffusion equations

    John, V., Knobloch, P., 2007. On spurious oscillations at layers diminishing (sold) methods for convection-diffusion equations. i: A review. Comput. Methods Appl. Mech. Eng. 196, 2197–2215

    work page 2007

  28. [28]

    On spurious oscillations at layers diminishing (sold) methods for convection-diffusion equations

    John, V., Knobloch, P., 2008. On spurious oscillations at layers diminishing (sold) methods for convection-diffusion equations. ii: Analysis for𝑝 1 and𝑞 1 finite elements. Comput. Methods Appl. Mech. Eng. 197, 1997–2014

    work page 2008

  29. [29]

    Finiteelementmethodsforlinearhyperbolicproblems

    Johnson,C.,Nävert,U.,Pitkäranta,J.,1984. Finiteelementmethodsforlinearhyperbolicproblems. Comput.MethodsAppl.Mech.Eng.45, 285–312

    work page 1984

  30. [30]

    A convergent FEM-DG method for the compressible Navier-Stokes equations

    Karper, T.K., 2013. A convergent FEM-DG method for the compressible Navier-Stokes equations. Numer. Math. 125, 441–510

    work page 2013

  31. [31]

    Monolithicconvexlimitingforcontinuousfiniteelementdiscretizationsofhyperbolicconservationlaws

    Kuzmin,D.,2020. Monolithicconvexlimitingforcontinuousfiniteelementdiscretizationsofhyperbolicconservationlaws. Comput.Methods Appl. Mech. Engrg. 361, 112804, 28

    work page 2020

  32. [32]

    Locallybound-preservingenrichedGalerkinmethodsforthelinearadvectionequation

    Kuzmin,D.,Hajduk,H.,Rupp,A.,2020. Locallybound-preservingenrichedGalerkinmethodsforthelinearadvectionequation. Comput.& Fluids 205, 104525, 15

    work page 2020

  33. [33]

    Kwon,Y.S.,Novotný,A.,2021.ConstructionofweaksolutionstocompressibleNavier-Stokesequationswithgeneralinflow/outflowboundary conditions via a numerical approximation. Numer. Math. 149, 717–778

    work page 2021

  34. [34]

    High-resolution conservative algorithms for advection in incompressible flow

    Leveque, R.J., 1996. High-resolution conservative algorithms for advection in incompressible flow. SIAM J. Numer. Anal. 33, 627–665

    work page 1996

  35. [35]

    NewerrorestimatesofnonconformingmixedfiniteelementmethodsfortheStokesproblem

    Li,M.,Mao,S.,Zhang,S.,2014. NewerrorestimatesofnonconformingmixedfiniteelementmethodsfortheStokesproblem. Math.Methods Appl. Sci. 37, 937–951

    work page 2014

  36. [36]

    Optimal error estimates for discontinuous galerkin methods based on upwind-biased fluxes for linear hyperbolic equations

    Meng, X., Shu, C.W., Wu, B., 2016. Optimal error estimates for discontinuous galerkin methods based on upwind-biased fluxes for linear hyperbolic equations. Math. Comput. 85, 1225–1261

    work page 2016

  37. [37]

    MFEM: Modular finite element methods [Software].mfem.org

    mfem, . MFEM: Modular finite element methods [Software].mfem.org. doi:10.11578/dc.20171025.1248

    work page doi:10.11578/dc.20171025.1248

  38. [38]

    A rational finite element basis

    Wachspress, E.L., 1975. A rational finite element basis. volume Vol. 114 ofMathematics in Science and Engineering. Academic Press, Inc. [Harcourt Brace Jovanovich, Publishers], New York-London

    work page 1975

  39. [39]

    TheL2-norm stability analysis of runge-kutta discontinuous galerkin methods for linear hyperbolic equations

    Xu, Y., Zhang, Q., Shu, C.W., Wang, H., 2019. TheL2-norm stability analysis of runge-kutta discontinuous galerkin methods for linear hyperbolic equations. SIAM J. Numer. Anal. 57, 1574–1601

    work page 2019

  40. [40]

    Fully multidimensional flux-corrected transport algorithms for fluids

    Zalesak, S.T., 1979. Fully multidimensional flux-corrected transport algorithms for fluids. J. Comput. Phys. 31, 335–362

    work page 1979

  41. [41]

    Theweakgalerkinmethodforlinearhyperbolicequation

    Zhai,Q.,Zhang,R.,Malluwawadu,N.,Hussain,S.,2018. Theweakgalerkinmethodforlinearhyperbolicequation. Commun.Comput.Phys. 24, 152–166. S. Mao, M. Zhang:Preprint submitted to ElsevierPage 19 of 19 CR FEM for Time-dependent transport equations −1 1 1 𝑥= 0 −1 1−1 −3 Figure 1.The CR interpolation attains a value of−3at the lower-left vertex of the blue triangle....

    work page 2018