pith. sign in

arxiv: 2510.18772 · v2 · submitted 2025-10-21 · 🧮 math.NA · cs.NA· physics.comp-ph

Adaptive hyperviscosity stabilisation for the RBF-FD method in solving advection-dominated transport equations

Pith reviewed 2026-05-18 04:32 UTC · model grok-4.3

classification 🧮 math.NA cs.NAphysics.comp-ph
keywords RBF-FDhyperviscosity stabilisationadaptive algorithmadvection-dominated transportspectral radiusnumerical stabilitymeshless methods
0
0 comments X

The pith

An adaptive algorithm sets the hyperviscosity constant for RBF-FD schemes from the spectral radius of the evolution matrix to stabilize advection-dominated equations without PDE-specific tuning.

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

The paper develops an adaptive stabilisation method for radial basis function finite difference discretizations of linear and nonlinear advection-dominated transport problems on unbounded domains. It computes the hyperviscosity coefficient directly from the spectral radius of the discrete evolution matrix, producing a choice that depends neither on the particular PDE nor on empirical constants or von Neumann analysis. The procedure works for arbitrary node layouts and reduces cost by allowing lower monomial augmentation in the hyperviscosity operator while still preserving stability. A hybrid choice of spline orders between the advection and hyperviscosity operators is shown to further improve robustness. Numerical tests on pure advection and the Burgers equation confirm stable evolution with limited added dissipation.

Core claim

The central claim is that the spectral radius of the RBF-FD evolution matrix supplies a reliable, problem-independent scale for the hyperviscosity constant, enabling a general stabilisation procedure that supports arbitrary node sets and both linear and nonlinear advection-dominated equations while avoiding empirical tuning and von Neumann estimates; lower monomial augmentation in the hyperviscosity operator and hybrid spline orders maintain consistency and efficiency.

What carries the argument

The adaptive hyperviscosity constant determined from the spectral radius of the RBF-FD evolution matrix, together with reduced monomial augmentation for the viscosity operator.

If this is right

  • The stabilisation procedure applies to general node layouts and is not restricted to specific equations.
  • Lower monomial augmentation in the hyperviscosity operator preserves stability while permitting smaller stencils and lower computational cost.
  • A hybrid strategy with different spline orders for advection and hyperviscosity operators improves overall stability.
  • Stable performance with limited numerical dissipation is obtained for both pure linear advection and the nonlinear Burgers equation.

Where Pith is reading between the lines

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

  • The spectral-radius approach could extend naturally to time-evolving node distributions where retuning would otherwise be needed at each step.
  • Similar spectral information might be used to adapt stabilisation parameters in other meshless or particle-based schemes for advection-dominated flows.
  • Application to bounded domains would require consistent treatment of boundary conditions alongside the adaptive hyperviscosity scaling.

Load-bearing premise

The spectral radius of the RBF-FD evolution matrix supplies a reliable, problem-independent measure for setting the hyperviscosity constant that stabilizes the scheme without introducing excessive numerical dissipation across varied node sets and both linear and nonlinear equations.

What would settle it

A test on an advection-dominated problem or node distribution not examined in the paper in which the adaptively chosen constant either permits growing oscillations or produces visibly larger smoothing than a well-tuned fixed value.

Figures

Figures reproduced from arXiv: 2510.18772 by Andrej Kolar-Po\v{z}un, Gregor Kosec, Miha Rot, \v{Z}iga Vaupoti\v{c}.

