Laplacian regularized eikonal equation with Soner boundary condition on polyhedral meshes

by   Jooyoung Hahn, et al.

In this paper, we propose a numerical algorithm based on a cell-centered finite volume method to compute a distance from given objects on a three-dimensional computational domain discretized by polyhedral cells. Inspired by the vanishing viscosity method, a Laplacian regularized eikonal equation is solved and the Soner boundary condition is applied to the boundary of the domain to avoid a non-viscosity solution. As the regularization parameter depending on a characteristic length of the discretized domain is reduced, a corresponding numerical solution is calculated. A convergence to the viscosity solution is verified numerically as the characteristic length becomes smaller and the regularization parameter accordingly becomes smaller. From the numerical experiments, the second experimental order of convergence in the L^1 norm error is confirmed for smooth solutions. Compared to solve a time-dependent form of eikonal equation, the Laplacian regularized eikonal equation has the advantage of reducing computational cost dramatically when a more significant number of cells is used or a region of interest is far away from the given objects. Moreover, the implementation of parallel computing using domain decomposition with 1-ring face neighborhood structure can be done straightforwardly by a standard cell-centered finite volume code.


page 12

page 13

page 14

page 15


Convergence analysis of a Crank-Nicolson Galerkin method for an inverse source problem for parabolic systems with boundary observations

This work is devoted to an inverse problem of identifying a source term ...

Convergence of a regularized finite element discretization of the two-dimensional Monge-Ampère equation

This paper proposes a regularization of the Monge-Ampère equation in pla...

Convergence analysis of a numerical scheme for a tumour growth model

We consider a one-spatial dimensional tumour growth model that consists ...

Numerical solution of a two dimensional tumour growth model with moving boundary

We consider a biphasic continuum model for avascular tumour growth in tw...

Iterative solution to the biharmonic equation in mixed form discretized by the Hybrid High-Order method

We consider the solution to the biharmonic equation in mixed form discre...

Upstream mobility Finite Volumes for the Richards equation in heterogenous domains

This paper is concerned with the Richards equation in a heterogeneous do...

4D Segmentation Algorithm with application to 3D+time Image Segmentation

In this paper, we introduce and study a novel segmentation method for 4D...

Please sign up or login with your details

Forgot password? Click here to reset