Fact-checked by Grok 2 weeks ago

Eikonal equation

The eikonal equation is a nonlinear first-order partial differential equation that models the high-frequency limit of wave propagation, approximating the phase function of waves in media where the wavelength is much smaller than variations in the propagation speed.^[1] It takes the general form |\nabla u(\mathbf{x})| = n(\mathbf{x}), where u is the eikonal function representing the optical path length or phase, and n(\mathbf{x}) is the refractive index or inverse speed function in the medium.^[2] This equation captures the geometry of ray paths, ensuring that solutions align with the principle of least time for wave travel.^[3] Originating from geometrical optics, the eikonal equation derives from the Helmholtz equation under the short-wavelength approximation, where the wave amplitude is assumed slowly varying compared to the rapidly oscillating phase.^[2] The term "eikonal" stems from the Greek word εἰκών (eikōn), meaning "image" or "reflection," reflecting its roots in optical imaging and ray tracing. It is closely tied to Fermat's principle, which states that light rays follow paths minimizing travel time, and the eikonal equation serves as its differential formulation, enabling derivations of phenomena like Snell's law at refractive interfaces.^[3] In homogeneous media, solutions yield straight-line rays, while in inhomogeneous cases, rays curve according to variations in n(\mathbf{x}).^[2] Solutions to the eikonal equation are often viscosity solutions, a weak formulation that handles discontinuities and ensures uniqueness for boundary value problems, such as the distance function in Euclidean space satisfying |\nabla d| = 1 with d = 0 on the boundary.^[1] Numerically, it is solved using methods like fast marching or sweeping algorithms, which propagate information along characteristics in O(N \log N) time for N grid points, making it practical for large-scale computations.^[1] Beyond optics, the eikonal equation finds applications in seismology for travel-time tomography, computer vision for shape-from-shading and segmentation, medical imaging for ultrasound and MRI reconstruction, and computational fluid dynamics for front-tracking.^[1] In heterogeneous media, it models effective distances weighted by local speeds, aiding homogenization theory for wave propagation in composites. Its hyperbolic nature and connection to Hamilton-Jacobi equations underscore its role in optimal control and variational problems across physics and engineering.^[3]

Mathematical Foundations

Definition and Formulation

The eikonal equation is a first-order nonlinear partial differential equation central to the high-frequency approximation of wave propagation problems. Its standard formulation is given by

|\nabla T(\mathbf{x})| = n(\mathbf{x}),

where T(\mathbf{x}) denotes the phase function or travel time at position \mathbf{x} \in \mathbb{R}^d, \nabla T is the gradient of T, and n(\mathbf{x}) represents the refractive index or slowness, defined as the reciprocal of the local wave speed.^[4] This equation is typically supplemented with boundary conditions, such as T = 0 on an initial surface or curve \Gamma, which specifies the starting point for propagation; here, T(\mathbf{x}) interprets as the optical path length or the arrival time of the wavefront at \mathbf{x}.^[4] Variations of the equation account for different medium properties. In the isotropic case with constant speed, n(\mathbf{x}) is uniform, yielding solutions where T scales linearly with Euclidean distance from the boundary. In the inhomogeneous case with position-dependent speed, n(\mathbf{x}) varies spatially, enabling modeling of heterogeneous media where propagation speed changes with location.^[4] Geometrically, the level sets of T(\mathbf{x}) = t for constant t describe wavefronts that advance normal to themselves at speed $1/n(\mathbf{x}).^[4]^[5]

Derivation from Wave Propagation

The eikonal equation emerges as the leading-order approximation in the high-frequency asymptotic analysis of the scalar Helmholtz equation, which describes time-harmonic wave propagation in inhomogeneous media. The Helmholtz equation takes the form

\nabla^2 u + k^2 n^2(\mathbf{x}) u = 0,

where u(\mathbf{x}) is the complex-valued wave amplitude, k = \omega / c_0 is the reference wavenumber with angular frequency \omega and speed c_0, and n(\mathbf{x}) is the position-dependent refractive index.^[6]^[7] In the Wentzel-Kramers-Brillouin (WKB) approximation, valid for large k \to \infty, the solution is postulated as u(\mathbf{x}) \approx A(\mathbf{x}) \exp(i k S(\mathbf{x})), separating the slowly varying amplitude A from the rapidly oscillating phase S.^[8] Substituting this ansatz into the Helmholtz equation, collecting terms by powers of k, and retaining only the dominant O(k^2) contribution yields the eikonal equation

|\nabla S|^2 = n^2(\mathbf{x}),

