Laplaciano

Operador laplaciano en coordenadas cartesianas y curvilíneas para ecuaciones de ondas y difusión.

Que es

El Laplaciano es un operador diferencial de segundo orden, denotado por $\nabla^2$ o $\Delta$, que surge de la divergencia del gradiente de un campo:

$\nabla^2 \phi = \vec{\nabla} \cdot (\vec{\nabla} \phi)$

En coordenadas cartesianas se expresa como la suma de las segundas derivadas parciales puras:

$\nabla^2 \phi = \frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} + \frac{\partial^2 \phi}{\partial z^2}$

Geométricamente y físicamente, el Laplaciano mide la curvatura espacial media o la discrepancia entre el valor central de una función en un punto y el promedio de la función en su entorno esférico inmediato:

$\nabla^2 \phi(P) = \lim_{R \to 0} \frac{6}{R^2} \left( \langle \phi \rangle_{S_R} - \phi(P) \right)$

• Si $\nabla^2 \phi(P) > 0$: El valor de $\phi(P)$ es menor que el promedio de sus vecinos (la función tiene concavidad positiva, como el fondo de un pozo).
• Si $\nabla^2 \phi(P) < 0$: El valor de $\phi(P)$ es mayor que el promedio de sus vecinos (la función tiene concavidad negativa, como la cima de una colina).
• Si $\nabla^2 \phi(P) = 0$: El valor coincide exactamente con el promedio de su entorno. Las funciones que satisfacen $\nabla^2 \phi = 0$ se denominan funciones armónicas.

Por qué aparece en física

El operador Laplaciano es omnipresente en todas las ramas de la física fundamental:

• Ecuación de Poisson: $\nabla^2 V = -\frac{\rho}{\varepsilon_0}$ (electrostática) y $\nabla^2 \Phi = 4\pi G \rho$ (gravitación). Relaciona la curvatura del potencial con la densidad de masa o carga.
• Ecuación de Laplace: $\nabla^2 \phi = 0$. Describe estados estacionarios de difusión, electrostática en el vacío y flujo potencial libre de vórtices.
• Ecuación de calor/difusión: $\frac{\partial T}{\partial t} = \alpha \nabla^2 T$. El calor fluye hacia donde la temperatura es inferior al promedio local.
• Ecuación de ondas: $\nabla^2 u = \frac{1}{c^2} \frac{\partial^2 u}{\partial t^2}$. Rige el sonido, la luz y las ondas elásticas.
• Ecuación de Schrödinger: $-\frac{\hbar^2}{2m} \nabla^2 \psi + V \psi = E \psi$. La energía cinética cuántica es proporcional al Laplaciano de la función de onda.

Idea visual o geométrica

Imagina una membrana de goma elástica tensada (como el parche de un tambor):

• Si la membrana está en reposo estático sin fuerzas externas que la empujen, su perfil de altura satisface la ecuación de Laplace $\nabla^2 z = 0$: cada punto está a la misma altura que el promedio de los puntos que lo rodean.
• Si colocas un peso pesado en un punto, la curvatura local cambia bruscamente ($\nabla^2 z \neq 0$), modelando la ecuación de Poisson con una fuente puntual.

Lenguaje simbólico mínimo

• Operador en cartesianas: $\nabla^2 = \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2}$.
• En cilíndricas $(r, \theta, z)$: $\nabla^2 \phi = \frac{1}{r} \frac{\partial}{\partial r}\left(r \frac{\partial \phi}{\partial r}\right) + \frac{1}{r^2} \frac{\partial^2 \phi}{\partial \theta^2} + \frac{\partial^2 \phi}{\partial z^2}$.
• En esféricas $(r, \theta, \varphi)$: $\nabla^2 \phi = \frac{1}{r^2} \frac{\partial}{\partial r}\left(r^2 \frac{\partial \phi}{\partial r}\right) + \frac{1}{r^2 \sin\theta} \frac{\partial}{\partial \theta}\left(\sin\theta \frac{\partial \phi}{\partial \theta}\right) + \frac{1}{r^2 \sin^2\theta} \frac{\partial^2 \phi}{\partial \varphi^2}$.
• Laplaciano vectorial: $\nabla^2 \vec{A} = \vec{\nabla}(\vec{\nabla}\cdot\vec{A}) - \vec{\nabla}\times(\vec{\nabla}\times\vec{A})$.

