Ecuación de Laplace

De Wikipedia, la enciclopedia libre
Saltar a: navegación, búsqueda
Pierre-Simon Laplace

En cálculo vectorial, la ecuación de Laplace es una ecuación en derivadas parciales de segundo orden de tipo elíptico, que recibe ese nombre en honor al físico y matemático Pierre-Simon Laplace.

Introducida por las necesidades de la mecánica newtoniana, la ecuación de Laplace aparece en muchas otras ramas de la física teórica como la astronomía, la electrostática, la mecánica de fluidos o la mecánica cuántica.

Definición[editar]

En tres dimensiones, el problema consiste en hallar funciones reales doblemente diferenciables, una función \scriptstyle u de variables reales x, y, y z, tal que

En coordenadas cartesianas


{\partial^2 u\over \partial x^2 } +
{\partial^2 u\over \partial y^2 } +
{\partial^2 u\over \partial z^2 } = 0.

En coordenadas cilíndricas,

  {1 \over \rho} {\partial \over \partial \rho}
  \left( \rho {\partial u \over \partial \rho} \right) 
+ {1 \over \rho^2} {\partial^2 u \over \partial \theta^2}
+ {\partial^2 u \over \partial z^2 } =0

En coordenadas esféricas,

  {1 \over r^2} {\partial \over \partial r}
  \left( r^2 {\partial u \over \partial r} \right) 
+ {1 \over r^2 \sin \varphi} {\partial \over \partial \varphi}
  \left( \sin \varphi {\partial u \over \partial \varphi} \right) 
+ {1 \over r^2 \sin^2 \varphi} {\partial^2 u \over \partial \theta^2} =0

Muchas veces se escribe de la siguiente manera:

\nabla^2 u = 0 \,

donde \nabla^2 es el operador de Laplace o "laplaciano"

que también se escribe como:

\nabla \cdot \nabla u = 0,

donde \scriptstyle \nabla \cdot es la divergencia, y \scriptstyle \nabla es el gradiente

o sino, algunas veces la notación puede ser:

\Delta u = 0,\,

donde \Delta también es el operador de Laplace.

Las soluciones de la ecuación de Laplace se denominan funciones armónicas.

Si del lado derecho de la igualdad se especifica una función, f(x, y, z), es decir, si la ecuación se escribe como:

\Delta u = f\,

entonces se tiene la "ecuación de Poisson", por lo que la ecuación de Laplace es un caso particular de esta. La ecuación de Laplace también es un caso particular de la ecuación de Helmholtz.

La ecuación de Laplace, así como también la ecuación de Poisson, son los ejemplos más simples de ecuaciones en derivadas parciales elípticas.


Condiciones de contorno o frontera[editar]

Ecuación de Laplace sobre una corona (r=2 y R=4) con condiciones de contorno de Dirichlet: u(r=2)=0 y u(r=4)=4sin(5*θ)

Problema de Dirichlet[editar]

El problema de Dirichlet para la ecuación de Laplace consiste de hallar una solución u en algún dominio D tal que u sobre su contorno o frontera \partial D es igual a una función determinada:

\begin{cases}
\triangle u = 0,\quad & x\in D\\
u=g, & x\in \partial D
\end{cases}

Como el operador de Laplace aparece en la ecuación del calor, una interpretación física de este problema es lo siguiente: fijar la temperatura sobre el contorno del dominio de acuerdo a una especificación determinada de la condición de contorno. La temperatura fluye hasta que alcanza un estado estacionario en el que dicha temperatura en cada punto del dominio no cambia más. La distribución de la temperatura en el interior será entonces la solución correspondiente al problema de Dirichlet.

Problema de Neumann[editar]

Las condiciones de contorno de Neumann para la ecuación de Laplace no especifica la función u en sí mismo sobre el contorno \partial D, pero sí su derivada normal. Físicamente, esto corresponde a la construcción de un potencial para un campo vectorial cuyo efecto es conocido en el contorno de D:

