pith. machine review for the scientific record. sign in

arxiv: 2604.02757 · v2 · submitted 2026-04-03 · ⚛️ physics.flu-dyn

Recognition: 2 theorem links

· Lean Theorem

Wave-appropriate reconstruction of compressible flows: physics-constrained acoustic dissipation and rank-1 entropy wave correction

Amareshwara Sainadh Chamarthi
Authors on Pith no claims yet

Pith reviewed 2026-05-13 18:25 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords wave-appropriate reconstructionacoustic dissipationentropy wave correctioncompressible flowscharacteristic decompositionDucros sensorTaylor-Green vortex
0
0 comments X

The pith

Optimizing one acoustic upwind parameter yields values that generalize without retuning from subsonic turbulence to hypersonic shocks.

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

The paper establishes a wave-appropriate reconstruction that decomposes the procedure into characteristic wave families, applying upwind bias only to acoustic waves while centralizing non-acoustic waves to reduce dissipation. It performs bounded scalar minimization treating the CFD solver as a black box, using an accuracy objective on subsonic inviscid Taylor-Green vortex subject to a stability constraint from the supersonic viscous case, to locate the minimal robust acoustic upwind parameter. The resulting values then apply directly to hypersonic flows with shocks and contacts; a separate rank-1 correction along the entropy eigenvector removes the need for explicit contact detectors and relies solely on the Ducros sensor, cutting wall time by 29-41 percent.

Core claim

Treating the solver as a black box in bounded minimization finds an optimal acoustic upwind parameter that generalizes across flow regimes from subsonic inviscid TGV to hypersonic flows with shocks, while a rank-1 entropy wave correction using the Ducros sensor removes the need for contact-discontinuity detectors and reduces wall time by 29-41% compared to full decomposition.

What carries the argument

Characteristic wave-family decomposition with a single tunable acoustic upwind parameter and a rank-1 update along the entropy right eigenvector.

If this is right

  • The optimal acoustic parameter requires no retuning when moving from subsonic turbulence to hypersonic flows with shocks and contacts.
  • The rank-1 entropy correction is limiter-agnostic and integrates directly into other schemes such as WENO.
  • Wall time drops by 29-41 percent relative to full characteristic decomposition.
  • Controlled acoustic bias applied only to normal momentum in KEP schemes removes spurious vortices in periodic shear layers.

Where Pith is reading between the lines

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

  • The framework independence from specific limiters or discretizations suggests the approach can be inserted into many existing CFD codes with limited code changes.
  • Similar black-box optimization on other reconstruction parameters could extend robustness to additional multi-physics or reacting-flow settings.

Load-bearing premise

The bounded optimization performed on two Taylor-Green vortex cases is sufficient to guarantee robustness and accuracy in all other regimes including those with strong shocks and arbitrary limiters.

What would settle it

A hypersonic simulation containing strong shocks and contact discontinuities that becomes unstable or loses accuracy when the optimized acoustic parameter and rank-1 entropy correction are used without further adjustment.

Figures

Figures reproduced from arXiv: 2604.02757 by Amareshwara Sainadh Chamarthi.

