Pith. sign in

REVIEW 3 major objections 6 minor 73 references

Discontinuous Galerkin methods for the complete stochastic Euler equations

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

Pith's one-line read Entropy-stable discontinuous Galerkin discretizations of the stochastically forced complete Euler equations are shown to converge in law to dissipative martingale solutions, with relative-energy rates $O(h^{1/2})$ and $O(h)$ near pathwise…

desk verdict First convergence result for entropy-stable DG on the complete stochastic Euler system, but every theorem rests on an unproved uniform-in-h bound and a bounded-domain strong-solution existence that is only cited. read the letter →

arxiv 2412.07613 v1 pith:LIJ263JQ submitted 2024-12-10 math.NA cs.NAmath.APmath.PR

classification math.NAcs.NAmath.APmath.PR MSC 65M6060H1535Q3165M1265M08
keywords stochasticEulerequationsdiscontinuousGalerkinentropystabilitydissipativemartingalesolutionsrelativefinitevolumecompactnessconvergencerates
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 establishes a first convergence theory for high-order entropy-stable discontinuous Galerkin discretizations of the complete compressible Euler equations with stochastic forcing in the momentum equation. The target object is a dissipative martingale solution, a measure-valued weak solution adapted to the driving Wiener process whose oscillations and concentrations are recorded by defect measures. The main result says that, as the mesh size $h\to 0$, the numerical laws converge, up to a subsequence, to such a solution, provided the discrete density stays away from vacuum and the discrete total energy stays bounded up to a stopping time. If a smooth pathwise strong solution exists, the expected $L^1$ relative energy error between the numerical and exact solution decays at least as $h^{1/2}$ for the finite-volume case and as $h$ for higher polynomial degrees, matching the deterministic theory. This matters because stochastic Euler-type models are used to represent turbulence and uncertainty, and until now no convergence result existed for such methods on the full stochastic system.

What carries the argument

The load-bearing object is the entropy-stable discontinuous Galerkin spectral element scheme (2.29): a nodal finite-element space on Gauss-Lobatto quadrature points written in summation-by-parts form, in which the volume integral is replaced by an entropy-conservative two-point numerical flux via flux differencing and interfaces are closed with a local Lax-Friedrichs flux, so the semidiscrete system dissipates the mathematical entropy. This scheme is applied to the total-energy formulation (2.22) of the stochastic Euler equations, with the noise entering only through the momentum and total-energy equations. Convergence is mediated by the dissipative martingale solution concept of Definition 3.1, where oscillations and concentrations are recorded by parametrised measures and defect measures, and the error analysis runs through the relative entropy, or ballistic free energy, inequality of Proposition 5.1, which compares the discrete state to the pathwise strong solution and reduces the error to a Gronwall-type estimate.

What would settle it

Run the $p=0$ local Lax-Friedrichs scheme on the one-dimensional SOD shock-tube initial data of Section 6.1.3 with the same noise protocol and 1000 samples, and record, for each mesh size $h=64,128,256,512$, the first time at which $\varrho_h$ falls below $1/K$ or $E_h$ exceeds $K$; if these hitting times converge to zero with positive probability as $h\to0$, then assumption (3.12) has no positive-time regime for that test case, whereas if they remain bounded away from zero, comparing the measured $E_2$ with the predicted $O(h^{1/2})$ bound at $T=0.2$ settles the rate claim of Theorem 3.8(a) empirically.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central claim is that the entropy-stable DGSEM scheme for the complete stochastic Euler system is consistent and convergent even though the underlying PDE is hyperbolic, stochastically forced, and may develop discontinuities. Under the a priori bound that the discrete density never approaches vacuum and the discrete total energy never blows up, uniformly in the mesh size and up to a stopping time with $\mathbb{P}(\tau>0)=1$, the laws of the numerical solutions are tight, and every limit point, after a Skorokhod-type change of probability space, is a dissipative martingale solution in the sense of the paper's Definition 3.1. When a local pathwise strong solution exists on the same stochastic basis, the expected $L^1$ relative energy between the numerical and exact states is $O(h^{1/2})$ for polynomial degree $p=0$ and $O(h)$ for $p\ge1$ under an additional unverified nodal Lipschitz condition, reproducing the deterministic convergence rates. The numerical experiments on smooth density waves in one and two dimensions show the scheme meeting or exceeding these rates, and the paper presents the first simulations for this stochastic Euler system in regimes where discontinuities may develop.

Load-bearing premise