\begin{cases}
\triangle u = 0,\quad & x\in D\\
\frac{\partial u}{\partial\eta} = g,& x\in \partial D
\end{cases}

Las soluciones de la ecuación de Laplace son funciones armónicas; son todas analíticas dentro del dominio donde la ecuación se satisface. Si cualesquiera de dos funciones son soluciones a la ecuación de Laplace (o de cualquier ecuación diferencial homogénea), su suma (o cualquier combinación lineal) es también una solución. Esta propiedad, llamada principio de superposición, es muy útil, por ejemplo, las soluciones de problemas complejos pueden construirse simplemente sumando las soluciones determinadas e variables.

Ecuación de Laplace en dos dimensiones[editar]

La ecuación de Laplace en dos variables independientes:

\varphi_{xx} + \varphi_{yy} = 0.\,

Funciones analíticas[editar]

Las partes reales e imaginarias de un función analítica en los complejos satisfacen la ecuación de Laplace. Es decir, si z=x + i.y, y si

f(z) = u(x,y) + i.v(x,y),\,

entonces la condición necesaria para que f(z) sea analítica es que se satisfagan las ecuaciones de Cauchy-Riemann:

u_x = v_y, \quad v_x = -u_y.\,

donde ux es la primera derivada parcial de u con respecto a x.

Entonces

u_{yy} = (-v_x)_y = -(v_y)_x = -(u_x)_x.\,

Por lo tanto u satisface la ecuación de Laplace. Un cálculo similar demuestra que v también satisface la ecuación de Laplace.

A la inversa, dada una función armónica, es la parte real de una función analítica, f(z) (al menos localmente). Una forma de probarlo es:

f(z) = \varphi(x,y) + i \psi(x,y),\,

entonces las ecuaciones de Cauchy-Riemann se satisfacen:

\psi_x = -\varphi_y, \quad \psi_y = \varphi_x.\,

Esta relación no determina \psi, solo sus incrementos:

d \psi = -\varphi_y\, dx + \varphi_x\, dy.\,

La ecuación de Laplace para \varphi implica que la condición de integrabilidad para \psi se satisface:

\psi_{xy} = \psi_{yx},\,

y así \psi puede definirse con una integral de línea. La condición de integrabilidad y el teorema de Stokes implica que el valor de la integral de línea que conecta dos puntos es independiente del camino. El par de soluciones resultante de la ecuación de Laplace se denominan funciones armónicas conjugadas. Esta construcción solo es válida localmente ,o siempre que el camino no esté rodeando una singularidad. Por ejemplo, si r y θ son coordenadas polares y

\varphi = \log r, \,

entonces una función analítica correspondiente es

f(z) = \log z = \log r + i\theta. \,

Sin embargo, el ángulo θ es univaluada solo en una región que no incluye el origen.

La estrecha relación entre la ecuación de Laplace y las funciones analíticas establece que cualquier solución de la ecuación de Laplace tiene derivadas en todos los órdenes, y puede expandirse en series de potencias, al menos dentro de un círculo que no incluya una singularidad. Esto está en contraste con las soluciones de la ecuación de onda, que por lo general tiene menor regularidad.

Hay una íntima conexión entre las series de potencias y las series de Fourier. Si expandimos una función f en series de potencias dentro de un círculo de radio R, esto significa que

f(z) = \sum_{n=0}^\infty c_n z^n,\,

con coeficientes definidos adecuadamente cuyas parte real e imaginaria están dadas por:

c_n = a_n + i b_n.\,

Entonces

f(z) = \sum_{n=0}^\infty \left[ a_n r^n \cos n \theta - b_n r^n \sin n \theta\right] + i \sum_{n=1}^\infty \left[ a_n r^n \sin n\theta + b_n r^n \cos n \theta\right],\,

la cual es una serie de fourier de f.

Flujo de fluido[editar]

Sean las cantidades u y v las componentes horizontal y vertical del campo de velocidad del flujo incompresible estacionario e irrotacional en dos dimensiones. LA condición de que el flujo es incompresible es que