Que significa el resultado

Un valor no nulo de $\nabla^2 \phi$ indica que la distribución espacial de $\phi$ no está en equilibrio difusivo o armónico. En problemas de difusión, $\nabla^2 T > 0$ significa que el punto se calentará espontáneamente; $\nabla^2 T < 0$ significa que se enfriará.

Cuando no debes usarlo

• No apliques las fórmulas de derivadas segundas cartesianas directas en coordenadas curvilíneas. El Laplaciano en polares no es $\partial_r^2 + \partial_\theta^2$, sino $\frac{1}{r}\partial_r(r\partial_r) + \frac{1}{r^2}\partial_\theta^2$.
• Para campos vectoriales $\vec{A}$, no apliques el Laplaciano escalar componente a componente si las coordenadas no son cartesianas (los vectores unitarios $\hat{e}_r, \hat{e}_\theta$ cambian con la posición). Usa la identidad vectorial $\nabla^2 \vec{A} = \vec{\nabla}(\vec{\nabla}\cdot\vec{A}) - \vec{\nabla}\times(\vec{\nabla}\times\vec{A})$.

Resumen operativo

Para calcular el Laplaciano: calcula las derivadas parciales segundas puras de la función escalar y súmalas (en cartesianas) o aplica la fórmula métrica de Sturm-Liouville en cilíndricas/esféricas. Si $\nabla^2 \phi = 0$, la función satisface el teorema del valor medio y no posee extremos locales en el interior del dominio (Teorema de Earnshaw).

Procedimiento

Antes de calcular

Identifica si la función es un campo escalar $\phi$ o un campo vectorial $\vec{A}$, y elige el sistema de coordenadas más simétrico.

Paso 1: identifica el objeto matemático

Escribe la función analítica y determina si se busca verificar la ecuación de Laplace ($\nabla^2 \phi = 0$), Poisson ($\nabla^2 \phi = f$), calor o Helmholtz.

Paso 2: elige la operación

Selecciona la expresión del Laplaciano adaptada al sistema de coordenadas:

• Cartesianas: $\frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} + \frac{\partial^2 \phi}{\partial z^2}$.
• Cilíndricas: $\frac{1}{r}\partial_r(r\partial_r \phi) + \frac{1}{r^2}\partial_\theta^2 \phi + \partial_z^2 \phi$.
• Esféricas: $\frac{1}{r^2}\partial_r(r^2\partial_r \phi) + \frac{1}{r^2\sin\theta}\partial_\theta(\sin\theta\partial_\theta \phi) + \frac{1}{r^2\sin^2\theta}\partial_\varphi^2 \phi$.

Paso 3: conserva unidades y signos

Diferencia paso a paso: primero las derivadas interiores de primer orden, aplica los factores métricos multiplicativos, y efectúa la segunda derivada exterior.

Paso 4: calcula

Suma todos los términos y simplifica algebraicamente. Si se evalúa en un punto numérico, sustituye las coordenadas para obtener el valor escalar.

Paso 5: interpreta el resultado

• Si $\nabla^2 \phi = 0$: la función es armónica (equilibrio de Laplace).
• Si $\nabla^2 \phi \neq 0$: el valor es proporcional a la densidad de fuentes (Poisson) o a la tasa de variación temporal (difusión).

Checklist de autocontrol

• ¿El resultado para un campo escalar es un número escalar?
• ¿Se respetó el orden de derivación en coordenadas cilíndricas/esféricas?
• ¿Las unidades resultantes son $[\phi]/[\text{m}^2]$?

Ejemplo guiado

Situación física

Un cable coaxial cilíndrico de longitud infinita tiene un conductor interno de radio $a = 1.0 \text{ mm}$ mantenido a potencial $V_a = 100 \text{ V}$ y un conductor externo de radio $b = 5.0 \text{ mm}$ conectado a tierra ($V_b = 0 \text{ V}$). El espacio intermedio ($a < r < b$) es vacío libre de cargas libres. Se desea:

