A Pressure-Robust Immersed Interface Method for Discrete Surfaces
Fuente:
arXiv
Enregistré dans:
| Auteurs principaux: | , , |
|---|---|
| Format: | Preprint |
| Publié: |
2026
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866908871677706240 |
|---|---|
| author | Facci, Michael J. Sun, Qi Griffith, Boyce E. |
| author_facet | Facci, Michael J. Sun, Qi Griffith, Boyce E. |
| contents | The immersed interface method (IIM) for fluid-structure interaction imposes discontinuities in the fluid stress along immersed boundaries that are generated by forces concentrated along those boundaries. For a viscous incompressible fluid, imposing these discontinuities requires decomposing the boundary force into its normal and tangential components, which determine jump conditions for the pressure and velocity gradient. Previously, we developed an IIM for C0 triangulated surfaces, with a focus on piecewise linear surface representations. In this setting, the normal and tangent vectors of the discrete surface are constant on each element, and that method uses those piecewise constant vectors to determine the normal and tangential force components and, ultimately, the jump condition. We demonstrated that this is substantially more accurate than immersed boundary methods that use regularized delta functions at corresponding grid resolutions for situations in which shear stresses dominate. However, this IIM formulation struggles to accurately capture pressure loads. Here, we identify that the primary cause of this limitation is the discontinuous surface normal inherent in C0 triangulated surfaces. We propose a procedure that uses approximations of the surface normals that more accurately account for the curvature and avoid discontinuities in the reconstructed normal vectors. We investigate two ways to construct the continuous surface normal approximation: a standard L2 projection of the discontinuous surface normal field into a continuous finite element space, and constructing vertex normal vectors using inverse centroid-distance weighting and applying linear interpolation across each element. Numerical experiments show that the use of jump conditions computed with reconstructed continuous normal vector fields reduce leakage by up to six orders of magnitude across a range of pressures. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2603_07248 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | A Pressure-Robust Immersed Interface Method for Discrete Surfaces Facci, Michael J. Sun, Qi Griffith, Boyce E. Numerical Analysis Fluid Dynamics The immersed interface method (IIM) for fluid-structure interaction imposes discontinuities in the fluid stress along immersed boundaries that are generated by forces concentrated along those boundaries. For a viscous incompressible fluid, imposing these discontinuities requires decomposing the boundary force into its normal and tangential components, which determine jump conditions for the pressure and velocity gradient. Previously, we developed an IIM for C0 triangulated surfaces, with a focus on piecewise linear surface representations. In this setting, the normal and tangent vectors of the discrete surface are constant on each element, and that method uses those piecewise constant vectors to determine the normal and tangential force components and, ultimately, the jump condition. We demonstrated that this is substantially more accurate than immersed boundary methods that use regularized delta functions at corresponding grid resolutions for situations in which shear stresses dominate. However, this IIM formulation struggles to accurately capture pressure loads. Here, we identify that the primary cause of this limitation is the discontinuous surface normal inherent in C0 triangulated surfaces. We propose a procedure that uses approximations of the surface normals that more accurately account for the curvature and avoid discontinuities in the reconstructed normal vectors. We investigate two ways to construct the continuous surface normal approximation: a standard L2 projection of the discontinuous surface normal field into a continuous finite element space, and constructing vertex normal vectors using inverse centroid-distance weighting and applying linear interpolation across each element. Numerical experiments show that the use of jump conditions computed with reconstructed continuous normal vector fields reduce leakage by up to six orders of magnitude across a range of pressures. |
| title | A Pressure-Robust Immersed Interface Method for Discrete Surfaces |
| topic | Numerical Analysis Fluid Dynamics |
| url | https://arxiv.org/abs/2603.07248 |