u_x + v_y=0,\,

y la condición de que el flujo es irrotatacional es que

\nabla \times \mathbf{V}=v_x - u_y =0. \,

Si definimos el diferencial de ψ como

d \psi = v\, dx - u\, dy,\,

entonces la condición de incompressibilidad es la condición de integrabilidad para este diferencial: la función resultante se llama función de corriente porque es constante a lo largo de las líneas de flujo. Las primeras derivadas de ψ son

\psi_x = v, \quad \psi_y=-u, \,

y la condición de irrotacionabilidad establece que ψ satisface la ecuación de Laplace. La función armónica φ que es el conjugado de ψ se denomina potencial de velocidad. Las ecuaciones de Cauchy-Riemann establecen que

\varphi_x=-u, \quad \varphi_y=-v. \,

Así, para cada función analítica corresponde a un flijo de fluido incompresible estacionario e irrotational en el plano. Las parte real es el potencial de velocidad, y la parte imaginaria es la función de corriente.

Electrostática[editar]

De acuerdo a las ecuaciones de Maxwell, un campo eléctrico (u,v) en un espacio de dos dimensiones que es independiente del tiempo satisface

\begin{cases}\nabla \times (u,v) = v_x -u_y =0,\,\\
\nabla \cdot (u,v) = \rho,\,\end{cases}

donde \rho es la densidad de carga. La primera ecuación de Maxwell es la condición de integrabilidad para el diferencial

d \varphi = -u\, dx -v\, dy,\,

así el potencial eléctrico φ puede construirse para satisfacer

\varphi_x = -u, \quad \varphi_y = -v.\,

La segunda ecuación de Maxwell establece que

\varphi_{xx} + \varphi_{yy} = -\rho,\,

la cual es la ecuación de Poisson.

Es importante observar que la ecuación de Laplace puede usarse en problemas de tres dimensiones en electroestática y flujo de fluido así como en dos dimensiones.

Ecuación de Laplace en tres dimensiones[editar]

Solución fundamental[editar]