The entire convergence and rate statements are conditional on (3.12): uniformly in the mesh size, the discrete density never drops below $1/K$ and the discrete total energy never exceeds $K$ until some stopping time, and this bound is assumed rather than derived from the scheme, with an extra unverified smoothness condition on the numerical nodal values needed for the $p\ge1$ rate.

Editorial extensions

If this is right

  • Under the stopping-time bound, the entropy-stable DGSEM scheme yields, along a subsequence, a dissipative martingale solution of the complete stochastic Euler system, providing the first convergence result for a high-order method on this system.
  • When a local pathwise strong solution exists, the expected $L^1$ relative energy at the final time is $O(h^{1/2})$ for $p=0$ and $O(h)$ for $p\ge1$ under the nodal Lipschitz assumption, giving concrete rates that can be checked numerically.
  • Because the scheme is entropy dissipative by construction, the limit inherits the entropy inequality and the total-energy balance, so the numerical method preserves the thermodynamic structure rather than merely conserving mass, momentum, and energy.
  • The finite-volume case $p=0$ is included as a special case, so the result also supplies convergence of local Lax-Friedrichs finite volume schemes for the stochastic Euler equations.
  • Outside the lifespan of a strong solution, the only guarantee is convergence in law to a dissipative martingale solution; the stopping-time formulation is what makes the a priori bounds and hence the whole argument possible.

Reading between the lines

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

  • A natural testable extension is to replace the unverified nodal Lipschitz condition (3.15) by an explicit limiter or slope-reconstruction step; if that provably enforces (3.15), the $O(h)$ rate for $p\ge1$ becomes a property of the algorithm rather than an assumption on its output.
  • The numerical tables show fourth-order relative-energy convergence on smooth data, far above the proved rates; this suggests the analysis is limited by the worst-case consistency estimates and that a regularity-aware error estimate could close the gap.
  • Because Gaussian noise makes uniform-in-time $L^\infty$ bounds impossible, the stopping time is not merely technical: any extension to deterministic times will likely need a noise-strength-dependent or path-dependent horizon, or a new a priori bound not present in this paper.
  • For genuinely infinite-dimensional noise, the treatment of truncating the Wiener process to finitely many modes adds a tail error controlled by the truncated basis; this yields a concrete recipe for an implementable method with a combined discretization-plus-truncation error estimate.
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 / 6 minor

Summary. The paper studies an entropy-dissipative discontinuous Galerkin spectral element method (DGSEM), including the finite-volume case, for the complete stochastic Euler equations with multiplicative noise in the momentum equation. The main theoretical results are conditional: under a uniform-in-h a priori bound on the numerical density and total energy up to a stopping time (3.12), the authors prove consistency and convergence in law, up to a subsequence, to a dissipative martingale solution in the sense of Definition 3.1 (Theorem 3.6). In the lifespan of a sufficiently regular pathwise strong solution, they prove relative-entropy error estimates of order h^{1/2} for polynomial degree p=0 and order h for p≥1, the latter under an additional Lipschitz condition on the nodal values (Theorem 3.8, Eqs. (3.14) and (3.16)). The paper also presents numerical experiments in one and two space dimensions using the Trixi.jl framework, reporting observed convergence rates for smooth and discontinuous test cases and studying the influence of noise strength.

Significance. If the results hold, this is the first convergence analysis for a high-order DG discretization of the complete stochastic Euler equations, and it extends the deterministic consistency framework of Lukáčová-Medvid'ová and Öffner to the stochastic setting using dissipative martingale solutions. The paper is carefully structured, states its hypotheses explicitly, and provides a concrete numerical implementation with a reproducible open-source software framework. The central convergence theorem is honest about its conditional character, and the relative-entropy argument gives a natural stochastic analogue of deterministic weak-strong uniqueness. However, the two main theorems rest entirely on a priori bounds that are not derived from the scheme, and one step in the proof of Theorem 3.6 (conversion of an energy inequality into the equality required by Definition 3.1) is not rigorously justified. These issues are load-bearing for the paper's central claims.