with subsequent orders providing a transport equation $2 \nabla S \cdot \nabla A + (\Delta S) A = 0 for the amplitude.^[6]^[9] Here, S(\mathbf{x}) represents the optical path length, and in the section's notation, the eikonal function T coincides with S (or travel time in units where c_0 = 1). The approximation holds by neglecting lower-order terms involving second derivatives of A and \Delta S, which become insignificant when the wavelength $2\pi / k is much smaller than the scale of medium variations.^[8]^[7] This derivation extends naturally to the time-dependent scalar wave equation

\frac{\partial^2 u}{\partial t^2} = c^2(\mathbf{x}) \nabla^2 u,

governing wave propagation with local speed c(\mathbf{x}). For high-frequency solutions, assume u(\mathbf{x}, t) \approx A(\mathbf{x}, t) \exp(i \omega (t - \tau(\mathbf{x}))), where \omega \gg 1 and \tau is the phase.^[9] Inserting this form and expanding asymptotically in powers of $1/\omega produces, at leading order, the eikonal equation

|\nabla \tau|^2 = \frac{1}{c^2(\mathbf{x})},

with \tau(\mathbf{x}) interpreting as the travel time along characteristics.^[6] The next-order terms yield an amplitude transport equation analogous to the Helmholtz case, ensuring conservation of wave energy flux along rays. This asymptotic approach underscores the eikonal equation's connection to geometric optics, where waves propagate along rays perpendicular to wavefronts defined by constant phase.^[9]^[7]

Historical Development

Origins in Optics

The origins of the eikonal equation trace back to the 17th century with Pierre de Fermat's principle of least time, which posits that light travels between two points along the path that minimizes the travel time.^[10] This principle provided a variational foundation for understanding light propagation in inhomogeneous media, serving as a precursor to the eikonal's role in describing optical paths.^[11] Fermat formulated this principle in the 1650s, using it to derive the law of refraction independently of earlier empirical observations.^[12] Building on this, René Descartes had earlier proposed a related formulation in the 1630s for the laws of reflection and refraction, though his approach assumed varying light speeds in media that contrasted with Fermat's temporal minimization.^[13] These 17th-century developments established the conceptual basis for ray paths as extrema of an optical length integral, central to the eikonal framework. In the late 19th century, the term "eikonal" was introduced by Heinrich Bruns in 1895, derived from the Greek word for "image," to denote a characteristic function in variational optics that maps wavefronts to ray paths.^[14] Bruns' work formalized the eikonal as a tool for analyzing light ray bundles in complex optical systems, emphasizing its geometrical interpretation without reliance on wave theory.^[15] Early 19th-century advancements by Augustin-Jean Fresnel further connected these ideas to wavefront propagation, integrating ray optics with emerging wave descriptions of light to explain how rays remain tangent to advancing wavefronts in refracting media.^[16] Fresnel's contributions, particularly through his elastic theory of light around 1818, highlighted the continuity between geometrical rays and wave surfaces, laying groundwork for the eikonal's interpretation as the phase function in high-frequency approximations.^[16]

Evolution in Hamiltonian Mechanics

The mathematical formalization of the eikonal equation advanced significantly in the 19th century through William Rowan Hamilton's development of characteristic functions in optics during the 1830s. Hamilton introduced the principal characteristic function, often denoted as V, which represents the optical path length between two points along a light ray in a medium with varying refractive index. This function satisfies a partial differential equation that is precisely the eikonal equation, |\nabla V| = n(\mathbf{r}), where n is the refractive index. Hamilton's framework bridged optics and dynamics by deriving equations of motion for light rays analogous to those in mechanics, with the eikonal serving as the generating function for ray trajectories. His approach emphasized the variational principle underlying ray paths, minimizing the optical path length, and laid the groundwork for treating light propagation as a Hamiltonian system.^[17]^[18] Building on Hamilton's ideas, Carl Gustav Jacob Jacobi provided a key generalization in the 1840s through his work on variational calculus and dynamics. In his lectures on dynamics delivered in 1842–1843 and later published, Jacobi extended the Hamilton-Jacobi equation to systems with multiple degrees of freedom, introducing the Jacobi integral as a conserved quantity in the context of constrained variational problems. This generalization transformed the eikonal into a more versatile tool within the broader framework of analytical mechanics, where the principal function—equivalent to the eikonal—facilitates the separation of variables and solution of complex ray paths. Jacobi's contributions clarified the mathematical structure, emphasizing the equation's role in optimizing paths via least action principles, and connected optical phenomena to general mechanical systems without relying on specific optical assumptions.^[19]^[20] The eikonal equation's scope expanded into relativistic contexts with David Hilbert's 1915 formulation of the foundations of physics, where variational methods yielded the Einstein field equations of general relativity. In this framework, the eikonal describes the propagation of light as null geodesics in curved spacetime, with the optical path corresponding to the invariant interval along light rays satisfying ds^2 = 0. Hilbert's axiomatic approach integrated electromagnetic and gravitational fields variationally, implicitly incorporating the eikonal as the phase function for high-frequency wave solutions in gravitational fields, thus linking geometrical optics to the geodesic motion of photons. This connection highlighted the equation's universality, treating light rays as extremals in spacetime metrics.^[21]^[22] Twentieth-century refinements synthesized these developments in authoritative treatments of geometrical optics. Notably, Max Born and Emil Wolf's 1959 textbook Principles of Optics provided a comprehensive analysis, deriving the eikonal equation from the scalar wave equation in the short-wavelength limit and exploring its implications for ray tracing, wavefront evolution, and intensity propagation. Their work emphasized practical computations of ray paths in inhomogeneous media and aberration theory, attributing the eikonal's foundational role to Hamilton's characteristic functions while extending applications to modern optical systems. This synthesis bridged 19th-century variational insights with wave-based derivations, solidifying the eikonal's place in both theoretical and applied physics.^[23]^[24]

