Stabilized isogeometric formulation of the Stokes problem on overlapping patches

Xiaodong Wei,Riccardo Puppi,Pablo Antolin,Annalisa Buffa
2023-09-21
Abstract:We present a novel stabilized isogeometric formulation for the Stokes problem, where the geometry of interest is obtained via overlapping NURBS (non-uniform rational B-spline) patches, i.e., one patch on top of another in an arbitrary but predefined hierarchical order. All the visible regions constitute the computational domain, whereas independent patches are coupled through visible interfaces using Nitsche's formulation. Such a geometric representation inevitably involves trimming, which may yield trimmed elements of extremely small measures (referred to as bad elements) and thus lead to the instability issue. Motivated by the minimal stabilization method that rigorously guarantees stability for trimmed geometries [1], in this work we generalize it to the Stokes problem on overlapping patches. Central to our method is the distinct treatments for the pressure and velocity spaces: Stabilization for velocity is carried out for the flux terms on interfaces, whereas pressure is stabilized in all the bad elements. We provide a priori error estimates with a comprehensive theoretical study. Through a suite of numerical tests, we first show that optimal convergence rates are achieved, which consistently agrees with our theoretical findings. Second, we show that the accuracy of pressure is significantly improved by several orders using the proposed stabilization method, compared to the results without stabilization. Finally, we also demonstrate the flexibility and efficiency of the proposed method in capturing local features in the solution field.
Numerical Analysis,Computational Engineering, Finance, and Science
What problem does this paper attempt to address?
Based on the provided text content, the main problem that this paper attempts to solve is: When using overlapping NURBS (Non - Uniform Rational B - Spline) patches to represent geometric shapes, how to deal with the instability problems caused by geometric trimming. Specifically, when a geometric shape is formed by one patch superimposed on another patch, this representation method inevitably involves trimming operations, which may lead to elements with extremely small measures (referred to as "bad elements"), thus causing numerical instability problems. These problems particularly affect the solution of the Stokes problem because the Stokes problem involves the coupling of the velocity field and the pressure field, and the performance of these two fields in "bad elements" is especially crucial. To solve these problems, the author proposes a new stabilized Isogeometric Analysis (IGA) method, which is especially suitable for the Stokes problem on overlapping patches. The core idea is to adopt different stabilization strategies for the velocity space and the pressure space: the velocity space stabilizes the flux term at the interface, while the pressure is stabilized in all "bad elements". This method not only ensures the stability of the numerical solution but also significantly improves the accuracy of the pressure field and shows the flexibility and efficiency in capturing the local characteristics of the solution field.