Figure 1
Figure 1. Figure 1: Wave-appropriate reconstruction algorithm. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Shock sensor components for underexpanded jet, reproduced from Sciacovelli et al. [ [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Optimizer convergence traces for the physics-constrained optimization of [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Subsonic inviscid Taylor-Green vortex (Ma = 0.1, p0 = 100, 643 grid, Sec. 7.1): time evolution of volume-averaged kinetic energy for the (a) third-order and (b) fifth-order schemes, and (c) the kinetic energy spectrum at t = 10. Three observations from the third-order results are worth noting. First, the instability at ηa = 0.5 is driven entirely by the acoustic waves: the entropy wave in WA-3 is always re… view at source ↗
Figure 5
Figure 5. Figure 5: Viscous Taylor-Green vortex (Re = 1600, Ma = 0.1, Sec. 7.1.1): time evolution of the volume-averaged dissipation rate ϵ on 643 and 963 grids for WA-3 scheme. DNS reference from [42]. Figures 6a and 6b show the fifth-order results. On both grids, WA-5 and WA-CR overlap in plotting accuracy and closely follow the DNS profile through the dissipation peak at t ≈ 9. The linear U-7 scheme is slightly more dissip… view at source ↗
Figure 6
Figure 6. Figure 6: Viscous Taylor-Green vortex (Re = 1600, Ma = 0.1, Sec. 7.1.1): time evolution of the volume-averaged dissipation rate ϵ on 643 and 963 grids for WA-5 scheme. DNS reference from [42]. Figure 7a shows a broader comparison on the 643 grid including WA-WENO-CR, Feng et al. [20] (TENO5DV), and ALDM (Adaptive Local Deconvolution Method) [43]. WA-5, WA-CR, and WA-WENO￾CR produce nearly identical results and outpe… view at source ↗
Figure 7
Figure 7. Figure 7: Viscous Taylor-Green vortex (Re = 1600, Ma = 0.1, 643 grid, Sec. 7.1.1): time evolution of volume-averaged kinetic energy dissipation rate ϵ(t). WA-3 and WA-5 use the optimized bias η ∗ a. WA-CR and WA-WENO-CR overlap WA-5 to plotting accuracy, confirming that the rank-1 correction introduces no accuracy penalty in smooth flows. DNS reference from Brachet et al. [42]. Figure 7c is reproduced with permissio… view at source ↗
Figure 8
Figure 8. Figure 8: Supersonic viscous Taylor-Green vortex ( [PITH_FULL_IMAGE:figures/full_fig_p023_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: presents the reference z-vorticity and compares three baseline variants on the coarse grid to clarify the specific behaviors of each scheme. The z-vorticity field is calculated using spectral derivatives via the fast Fourier transform, eliminating numerical differentiation error from the vorticity calculation regardless of the flow solver’s order. When using the KEP scheme, which conserves energy exactly, … view at source ↗
Figure 10
Figure 10. Figure 10: Double periodic shear layer (θ = 80, inviscid, t = 1, 320 × 320 grid, Sec. 7.2): z-vorticity contours for the optimized schemes. WA-3 at η ∗ a = 0.54 and WA-5 at η ∗ a = 0.6010 both reproduce two clean primary vortices. WA-CR and WA-WENO￾CR match WA-5. TENO5 produces spurious braid vortices on this grid. The WA-KEP scheme, formulated in Section 6, is tested here on the double periodic shear layer. The unm… view at source ↗
Figure 11
Figure 11. Figure 11: Double periodic shear layer (θ = 80, inviscid, t = 1, 320 × 320 grid, Sec. 7.2): z-vorticity contours for the WA-KEP approach and that of Feng et al. [20]. The value ηa = 0.56 is slightly above the WA-3 threshold of η ∗ a = 0.54. This is expected: WA-3 applies either a third-order upwind or fourth-order central scheme to the flow variables, whereas KEP is purely second-order central apart from the normal … view at source ↗
Figure 12
Figure 12. Figure 12: Rayleigh-Taylor instability (t = 1.95, Sec. 7.3): density contours at 128 × 512 (a,b) and 512 × 2048 (c,d,e). WA-CR matches WA-5 at both resolutions. WA-WENO-CR matches WA-CR on the fine grid. On the coarse 128 × 512 grid, WA-5 and WA-CR produce nearly identical density contours. The mush￾room cap structure, rolled-up vortex sheets, and interface symmetry are preserved equally well. This result confirms t… view at source ↗
Figure 13
Figure 13. Figure 13: Rayleigh-Taylor instability (t = 1.95, Sec. 7.3): density contours. Left: reference results from Fleischmann et al. [47] on the same configuration. Right: WA-CR at 1024 × 4096 resolution. The present scheme resolves finer secondary instabilities and filaments along the stem while remaining free of spurious density oscillations at the interface. 7.4. Explosion problem In this example, the initial condition… view at source ↗
Figure 14
Figure 14. Figure 14: Explosion problem (400 × 400 grid, t = 0.25, Sec. 7.4): density contours with the cross-sectional density profile along y = 0 shown in blue. The shock and contact discontinuity are captured cleanly by all schemes. 7.5. 2-D shock-entropy wave interaction The two-dimensional shock-entropy wave interaction of Acker et al. [49] tests the ability of a scheme to resolve fine-scale entropy waves. This occurs in … view at source ↗
Figure 15
Figure 15. Figure 15: 2-D shock-entropy wave interaction (t = 1.8, 400 × 80 grid, Sec. 7.5): density contours and local density profile along y = 0 for WA-CR. The exact solution is shown for reference in (b). 7.6. Double Mach reflection The double Mach reflection [51] involves a Mach 10 shock impinging on a 30◦ wedge, testing both shock￾capturing fidelity and the resolution of slip-layer vortices formed behind the Mach stem. F… view at source ↗
Figure 16
Figure 16. Figure 16: Double Mach reflection (Ma = 10, 768 × 256 grid, t = 0.3, Sec. 7.6): density contours zoomed into the Mach stem region. WA-5, WA-CR, and WA-WENO-CR all resolve the slip-layer roll-up structures. WA-3 resolves fewer vortices due to its higher background dissipation. 7.7. Two-dimensional Riemann problem The 2-D Riemann problem of configuration 3 [52] initiates four shocks at the quadrant boundaries and deve… view at source ↗
Figure 17
Figure 17. Figure 17: Two-dimensional Riemann problem, configuration 3 ( [PITH_FULL_IMAGE:figures/full_fig_p033_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Shock-bubble interaction (Ma = 6 shock, helium bubble, 400 × 400 grid, Sec. 7.8): density contours at the final time. WA-WENO-CR confirms the limiter-agnostic property of the rank-1 correction. This configuration tests the algorithm for both reconstruction paths simultaneously: the Ducros sensor activates in the shock region, where the full characteristic path is used, while the rank-1 entropy wave correc… view at source ↗
Figure 19
Figure 19. Figure 19: Viscous shock tube (Re = 1000, 1280 × 640 grid, t = 1, Sec. 7.9): flow-field density contours and density profile along the wall (y = 0). WA-CR matches WA-5. WA-3 captures the primary shock structure and the lambda-shock pattern correctly, but provides less fine-scale vortical detail in the separation region compared to the fifth-order schemes. WA-5 and WA￾CR yield nearly identical contours, both offering… view at source ↗
Figure 20
Figure 20. Figure 20: shows results at the higher Reynolds number Re = 2500. Under these conditions, the separation region contains finer vortical structures that are more sensitive to numerical dissipation. WA-3 captures the primary lambda shock and the broad separation bubble correctly, but it resolves the fine-scale Kelvin￾Helmholtz roll-up in the shear layer near x = 0.8–1.0 less clearly than WA-5 or WA-CR. Both WA-5 and W… view at source ↗
Figure 21
Figure 21. Figure 21: (a) Ducros sensor, x-direction (b) Ducros sensor, y-direction (c) Fine-grid reference, WA-CR [PITH_FULL_IMAGE:figures/full_fig_p037_21.png] view at source ↗
read the original abstract

The wave-appropriate reconstruction approach decomposes the reconstruction procedure into characteristic wave families, centralizing non-acoustic waves to minimize dissipation while retaining an upwind bias for acoustic waves. In previous implementations, the acoustic upwind parameter $\eta_a$ was fixed at its maximum value of $1.0$; however, this choice is conservative and motivated a systematic search for the minimum value that is robust across flow regimes. To this end, the CFD solver is treated as a black box within a bounded scalar minimization, which minimizes an accuracy objective for the subsonic inviscid TGV subject to a stability constraint enforced by the supersonic viscous TGV. Because the wave-appropriate framework leaves $\eta_a$ as the sole degree of freedom, the optimization converges in approximately 25 evaluations. The resulting optimal values generalize without retuning across a wide range from subsonic turbulence to hypersonic flows with shocks and contact discontinuities. The second contribution focuses on eliminating the need for an explicit contact-discontinuity detector, which is commonly required in flows involving both shock waves and contact discontinuities. In such cases, the reconstruction deficiency appears solely within the entropy characteristic wave and can be corrected by a rank-1 update along the entropy right eigenvector. The proposed algorithm relies only on the Ducros sensor and is limiter-agnostic, facilitating direct use in other schemes, such as WENO. This approach reduces wall time by $29$--$41\%$ compared to full characteristic decomposition. To further demonstrate the method's generality, introducing a controlled acoustic bias exclusively to the normal momentum in a KEP scheme eliminates spurious vortices in periodic shear layers, confirming that the acoustic stability mechanism operates independently of the discretization framework.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript introduces a wave-appropriate reconstruction for compressible flows by decomposing into characteristic waves, optimizing the acoustic parameter η_a via bounded minimization on two Taylor-Green vortex (TGV) cases (accuracy on subsonic inviscid, stability on supersonic viscous), and proposing a rank-1 entropy wave correction using only the Ducros sensor to handle contacts without explicit detectors. It claims the optimized η_a generalizes without retuning from subsonic turbulence to hypersonic shocks, reducing wall time by 29-41%, and demonstrates acoustic bias in KEP schemes.

Significance. If the generalization of the optimized η_a holds, the method offers a practical way to reduce dissipation on non-acoustic waves and computational cost in multi-regime compressible simulations, with the rank-1 correction simplifying implementation across schemes like WENO. The black-box optimization approach is efficient, converging in ~25 evaluations.

major comments (2)
  1. [Results / generalization claim] The central claim that the optimal η_a obtained from the two TGV cases generalizes without retuning to hypersonic flows with strong shocks and contacts is load-bearing but unsupported by quantitative evidence; the abstract reports time savings and the optimization procedure but provides no error tables, L2 norms, or stability metrics for hypersonic test cases, leaving open whether the acoustic dissipation requirement grows with shock strength as the skeptic analysis indicates.
  2. [Optimization procedure] The bounded minimization enforces stability only via the supersonic viscous TGV; § on optimization does not report the achieved stability margin, the objective function values at convergence, or a sensitivity study showing that the minimal η_a remains sufficient when acoustic amplitudes increase in untested hypersonic regimes.
minor comments (2)
  1. [Abstract] The 29-41% wall-time reduction is stated relative to full characteristic decomposition but the specific test cases, grid sizes, and timing methodology are not detailed in the abstract or results summary.
  2. [Methods / rank-1 correction] Clarify the precise form of the rank-1 entropy correction (right eigenvector update) and its interaction with the Ducros sensor when limiters are present; this would aid reproducibility in other schemes such as WENO.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback on our manuscript. We address each major comment below and will revise the paper to strengthen the quantitative support for our claims where appropriate.

read point-by-point responses

  1. Referee: [Results / generalization claim] The central claim that the optimal η_a obtained from the two TGV cases generalizes without retuning to hypersonic flows with strong shocks and contacts is load-bearing but unsupported by quantitative evidence; the abstract reports time savings and the optimization procedure but provides no error tables, L2 norms, or stability metrics for hypersonic test cases, leaving open whether the acoustic dissipation requirement grows with shock strength as the skeptic analysis indicates.

    Authors: We agree that explicit quantitative metrics for the hypersonic cases would strengthen the generalization claim. While the results section demonstrates stable and accurate performance on hypersonic flows with shocks and contacts (including reduced dissipation on non-acoustic waves), we did not include tabulated L2 norms or direct comparisons of acoustic dissipation growth. We will add a new table in the results section reporting L2 errors, stability indicators, and wall-time savings for the hypersonic test cases, along with a brief discussion addressing whether acoustic bias requirements increase with shock strength. revision: yes

  2. Referee: [Optimization procedure] The bounded minimization enforces stability only via the supersonic viscous TGV; § on optimization does not report the achieved stability margin, the objective function values at convergence, or a sensitivity study showing that the minimal η_a remains sufficient when acoustic amplitudes increase in untested hypersonic regimes.

    Authors: We acknowledge that the optimization section focuses on the procedure and convergence in ~25 evaluations but omits explicit reporting of the stability margin, final objective values, and sensitivity to higher acoustic amplitudes. We will revise the optimization section to include the converged objective function value, the achieved stability margin from the supersonic TGV constraint, and a short sensitivity study (e.g., perturbing acoustic amplitudes by factors of 2–5) confirming that the minimal η_a remains robust for hypersonic regimes. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation remains self-contained

full rationale

The optimization of η_a is performed explicitly via bounded scalar minimization on two TGV cases (subsonic inviscid for accuracy, supersonic viscous for stability constraint), after which the resulting scalar is applied to other regimes. This constitutes an empirical generalization claim rather than a mathematical reduction in which the target performance is forced by construction from the fitted value. The rank-1 entropy correction is obtained directly from the characteristic decomposition of the Euler equations (entropy right eigenvector) and the Ducros sensor; it does not rely on any fitted parameter or self-citation chain. No load-bearing step equates a derived prediction to its own input by definition, and the framework treats the CFD solver as a black box without smuggling ansatzes or uniqueness theorems from prior author work. The central claims are therefore independent of the optimization inputs and can be falsified by external test cases.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claims rest on the assumption that the two TGV configurations adequately sample the stability-accuracy trade-off for all compressible regimes and that the Ducros sensor plus rank-1 update fully compensates for contact discontinuities without side effects.

free parameters (1)
  • η_a
    Acoustic upwind bias parameter whose value is obtained by black-box minimization on TGV test cases.
axioms (1)
  • domain assumption The CFD solver can be treated as a black-box function for scalar minimization subject to stability constraints.
    Invoked to justify the 25-evaluation optimization procedure.

pith-pipeline@v0.9.0 · 5611 in / 1211 out tokens · 50891 ms · 2026-05-13T18:25:53.420254+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. On Reconstructing Conservative and Primitive Variables: An Eigenvector Analysis on Curvilinear Grids

    physics.comp-ph 2026-04 conditional novelty 7.0

    Conservative curvilinear eigenvectors isolate contact discontinuities in a single metric-independent entropy direction while primitive ones mix metric-dependent density and velocity terms, enabling exact rank-one corr...

  2. Wave-Appropriate Reconstruction of Compressible Multiphase and Multicomponent Flows: Fully Conservative and Semi-Conservative Eigenstructures

    physics.flu-dyn 2026-04 unverdicted novelty 7.0

    Derives explicit left and right eigenvectors for the Allaire model that enforce pressure-velocity equilibrium at interfaces via a thermodynamic jump term in the conservative form and a structural zero in the semi-cons...

Reference graph

Works this paper leans on

56 extracted references · 56 canonical work pages · cited by 2 Pith papers

  1. [1]

    A. S. Chamarthi, N. Hoffmann, S. Frankel, A wave appropriate discontinuity sensor approach for compressible flows, Physics of Fluids 35 (6) (2023)

    work page 2023

  2. [2]

    Hoffmann, A

    N. Hoffmann, A. S. Chamarthi, S. H. Frankel, Centralized gradient-based reconstruction for wall mod- elled large eddy simulations of hypersonic boundary layer transition, Journal of Computational Physics (2024) 113128

    work page 2024

  3. [3]

    A. S. Chamarthi, Wave-appropriate multidimensional upwinding approach for compressible multiphase flows, Journal of Computational Physics 538 (2025) 114157

    work page 2025

  4. [4]

    A. S. Chamarthi, A generalized adaptive central-upwind scheme for compressible flow simulations and preventing spurious vortices, arXiv preprint arXiv:2409.02340 (2024)

    work page arXiv 2024

  5. [5]

    A. S. Chamarthi, Physics appropriate interface capturing reconstruction approach for viscous compress- ible multicomponent flows, Computers & Fluids 303 (2025) 106858

    work page 2025

  6. [6]

    Ducros, V

    F. Ducros, V. Ferrand, F. Nicoud, C. Weber, D. Darracq, C. Gacherieu, T. Poinsot, Large-eddy simu- lation of the shock/turbulence interaction, Journal of Computational Physics 152 (2) (1999) 517–549

    work page 1999

  7. [7]

    Sandham, E

    N. Sandham, E. Schuelein, A. Wagner, S. Willems, J. Steelant, Transitional shock-wave/boundary-layer interactions in hypersonic flow, Journal of Fluid Mechanics 752 (2014) 1–33

    work page 2014

  8. [8]

    G. K. Batchelor, An introduction to fluid dynamics, Cambridge university press, 1967

    work page 1967

  9. [9]

    J. C. Meng, T. Colonius, Numerical simulation of the aerobreakup of a water droplet, Journal of Fluid Mechanics 835 (2018) 1108–1135

    work page 2018

  10. [10]

    De Vanna, F

    F. De Vanna, F. Picano, E. Benini, A sharp-interface immersed boundary method for moving objects in compressible viscous flows, Computers & Fluids 201 (2020) 104415

    work page 2020

  11. [11]

    P. K. Subbareddy, G. V. Candler, A fully discrete, kinetic energy consistent finite-volume scheme for compressible flows, Journal of Computational Physics 228 (5) (2009) 1347–1364. 39

    work page 2009

  12. [12]

    van Noordt, S

    W. van Noordt, S. Ganju, C. Brehm, An immersed boundary method for wall-modeled large-eddy simulation of turbulent high-mach-number flows, Journal of Computational Physics 470 (2022) 111583

    work page 2022

  13. [13]

    De Vanna, F

    F. De Vanna, F. Avanzi, M. Cogo, S. Sandrin, M. Bettencourt, F. Picano, E. Benini, Uranos: A gpu ac- celerated navier-stokes solver for compressible wall-bounded flows, Computer Physics Communications 287 (2023) 108717

    work page 2023

  14. [14]

    Ghate, S

    A. Ghate, S. K. Lele, Finite difference methods for turbulence simulations, in: Numerical Methods in Turbulence Simulation, Elsevier, 2023, pp. 235–284

    work page 2023

  15. [15]

    Sciacovelli, D

    L. Sciacovelli, D. Passiatore, P. Cinnella, G. Pascazio, Assessment of a high-order shock-capturing central-difference scheme for hypersonic turbulent flow simulations, Computers & Fluids 230 (2021) 105134

    work page 2021

  16. [16]

    F. S. Schranner, X. Y. Hu, N. A. Adams, A physically consistent weakly compressible high-resolution approach to underresolved simulations of incompressible flows, Computers & Fluids 86 (2013) 109–124

    work page 2013

  17. [17]

    F. S. Schranner, V. Rozov, N. A. Adams, Optimization of an implicit large-eddy simulation method for underresolved incompressible flow simulations, AIAA Journal 54 (5) (2016) 1567–1577

    work page 2016

  18. [18]

    J. M. Winter, F. S. Schranner, N. A. Adams, Iterative bayesian optimization of an implicit les method for under-resolved simulations of incompressible flows, in: en. In: 10th International Symposium on Turbulence and Shear Flow Phenomena, TSFP, Vol. 10, 2016

    work page 2016

  19. [19]

    Y. Feng, F. S. Schranner, J. Winter, N. A. Adams, A multi-objective bayesian optimization environment for systematic design of numerical schemes for compressible flow, Journal of Computational Physics 468 (2022) 111477

    work page 2022

  20. [20]

    Y. Feng, J. Winter, N. A. Adams, F. S. Schranner, A general multi-objective bayesian optimization framework for the design of hybrid schemes towards adaptive complex flow simulations, Journal of Computational Physics 510 (2024) 113088

    work page 2024

  21. [21]

    Y.Feng, F.S.Schranner, J.Winter, N.A.Adams, Adeepreinforcementlearningframeworkfordynamic optimization of numerical schemes for compressible flow simulations, Journal of Computational Physics 493 (2023) 112436

    work page 2023

  22. [22]

    D. A. Bezgin, A. B. Buhendwa, S. J. Schmidt, N. A. Adams, Ml-iles: End-to-end optimization of data-driven high-order godunov-type finite-volume schemes for compressible homogeneous isotropic turbulence, Journal of Computational Physics 522 (2025) 113560

    work page 2025

  23. [23]

    B. Van Leer, Upwind and high-resolution methods for compressible flow: From donor cell to residual- distribution schemes, in: 16th AIAA Computational Fluid Dynamics Conference, 2003, p. 3559

    work page 2003

  24. [24]

    P. L. Roe, Discrete models for the numerical analysis of time-dependent multidimensional gas dynamics, Journal of Computational Physics 63 (2) (1986) 458–476

    work page 1986

  25. [25]

    P. Chandrashekar, Kinetic energy preserving and entropy stable finite volume schemes for compressible euler and navier-stokes equations, Communications in Computational Physics 14 (5) (2013) 1252–1286

    work page 2013

  26. [26]

    Nishikawa, Beyond Interface Gradient: A General Principle for Constructing Diffusion Schemes, 40th Fluid Dynamics Conference and Exhibit (2010)

    H. Nishikawa, Beyond Interface Gradient: A General Principle for Constructing Diffusion Schemes, 40th Fluid Dynamics Conference and Exhibit (2010)

    work page 2010

  27. [27]

    E. F. Toro, M. Spruce, W. Speares, Restoration of the contact surface in the hll-riemann solver, Shock waves 4 (1) (1994) 25–34

    work page 1994

  28. [28]

    Harten, P

    A. Harten, P. D. Lax, B. v. Leer, On upstream differencing and godunov-type schemes for hyperbolic conservation laws, SIAM review 25 (1) (1983) 35–61

    work page 1983

  29. [29]

    A. S. Chamarthi, Gradient based reconstruction: Inviscid and viscous flux discretizations, shock cap- turing, and its application to single and multicomponent flows, Computers & Fluids 250 (2023) 105706. 40

    work page 2023

  30. [30]

    A. S. Chamarthi, Efficient high-order gradient-based reconstruction for compressible flows, Journal of Computational Physics 486 (2023) 112119

    work page 2023

  31. [31]

    A. S. Chamarthi, N. Hoffmann, H. Nishikawa, S. H. Frankel, Implicit gradients based conservative numerical scheme for compressible flows, Journal of Scientific Computing 95 (1) (2023) 17

    work page 2023

  32. [32]

    C.-W. Shu, Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, in: Advanced Numerical Approximation of Nonlinear Hyperbolic Equations: Lec- tures given at the 2nd Session of the Centro Internazionale Matematico Estivo (CIME) held in Cetraro, Italy, June 23–28, 1997, Springer, 2006, pp. 325–432

    work page 1997

  33. [33]

    S.K.Godunov, Adifferencemethodfornumericalcalculationofdiscontinuoussolutionsoftheequations of hydrodynamics, Matematicheskii Sbornik 89 (3) (1959) 271–306

    work page 1959

  34. [34]

    Van Leer, Towards the ultimate conservative difference scheme

    B. Van Leer, Towards the ultimate conservative difference scheme. iv. a new approach to numerical convection, Journal of Computational Physics 23 (3) (1977) 276–299

    work page 1977

  35. [35]

    Van Leer, Towards the ultimate conservative difference scheme

    B. Van Leer, Towards the ultimate conservative difference scheme. v. a second-order sequel to godunov’s method, Journal of computational Physics 32 (1) (1979) 101–136

    work page 1979

  36. [36]

    Suresh, H

    A. Suresh, H. Huynh, Accurate monotonicity-preserving schemes with runge-kutta time stepping, Jour- nal of Computational Physics 136 (1) (1997) 83–99

    work page 1997

  37. [37]

    Pirozzoli, On the spectral properties of shock-capturing schemes, Journal of Computational Physics 219 (2) (2006) 489–497

    S. Pirozzoli, On the spectral properties of shock-capturing schemes, Journal of Computational Physics 219 (2) (2006) 489–497

    work page 2006

  38. [38]

    D. J. Lusher, N. D. Sandham, Assessment of low-dissipative shock-capturing schemes for the compress- ible taylor–green vortex, AIAA Journal 59 (2) (2021) 533–545

    work page 2021

  39. [39]

    Virtanen, R

    P. Virtanen, R. Gommers, T. E. Oliphant, M. Haberland, T. Reddy, D. Cournapeau, E. Burovski, P. Peterson, W. Weckesser, J. Bright, et al., Scipy 1.0: fundamental algorithms for scientific computing in python, Nature methods 17 (3) (2020) 261–272

    work page 2020

  40. [40]

    Borges, M

    R. Borges, M. Carmona, B. Costa, W. S. Don, An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws, Journal of Computational Physics 227 (6) (2008) 3191–3211

    work page 2008

  41. [41]

    Jiang, C.-W

    G.-S. Jiang, C.-W. Shu, Efficient Implementation of Weighted ENO Schemes, Journal of Computational Physics 126 (126) (1995) 202–228

    work page 1995

  42. [42]

    M. E. Brachet, D. I. Meiron, S. A. Orszag, B. Nickel, R. H. Morf, U. Frisch, Small-scale structure of the taylor–green vortex, Journal of Fluid Mechanics 130 (1983) 411–452

    work page 1983

  43. [43]

    Hickel, N

    S. Hickel, N. A. Adams, J. A. Domaradzki, An adaptive local deconvolution method for implicit les, Journal of Computational Physics 213 (1) (2006) 413–436

    work page 2006

  44. [44]

    T.Liang, L.Fu, Anewhigh-ordershock-capturingtenoschemecombinedwithskew-symmetric-splitting method for compressible gas dynamics and turbulence simulation, Computer Physics Communications 302 (2024) 109236

    work page 2024

  45. [45]

    F. F. Grinstein, L. G. Margolin, W. J. Rider, Implicit large eddy simulation, Vol. 10, Cambridge university press Cambridge, 2007

    work page 2007

  46. [46]

    M. L. Minion, D. L. Brown, Performance of under-resolved two-dimensional incompressible flow simu- lations, ii, Journal of Computational Physics 138 (2) (1997) 734–765

    work page 1997

  47. [47]

    Fleischmann, S

    N. Fleischmann, S. Adami, N. A. Adams, Numerical symmetry-preserving techniques for low-dissipation shock-capturing schemes, Computers & Fluids 189 (2019) 94–107

    work page 2019

  48. [48]

    Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction, Springer Berlin Heidelberg, 2009

    E. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction, Springer Berlin Heidelberg, 2009. 41

    work page 2009

  49. [49]

    Acker, R

    F. Acker, R. d. R. Borges, B. Costa, An improved weno-z scheme, Journal of Computational Physics 313 (2016) 726–753

    work page 2016

  50. [50]

    X. Deng, Y. Shimizu, F. Xiao, A fifth-order shock capturing scheme with two-stage boundary variation diminishing algorithm, Journal of Computational Physics 386 (2019) 323–349

    work page 2019

  51. [51]

    Woodward, P

    P. Woodward, P. Colella, The numerical simulation of two-dimensional fluid flow with strong shocks, Journal of Computational Physics 54 (1) (1984) 115–173

    work page 1984

  52. [52]

    C. W. Schulz-Rinne, J. P. Collins, H. M. Glaz, Numerical solution of the riemann problem for two- dimensional gas dynamics, SIAM Journal on Scientific Computing 14 (6) (1993) 1394–1414

    work page 1993

  53. [53]

    V. Daru, C. Tenaud, Numerical simulation of the viscous shock tube problem by using a high resolution monotonicity-preserving scheme, Computers & Fluids 38 (3) (2009) 664–676

    work page 2009

  54. [54]

    A. S. Chamarthi, S. Bokor, S. H. Frankel, On the importance of high-frequency damping in high- order conservative finite-difference schemes for viscous fluxes, Journal of Computational Physics (2022) 111195

    work page 2022

  55. [55]

    Kundu, M

    A. Kundu, M. Thangadurai, G. Biswas, Investigation on shear layer instabilities and generation of vortices during shock wave and boundary layer interaction, Computers & Fluids 224 (2021) 104966

    work page 2021

  56. [56]

    M.Macklin, Warp: Ahigh-performancepythonframeworkforgpusimulationandgraphics, in: NVIDIA GPU Technology Conference (GTC), Vol. 3, 2022. 42

    work page 2022