pith. sign in

arxiv: 2512.17666 · v2 · submitted 2025-12-19 · 🧮 math.NA · cs.NA

Local h-, p-, and k-Refinement Strategies for the Isogeometric Shifted Boundary Method Using THB-Splines

Pith reviewed 2026-05-16 20:47 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords isogeometric analysisshifted boundary methodTHB-splineslocal refinementtrimmed domainsNeumann boundary conditionsconvergence ratesh-p-k refinement
0
0 comments X

The pith

Targeted local degree elevation in THB-splines mitigates the convergence reduction for Neumann boundaries in the Shifted Boundary Method.

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

This paper examines how local h-, p-, and k-refinement strategies using Truncated Hierarchical B-splines interact with the Shifted Boundary Method on trimmed geometries in isogeometric analysis. The standard SBM represents the domain through uncut elements and enforces boundary conditions with a Taylor expansion from a surrogate boundary, but Neumann conditions require extra derivatives that drop the convergence order by one. The authors introduce local p- and k-refinement schemes for THB-splines plus an enhanced shift operator that includes mixed partial derivatives, then compare their effects on accuracy, stability, and efficiency against plain h-refinement and unmodified SBM. A sympathetic reader would care because the approach avoids ill-conditioning from small cut elements and body-fitted mesh generation while offering a practical route to higher accuracy on complex immersed domains.

Core claim

The central claim is that local p- and k-refinement schemes for THB-splines, paired with an enhanced shift operator that incorporates mixed partial derivatives, change convergence behavior in SBM-based trimmed IGA formulations and that targeted degree elevation can mitigate the order reduction that standard SBM exhibits for Neumann boundary conditions, while stability is maintained across benchmark trimmed problems.

What carries the argument

The enhanced mixed-derivative shift operator applied within the Shifted Boundary Method to locally refined THB-splines.

Load-bearing premise

The local p- and k-refinement schemes together with the mixed-derivative shift operator preserve stability and optimal convergence rates across trimmed configurations without introducing new ill-conditioning.

What would settle it

A numerical test on a trimmed domain with Neumann boundary conditions in which p-refinement and the mixed-derivative operator still produce the same one-order drop in convergence rate as unmodified SBM.

Figures

Figures reproduced from arXiv: 2512.17666 by Andrea Gorgi, Christoph Hollweck, Marcus Wagner, Nicolo Antonelli, Roland W\"uchner.