major comments (3)
  1. [§3.3, Eq. (3.12); Theorems 3.6 and 3.8] Theorems 3.6 and 3.8 are both conditioned on the uniform-in-h bound (3.12), but this bound is neither proved for the scheme nor shown to hold on a probability-one set. The discrete entropy inequality (4.4) and the energy balance (4.5) provide only L^p-type and L^2(U) control, not the uniform L^∞ upper bound on E_h and the uniform positive lower bound on ϱ_h required by (3.12); in particular, the scheme is not shown to be positivity preserving or to preclude blow-up before the common stopping time τ. Section 6 in fact reports stability failures for the Kelvin–Helmholtz test once the noise strength exceeds 1. Because τ must be common to all mesh refinements, the theorem has no content on noise paths for which some refinement violates (3.12) arbitrarily early. The authors should either prove (3.12) for a stabilized (e.g., limited) version of the scheme, or state the results as strictly conditional and discuss the verifiability and plausibility of the assumption in detail.
  2. [§3.3, Theorem 3.8(b), Eq. (3.15); Remark 4.1] The O(h) rate for p≥1 is claimed only under the additional Lipschitz condition (3.15) on the nodal values of the numerical solution within each element. This is an unverified a priori assumption about the numerical solution: the entropy-dissipative DGSEM is not shown to satisfy it, and Remark 4.1 explicitly notes that limiters are disregarded and that with limiters only finite-volume-type rates are expected. Since (3.15) is needed to obtain the strong consistency estimate (4.7) from [52], the O(h) rate is not a property of the scheme as analyzed unless (3.15) is established or at least numerically verified. The authors should prove (3.15) under suitable assumptions on the mesh, time step, and data, or downgrade the claim to the p=0 result.
  3. [§4.3, after Eq. (4.23)] The passage from the energy inequality (4.23) to the energy equality required by Definition 3.1(l) is not rigorously justified. The text first says "Performing the limit εm → 0 yields an energy inequality" (no εm has been introduced at that point), and then states that augmenting Rpress by a spatially homogeneous h(t)dx converts the inequality into equality. The authors do not verify that the augmented measure has the required measurability, adaptedness, and weak-* compactness properties, nor that the energy deficit can always be represented in this way while preserving the momentum equation (3.6). This step is load-bearing because Definition 3.1(l) demands equality. A detailed construction of the augmented measure, or a reformulation of the theorem with an energy inequality in the definition of the limit object, is needed.
minor comments (6)
  1. [Abstract and Section 1] The abstract contains a typo: "we proof" should be "we prove", and "The results built" should be "The results build".
  2. [§2.3, Eq. (2.20)] The solution space V_h is defined using L^1(O), although its elements are piecewise polynomials and hence bounded; L^∞(O) or a more standard broken polynomial space notation would be clearer.
  3. [§2.3, after Eq. (2.28)] The text reads "I4 is the 4 × 4-identy matrix"; this should be "identity matrix".
  4. [§3.1, Definition 3.1] The dummy-variable space A in (3.1) is defined with m' ∈ R^3, while the paper explicitly focuses on two space dimensions; the dimension should be made consistent, or the notation should be explained as generic.
  5. [§4.3, last paragraph] The phrase "Performing the limit εm → 0" appears without any prior definition of εm; if the intended limit is m → ∞, it should be stated correctly.
  6. [§6.1.1 and Table 5] The text contains a German fragment "f¨ urt" instead of "for", and Table 5 writes "162", "322", etc., where the intended entries are 16^2, 32^2, and so on; these should be formatted correctly.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the convergence and rate theorems are conditional on an explicit a priori bound, but they do not reduce to their inputs by construction.

full rationale

The paper's central claim is that, under the explicit uniform-in-h bound (3.12) (positive density and bounded total energy up to a stopping time), the entropy-stable DGSEM solutions converge in law to a dissipative martingale solution of the stochastic Euler system, and that the relative-energy error against a strong solution is O(h^{1/2}) for p=0 or O(h) for p>=1 under the additional nodal-Lipschitz assumption (3.15). This is a conditional theorem: (3.12) and (3.15) are hypotheses, not consequences, and the paper says so explicitly ('We work under the following hypothesis...'). A conditional statement is not circular merely because its hypothesis is unverified; the convergence proof supplies consistency estimates (4.2)-(4.8), compactness via Jakubowski-Skorokhod, and a relative-entropy inequality (Proposition 5.1) that are not the same as the assumed bounds. The consistency estimates and the dissipative-solution concept are imported from the authors' own prior works [52] and [56], but those are published, parameter-free results with stated assumptions that do not include the stochastic convergence target; therefore, by the independence rule, they are real evidence rather than circularity. The step of 'augmenting' Rpress by a spatially homogeneous measure to convert the energy inequality (4.23) into the equality required by Definition 3.1 is a standard defect-measure completion and does not alter the momentum or entropy equations, so it is not a self-definitional reduction. The paper also flags genuine scope limitations: (3.12) is not shown to hold for any fixed deterministic horizon, the scheme is not shown to be positivity-preserving, and the numerical section reports instability for large noise in the Kelvin-Helmholtz test. These are correctness and robustness gaps, not circular derivations; no equation is defined in terms of the result it is used to prove.

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

