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 desired Neumann boundary conditions. As such, the flux-based VP method can be easily incorporated into existing fictitious domain codes. Sakurai et al. relied on 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 along simple interfaces were considered. In this paper, we present a numerical approach to construct 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 method even more general. We consider several two- and three-dimensional test examples to access the spatial accuracy of the numerical solutions. We formally derive 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, an 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.
翻译:Sakurai 等人(J Comput Phys, 2019年) 展示了一种基于通量的内流量惩罚法(VP), 用于在嵌入界面上强制设置不相容的 Neumann 边界条件。 基于通量的VP 方法改变原椭圆(Poisson) 方程式的传播系数, 并使用通量强制函数作为方程式的源术语, 以强制实施所期望的Neumann边界条件。 因此, 基于通量的VP 方法可以很容易地融入现有的虚构域域代码中。 Sakurai 等人(VP 等人) 依靠对通量配置功能的分析构建, 从而限制该方法的实用性。 由于在先前的工作中采用的分析方法, 仅( spacially) 常量通量的通量通量计算值在简单界面中修改 。 我们把基于通量基的VPPP- 的内值计算方法的进度/ dismainal 运算法的内, 的进量的进度的进度的进度 。 我们认为, 以两种进度的进度的进度的进度的进度的进度的进度的进度/进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度的进度功能。