The asymptotic analysis developed in this paper guarantees that the following relaxed equation with a space-dependent relaxation parameter \(\varepsilon a(x)\partial u_ \varepsilon/\partial t=\varepsilon\text{div}(a(x)\nabla_ x u_ \varepsilon)-\psi(u_ \varepsilon)/(2\varepsilon a(x))\) in \(\mathbb{R}^ n\times (0,T)\), coupled with the initial condition \(u_ \varepsilon(.,0)=\chi_ \varepsilon(.)\) in \(\mathbb{R}^ n\), approximates a flow by mean curvature, as \(\varepsilon\to 0\).
It is shown that the zero-level surface of the relaxed solution approximates, with an error of order \(O(\varepsilon^ 2)\), a surface evolving according to the mean curvature motion. The mean curvature evolution of the boundary of various 2D and axisymmetric 3D sets is simulated numerically. The numerous numerical experiments confirm the reliability of the asymptotic result.