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

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