A Pressure-Robust Immersed Interface Method for Discrete Surfaces

TL;DR

This paper introduces a pressure-robust immersed interface method that uses continuous surface normal approximations to improve pressure load accuracy in fluid-structure interaction simulations involving discrete surfaces.

Contribution

It proposes a novel procedure for constructing continuous surface normals, significantly enhancing pressure load accuracy compared to previous methods with discontinuous normals.

Findings

Reduced pressure leakage by up to six orders of magnitude

Improved accuracy in capturing pressure loads on discrete surfaces

Demonstrated effectiveness across various pressure ranges

Abstract

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.