1. Resolver la ecuación de Laplace $\nabla^2 V = 0$ para encontrar el potencial $V(r)$.
2. Calcular el campo eléctrico $\vec{E}(r) = -\vec{\nabla} V$.

Traducción matemática

Por simetría cilíndrica y longitud infinita, el potencial depende únicamente del radio $r$ ($V = V(r)$). En coordenadas cilíndricas, el Laplaciano se reduce a: $\nabla^2 V = \frac{1}{r} \frac{d}{dr}\left( r \frac{dV}{dr} \right) = 0$

Calculo paso a paso

1. Integración de la ecuación diferencial: $\frac{d}{dr}\left( r \frac{dV}{dr} \right) = 0 \implies r \frac{dV}{dr} = C_1 \implies \frac{dV}{dr} = \frac{C_1}{r}$ Integrando nuevamente respecto a $r$: $V(r) = C_1 \ln(r) + C_2$

2. Aplicación de condiciones de contorno:

   • En $r = a$: $V(a) = C_1 \ln(a) + C_2 = V_a$
   • En $r = b$: $V(b) = C_1 \ln(b) + C_2 = 0 \implies C_2 = -C_1 \ln(b)$ Sustituyendo $C_2$: $C_1 \ln(a) - C_1 \ln(b) = V_a \implies C_1 \ln\left(\frac{a}{b}\right) = V_a \implies C_1 = \frac{V_a}{\ln(a/b)} = -\frac{V_a}{\ln(b/a)}$ Potencial analítico: $V(r) = V_a \frac{\ln(b/r)}{\ln(b/a)}$

3. Evaluación numérica con $a = 1\text{ mm}$, $b = 5\text{ mm}$, $V_a = 100\text{ V}$: $\ln(b/a) = \ln(5) \approx 1.60944$ $V(r) = 100 \frac{\ln(0.005/r)}{1.60944} \approx 62.133 \ln(0.005/r) \quad [\text{V}]$

4. Cálculo del campo eléctrico: $\vec{E}(r) = -\frac{dV}{dr} \hat{e}_r = -\left( \frac{C_1}{r} \right) \hat{e}_r = \frac{V_a}{r \ln(b/a)} \hat{e}_r = \frac{62.133}{r} \hat{e}_r \quad [\text{V/m}]$ En la superficie del conductor interno ($r = a = 0.001\text{ m}$): $E(a) = \frac{62.133}{0.001} \approx 62133 \text{ V/m} = 62.13 \text{ kV/m}$

Lectura del resultado

El potencial disminuye logarítmicamente entre los conductores, y el campo eléctrico decrece como $1/r$. El Laplaciano es idénticamente cero en todo el dieléctrico, lo que confirma que no hay acumulación de carga espacial entre las placas cilíndricas.

Error común evitado

Plantear $\frac{d^2 V}{dr^2} = 0$, lo que llevaría a un falso perfil lineal $V(r) = C_1 r + C_2$. En coordenadas cilíndricas, el factor métrico $r$ dentro de $\frac{1}{r}\frac{d}{dr}(r\frac{dV}{dr})$ es obligatorio.

Ejemplo inverso

Resultado dado

Una nube esférica de carga espacial produce un potencial eléctrico con Laplaciano uniforme en su interior ($r < R$): $\nabla^2 V = -K \quad (K > 0)$

Que pregunta física permite responder

¿Cuál es la distribución de densidad de carga $\rho$ en el interior de la nube? Por la ecuación de Poisson electrostática: $\nabla^2 V = -\frac{\rho}{\varepsilon_0}$ Igualando con el resultado dado: $-\frac{\rho}{\varepsilon_0} = -K \implies \rho = K \varepsilon_0$ El Laplaciano constante revela inmediatamente que la densidad de carga es estrictamente homogénea y constante en toda la esfera: $\rho(r) = \rho_0 = K \varepsilon_0$.

Errores frecuentes

Error 1: Confundir el Laplaciano con el gradiente al cuadrado

