Abstract:This paper introduces an approach to improve volume conservation in the immersed boundary (IB) method using regularized delta functions derived from composite B-splines. These delta functions employ tensor product kernels using B-splines, whose polynomial degrees vary in normal and tangential directions based on the corresponding velocity component. Our method addresses the long-standing volume conservation issues in the conventional IB method, particularly evident in simulations of pressurized, closed membranes. We demonstrate that our approach significantly enhances volume conservation, rivaling the performance of the non-local Divergence-Free Immersed Boundary (DFIB) method introduced by Bao et al. while maintaining the local nature of the classical IB method. This avoids the computational overhead associated with the DFIB method's construction of an explicit velocity potential which requires additional Poisson solves. Numerical experiments show that sufficiently regular composite B-spline kernels can maintain initial volumes to within machine precision. We analyze the relationship between kernel regularity and the accuracy of force spreading and velocity interpolation operations. Our findings indicate that composite B-splines of at least $C^1$ regularity produce results comparable to the DFIB method in dynamic simulations, with volume conservation errors primarily dominated by the time-stepping scheme's truncation error. This work offers a computationally efficient alternative for improving volume conservation in IB methods, particularly beneficial for large-scale, three-dimensional simulations. The proposed approach requires minimal modifications to an existing IB code, making it an accessible improvement for a wide range of applications in computational fluid dynamics and fluid-structure interaction.

Higher order divergence-free and curl-free interpolation on MAC grids

Surface Interpolation Scheme by Distance Blending over Convex Sets

Curl-Flow: Boundary-Respecting Pointwise Incompressible Velocity Interpolation for Grid-Based Fluids

A Partition of Unity Method for Divergence-free or Curl-free Radial Basis Function Approximation

Local decay rates of best-approximation errors using vector-valued finite elements for fields with low regularity and integrable curl or divergence

Higher-order conservative discretizations on arbitrarily varying non-uniform grids

Implementation Of The Divergence-Free And Pressure-Oscillation-Free Projection Method For Solving The Incompressible Navier-Stokes Equations On The Collocated Grids

An Efficient Second-Order Accurate and Continuous Interpolation for Block-Adaptive Grids

A highly efficient and accurate divergence-free spectral method for curl-curl equation in two and three dimensions

A new adaptive Cartesian-grid CIP method for computation of violent free-surface flows

Locally Structure-Preserving div-curl operators for high order Discontinuous Galerkin schemes

Quasi-interpolation for the Helmholtz-Hodge decomposition

Local divergence-free velocity interpolation for the immersed boundary method using composite B-splines

An Efficient Quadratic Interpolation Scheme for a Third-Order Cell-Centered Finite-Volume Method on Tetrahedral Grids

A comparison of interpolation techniques for non-conformal high-order discontinuous Galerkin methods

A Conservative Solver for Surface-Tension-driven Multiphase Flows on Collocated Unstructured Grids

High order Finite Difference/Discontinuous Galerkin schemes for the incompressible Navier-Stokes equations with implicit viscosity

Moving Mesh Finite Element Methods for the Incompressible Navier-Stokes Equations.

Finite difference methods for 3D viscous incompressible flows in the vorticity-vector potential formulation on nonstaggered grids

Exponentially-fitted finite elements for $H({\rm curl})$ and $H({\rm div})$ convection-diffusion problems

A variational discretization method for mean curvature flows by the Onsager principle