next up previous
Next: Comentários Up: Cálculo Vetorial Prático Previous: Aplicações

Aplicações ``peso-pesado''

Aqui veremos resultados que são difíceis, ou muito trabalhosos, de obter por outros métodos. Em primeiro lugar, vamos mostrar (como aquecimento) que $div\;rot=0$.
\begin{displaymath}
div\;rot\;\vec{V}=\vec{\nabla}.(\vec{\nabla}\times\vec{V})=
...
...{ijk}
\partial_j V_k= \epsilon_{ijk} \partial_i \partial_j V_k
\end{displaymath} (39)

Ora,
\begin{displaymath}
\epsilon_{ijk} \partial_i \partial_j =0
\end{displaymath} (40)

pois $\epsilon _{ijk}$ é antissimétrico pela troca de índices $(ij)$, e $\partial_i \partial_j = \partial_j \partial_i$. Segue que
\begin{displaymath}
div\; rot\; \vec{V}=0
\end{displaymath} (41)

qualquer que seja $\vec{V}$.

Seja $\vec{E}$ um campo vetorial. Por exemplo, o campo elétrico. Na teoria de Maxwell precisamos calcular $rot\;rot\;\vec{E}$. Com o nosso método, isso é simples:

$\displaystyle \left(\vec{\nabla}\times(\vec{\nabla}\times\vec{E})\right)_i$ $\textstyle =$ $\displaystyle \epsilon_{ijk}\partial_j(rot\;\vec{E})_k$ (42)
  $\textstyle =$ $\displaystyle \epsilon_{ijk}\partial_j \epsilon_{klm}\partial_lE_m$  
  $\textstyle =$ $\displaystyle \epsilon_{ijk} \epsilon_{klm}\partial_j\partial_lE_m$  
  $\textstyle =$ $\displaystyle \epsilon_{kij}\epsilon_{klm}\partial_j\partial_lE_m$  
  $\textstyle =$ $\displaystyle (\delta_{il}\delta_{jm}-\delta_{im}\delta_{lj})\partial_j\partial_l
E_m$  
  $\textstyle =$ $\displaystyle \partial_m\partial_iE_m-\partial_j\partial_j E_i$  

Resta interpretar o resultado. A última linha pode ser escrita:
\begin{displaymath}
\partial_i\partial_mE_m-\partial_j\partial_j
E_i=\partial_i(...
...2E_i=\left(grad\;div\;\vec{E}
-\vec{\nabla}^2 \vec{E}\right)_i
\end{displaymath} (43)

Logo,
\begin{displaymath}
\vec{\nabla}\times(\vec{\nabla}\times\vec{E})=
grad\;div\;\vec{E}-\vec{\nabla}^2\vec{E}
\end{displaymath} (44)

Outra relação de grande importância na eletrodinâmica envolve o cálculo de $div(\vec{E}\times\vec{B})$.

$\displaystyle div(\vec{E}\times\vec{B})=\partial_i \epsilon_{ijk} E_j B_k$ $\textstyle =$ $\displaystyle \epsilon_{ijk} (\partial_i E_j)B_k+\epsilon_{ijk} E_j(\partial_i
B_k)$ (45)
  $\textstyle =$ $\displaystyle B_k \epsilon_{kij}\partial_i E_j - E_j
\epsilon_{jik}\partial_iB_k$  
  $\textstyle =$ $\displaystyle B_k(rot\;\vec{E})_k-E_j(rot\;\vec{B})_j$  
  $\textstyle =$ $\displaystyle \vec{B}.rot\;\vec{E}-\vec{E}.rot\;\vec{B}$  

Finalmente, vamos ao nosso tour de force: calcular $\vec{\nabla}.
(\vec{A}.\vec{B})$.

\begin{displaymath}
\left(\vec{\nabla}(\vec{A}.\vec{B})\right)_i=
\partial_i(A_j B_j)=(\partial_i A_j)B_j+A_j(\partial_i B_j)
\end{displaymath} (46)


\begin{displaymath}
=(\partial_i A_j-\partial_j A_i)B_j+(\partial_j A_i)B_j+
A_j(\partial_i B_j-\partial_j B_i)+A_j\partial_j B_i
\end{displaymath}

Temos, preliminarmente, que
$\displaystyle (\partial_j A_i)B_j$ $\textstyle =$ $\displaystyle B_j\partial_j A_i
=(\vec{B}.\vec{\nabla})A_i$ (47)
$\displaystyle (A_j\partial_j)B_i$ $\textstyle =$ $\displaystyle (\vec{A}.\vec{\nabla})B_i$ (48)

de maneira que, também preliminarmente,
\begin{displaymath}
\vec{\nabla}(\vec{A}.\vec{B})= (\vec{B}.\vec{\nabla})\vec{A}
+(\vec{A}.\vec{\nabla})\vec{B}+ ...
\end{displaymath} (49)

Para calcular os termos adicionais, notemos que
\begin{displaymath}
\partial_iA_j-\partial_j A_i=\epsilon_{kij}\epsilon_{klm}
\partial_l A_m
\end{displaymath} (50)

como o leitor, a esta altura, poderá facilmente verificar. Logo,
$\displaystyle (\partial_i A_j-\partial_j A_i)B_j=B_j \epsilon_{kij}\epsilon_{klm}
\partial_l A_m$ $\textstyle =$ $\displaystyle B_j \epsilon_{kij}(rot\; \vec{A})_k$ (51)
  $\textstyle =$ $\displaystyle \epsilon_{ijk} B_j(rot\vec{A})_k=(\vec{B}\times rot\;\vec{A})_i$  

Analogamente,
\begin{displaymath}
A_j(\partial_i B_j-\partial_j B_i)=(\vec{A}\times rot\;\vec{B})_i
\end{displaymath} (52)

Juntando estes termos à eq.(49), temos
\begin{displaymath}
\vec{\nabla}(\vec{A}.\vec{B})= (\vec{B}.\vec{\nabla})\vec{A}...
...B}+ (\vec{B}\times rot\;\vec{A})
+(\vec{A}\times rot\;\vec{B})
\end{displaymath} (53)


next up previous
Next: Comentários Up: Cálculo Vetorial Prático Previous: Aplicações
Henrique Fleming 2001-12-18