We present a compact upwind second order scheme for computing the viscosity solution of the Eikonal equation. This new scheme is based on: 1. the numerical observation that classical first order monotone upwind sche...We present a compact upwind second order scheme for computing the viscosity solution of the Eikonal equation. This new scheme is based on: 1. the numerical observation that classical first order monotone upwind schemes for the Eikonal equation yield numerical upwind gradient which is also first order accurate up to singularities; 2. a remark that partial information on the second derivatives of the solution is known and given in the structure of the Eikonal equation and can be used to reduce the size of the stencil. We implement the second order scheme as a correction to the well known sweeping method but it should be applicable to any first order monotone upwind scheme. Care is needed to choose the appropriate stencils to avoid instabilities.展开更多
基金partially supported by ONR Grant N00014-02-1-0090ARO MURI Grant W911NF-07-1-0185NSF Grant DMS0811254
文摘We present a compact upwind second order scheme for computing the viscosity solution of the Eikonal equation. This new scheme is based on: 1. the numerical observation that classical first order monotone upwind schemes for the Eikonal equation yield numerical upwind gradient which is also first order accurate up to singularities; 2. a remark that partial information on the second derivatives of the solution is known and given in the structure of the Eikonal equation and can be used to reduce the size of the stencil. We implement the second order scheme as a correction to the well known sweeping method but it should be applicable to any first order monotone upwind scheme. Care is needed to choose the appropriate stencils to avoid instabilities.