pith. sign in

arxiv: 2605.20583 · v3 · pith:WREYNZMXnew · submitted 2026-05-20 · 🧮 math.NA · cs.NA

Multilevel Isogeometric Projection Stabilization via Quasi-Interpolation for Advection-Dominated Problems

Zakaria El Hasnaoui , Ahmed Ratnani This is my paper

Pith reviewed 2026-05-22 09:19 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords isogeometric analysisstabilizationadvection-dominated problemsquasi-interpolationmultilevel methodsconvection-diffusion equationsB-splineserror estimates
0
0 comments X

The pith

A linear multilevel quasi-interpolation method stabilizes advection-dominated isogeometric problems with reduced parameter sensitivity.

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

The paper introduces a stabilization technique for convection-diffusion problems in isogeometric analysis that uses continuous B-spline quasi-interpolants to extract and penalize fine-scale fluctuations at multiple levels. This hierarchical approach with mesh-dependent weights avoids the tuning issues of residual-based methods and the discontinuous spaces of local projection stabilization. It provides a priori error estimates and a discrete inf-sup condition in one dimension under a stability hypothesis. Numerical results on benchmarks indicate substantial reduction of oscillations, including in pure advection, and performance on par with nonlinear schemes despite the linear nature of the formulation.

Core claim

The central discovery is a multilevel projection stabilization method that penalizes fine-scale fluctuations extracted via continuous B-spline quasi-interpolants applied hierarchically across nested discrete spaces, yielding robust control of spurious oscillations in advection-dominated regimes while maintaining a fully linear formulation and using a global parameter scaling that reduces problem-dependent tuning.

What carries the argument

Multilevel projection-based stabilization via continuous B-spline quasi-interpolants that hierarchically penalizes fine-scale fluctuations using explicit mesh-dependent weights.

Load-bearing premise

The discrete inf-sup condition holds under a numerically validated stability hypothesis ensuring robust streamline derivative control in the one-dimensional setting with constant advection.

What would settle it

Observing failure to satisfy the inf-sup condition or excessive oscillations in a one-dimensional constant advection test case with the proposed method would falsify the stability claim.

Figures

Figures reproduced from arXiv: 2605.20583 by Ahmed Ratnani, Zakaria El Hasnaoui.

