phi-FD, a second order finite difference scheme for geometries defined by a level-set function: the Neumann case

Abstract

We extend the phi-FD finite-difference scheme, which handles domains described by a level-set function on a Cartesian grid, to Neumann boundary conditions. The condition is relaxed near the boundary through the level-set and an auxiliary field, which is eliminated by writing the relaxed condition at a boundary node and at its nearest interior node. The scheme uses only interior nodes and the level-set: no ghost values, flux reconstruction or cut-cell volumes. We prove optimal consistency and show second-order convergence in the L2, H1 and Linf norms on a non-convex domain. Two elementary precautions, an individual normalization of the boundary equations and a rule for the choice of the coupled interior node, make the condition number O(h^-2) and independent of the smallest cut cell, without any stabilization term.