Figure 1
Figure 1. Figure 1: Two-scale relation for different p-refined basis functions. The first row shows basis functions on level 0. The second row shows their corresponding children after p-refinement. The third row illustrates the scaled children whose weighted sum reproduces the parent basis function. With the two-scale relation at hand, the conceptual idea of THB-splines is identical for all refinement types. In contrast to cl… view at source ↗
Figure 2
Figure 2. Figure 2: Comparison of refinement strategies: h-, p-, and k-refinement. The first row shows the initial basis with p “ 2 and respective control points for all strategies. The second row illustrates the basis functions after global refinement for each method. The third row displays the final THB-spline basis after local refinement of the leftmost three basis functions (grey), which are replaced by finer-level child … view at source ↗
Figure 3
Figure 3. Figure 3: Refinement domains and final mesh Mˆ for local k-refinement [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: Visualization of the SBM for B-splines of order [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Surrogate domain and boundary for an immersed circle with standard B-splines of degree [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Surrogate domain for an immersed circle with THB-splines of degree [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Manufactured solution upx, yq “ xy p1 ´ xqp1 ´ yq rsinh p10p1 ´ xqp1 ´ yqq ` sinh p10xyqs ` 10 “ p1 ´ xq 2 p1 ´ yq 2 ` x 2 y 2 ‰ . 0 0.2 0.4 0.6 0.8 1 x 0 0.2 0.4 0.6 0.8 1 y |u−uh| kukL2 0 0.5 1 1.5 ×10−1 (a) no refinement, 144 DOFs, 100 Elements 0 0.2 0.4 0.6 0.8 1 x 0 0.2 0.4 0.6 0.8 1 y |u−uh| kukL2 0 0.5 1 1.5 2 ×10−2 (b) h-refinement (1 step), 294 DOFs, 250 Elements 0 0.2 0.4 0.6 0.8 1 x 0 0.2 0.4 0.… view at source ↗
Figure 9
Figure 9. Figure 9: Error distribution for various local refinement strategies. [PITH_FULL_IMAGE:figures/full_fig_p015_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Effect of the enhanced shift operator rS p d on the convergence for the manufactured solution in Eq. 30 with Dirichlet boundary conditions applied to the geometry shown in [PITH_FULL_IMAGE:figures/full_fig_p017_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Effect of the enhanced shift operator Srp d on convergence for the manufactured solution in Eq. 30, with Neumann conditions on the exterior boundary and Dirichlet conditions on the interior boundary of the geometry shown in [PITH_FULL_IMAGE:figures/full_fig_p018_11.png] view at source ↗
Figure 14
Figure 14. Figure 14: The color legend is the same as in the previous example. [PITH_FULL_IMAGE:figures/full_fig_p018_14.png] view at source ↗
Figure 12
Figure 12. Figure 12: Convergence study for the geometry in Fig. 6 with di [PITH_FULL_IMAGE:figures/full_fig_p019_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Head-to-head convergence comparison, from 20 [PITH_FULL_IMAGE:figures/full_fig_p020_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Convergence study for the geometry 6 with di [PITH_FULL_IMAGE:figures/full_fig_p021_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Head-to-head convergence comparison, from 20 [PITH_FULL_IMAGE:figures/full_fig_p022_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Convergence study for the geometry 5 with immersed exterior Neumann boundary and immersed interior Dirichlet boundary, di [PITH_FULL_IMAGE:figures/full_fig_p023_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Head-to-head convergence comparison, from 20 [PITH_FULL_IMAGE:figures/full_fig_p024_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Example visualizations of the geometries with 30 [PITH_FULL_IMAGE:figures/full_fig_p024_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: Step-by-step convergence study for the geometry in Fig. 18a, using di [PITH_FULL_IMAGE:figures/full_fig_p025_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: Step-by-step convergence study for the geometry in Fig. 18b, using di [PITH_FULL_IMAGE:figures/full_fig_p026_20.png] view at source ↗
read the original abstract

The concept of trimming, embedding, or immersing geometries into a computational background mesh has gained considerable attention in recent years, particularly in isogeometric analysis (IGA). In this approach, the physical domain is represented independently from the computational mesh, allowing the latter to be generated more easily compared with body-fitted meshes. While this facilitates the treatment of complex geometries, it also introduces challenges, such as ill-conditioning of the stiffness matrix caused by small cut elements and difficulties in accurately enforcing boundary conditions. A recently proposed technique to address these issues is the Shifted Boundary Method (SBM), which represents the computational domain solely through uncut elements and enforces boundary conditions via a Taylor expansion from a surrogate boundary to the true boundary. Previous studies have shown that, for Neumann boundary conditions, the flux evaluation requires additional derivatives in the Taylor expansion, effectively reducing the order of convergence by one. In this work, we investigate for the first time the performance of SBM combined with Truncated Hierarchical B-splines (THB-splines) under various local refinement strategies. In particular, we propose local p- and k-refinement schemes for THB-splines and compare them with local h-refinement and the unmodified SBM. Furthermore, we propose an enhanced shift operator that incorporates mixed partial derivatives, in contrast to the standard operator. The study assesses accuracy, stability, and computational efficiency for benchmark problems on trimmed domains. The results highlight how different refinement strategies affect convergence behavior in trimmed IGA formulations using SBM and demonstrate that targeted degree elevation can mitigate the Neumann boundary limitations of the standard method.

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 claims to investigate local h-, p-, and k-refinement strategies for the Isogeometric Shifted Boundary Method (SBM) using THB-splines on trimmed domains. It proposes an enhanced shift operator incorporating mixed partial derivatives to overcome the order reduction in Neumann boundary conditions associated with the standard SBM. Through numerical benchmarks on trimmed domains, it demonstrates the impact of different refinement strategies on convergence behavior and shows that targeted degree elevation can mitigate the Neumann limitations.

Significance. If the results hold, this work advances the field by providing practical refinement strategies for immersed IGA methods, improving accuracy and efficiency for complex geometries. The enhanced operator and THB-spline refinements could lead to more stable and optimal convergent methods for problems involving Neumann conditions, building on prior SBM and THB literature.

major comments (2)
  1. [Enhanced shift operator description] The central claim that the mixed-derivative shift operator preserves consistency and stability when paired with local THB-spline refinements lacks supporting a-priori error estimates. The Taylor expansion accuracy under variable knot multiplicities and hierarchical truncation is not analyzed, placing the full burden on numerical benchmarks which may be configuration-specific.
  2. [Numerical benchmarks] The abstract asserts that benchmark results support the mitigation of Neumann order reduction via targeted p-refinement, but no quantitative error tables, convergence rates, or stability metrics are provided in the available text. This undermines the ability to assess the strength of the claims regarding generality across trimmed configurations.
minor comments (2)
  1. [Abstract] The abstract could more clearly specify the benchmark problems used and the observed convergence orders to strengthen the summary of results.
  2. [Refinement strategies] Clarify the distinction between p- and k-refinement in the context of THB-splines, as the notation may be non-standard for some readers.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and the opportunity to clarify and strengthen our manuscript. We address each major point below and indicate the planned revisions.

read point-by-point responses
  1. Referee: The central claim that the mixed-derivative shift operator preserves consistency and stability when paired with local THB-spline refinements lacks supporting a-priori error estimates. The Taylor expansion accuracy under variable knot multiplicities and hierarchical truncation is not analyzed, placing the full burden on numerical benchmarks which may be configuration-specific.

    Authors: We acknowledge that a full a-priori error analysis for the enhanced shift operator under THB-spline hierarchies would provide stronger theoretical backing. The manuscript focuses on numerical investigation of local refinement strategies; the mixed-derivative operator is derived directly from the Taylor expansion to restore consistency for Neumann data. We will add a dedicated subsection deriving the consistency of the operator, explicitly addressing the effect of knot multiplicity changes and truncation on the expansion remainder, and include additional numerical tests on configurations with varying hierarchical levels to support generality. A complete a-priori proof is beyond the present scope. revision: partial

  2. Referee: The abstract asserts that benchmark results support the mitigation of Neumann order reduction via targeted p-refinement, but no quantitative error tables, convergence rates, or stability metrics are provided in the available text. This undermines the ability to assess the strength of the claims regarding generality across trimmed configurations.

    Authors: The full manuscript contains detailed numerical results, including L2 and H1 error tables, observed convergence rates, and condition-number metrics for several trimmed geometries. We will revise the text to ensure these quantitative data are explicitly referenced from the abstract onward, add a summary table of convergence rates across all refinement strategies, and include two additional trimmed-domain examples to further illustrate generality. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation or claims.

full rationale

The paper proposes local h-, p-, and k-refinement strategies for THB-splines in the context of the Shifted Boundary Method and introduces an enhanced shift operator incorporating mixed partial derivatives. These are evaluated via numerical benchmarks on trimmed domains, with convergence behavior compared to the unmodified SBM. No equations, parameters, or central claims reduce by construction to fitted inputs, self-definitions, or load-bearing self-citations; the work builds on prior SBM and THB literature as external foundations while presenting independent numerical assessments of accuracy, stability, and efficiency.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available; no explicit free parameters, axioms, or invented entities are identifiable. The approach inherits standard assumptions from SBM (Taylor expansion validity near the boundary) and THB-splines (truncation properties) without additional postulates stated.

pith-pipeline@v0.9.0 · 5612 in / 1066 out tokens · 26346 ms · 2026-05-16T20:47:29.945373+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

57 extracted references · 57 canonical work pages

  1. [1]

    Bayesian Segmentation of Atrium Wall Using Globally-Optimal Graph Cuts on 3D Meshes

    L. Piegl, W. Tiller, The NURBS Book, Monographs in Visual Communication, Springer Berlin Heidelberg, 1997. doi:10.1007/978-3-642- 59223-2

  2. [2]

    Bathe, Finite Element Procedures, Prentice Hall, Pearson Education Inc., 2014

    K.-J. Bathe, Finite Element Procedures, Prentice Hall, Pearson Education Inc., 2014

  3. [3]

    T. J. Hughes, J. A. Cottrell, Y . Bazilevs, Isogeometric analysis: Cad, finite elements, nurbs, exact geometry and mesh refinement, Computer methods in applied mechanics and engineering 194 (39-41) (2005) 4135–4195

  4. [4]

    J. A. Cottrell, T. J. Hughes, Y . Bazilevs, Isogeometric analysis: toward integration of CAD and FEA, John Wiley & Sons, 2009

  5. [5]

    J. A. Cottrell, A. Reali, Y . Bazilevs, T. J. Hughes, Isogeometric Analysis of structural vibrations, Computer methods in applied mechanics and engineering 195 (41-43) (2006) 5257–5296

  6. [6]

    Apostolatos, R

    A. Apostolatos, R. Schmidt, R. W ¨uchner, K.-U. Bletzinger, A Nitsche-type formulation and comparison of the most common domain de- composition methods in isogeometric analysis, International Journal for Numerical Methods in Engineering 97 (7) (2014-02-17) 473–504. doi:10.1002/nme.4568

  7. [7]

    Bauer, M

    A. Bauer, M. Breitenberger, B. Philipp, R. W ¨uchner, K.-U. Bletzinger, Nonlinear isogeometric spatial Bernoulli beam, Computer Methods in Applied Mechanics and Engineering 303 (2016-05) 101–127. doi:10.1016/j.cma.2015.12.027

  8. [8]

    A. M. Bauer, R. W ¨uchner, K.-U. Bletzinger, Innovative CAD-integrated isogeometric simulation of sliding edge cables in lightweight structures, Journal of the International Association for Shell and Spatial Structures 59 (4) (2018-12-01) 251–258. doi:10.20898/j.iass.2018.198.039

  9. [9]

    Bauer, R

    A. Bauer, R. W ¨uchner, K.-U. Bletzinger, Weak coupling of nonlinear isogeometric spatial Bernoulli beams, Computer Methods in Applied Mechanics and Engineering 361 (2020) 112747. doi:10.1016/j.cma.2019.112747

  10. [10]

    Leidinger, M

    L. Leidinger, M. Breitenberger, A. Bauer, S. Hartmann, R. W ¨uchner, K.-U. Bletzinger, F. Duddeck, L. Song, Explicit dynamic isogeometric B-Rep analysis of penalty-coupled trimmed NURBS shells, Computer Methods in Applied Mechanics and Engineering 351 (2019-07) 891–

  11. [11]

    doi:10.1016/j.cma.2019.04.016

  12. [12]

    Meßmer, T

    M. Meßmer, T. Teschemacher, L. F. Leidinger, R. W ¨uchner, K.-U. Bletzinger, Efficient CAD-integrated isogeometric analysis of trimmed solids, Computer Methods in Applied Mechanics and Engineering 400 (2022-10) 115584. doi:10.1016/j.cma.2022.115584

  13. [13]

    L. Chen, W. Zhao, C. Liu, H. Chen, S. Marburg, Isogeometric fast multipole boundary element method based on Burton-Miller formulation for 3d acoustic problems, Archives of Acoustics 44 (3) (2019) 475–492. doi:10.24425/aoa.2019.129263

  14. [14]

    L. Chen, S. Marburg, W. Zhao, C. Liu, H. Chen, Implementation of isogeometric fast multipole boundary element methods for 2d half-space acoustic scattering problems with absorbing boundary condition, Journal of Theoretical and Computational Acoustics 27 (2018) 1850024. doi:10.1142/S259172851850024X

  15. [15]

    C. Ding, K. K. Tamma, H. Lian, Y . Ding, T. J. Dodwell, S. P. Bordas, Uncertainty quantification of spatially uncorrelated loads with a reduced-order stochastic isogeometric method, Computational Mechanics 67 (5) (2021) 1255–1271

  16. [16]

    N. A. Manque, J. Liedmann, F.-J. Barthold, M. A. Valdebenito, M. G. Faes, Interval isogeometric analysis for coping with geometric uncer- tainty, Computer Methods in Applied Mechanics and Engineering 437 (2025) 117773

  17. [17]

    M. E. Matzen, Isogeometrische Modellierung und Diskretisierung von Kontaktproblemen, Ph.D. thesis, Universit ¨at Stuttgart (2015). URLhttps://elib.uni-stuttgart.de/handle/11682/659

  18. [18]

    De Lorenzis, P

    L. De Lorenzis, P. Wriggers, T. J. Hughes, Isogeometric contact: a review, GAMM-Mitteilungen 37 (1) (2014) 85–123. doi:10.1002/gamm.201410005

  19. [19]

    Bazilevs, M.-C

    Y . Bazilevs, M.-C. Hsu, M. Scott, Isogeometric fluid–structure interaction Analysis with emphasis on non-matching discretiza- tions, and with application to wind turbines, Computer Methods in Applied Mechanics and Engineering 249-252 (2012) 28–41. doi:10.1016/j.cma.2012.03.028

  20. [20]

    Bazilevs, V

    Y . Bazilevs, V . M. Calo, Y . Zhang, T. J. R. Hughes, Isogeometric fluid–structure interaction Analysis with applications to arterial, Computa- tional Mechanics (2006). doi:10.1007/s00466-006-0084-3

  21. [21]

    M.-C. Hsu, I. Akkerman, Y . Bazilevs, High-performance computing of wind turbine aerodynamics using isogeometric Analysis, Computers & Fluids 49 (1) (2011) 93–100. doi:10.1016/j.compfluid.2011.05.002

  22. [22]

    Bazilevs, V

    Y . Bazilevs, V . Calo, J. Cottrell, T. Hughes, A. Reali, G. Scovazzi, Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows, Computer Methods in Applied Mechanics and Engineering 197 (1) (2007) 173–201. doi:10.1016/j.cma.2007.07.016

  23. [23]

    X. Qian, O. Sigmund, Isogeometric shape optimization of photonic crystals via coons patches, Computer Methods in Applied Mechanics and Engineering 200 (25) (2011) 2237–2255. doi:10.1016/j.cma.2011.03.007

  24. [24]

    Kang, S.-K

    P. Kang, S.-K. Youn, Isogeometric shape optimization of trimmed shell structures, Structural and Multidisciplinary Optimization (2016). doi:10.1007/s00158-015-1361-6

  25. [25]

    Kiendl, K.-U

    J. Kiendl, K.-U. Bletzinger, J. Linhard, R. W ¨uchner, Isogeometric shell analysis with Kirchhoff–Love elements, Computer Methods in Applied Mechanics and Engineering 198 (49) (2009-11) 3902–3914. doi:10.1016/j.cma.2009.08.013

  26. [26]

    Kiendl, Y

    J. Kiendl, Y . Bazilevs, M.-C. Hsu, R. W ¨uchner, K.-U. Bletzinger, The bending strip method for isogeometric analysis of Kirchhoff–Love shell structures comprised of multiple patches, Computer Methods in Applied Mechanics and Engineering 199 (37) (2010-08) 2403–2416. doi:10.1016/j.cma.2010.03.029

  27. [27]

    Kiendl, R

    J. Kiendl, R. Schmidt, R. W ¨uchner, K.-U. Bletzinger, Isogeometric shape optimization of shells using semi-analytical sensitivity analysis and sensitivity weighting, Computer Methods in Applied Mechanics and Engineering 274 (2014-06) 148–167. doi:10.1016/j.cma.2014.02.001

  28. [28]

    J. M. Kiendl, Isogeometric analysis and shape optimal design of shell structures, Ph.D. thesis, Technische Universit ¨at M¨unchen (2011). URLhttps://mediatum.ub.tum.de/doc/1002634/464162.pdf

  29. [29]

    Goswami, C

    S. Goswami, C. Anitescu, T. Rabczuk, Adaptive fourth-order phase field analysis for brittle fracture, Computer Methods in Applied Mechanics and Engineering 361 (2020) 112808. doi:https://doi.org/10.1016/j.cma.2019.112808. URLhttps://www.sciencedirect.com/science/article/pii/S0045782519307005

  30. [30]

    Rorris, A

    L. Rorris, A. Nagy, S. Hartmann, I. Chalkidis, A. Vafeidis, The ansa/ls-dyna approach for iga simulations, in: 12th European LS-DYNA 28 Conference, BETA CAE Systems International AG, Livermore Software Technology Corporation, DYNAmore GmbH, BETA CAE Systems SA, Koblenz, Germany, 2019

  31. [31]

    Rorris, I

    L. Rorris, I. Chalkidis, The latest developments of the ansa preprocessor for iga applications of ls-dyna, in: 16th International LS-DYNA® Users Conference Isogeometric Analysis, BETA CAE Systems International AG and BETA CAE Systems SA, Switzerland and Greece, 2020, pp. 1–XX

  32. [32]

    Breitenberger, A

    M. Breitenberger, A. Apostolatos, B. Philipp, R. W ¨uchner, K.-U. Bletzinger, Analysis in computer aided design: Nonlinear iso- geometric B-Rep analysis of shell structures, Computer Methods in Applied Mechanics and Engineering 284 (2015-02) 401–457. doi:10.1016/j.cma.2014.09.033

  33. [33]

    Teschemacher, A

    T. Teschemacher, A. M. Bauer, T. Oberbichler, M. Breitenberger, R. Rossi, R. W ¨uchner, K.-U. Bletzinger, Realization of CAD-integrated shell simulation based on isogeometric B-Rep analysis, Advanced Modeling and Simulation in Engineering Sciences 5 (1) (2018-12) 19. doi:10.1186/s40323-018-0109-4

  34. [34]

    Teschemacher, A

    T. Teschemacher, A. M. Bauer, R. Aristio, M. Meßmer, R. W ¨uchner, K.-U. Bletzinger, Concepts of data collection for the cad-integrated isogeometric analysis, Engineering with Computers 38 (6) (2022) 5675–5693

  35. [35]

    Parvizian, A

    J. Parvizian, A. D ¨uster, E. Rank, Finite cell method: h-and p-extension for embedded domain problems in solid mechanics, Computational Mechanics 41 (1) (2007) 121–133

  36. [36]

    E. Rank, M. Ruess, S. Kollmannsberger, D. Schillinger, A. D ¨uster, Geometric modeling, isogeometric analysis and the finite cell method, Computer Methods in Applied Mechanics and Engineering 249-252 (2012) 104–115

  37. [37]

    Leidinger, Explicit isogeometric B-Rep analysis for nonlinear dynamic crash simulations, Ph.D

    L. Leidinger, Explicit isogeometric B-Rep analysis for nonlinear dynamic crash simulations, Ph.D. thesis, Technische Universit ¨at M¨unchen (2020). URLhttps://mediatum.ub.tum.de/1542623

  38. [38]

    de Prenter, C

    F. de Prenter, C. V . Verhoosel, E. H. van Brummelen, M. G. Larson, S. Badia, Stability and conditioning of immersed finite element methods: Analysis and remedies, Archives of Computational Methods in Engineering 30 (6) (2023) 3617–3656. doi:10.1007/s11831-023-09913-0. URLhttps://doi.org/10.1007/s11831-023-09913-0

  39. [39]

    A. Main, G. Scovazzi, The shifted boundary method for embedded domain computations. part i: Poisson and stokes problems, Journal of Computational Physics 372 (2018) 972–995

  40. [40]

    N. M. Atallah, C. Canuto, G. Scovazzi, The shifted boundary method for solid mechanics, International Journal for Numerical Methods in Engineering 122 (20) (2021) 5935–5970

  41. [41]

    K. Li, A. Gorgi, R. Rossi, G. Scovazzi, The shifted boundary method for contact problems, Computer Methods in Applied Mechanics and Engineering 440 (2025) 117940

  42. [42]

    Antonelli, R

    N. Antonelli, R. Aristio, A. Gorgi, R. Zorrilla, R. Rossi, G. Scovazzi, R. W ¨uchner, The shifted boundary method in isogeometric analysis, Computer Methods in Applied Mechanics and Engineering 430 (2024) 117228

  43. [43]

    Antonelli, A

    N. Antonelli, A. Gorgi, R. Zorrilla, R. Rossi, Isogeometric analysis for non-newtonian viscoplastic fluids: challenges for non-smooth solu- tions, Computer Methods in Applied Mechanics and Engineering 447 (2025) 118386

  44. [44]

    J. H. Collins, K. Li, A. Lozinski, G. Scovazzi, Gap-sbm: A new conceptualization of the shifted boundary method with optimal convergence for the neumann and dirichlet problems, arXiv preprint arXiv:2508.09613 (2025)

  45. [45]

    T. J. R. Hughes, Y . Bazilevs, Isogeometric analysis: toward integration of CAD and FEA, Wiley, 2009

  46. [46]

    Giannelli, B

    C. Giannelli, B. J ¨uttler, H. Speleers, THB-splines: The truncated basis for hierarchical splines, Computer Aided Geometric Design 29 (7) (2012-10) 485–498. doi:10.1016/j.cagd.2012.03.025

  47. [47]

    Giannelli, B

    C. Giannelli, B. J ¨uttler, H. Speleers, Strongly stable bases for adaptively refined multilevel spline spaces, Advances in Computational Math- ematics 40 (2) (2014-04) 459–490. doi:10.1007/s10444-013-9315-2

  48. [48]

    Buffa, G

    A. Buffa, G. Gantner, C. Giannelli, D. Praetorius, R. V ´azquez, Mathematical foundations of adaptive isogeometric Analysis, Archives of Computational Methods in Engineering (Nov 2022). doi:10.1007/s11831-022-09752-5

  49. [49]

    E. M. Garau, R. V ´azquez, Algorithms for the implementation of adaptive isogeometric methods using hierarchical B-splines, Applied Nu- merical Mathematics 123 (2018) 58–87. doi:10.1016/j.apnum.2017.08.006

  50. [50]

    D. Schillinger, The p- and b-spline versions of the geometrically nonlinear finite cell method and hierarchical refinement strategies for adaptive isogeometric and embedded domain analysis, Ph.D. thesis, Technische Universit¨at M¨unchen (2012). URLhttps://mediatum.ub.tum.de/1093297

  51. [51]

    J. A. Cottrell, T. J. Hughes, A. Reali, Studies of refinement and continuity in isogeometric structural analysis, Computer methods in applied mechanics and engineering 196 (41) (2007) 4160–4183

  52. [52]

    Burman, A penalty-free nonsymmetric nitsche-type method for the weak imposition of boundary conditions, SIAM Journal on Numerical Analysis 50 (4) (2012) 1959–1981

    E. Burman, A penalty-free nonsymmetric nitsche-type method for the weak imposition of boundary conditions, SIAM Journal on Numerical Analysis 50 (4) (2012) 1959–1981

  53. [53]

    Q. Hu, F. Chouly, P. Hu, G. Cheng, S. P. Bordas, Skew-symmetric nitsche’s formulation in isogeometric analysis: Dirichlet and symmetry conditions, patch coupling and frictionless contact, Computer Methods in Applied Mechanics and Engineering 341 (2018) 188–220

  54. [54]

    Giannelli, B

    C. Giannelli, B. J ¨uttler, S. K. Kleiss, A. Mantzaflaris, B. Simeon, J. ˇSpeh, THB-splines: An effective mathematical technology for adaptive refinement in geometric design and isogeometric analysis, Computer Methods in Applied Mechanics and Engineering 299 (2016) 337–365. doi:10.1016/j.cma.2015.11.002

  55. [55]

    J. Gu, T. Yu, L. Van Lich, T.-T. Nguyen, T. Q. Bui, Adaptive multi-patch isogeometric analysis based on locally refined B-splines, Computer Methods in Applied Mechanics and Engineering 339 (2018-09) 704–738. doi:10.1016/j.cma.2018.04.013

  56. [56]

    Patrizi, C

    F. Patrizi, C. Manni, F. Pelosi, H. Speleers, Adaptive refinement with locally linearly independent LR B-splines: Theory and applications, Computer Methods in Applied Mechanics and Engineering 369 (2020-09) 113230. doi:10.1016/j.cma.2020.113230

  57. [57]

    Kumar, T

    M. Kumar, T. Kvamsdal, K. A. Johannessen, Simple a posteriori error estimators in adaptive isogeometric Analysis, Computers & Mathe- matics with Applications 70 (7) (2015-10) 1555–1582. doi:10.1016/j.camwa.2015.05.031. 29