La teoría general de las soluciones de la ecuación de Laplace se conoce como teoría del potencial . Las soluciones dos veces continuamente diferenciables de la ecuación de Laplace son las funciones armónicas , [1] que son importantes en múltiples ramas de la física, en particular la electrostática, la gravitación y la dinámica de fluidos . En el estudio de la conducción de calor , la ecuación de Laplace es la ecuación del calor en estado estacionario . [2] En general, la ecuación de Laplace describe situaciones de equilibrio, o aquellas que no dependen explícitamente del tiempo.
Ecuación de Laplace en un anillo (radio interior r = 2 y radio exterior R = 4 ) con condiciones de contorno de Dirichlet u ( r =2) = 0 y u ( R =4) = 4 sin(5 θ )
El problema de Dirichlet para la ecuación de Laplace consiste en encontrar una solución φ en algún dominio D tal que φ en la frontera de D sea igual a alguna función dada. Dado que el operador de Laplace aparece en la ecuación del calor , una interpretación física de este problema es la siguiente: fijar la temperatura en el límite del dominio de acuerdo con la especificación dada de la condición de contorno. Permita que el calor fluya hasta que se alcance un estado estacionario en el que la temperatura en cada punto del dominio ya no cambie. La distribución de la temperatura en el interior vendrá dada entonces por la solución al correspondiente problema de Dirichlet.
Las condiciones de frontera de Neumann para la ecuación de Laplace no especifican la función φ en sí en la frontera de D sino su derivada normal . Físicamente, esto corresponde a la construcción de un potencial para un campo vectorial cuyo efecto se conoce sólo en la frontera de D. En el ejemplo de la ecuación del calor, equivale a prescribir el flujo de calor a través de la frontera. En particular, en una frontera adiabática, la derivada normal de φ es cero.
Las soluciones de la ecuación de Laplace se llaman funciones armónicas ; todos son analíticos dentro del dominio donde se satisface la ecuación. Si dos funciones cualesquiera son soluciones de la ecuación de Laplace (o cualquier ecuación diferencial lineal homogénea), su suma (o cualquier combinación lineal) también es una solución. Esta propiedad, llamada principio de superposición , es muy útil. Por ejemplo, las soluciones a problemas complejos se pueden construir sumando soluciones simples.
En dos dimensiones
La ecuación de Laplace en dos variables independientes en coordenadas rectangulares tiene la forma
Funciones analíticas
Las partes real e imaginaria de una función analítica compleja satisfacen la ecuación de Laplace. Es decir, si z = x + iy , y si
Sin embargo, el ángulo θ tiene un solo valor solo en una región que no encierra el origen.
La estrecha conexión entre la ecuación de Laplace y las funciones analíticas implica que cualquier solución de la ecuación de Laplace tiene derivadas de todos los órdenes y puede expandirse en una serie de potencias , al menos dentro de un círculo que no encierre una singularidad. Esto contrasta marcadamente con las soluciones de la ecuación de onda , que generalmente tienen menos regularidad [ cita requerida ] .
Existe una conexión íntima entre las series de potencias y las series de Fourier . Si desarrollamos una función f en una serie de potencias dentro de un círculo de radio R , esto significa que
Sean las cantidades u y v las componentes horizontal y vertical del campo de velocidades de un flujo irrotacional, estable e incompresible en dos dimensiones. La condición de continuidad para un flujo incompresible es que
La ecuación de Laplace no cambia bajo una rotación de coordenadas y, por tanto, podemos esperar que se pueda obtener una solución fundamental entre soluciones que sólo dependen de la distancia r desde el punto fuente. Si elegimos que el volumen sea una bola de radio a alrededor del punto fuente, entonces el teorema de divergencia de Gauss implica que
Una función de Green es una solución fundamental que también satisface una condición adecuada en la frontera S de un volumen V. Por ejemplo,
Ahora bien, si u es cualquier solución de la ecuación de Poisson en V :
y u asume los valores límite g en S , entonces podemos aplicar la identidad de Green (una consecuencia del teorema de la divergencia) que establece que
Las notaciones u n y G n denotan derivadas normales en S . En vista de las condiciones satisfechas por u y G , este resultado se simplifica a
Así, la función de Green describe la influencia en ( x ′, y ′, z ′) de los datos f y g . Para el caso del interior de una esfera de radio a , la función de Green se puede obtener mediante una reflexión (Sommerfeld 1949): el punto fuente P a una distancia ρ del centro de la esfera se refleja a lo largo de su línea radial hasta una punto P' que está a una distancia
Tenga en cuenta que si P está dentro de la esfera, entonces P′ estará fuera de la esfera. La función de Green viene dada entonces por
Armónicos esféricos reales (Laplace) Y ℓ m para ℓ = 0, ..., 4 (de arriba a abajo) y m = 0, ..., ℓ (de izquierda a derecha). Los armónicos zonales, sectoriales y teselares se representan a lo largo de la columna de la izquierda, la diagonal principal y en otros lugares, respectivamente. (Los armónicos de orden negativo se mostrarían girados alrededor del eje z con respecto a los de orden positivo).
Considere el problema de encontrar soluciones de la forma f ( r , θ , φ ) = R ( r ) Y ( θ , φ ) . Por separación de variables , resultan dos ecuaciones diferenciales al imponer la ecuación de Laplace:
La segunda ecuación se puede simplificar suponiendo que Y tiene la forma Y ( θ , φ ) = Θ( θ ) Φ( φ ) . Aplicar nuevamente la separación de variables a la segunda ecuación da paso al par de ecuaciones diferenciales
para algún número m . A priori, m es una constante compleja, pero debido a que Φ debe ser una función periódica cuyo período divide uniformemente a 2 π , m es necesariamente un número entero y Φ es una combinación lineal de las exponenciales complejas e ± imφ . La función solución Y ( θ , φ ) es regular en los polos de la esfera, donde θ = 0, π . Imponer esta regularidad en la solución Θ de la segunda ecuación en los puntos límite del dominio es un problema de Sturm-Liouville que obliga al parámetro λ a tener la forma λ = ℓ ( ℓ + 1) para algún entero no negativo con ℓ ≥ | metro | ; esto también se explica a continuación en términos del momento angular orbital . Además, un cambio de variables t = cos θ transforma esta ecuación en la ecuación de Legendre , cuya solución es un múltiplo del polinomio de Legendre asociado P ℓ m (cos θ ) . Finalmente, la ecuación de R tiene soluciones de la forma R ( r ) = A r ℓ + B r − ℓ − 1 ; requerir que la solución sea regular en todo R 3 fuerza a B = 0 . [nota 2]
Aquí se supuso que la solución tenía la forma especial Y ( θ , φ ) = Θ( θ ) Φ( φ ) . Para un valor dado de ℓ , hay 2 ℓ + 1 soluciones independientes de esta forma, una para cada entero m con − ℓ ≤ m ≤ ℓ . Estas soluciones angulares son producto de funciones trigonométricas , aquí representadas como una exponencial compleja , y polinomios de Legendre asociados:
Aquí Y ℓ m se llama función armónica esférica de grado ℓ y orden m , P ℓ m es un polinomio de Legendre asociado , N es una constante de normalización y θ y φ representan colatitud y longitud, respectivamente. En particular, la colatitud θ , o ángulo polar, varía desde 0 en el Polo Norte, hasta π /2 en el Ecuador, y π en el Polo Sur, y la longitud φ , o azimut , puede asumir todos los valores con 0 ≤ φ. < 2π . Para un entero fijo ℓ , cada solución Y ( θ , φ ) del problema de valores propios
La solución general a la ecuación de Laplace en una bola centrada en el origen es una combinación lineal de las funciones armónicas esféricas multiplicadas por el factor de escala apropiado r ℓ ,
Para , en su lugar se eligen los armónicos sólidos con potencias negativas de . En ese caso, es necesario expandir la solución de regiones conocidas en la serie de Laurent (aproximadamente ), en lugar de la serie de Taylor (aproximadamente ), para hacer coincidir los términos y encontrar .
Electrostática
Sea el campo eléctrico, la densidad de carga eléctrica y la permitividad del espacio libre. Entonces la ley de Gauss para la electricidad (primera ecuación de Maxwell) en forma diferencial establece [5]
Ahora, el campo eléctrico se puede expresar como el gradiente negativo del potencial eléctrico ,
[5]
Al conectar esta relación a la ley de Gauss, obtenemos la ecuación de Poisson para la electricidad, [5]
En el caso particular de una región sin fuente, la ecuación de Poisson se reduce a la ecuación de Laplace para el potencial eléctrico. [5]
Si el potencial electrostático se especifica en el límite de una región , entonces se determina de forma única. Si está rodeado por un material conductor con una densidad de carga específica y si se conoce la carga total, entonces también es única. [6]
Un potencial que no satisface la ecuación de Laplace junto con la condición de contorno es un potencial electrostático no válido.
Gravitación
Sean el campo gravitacional, la densidad de masa y la constante gravitacional. Entonces la ley de Gauss para la gravitación en forma diferencial es [7]
El campo gravitacional es conservativo y, por tanto, puede expresarse como el gradiente negativo del potencial gravitacional:
Usando la forma diferencial de la ley de gravitación de Gauss, tenemos
[7]
En el espacio vacío, y tenemos
En la métrica de Schwarzschild
S. Persides [8] resolvió la ecuación de Laplace en el espaciotiempo de Schwarzschild en hipersuperficies de t constante . Usando las variables canónicas r , θ , φ la solución es
Aquí P l y Q l son funciones de Legendre de primer y segundo tipo, respectivamente, mientras que r s es el radio de Schwarzschild . El parámetro l es un número entero arbitrario no negativo.
^ El símbolo delta, Δ, también se usa comúnmente para representar un cambio finito en alguna cantidad, por ejemplo . Su uso para representar al laplaciano no debe confundirse con este uso.
^ Las aplicaciones físicas suelen tomar la solución que desaparece en el infinito, haciendo A = 0 . Esto no afecta la porción angular de los armónicos esféricos.
^ Zill, Dennis G y Michael R Cullen. Ecuaciones diferenciales con problemas de valores en la frontera . 8.ª edición / ed., Brooks/Cole, Cengage Learning, 2013. Capítulo 12: Problemas de valores en la frontera en coordenadas rectangulares. pag. 462. ISBN 978-1-111-82706-9 .
^ abc Griffiths, David J. Introducción a la electrodinámica . 4ª ed., Pearson, 2013. Portada interior. ISBN 978-1-108-42041-9 .
^ El enfoque de los armónicos esféricos adoptado aquí se encuentra en (Courant & Hilbert 1962, §V.8, §VII.5).
^ abcd Griffiths, David J. Introducción a la electrodinámica . 4ª ed., Pearson, 2013. Capítulo 2: Electrostática. pag. 83-4. ISBN 978-1-108-42041-9 .
^ Griffiths, David J. Introducción a la electrodinámica . 4ª ed., Pearson, 2013. Capítulo 3: Potenciales. pag. 119-121. ISBN 978-1-108-42041-9 .
^ ab Chicone, C.; Mashhoon, B. (20 de noviembre de 2011). "Gravedad no local: ecuación de Poisson modificada". Revista de Física Matemática . 53 (4): 042501. arXiv : 1111.4702 . doi : 10.1063/1.3702449. S2CID 118707082.
^ Perside, S. (1973). "Las ecuaciones de Laplace y Poisson en el espacio-tiempo de Schwarzschild". Revista de Análisis y Aplicaciones Matemáticas . 43 (3): 571–578. Código bibliográfico : 1973JMAA...43..571P. doi : 10.1016/0022-247X(73)90277-1 .
Sommerfeld, A. (1949). Ecuaciones diferenciales parciales en Física . Nueva York: Academic Press.
Zachmanoglou, CE; Thoe, Dale W. (1986). Introducción a las ecuaciones diferenciales parciales con aplicaciones . Nueva York: Dover. ISBN 9780486652511.
Polianina, AD (2002). Manual de ecuaciones diferenciales parciales lineales para ingenieros y científicos . Boca Ratón: Chapman & Hall/CRC Press. ISBN 978-1-58488-299-2.
Descubra cómo los problemas de valores en la frontera regidos por la ecuación de Laplace se pueden resolver numéricamente mediante el método del elemento en la frontera