Physical Interpretations

Geometrical Optics and Ray Tracing

In geometrical optics, the eikonal equation governs the propagation of light rays in inhomogeneous media by describing the evolution of wavefronts, which are level sets of the eikonal function T(\mathbf{x}), the optical path length from a reference wavefront. Rays represent the orthogonal trajectories to these level sets, indicating the direction of energy flow perpendicular to the wavefronts. The path of a ray parameterized by arc length s follows the differential equation

\frac{d\mathbf{x}}{ds} = \frac{\nabla T}{|\nabla T|},

where \nabla T points normal to the wavefront and |\nabla T| = n(\mathbf{x}), with n(\mathbf{x}) denoting the refractive index. This formulation arises from Fermat's principle, which posits that rays minimize the optical path \int n \, ds, ensuring stationary travel time along the ray.^[4]^[25] A key consequence of the eikonal equation is Snell's law of refraction at interfaces between media with different refractive indices. Consider a planar interface where n jumps from n_1 to n_2; the eikonal T remains continuous across the boundary, while its tangential gradient component is preserved to maintain phase matching. For incident and refracted rays with angles \theta_1 and \theta_2 relative to the normal, the condition n_1 \sin \theta_1 = n_2 \sin \theta_2 emerges directly from the constancy of the tangential slowness p = \frac{\sin \theta}{c} = \frac{n \sin \theta}{c_0}, where c_0 is the speed in vacuum. This ensures the ray bends such that the optical path is stationary, validating the eikonal approximation for high-frequency waves.^[26]^[27] The eikonal framework breaks down at caustics, regions where rays converge and intersect, forming singularities in the ray density. These occur when the Jacobian of the ray mapping vanishes, leading to focusing points or lines where the amplitude diverges as $1/\sqrt{\delta} near the caustic, with \delta the distance from it. Such singularities, exemplified by fold or cusp caustics in lens systems, signal the limits of geometrical optics, as diffraction effects from the full wave equation become prominent beyond the high-frequency approximation.^[4]^[28] The eikonal equation aligns with Huygens' principle in the geometrical optics limit, where wavefronts advance as the envelope of secondary wavelets emanating from points on the prior wavefront, each expanding at local speed c = 1/n. Rays, as normals to these evolving level sets of T, trace the direction of wavelet propagation, capturing the normal advancement without diffraction in homogeneous or slowly varying media. This geometric interpretation unifies ray and wavefront descriptions for light propagation.^[4]

Continuous Shortest-Path Formulations

The eikonal equation arises as a Hamilton-Jacobi equation characterizing the geodesic distance function in Riemannian or Finsler geometry, providing a continuous formulation of shortest-path problems beyond wave propagation contexts. Specifically, the solution T(\mathbf{x}) to the eikonal equation |\nabla T(\mathbf{x})| = n(\mathbf{x}) with boundary condition T = 0 on a source set represents the infimum travel time from the source to point \mathbf{x}, defined variationally as

T(\mathbf{x}) = \inf_{\gamma} \int_{\gamma} n(\mathbf{r}) \, ds,

