2023年10月20日金曜日

円筒座標のベクトル解析

円筒座標系の基本ベクトルは,$\bm{e}_\rho, \bm{e}_\phi, \bm{e}_z$であり互いに直交している。
$\bm{e}_\rho = \bm{e}_x \cos\phi + \bm{e}_y \sin\phi $, $\bm{e}_\phi = -\bm{e}_x \sin\phi + \bm{e}_y \cos\phi$から,$\dfrac{\partial \bm{e}_\rho }{\partial \phi} = \bm{e}_\phi,\ \  \dfrac{\partial \bm{e}_\phi }{\partial \phi} = - \bm{e}_\rho$ が成り立つ。その他の基本ベクトルの各変数での微分はゼロ。

ここで,$\nabla = \bm{e}_\rho \dfrac{\partial}{\partial \rho} + \bm{e}_\phi \dfrac{1}{\rho} \dfrac{\partial}{\partial \phi}  + \bm{e}_z \dfrac{\partial}{\partial z} $であり,任意のベクトルは,$\bm{A}= \bm{e}_\rho A_\rho + \bm{e}_\phi A_\phi  + \bm{e}_z A_z$とかける。

したがって,円筒座標系での発散や回転は,この演算子とベクトルの内積や外積を機械的に計算すれば良い。ただし,微分演算が基本ベクトルに作用して現れる項があることだけ注意が必要となる。今回は,$\bm{e}_\rho,\ \bm{e}_\phi$を$\phi$で微分する項の存在に気をつける。

発散は,$\nabla\cdot\bm{A} = \Bigl(  \bm{e}_\rho \dfrac{\partial}{\partial \rho} + \bm{e}_\phi \dfrac{1}{r} \dfrac{\partial}{\partial \phi}  + \bm{e}_z \dfrac{\partial}{\partial z}\Bigr) \cdot \Bigl( \bm{e}_\rho A_\rho + \bm{e}_\phi A_\phi  + \bm{e}_z A_z \Bigr)$
余分の項は,$\bm{e}_\phi \dfrac{1}{\rho} \dfrac{\partial}{\partial \phi} \cdot  \bm{e}_\rho A_\rho = \bm{e}_\phi \dfrac{1}{\rho} \cdot  \bm{e}_\phi A_\rho = \dfrac{1}{\rho} A_\rho$である。
$\therefore \  \nabla\cdot\bm{A} = \dfrac{1}{\rho}\dfrac{\partial}{\partial \rho} \bigl(\rho A_\rho \bigr) + \dfrac{1}{\rho}\dfrac{\partial}{\partial \phi} A_\phi + \dfrac{\partial}{\partial z} A_z$

回転は,$\nabla \times \bm{A} = \Bigl(  \bm{e}_\rho \dfrac{\partial}{\partial \rho} + \bm{e}_\phi \dfrac{1}{\rho} \dfrac{\partial}{\partial \phi}  + \bm{e}_z \dfrac{\partial}{\partial z}\Bigr) \times \Bigl( \bm{e}_\rho A_\rho + \bm{e}_\phi A_\phi  + \bm{e}_z A_z \Bigr)$
余分の項は,$\bm{e}_\phi \dfrac{1}{\rho} \dfrac{\partial}{\partial \phi} \times  \bm{e}_\phi A_\phi = \bm{e}_\phi \dfrac{1}{\rho}  \times  (-\bm{e}_\rho ) A_\phi  =  \dfrac{1}{\rho} A_\phi \bm{e}_z$である。
$\therefore \  \nabla \times \bm{A} = \Bigl( \dfrac{1}{\rho} \dfrac{\partial A_z}{\partial\phi}- \dfrac{\partial A_\phi}{\partial z}\Bigr) \bm{e}_r + \Bigl(  \dfrac{\partial A_\rho}{\partial z}- \dfrac{\partial A_z}{\partial \rho}\Bigr) \bm{e}_\phi +  \dfrac{1}{\rho} \Bigl(  \dfrac{\partial (\rho A_\phi )}{\partial \rho}- \dfrac{\partial A_\rho}{\partial \phi}\Bigr) \bm{e}_z $

0 件のコメント: