pith. sign in

arxiv: 2606.18026 · v1 · pith:AJJOX4DCnew · submitted 2026-06-16 · 🧮 math.NA · cs.NA· physics.comp-ph

A Fourth-order Conservative Adaptive Multiresolution Wavelet Upwind Scheme for Compressible Flows

Pith reviewed 2026-06-26 23:56 UTC · model grok-4.3

classification 🧮 math.NA cs.NAphysics.comp-ph
keywords adaptive multiresolutionwavelet schemeconservative finite volumeupwind discretizationshock capturingcompressible flowsfourth-order accuracy
0
0 comments X

The pith

The scheme achieves fourth-order accuracy and machine-precision conservation for adaptive compressible flow simulations by operating entirely on cell averages.

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

The paper introduces a fourth-order conservative adaptive multiresolution wavelet upwind scheme for solving hyperbolic conservation laws that govern compressible flows. Asymmetric average-interpolating wavelets supply upwind properties for the finite-volume discretization, while symmetric wavelets handle the multiresolution decomposition that drives mesh adaptation. Because both operations use cell-average quantities, strict conservation holds during time evolution and during dynamic redistribution of cells. Interface values for flux computation are obtained directly from the wavelet reconstruction, eliminating separate ghost-cell procedures at coarse-fine boundaries. Boundary variation diminishing reconstruction at the finest level supplies non-oscillatory shock capture.

Core claim

By constructing a family of asymmetric average-interpolating wavelets that possess upwind bias and by employing symmetric counterparts for adaptation, the method performs both conservative finite-volume updates and adaptive mesh refinement on cell averages alone, thereby preserving conservation to machine precision while attaining fourth-order accuracy in smooth regions and sharp, oscillation-free capture of discontinuities.

What carries the argument

Asymmetric average-interpolating wavelets with upwind properties that reconstruct interface values directly for numerical flux evaluation inside a cell-average finite-volume framework.

If this is right

  • The scheme attains the design fourth-order accuracy on smooth problems.
  • Conservation errors remain near machine precision throughout both evolution and adaptation.
  • Numerical errors stay controlled near the user-specified threshold.
  • Shock waves and contact discontinuities are captured sharply without spurious oscillations.
  • Multiscale smooth structures are resolved with a sparse adaptive representation.

Where Pith is reading between the lines

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

  • The direct wavelet reconstruction of interface values may reduce implementation complexity relative to standard adaptive-mesh-refinement codes that require special interface bookkeeping.
  • Similar wavelet constructions could be explored for other hyperbolic systems or for extension to three space dimensions.
  • The error-threshold control built into the adaptation criterion suggests the method could be applied to long-time integrations where accumulated conservation drift must be avoided.

Load-bearing premise

Such a family of asymmetric average-interpolating wavelets with the required upwind properties can be built so that both discretization and adaptation remain strictly on cell averages.

What would settle it

A convergence study on a smooth isentropic vortex or similar problem that yields an observed order of accuracy below four, or a long adaptive run whose conservation error grows beyond round-off level, would refute the central claim.

Figures

Figures reproduced from arXiv: 2606.18026 by Bing Yang, Feng Xiao, Jizeng Wang, Xiaojing Liu, Youhe Zhou.