The central claims are conditional on several domain assumptions rather than on fitted parameters. There are no newly invented physical entities: the dissipative martingale solution with defect measures Rconv, Rpress and Young measure V was introduced in Moyo [56]. The scheme itself is a standard entropy-stable DGSEM with a stochastic forcing term. The free parameter list is empty because the constant K and the stopping time t are assumed inputs, not fitted values.

assumptions (6)
  • domain assumption A priori bound (3.12): inf_{[0,t]×O} ϱh ≥ 1/K and sup_{[0,t]×O} E_h ≤ K P-a.s., uniformly in h, for some stopping time t with P(t>0)=1.
    Central to Theorem 3.6 and Theorem 3.8; not derived from the scheme. The paper argues it is plausible only for short random times, but the stopping time is not quantified.
  • domain assumption Existence of a local strong pathwise solution [r, Θ, v] on bounded domain O with no-flux boundary conditions up to stopping time s (Definition 3.3).
    Required in Theorem 3.8. The paper says 'we expect' existence and cites [47] for the whole space, but does not provide a proof for the bounded-domain case.
  • ad hoc to paper Lipschitz continuity of nodal values within each element, (3.15): |U_h(x_j,t∧t) - U_h(x_i,t∧t)| ≤ ch for all x_i, x_j in each element.
    Introduced specifically to obtain the O(h) rate for p≥1 in Theorem 3.8(b). No justification is given that realistic DG solutions satisfy it; the authors note limiters are not considered.
  • domain assumption Deterministic consistency estimates (4.7) and (4.8) from [52] and [53] hold for the stochastic scheme.
    The paper asserts 'line by line' following [52] and [53] without re-deriving them. These estimates control the consistency error terms e_h in the stochastic setting.
  • domain assumption Weak-strong uniqueness principle for dissipative martingale solutions, Theorem 3.4 from [56].
    Used in the proof of Theorem 3.8 to control the relative entropy production Q. The result is from an author's prior paper and is not re-proved here.
  • standard math Itô formula, martingale representation, and Jakubowski's extension of the Skorokhod representation theorem.
    Standard tools in stochastic PDE analysis, used without proof in Section 4.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Discontinuous Galerkin methods for the complete stochastic Euler equations." pith.science (2026). https://pith.science/paper/LIJ263JQ

@misc{pith2026241207613,
  author       = {Pith},
  title        = {Pith review of: Discontinuous Galerkin methods for the complete stochastic Euler equations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LIJ263JQ}},
  note         = {Machine review of arXiv:2412.07613}
}
abstract

In recent years, stochastic effects have become increasingly relevant for describing fluid behaviour, particularly in the context of turbulence. The most important model for inviscid fluids in computational fluid dynamics are the Euler equations of gas dynamics which we focus on in this paper. To take stochastic effects into account, we incorporate a stochastic forcing term in the momentum equation of the Euler system. To solve the extended system, we apply an entropy dissipative discontinuous Galerkin spectral element method including the Finite Volume setting, adjust it to the stochastic Euler equations and analyze its convergence properties. Our analysis is grounded in the concept of dissipative martingale solutions, as recently introduced by Moyo (J. Diff. Equ. 365, 408-464, 2023). Assuming no vacuum formation and bounded total energy, we proof that our scheme converges in law to a dissipative martingale solution. During the lifespan of a pathwise strong solution, we achieve convergence of at least order 1/2, measured by the expected $L^1$ norm of the relative energy. The results built a counterpart of corresponding results in the deterministic case. In numerical simulations, we show the robustness of our scheme, visualise different stochastic realizations and analyze our theoretical findings.

Figures

Figures reproduced from arXiv: 2412.07613 by the authors.

