Frequency comb: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>Krishnavedala
m vector image
 
en>Fgnievinski
No edit summary
Line 1: Line 1:
Customer Service Manager Brzozowski from Carlisle, loves to spend time computers, como ganhar dinheiro na internet and actions. Loves to travel and have been motivated after going to Anjar.<br><br>Feel free to visit my web page; [http://www.comoganhardinheiro101.com ganhando dinheiro na internet]
{{Use dmy dates|date=July 2013}}
The '''eikonal equation''' (from German ''Eikonal'', which is from [[Ancient Greek|Greek]] εἰκών, image<ref>The Oxford English Dictionary. 2nd ed. 1989. OED Online. Oxford University Press. 4 April 2000 http://dictionary.oed.com/cgi/entry/00292404</ref><ref>Evans, L. C., ''Partial Differential Equations'', AMS Graduate Texts in Mathematics, Vol. 19, pg. 93.</ref>) is a [[non-linear]] [[partial differential equation]] encountered in problems of [[wave propagation]], when the [[wave equation]] is approximated using the [[WKB theory]]. It is derivable from [[Maxwell's equations]] of electromagnetics, and provides a link between [[physical optics|physical (wave) optics]] and [[geometric optics|geometric (ray) optics]].
 
The eikonal equation is of the form
: <math> | \nabla u(x)|=F(x), \ x\in \Omega</math>
subject to <math>u|_{\partial \Omega}=0</math>, where <math>\Omega</math> is an open set in <math>\mathbb{R}^n</math> with [[well-behaved]] boundary, <math>F(x)</math> is a function with positive values, <math>\nabla</math>  denotes the [[gradient]] and |·| is the [[Euclidean_norm#Euclidean_norm|Euclidean norm]]. Here, the right-hand side <math>F(x)</math> is typically supplied as known input. Physically, the solution <math>u(x)</math> is the shortest time needed to travel from the [[Boundary (topology)|boundary]] <math>\partial \Omega</math> to <math>x</math> inside <math>\Omega,</math> with <math>F(x)</math> being the time cost (not speed) at <math>x</math>.
 
In the special case when <math>F=1</math>, the solution gives the [[signed distance]] from  <math>\partial \Omega</math>.
 
One fast computational algorithm to approximate the solution to the eikonal equation is the [[fast marching method]].
 
==Physical interpretation==
The physical meaning of the eikonal equation is related to the formula
 
: <math> E = -\nabla V </math>
 
where <math>E</math> is the electric field strength and <math>V</math> 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 in the region is pushed to move at right angles to the lines of constant potential, and along lines of force determined by the field of the E vector and the sign of the charge.  Corresponding variables occur in fluid flow and thermodynamics.  Ray optics and electromagnetism 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.
 
==Mathematical description==
An eikonal equation is one of the form
:<math>H(x,\nabla u(x)) = 0</math>
:<math>u(0,x') = u_0(x'),\text{ for } x = (x_1,x')</math>
The plane <math>x = (0,x')</math> can be thought of as the initial condition, by thinking of <math>x_1</math> as <math>t.</math>  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 <math>c(x)^2 |\nabla_x u(x,t)|^2 = |\partial_t u(x,t)|^2</math>.  There it is an equation describing the phase fronts of waves.  Under reasonable hypothesis on the "initial" data, the eikonal equation admits a local solution, but a global solution (e.g. a solution for all time in the geometrical optics case) is not possible.  The reason is that [[caustic (mathematics)|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 <math>\partial_{x_0} H(x,\nabla u(x)) \neq 0</math> for <math>x = (0,x').</math>  We must also assume <math>H(x,\nabla u(x)) = 0</math>, for <math>x = (0,x').</math>
 
First, solve the problem <math>H(x,\xi(x)) = 0</math>, <math>\xi(x) = \nabla u(x), x\in H</math>.  This is done by defining curves (and values of <math>\xi</math> on those curves) as
:<math>\dot x(s) = \nabla_\xi H(x(s),\xi(s)), \;\;\;\; \dot \xi(s) = -\nabla_x H(x(s),\xi(s)).</math>
:<math>x(0) = x_0, \;\;\;\; \xi(x(0)) = \nabla u(x(0)).</math>  Note that even before we have a solution <math>u</math>, we know <math>\nabla u(x)</math> for <math>x = (0,x')</math> due to our equation for <math>H</math>.
That these equations have a solution for some interval <math>0 \leq s < s_1</math> follows from standard ODE theorems (using the non-characteristic hypothesis).  These curves fill out an [[open set]] around the plane <math>x = (0,x')</math>.  Thus the curves define the value of <math>\xi</math> in an open set about our initial plane.  Once defined as such it is easy to see using the chain rule that <math>\partial_s H(x(s), \xi(s)) = 0</math>, and therefore <math>H = 0</math> along these curves.
 
We want our solution <math>u</math> to satisfy <math>\nabla u = \xi</math>, or more specifically, for every <math>s</math>, <math>(\nabla u)(x(s)) = \xi(x(s)).</math>  Assuming for a minute that this is possible, for any solution <math>u(x)</math> we must have
:<math>\frac{d}{d s} u(x(s)) = \nabla u(x(s)) \cdot \dot x(s) = \xi \cdot \frac{\partial H}{\partial \xi}</math>,
and therefore
:<math>u(x(t)) = u(x(0)) + \int_0^t \xi(x(s))\cdot \dot x(s)\, ds.</math>
In other words, the solution <math>u</math> will be given in a neighborhood of the initial plane by an explicit equation.  However, since the different paths <math>x(t)</math>, 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 <math>u</math> is a solution)
:<math>\xi(x(t)) = \xi(x(0)) - \int_0^s \nabla_x H(x(s),\xi(x(s))).</math>
 
It remains to show that <math>\xi</math>, which we have defined in a neighborhood of our initial plane, is the gradient of some function <math>u</math>. This will follow if we show that the vector field <math>\xi</math> is curl free.  Consider the first term in the definition of <math>\xi</math>.  This term, <math>\xi(x(0)) = \nabla u(x(0))</math> is curl free as it is the gradient of a function.  As for the other term, we note
:<math>\frac{\partial^2}{\partial x_k \, \partial x_j} H = \frac{\partial^2}{\partial x_j \, \partial x_k} H. </math>
The result follows
 
==Applications==
*A concrete application is the [[Computation of radiowave attenuation in the atmosphere]].
*Finding the [[Shape from Shading]] in computer vision.
 
==See also==
* [[Hamilton–Jacobi equation]]
* [[Fermat's principle]]
 
==References==
* Paris, D. T. and Hurd F.K., ''Basic Electromagnetic Theory'', McGraw-Hill 1969, pg. 383–385.
* Arnold, V. I., ''Lectures on Partial Differential Equations'', Springer 2004, 2nd Edition, pg. 2–3.
 
==External links==
*[http://www.reproducibility.org/RSF/book/sep/lineiko/paper_html/node2.html The linearized 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
*[http://neo-classical-physics.info/uploads/3/0/6/5/3065888/bruns_-_the_eikonal.pdf  English translation of "Das Eikonal" by Heinrich Bruns]
 
==Notes==
{{Reflist}}
 
{{DEFAULTSORT:Eikonal Equation}}
[[Category:Partial differential equations]]

Revision as of 23:03, 29 October 2013

30 year-old Entertainer or Range Artist Wesley from Drumheller, really loves vehicle, property developers properties for sale in singapore singapore and horse racing. Finds inspiration by traveling to Works of Antoni Gaudí. The eikonal equation (from German Eikonal, which is from Greek εἰκών, image[1][2]) is a non-linear partial differential equation encountered in problems of wave propagation, when the wave equation is approximated using the WKB theory. It is derivable from Maxwell's equations of electromagnetics, and provides a link between physical (wave) optics and geometric (ray) optics.

The eikonal equation is of the form

|u(x)|=F(x),xΩ

subject to u|Ω=0, where Ω is an open set in n with well-behaved boundary, F(x) is a function with positive values, 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 Ω to x inside Ω, with F(x) being the time cost (not speed) at x.

In the special case when F=1, the solution gives the signed distance from Ω.

One fast computational algorithm to approximate the solution to the eikonal equation is the fast marching method.

Physical interpretation

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

E=V

where E is the electric field strength and V 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 in the region is pushed to move at right angles to the lines of constant potential, and along lines of force determined by the field of the E vector and the sign of the charge. Corresponding variables occur in fluid flow and thermodynamics. Ray optics and electromagnetism 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.

Mathematical description

An eikonal equation is one of the form

H(x,u(x))=0
u(0,x)=u0(x), for x=(x1,x)

The plane x=(0,x) can be thought of as the initial condition, by thinking of x1 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|xu(x,t)|2=|tu(x,t)|2. There it is an equation describing the phase fronts of waves. Under reasonable hypothesis on the "initial" data, the eikonal equation admits a local solution, but 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 x0H(x,u(x))0 for x=(0,x). We must also assume H(x,u(x))=0, for x=(0,x).

First, solve the problem H(x,ξ(x))=0, ξ(x)=u(x),xH. This is done by defining curves (and values of ξ on those curves) as

x˙(s)=ξH(x(s),ξ(s)),ξ˙(s)=xH(x(s),ξ(s)).
x(0)=x0,ξ(x(0))=u(x(0)). Note that even before we have a solution u, we know u(x) for x=(0,x) due to our equation for H.

That these equations have a solution for some interval 0s<s1 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 ξ in an open set about our initial plane. Once defined as such it is easy to see using the chain rule that sH(x(s),ξ(s))=0, and therefore H=0 along these curves.

We want our solution u to satisfy u=ξ, or more specifically, for every s, (u)(x(s))=ξ(x(s)). Assuming for a minute that this is possible, for any solution u(x) we must have

ddsu(x(s))=u(x(s))x˙(s)=ξHξ,

and therefore

u(x(t))=u(x(0))+0tξ(x(s))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)

ξ(x(t))=ξ(x(0))0sxH(x(s),ξ(x(s))).

It remains to show that ξ, 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 ξ is curl free. Consider the first term in the definition of ξ. This term, ξ(x(0))=u(x(0)) is curl free as it is the gradient of a function. As for the other term, we note

2xkxjH=2xjxkH.

The result follows

Applications

See also

References

  • Paris, D. T. and Hurd F.K., Basic Electromagnetic Theory, McGraw-Hill 1969, pg. 383–385.
  • Arnold, V. I., Lectures on Partial Differential Equations, Springer 2004, 2nd Edition, pg. 2–3.

External links

Notes

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

  1. The Oxford English Dictionary. 2nd ed. 1989. OED Online. Oxford University Press. 4 April 2000 http://dictionary.oed.com/cgi/entry/00292404
  2. Evans, L. C., Partial Differential Equations, AMS Graduate Texts in Mathematics, Vol. 19, pg. 93.