Eikonal equation

Eikonal equation

The eikonal equation is a non-linear partial differential equation of the form

: | abla u(x)|=F(x), xin Omega

subject to u|_{partial Omega}=0, where Omega is an open set in mathbb{R}^n with well-behaved boundary, F(x) is a function with positive values, abla denotes the gradient and |·| is the Euclidean norm. Here, the right-hand side F(x) is typically supplied as known input. Physically, the solution u(x) is the shortest time needed to travel from the boundary partial Omega to x inside Omega, with F(x) being the time cost (not speed) at x.

A fast computational algorithm to approximate the solution to the eikonal equation is the fast marching method. In the special case when F=1, the solution gives the signed distance from partial Omega.

The eikonal equation is encountered in problems of wave propagation, when the wave equation is approximated using the WKB theory.

The eikonal equation is derivable from Maxwell's equations of electromagnetics. It is the link between physical (wave) optics and geometric (ray) optics.

Physical interpretation

The physical meaning of the eikonal equation is related to the formula

: E = - abla Omega

where E is the electric field intensity and Omega is the electric potential. There is a similar equation for velocity potential in fluid flow and temperature in heat transfer. The physical meaning of this equation in the electromagnetic example is that any charge occurring in the region is pushed outward at a right angle from lines of constant potential and this charge travels along lines of constant force given by the field of the E vector. Corresponding variables occur in thermodynamics and fluid flow. Ray optics and electromagnetics are related by the fact that the eikonal equation gives a second electromagnetic formula of the same form as the potential equation above where the line of constant potential has been replaced by a line of constant phase and the force lines have been replaced by normal vectors coming out of the constant phase line at right angles. The magnitude of these normal vectors is given by the square root of the relative permittivity. The line of constant phase can be considered the edge of one of the advancing light waves. The normal vectors are the rays the light is traveling down in ray optics. This explanation is in the RMKS system of units used by electrical engineers.

Mathematical description

An eikonal equation is one of the form:H(x, abla u(x)) = 0:u(0,x') = u_0(x'), mbox{ for } x = (x_1,x')The plane x = (0,x') can be thought of as the initial condition, by thinking of x_1 as t. We could also solve the equation on a subset of this plane, or on a curved surface, with obvious modifications.This shows up in geometrical optics for example, where the equation is c(x)^2 | abla_x u(x,t)|^2 = |partial_t u(x,t)|^2. There it is an equation describing the phase fronts of waves. The lucky thing is that, under reasonable hypothesis on the "initial" data, the eikonal equation admits a local solution. Unfortunately, a global solution (e.g. a solution for all time in the geometrical optics case) is not possible. The reason is that caustics may develop. In the geometrical optics case, this means that wavefronts cross.

We can solve the eikonal equation using the method of characteristics. Note though that one must make the "non-characteristic" hypothesis partial_{x_0} H(x, abla u(x)) eq 0 for x = (0,x'). We must also clearly assume H(x, abla u(x)) = 0, for x = (0,x').

First, solve the problem H(x,xi(x)) = 0, xi(x) = abla u(x), xin H. This is done by defining curves (and values of xi on those curves) as:dot x(s) = abla_xi H(x(s),xi(s)), ;;;; dot xi(s) = - abla_x H(x(s),xi(s)).:x(0) = x_0, ;;;; xi(x(0)) = abla u(x(0)). Note that even before we have a solution u, we know abla u(x) for x = (0,x') due to our equation for H.That these equations have a solution for some interval 0 leq s < s_1 follows from standard ODE theorems (using the non-characteristic hypothesis). These curves fill out an open set around the plane x = (0,x'). Thus the curves define the value of xi in an open set about our initial plane. Once defined as such it is easy to see using the chain rule that partial_s H(x(s), xi(s)) = 0, and therefore H = 0 along these curves.