Figure 1
Figure 1. Figure 1: Stencils for symmetric and asymmetric wavelets. Eq. (20) to determine the coefficients of the polynomial 𝑝𝑙 (𝑥) and subsequently compute ℎ2𝑙 via integration over the half cell. However, this approach requires solving a system of linear equations and is therefore computationally implicit. Here, we propose an explicit formulation that directly evaluates ℎ2𝑙 using the primitive function of the polynomial 𝑝𝑙 (… view at source ↗
Figure 2
Figure 2. Figure 2: Scaling and wavelet functions of the 𝑁 = 4, 𝑏 = ±1 average-interpolating wavelets. When using wavelet approximation at a single resolution level, the corresponding cells are distributed uniformly. Since the average-interpolating wavelet is constructed in the dyadic framework, the mesh size only depends on the resolution level. Hence, we utilize asymmetric average-interpolating wavelet approximation at the … view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of the spectral properties of different schemes. 2.3. Conservative adaptive wavelet upwind schemes on multiresolution non-uniform cells 2.3.1. Multiresolution cell adaptation First, we introduce the cell structure of the average-interpolating wavelet multiresolution analysis, which corre￾sponds to a graded tree structure in the dyadic framework [29, 51]. It should be noted that the connectivity … view at source ↗
Figure 4
Figure 4. Figure 4: Schematic illustration of the graded-tree structure used in the adaptive average-interpolating wavelet method. As presented in subsection. 2.1.3, the physical variable 𝑢 defined on a domain Ω can be decomposed into the following average-interpolating wavelet multiresolution approximation:  𝐽max 𝐽0 𝑢(𝑥) = ∑ 𝑘∈𝐽0 ̄𝑢𝐽0 ,𝑘𝜙𝐽0 ,𝑘(𝑥) + 𝐽max ∑ −1 𝐽=𝐽0 ∑ 𝑚∈Λ𝐽 𝑑𝐽,𝑚𝜓𝐽,𝑚(𝑥), 𝐽0 = { 𝑘 ∈ ℤ ∶ supp 𝜙𝐽0 ,𝑘 ∩ Ω ≠ ∅ } , … view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of conservation errors between the adaptive AMAIW4T2 -BVD and uniform-cell AIW4T2 -BVD schemes for the sine wave advection.      0. 0 0. 5 1 . 0    0. 0 0. 2 0. 4 0. 6 0. 8 1 . 0 1 . 2 E xact AI W4T2 -BVD WE N O5-Z u x (a) 𝑡 = 2.0, 𝑁1 = 256.      0. 0 0. 5 1 . 0    0. 0 0. 2 0. 4 0. 6 0. 8 1 . 0 1 . 2 E xact AI W4T2 -BVD WE N O5-Z u x (b) 𝑡 = 500.0, 𝑁1 = 512 [PITH_FULL_IMAGE:fig… view at source ↗
Figure 6
Figure 6. Figure 6: Numerical results of the Jiang and Shu’s problem computed by the AIW4T2 -BVD and WENO5-Z schemes at 𝑡 = 2.0 and 𝑡 = 500.0. Page 29 of 54 [PITH_FULL_IMAGE:figures/full_fig_p029_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Comparison of the numerical results at 𝑡 = 2.0 of the Jiang and Shu’s problem obtained by the proposed AMAIW4T2 -BVD and the AIW4T2 -BVD schemes for different solution profiles, 𝐽0 = 6, 𝜖0 = 1.0 × 10−5 . Page 30 of 54 [PITH_FULL_IMAGE:figures/full_fig_p030_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Numerical solutions at 𝑡 = 2.0 of the Jiang and Shu’s problem obtained by the proposed adaptive AMAIW4T2 -BVD scheme with different finest resolution levels, 𝐽0 = 6, 𝜖0 = 1.0 × 10−5 . Page 31 of 54 [PITH_FULL_IMAGE:figures/full_fig_p031_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Numerical solutions at 𝑡 = 2.0 of the Jiang and Shu’s problem obtained by the proposed adaptive AMAIW4T2 -BVD scheme with different thresholding parameters, 𝐽0 = 6, and 𝐽max = 10. 0. 0 0. 5 1 . 0 1 . 5 2. 0 0 20 40 60 80 1 00  act= 71 . 4%  act(%) t J0 = 6 Jmax = 8 J0 = 6 Jmax = 1 0 J0 = 6 Jmax = 1 2  act= 30. 1 %  act= 1 1 . 2% (a) Different finest resolution levels 𝐽max, 𝜖0 = 1.0 × 10−5 . 0. 0 0. 5 1… view at source ↗
Figure 10
Figure 10. Figure 10: Active-cell percentage of the adaptive AMAIW4T2 -BVD scheme for the Jiang and Shu’s problem with different 𝐽max and 𝜖0 . Page 32 of 54 [PITH_FULL_IMAGE:figures/full_fig_p032_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Conservation errors of the AMAIW4T2 -BVD scheme for the Jiang and Shu’s problem with different 𝐽max and 𝜖0 . Page 33 of 54 [PITH_FULL_IMAGE:figures/full_fig_p033_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Numerical results at 𝑡 = 0.1∕𝜋 obtained by the proposed adaptive AMAIW4T2 -BVD scheme for the Burgers’ equation with different finest resolution levels, 𝐽0 = 4, and 𝜖0 = 1.0 × 10−5 .      0. 0 0. 5 1 . 0 1 0  1 0   1 0   1 0   1 0   1 0  1 0 0 1 0 2 E rror x  0 = 1 . 0×1 0   N1 = 1 60  0 = 1 . 0×1 0  N1 = 1 98  0 = 1 . 0×1 0  N1 = 261  0 = 1 . 0×1 0   N1 = 379 J = 4 N1 = 32 J = 1 … view at source ↗
Figure 13
Figure 13. Figure 13: Numerical errors at 𝑡 = 0.1∕𝜋 and variations in active-cell percentages with time obtained by the adaptive AMAIW4T2 -BVD scheme for the Burgers’ equation with different thresholding parameters, 𝐽0 = 4, and 𝐽max = 10. Page 36 of 54 [PITH_FULL_IMAGE:figures/full_fig_p036_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Numerical solutions at 𝑡 = 24 for different finest resolution levels and the corresponding temporal evolutions of the active-cell percentages obtained by the adaptive AMAIW4T2 -BVD scheme for the Burgers’ equation with 𝜖0 = 1.0 × 10−5 . 0 1 ×1 0 4 2×1 0 4 3×1 0 4 4×1 0 4 1 0   1 0   1 0   1 0  1 0  1 0   1 0   Econ s steps J = 6 J0 = 4 Jmax = 6 J = 8 J0 = 4 Jmax = 8 J = 1 0 J0 = 4 Jmax = 1 0 … view at source ↗
Figure 15
Figure 15. Figure 15: Conservation errors obtained by the adaptive AMAIW4T2 -BVD scheme and the corresponding uniform-cell scheme for the Burgers’ equation with different finest resolution levels up to 4.0 × 104 time steps, 𝜖0 = 1.0 × 10−5 . Page 37 of 54 [PITH_FULL_IMAGE:figures/full_fig_p037_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Comparison of the density, velocity, pressure, and zoomed density profile near the discontinuities at 𝑡 = 0.2 obtained by the AMAIW4T2 -BVD scheme (𝜖0 = 1.0 × 10−4), and the AIW4T2 -BVD scheme for Sod’s problem. Page 40 of 54 [PITH_FULL_IMAGE:figures/full_fig_p040_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Numerical solutions at 𝑡 = 0.2 of Sod’s problem obtained by the adaptive AMAIW4T2 -BVD scheme with different finest resolution levels and 𝜖0 = 1.0 × 10−4. The last panel shows a zoomed-in view near the contact discontinuity. 0. 00 0. 05 0. 1 0 0. 1 5 0. 20 0 1 0 20 30 40 50 60  act(%) t J0 = 6 Jmax = 8 J0 = 6 Jmax = 1 0 J0 = 6 Jmax = 1 2 (a) Active-cell percentage. 0. 0 0. 2 0. 4 0. 6 0. 8 1 . 0 5 6 7 8 … view at source ↗
Figure 18
Figure 18. Figure 18: Temporal evolutions of the active-cell percentages for different finest resolution levels and the active leaf cell distribution at 𝑡 = 0.2 with 𝐽max = 12 obtained by the adaptive AMAIW4T2 -BVD scheme for Sod’s problem with 𝜖0 = 1.0 × 10−4 . Page 41 of 54 [PITH_FULL_IMAGE:figures/full_fig_p041_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: Conservation errors of mass, momentum, and total energy up to 1000 time steps for Sod’s problem. The adaptive AMAIW4T2 -BVD scheme with 𝜖0 = 1.0 × 10−4, and different finest resolution levels is compared with the corresponding uniform-cell AIW4T2 -BVD scheme at 𝐽 = 𝐽max. Page 42 of 54 [PITH_FULL_IMAGE:figures/full_fig_p042_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: Numerical solutions at 𝑡 = 0.13 of the Lax’s problem obtained by the adaptive AMAIW4T2 -BVD scheme with 𝜖0 = 1.0 × 10−4 and different finest resolution levels, and WENO5-Z scheme. The right panels show the corresponding zoomed-in views near the strong discontinuities. Page 46 of 54 [PITH_FULL_IMAGE:figures/full_fig_p046_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: Active leaf cell distributions at 𝑡 = 0.13 obtained by the adaptive AMAIW4T2 -BVD scheme for the Lax’s problem with 𝜖0 = 1.0 × 10−4, and different finest resolution levels.         0 1 2 3 4 5 0. 5 1 . 0 1 . 5 2. 0 2. 5 3. 0 3. 5 4. 0 4. 5 5. 0 J0 = 3 Jm ax = 6 N1 = 31 7 E xact  x (a) 𝐽max = 6.         0 1 2 3 4 5 0. 5 1 . 0 1 . 5 2. 0 2. 5 3. 0 3. 5 4. 0 4. 5 5. 0 J0 = 3 Jm ax = 8 N1 = 4… view at source ↗
Figure 22
Figure 22. Figure 22: Numerical solutions at 𝑡 = 1.8 obtained by the adaptive AMAIW4T2 -BVD scheme for the shock–turbulence interaction problem with 𝜖0 = 1.0 × 10−3 and different finest resolution levels. Page 47 of 54 [PITH_FULL_IMAGE:figures/full_fig_p047_22.png] view at source ↗
Figure 23
Figure 23. Figure 23: Active leaf cell distributions at 𝑡 = 1.8 obtained by the adaptive AMAIW4T2 -BVD scheme for the shock–turbulence interaction problem with 𝜖0 = 1.0 × 10−3 and different finest resolution levels. Page 48 of 54 [PITH_FULL_IMAGE:figures/full_fig_p048_23.png] view at source ↗
Figure 24
Figure 24. Figure 24: Numerical solutions at 𝑡 = 0.038 of the two interacting blast waves problem obtained by the uniform-cell AIW4T2 -BVD scheme at 𝐽 = 8, 9, 10 and the adaptive AMAIW4T2 -BVD scheme with 𝜖0 = 1.0 × 10−4, and 𝐽max = 8, 9, 10. Page 49 of 54 [PITH_FULL_IMAGE:figures/full_fig_p049_24.png] view at source ↗
Figure 25
Figure 25. Figure 25: Comparison of the numerical solutions at 𝑡 = 0.038 of the two interacting blast waves problem obtained by the adaptive AMAIW4T2 -BVD scheme with 𝜖0 = 1.0 × 10−4 , and the WENO5-Z scheme. Page 50 of 54 [PITH_FULL_IMAGE:figures/full_fig_p050_25.png] view at source ↗
read the original abstract

A fourth-order conservative adaptive multiresolution average-interpolating wavelet upwind scheme is proposed for compressible flows governed by hyperbolic conservation laws. A family of asymmetric average-interpolating wavelets with upwind properties is constructed for conservative finite volume discretization, while symmetric average-interpolating wavelets are employed for multiresolution decomposition and reconstruction of physical variables in the adaptive procedure. Since both the conservative discretization and the adaptive multiresolution representation are constructed from cell-average quantities, the proposed scheme preserves strict conservation during both numerical evolution and adaptive cell redistribution. Unlike hybrid adaptive wavelet methods that use wavelets mainly for data compression and mesh adaptation, the present adaptive wavelet upwind scheme utilizes average-interpolating wavelet multiresolution approximation to reconstruct the interface values directly for numerical flux evaluation, thereby avoiding additional ghost-cell marking and reconstruction near coarse--fine mesh interfaces. The boundary variation diminishing reconstruction is incorporated at the finest resolution level to achieve non-oscillatory shock-capturing capability. Numerical tests demonstrate that the proposed scheme achieves the expected fourth-order accuracy, maintains conservation errors close to machine precision, and controls numerical errors around the prescribed threshold. The proposed method also sharply captures shock waves and contact discontinuities without spurious oscillations and resolves multiscale smooth structures through a sparse adaptive representation. These results indicate that the proposed scheme provides an efficient, conservative, and reliable approach for high-resolution simulations of compressible flows.

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

0 major / 3 minor

Summary. The manuscript proposes a fourth-order conservative adaptive multiresolution average-interpolating wavelet upwind scheme for compressible flows governed by hyperbolic conservation laws. It constructs a family of asymmetric average-interpolating wavelets with upwind bias for the finite-volume discretization and employs symmetric average-interpolating wavelets for the multiresolution decomposition and reconstruction, with both operating exclusively on cell-average quantities to ensure strict conservation during evolution and adaptive redistribution without ghost-cell procedures. Boundary variation diminishing reconstruction is incorporated at the finest level for non-oscillatory shock capturing. Numerical experiments are reported to confirm fourth-order accuracy on smooth problems, conservation errors at machine precision, and sharp capture of shocks and contacts without oscillations while using a sparse adaptive representation.

Significance. If the central claims hold, the work offers a direct integration of wavelet-based multiresolution adaptation into the flux evaluation step of a conservative finite-volume scheme, avoiding auxiliary ghost-cell marking near coarse-fine interfaces. The explicit construction of the asymmetric wavelet family (Section 3), derivation of filter coefficients enforcing cell-average reproduction and upwind bias, and proof that conservation follows from telescoping flux differencing without fitted parameters constitute clear strengths. The numerical validation in Section 5 further supports the approach for efficient high-resolution simulation of multiscale compressible flows.

minor comments (3)
  1. Abstract: the statement that the scheme 'controls numerical errors around the prescribed threshold' is imprecise; the manuscript should state explicitly what the threshold value is and how it is enforced in the adaptation criterion.
  2. Section 5: while the experiments are described as confirming the expected properties, the manuscript would benefit from a brief table summarizing the specific test cases, grid sizes, and computed error norms (L1, L2, L∞) used to verify fourth-order convergence.
  3. Notation: ensure consistent use of symbols for the adaptation threshold and the wavelet filter coefficients across Sections 3 and 4 to avoid minor ambiguity for readers.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the positive assessment leading to the recommendation to accept. There are no major comments requiring a point-by-point response.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The manuscript constructs the asymmetric average-interpolating wavelets and their filter coefficients explicitly in Section 3 from cell-average reproduction and upwind-bias requirements, then shows that both the finite-volume update and multiresolution operators act only on cell averages, so conservation follows directly from the telescoping property of flux differencing. No parameter is fitted to data and then relabeled a prediction, no uniqueness theorem is imported from prior self-work, and no ansatz is smuggled via citation. Numerical experiments in Section 5 serve only as verification, not as load-bearing inputs to the claimed fourth-order accuracy or machine-precision conservation. The central argument therefore reduces to standard finite-volume and wavelet principles without internal reduction to its own fitted quantities.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard domain assumptions for hyperbolic conservation laws and wavelet constructions, with the main addition being the specific wavelet family and its conservative integration; no invented entities are introduced.

free parameters (1)
  • adaptation threshold
    The threshold used to control numerical errors and decide mesh adaptation is a user-chosen parameter typical in adaptive methods.
axioms (2)
  • domain assumption The target problems are governed by hyperbolic conservation laws.
    Invoked in the first sentence of the abstract as the governing equations.
  • domain assumption Both conservative discretization and adaptive multiresolution representation can be built from cell-average quantities to preserve strict conservation.
    Central premise stated in the abstract for avoiding ghost-cell issues and maintaining conservation during adaptation.

pith-pipeline@v0.9.1-grok · 5795 in / 1499 out tokens · 34850 ms · 2026-06-26T23:56:26.641936+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

83 extracted references · 60 canonical work pages

  1. [1]

    C.-W.Shu,Essentiallynon-oscillatoryandweightedessentiallynon-oscillatoryschemes,ActaNumerica29(2020)701–762.doi:10.1017/ S0962492920000057

  2. [2]

    Z. J. Wang, K. Fidkowski, R. Abgrall, F. Bassi, D. Caraeni, A. Cary, H. Deconinck, R. Hartmann, K. Hillewaert, H. T. Huynh, N. Kroll, G. May,P.-O. Persson,B. VanLeer, M.Visbal, High-order CFDmethods: currentstatus andperspective, InternationalJournal for Numerical Methods in Fluids 72 (8) (2013) 811–845.doi:10.1002/fld.3767

  3. [3]

    R. J. LeVeque, Numerical methods for conservation laws, Vol. 214, Birkhäuser, Basel, 1992

  4. [4]

    J. Tang, P. C. Cui, J. Zhang, N. C. Zhou, X. J. Wu, X. Q. Gong, Y. B. Zhang, Review of mesh adaptation for fluid numerical simulation, Advances in Mechanics 53 (3) (2023) 661–692, in Chinese.doi:10.6052/1000-0992-23-013

  5. [5]

    Cockburn, C.-W

    B. Cockburn, C.-W. Shu, TVB runge-kutta local projection discontinuous galerkin finite element method for conservation laws II: General framework, Mathematics of Computation 52 (186) (1989) 411–435

  6. [6]

    S.Ii,F.Xiao,Highordermulti-momentconstrainedfinitevolumemethod.partI:Basicformulation,JournalofComputationalPhysics228(10) (2009) 3669–3707.doi:10.1016/j.jcp.2009.02.009

  7. [7]

    C.-W. Shu, S. Osher, Efficient implementation of essentially non-oscillatory shock-capturing schemes, II, Journal of computational physics 83 (1) (1989) 32–78

  8. [8]

    Merriman, Understanding the Shu-Osher conservative finite difference form, Journal of Scientific Computing 19 (1) (2003) 309–322

    B. Merriman, Understanding the Shu-Osher conservative finite difference form, Journal of Scientific Computing 19 (1) (2003) 309–322

  9. [9]

    Sebastian, C.-W

    K. Sebastian, C.-W. Shu, Multidomain WENO finite difference method with interpolation at subdomain interfaces, Journal of Scientific Computing 19 (1) (2003) 405–438.doi:10.1023/A:1025372429380. Page 51 of 54 Preprint June 17, 2026 Submitted to Journal

  10. [10]

    F. Alauzet, A parallel matrix-free conservative solution interpolation on unstructured tetrahedral meshes, Computer Methods in Applied Mechanics and Engineering 299 (2016) 116–142.doi:10.1016/j.cma.2015.10.012

  11. [11]

    J.-F.Remacle,J.E.Flaherty,M.S.Shephard,Anadaptivediscontinuousgalerkintechniquewithanorthogonalbasisappliedtocompressible flow problems, SIAM Review 45 (1) (2003) 53–72.doi:10.1137/S00361445023830

  12. [12]

    Y. Chen, G. Tóth, T. I. Gombosi, A fifth-order finite difference scheme for hyperbolic equations on block-adaptive curvilinear grids, Journal of Computational Physics 305 (2016) 604–621.doi:10.1016/j.jcp.2015.11.003

  13. [13]

    H. Li, S. Do, M. Kang, A wavelet-based adaptive WENO algorithm for Euler equations, Computers & Fluids 123 (2015) 10–22.doi: 10.1016/j.compfluid.2015.09.005

  14. [14]

    Z. Tang, H. Chen, L. Bi, X. Yuan, Numerical simulation of supersonic viscous flow based on adaptive cartesian grid, Acta Aeronautica et Astronautica Sinica 39 (5) (2018) 121697, in Chinese

  15. [15]

    Alauzet, A

    F. Alauzet, A. Loseille, A decade of progress on anisotropic mesh adaptation for computational fluid dynamics, Computer-Aided Design 72 (2016) 13–39.doi:10.1016/j.cad.2015.09.005

  16. [16]

    K.Schaal,A.Bauer,P.Chandrashekar,R.Pakmor,C.Klingenberg,V.Springel,Astrophysicalhydrodynamicswithahigh-orderdiscontinuous galerkin scheme and adaptive mesh refinement, Monthly Notices of the Royal Astronomical Society 453 (4) (2015) 4278–4300.doi: 10.1093/mnras/stv1859

  17. [17]

    Z.Zhang,C.Groth,Parallelhigh-orderanisotropicblock-basedadaptivemeshrefinementfinite-volumescheme,in:20thAIAAComputational Fluid Dynamics Conference, p. 3695

  18. [18]

    X. Qi, Z. Wang, J. Zhu, L. Tian, N. Zhao, High-order discontinuous Galerkin method with immersed boundary treatment for compressible flows on parallel adaptive cartesian grids, Physics of Fluids 36 (11) (2024) 116139.doi:10.1063/5.0238605

  19. [19]

    J. Tang, P. Cui, B. Li, Y. Zhang, H. Si, Parallel hybrid mesh adaptation by refinement and coarsening, Graphical Models 111 (2020) 101084. doi:10.1016/j.gmod.2020.101084

  20. [20]

    M. A. Park, Adjoint-based, three-dimensional error prediction and grid adaptation, AIAA Journal 42 (9) (2004) 1854–1862.doi:10.2514/ 1.10051

  21. [21]

    W.G.Habashi,J.Dompierre,Y.Bourgault,D.Ait-Ali-Yahia,M.Fortin,M.-G.Vallet,Anisotropicmeshadaptation:towardsuser-independent, mesh-independentandsolver-independentCFD.partI:generalprinciples,InternationalJournalforNumericalMethodsinFluids32(6)(2000) 725–744

  22. [22]

    S.J.Kamkar,A.M.Wissink,V.Sankaran,A.Jameson,Feature-drivencartesianadaptivemeshrefinementforvortex-dominatedflows,Journal of Computational Physics 230 (16) (2011) 6271–6298.doi:10.1016/j.jcp.2011.04.024

  23. [23]

    Roy, Strategies for driving mesh adaptation in CFD, in: 47th AIAA aerospace sciences meeting including the new horizons forum and aerospace exposition, 2009, p

    C. Roy, Strategies for driving mesh adaptation in CFD, in: 47th AIAA aerospace sciences meeting including the new horizons forum and aerospace exposition, 2009, p. 1302

  24. [24]

    K.J.Fidkowski,D.L.Darmofal,Reviewofoutput-basederrorestimationandmeshadaptationincomputationalfluiddynamics,AIAAJournal 49 (4) (2011) 673–694.doi:10.2514/1.J050073

  25. [25]

    Alauzet, O

    F. Alauzet, O. Pironneau, Continuous and discrete adjoints to the Euler equations for fluids, International Journal for Numerical Methods in Fluids 70 (2) (2012) 135–157.doi:10.1002/fld.2681

  26. [26]

    P.J.Frey,F.Alauzet,AnisotropicmeshadaptationforCFDcomputations,ComputerMethodsinAppliedMechanicsandEngineering194(48) (2005) 5068–5082.doi:10.1016/j.cma.2004.11.025

  27. [27]

    B. Yang, J. Wang, X. Liu, Y. Zhou, Y. Feng, Wavelet-based numerical methods and their applications in computational mechanics, Advances in Mechanics 54 (3) (2024) 427–476.doi:10.6052/1000-0992-24-009

  28. [28]

    A. Harten, Multiresolution algorithms for the numerical solution of hyperbolic conservation laws, Communications on Pure and Applied Mathematics 48 (12) (1995) 1305–1342.doi:10.1002/cpa.3160481201

  29. [29]

    Cohen, S

    A. Cohen, S. Kaber, S. M?ller, M. Postel, Fully adaptive multiresolution finite volume schemes for conservation laws, Mathematics of Computation 72 (241) (2003) 183–225

  30. [30]

    O. V. Vasilyev, C. Bowman, Second-generation wavelet collocation method for the solution of partial differential equations, Journal of Computational Physics 165 (2) (2000) 660–693.doi:10.1006/jcph.2000.6638

  31. [31]

    O. V. Vasilyev, S. Paolucci, A dynamically adaptive multilevel wavelet collocation method for solving partial differential equations in a finite domain, Journal of Computational Physics 125 (2) (1996) 498–512.doi:10.1006/jcph.1996.0111

  32. [32]

    Schneider, O

    K. Schneider, O. V. Vasilyev, Wavelet methods in computational fluid dynamics, Annual review of fluid mechanics 42 (2010) 473–503

  33. [33]

    Harten, Discrete multi-resolution analysis and generalized wavelets, Applied Numerical Mathematics 12 (1) (1993) 153–192.doi: 10.1016/0168-9274(93)90117-A

    A. Harten, Discrete multi-resolution analysis and generalized wavelets, Applied Numerical Mathematics 12 (1) (1993) 153–192.doi: 10.1016/0168-9274(93)90117-A

  34. [34]

    R. A. DeVore, Nonlinear approximation, Acta Numerica 7 (1998) 51–150.doi:10.1017/S0962492900002816

  35. [35]

    Deiterding, M

    R. Deiterding, M. O. Domingues, S. M. Gomes, K. Schneider, Comparison of adaptive multiresolution and adaptive mesh refinement applied to simulations of the compressible Euler equations, SIAM Journal on Scientific Computing 38 (5) (2016) S173–S193.doi: 10.1137/15m1026043

  36. [36]

    B.Yang,J.Wang,X.Liu,Y.Zhou,High-orderadaptivemultiresolutionwaveletupwindschemesforhyperbolicconservationlaws,Computers & Fluids 269 (2024) 106111.doi:10.1016/j.compfluid.2023.106111

  37. [37]

    Minbashian, H

    H. Minbashian, H. Adibi, M. Dehghan, An adaptive wavelet space-time SUPG method for hyperbolic conservation laws, Numerical Methods for Partial Differential Equations 33 (6) (2017) 2062–2089.doi:10.1002/num.22180

  38. [38]

    doi:10.1137/23M1588627

    R.M.Pereira,N.Nguyenvanyen,K.Schneider,M.Farge,Areadaptivegalerkinschemesdissipative?,SIAMReview65(4)(2023)1109–1134. doi:10.1137/23M1588627

  39. [39]

    Bertoluzza, Adaptive wavelet collocation method for the solution of Burgers equation, Transport Theory and Statistical Physics 25 (3-5) (1996) 339–352.doi:10.1080/00411459608220705

    S. Bertoluzza, Adaptive wavelet collocation method for the solution of Burgers equation, Transport Theory and Statistical Physics 25 (3-5) (1996) 339–352.doi:10.1080/00411459608220705

  40. [40]

    Page 52 of 54 Preprint June 17, 2026 Submitted to Journal

    A.Harten,P.D.Lax,B.V.Leer,OnupstreamdifferencingandGodunov-typeschemesforhyperbolicconservationlaws,SIAMReview25(1) (1983) 35–61.doi:10.1137/1025002. Page 52 of 54 Preprint June 17, 2026 Submitted to Journal

  41. [41]

    B. Yang, Y. Zhou, J. Wang, High-order second-generation wavelet upwind schemes with multiresolution self-adaptive capabilities for hyperbolic conservation laws, Extreme Mechanics Letters 70 (2024) 102192.doi:10.1016/j.eml.2024.102192

  42. [42]

    B. Yang, X. Liu, Y. Zhou, J. Wang, A conservative wavelet upwind scheme for compressible flows, Acta Mechanica Sinica 41 (7) (2025) 325178.doi:10.1007/s10409-025-25178-x

  43. [43]

    J.M.Restrepo,G.K.Leaf,Wavelet-galerkindiscretizationofhyperbolicequations,JournalofComputationalPhysics122(1)(1995)118–128

  44. [44]

    B.K.Alpert,AclassofbasesinL2forthesparserepresentationofintegraloperators,SIAMJournalonMathematicalAnalysis24(1)(1993) 246–262.doi:10.1137/0524016

  45. [45]

    Hovhannisyan, S

    N. Hovhannisyan, S. Müller, R. Schäfer, Adaptive multiresolution discontinuous galerkin schemes for conservation laws, Mathematics of computation 83 (285) (2014) 113–151

  46. [46]

    N.Gerhard,S.Müller,Adaptivemultiresolutiondiscontinuousgalerkinschemesforconservationlaws:multi-dimensionalcase,Computational and Applied Mathematics 35 (2) (2016) 321–349.doi:10.1007/s40314-014-0134-y

  47. [47]

    Gerhard, F

    N. Gerhard, F. Iacono, G. May, S. Müller, R. Schäfer, A high-order discontinuous galerkin discretization with multiwavelet-based grid adaptation for compressible flows, Journal of Scientific Computing 62 (1) (2015) 25–52.doi:10.1007/s10915-014-9846-9

  48. [48]

    Harten, Adaptive multiresolution schemes for shock computations, Journal of Computational Physics 115 (2) (1994) 319–338.doi: 10.1006/jcph.1994.1199

    A. Harten, Adaptive multiresolution schemes for shock computations, Journal of Computational Physics 115 (2) (1994) 319–338.doi: 10.1006/jcph.1994.1199

  49. [49]

    B.L.Bihari,A.Harten,Multiresolutionschemesforthenumericalsolutionof2-DconservationlawsI,SIAMJournalonScientificComputing 18 (2) (1997) 315–354.doi:10.1137/S1064827594278848

  50. [50]

    Harten, Multiresolution representation of data: A general framework, SIAM Journal on Numerical Analysis 33 (3) (1996) 1205–1256

    A. Harten, Multiresolution representation of data: A general framework, SIAM Journal on Numerical Analysis 33 (3) (1996) 1205–1256

  51. [51]

    Roussel, K

    O. Roussel, K. Schneider, A. Tsigulin, H. Bockhorn, A conservative fully adaptive multiresolution algorithm for parabolic pdes, Journal of Computational Physics 188 (2) (2003) 493–523.doi:10.1016/S0021-9991(03)00189-X

  52. [52]

    D. A. Castro, S. M. Gomes, J. Stolfi, High-order adaptive finite-volume schemes in the context of multiresolution analysis for dyadic grids, Computational and Applied Mathematics 35 (1) (2016) 1–16.doi:10.1007/s40314-014-0159-2

  53. [53]

    D. L. Donoho, Smooth wavelet decompositions with blocky coefficient kernels, in: Recent Advances in Wavelet Analysis, Academic Press, 1994

  54. [54]

    P. G. Ciarlet, Linear and Nonlinear Functional Analysis with Applications, Society for Industrial and Applied Mathematics, Philadelphia, 2013

  55. [55]

    Daubechies, Ten Lectures on Wavelets, Society for industrial and applied mathematics, Philadelphia, 1992

    I. Daubechies, Ten Lectures on Wavelets, Society for industrial and applied mathematics, Philadelphia, 1992

  56. [56]

    Sweldens, P

    W. Sweldens, P. Schroder, Building your own wavelets at home, Vol. Springer, Springer, Heidelberg, 2005

  57. [57]

    X. Liu, G. Liu, J. Wang, Y. Zhou, A wavelet multiresolution interpolation galerkin method for targeted local solution enrichment, Computational Mechanics 64 (2019) 986–1096

  58. [58]

    B. Yang, J. Wang, X. Liu, Y. Zhou, Stability and resolution analysis of the wavelet collocation upwind schemes for hyperbolic conservation laws, Fluids 8 (2) (2023) 65.doi:10.3390/fluids8020065

  59. [59]

    Vichnevetsky, J

    R. Vichnevetsky, J. B. Bowles, Fourier Analysis of Numerical Approximations of Hyperbolic Equations, SIAM, Philadelphia, 1982.doi: 10.1137/1.9781611970876

  60. [60]

    Gottlieb, C.-W

    S. Gottlieb, C.-W. Shu, E. Tadmor, Strong stability-preserving high-order time discretization methods, SIAM Review 43 (1) (2001) 89–112. doi:10.1137/S003614450036757X

  61. [61]

    Z. Sun, S. Inaba, F. Xiao, Boundary variation diminishing (BVD) reconstruction: A new approach to improve Godunov schemes, Journal of Computational Physics 322 (2016) 309–325.doi:10.1016/j.jcp.2016.06.051

  62. [62]

    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.doi:10.1016/j.jcp.2019.02.024

  63. [63]

    X. Deng, Y. Shimizu, B. Xie, F. Xiao, Constructing higher order discontinuity-capturing schemes with upwind-biased interpolations and boundary variation diminishing algorithm, Computers & Fluids 200 (2020) 104433.doi:10.1016/j.compfluid.2020.104433

  64. [64]

    Deng, Z.-H

    X. Deng, Z.-H. Jiang, P. Vincent, F. Xiao, C. Yan, A new paradigm of dissipation-adjustable, multi-scale resolving schemes for compressible flows, Journal of Computational Physics 466 (2022) 111287.doi:10.1016/j.jcp.2022.111287

  65. [65]

    Wakimura, S

    H. Wakimura, S. Takagi, F. Xiao, Symmetry-preserving enforcement of low-dissipation method based on boundary variation diminishing principle, Computers & Fluids 233 (2022) 105227.doi:10.1016/j.compfluid.2021.105227

  66. [66]

    C. Pan, S. Song, C. Chen, X. Li, X. Shen, F. Xiao, A high-fidelity finite volume scheme for ideal magnetohydrodynamics equations using boundary variation diminishing algorithm, Physics of Fluids 36 (11) (2024) 116136.doi:10.1063/5.0237231

  67. [67]

    Z.He,Y.Ruan,Y.Yu,B.Tian,F.Xiao,Self-adjustingsteepness-basedschemesthatpreservediscontinuousstructuresincompressibleflows, Journal of Computational Physics 463 (2022) 111268.doi:10.1016/j.jcp.2022.111268

  68. [68]

    P. L. Roe, Approximate Riemann solvers, parameter vectors, and difference schemes, Journal of Computational Physics 135 (2) (1997) 250– 258.doi:10.1006/jcph.1997.5705

  69. [69]

    Pirozzoli, On the spectral properties of shock-capturing schemes, Journal of Computational Physics 219 (2) (2006) 489–497.doi: 10.1016/j.jcp.2006.07.009

    S. Pirozzoli, On the spectral properties of shock-capturing schemes, Journal of Computational Physics 219 (2) (2006) 489–497.doi: 10.1016/j.jcp.2006.07.009

  70. [70]

    Liandrat, P

    J. Liandrat, P. H. Tchamitchian, Resolution of the 1D regularized burgers equation using a spatial wavelet approximation, Report, NASA Langley Research Center (01/01 1990)

  71. [71]

    Jiang, C.-W

    G.-S. Jiang, C.-W. Shu, Efficient implementation of weighted ENO schemes, Journal of Computational Physics 126 (1) (1996) 202–228. doi:10.1006/jcph.1996.0130

  72. [72]

    Holmström, Solving hyperbolic PDEs using interpolating wavelets, SIAM Journal on Scientific Computing 21 (2) (1999) 405–420

    M. Holmström, Solving hyperbolic PDEs using interpolating wavelets, SIAM Journal on Scientific Computing 21 (2) (1999) 405–420. doi:10.1137/s1064827597316278

  73. [73]

    L. Tang, S. Song, A multiresolution finite volume scheme for two-dimensional hyperbolic conservation laws, Journal of Computational and Applied Mathematics 214 (2) (2008) 583–595.doi:10.1016/j.cam.2007.03.014. Page 53 of 54 Preprint June 17, 2026 Submitted to Journal

  74. [74]

    H. T. Huynh, A Flux Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin Methods, Fluid Dynamics and Co-located Conferences, American Institute of Aeronautics and Astronautics, 2007.doi:10.2514/6.2007-4079

  75. [75]

    Abgrall, Apersonal discussion on conservation,and how to formulate it,in: E

    R. Abgrall, Apersonal discussion on conservation,and how to formulate it,in: E. Franck, J. Fuhrmann,V. Michel-Dansac, L. Navoret (Eds.), Finite Volumes for Complex Applications X–Volume 1, Elliptic and Parabolic Problems, Springer Nature Switzerland, 2023, pp. 3–19

  76. [76]

    E. F. Toro, Riemann solvers and numerical methods for fluid dynamics: a practical introduction, Springer Science & Business Media, Heidelberg, 2013

  77. [77]

    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.doi:10.1016/j.jcp.2007.11.038

  78. [78]

    W.-S. Don, R. Borges, Accuracy of the weighted essentially non-oscillatory conservative finite difference schemes, Journal of Computational Physics 250 (2013) 347–372.doi:10.1016/j.jcp.2013.05.018

  79. [79]

    Harten, B

    A. Harten, B. Engquist, S. Osher, S. R. Chakravarthy, Uniformly high order accurate essentially non-oscillatory schemes, III, Journal of Computational Physics 71 (2) (1987) 231–303.doi:10.1016/0021-9991(87)90031-3

  80. [80]

    Cockburn, M

    B. Cockburn, M. Luskin, C.-W. Shu, E. Süli, Enhanced accuracy by post-processing for finite element methods for hyperbolic equations, Mathematics of Computation 72 (242) (2003) 577–606.doi:10.1090/S0025-5718-02-01464-3

Showing first 80 references.