Interpretar erróneamente $\nabla^2 \phi$ como $\|\vec{\nabla}\phi\|^2$. El Laplaciano $\nabla^2 \phi = \vec{\nabla}\cdot(\vec{\nabla}\phi)$ es una suma de segundas derivadas (lineal), mientras que $\|\vec{\nabla}\phi\|^2 = (\partial_x \phi)^2 + (\partial_y \phi)^2 + (\partial_z \phi)^2$ es un producto no lineal de primeras derivadas.

Error 2: Omitir los factores métricos en coordenadas no cartesianas

Calcular el Laplaciano esférico como $\frac{\partial^2 \phi}{\partial r^2}$. El término radial correcto es $\frac{1}{r^2}\frac{\partial}{\partial r}(r^2\frac{\partial \phi}{\partial r}) = \frac{\partial^2 \phi}{\partial r^2} + \frac{2}{r}\frac{\partial \phi}{\partial r}$. Omitir el término $\frac{2}{r}\partial_r$ destruye la solución física.

Error 3: Aplicar el Laplaciano escalar a un vector componente a componente en coordenadas esféricas

Intentar calcular $\nabla^2 \vec{A}$ en esféricas calculando $\nabla^2 A_r \hat{e}_r + \nabla^2 A_\theta \hat{e}_\theta + \nabla^2 A_\varphi \hat{e}_\varphi$. Los vectores unitarios esféricos varían con el ángulo, por lo que este procedimiento es falso. Debe usarse la identidad $\nabla^2 \vec{A} = \vec{\nabla}(\vec{\nabla}\cdot\vec{A}) - \vec{\nabla}\times(\vec{\nabla}\times\vec{A})$.

Error 4: Suponer que una función armónica puede tener un máximo interior

Buscar trampas electrostáticas estables para partículas cargadas en el vacío. El Teorema de Earnshaw demuestra que al ser $\nabla^2 V = 0$, no pueden existir mínimos ni máximos de potencial locales estables en el vacío.

Error 5: Olvidar que el Laplaciano de una función lineal es cero

Sorprenderse de que $\nabla^2(a x + b y + c z + d) = 0$. Cualquier plano o función lineal satisface trivialmente la ecuación de Laplace.

Regla rápida de detección

Si la expresión de un Laplaciano contiene derivadas cruzadas del tipo $\partial_{xy}^2$ en cartesianas, o si las unidades no son $[\phi]/[\text{m}^2]$, el cálculo es incorrecto.

Aplicaciones físicas

1. Electrostática y gravitación: Ecuaciones de Poisson y Laplace

Gobierna los potenciales fundamentales de la naturaleza:

• $\nabla^2 V = -\rho/\varepsilon_0$ (potencial de Coulomb).
• $\nabla^2 \Phi = 4\pi G \rho$ (potencial gravitatorio newtoniano).
• $\nabla^2 V = 0$ en el vacío (Teorema de Earnshaw e imposibilidad de trampas electrostáticas puramente estáticas).

2. Mecánica cuántica no relativista: Ecuación de Schrödinger

La energía cinética del operador Hamiltoniano cuántico es $\hat{T} = -\frac{\hbar^2}{2m} \nabla^2$. El Laplaciano determina los orbitales atómicos del hidrógeno, bandas electrónicas en sólidos y niveles de osciladores armónicos.

3. Conducción térmica y procesos de difusión molecular

La ecuación del calor de Fourier $\partial_t T = \alpha \nabla^2 T$ y la segunda ley de Fick $\partial_t c = D \nabla^2 c$ describen cómo los gradientes espaciales se homogenizan irreversiblemente hacia el equilibrio termodinámico.

4. Acústica, óptica y radiación electromagnética

La ecuación de ondas d'Alembertiana $\nabla^2 u - \frac{1}{c^2} \partial_t^2 u = 0$ predice la propagación de ondas de sonido en fluidos, ondas sísmicas en la corteza terrestre y ondas de luz en el vacío.

5. Elasticidad y deformación de membranas y sólidos

En mecánica de sólidos, las tensiones de corte y deformaciones elásticas estáticas en placas delgadas y vigas bajo flexión se modelan mediante la ecuación biarmónica $\nabla^4 w = \nabla^2(\nabla^2 w) = q/D$.