Figure 1
Figure 1. Figure 1: Results for the 1D advection-reaction problem [PITH_FULL_IMAGE:figures/full_fig_p016_1.png] view at source ↗
Figure 1
Figure 1. Figure 1: Results for the 1D advection-reaction problem [PITH_FULL_IMAGE:figures/full_fig_p017_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Standard Galerkin solution for the 1D advection-diffusion problem [PITH_FULL_IMAGE:figures/full_fig_p017_2.png] view at source ↗
Figure 2
Figure 2. Figure 2: Standard Galerkin solution for the 1D advection-diffusion problem [PITH_FULL_IMAGE:figures/full_fig_p018_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of stabilized solutions for the 1D advection-diffusion problem [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Relative L 2 -error convergence plot for the 1D advection-diffusion problem (38) (Test 2) with polynomial degree p = 5 [PITH_FULL_IMAGE:figures/full_fig_p018_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Zoom near the boundary layer showing solutions at different levels for [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 5
Figure 5. Figure 5: Zoom near the boundary layer showing solutions at different levels for [PITH_FULL_IMAGE:figures/full_fig_p019_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Condition number comparison for the 1D advection-diffusion problem [PITH_FULL_IMAGE:figures/full_fig_p019_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Standard Galerkin surface plot for the parabolic boundary layers problem [PITH_FULL_IMAGE:figures/full_fig_p020_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Multilevel stabilized surface plot for the parabolic boundary layers problem [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Comparison of stabilized solutions for the parabolic boundary layers problem [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Surface plots of the numerical solutions for the two internal layers problem [PITH_FULL_IMAGE:figures/full_fig_p021_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Cross-sectional profiles at x = 0.5 for the two internal layers problem (1) (Test 4) with p = 3, ne = 64 × 64, and L = 5. (a) Comparison between Multilevel and SUPG methods. (b) Comparison between Multilevel method and the reference solution [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: SUPG surface plot for the rotational problem (1) (Test 5) with [PITH_FULL_IMAGE:figures/full_fig_p023_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Multilevel stabilized surface plot for the rotational problem (1) (Test 5) with [PITH_FULL_IMAGE:figures/full_fig_p023_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Surface plots for the 2D pure advection problem [PITH_FULL_IMAGE:figures/full_fig_p024_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Cross-sectional profiles for the 2D pure advection problem [PITH_FULL_IMAGE:figures/full_fig_p024_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Relative error convergence plots for the 2D pure advection problem [PITH_FULL_IMAGE:figures/full_fig_p025_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Condition number comparison for the pure advection problem across different polynomial degrees. [PITH_FULL_IMAGE:figures/full_fig_p025_17.png] view at source ↗
read the original abstract

This paper presents a novel multilevel projection-based stabilization method for advection-dominated convection--diffusion problems within the framework of Isogeometric Analysis. The proposed approach extracts and penalizes fine-scale fluctuations using continuous B-spline quasi-interpolants, avoiding both the highly sensitive parameters used in residual-based stabilization methods and the discontinuous auxiliary spaces required by classical Local Projection Stabilization. Stabilization is applied hierarchically across nested levels of the discrete space via explicit mesh-dependent weights. We establish the theoretical foundation of the method by deriving a priori error estimates, supplemented by a discrete inf-sup condition established for the one-dimensional setting with constant advection under a numerically validated stability hypothesis that ensures robust streamline derivative control. Numerical experiments on stringent benchmarks demonstrate the method's ability to significantly reduce spurious oscillations across a variety of regimes, including the limiting cases of pure advection and advection--reaction. Notably, despite being a fully linear formulation, the method achieves significant reduction of undershoots near sharp layers, providing an effective stabilization mechanism without requiring complex nonlinear terms. Furthermore, by utilizing a robust global parameter scaling, the proposed approach significantly alleviates the parameter sensitivity that typically affects residual-based alternatives, reducing the strong dependence on problem-dependent tuning.

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 multilevel projection-based stabilization method for advection-dominated convection-diffusion problems within Isogeometric Analysis. It employs continuous B-spline quasi-interpolants to extract and penalize fine-scale fluctuations hierarchically across nested levels using explicit mesh-dependent weights. The authors derive a priori error estimates and establish a discrete inf-sup condition under a stability hypothesis for the one-dimensional constant-advection case that is numerically validated to ensure robust streamline derivative control. Numerical experiments on benchmarks, including pure advection and advection-reaction limits, demonstrate significant reduction of spurious oscillations and undershoots near sharp layers, with performance comparable to nonlinear shock-capturing schemes and reduced sensitivity via robust global parameter scaling.

Significance. If the stability hypothesis extends to the multidimensional and variable-coefficient regimes in the benchmarks, this would represent a useful contribution by providing a fully linear stabilization approach in IGA that achieves oscillation control comparable to nonlinear methods while avoiding highly sensitive parameters and discontinuous auxiliary spaces. The emphasis on a robust global parameter scaling and the numerical results on stringent benchmarks are strengths that could support practical adoption in advection-dominated problems.

major comments (2)
  1. [Abstract and theoretical foundation section on discrete inf-sup condition] The discrete inf-sup condition and a priori error estimates are derived under a stability hypothesis that ensures robust streamline derivative control and is numerically validated only for the one-dimensional setting with constant advection (as stated in the abstract and the theoretical foundation section). This hypothesis directly supports the central theoretical claims, yet the numerical benchmarks include multidimensional problems; without additional validation or proof that the hypothesis carries over, the justification for robustness in the claimed regimes is incomplete.
  2. [Numerical experiments section] The numerical experiments section reports performance on benchmarks but does not include explicit error bars, multiple-run statistics, or details on data exclusion criteria. This weakens the strength of the claims regarding significant reduction of undershoots and comparability to nonlinear schemes, as quantitative reliability measures are needed to support the central numerical validation.
minor comments (2)
  1. [Method description] The notation for the global parameter scaling and mesh-dependent weights should be defined with explicit formulas early in the method section to enhance clarity and reproducibility.
  2. [Figures in numerical section] Figure captions for the benchmark results could be expanded to specify the exact parameter values used in the global scaling for each test case.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address the major comments point by point below, indicating the revisions we intend to make.

read point-by-point responses

  1. Referee: [Abstract and theoretical foundation section on discrete inf-sup condition] The discrete inf-sup condition and a priori error estimates are derived under a stability hypothesis that ensures robust streamline derivative control and is numerically validated only for the one-dimensional setting with constant advection (as stated in the abstract and the theoretical foundation section). This hypothesis directly supports the central theoretical claims, yet the numerical benchmarks include multidimensional problems; without additional validation or proof that the hypothesis carries over, the justification for robustness in the claimed regimes is incomplete.

    Authors: We agree that the discrete inf-sup condition and the supporting stability hypothesis are established and numerically validated only for the one-dimensional constant-advection case, as explicitly stated in the manuscript. The a priori error estimates are derived under this hypothesis. For the multidimensional and variable-coefficient benchmarks, robustness is demonstrated through the numerical results, which show effective control of oscillations and undershoots. While a complete theoretical extension to higher dimensions would strengthen the claims, the existing numerical evidence across the reported benchmarks supports the method's performance in those regimes. In the revised version, we will clarify the scope of the theoretical analysis in the relevant sections and emphasize that multidimensional robustness relies on numerical validation. revision: partial

  2. Referee: [Numerical experiments section] The numerical experiments section reports performance on benchmarks but does not include explicit error bars, multiple-run statistics, or details on data exclusion criteria. This weakens the strength of the claims regarding significant reduction of undershoots and comparability to nonlinear schemes, as quantitative reliability measures are needed to support the central numerical validation.

    Authors: The numerical experiments consist of deterministic finite-element simulations of the underlying PDEs. Consequently, there is no stochastic variability, and statistics such as error bars from multiple independent runs or data exclusion criteria are not applicable. All simulation results are reported without selective omission. We will add an explicit statement in the numerical experiments section of the revised manuscript clarifying the deterministic character of the computations to address this point. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected in derivation chain

full rationale

The paper derives a priori error estimates after establishing a discrete inf-sup condition under a numerically validated stability hypothesis for the 1D constant-advection case. This is a conditional theoretical result resting on an external numerical check rather than any self-referential reduction where a claimed prediction or first-principles result equals its own inputs by construction. No self-definitional steps, fitted inputs renamed as predictions, load-bearing self-citations, uniqueness theorems imported from the authors' prior work, or ansatzes smuggled via citation appear in the provided text. The central claims are supported by independent numerical benchmarks on stringent test cases, keeping the derivation self-contained against external validation.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claims rest on a stability hypothesis for the inf-sup condition and mesh-dependent weights plus a global scaling parameter whose selection is not fully derived from first principles.

free parameters (2)
  • global parameter scaling
    Used for robust stabilization across regimes; chosen to alleviate problem-dependent tuning but still requires selection.
  • mesh-dependent weights
    Explicit weights applied hierarchically across nested levels; depend on mesh size and are part of the stabilization definition.
axioms (1)
  • domain assumption Stability hypothesis ensuring robust streamline derivative control in 1D constant advection case
    Invoked to establish the discrete inf-sup condition and support a priori error estimates; numerically validated but not proven in general.

pith-pipeline@v0.9.0 · 5744 in / 1494 out tokens · 37022 ms · 2026-05-22T09:19:22.117830+00:00 · methodology

Review history (2 revisions) →

discussion (0)

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