where the infimum is taken over all admissible paths \gamma connecting the source to \mathbf{x}, and ds denotes the arc length element along \gamma. This variational principle encodes the minimal "optical length" or time in media with position-dependent speed $1/n(\mathbf{x}), generalizing classical distance measures to curved or anisotropic spaces.^[29]^[30] In uniform media, where the refractive index n(\mathbf{x}) is constant, the eikonal equation simplifies to |\nabla T| = n, yielding T(\mathbf{x}) = n \cdot d_E(\mathbf{x}), with d_E the standard Euclidean distance from the source; this reflects the isotropic propagation at constant speed. For spatially varying n(\mathbf{x}), the formulation induces a Finsler or Riemannian metric on the domain, where the line element becomes ds = n(\mathbf{r}) \sqrt{d\mathbf{r} \cdot d\mathbf{r}} in the Riemannian case, defining geodesics as curves minimizing the integral and enabling shortest paths in heterogeneous environments. Such metrics capture anisotropy and inhomogeneity, with Finsler structures allowing direction-dependent lengths while preserving positive homogeneity.^[31]^[30] To handle discontinuities in n(\mathbf{x}) or non-smooth boundaries, the appropriate notion of solution is the viscosity solution, which ensures uniqueness and stability for the eikonal equation as a first-order Hamilton-Jacobi partial differential equation. Viscosity solutions guarantee that T(\mathbf{x}) coincides with the minimal arrival time, even across interfaces, by satisfying the equation in a weak sense via test functions at maxima and minima; this framework, developed by Crandall and Lions, resolves ambiguities in classical derivatives and confirms that the distance function is the unique continuous viscosity solution.^[32]^[33] This continuous shortest-path perspective links the eikonal equation to optimal control theory, where the characteristics—or rays—emerge as optimal trajectories minimizing the time integral under velocity constraints. In the Hamilton-Jacobi-Bellman framework, the eikonal corresponds to a static case of the value function for a control problem with unit-speed dynamics and cost n(\mathbf{x}), with gradient lines tracing the minimizing paths.^[34]

Applications

Geophysics and Seismology

In geophysics and seismology, the eikonal equation is fundamental for modeling seismic wave propagation, particularly in travel-time tomography, where it enables the inversion of observed first-arrival times T to recover the subsurface velocity model v(\mathbf{x}) through the relation |\nabla T| = 1/v. This approach avoids explicit ray tracing by treating the traveltime field as a solution to a Hamilton-Jacobi equation, allowing efficient computation of wavepaths that broaden to account for finite-frequency effects and off-ray sensitivity. A key advancement is the wavepath eikonal traveltime inversion (WET) method, which unifies various tomographic schemes by back-projecting residuals along curved wavepaths derived from the eikonal solution, offering improved resolution over traditional straight-ray approximations while remaining computationally efficient—nearly an order of magnitude faster than full wave-equation methods. The eikonal equation also underpins ray-based migration techniques in seismic imaging, where rays are traced as characteristics of the equation through heterogeneous, layered Earth models to position reflections accurately. In prestack depth migration, solving the eikonal provides traveltimes for ray propagation, enabling the construction of imaging operators that handle complex velocity structures like salt domes or faults in sedimentary basins. This ray-theoretic framework, rooted in high-frequency asymptotics, facilitates the simulation of wave fronts and amplitudes, essential for constructing velocity models in exploration seismology. To address anisotropy prevalent in crustal rocks, such as shales or fractured reservoirs, the eikonal equation is extended to transversely isotropic (TI) or elliptically anisotropic media, where the phase velocity depends on direction via parameters like Thomsen's \epsilon and \delta. In vertically transverse isotropy (VTI), an acoustic approximation sets the vertical shear-wave velocity to zero, simplifying the eikonal to depend primarily on the normal-moveout velocity and the \eta parameter, which captures anellipticity and improves kinematic accuracy for P-wave propagation with minimal error. This formulation is crucial for migration and tomography in anisotropic overburdens, enhancing image quality by correcting for velocity stretching in tilted or elliptical wavefronts. For more general transverse isotropy, factored forms of the eikonal allow numerical solutions that handle quasi-P and quasi-SV waves, vital for crustal velocity inversion. In earthquake studies, the eikonal equation supports first-arrival picking by computing theoretical traveltimes for P- and S-phases across networks, aiding automated detection and hypocenter relocation in heterogeneous media. For instance, during the 2021 M_S 6.4 Yangbi earthquake in China, eikonal-based inversion of multiple arrivals (P, sPg, PmP) refined locations by incorporating 3D velocity variations, reducing uncertainties compared to 1D layered models.^[35] Additionally, velocity inversion using eikonal-derived finite-frequency kernels—such as those from phase-front tracking of surface waves—maps crustal structure by weighting sensitivities to off-great-circle paths. These applications highlight the eikonal's role in bridging kinematic modeling and structural imaging for seismic hazard assessment.^[35]

Computer Vision and Graphics

