The linearized eikonal equation stanford university. In 1834, william rowan hamilton published a paper in which he displayed an understanding between the mathematics particle motion and the propagation of linear waves in the eikonal approximation. It is possible to transform this equation into a system of firstorder ordinary. The methods of group analysis are applied to establish a classification of eikonal equations for anisotropic stationary media, g. The scheme relies on an upwind finite difference approximation to the gradient and a resulting causality relationship that lends itself to a dijkstralike programming approach. Computational study of fast methods for the eikonal equation. Lemma applied to eikonal equation eikonal equation. It is the basis for continuous versions of mathematical m.
In this appendix, i remind the reader how the eikonal equation is derived from the wave equation. A fast eikonal equation solver using the schrodinger wave equation. R2, 1 with various boundary conditions and choices of n index of refraction and f source term. In this paper we propose a novel numerical algorithm, which we call the fast iterative method. This equation, like the eikonal equation, is an equation in the first partial derivatives and is of second degree. We show that a scheme to solve the 2d eikonal equation by a finitedifference method can violate causality for moderate to large velocity contrasts v 2 v 1 2. The algorithm takes advantage of the upwind nature of the eikonal equation in order to update the value of each grid point by a number of times bounded by the number of its neighbors. In there, the authors use the approximation of low wavelengths high wavenumbers, to neglect some terms in the wave equation and arrive at the eikonal equation.
The eikonal equation, which was transformed to the ode system by the method of characteristics, is considered. Finitedifferencesolution of the eikonal equation for transversely isotropic media david w. The only assumption on the weak solutions is that the. The eikonal equation from german eikonal, which is from greek, image is a nonlinear partial differential equation encountered in problems of wave propagation, when the wave equation is approximated using the wkb theory. On the solution of the complex eikonal equation in. Eikonal equation eikon image greek consider the following solution to the wave equation. The solution of the eikonal equation may have singularities. The proposed method manages a list of active nodes and iteratively updates the solutions on those grid points until they. Alternatively, linearized ray equations can be used to obtain perturbed rays for which traveltimes can be found. Progress in nonlinear differential equations and their applications, vol 37.
The crucial idea is that each sweeping ordering follows a family of. It describes simd, shared memory, and distributed memory machine models. Nov 11, 2016 the paper present principles and derivation of the iterative method for solving the eikonal equation. On the solution of the complex eikonal equation in acoustic. A linearized eikonal equation is developed for transversely isotropic ti media with a vertical symmetry axis vti. It is stressed that both algorithms should be considered. Since there is the index in independent of z, we may assume the slope of phase change in z direction is linear. Pdf numerical solution of eikonal equation using finite. Abstract we introduce a new fast marching method for the eikonal equation. Journal of advances in mathematics vol 15 2018 issn. All forms of the eikonal equation are firstorder partial differential equations. Ordered line integral methods for solving the eikonal equation.
Here, we propose eikonet, a deep learning approach to solving the eikonal equation, which characterizes the firstarrivaltime field in heterogeneous 3d velocity structures. A highly accurate solution to the eikonal equation on cartesian domains. This paper presents a novel solver for the eikonal equation that is designed to run e. Fim, to solve eikonal equations on massively parallel architectures. I think i came up with a way but i wonder if there is a simpler way to argue this. Eikonal equation, a nonlinear partial differential equation encountered in problems of wave propagation. The theory of viscosity solutions for eikonal equations or more general hamiltonjacobi equations is welldeveloped when x is a euclidean space. Raytracing and eikonal solutions for lowfrequency wavefields.
In the papers 16,17,18,19,20,21,22, we studied the eikonal equation using the subgroup structure of the proper orthochronous group p 1, 4. It is linear with respect to perturbations in the horizontal velocity or the anisotropy parameter. In other words, the derivatives are of the first order, while the degree of the exponent is equal to 2. Abstract pdf 327 kb 2007 multistencils fast marching methods. To formulate a wellposed initialvalue problem on equation, it is sufficient to specify at some closed surface and to choose one of the two branches. A new front marching algorithm for solving the eikonal equation is presented. We would like to devise a similar algorithm for equation 2. Specifically, it is a firstorder and seconddegree partial differential equation. Pdf a new finite difference eikonal equation solver for. As an alternative, we present a finitedifference scheme in which the solution region progresses outward from an expanding wavefront rather than an expanding square, and therefore honors causality.
That is, its derivation does not involve spatial approximation of derivatives and thus does not require that the nodes be. Approximate solutions of the eikonal equation can also be obtained. It is a raytheoretical approximation to the scalar wave equation. A simple derivation of the eikonal and transport equations. This section will discuss the derivation and the application conditions of wkbj method together with the eikonal equation and. A generalized front marching algorithm for the solution of the eikonal. On classification of symmetry reductions for the eikonal equation. Finitedifference solution to the eikonal equation for. The method is an extension of other recently developed, isotropic finitedifference algorithms. In the derivation of the complex eikonal equation, we use the acoustic approximation, which is a useful tool in the. The eikonal equation, which defines the relationship between the phase of the optical wave. A method of solving the eikonal equation, in either flat or curved space times, with arbitrary cauchy data, is extended to the case of data given on a characteristic. The method is an iterative method which uses upwind difference for discretization and uses gaussseidel iterations with alternating sweeping ordering to solve the discretized system.
The eikonal approximation was not born in the study of quantum mechanics. This paper shows the application of geometric methods of electrodynamics to the calculation of optical devices, such as maxwell and luneburg lenses. By using nonconjugate subalgebras of ranks 1, 2, and 3, we performed the symmetry reduction of the equation 1 to differential equations with less a number of independent variables. An ecient solution to the eikonal equation on parametric. That is, the eikonal equation links neutron wave optics to neutron ray optics. From the form of the eikonal equation there results a remarkable analogy between geometrical optics and the mechanics of material particles.
Margrave abstract highfrequency approximations to the wave. Structure of entropy solutions to the eikonal equation ias school of. In general, the solution of the eikonal equation requires numerical methods. On the eikonal function cornbleet 1996 radio science. Unstructured grid solution approach for eikonal equation with. A characteristics driven fast marching method for the eikonal. From fermats principle, firstorder changes in the traveltime can be found by integrating along initial ray trajectories. The eikonal equation is a nonlinear partial differential equation encountered in problems of. A similar derivation is applicable for color images. The eikonal equation is a firstorder nonlinear partial differential equation. Eikonal approximation university of oklahoma physics. Blokhintsev 8 notes that the ray velocity is equal to the sum of the fluid. Eikonal is the german form of the greek word, meaning likeness, icon or image.
It is shown that the ansatz is a form of bateman potential, and the derivation explicitly defines rays, phase fronts, and the field polarization. Fo r a n m dimensional riemannian sp ace of c onstant curvature k with the metric ds 2. Singular sets regular sets solarity points eikonal equation. The eikonal equation is a nonlinear partial differential equation. Derivation of the eikonal equation physics stack exchange. Paraxial eikonal equations and applications william w. It describes the evolution of the wavefront, which is given by the equation. The eikonal equation, describing the traveltime propagation in an isotropic medium, has the form 1 where is the traveltime eikonal from the source to the point with the coordinates x, y, z, and n is the slowness at that point the velocity v equals 1 n. Homework statement a duck is walking in a flat desert yes, the proff. Eikonal equation, a nonlinear partial differential equation encountered in problems of wave propagation eikonal approximation, a method of approximation useful in wave scattering equations operation eikonal, an operation conducted by the national security agency in collaboration.
Nov, 2009 homework statement a duck is walking in a flat desert yes, the proff. Pdf iterative method for solving the eikonal equation. May 23, 2000 the fast marching method is a numerical algorithm for solving the eikonal equation on a rectangular orthogonal mesh in o m log m steps, where m is the total number of grid points. It is derivable from maxwells equations of electromagnetics, and provides a link between physical wave optics and geometric ray optics. In the eikonal equation, the rays of neutron optics are shown to be the trajectories of the particle. While the solution of the scalar wave equation represents the wavefield px, y, z. On classification of symmetry reductions for the eikonal. The application of the solutions of eikonal equation is numerous, e. According to atmospheric conditions the refractive index is ny1. Eikonal equations are fundamental to describe propagation of a wave front or interfaces in various disciplines of sciences and technology. Jun 22, 2019 thus, the eikonal equation relates the wave and geometric optics. Raytracing and eikonal solutions for lowfrequency wave. A highly accurate solution to the eikonal equation on cartesian domains a wide range of computer vision applications require an accurate solution of a particular hamiltonjacobi hj equation known as the eikonal equation. We show that a solution to the eikonal equation is obtained in the limit as 0 of the solution to the corresponding linear schrodinger.
The mathematical theory of geometrical optics can be regarded as the theory of the eikonal equation. For acoustics the eikonal equation will directly give wave fronts and thus the normal component of acoustic rays. Pdf on the solution of the complex eikonal equation in acoustic. We present a family of fast and accurate dijkstralike solvers for the eikonal equation and factored eikonal equation which compute solutions on a regular grid by solving local variational minimization problems. In terms of maxwells equations light is understood as a wave obeying a wave equation. Eikonal equation article about eikonal equation by the free.
However, finding solutions of maxwells equations is an extremely difficult task. Then, step 8 is used to obtain the final numerical solutions of the complex eikonal equation in. The derivation is classic and can be found in many popular textbooks. Then, step 8 is used to obtain the final numerical solutions of the complex eikonal equation in inhomogeneous vti media. Finitedifference solution of the eikonal equation along. A number of different numerical strategies have been proposed to efficiently solve the eikonal equation. The resulting system is written for the case of maxwell and luneburg lenses. Eaton abstract a new finitedifference technique is presented for solving the eikonal equation for inhomogeneous, transversely isotropic media. The fast marching method is a numerical algorithm for solving the eikonal equation on a rectangular orthogonal mesh in o m log m steps, where m is the total number of grid points. Fast methods for the eikonal and related hamilton jacobi. In principles of optics by max born and emil wolf, the eikonal equation is derived from maxwells equations. Here, we propose eikonet, a deep learning approach to solving the eikonal equation, which characterizes the firstarrivaltime field in heterogeneous 3d velocity. The eikonal approximation requires that the wavepacket is large compared to one wavelength, yet small compared with the distance one must travel before the coefficients a, b, c change significantly. The paper present principles and derivation of the iterative method for solving the eikonal equation.
In fact, steps 16 are used to solve the complex eikonal equation in the background medium, but step 7 is used to compute the coefficients of the taylor expansions 9 and 10. In this paper a fast sweeping method for computing the numerical solution of eikonal equations on a rectangular grid is presented. There was proli c activity in the application of eikonal approximation in high energy physics, especially in qcd. Over a century later, herbert goldstein would remark that hamilton would have postulated schrodingers equation had there been experimental evidence that particles were waves. A computational study of the fast marching and the fast sweeping methods for the eikonal equation is given. A new finite difference eikonal equation solver for anisotropic medium. The motion of a material particle is determined by the hamiltonjacobi equation 16. In 3d, we present two different families of algorithms which. Eikonal equation article about eikonal equation by the. Advanced classical mechanicsthe eikonal approximation and. In fact, the eikonal equation is a quasiclassical approximation the wentzelkramersbrillouin method of wave optics.
We use a schrodinger wave equation formalism to solve the eikonal equation. Our methods converge linearly but compute significantly more accurate solutions than competing first order methods. The eikonal equation and variants of it are of significant interest for problems in computer vision and image processing. In terms of maxwells equations light is understood as a. In order to overcome the problem of singularities and nonuniform grids arising when solving eikonal equation in spherical coordinate systems, a spherical cartesian coordinate system is defined and the hamiltonian form of the eikonal equation according to this coordinate system is given. Unstructured grid solution approach for eikonal equation. The recent deep learning revolution has created an enormous opportunity for accelerating compute capabilities in the context of physicsbased simulations. It is derivable from maxwells equations of electromagnetics, and provides a link between physical wave. The timefield tx again satisfies the eikonal equation, but now is complex with. Traveltime computation with the linearized eikonal equation. Nowack department of earth and atmospheric sciences, purdue university, west lafayette, in 47907, usa accepted 1992 january 15. A new method for solving the eikonal equation in the. Traveltime computation with the linearized eikonal. Alternatively, linearized ray equations can be used to obtain.
1337 14 419 99 1448 131 188 1297 1025 1119 420 228 401 149 503 1560 124 1260 1378 1426 274 845 1450 1510 1519 797 1486 1544 145 605 1440 1152 1499 102 1260 565 1373 133 1203 1296