Practica

Problema 1: Teorema del valor medio y funciones armónicas

Demuestre si la función potencial bidimensional $V(x,y) = e^x \cos(y) + 3xy$ es armónica calculando su Laplaciano $\nabla^2 V$.

Solución:

1. Derivadas respecto a $x$: $\frac{\partial V}{\partial x} = e^x \cos(y) + 3y \implies \frac{\partial^2 V}{\partial x^2} = e^x \cos(y)$
2. Derivadas respecto a $y$: $\frac{\partial V}{\partial y} = -e^x \sin(y) + 3x \implies \frac{\partial^2 V}{\partial y^2} = -e^x \cos(y)$
3. Sumando las derivadas segundas: $\nabla^2 V = \frac{\partial^2 V}{\partial x^2} + \frac{\partial^2 V}{\partial y^2} = e^x \cos(y) - e^x \cos(y) = 0$ La función satisface idénticamente la ecuación de Laplace en todo el plano $\mathbb{R}^2$; es una función armónica. Por el Teorema de Earnshaw, no posee máximos ni mínimos locales aislados en ninguna región acotada.

Problema 2: Ecuación de Poisson y densidad de carga en el átomo de hidrógeno

El potencial electrostático del estado fundamental del átomo de hidrógeno viene modelado por: $V(r) = \frac{e}{4\pi\varepsilon_0} \frac{e^{-2r/a_0}}{r}$ donde $a_0$ es el radio de Bohr. Utilizando la expresión del Laplaciano esférico con simetría radial, calcule $\nabla^2 V$ para $r > 0$ y determine la densidad de carga electrónica $\rho(r)$.

Solución:

1. Laplaciano radial esférico: $\nabla^2 V = \frac{1}{r} \frac{d^2}{dr^2}(r V)$. Sea $u(r) = r V(r) = \frac{e}{4\pi\varepsilon_0} e^{-2r/a_0}$.
2. Derivadas de $u(r)$: $\frac{du}{dr} = -\frac{2}{a_0} \frac{e}{4\pi\varepsilon_0} e^{-2r/a_0}$ $\frac{d^2 u}{dr^2} = \frac{4}{a_0^2} \frac{e}{4\pi\varepsilon_0} e^{-2r/a_0}$
3. Laplaciano: $\nabla^2 V = \frac{1}{r} \frac{d^2 u}{dr^2} = \frac{e}{\pi\varepsilon_0 a_0^2} \frac{e^{-2r/a_0}}{r}$
4. Densidad de carga por Poisson $\rho(r) = -\varepsilon_0 \nabla^2 V$: $\rho(r) = -\frac{e}{\pi a_0^3} \left(\frac{a_0}{r}\right) e^{-2r/a_0}$

Problema 3: Conducción térmica estacionaria en una esfera hueca

Una cáscara esférica de radios interior $R_1 = 0.10 \text{ m}$ y exterior $R_2 = 0.20 \text{ m}$ tiene temperaturas fijas $T(R_1) = 100^\circ\text{C}$ y $T(R_2) = 20^\circ\text{C}$. Sin fuentes internas, la temperatura satisface $\nabla^2 T = 0$. Determine el perfil radial $T(r)$ y calcule la temperatura a mitad de camino ($r = 0.15 \text{ m}$).

Solución:

1. Ecuación de Laplace radial: $\frac{1}{r^2} \frac{d}{dr}\left( r^2 \frac{dT}{dr} \right) = 0 \implies r^2 \frac{dT}{dr} = A \implies T(r) = -\frac{A}{r} + B$
2. Condiciones de contorno: $-\frac{A}{0.1} + B = 100, \quad -\frac{A}{0.2} + B = 20$ Restando: $-10A - (-5A) = 80 \implies -5A = 80 \implies A = -16 \text{ K}\cdot\text{m}$ $B = 20 + \frac{-16}{0.2} = 20 - 80 = -60^\circ\text{C}$ $T(r) = \frac{16}{r} - 60$
3. En $r = 0.15 \text{ m}$: $T(0.15) = \frac{16}{0.15} - 60 = 106.67 - 60 = 46.67^\circ\text{C}$