In computer vision, the eikonal equation facilitates the fast computation of distance fields on 2D and 3D voxel grids, which are essential for shape analysis by yielding Euclidean or weighted distances from arbitrary points to object boundaries. These fields, obtained by solving |\nabla d(\mathbf{x})| = 1/F(\mathbf{x}) where F denotes the local speed (often constant for Euclidean distances), enable applications such as silhouette extraction and surface reconstruction from point clouds, with errors as low as 3.89% compared to exact methods in 3D benchmarks like the Stanford Dragon model. The Fast Marching Method (FMM), a seminal discretization technique, achieves this in O(N \log N) time for N grid points by propagating a monotone front, supporting weighted variants for anisotropic media in shape modeling and visualization.^[36]^[37] For path planning and collision avoidance, the eikonal equation models optimal trajectories by treating obstacles as regions of zero speed (corresponding to infinite refractive index n = \infty), computing a travel-time field that guides agents around hazards via gradient descent on the solution. In robotic navigation, such as marine vessel routing, FMM-based solvers generate collision-free paths by expanding wavefronts from start points, incorporating dynamic obstacles through anisotropic speed functions that enforce safety distances (e.g., 200 m at the bow), with real-time replanning in complex environments like crowded ports. This approach avoids local minima by leveraging the eikonal's viscosity solution properties, ensuring globally optimal paths under spatial-temporal constraints.^[38]^[39] In computer graphics, the eikonal equation drives refraction simulation by deriving ray equations from |\nabla S| = n(\mathbf{x}), where n is the spatially varying refractive index, to trace light paths accurately in non-homogeneous media like glass or scattering fluids. The Eikonal Rendering technique, implemented on GPUs, approximates global illumination by adaptively marching wavefronts for caustics and internal reflections, achieving interactive frame rates (e.g., 24.8 fps for refractive wine glasses) while handling attenuation and anisotropy for realistic rendering of transparent objects. This method extends to broader light transport approximations, reducing computational overhead in scenes with complex refractive elements compared to full Monte Carlo ray tracing.^[40] The eikonal equation integrates seamlessly with level-set methods for evolving interfaces in image segmentation, where solutions provide signed distance functions to reinitialize the level-set \phi and propagate fronts at speeds modulated by image gradients (e.g., F = 1 - |\nabla I|). In applications like medical imaging, competing eikonal-derived distance fields from interior and exterior seeds delineate regions such as lymph nodes, enabling topology changes and robust boundary detection without explicit parameterization, often via fast marching for O(N \log N) efficiency. This front propagation framework, as in the Chan-Vese model, minimizes energy functionals for accurate segmentation of noisy data like MRI tumors.^[41]^[42]

Numerical Solution Techniques

Grid-Based Discretization Methods

Grid-based discretization methods approximate solutions to the eikonal equation on structured Cartesian grids by replacing the continuous gradient operator with finite difference approximations, typically ensuring consistency and stability through upwind schemes that respect the directional propagation of wavefronts. These methods discretize the domain into a uniform grid with spacing h in each direction, where the arrival time T is computed at grid points (x_i, y_j) = (ih, jh), and the slowness n(x) is evaluated at the same points. The core challenge is to approximate |\nabla T| in a way that maintains the monotonicity of the solution, avoiding oscillations and ensuring convergence to the physically relevant viscosity solution. In two dimensions, a basic upwind finite difference approximation uses backward differences to estimate the gradient components, assuming the wavefront advances from known values with smaller T. For instance, at grid point (i,j), the scheme solves

\sqrt{ \left( \frac{T_{i,j} - T_{i-1,j}}{h} \right)^2 + \left( \frac{T_{i,j} - T_{i,j-1}}{h} \right)^2 } = n_{i,j}

for T_{i,j}, which corresponds to an upwind stencil in the first quadrant and enforces causality by depending only on previously computed upwind values. This formulation arises from the continuous eikonal equation |\nabla T| = n(x), where the square root structure preserves the Euclidean norm of the gradient. More general upwind schemes consider all four quadrants by taking the minimum over possible stencil orientations to handle arbitrary propagation directions. Godunov's scheme enhances this by providing a monotone approximation through a Godunov-type flux that selects upwind differences based on the local gradient signs, ensuring the numerical Hamiltonian is non-decreasing and the solution remains causal with respect to wavefront propagation. Specifically, the scheme approximates the gradient as

|\nabla T_{i,j}| \approx \sqrt{ \left( \max\left( D_x^- T_{i,j}, 0 \right) - \min\left( D_x^+ T_{i,j}, 0 \right) \right)^2 + \left( \max\left( D_y^- T_{i,j}, 0 \right) - \min\left( D_y^+ T_{i,j}, 0 \right) \right)^2 },

where D_x^\pm and D_y^\pm are the backward and forward differences in the x and y directions, respectively, set equal to n_{i,j}. This construction guarantees the scheme is monotone, as the update for T_{i,j} takes the minimum of the current value and the solved value from the scheme, ensuring non-increasing arrival times and preventing non-physical values.^[43] Extensions to n dimensions generalize the upwind approximation by including squared differences along all coordinate axes, yielding

\sqrt{ \sum_{k=1}^n \left( \max\left( D_k^- T, 0 \right) - \min\left( D_k^+ T, 0 \right) \right)^2 } = n,

