Handling Neumann and Robin boundary conditions in a fictitious domain volume penalization framework

Abstract

Sakurai et al. (J Comput Phys, 2019) presented a flux-based volume penalization (VP) approach for imposing inhomogeneous Neumann boundary conditions on embedded interfaces. The flux-based VP method modifies the diffusion coefficient of the original elliptic (Poisson) equation and uses a flux-forcing function as a source term in the equation to impose the Neumann boundary conditions. As such, the flux-based VP method can be easily incorporated into existing fictitious domain codes. Sakurai et al. relied on an analytical construction of flux-forcing functions, which limits the practicality of the approach. Because of the analytical approach taken in the prior work, only (spatially) constant flux values on simple interfaces were considered. In this paper, we present a numerical technique for constructing flux-forcing functions for arbitrarily complex boundaries. The imposed flux values are also allowed to vary spatially in our approach. Furthermore, the flux-based VP method is extended to include (spatially varying) Robin boundary conditions, which makes the flux-based VP method even more general. We consider several two- and three-dimensional test examples to access the spatial accuracy of the numerical solutions. The method is also used to simulate flux-driven thermal convection in a concentric annular domain. We formally derive the flux-based volume penalized Poisson equation satisfying Neumann/Robin boundary condition in strong form; such a derivation was not presented in Sakurai et al., where the equation first appeared for the Neumann problem. The derivation reveals that the flux-based VP approach relies on a surface delta function to impose inhomogeneous Neumann/Robin boundary conditions. However, explicit construction of the delta function is not necessary for the flux-based VP method, which makes it different from other diffuse domain equations presented in the literature.

References in corpus (11)

• Immersed boundary model of aortic heart valve dynamics with physiological driving and loading conditions
• A robust incompressible Navier-Stokes solver for high density ratio multiphase flows
• A DLM immersed boundary method based wave-structure interaction solver for high density ratio multiphase flows
• Simulating water-entry/exit problems using Eulerian-Lagrangian and fully-Eulerian fictitious domain methods within the open-source IBAMR library
• The inertial sea wave energy converter (ISWEC) technology: device-physics, multiphase modeling and simulations
• Comparison of wave-structure interaction dynamics of a submerged cylindrical point absorber with three degrees of freedom using potential flow and computational fluid dynamics models
• A fully resolved active musculo-mechanical model for esophageal transport
• Immersed boundary methods for numerical simulation of confined fluid and plasma turbulence in complex geometries: a review
• Volume penalization for inhomogeneous Neumann boundary conditions modeling scalar flux in complicated geometry
• Analysis and discretization of the volume penalized Laplace operator with Neumann boundary conditions
• Critique on "Volume penalization for inhomogeneous Neumann boundary conditions modeling scalar flux in complicated geometry"