Figure 1
Figure 1. Figure 1: Eigenvalue spectra of the stabilised advection operator [PITH_FULL_IMAGE:figures/full_fig_p012_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Spectral radius ρ of the evolution matrix Gh (19) at ∆t = 10−4 with respect to c ∈ [10−4 , 102 ] for internodal distances h ∈ {0.01, 0.02, 0.03, 0.04} and orders of hyper￾viscosity α ∈ {2, 3, 4}. The solid line shows ρ calculated with the full order monomial augmentation as suggested for the RBF-FD approximation while the dotted line displays effects of using a reduced monomial order of m = 2 with stencil … view at source ↗
Figure 3
Figure 3. Figure 3: Eigenvalue spectra of the stabilised advection operator [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Eigenvalue spectra of the advection operator [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Solution of the linear advection equation at two different times [PITH_FULL_IMAGE:figures/full_fig_p019_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Relative energy and error of numerical simulation of the pure advection with [PITH_FULL_IMAGE:figures/full_fig_p020_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Relative error of pure advection equation at [PITH_FULL_IMAGE:figures/full_fig_p021_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Relative energy of the hyperviscosity stabilized solution discretised with [PITH_FULL_IMAGE:figures/full_fig_p021_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Error convergence under h-refinement for linear advection equation at t = 5 with ∆t = 10−5 . Different lines represent different orders of monomial augmentation, while columns represent the order of the hyperviscosity operator. The black lines are fitted based on the ending 8 markers without the last 2 markers with h k + d. Next, we take a look at how the stencil size and different orders of poly￾nomial au… view at source ↗
Figure 10
Figure 10. Figure 10: Error of pure advection numerical simulation with respect to stencil size of the [PITH_FULL_IMAGE:figures/full_fig_p023_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Error of pure advection numerical simulation with respect to [PITH_FULL_IMAGE:figures/full_fig_p023_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Comparison between stabilised (centre) and non-stabilised (right) Burgers’ [PITH_FULL_IMAGE:figures/full_fig_p025_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Error for the case of Burgers’ equation at [PITH_FULL_IMAGE:figures/full_fig_p026_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Relative energy with respect to c at h = 0.02 and ∆t = 0.01. Different lines represent sampling at different times of the simulation. The columns show simulation at different Reynolds numbers. As the linearisation of the advection term causes the matrix Gh(tn) to become time-dependent, a natural question arises regarding the required fre￾quency of copt re-computation. We can assume that the number of time… view at source ↗
Figure 15
Figure 15. Figure 15: Time evolution of relative error of Burgers’ equation at [PITH_FULL_IMAGE:figures/full_fig_p028_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Convergence of the approximation error as a function of [PITH_FULL_IMAGE:figures/full_fig_p029_16.png] view at source ↗
read the original abstract

This paper presents an adaptive hyperviscosity stabilisation procedure for the Radial Basis Function-generated Finite Difference (RBF-FD) method, aimed at solving linear and non-linear advection-dominated transport equations on domains without a boundary. The approach employs a PDE-independent algorithm that adaptively determines the hyperviscosity constant based on the spectral radius of the RBF-FD evolution matrix. The proposed procedure supports general node layouts and is not tailored for specific equations, avoiding the limitations of empirical tuning and von Neumann-based estimates. To reduce computational cost, it is shown that lower monomial augmentation in the approximation of the hyperviscosity operator can still ensure consistent stabilisation, enabling the use of smaller stencils and improving overall efficiency. A hybrid strategy employing different spline orders for the advection and hyperviscosity operators is also implemented to enhance stability. The method is evaluated on pure linear advection and non-linear Burgers' equation, demonstrating stable performance with limited numerical dissipation. The two main contributions are: (1) a general hyperviscosity RBF-FD solution procedure demonstrated on both linear and non-linear advection-dominated problems, and (2) an in-depth analysis of the behaviour of hyperviscosity within the RBF-FD framework, addressing the interplay between key free parameters and their influence on numerical results.

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 paper proposes an adaptive hyperviscosity stabilization procedure for RBF-FD discretizations of advection-dominated transport equations on domains without boundaries. The hyperviscosity constant is chosen adaptively from the spectral radius of the RBF-FD evolution matrix in a claimed PDE-independent manner that supports general node sets. Additional contributions include showing that lower-order monomial augmentation suffices for the hyperviscosity operator and implementing a hybrid spline strategy; numerical tests on linear advection and nonlinear Burgers' equation are reported to yield stable solutions with limited dissipation.

Significance. If the adaptive spectral-radius procedure can be shown to deliver stable, low-dissipation solutions for both linear and nonlinear problems on irregular nodes without problem-specific tuning, the work would provide a practical advance in meshless methods for advection-dominated flows and reduce reliance on empirical or von Neumann estimates.

major comments (2)
  1. [Abstract and numerical experiments] The central claim that the spectral radius of the (time-dependent) RBF-FD evolution matrix supplies a reliable, problem-independent hyperviscosity coefficient is not supported by any derivation, stability bound, or quantitative dissipation measure. For Burgers' equation the matrix changes at each step, yet the manuscript provides no analysis of recomputation frequency or uniform application of the coefficient, leaving open the possibility that local truncation error or stability is violated on irregular nodes.
  2. [Abstract] The assertion that lower monomial augmentation for the hyperviscosity operator still ensures 'consistent stabilisation' is presented without supporting analysis or comparison of truncation errors; this choice is load-bearing for the efficiency claim but lacks justification that stability properties are preserved.
minor comments (2)
  1. The interplay between free parameters (augmentation order, spline orders, radius-to-constant mapping) is discussed qualitatively; explicit equations or pseudocode for the adaptive algorithm would improve reproducibility.
  2. [Numerical results] No tables or figures quantify dissipation (e.g., L2 error growth or energy decay rates) across node sets, making the 'limited numerical dissipation' statement difficult to verify.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their insightful comments and the opportunity to clarify and strengthen our manuscript. We respond to each major comment below, providing explanations based on the content of our work and indicating planned revisions.

read point-by-point responses
  1. Referee: The central claim that the spectral radius of the (time-dependent) RBF-FD evolution matrix supplies a reliable, problem-independent hyperviscosity coefficient is not supported by any derivation, stability bound, or quantitative dissipation measure. For Burgers' equation the matrix changes at each step, yet the manuscript provides no analysis of recomputation frequency or uniform application of the coefficient, leaving open the possibility that local truncation error or stability is violated on irregular nodes.

    Authors: While the manuscript does not include a formal derivation or stability bound, the adaptive spectral-radius procedure is motivated by ensuring the hyperviscosity term counters unstable modes in the discrete operator spectrum and is shown to be effective and PDE-independent through numerical tests on general node sets for both linear advection and Burgers' equation. For the nonlinear case, the evolution matrix is recomputed at each time step to adapt the coefficient to the changing dynamics. We will revise the relevant sections to provide a clearer description of the recomputation strategy and include additional quantitative dissipation measures, such as energy norms and error comparisons, to better address potential concerns on irregular nodes. revision: yes

  2. Referee: The assertion that lower monomial augmentation for the hyperviscosity operator still ensures 'consistent stabilisation' is presented without supporting analysis or comparison of truncation errors; this choice is load-bearing for the efficiency claim but lacks justification that stability properties are preserved.

    Authors: Numerical results in the manuscript demonstrate that lower monomial augmentation for the hyperviscosity operator maintains consistent stabilization while allowing smaller stencils and reduced cost. We acknowledge the absence of explicit truncation error comparisons or formal analysis. In the revised manuscript we will add targeted numerical comparisons of truncation errors and stability metrics across augmentation orders to provide stronger empirical justification for the efficiency claim. revision: yes

Circularity Check

0 steps flagged

No circularity: adaptive hyperviscosity scale uses external spectral property of the base operator

full rationale

The paper's central procedure computes the hyperviscosity coefficient from the spectral radius of the RBF-FD evolution matrix for the underlying advection operator. This is a direct functional mapping from a property of the discretization (spectral radius) to a stabilization parameter; it does not redefine the target quantity in terms of itself, fit the constant to output data, or rely on a self-citation chain for its justification. The abstract and claimed contributions present the mapping as an algorithmic choice that remains independent of specific PDEs and node sets, with no equations shown that reduce the output to a tautological restatement of the input. The method is evaluated on both linear and nonlinear test cases, confirming that the derivation chain does not collapse by construction. This is the normal, self-contained case for a stabilization heuristic.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard RBF-FD approximation properties plus an ad-hoc rule linking spectral radius to hyperviscosity scale. Several implementation parameters remain free and are explored rather than derived.

free parameters (3)
  • mapping from spectral radius to hyperviscosity constant
    The exact functional form or scaling factor that converts the computed radius into the stabilization coefficient is not derived from first principles and must be chosen or calibrated.
  • monomial augmentation order for hyperviscosity operator
    Lower orders are selected for efficiency; the choice affects both cost and the quality of stabilization and is tested rather than fixed by theory.
  • spline orders in hybrid strategy
    Different orders are assigned to advection versus hyperviscosity operators; the specific pairing is an implementation choice explored for stability.
axioms (2)
  • domain assumption The RBF-FD differentiation matrix accurately approximates the underlying differential operators on the chosen node set.
    Invoked throughout as the foundation for both the evolution matrix and the hyperviscosity operator.
  • ad hoc to paper The spectral radius of the evolution matrix is a sufficient indicator of the instability that hyperviscosity must counteract.
    This is the core premise of the adaptive algorithm and is not proven but asserted as PDE-independent.

pith-pipeline@v0.9.0 · 5788 in / 1722 out tokens · 58372 ms · 2026-05-18T04:32:27.219027+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.

Reference graph

Works this paper leans on

54 extracted references · 54 canonical work pages

  1. [1]

    1007/978-3-319-78619-3

    S.Širca, M.Horvat, ComputationalMethodsinPhysics, GraduateTexts in Physics, Springer International Publishing, Cham, 2018.doi:10. 1007/978-3-319-78619-3. URLhttp://link.springer.com/10.1007/978-3-319-78619-3

  2. [2]

    Liu, Mesh free methods: moving beyond the finite element method, CRC Press, 2002.doi:10.1201/9781420040586

    G.-R. Liu, Mesh free methods: moving beyond the finite element method, CRC Press, 2002.doi:10.1201/9781420040586

  3. [3]

    V. P. Nguyen, T. Rabczuk, S. Bordas, M. Duflot, Meshless methods: a review and computer implementation aspects, Math. Comput. Simul 79 (3) (2008) 763–813.doi:10.1016/j.matcom.2008.01.003. URLhttps://sourceforge.net/projects/elemfregalerkin/

  4. [4]

    J. Slak, G. Kosec, On generation of node distributions for meshless pde discretizations, SIAM Journal on Scientific Computing 41 (5) (2019) A3202–A3229.doi:10.1137/18m1231456. URLhttp://dx.doi.org/10.1137/18M1231456

  5. [5]

    van der Sande, B

    K. van der Sande, B. Fornberg, Fast variable density 3-D node genera- tion, SIAM Journal on Scientific Computing 43 (1) (2021) A242–A257, arXiv:1906.00636 [cs, math].doi:10.1137/20M1337016. URLhttp://arxiv.org/abs/1906.00636

  6. [6]

    O. C. Zienkiewicz, R. L. Taylor, J. Z. Zhu, The finite element method: its basis and fundamentals, Elsevier, 2005

  7. [7]

    S. M. Mirfatah, B. Boroomand, E. Soleimanifar, On the solution of 3d problems in physics: from the geometry definition in cad to the solution by a meshless method, Journal of Computational Physics 393 (2019) 351–374

  8. [8]

    U. Duh, V. Shankar, G. Kosec, Discretization of Non-uniform Ra- tional B-Spline (NURBS) Models for Meshless Isogeometric Analysis, Journal of Scientific Computing 100 (2) (2024) 51.doi:10.1007/ s10915-024-02597-z. URLhttps://link.springer.com/10.1007/s10915-024-02597-z

  9. [9]

    URLhttps://doi.org/10.1137/23M1598052

    J.A.Reeger, Adaptivityinlocalkernelbasedmethodsforapproximating the action of linear operators, SIAM Journal on Scientific Computing 32 46 (4) (2024) A2683–A2708.doi:10.1137/23M1598052. URLhttps://doi.org/10.1137/23M1598052

  10. [10]

    D. T. Oanh, O. Davydov, H. X. Phu, Adaptive rbf-fd method for ellip- tic problems with point singularities in 2d, Applied Mathematics and Computation 313 (2017) 474–497

  11. [11]

    J. Slak, G. Kosec, Adaptive radial basis function–generated finite differ- ences method for contact problems, International Journal for Numerical Methods in Engineering 119 (7) (2019) 661–686

  12. [12]

    Ortega, E

    E. Ortega, E. Oñate, S. Idelsohn, R. Flores, A meshless finite point method for three-dimensional analysis of compressible flow problems involving moving boundaries and adaptivity, International Journal for Numerical Methods in Fluids 73 (4) (2013) 323–343

  13. [13]

    Jančič, M

    M. Jančič, M. Založnik, G. Kosec, Meshless interface tracking for the simulation of dendrite envelope growth, Journal of Computational Physics 507 (2024) 112973.doi:10.1016/j.jcp.2024.112973. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0021999124002225

  14. [14]

    Shankar, R

    V. Shankar, R. M. Kirby, A. L. Fogelson, Robust node generation for mesh-free discretizations on irregular domains and surfaces, SIAM Jour- nal on Scientific Computing 40 (4) (2018) A2584–A2608

  15. [15]

    Jančič, J

    M. Jančič, J. Slak, G. Kosec, Monomial augmentation guidelines for rbf- fd from accuracy versus computational time perspective, Journal of Sci- entific Computing 87 (1) (2021) 9.doi:10.1007/s10915-020-01401-y

  16. [16]

    fi- nite difference mode

    A. Tolstykh, D. Shirobokov, On using radial basis functions in a “fi- nite difference mode” with applications to elasticity problems, Compu- tational Mechanics 33 (1) (2003) 68–79

  17. [17]

    Bayona, N

    V. Bayona, N. Flyer, B. Fornberg, G. A. Barnett, On the role of polynomials in rbf-fd approximations: Ii. numerical solution of ellip- tic pdes, Journal of Computational Physics 332 (2017) 257–273.doi: 10.1016/j.jcp.2016.12.008. 33

  18. [18]

    Bayona, An insight into rbf-fd approximations augmented with poly- nomials, Computers & Mathematics with Applications 77 (9) (2019) 2337–2353

    V. Bayona, An insight into rbf-fd approximations augmented with poly- nomials, Computers & Mathematics with Applications 77 (9) (2019) 2337–2353

  19. [19]

    Jančič, G

    M. Jančič, G. Kosec, Strong form mesh-free hp-adaptive solution of linear elasticity problem, Engineering with Computers (May 2024).doi: 10.1007/s00366-023-01843-6. URLhttps://link.springer.com/10.1007/s00366-023-01843-6

  20. [20]

    S.Shahane, A.Radhakrishnan, S.P.Vanka, Ahigh-orderaccuratemesh- less method for solution of incompressible fluid flow problems, Journal of Computational Physics 445 (2021) 110623

  21. [21]

    M. E. Nielsen, B. Fornberg, High-order numerical method for solving ellipticpartialdifferentialequationsonunfittednodesets, arXivpreprint arXiv:2407.15825 (2024)

  22. [22]

    M. Rot, M. Jančič, G. Kosec, Spatially dependent node regularity in meshlessapproximationofpartialdifferentialequations, JournalofCom- putational Science 79 (2024) 102306

  23. [23]

    Platte, T

    R. Platte, T. Driscoll, Eigenvalue stability of radial basis function dis- cretizations for time-dependent problems, Mathematics with Applica- tions 51 (8) (2006) 1251–1268.doi:10.1016/j.camwa.2006.04.007. URLhttp://dx.doi.org/10.1016/j.camwa.2006.04.007

  24. [24]

    Fornberg, E

    B. Fornberg, E. Lehto, Stabilization of rbf-generated finite difference methods for convective pdes, Journal of Computational Physics 230 (6) (2011) 2270–2285.doi:10.1016/j.jcp.2010.12.014. URLhttp://dx.doi.org/10.1016/j.jcp.2010.12.014

  25. [26]

    Tominec, E

    I. Tominec, E. Larsson, A. Heryudono, A least squares radial basis function finite difference method with improved stability properties, SIAM Journal on Scientific Computing 43 (2) (2021) A1441–A1471. 34 doi:10.1137/20m1320079. URLhttp://dx.doi.org/10.1137/20M1320079

  26. [27]

    Larsson, V

    E. Larsson, V. Shcherbakov, A. Heryudono, A least squares radial basis function partition of unity method for solving pdes, SIAM Journal on Scientific Computing 39 (6) (2017) A2538–A2563.arXiv:https://doi. org/10.1137/17M1118087,doi:10.1137/17M1118087. URLhttps://doi.org/10.1137/17M1118087

  27. [28]

    J. M. Martel, R. B. Platte, Stability of radial basis function methods for convection problems on the circle and sphere, Journal of Scientific Computing 69 (2) (2016) 487–505.doi:10.1007/s10915-016-0206-9. URLhttp://dx.doi.org/10.1007/s10915-016-0206-9

  28. [29]

    A. Golbabai, et al., Analysis on the upwind local radial basis functions method to solve convection dominated problems and it’s application for mhd flow, Engineering Analysis with Boundary Elements 100 (2019) 59–67

  29. [30]

    P. Lyra, K. Morgan, A review and comparative study of upwind biased schemes for compressible flow computation. part i: 1—d first—order schemes, Archives of Computational Methods in Engineering 7 (2000) 19–55

  30. [31]

    Kosec, B

    G. Kosec, B. Šarler, Simulation of macrosegregation with mesoseg- regates in binary metallic casts by a meshless method, En- gineering Analysis with Boundary Elements 45 (2014) 36–44. doi:10.1016/j.enganabound.2014.01.016. URLhttps://linkinghub.elsevier.com/retrieve/pii/ S0955799714000290

  31. [32]

    Ma, Chebyshev–legendre spectral viscosity method for nonlinear con- servation laws, SIAM Journal on Numerical Analysis 35 (3) (1998) 869–892.doi:10.1137/s0036142995293900

    H. Ma, Chebyshev–legendre spectral viscosity method for nonlinear con- servation laws, SIAM Journal on Numerical Analysis 35 (3) (1998) 869–892.doi:10.1137/s0036142995293900. URLhttp://dx.doi.org/10.1137/S0036142995293900

  32. [33]

    J. P. Boyd, Hyperviscous shock layers and diffusion zones: Monotonicity, spectral viscosity, and pseudospectral methods for very high order dif- ferential equations, Journal of Scientific Computing 9 (1) (1994) 81–106. doi:10.1007/BF01573179. 35

  33. [34]

    Shankar, A

    V. Shankar, A. L. Fogelson, Hyperviscosity-based stabilization for ra- dial basis function-finite difference (rbf-fd) discretizations of advec- tion–diffusion equations, Journal of Computational Physics 372 (2018) 616–639.doi:10.1016/j.jcp.2018.06.036

  34. [35]

    Shankar, G

    V. Shankar, G. B. Wright, A. Narayan, A robust hyperviscosity for- mulation for stable rbf-fd discretizations of advection-diffusion-reaction equations on manifolds, SIAM Journal on Scientific Computing 42 (4) (2020) A2371–A2401.doi:10.1137/19m1288747

  35. [36]

    Physica E: Low-dimensional Systems and Nanostructures106, 208–238 (2019) https://doi.org/10.1016/j

    N. Flyer, G. A. Barnett, L. J. Wicker, Enhancing finite differences with radial basis functions: Experiments on the navier–stokes equations, Journal of Computational Physics 316 (2016) 39–62.doi:10.1016/j. jcp.2016.02.078. URLhttp://dx.doi.org/10.1016/j.jcp.2016.02.078

  36. [37]

    Flyer, E

    N. Flyer, E. Lehto, S. Blaise, G. B. Wright, A. St-Cyr, A guide to rbf-generated finite differences for nonlinear transport: Shallow water simulations on a sphere, Journal of Computational Physics 231 (11) (2012) 4078–4095.doi:10.1016/j.jcp.2012.01.028. URLhttp://dx.doi.org/10.1016/j.jcp.2012.01.028

  37. [38]

    Tominec, M

    I. Tominec, M. Nazarov, Residual viscosity stabilized rbf-fd methods for solving nonlinear conservation laws, Journal of Scientific Computing 94 (1) (Dec. 2022).doi:10.1007/s10915-022-02055-8. URLhttp://dx.doi.org/10.1007/s10915-022-02055-8

  38. [39]

    C. A. J. Fletcher, Generating exact solutions of the two-dimensional burgers’ equations, International Journal for Numerical Methods in Flu- ids 3 (3) (1983) 213–216.doi:10.1002/fld.1650030302. URLhttp://dx.doi.org/10.1002/fld.1650030302

  39. [40]

    H. Zhu, H. Shu, M. Ding, Numerical solutions of two-dimensional burg- ers’ equations by discrete adomian decomposition method, Mathematics with Applications 60 (3) (2010) 840–848.doi:10.1016/j.camwa.2010. 05.031. URLhttp://dx.doi.org/10.1016/j.camwa.2010.05.031

  40. [41]

    Bayona, An insight into rbf-fd approximations augmented with polynomials, Mathematics with Applications 77 (9) (2019) 2337–2353

    V. Bayona, An insight into rbf-fd approximations augmented with polynomials, Mathematics with Applications 77 (9) (2019) 2337–2353. 36 doi:10.1016/j.camwa.2018.12.029. URLhttp://dx.doi.org/10.1016/j.camwa.2018.12.029

  41. [42]

    G. E. Fasshauer, Meshfree Approximation Methods with Matlab: (With CD-ROM), WORLD SCIENTIFIC, 2007.doi:10.1142/6437. URLhttp://dx.doi.org/10.1142/6437

  42. [43]

    Le Borne, W

    S. Le Borne, W. Leinen, Guidelines for rbf-fd discretization: Nu- merical experiments on the interplay of a multitude of parameter choices, Journal of Scientific Computing 95 (1) (2023) 8.doi:10.1007/ s10915-023-02123-7. URLhttps://doi.org/10.1007/s10915-023-02123-7

  43. [44]

    Vaupotič, M

    Ž. Vaupotič, M. Rot, G. Kosec, Hyperviscosity stabilisation of the rbf- fd solution to natural convection, Journal of Physics: Conference Series 2766 (1) (2024) 012160.doi:10.1088/1742-6596/2766/1/012160. URLhttp://dx.doi.org/10.1088/1742-6596/2766/1/012160

  44. [45]

    Flyer, B

    N. Flyer, B. Fornberg, V. Bayona, G. A. Barnett, On the role of polyno- mials in rbf-fd approximations: I. interpolation and accuracy, Journal of Computational Physics 321 (2016) 21–38.doi:10.1016/j.jcp.2016. 05.026. URLhttp://dx.doi.org/10.1016/j.jcp.2016.05.026

  45. [46]

    Ma, H.-y

    H.-p. Ma, H.-y. Li, Super spectral viscosity method for nonlinear con- servation laws, Journal of Shanghai University (English Edition) 10 (1) (2006) 9–14.doi:10.1007/s11741-006-0098-2. URLhttp://dx.doi.org/10.1007/s11741-006-0098-2

  46. [47]

    R. Horn, C. Johnson, Matrix Analysis, Cambridge University Press, 2012

  47. [48]

    Y. Saad, Variations on arnoldi’s method for computing eigenelements of large unsymmetric matrices, Linear Algebra and its Applications 34 (1980) 269–295.doi:10.1016/0024-3795(80)90169-x. URLhttp://dx.doi.org/10.1016/0024-3795(80)90169-X

  48. [49]

    Tadmor, Burgers’ equation with vanishing hyper-viscosity, Com- munications in Mathematical Sciences 2 (2) (2004) 317–324.doi: 10.4310/cms.2004.v2.n2.a9

    E. Tadmor, Burgers’ equation with vanishing hyper-viscosity, Com- munications in Mathematical Sciences 2 (2) (2004) 317–324.doi: 10.4310/cms.2004.v2.n2.a9. URLhttp://dx.doi.org/10.4310/CMS.2004.v2.n2.a9 37

  49. [50]

    Hussein, Partial and full hyper-viscosity for navier-stokes and primitive equations, Journal of Differential Equations 269 (4) (2020) 3003–3030.doi:10.1016/j.jde.2020.02.019

    A. Hussein, Partial and full hyper-viscosity for navier-stokes and primitive equations, Journal of Differential Equations 269 (4) (2020) 3003–3030.doi:10.1016/j.jde.2020.02.019. URLhttp://dx.doi.org/10.1016/j.jde.2020.02.019

  50. [51]

    S. C. Brenner, L. R. Scott, The Mathematical Theory of Fi- nite Element Methods, Springer New York, 2008.doi:10.1007/ 978-0-387-75934-0. URLhttp://dx.doi.org/10.1007/978-0-387-75934-0

  51. [52]

    Kolar-Požun, M

    A. Kolar-Požun, M. Jančič, G. Kosec, A superconvergence result in the rbf-fd method, Journal of Physics: Conference Series 2766 (1) (2024) 012161.doi:10.1088/1742-6596/2766/1/012161. URLhttp://dx.doi.org/10.1088/1742-6596/2766/1/012161

  52. [53]

    J. Slak, G. Kosec, Medusa: A C++ Library for Solving PDEs Using Strong Form Mesh-free Methods, ACM Transactions on Mathematical Software 47 (3) (2021) 1–25.doi:10.1145/3450966. URLhttps://dl.acm.org/doi/10.1145/3450966

  53. [54]

    Qiu, Spectra 1.1.0, a c++ library for large scale eigenvalue problems, https://spectralib.org/ (2025)

    Y. Qiu, Spectra 1.1.0, a c++ library for large scale eigenvalue problems, https://spectralib.org/ (2025)

  54. [55]

    Gunnebaud, B

    G. Gunnebaud, B. Jacob, et al., Eigen v3, http://eigen.tuxfamily.org (2010). 38