which maintains consistency for the full tensor product grid. For computational efficiency in high dimensions, dimensional splitting decomposes the operator into sequential one-dimensional solves along each axis, reducing the nonlinearity while preserving first-order accuracy. Error analysis for these schemes demonstrates first-order accuracy, with the local truncation error bounded by O(h) due to the one-sided differences, leading to global convergence at rate O(h) under uniform grid refinement. Monotonicity and consistency ensure that, as h \to 0, the numerical solution converges uniformly to the unique viscosity solution of the eikonal equation, as established by the Crandall-Lions theory for Hamilton-Jacobi equations.^[44]

Upwind and Monotone Schemes

Upwind and monotone schemes represent a class of efficient numerical methods for solving the eikonal equation, particularly suited for large-scale problems where causality and solution monotonicity must be preserved to ensure accurate front propagation without spurious oscillations. These schemes build on finite difference discretizations by enforcing an upwind ordering that respects the direction of information flow, akin to characteristics in hyperbolic equations, thereby achieving high speed and stability for the static Hamilton-Jacobi formulation of the eikonal.^[45] The fast marching method, introduced by Sethian in 1996, is a seminal Dijkstra-like algorithm that solves the eikonal equation by sorting grid points according to tentative arrival times T and updating neighboring points using the solved eikonal values in a single forward pass. This approach treats the discrete domain as a graph where each grid point is a node, and edges connect to upwind neighbors; a priority queue efficiently selects the next point with the smallest tentative T, freezing it once accepted to maintain causality. The method's efficiency stems from its O(N log N) complexity in d dimensions, where N is the number of grid points, making it particularly effective for 3D domains by avoiding iterative relaxations.^[46]^[45] Fast sweeping methods, developed by Zhao in 2005, provide an alternative iterative approach using the same upwind discretizations, such as Godunov schemes, combined with Gauss-Seidel relaxations and multiple directional sweeps across the grid to propagate information along characteristics. Unlike fast marching, fast sweeping does not require sorting and achieves practical O(N) complexity with a fixed number of sweeps (typically 4 in 2D, more in higher dimensions), making it simpler to parallelize and suitable for inhomogeneous media. Convergence is ensured by the monotonicity of the updates, where each grid point's value is non-increasing until the viscosity solution is reached.^[47] Ordered upwind methods extend the fast marching paradigm to more general static Hamilton-Jacobi equations, including anisotropic eikonal variants, by using structured priority queues to perform multi-dimensional updates in a causal order. Developed by Sethian and Vladimirsky in 2000, these methods divide the grid into ordered simplices and solve local systems along characteristic directions, ensuring single-pass computation with O(N log N) complexity even for non-isotropic speeds. This allows for accurate handling of complex geometries and varying refraction indices without repeated sweeps.^[48] A key feature of both fast marching and ordered upwind schemes is their monotonicity preservation, which guarantees convergence to the unique viscosity solution of the eikonal equation as a static Hamilton-Jacobi equation by ensuring that the numerical operator is non-decreasing in the solution values. This property prevents oscillations and ensures stability, as upwind differencing selects only known (upwind) information, aligning with the maximum principle for the eikonal. Such monotonicity is crucial for applications requiring precise travel-time fields, like seismic modeling.^[45]^[48] In terms of runtime comparisons for 3D domains, fast marching methods typically outperform direct sparse linear solvers, which discretize the eikonal into a nonlinear system solvable via methods like multigrid (O(N)) or direct factorization (O(N^{1.5})), especially for large N exceeding 10^6 points; this efficiency gap widens in heterogeneous media, where iterative solvers may require hundreds of relaxations, while fast marching achieves optimal single-pass performance.^[49]

References

