数式で独楽する

数式を使って楽しむブログです

[tex: ]

3次元円柱座標系の回転 ~ 外積のように導く

本稿では、
\begin{eqnarray}
x &=& r \cos \theta \\
y &=& r \sin \theta \\
z &=& z \tag{1}
\end{eqnarray}で表される3次元の円柱座標系$(r, \theta, z)$の回転について述べます。
ここでは、ベクトルの外積のように導きます。
極座標 - 数式で独楽する
ベクトルの外積 - 数式で独楽する
ベクトルの回転 - 数式で独楽する

円柱座標系のナブラとベクトルはそれぞれ
\begin{eqnarray}
\nabla &=& \boldsymbol{e}_r \, \frac{\partial}{\partial r} + \boldsymbol{e}_\theta \, \frac{1}{r} \frac{\partial}{\partial \theta} + \boldsymbol{k} \, \frac{\partial}{\partial z} \tag{2} \\
\boldsymbol{A} &=& A_r \,\boldsymbol{e}_r + A_\theta \, \boldsymbol{e}_\theta + A_z \, \boldsymbol{k} \tag{3}
\end{eqnarray}で表します。
3次元円柱座標系の勾配 - 数式で独楽する
3次元円柱座標系の発散 - 数式で独楽する

これより、直接、外積を求めるようにしていきます。ただし、単位ベクトルの偏微分外積に注意して計算していきます。
3次元円柱座標系の偏微分 - 数式で独楽する
\begin{eqnarray}
\nabla \times \boldsymbol{A}&=& \left( \boldsymbol{e}_r \, \frac{\partial}{\partial r} + \boldsymbol{e}_\theta \, \frac{1}{r} \frac{\partial}{\partial \theta} + \boldsymbol{k} \, \frac{\partial}{\partial z} \right) \times (A_r \, \boldsymbol{e}_r + A_\theta \, \boldsymbol{e}_\theta + A_z \, \boldsymbol{k} ) \\
&=& \left( \frac{1}{r} \frac{\partial A_z}{\partial \theta} - \frac{\partial A_\theta}{\partial z} \right) \, \boldsymbol{e}_r + \left( \frac{\partial A_r}{\partial z} - \frac{\partial A_z}{\partial r} \right) \, \boldsymbol{e}_\theta + \left[ \frac{1}{r} \frac{\partial}{\partial \theta} (r A_\theta) - \frac{1}{r} \frac{\partial A_r}{\partial \theta} \right] \boldsymbol{k}
\end{eqnarray}
こちらの方が計算量は少ないです。
f:id:toy1972:20200615223154p:plain:w400
f:id:toy1972:20200523064844p:plain:w300

f:id:toy1972:20200116231500g:plain:w300