Cálculo de divergencias y rotacionales
De Laplace
Contenido |
1 Enunciado
Para los campos vectoriales
calcule su divergencia y su rotacional, empleando en cada caso, coordenadas cartesianas, cilíndricas y esféricas. ¿Cuáles son irrotacionales y cuáles solenoidales?
2 Campo A
2.1 Divergencia
La divergencia, calculada en cartesianas, del vector de posición, es
![\nabla{\cdot}\mathbf{A} = \frac{\partial x}{\partial x}+\frac{\partial y}{\partial y}+\frac{\partial z}{\partial z} = 3](/wiki/images/math/1/3/a/13a88fd5baacc5594760e6a5e66ca4da.png)
Para este mismo campo, en cilíndricas, sustituyendo la expresión de dada en otro problema
![\nabla{\cdot}\mathbf{A} = \frac{1}{\rho}\frac{\partial \rho^2}{\partial \rho}+\frac{\partial z}{\partial z} =
\frac{2\rho}{\rho}+1 = 3](/wiki/images/math/4/e/3/4e32b4218114696b4f72105ba358f805.png)
y, en esféricas,
![\nabla{\cdot}\mathbf{A} = \frac{1}{r^2}\frac{\partial (r^3)}{\partial r}= 3](/wiki/images/math/b/2/c/b2c2c5d9211b7cef39097c817f103eb3.png)
2.2 Rotacional
Para el rotacional de este mismo campo, empleando coordenadas cartesianas
![\nabla\times\mathbf{A} = \left|\begin{matrix}
\mathbf{u}_{x} & \mathbf{u}_{y} & \mathbf{u}_{z} \\ & & \\
\displaystyle\frac{\partial \ }{\partial x} & \displaystyle\frac{\partial \ }{\partial y} & \displaystyle\frac{\partial \ }{\partial z} \\ && \\ x & y & z
\end{matrix}\right| = 0\mathbf{u}_{x}+0\mathbf{u}_{y}+0\mathbf{u}_{z}=\mathbf{0}](/wiki/images/math/1/1/8/118899cc72356054376e5a8b90759dc3.png)
en cilíndricas
![\nabla\times\mathbf{A} = \frac{1}{\rho}\left|\begin{matrix}\mathbf{u}_{\rho} & \rho\mathbf{u}_{\varphi} & \mathbf{u}_{z} \\ & & \\
\displaystyle\frac{\partial \ }{\partial \rho} &\displaystyle\frac{\partial \ }{\partial \varphi} & \displaystyle\frac{\partial \ }{\partial z} \\ && \\ \rho & 0 & z\end{matrix}\right|= \mathbf{0}](/wiki/images/math/a/e/7/ae7befbab33803f9ec30404158449c1f.png)
y en esféricas
![\nabla\times\mathbf{A} = \frac{1}{r^2\mathrm{sen}\,\theta}\left|\begin{matrix}\mathbf{u}_{r} & r\mathbf{u}_{\theta} & r\,\mathrm{sen}\,\theta\mathbf{u}_{\varphi} \\ & & \\
\displaystyle \frac{\partial \ }{\partial r} &\displaystyle\frac{\partial \ }{\partial \theta} & \displaystyle\frac{\partial \ }{\partial \varphi} \\ && \\ r & 0 & 0\end{matrix}\right| =
\mathbf{0}](/wiki/images/math/a/6/8/a68e146211ce3722e58d0a8d04681d55.png)
Naturalmente los resultados son los mismos independientemente del sistema empleado para calcularlos.
3 Campo B
3.1 Divergencia
Para el segundo campo, su divergencia, calculada en cartesianas,
![\nabla{\cdot}\mathbf{B} = \frac{\partial (-y)}{\partial x}+\frac{\partial x}{\partial y}+\frac{\partial 0}{\partial z}=0](/wiki/images/math/9/8/f/98f1b16d681e33bff68ffb1865508022.png)
En cilíndricas este campo se escribe
![\mathbf{B} = -\rho\,\mathrm{sen}\varphi\mathbf{u}_{x}+\rho\cos\varphi\mathbf{u}_{y}=\rho\mathbf{u}_{\varphi}](/wiki/images/math/5/9/b/59b5fc8367fad6d507a438bf66c73b73.png)
y la divergencia
![\nabla{\cdot}\mathbf{B} = 0+ \frac{1}{\rho}\frac{\partial \rho}{\partial \varphi}+0 = 0](/wiki/images/math/d/a/2/da22759c77e8c4313e7d5a43c2dc83be.png)
En esféricas el campo es
![\mathbf{B} = r \,\mathrm{sen}\,\theta\mathbf{u}_{\varphi}](/wiki/images/math/9/d/c/9dc7a52f94fb882866b206eba16ffb29.png)
y la divergencia
![\nabla{\cdot}\mathbf{B} = 0+0+\frac{1}{r\,\mathrm{sen}\,\theta}\frac{\partial \left(r\,\mathrm{sen}\,\theta\right)}{\partial \varphi} = 0](/wiki/images/math/d/a/e/dae0565627356f699d8b63fd0f9d1f85.png)
3.2 Rotacional
Para el rotacional, en cartesianas,
![\nabla\times\mathbf{B} = \left|\begin{matrix}
\mathbf{u}_{x} & \mathbf{u}_{y} & \mathbf{u}_{z} \\ & & \\
\displaystyle\frac{\partial \ }{\partial x} & \displaystyle\frac{\partial \ }{\partial y} & \displaystyle\frac{\partial \ }{\partial z} \\ && \\ -y & x & 0\end{matrix}\right|=
0\mathbf{u}_{x}+0\mathbf{u}_{y}+2\mathbf{u}_{z}=2\mathbf{u}_{z}](/wiki/images/math/8/7/b/87bb98d27ecf9fe1ebc232d85eda8f2d.png)
En cilíndricas
![\nabla\times\mathbf{B} = \frac{1}{\rho}\left|\begin{matrix}\mathbf{u}_{\rho} & \rho\mathbf{u}_{\varphi} & \mathbf{u}_{z} \\ & & \\
\displaystyle\frac{\partial \ }{\partial \rho} &\displaystyle\frac{\partial \ }{\partial \varphi} & \displaystyle\frac{\partial \ }{\partial z} \\ && \\ 0 & \rho^2 & 0\end{matrix}\right|=
2\mathbf{u}_{z}](/wiki/images/math/0/0/e/00ee3bff04a3ab357d7bf470ba939ec7.png)
y en esféricas
![\nabla\times\mathbf{B} = \frac{1}{r^2\mathrm{sen}\,\theta}\left|\begin{matrix}\mathbf{u}_{r} & r\mathbf{u}_{\theta} & r\,\mathrm{sen}\,\theta\mathbf{u}_{\varphi} \\ & & \\
\displaystyle \frac{\partial \ }{\partial r} &\displaystyle\frac{\partial \ }{\partial \theta} & \displaystyle\frac{\partial \ }{\partial \varphi} \\ && \\ 0 & 0 & r^2\mathrm{sen}^2\theta\end{matrix}\right|=
2\cos\theta\mathbf{u}_{r}-2\,\mathrm{sen}\,\theta\mathbf{u}_{\theta}](/wiki/images/math/d/b/4/db44b94ab2a38fc5edf024a05d2af9fa.png)
De nuevo el resultado es el mismo aunque, al estar expresado en base diferentes, parece formalmente distinto.
4 Campo C
4.1 Divergencia
Para el tercer campo, la divergencia en cartesianas
![\nabla{\cdot}\mathbf{C} = -1-1+2 = 0](/wiki/images/math/c/9/f/c9ff0e65e4ad977c5c8f148191192a68.png)
En cilíndricas, aplicando los resultados del problema de cálculo de gradientes
![\mathbf{C} = -\rho\mathbf{u}_{\rho}+2z\mathbf{u}_{z}](/wiki/images/math/2/6/9/2694580501c678234edc4c060c76219d.png)
y el cálculo de la divergencia da
![\nabla{\cdot}\mathbf{C} =
-\frac{1}{\rho}\frac{\partial
\left(\rho^2\right)}{\partial \rho}+\frac{\partial (2z)}{\partial z}= -2 + 2 = 0](/wiki/images/math/4/4/5/4451bec49b4fbf6ec325550ec03f8182.png)
y en esféricas
![\mathbf{C} = r(3\cos^2\theta-1)\mathbf{u}_{r}-3r\,\mathrm{sen}\,\theta\cos\theta\mathbf{u}_{\theta}](/wiki/images/math/f/b/e/fbe49a07097085615ef8e5d3455d56fa.png)
![\nabla{\cdot}\mathbf{C} = 3(3\cos^2\theta-1)-3(\mathrm{sen}^2\theta-2\cos^2\theta)=0](/wiki/images/math/5/2/c/52cb1906714bdae021c2a7389aa1f476.png)
4.2 Rotacional
Para el rotacional, en cartesianas
![\nabla\times\mathbf{C} = \left|\begin{matrix}
\mathbf{u}_{x} & \mathbf{u}_{y} & \mathbf{u}_{z} \\ & & \\
\displaystyle\frac{\partial \ }{\partial x} & \displaystyle\frac{\partial \ }{\partial y} & \displaystyle\frac{\partial \ }{\partial z} \\ && \\ -x & -y & 2z\end{matrix}\right|=
\mathbf{0}](/wiki/images/math/c/5/0/c5070dfac4b30068e9137b4ca9fb88ba.png)
En cilíndricas
![\nabla\times\mathbf{C} = \frac{1}{\rho}\left|\begin{matrix}\mathbf{u}_{\rho} & \rho\mathbf{u}_{\varphi} & \mathbf{u}_{z} \\ & & \\
\displaystyle\frac{\partial \ }{\partial \rho} &\displaystyle\frac{\partial \ }{\partial \varphi} & \displaystyle\frac{\partial \ }{\partial z} \\ && \\ -\rho & 0 & 2z\end{matrix}\right|=
\mathbf{0}](/wiki/images/math/2/b/8/2b873ac3e86f0933645918ebd1a668b8.png)
y en esféricas
![\nabla\times\mathbf{C} = \frac{1}{r^2\mathrm{sen}\,\theta}\left|\begin{matrix}\mathbf{u}_{r} & r\mathbf{u}_{\theta} & r\,\mathrm{sen}\,\theta\mathbf{u}_{\varphi} \\ & & \\
\displaystyle \frac{\partial \ }{\partial r} &\displaystyle\frac{\partial \ }{\partial \theta} & \displaystyle\frac{\partial \ }{\partial \varphi} \\ && \\ r(3\cos^2\theta-1) & -3r^2\mathrm{sen}\,\theta\cos\theta & 0\end{matrix}\right|=
\mathbf{0}](/wiki/images/math/f/e/8/fe8a3b7af05bf985c18004249c6d0175.png)
5 Campo D
Por último, para el campo
![\mathbf{D} =
\rho^2\cos\varphi\,\mathbf{u}_{\rho}+\rho^2\mathrm{sen}\,\varphi\,\mathbf{u}_{\varphi}](/wiki/images/math/2/6/d/26d26ec4ad1930ce6b95fed526230a32.png)
calculamos en primer lugar su divergencia y su rotacional en cilíndricas, ya que en estas coordenadas viene expresado el campo.
![\nabla{\cdot}\mathbf{D} = 3\rho\cos\varphi +\rho\cos\varphi = 4\rho\cos\varphi](/wiki/images/math/4/e/9/4e94780680d12af26bbf201a6e380e46.png)
![\nabla\times\mathbf{D} = \frac{1}{\rho}\left|\begin{matrix}\mathbf{u}_{\rho} & \rho\mathbf{u}_{\varphi} & \mathbf{u}_{z} \\ & & \\
\displaystyle\frac{\partial \ }{\partial \rho} &\displaystyle\frac{\partial \ }{\partial \varphi} & \displaystyle\frac{\partial \ }{\partial z} \\ && \\
\rho^2\cos\varphi & \rho^3\mathrm{sen}\varphi & 0\end{matrix}\right| = 4\rho\mathrm{sen}\varphi\mathbf{u}_{z}](/wiki/images/math/2/b/9/2b9b1271fbcb047a68556449dbd7e13d.png)
Para calcular estas cantidades en cartesianas, pasamos el campo a este sistema
![\mathbf{D} =
\rho^2\cos\varphi(\cos\varphi\mathbf{u}_{x}+\mathrm{sen}\,\varphi\mathbf{u}_{y})+\rho^2\mathrm{sen}\,\varphi\left(-\mathrm{sen}\,\varphi\mathbf{u}_{x}+\cos\varphi\mathbf{u}_{y}\right)
=\left(x^2-y^2\right)\mathbf{u}_{x}+2xy\mathbf{u}_{y}](/wiki/images/math/9/8/3/983f5546b68301006545e1ce2849f659.png)
y calculamos su divergencia
![\nabla{\cdot}\mathbf{D} = 2x+2x = 4x](/wiki/images/math/9/c/e/9ce3d031290db17526f4b0948e11eb5e.png)
y su rotacional
![\nabla\times\mathbf{D} = \left|\begin{matrix}
\mathbf{u}_{x} & \mathbf{u}_{y} & \mathbf{u}_{z} \\ & & \\
\displaystyle\frac{\partial \ }{\partial x} & \displaystyle\frac{\partial \ }{\partial y} & \displaystyle\frac{\partial \ }{\partial z} \\ && \\ x^2-y^2 & 2xy & 0\end{matrix}\right|= 4y\mathbf{u}_{z}](/wiki/images/math/a/0/2/a023b4f62554a900a75a81cb1754e1c6.png)
Para pasar a esféricas, primero expresamos en sus componentes cartesianas
![\mathbf{D}=(x^2-y^2)\mathbf{u}_{x}+2xy\mathbf{u}_{y}
=r^2\mathrm{sen}^2\theta\cos(2\varphi)\mathbf{u}_{x}+ r^2\mathrm{sen}^2\theta\,\mathrm{sen}(2\varphi)\mathbf{u}_{y}](/wiki/images/math/e/4/e/e4ef7c99f5c6a243f20ab44efdaa22db.png)
A continuación hallamos las diferentes componentes en esféricas
![D_r = \mathbf{D}{\cdot}\mathbf{u}_{r}=r^2\mathrm{sen}^3\theta\cos\varphi](/wiki/images/math/f/a/4/fa4980c90fb4305722d378aaa0864732.png)
![D_\theta = \mathbf{D}{\cdot}\mathbf{u}_{\theta}=r^2\mathrm{sen}^2\theta\cos\theta\cos\varphi](/wiki/images/math/1/c/f/1cf219de82098a8d8b626f5a7f3cf713.png)
![D_\varphi = \mathbf{D}{\cdot}\mathbf{u}_{\varphi}=r^2\mathrm{sen}^2\theta\,\mathrm{sen}\varphi](/wiki/images/math/e/9/7/e9746366cfb0de2ba3421480fe7ee45f.png)
Luego, calculamos los diferentes sumandos que constituyen la divergencia
![\frac{1}{r^2}\frac{\partial \ }{\partial r}\left(r^2 D_r\right)=4r\,\mathrm{sen}^3\theta\cos\varphi](/wiki/images/math/7/7/b/77bad807ea639a29f89955277c2df96b.png)
![\frac{1}{r\,\mathrm{sen}\,\theta}\frac{\partial \ }{\partial \varphi}(D_\varphi)=r\,\mathrm{sen}\,\theta\cos\varphi](/wiki/images/math/e/2/b/e2b37b4dcd9119c51f4e0fced8f9fddd.png)
y, por último, sumamos
![\nabla{\cdot}\mathbf{D}=4r\,\mathrm{sen}\,\theta\cos\varphi](/wiki/images/math/3/0/a/30adf0604349ce114d5593f6df6951cc.png)
Para el rotacional
![\nabla\times\mathbf{D} = \frac{1}{r^2\mathrm{sen}\,\theta}\left|\begin{matrix}\mathbf{u}_{r} & r\mathbf{u}_{\theta} & r\,\mathrm{sen}\,\theta\mathbf{u}_{\varphi} \\ & & \\
\displaystyle \frac{\partial \ }{\partial r} &\displaystyle\frac{\partial \ }{\partial \theta} & \displaystyle\frac{\partial \ }{\partial \varphi} \\ && \\
r^2\mathrm{sen}^3\theta\cos\varphi & r^3\mathrm{sen}^2\theta\cos\theta\cos\varphi &
r^3\mathrm{sen}^3\theta\,\mathrm{sen}\varphi\end{matrix}\right|=
4r\,\mathrm{sen}\,\theta\cos\theta\,\mathrm{sen}\,\varphi\mathbf{u}_{r}-4r\,\mathrm{sen}^2\theta\,\mathrm{sen}\,\varphi\mathbf{u}_{\theta}](/wiki/images/math/6/d/b/6dbe7c097c0ef4cef477471698722df3.png)
Teniendo en cuenta que
![y=\rho\,\mathrm{sen}\,\varphi=r\,\mathrm{sen}\theta\,\mathrm{sen}\varphi](/wiki/images/math/6/e/0/6e04e0be61a43392b862c4ac8ceb4ae4.png)
![\mathbf{u}_{z}=\cos\theta\mathbf{u}_{r}-\mathrm{sen}\,\theta\mathbf{u}_{\theta}](/wiki/images/math/c/2/6/c2608411f172476ccb3348b6eef04d33.png)
puede verse que los tres resultados son coincidentes.