[1]
[PDF] Math 5587 – Lecture 21
Nov 16, 2016 · This nonlinear partial differential equation is called the eikonal equation. The eikonal equation also shows up in wave propagation ...
[2]
[PDF] GEOMETRICAL OPTICS
Equation (1.1.5) is often called the eikonal equation or the ray-tracing equation. Many of the waves that we have already study satisfy this equation, such as ...
[3]
The eikonal equation - Book chapter - IOPscience - Institute of Physics
The eikonal equation is a partial differential equation with non-linearity found in wave propagation. It is an approximated version of the wave equation.
[4]
[PDF] Chapter 2 Geometrical optics - MIT OpenCourseWare
We now give a general solution of the eikonal equation, albeit in a somewhat implicit form, in terms of rays. The rays are the characteristic curves for the.
[5]
[PDF] Mathematical Morphology and Eikonal Equations on Graphs ... - HAL
Mar 3, 2009 · Our formulation of the eikonal leads to novel applications of this equation such as nonlocal image segmentation and weighted distances or data ...<|control11|><|separator|>
[6]
[PDF] A level set-based Eulerian approach for anisotropic wave propagation
For isotropic wave propagations, the eikonal equations, which are of Hamilton–Jacobi form, are usually solved using viscosity solution-based solvers over ...
[7]
[PDF] Wave propagation and ray methods - Arizona Math
WKB approximation u ∼ exp(i(ωt − θ(x))/)[u0(x, y) + u1(x, y) + ...] leads ... The eikonal equation. The equation |∇θ|2 = µ2(x) is nonlinear and ...
[8]
Principles of Optics. Second (Revised) Edition. : Born, Max, & Emil ...
Oct 17, 2023 · Principles of Optics. Second (Revised) Edition. by: Born, Max, & Emil Wolf,. Publication date: 1964-01-01. Publisher: Pergamon Press. Collection ...Missing: official | Show results with:official
[9]
[PDF] Lecture 10 Approximation methods: WKB method - TCM
In optics, WKB is known as eikonal method, and in general referred to by mathematicians as short wavelength asymptotics. In quantum mechanics, it provides ...
[10]
[PDF] The WKB approximation and Ray Theory - Biello
We will use the Eikonal equation to construct the ODEs for the rays and we will use the Amplitude equation to describe how the amplitude of the waves changes as ...
[11]
3.Fermat's Principle of Least Time - Galileo and Einstein
Fermat famously stated in the 1630's that a ray of light going from point A to point B always takes the route of least time.
[12]
The principle of least time - SEG Wiki
Apr 26, 2021 · Pierre de Fermat (1601–1665) formulated the rule known as Fermat's principle of least time. In his original statement, Fermat asserted that ...
[13]
Descartes on the Refraction and the Velocity of Light
Five months later, Clerselier dispatched two letters to Fermat in which the latter's assumption-the principle of least time-was attacked and torn to shreds. The ...<|separator|>
[14]
Detailed derivation of the generalized Snell–Descartes laws from ...
Mar 6, 2023 · Beginning with Fermat's principle, we provide a detailed derivation of the generalized laws of refraction and reflection for a geometry ...
[15]
[PDF] The Eikonal - Neo-classical physics
HEINRICH BRUNS. REGULAR MEMBER OF THE SAXON SOCIETY OF SCIENCE ... will use the term “eikonal” in order to have a concise expression. Each map ...
[16]
Applications of the eikonal approximation in quantum mechanical ...
May 1, 2023 · The term “eikonal” was coined by Heinrich Bruns in his 1895 manuscript “Das Eikonal” published in Leipzig by S. Hirzel, an English ...
[17]
(PDF) Eikonal Functions: Old and New - ResearchGate
Jan 1, 2025 · ... Bruns, unaware of Hamilton's work, proposed a similar. idea [Bruns, 1895]. Bruns also coined the name 'eikonal' to the same. function called ...
[18]
[PDF] On some Results of the View of a Characteristic Function in Optics ...
Hamilton divides mathematical optics into two principal parts: one part proposing to find in every particular case the form of the characteristic function V , ...Missing: eikonal | Show results with:eikonal
[19]
Hamilton's Method in Geometrical Optics - Optica Publishing Group
Hamilton's method. In 18325 Hamilton defined the characteristic function V of an optical system as the optical length of the ray joining ...Missing: primary source
[20]
Hamilton-Jacobi equation - Scholarpedia
Oct 21, 2011 · The Hamilton-Jacobi equation is used to generate particular canonical transformations that simplify the equations of motion.
[21]
[PDF] On the Geometry of the Hamilton-Jacobi Equation - ICMAT
With an emphasis on mechanics, Carl Gustav Jacob Jacobi (1804 − 1851) deepened. Hamilton's formulation, by clarifying some mathematical issues, and ...<|separator|>
[22]
Einstein and Hilbert: The Creation of General Relativity - arXiv
Apr 25, 2005 · Recent controversy, raised by a much publicized 1997 reading of Hilbert's proof-sheets of his article of November 1915, is also discussed.
[23]
[PDF] Einstein and Hilbert: Two Months in the History of General Relativity
In his postcard of November 7, 1915, EINSTEIN tells HILBERT that he is ... HILBERT-EINSTEIN equations, it seems to us proper to call the general theory.Missing: eikonal | Show results with:eikonal<|separator|>
[24]
Geometrical theory of optical imaging (IV) - Principles of Optics
Since S(r) satisfies the eikonal equation §3.1 (15), this function is fully specified by the refractive index function (r) alone, together with the appropriate ...
[25]
Eikonal Functions: Old and New - SpringerLink
Eikonal functions are among the oldest and most useful tools in optics. They form the foundations of geometrical optics and optical aberrations theory.
[26]
None
### Summary of Sections from Notes on Geometrical Optics
[27]
[PDF] Introduction to Partial Differential Equations, Math 463/513, Spring ...
Apr 10, 2015 · 1.7.1 Derivation of the Eikonal Equation . ... Equation (13.10) is called Snell's law. 13.3 Electromagnetic Waves in the Atmosphere. The ...
[28]
[PDF] The Wave Equation - MIT OpenCourseWare
Mar 3, 2007 · Fermat's Principle implies Snell's Law. The travel time curve, plotted as a function of offset, is typically a hyperbolic function. Near ...Missing: derivation | Show results with:derivation
[29]
[PDF] 1207.1.K.pdf - Caltech PMA
This equation is known in optics as the eikonal equation. It is formally the same as the. Hamilton-Jacobi equation of classical mechanics3 if we identify Ω ...
[30]
[PDF] Geodesic distances
The eikonal equation is simply a reflection of the fact that distance functions grow with constant speed radially away from the center point. It is a ...
[31]
[PDF] fast sweeping methods for static Hamilton-Jacobi equations
A special case of this type of equations is the eikonal equation: (3) ... An important application for (2) is obtaining geodesic distance on a manifold.
[32]
[PDF] Efficient Fast Marching with Finsler metrics. - HAL
Sep 28, 2012 · The Riemannian metric is equal to the euclidean norm, except on a band of width 1/100 along a spiraling curve Γ where it has eigenvalues 1 ...Missing: media | Show results with:media
[33]
viscosity solutions of hamilton-jacobi equations1
VISCOSITY SOLUTIONS OF HAMILTON-JACOBI EQUATIONS1. BY. MICHAEL G. CRANDALL AND PIERRE-LOUIS LIONS. Abstract. Problems involving Hamilton-Jacobi equations—which ...
[34]
[PDF] viscosity solutions of the eikonal equations - UChicago Math
The eikonal equation is a nonlinear PDE related to wave propagation. This paper will investigate how distance functions are viscosity solutions of eikonal ...
[35]
[PDF] An overview of static Hamilton-Jacobi equations - UCSB Math
This paper reviews some of the existence and uniqueness techniques for the time independent cases with a partic- ular interest in Eikonal-like Equations. In ...
[36]
https://www2.imm.dtu.dk/pubdb/edoc/imm4145.pdf
[37]
[PDF] 3D Distance Fields: A Survey of Techniques and Applications
The Eikonal equation states the obvious inverse relationship between the speed of the front and the gradient of the arrival time. Since F does not have to be ...
[38]
[PDF] A Schrödinger wave equation approach to the eikonal ... - UF CISE
The eikonal equation has several applications in image analysis, viz. signed distance functions for shape silhouettes, surface reconstruction from point clouds ...<|separator|>
[39]
Marine Applications of the Fast Marching Method - Frontiers
Jan 27, 2020 · In this paper, several solutions based on the Fast Marching Method are proposed. The basic method focus on collision avoidance and optimal planning.
[40]
A Path Planning Method for Ship Collision Avoidance Considering ...
This research proposes a novel path planning method based on the fast marching method to specifically assist with safe navigation for autonomous ships.A Path Planning Method For... · 3. A Proposed Approach · 3.3. 1. Trajectory...
[41]
[PDF] Eikonal Rendering: Efficient Light Transport in Refractive Objects
To summarize, we present a new fast and versatile framework de- rived from the eikonal equation that can jointly reproduce many lighting effects around complex ...
[42]
[PDF] Efficient Segmentation Based on Eikonal and Diffusion Equations
Segmentation uses Eikonal and diffusion equations to solve PDEs, computing distance functions for each region. A competition between these functions determines ...
[43]
[PDF] level set methods and their applications in image science
Level set methods are used in image science, especially for segmentation, processing level lines, and representing boundaries, often with PDE techniques.
[44]
A second-order distributed memory parallel fast sweeping method ...
Feb 1, 2023 · This is called a Godunov upwind difference scheme [4]. Using this discretization scheme in the Eikonal equation results in a large nonlinear ...
[45]
[PDF] Lecture notes on viscosity solutions
These notes are concerned with viscosity solutions for fully nonlinear equa- tions. A majority of the notes are concerned with Hamilton-Jacobi equations.
[46]
https://www.pnas.org/doi/10.1073/pnas.93.4.1591
[47]
A fast marching level set method for monotonically advancing fronts.
A fast marching level set method is presented for monotonically advancing fronts, which leads to an extremely fast scheme for solving the Eikonal equation.
[48]
Ordered upwind methods for static Hamilton–Jacobi equations | PNAS
We introduce a family of fast ordered upwind methods for approximating solutions to a wide class of static Hamilton–Jacobi equations with Dirichlet boundary ...
[49]
[PDF] arXiv:1403.1937v2 [math.NA] 8 Feb 2015
Feb 8, 2015 · Our approach is straightforward to implement (by deploying any sparse linear solver) and holds its own against contemporary fast marching ...