Figure 1
Figure 1. Different noise (left) and different noise strength(right) 6.2. Two dimension. We consider the two-dimensional stochastic compressible Euler equations (2.23) with the right hand side h(U) = µ 2 ϱ   0 0 0 1   , g(U)dW = µ   0 ϱ dβ1 ϱ dβ2 m1 dβ1 + m2 dβ2   . where µ denotes again a stochastic forcing factor. The domain is divided into n × n uniform quads for the elements (n ∈ N) [PITH_FULL_IMAGE:figur… view at source ↗
Figure 2
Figure 2. Kelvin-Helmholtz and influence of the noise references for instance in [44, 46, 51, 60]. However, our analysis is relying on the consistency estimation where no regularity have been assumed. In the future, we plan to derive new error estimates using the consistency analysis and transfer it also from the purely deterministic setting to the stochastic one. Meanwhile, we have numerically analyzed the stochastic Euler e… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

73 extracted references · 69 canonical work pages

  1. [52]

    Luk´ aˇ cov´ a-Medvid’ov´ a and P.¨Offner

    M. Luk´ aˇ cov´ a-Medvid’ov´ a and P.¨Offner. Convergence of discontinuous Galerkin schemes for the Euler equations via dissi- pative weak solutions. Appl. Math. Comput. , 436:22, 2023. Id/No 127508

  2. [53]

    Luk´ aˇ cov´ a-Medvid’ov´ a, B

    M. Luk´ aˇ cov´ a-Medvid’ov´ a, B. She, and Y. Yuan. Error estimates of the Godunov method for the multidimensional com- pressible Euler system. J. Sci. Comput. , 91(3):27, 2022. Id/No 71

  3. [1]

    Abgrall, P

    R. Abgrall, P. Bacigaluppi, and S. Tokareva. A high-order nonconservative approach for hyperbolic equations in fluid dynamics. Comput. Fluids , 169:10–22, 2018

  4. [2]

    Abgrall and S

    R. Abgrall and S. Karni. A comment on the computation of non-conservative products. J. Comput. Phys. , 229(8):2759– 2763, 2010

  5. [3]

    Abgrall, M

    R. Abgrall, M. Luk´ aˇ cova-Medvid’ov´ a, and P.¨Offner. On the convergence of residual distribution schemes for the compress- ible Euler equations via dissipative weak solutions. Math. Models Methods Appl. Sci. , 33(1):139–173, 2023

  6. [4]

    Abgrall and S

    R. Abgrall and S. Mishra. Uncertainty qualification for hyperbolic systems of conservation laws. In Handbook on numerical methods for hyperbolic problems. Applied and modern issues , pages 507–544. Amsterdam: Elsevier/North Holland, 2017

  7. [5]

    Abgrall, P

    R. Abgrall, P. ¨Offner, and H. Ranocha. Reinterpretation and extension of entropy correction terms for residual distribution and discontinuous Galerkin schemes: application to structure preserving discretization. J. Comput. Phys. , 453:24, 2022. Id/No 110955

  8. [6]

    Benzoni-Gavage and D

    S. Benzoni-Gavage and D. Serre. Multi-dimensional hyperbolic partial differential equations. First-order systems and ap- plications. Oxford Math. Monogr. Oxford: Oxford University Press, 2007

Show all 73 references
  1. [7]

    Bezanson, A

    J. Bezanson, A. Edelman, S. Karpinski, and V. B. Shah. Julia: A fresh approach to numerical computing. SIAM Review, 59(1):65–98, 2017

  2. [8]

    B. Birnir. The Kolmogorov-Obukhov statistical theory of turbulence. J. Nonlinear Sci. , 23(4):657–688, 2013

  3. [9]

    D. Breit. Existence theory for generalized Newtonian fluids . Math. Sci. Eng. Amsterdam: Elsevier/Academic Press, 2017

  4. [10]

    Breit and A

    D. Breit and A. Dodgson. Space-time approximation of local strong solutions to the 3d stochastic Navier-Stokes equations. Comput. Methods Appl. Math. , 24(3):577–597, 2024

  5. [11]

    Breit, E

    D. Breit, E. Feireisl, and M. Hofmanov´ a.Stochastically forced compressible fluid flows, volume 3 of De Gruyter Ser. Appl. Numer. Math. Berlin: De Gruyter, 2018

  6. [12]

    Breit, E

    D. Breit, E. Feireisl, and M. Hofmanov´ a. Dissipative solutions and semiflow selection for the complete Euler system. Commun. Math. Phys. , 376(2):1471–1497, 2020

  7. [13]

    Breit, E

    D. Breit, E. Feireisl, and M. Hofmanov´ a. On solvability and ill-posedness of the compressible Euler system subject to stochastic forces. Anal. PDE, 13(2):371–402, 2020

  8. [14]

    Breit, T

    D. Breit, T. C. Moyo, A. Prohl, and J. Wichmann. Mean square temporal error estimates for the 2d stochastic navier-stokes equations with transport noise. arXiv preprint arXiv:2305.10999 , 2023

  9. [15]

    Brenier, C

    Y. Brenier, C. De Lellis, and L. j. Sz´ ekelyhidi. Weak-strong uniqueness for measure-valued solutions. Commun. Math. Phys., 305(2):351–361, 2011

  10. [16]

    Bˇ rezina and E

    J. Bˇ rezina and E. Feireisl. Measure-valued solutions to the complete Euler system. J. Math. Soc. Japan , 70(4):1227–1245, 2018

  11. [17]

    Brze´ zniak, E

    Z. Brze´ zniak, E. Hausenblas, and P. A. Razafimandimby. Stochastic reaction-diffusion equations driven by jump processes. Potential Anal., 49(1):131–201, 2018

  12. [18]

    Brze´ zniak and M

    Z. Brze´ zniak and M. Ondrej´ at. Strong solutions to stochastic wave equations with values in Riemannian manifolds. J. Funct. Anal., 253(2):449–481, 2007

  13. [19]

    Buckmaster, S

    T. Buckmaster, S. Shkoller, and V. Vicol. Shock formation and vorticity creation for 3d Euler. Commun. Pure Appl. Math., 76(9):1965–2072, 2023

  14. [20]

    Chaudhary and U

    A. Chaudhary and U. Koley. A convergent finite volume scheme for the stochastic barotropic compressible Euler equations. ESAIM, Math. Model. Numer. Anal. , 57(6):3403–3437, 2023

  15. [21]

    Chen and C.-W

    T. Chen and C.-W. Shu. Entropy stable high order discontinuous Galerkin methods with suitable quadrature rules for hyperbolic conservation laws. J. Comput. Phys. , 345:427–461, 2017

  16. [22]

    Chiodaroli, C

    E. Chiodaroli, C. De Lellis, and O. Kreml. Global ill-posedness of the isentropic system of gas dynamics. Commun. Pure Appl. Math. , 68(7):1157–1190, 2015

  17. [23]

    Chiodaroli, E

    E. Chiodaroli, E. Feireisl, and F. Flandoli. Ill posedness for the full Euler system driven by multiplicative white noise. Indiana Univ. Math. J. , 70(4):1267–1282, 2021

  18. [24]

    Christ, D

    S. Christ, D. Schwabeneder, C. Rackauckas, M. K. Borregaard, and T. Breloff. Plots.jl — a user extendable plotting API for the Julia programming language. Journal of Open Research Software , 2023

  19. [25]

    Cockburn, G

    B. Cockburn, G. E. Karniadakis, and C.-W. Shu. The development of discontinuous Galerkin methods. In Discontinuous Galerkin methods. Theory, computation and applications. 1st international symposium on DGM, Newport, RI, USA, May 24–26, 1999 , pages 3–50. Berlin: Springer, 2000

  20. [26]

    Da Prato and J

    G. Da Prato and J. Zabczyk. Stochastic equations in infinite dimensions , volume 44 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, 1992

  21. [27]

    Da Prato and J

    G. Da Prato and J. Zabczyk. Stochastic equations in infinite dimensions , volume 44 of Encycl. Math. Appl. Cambridge etc.: Cambridge University Press, 1992

  22. [28]

    De Lellis and L

    C. De Lellis and L. j. Sz´ ekelyhidi. On admissibility criteria for weak solutions of the Euler equations. Arch. Ration. Mech. Anal., 195(1):225–260, 2010

  23. [29]

    Debussche, N

    A. Debussche, N. Glatt-Holtz, and R. Temam. Local martingale and pathwise solutions for an abstract fluids model. Physica D, 240(14-15):1123–1144, 2011

  24. [30]

    R. J. Diperna. Measure-valued solutions to conservation laws. Arch. Ration. Mech. Anal. , 88:223–270, 1985

  25. [31]

    R. J. DiPerna and A. J. Majda. Oscillations and concentrations in weak solutions of the incompressible fluid equations. Commun. Math. Phys. , 108:667–689, 1987

  26. [32]

    Feireisl, C

    E. Feireisl, C. Klingenberg, O. Kreml, and S. Markfelder. On oscillatory solutions to the complete Euler system. J. Differ. Equations, 269(2):1521–1543, 2020. 28 DISCONTINUOUS GALERKIN METHODS FOR THE COMPLETE STOCHASTIC EULER EQUATIONS

  27. [33]

    Feireisl, M

    E. Feireisl, M. Luk´ aˇ cov´ a-Medvid’ov´ a, and H. Mizerov´ a. Convergence of finite volume schemes for the Euler equations via dissipative measure-valued solutions. Found. Comput. Math. , 20(4):923–966, 2020

  28. [34]

    Feireisl, M

    E. Feireisl, M. Luk´ aˇ cov´ a-Medviˇdov´ a, H. Mizerov´ a, and B. She.Numerical analysis of compressible fluid flows , volume 20 of MS&A, Model. Simul. Appl. Cham: Springer, 2021

  29. [35]

    T. C. Fisher, M. H. Carpenter, J. Nordstr¨ om, N. K. Yamaleev, and C. Swanson. Discretely conservative finite-difference formulations for nonlinear conservation laws in split form: theory and boundary conditions.J. Comput. Phys., 234:353–375, 2013

  30. [36]

    U. S. Fjordholm, S. Mishra, and E. Tadmor. On the computation of measure-valued solutions. Acta Numerica, 25:567–679, 2016

  31. [37]

    G. J. Gassner, A. R. Winters, and D. A. Kopriva. Split form nodal discontinuous Galerkin schemes with summation-by-parts property for the compressible Euler equations. J. Comput. Phys. , 327:39–66, 2016

  32. [38]

    Glaubitz, H

    J. Glaubitz, H. Ranocha, A. R. Winters, M. Schlottke-Lakemper, P. ¨Offner, and G. Gassner. Generalized up- wind summation-by-parts operators and their application to nodal discontinuous galerkin methods. arXiv preprint arXiv:2406.14557, 2024

  33. [39]

    S. K. Godunov. A difference method for numerical calculation of discontinuous solutions of the equations of hydrodynamics. Mat. Sb., Nov. Ser. , 47:271–306, 1959

  34. [40]

    Harten, P

    A. Harten, P. D. Lax, and B. van Leer. On upstream differencing and Godunov-type schemes for hyperbolic conservation laws. SIAM Rev., 25:35–61, 1983

  35. [41]

    Hennemann, A

    S. Hennemann, A. M. Rueda-Ram ´ ırez, F. J. Hindenlang, and G. J. Gassner. A provably entropy stable subcell shock capturing approach for high order split form DG for the compressible Euler equations. J. Comput. Phys. , 426:31, 2021. Id/No 109935

  36. [42]

    J. S. Hesthaven and T. Warburton. Nodal discontinuous Galerkin methods. Algorithms, analysis, and applications , vol- ume 54 of Texts Appl. Math. New York, NY: Springer, 2008

  37. [43]

    Hofmanov´ a

    M. Hofmanov´ a. Degenerate parabolic stochastic partial differential equations. Stochastic Processes Appl. , 123(12):4294– 4336, 2013

  38. [44]

    Huang and C.-W

    J. Huang and C.-W. Shu. Error estimates to smooth solutions of semi-discrete discontinuous Galerkin methods with quadrature rules for scalar conservation laws. Numer. Methods Partial Differ. Equations , 33(2):467–488, 2017

  39. [45]

    Jakubowski

    A. Jakubowski. The almost sure Skorokhod representation for subsequences in nonmetric spaces. Teor. Veroyatn. Primen., 42(1):209–216, 1997

  40. [46]

    M. Jiao, Y. Jiang, C.-W. Shu, and M. Zhang. Optimal error estimates to smooth solutions of the central discontinuous Galerkin methods for nonlinear scalar conservation laws. ESAIM, Math. Model. Numer. Anal. , 56(4):1401–1435, 2022

  41. [47]

    J. U. Kim. On the stochastic quasi-linear symmetric hyperbolic system. J. Differ. Equations , 250(3):1650–1684, 2011

  42. [48]

    Kr¨ oner and W

    D. Kr¨ oner and W. M. Zajaczkowski. Measure-valued solutions of the Euler equations for ideal compressible polytropic fluids. Math. Methods Appl. Sci. , 19(3):235–252, 1996

  43. [49]

    Kuzmin and H

    D. Kuzmin and H. Hajduk. Property-Preserving Numerical Schemes for Conservation Laws . World Scientific, 2023

  44. [50]

    Kuzmin, M

    D. Kuzmin, M. Luk´ acova-Medvid’ov´ a, and P.¨Offner. Consistency and convergence of flux-corrected finite element methods for nonlinear hyperbolic problems. arXiv preprint arXiv:2308.14872 , 2023

  45. [51]

    Liu, C.-W

    Y. Liu, C.-W. Shu, and M. Zhang. Optimal error estimates of the semidiscrete discontinuous Galerkin methods for two dimensional hyperbolic equations on Cartesian meshes using P k elements. ESAIM, Math. Model. Numer. Anal., 54(2):705– 726, 2020

  46. [54]

    Luk´ aˇ cov´ a-Medvid’ov´ a and Y

    M. Luk´ aˇ cov´ a-Medvid’ov´ a and Y. Yuan. Convergence of first-order finite volume method based on exact Riemann solver for the complete compressible Euler equations. Numer. Methods Partial Differ. Equations , 39(5):3777–3810, 2023

  47. [55]

    Mikulevicius and B

    R. Mikulevicius and B. L. Rozovskii. Stochastic Navier-Stokes equations for turbulent flows. SIAM J. Math. Anal. , 35(5):1250–1310, 2004

  48. [56]

    T. C. Moyo. Dissipative solutions and Markov selection to the complete stochastic Euler system. J. Differ. Equations , 365:408–464, 2023

  49. [57]

    Neustupa

    J. Neustupa. Measure-valued solutions of the Euler and Navier-Stokes equations for compressible barotropic fluids. Math. Nachr., 163:217–227, 1993

  50. [58]

    P. ¨Offner. Zweidimensionale klassische und diskrete orthogonale polynome und ihre anwendung auf spektrale methoden zur l¨ osung von hyperbolischen erhaltungsgleichungen. PhD thesis, 2015

  51. [59]

    P. ¨Offner. Approximation and stability properties of numerical methods for hyperbolic conservation laws . Wiesbaden: Springer Spektrum, 2023

  52. [60]

    ¨Offner and H

    P. ¨Offner and H. Ranocha. Error boundedness of discontinuous Galerkin methods with variable coefficients.J. Sci. Comput., 79(3):1572–1607, 2019

  53. [61]

    Ondrej´ at, A

    M. Ondrej´ at, A. Prohl, and N. J. Walkington. Numerical approximation of nonlinear SPDE’s. Stoch. Partial Differ. Equ., Anal. Comput., 11(4):1553–1634, 2023

  54. [62]

    Rackauckas and Q

    C. Rackauckas and Q. Nie. DifferentialEquations.jl – A performant and feature-rich ecosystem for solving differential equations in Julia. Journal of Open Research Software , 5(1):15, 2017

  55. [63]

    H. Ranocha. Comparison of some entropy conservative numerical fluxes for the Euler equations.J. Sci. Comput., 76(1):216– 242, 2018. DISCONTINUOUS GALERKIN METHODS FOR THE COMPLETE STOCHASTIC EULER EQUATIONS 29

  56. [64]

    Ranocha, M

    H. Ranocha, M. Schlottke-Lakemper, A. R. Winters, E. Faulhaber, J. Chan, and G. J. Gassner. Adaptive numerical simulations with Trixi.jl: A case study of Julia for scientific computing. Proceedings of the JuliaCon Conferences, 1(1):77, 01 2022

  57. [65]

    Ranocha, A

    H. Ranocha, A. R. Winters, M. Schlottke-Lakemper, P. ¨Offner, J. Glaubitz, and G. J. Gassner. On the robustness of high- order upwind summation-by-parts methods for nonlinear conservation laws.Journal of Computational Physics , 520:113471, 2025

  58. [66]

    W. H. Reed and T. R. Hill. Triangular mesh methods for the neutron transport equation. Technical report, Los Alamos Scientific Lab., N. Mex.(USA), 1973

  59. [67]

    F. Renac. Entropy stable DGSEM for nonlinear hyperbolic systems in nonconservative form with application to two-phase flows. J. Comput. Phys. , 382:1–26, 2019

  60. [68]

    P. L. Roe. Approximate Riemann solvers, parameter vectors, and difference schemes. J. Comput. Phys. , 43:357–372, 1981

  61. [69]

    Schlottke-Lakemper, A

    M. Schlottke-Lakemper, A. R. Winters, H. Ranocha, and G. J. Gassner. A purely hyperbolic discontinuous Galerkin approach for self-gravitating gas dynamics. J. Comput. Phys. , 442:20, 2021. Id/No 110467

  62. [70]

    J. Smoller. Shock waves and reaction-diffusion equations. , volume 258 of Grundlehren Math. Wiss. New York: Springer- Verlag, 2nd ed. edition, 1994

  63. [71]

    E. F. Toro. Riemann solvers and numerical methods for fluid dynamics. A practical introduction . Berlin: Springer, 3rd ed. edition, 2009

  64. [72]

    Vandeven

    H. Vandeven. Family of spectral filters for discontinuous problems. J. Sci. Comput. , 6(2):159–192, 1991

  65. [73]

    Woodfield

    J. Woodfield. Strong stability preservation for stochastic partial differential equations. arXiv preprint arXiv:2411.11172 , 2024

Pith tools

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