Una solución fundamental de la ecuación de Laplace satisface:

 \Delta u = u_{xx} + u_{yy} + u_{zz} = -\delta(x-x',y-y',z-z'), \,

donde la función delta de Dirac \delta es una fuente unitaria concentrada en un punto  (x',\, y', \, z'). No es una función en sí, sin embargo puede pensarse como el límite de funciones cuya integral sobre todo el espacio es unitaria, y cuya región donde la función es distinta de cero es solo en un punto (ver solución débil). Es común elegir una convención de signos diferente para esta ecuación, esto se hace cuando se define la solución fundamental. Frecuentemente la elección de este signo es conveniente para trabajar con un -\Delta que es un operador positivo. Así la definición de la solución fundamental implica que, si el laplaciano de u es integrado sobre cualquier volumen que encierra el punto de la fuente, entonces

 \iiint_V \nabla \cdot \nabla u \, dV =-1. \,

La ecuación de Laplace no cambia bajo un cambio de coordenadas, y entonces podemos esperar que la solución fundamental puede obtenerse entre soluciones que solo dependen de la distancia r del punto de la fuente. Si elegimos el volumen de una bola de radio a alrededor del punto de la fuente, entonce por el teorema de la divergencia de Gauss:

 -1= \iiint_V \nabla \cdot \nabla u \, dV = \iint_S u_r \, dS = 4\pi a^2 u_r(a).\,

Entonces

  u_r(r) = -\frac{1}{4\pi r^2},\,

sobre una esfera de radio r que tiene como centro al punto de la fuente y por lo tanto

 u = \frac{1}{4\pi r}.\,

Un argumento similar muestra que en dos dimensiones:

 u = \frac{-\log r}{2\pi}. \,

Función de Green[editar]

Una función de Green es una solución fundamental que también satisface una condición adecuada de contorno S de un volumen V. Por ejemplo, G(x,y,z;x',y',z')\, satisface

\begin{cases} \nabla \cdot \nabla G = -\delta(x-x',y-y',z-z') \quad \hbox{en} \quad V, \,\\
G = 0 \quad \hbox{si} \quad (x,y,z) \quad \hbox{en} \quad S. \,\end{cases}

Ahora si u es cualquier solución de la ecuación de Poisson en V:

 \nabla \cdot \nabla u = -f, \,

y u toma valores de contorno g sobre S, entonces podemos aplicar la identidad de Green, (una consecuencia del teorema de la divergencia) el cual satisface

 \iiint_V \left[ G \, \nabla \cdot \nabla u - u \, \nabla \cdot \nabla G \right]\, dV  = \iiint_V \nabla \cdot \left[ G \nabla u - u \nabla G \right]\, dV = \iint_S \left[ G u_n -u G_n \right] \, dS. \,

Las notaciones un y Gn se refieren a derivadas normales a S. En vista de que las condiciones satisfacen u y G, este resultado simplifica a

 u(x',y',z') =  \iiint_V G f \, dV + \iint_S G_n g \, dS. \,

Así la función de Green describe la influencia en \scriptstyle (x',y',z')\, de f y g. Para el caso del interior de una esfera de radio a, la función de Green puede obtenerse por medio de la reflexión:[1] el punto de la fuente P a distancia ρ del centro de la esfera se refleja a lo largo de la línea radial al punto P' que es en una distancia

 \rho' = \frac{a^2}{\rho}. \,

Se observa que si P está dentro de la esfera, entonces P' estará fuera de la esfera. La función de Green está dada entonces por

 \frac{1}{4 \pi R} - \frac{a}{4 \pi \rho R'}, \,

donde R es la distancia al punto de la fuente P y R' es la distancia al punto reflejado P'. Una consecuencia de esta expresión para la función de Green es la fórmula integral de Poisson. Sea ρ, θ, y φ las componentes de coordenadas esféricas del punto de la fuente P. Aquí θ es el ángulo con el eje vertical, la cual es contraria a la notación matemática estadounidense, pero cumple con el estándar europeo y la práctica de la física . Entonces la solución de la ecuación de Laplace equation dentro de la esfera está dado por

 u(P) = \frac{1}{4\pi} a^3\left( 1 - \frac{\rho^2}{a^2} \right) \iint \frac{g(\theta',\varphi') \sin \varphi' \, d\theta' \, d\varphi'}{(a^2 + \rho^2 - 2 a \rho \cos \Theta)^{3/2} },  \,

donde

 \cos \Theta = \cos \varphi \cos \varphi' + \sin\varphi \sin\varphi'\cos(\theta -\theta'). \,

Una consecuencia simple de esta fórmula es que si u es una función armónica, enotnces el valor de u en el dentro de la esfera es el valor medio de los valores sobre la esfera. Esta propiedad de valor medio implica inmediatamente que funciones armónicas no constantes no pueden tomar su valor máximo en un punto interior.

Electroestática[editar]

En el espacio libre la ecuación de laplace de cualquier potencial electroestático debe ser igual a cero ya que \rho (densidad de carga volumétrica) es cero en el espacio libre.

A partir del gradiente del potencial se obtiene el campo eléctrico

\mathbf{E}=-\nabla V

Tomando la divergencia del campo eléctrico se obtiene la ecuación de Poisson, que relaciona el potencial eléctrico con la densidad de carga

\nabla^2V = -\frac{\rho}{\varepsilon_0}

En el caso particular del espacio libre (\rho=0) la ecuación de Poisson se reduce a la de Laplace.

Usando el teorema de la unicidad y mostrando que un potencial satisface la ecuación de Laplace (la segunda derivada de V debería ser cero en el espacio libre) y el potencial tiene los valores correctos en el contorno, el potencial entonces está unívocamente definido.

Un potencial que no satisface la ecuación de Laplace junto con la condición de contorno es un potencial electroestático inválido.

Véase también[editar]

Notas[editar]

  1. Sommerfeld, 1949

Referencias[editar]