We want our solution u to satisfy abla u = xi, or more specifically, for every s, ( abla u)(x(s)) = xi(x(s)). Assuming for a minute that this is possible, for any solution u(x) we must have :frac{d}{d s} u(x(s)) = abla u(x(s)) cdot dot x(s) = xi cdot frac{partial H}{partial xi},and therefore:u(x(t)) = u(x(0)) + int_0^t xi(x(s))cdot dot x(s), ds.In other words, the solution u will be given in a neighborhood of the initial plane by an explicit equation. However, since the different paths x(t), starting from different initial points may cross, the solution may become multi-valued, at which point we have developed caustics. We also have (even before showing that u is a solution):xi(x(t)) = xi(x(0)) - int_0^s abla_x H(x(s),xi(x(s))).

It remains to show that xi, which we have defined in a neighborhood of our initial plane, is the gradient of some function u. This will follow if we show that the vector field xi is curl free. Consider the first term in the definition of xi. This term, xi(x(0)) = abla u(x(0)) is curl free as it is the gradient of a function. As for the other term, we note:frac{partial^2}{partial x_k partial x_j} H = frac{partial^2}{partial x_j partial x_k} H .The result follows

References

* Paris, D. T. and Hurd F.K., "Basic Electromagnetic Theory", McGraw-Hill 1969, pg. 383-385.

External links

* [http://sepwww.stanford.edu/public/docs/sep94/sergey3/paper_html/node2.html The linearized eikonal equation]
* [http://mathworld.wolfram.com/EikonalEquation.html Mathworld entry on eikonal equation]
*The eikonal equation was used for continuum crowd simulation by [http://grail.cs.washington.edu/projects/crowd-flows/ Treuille, Cooper, and Popović] at the University of Washington Animation Research Labs


Wikimedia Foundation. 2010.

Игры ⚽ Нужен реферат?

Look at other dictionaries:

  • Eikonal approximation — The eikonal approximation is a method of approximation useful in wave scattering equations within the realms of quantum mechanics, optics, quantum electrodynamics, and partial wave expansionInformal descriptionThe main advantage the eikonal… …   Wikipedia

  • eikonal — adjective Describing the properties of a class of partial differential equation encountered in descriptions of wave propagation …   Wiktionary

  • Ecuación de la Eikonal — Saltar a navegación, búsqueda La Ecuación Eikonal es una Ecuación en derivadas parciales con no linealidad encontrada en propagación de ondas, cuando la ecuacion de onda es aproximada usando la Teoría WKB. Esto es derivable desde las Ecuaciones… …   Wikipedia Español

  • Mild-slope equation — Simulation of wave penetration involving diffraction and refraction into Tedious Creek, Maryland, using CGWAVE (which solves the mild slope equation). In fluid dynamics, the mild slope equation describes the combined effects of diffraction and… …   Wikipedia

  • Hamilton–Jacobi equation — In physics, the Hamilton–Jacobi equation (HJE) is a reformulation of classical mechanics and, thus, equivalent to other formulations such as Newton s laws of motion, Lagrangian mechanics and Hamiltonian mechanics. The Hamilton–Jacobi equation is… …   Wikipedia

  • Computation of radiowave attenuation in the atmosphere — One of the causes of attenuation of radio propagation is the absorption by the atmosphere. There are many well known facts on the phenomenon and qualitative treatments in textbooks.[1] A document published by the International Telecommunication… …   Wikipedia

  • List of mathematics articles (E) — NOTOC E E₇ E (mathematical constant) E function E₈ lattice E₈ manifold E∞ operad E7½ E8 investigation tool Earley parser Early stopping Earnshaw s theorem Earth mover s distance East Journal on Approximations Eastern Arabic numerals Easton s… …   Wikipedia

  • Geometrical optics — As a mathematical study, geometrical optics emerges as a short wavelength limit for solutions to hyperbolic partial differential equations. For a less mathematical introduction, please see optics. In this short wavelength limit, it is possible to …   Wikipedia

  • Monge cone — In the mathematical theory of partial differential equations (PDE), the Monge cone is a geometrical object associated with a first order equation. It is named for Gaspard Monge. In two dimensions, let be a PDE for an unknown real valued function… …   Wikipedia

  • List of nonlinear partial differential equations — In mathematics and physics, nonlinear partial differential equations are (as their name suggests) partial differential equations with nonlinear terms. They describe many different physical systems, ranging from gravitation to fluid dynamics, and… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”