pith. sign in

arxiv: 2605.20583 · v2 · 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, delivering performance comparable to complex nonlinear shock-capturing schemes. 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.

Reference graph

Works this paper leans on

40 extracted references · 40 canonical work pages

  1. [1]

    V . John, P. Knobloch, and J. Novo. Finite elements for scalar convection-dominated equations and incompressible flow problems: a never ending story?Computing and Visualization in Science, 19(5):47–63, December 2018

    work page 2018

  2. [2]

    Roos and M

    H.-G. Roos and M. Stynes. Some open questions in the numerical analysis of singularly perturbed differential equations.Computational Methods in Applied Mathematics, 15(4):531–550, 2015

    work page 2015

  3. [3]

    Harari and T.J.R

    I. Harari and T.J.R. Hughes. Stabilized finite element methods for steady advection—diffusion with production. Computer Methods in Applied Mechanics and Engineering, 115(1):165–191, 1994

    work page 1994

  4. [4]

    Hughes, G

    T.J.R. Hughes, G. Scovazzi, and L.P. Franca. Multiscale and Stabilized Methods. InEncyclopedia of Computational Mechanics. John Wiley & Sons, Ltd, Chichester, 2007

    work page 2007

  5. [5]

    Augustin, A

    M. Augustin, A. Caiazzo, A. Fiebach, J. Fuhrmann, V . John, A. Linke, and R. Umla. An assessment of discretiza- tions for convection-dominated convection–diffusion equations.Computer Methods in Applied Mechanics and Engineering, 200(47):3395–3409, 2011

    work page 2011

  6. [6]

    Hauke, G

    G. Hauke, G. Sangalli, and M.H. Doweidar. Combining adjoint stabilized methods for the advection-diffusion- reaction problem.Mathematical Models and Methods in Applied Sciences, 17(02):305–326, 2007

    work page 2007

  7. [7]

    ten Eikelder and I

    M.F.P. ten Eikelder and I. Akkerman. Correct energy evolution of stabilized formulations: The relation between vms, supg and gls via dynamic orthogonal small-scales and isogeometric analysis. i: The convective–diffusive context.Computer Methods in Applied Mechanics and Engineering, 331:259–280, 2018

    work page 2018

  8. [8]

    Braack and E

    M. Braack and E. Burman. Local projection stabilization for the oseen problem and its interpretation as a variational multiscale method.SIAM Journal on Numerical Analysis, 43(6):2544–2566, 2006

    work page 2006

  9. [9]

    Guermond

    J.-L. Guermond. Subgrid stabilization of Galerkin approximations of monotone operators.ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik, 79(S1):29–32, 1999

    work page 1999

  10. [10]

    R. Codina. On stabilized finite element methods for linear systems of convection–diffusion-reaction equations. Computer Methods in Applied Mechanics and Engineering, 188(1):61–82, July 2000

    work page 2000

  11. [11]

    Knobloch and G

    P. Knobloch and G. Lube. Local projection stabilization for advection–diffusion–reaction problems: One-level vs. two-level approach.Applied Numerical Mathematics, 59(12):2891–2907, December 2009

    work page 2009

  12. [12]

    Knobloch, P

    P. Knobloch, P. Lukáš, and P. Solin. On error indicators for optimizing parameters in stabilized methods.Advances in Computational Mathematics, 45(4):1853–1862, August 2019. 26 Stabilized Multilevel Isogeometric Quasi-Interpolant Projection Method

    work page 2019

  13. [13]

    Hughes, J.A

    T.J.R. Hughes, J.A. Cottrell, and Y . Bazilevs. Isogeometric analysis: Cad, finite elements, nurbs, exact geometry and mesh refinement.Computer Methods in Applied Mechanics and Engineering, 194(39):4135–4195, 2005

    work page 2005

  14. [14]

    Bazilevs, L

    Y . Bazilevs, L. Beirão Da Veiga, J.A. Cottrell, T.J.R. Hughes, and G. Sangalli. Isogeometric analysis: approxima- tion, stability and error estimates for h-refined meshes.Mathematical Models and Methods in Applied Sciences, 16(07):1031–1090, July 2006

    work page 2006

  15. [15]

    Manni, F

    C. Manni, F. Pelosi, and M.L. Sampoli. Isogeometric analysis in advection–diffusion problems: Tension splines approximation.Journal of Computational and Applied Mathematics, 236(4):511–528, 2011. International Workshop on Multivariate Approximation and Interpolation with Applications (MAIA 2010)

    work page 2011

  16. [16]

    John and L

    V . John and L. Schumacher. A study of isogeometric analysis for scalar convection–diffusion equations.Applied Mathematics Letters, 27:43–48, 2014

    work page 2014

  17. [17]

    Key, M.R.A

    K. Key, M.R.A. Abdelmalik, S. Elgeti, T.J.R. Hughes, and F.A. Baidoo. Finite element and isogeometric stabilized methods for the advection-diffusion-reaction equation.Computer Methods in Applied Mechanics and Engineering, 417:116354, December 2023

    work page 2023

  18. [18]

    Brooks and T.J.R

    A.N. Brooks and T.J.R. Hughes. Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations.Computer Methods in Applied Mechanics and Engineering, 32(1):199–259, September 1982

    work page 1982

  19. [19]

    Hughes, G.R

    T.J.R. Hughes, G.R. Feijóo, L. Mazzei, and J.-B. Quincy. The variational multiscale method—a paradigm for computational mechanics.Computer Methods in Applied Mechanics and Engineering, 166(1):3–24, 1998. Advances in Stabilized Methods in Computational Mechanics

    work page 1998

  20. [20]

    Matthies, P

    G. Matthies, P. Skrzypacz, and L. Tobiska. A unified convergence analysis for local projection stabilisations applied to the Oseen problem.ESAIM : Modélisation mathématique et analyse numérique, 41(4):713–742, 2007

    work page 2007

  21. [21]

    Hauke and M.H

    G. Hauke and M.H. Doweidar. Exploring new subgrid scale stabilized methods for advection-diffusion-reaction. InECCOMAS 2004 – European Congress on Computational Methods in Applied Sciences and Engineering, Jyväskylä, Finland, July 2004

    work page 2004

  22. [22]

    Lyche, C

    T. Lyche, C. Manni, and H. Speleers. B-Splines and Spline Approximation. InSplines and PDEs: From Approximation Theory to Numerical Linear Algebra, volume 2219 ofLecture Notes in Mathematics, pages 1–57. Springer, Cham, 2018

    work page 2018

  23. [23]

    B.-G. Lee, T. Lyche, and K. Mørken. Some examples of quasi-interpolants constructed from local spline projectors. InMathematical Methods for Curves and Surfaces: Oslo 2000, pages 243–252. Vanderbilt University, USA, 1 edition, January 2001

    work page 2000

  24. [24]

    He and L

    L. He and L. Tobiska. The two-level local projection stabilization as an enriched one-level approach.Advances in Computational Mathematics, 36(4):503–523, September 2011

    work page 2011

  25. [25]

    Knobloch

    P. Knobloch. A Generalization of the Local Projection Stabilization for Convection-Diffusion-Reaction Equations. SIAM Journal on Numerical Analysis, 48(2):659–680, 2010

    work page 2010

  26. [26]

    Barrenechea, E

    G.R. Barrenechea, E. Burman, and F. Karakatsani. Blending low-order stabilised finite element methods: A positivity-preserving local projection method for the convection–diffusion equation.Computer Methods in Applied Mechanics and Engineering, 317:1169–1193, April 2017

    work page 2017

  27. [27]

    Matthies and L

    G. Matthies and L. Tobiska. Local projection type stabilization applied to inf–sup stable discretizations of the Oseen problem.IMA Journal of Numerical Analysis, 35(1):239–269, January 2015

    work page 2015

  28. [28]

    Schoenberg

    I.J. Schoenberg. Contributions to the Problem of Approximation of Equidistant Data by Analytic Functions. Quarterly of Applied Mathematics, 4:45–99 and 112–141, 1946

    work page 1946

  29. [29]

    de Boor.A Practical Guide to Splines

    C. de Boor.A Practical Guide to Splines. Springer, New York, revised edition, 2001

    work page 2001

  30. [30]

    de Boor and G.J

    C. de Boor and G.J. Fix. Spline approximation by quasiinterpolants.Journal of Approximation Theory, 8(1):19–45, May 1973

    work page 1973

  31. [31]

    Lyche and L.L

    T. Lyche and L.L. Schumaker. Local spline approximation methods.Journal of Approximation Theory, 15(4):294– 325, December 1975

    work page 1975

  32. [32]

    Beirão da Veiga, A

    L. Beirão da Veiga, A. Buffa, G. Sangalli, and R. Vázquez. Mathematical analysis of variational isogeometric methods.Acta Numerica, 23:157–287, 2014

    work page 2014

  33. [33]

    Schumaker.Spline Functions: Basic Theory

    L.L. Schumaker.Spline Functions: Basic Theory. Cambridge University Press, Cambridge, 3rd edition, 2007

    work page 2007

  34. [34]

    Cottrell, T.J.R

    J.A. Cottrell, T.J.R. Hughes, and Y . Bazilevs.Isogeometric Analysis: Toward Integration of CAD and FEA. John Wiley & Sons, Chichester, 2009. 27 Stabilized Multilevel Isogeometric Quasi-Interpolant Projection Method

    work page 2009

  35. [35]

    Bressan and G

    A. Bressan and G. Sangalli. Isogeometric discretizations of the Stokes problem: stability analysis by the macroelement technique.IMA Journal of Numerical Analysis, 33(2):629–651, April 2013

    work page 2013

  36. [36]

    Matthies, P

    G. Matthies, P. Skrzypacz, and L. Tobiska. Stabilization of local projection type applied to convection-diffusion problems with mixed boundary conditions.Electronic Transactions on Numerical Analysis, 2008

    work page 2008

  37. [37]

    John and P

    V . John and P. Knobloch. A computational comparison of methods diminishing spurious oscillations in finite element solutions of convection–diffusion equations. InPrograms and Algorithms of Numerical Mathematics 13, pages 122–136, Prague, 2006. Institute of Mathematics AS CR

    work page 2006

  38. [38]

    John and P

    V . John and P. Knobloch. On spurious oscillations at layers diminishing (SOLD) methods for convection–diffusion equations: Part II – Analysis for P1 and Q1 finite elements.Computer Methods in Applied Mechanics and Engineering, 197(21):1997–2014, April 2008

    work page 1997

  39. [39]

    Guermond

    J.-L. Guermond. Stabilization of Galerkin approximations of transport equations by subgrid modeling.ESAIM: Mathematical Modelling and Numerical Analysis, 33:1293–1316, 1999

    work page 1999

  40. [40]

    Garg and S

    D. Garg and S. Ganesan. Local projection stabilized finite element methods for advection-reaction problems. Calcolo, 60(4):45, 2023. 28

    work page 2023