Operador Laplaciano en coordenadas curvilíneas

El operador Laplaciano en coordenadas cartesianas es:

$latex Delta u = frac{partial^2}{partial x^2} u + frac{partial^2}{partial y^2} u + frac{partial^2}{partial z^2} u = frac{partial^2}{(partial x^i)^2} u$.

¿Qué pasa cuando queremos expresarlo en otro sistema de coordenadas curvilineas tal que

$latex x=x(q^1,q^2,q^3) = x(q^i)$, $latex y=y(q^i)$, $latex z = z(q^i)$

cualesquiera? Pues despues de un poco de trabajo, se puede llegar a la expresión:

$latex boxed{frac{1}{h_1 h_2 h_3} frac{partial}{partial q^i} (frac{h_1 h_2 h_3}{h_i^2} frac{partial}{partial q^i} u)}$,

donde, si definimos el cambio de coordenadas

$latex boxed{phi(q^i) := (x(q^i),y(q^i),z(q^i))}$,

tenemos que

$latex boxed{h_i = |frac{partial}{partial q^i} phi|}$.

Para ver como funciona la formula, vamos a calcular el Laplaciano en coordenadas cilíndricas y en esféricas.

En el primer caso, tenemos que

$latex phi(r,theta,z) = (r cos theta, r sin theta, z)$ con

$latex r in mathbb{R}^+, theta in [0,2 pi], z in mathbb{R}$,

por lo que tenemos

$latex h_r = | frac{partial}{partial r} phi| = sqrt{(cos theta, sin theta, 0) cdot (cos theta, sin theta, 0)} = 1$,

$latex h_theta = |frac{partial}{partial theta} phi| = sqrt{(-r sin theta, r cos theta, 0) cdot (-r sin theta, r cos theta, 0)} = r$,

$latex h_z = |frac{partial}{partial z} phi| = sqrt{(0,0,1) cdot (0,0,1)} = 1$.

y entonces, al aplicar nuestra fórmula, obtenemos:

$latex frac{1}{r} bigg [ frac{partial}{partial r} ( r frac{partial}{partial r} u ) + frac{partial}{partial theta} ( frac{1}{r} frac{partial}{partial theta} u ) + frac{partial}{partial z} ( r frac{partial}{partial z} u ) bigg ] =$

$latex = frac{1}{r} (1 frac{partial}{partial r} + r frac{partial^2}{partial r^2}) u + frac{1}{r^2} frac{partial^2}{partial theta^2} u + frac{partial^2}{partial z^2} u =$

$latex = boxed{frac{partial^2}{partial r^2} u + frac{1}{r} frac{partial}{partial r} u + frac{1}{r^2} frac{partial^2}{partial theta^2} u + frac{partial^2}{partial z^2} u}$.

Ahora, en el caso de esféricas, tenemos:

$latex phi(r,theta,varphi) = (r sin theta cos varphi, r sin theta sin varphi, r cos theta)$ con

$latex r in mathbb{R}^+$, $latex theta in [0,pi]$ y $latex varphi in [0,2pi]$,

de manera que (el cuadrado hace referencia al producto escalar):

$latex h_r = |frac{partial}{partial r} phi| = sqrt{(sin theta cos varphi, sin theta sin varphi, cos theta)^2} = $

$latex = sqrt{sin^2 theta cos^2 varphi + sin^2 theta sin^2 varphi + cos^2 theta} = 1$

$latex h_theta = |frac{partial}{partial theta} phi| = sqrt{(r cos theta cos varphi, r cos theta sin varphi, – r sin theta)^2} = $

$latex = sqrt{r^2 cos^2 theta cos^2 varphi + r^2 cos^2 theta sin^2 varphi + r^2 sin^2 theta} = r$

$latex h_varphi = |frac{partial}{partial varphi} phi| = sqrt{(-r sin theta sin varphi,r sin theta cos varphi,0)^2} = $

$latex = sqrt{r^2 sin^2 theta sin^2 varphi + r^2 sin^2 theta cos^2 varphi + 0} = r sin theta$,

por lo que, con la fórmula, tenemos:

$latex frac{1}{r^2 sin theta} bigg [ frac{partial}{partial r} (r^2 sin theta frac{partial}{partial r}u) + frac{partial}{partial theta} (sin theta frac{partial}{partial theta}u) + frac{partial}{partial varphi} (frac{1}{sin theta} frac{partial}{partial varphi}u) bigg ] =$

$latex = frac{1}{r^2} frac{partial}{partial r} (r^2 frac{partial}{partial r}u) + frac{1}{r^2 sin theta} frac{partial}{partial theta} (sin theta frac{partial}{partial theta}u) + frac{1}{r^2 sin^2 theta} frac{partial^2}{partial varphi^2} u =$

$latex = frac{1}{r^2}(2r frac{partial}{partial r} + r^2 frac{partial^2}{partial r^2})u + frac{1}{r^2 sin theta}(cos theta frac{partial}{partial theta} + sin theta frac{partial^2}{partial theta^2})u + frac{1}{r^2 sin^2 theta} frac{partial^2}{partial varphi^2}u = $

$latex =boxed{ frac{partial^2}{partial r^2}u + frac{2}{r} frac{partial}{partial r}u + frac{1}{r^2} frac{partial^2}{partial theta^2}u + frac{cot theta}{r^2} frac{partial}{partial theta}u + frac{csc^2 theta}{r^2} frac{partial^2}{partial varphi^2}u}$

Únete a la conversación

2 comentarios

Dejar un comentario

Your email address will not be published. Required fields are marked *


¡IMPORTANTE! Responde a la pregunta: ¿Cuál es el valor de 8 6 ?