logo

曲線座標系でのスカラー関数の勾配 📂数理物理学

曲線座標系でのスカラー関数の勾配

概要

曲線座標系でスカラー関数のf=f(q1,q2,q3)f=f(q_{1},q_{2},q_{3})の勾配は次の通りだ。

∇f=1h1∂f∂q1q^1+1h2∂f∂q2q^2+1h3∂f∂q3q^3=∑i=131hi∂f∂qiq^i \nabla f= \frac{1}{h_{1}}\frac{ \partial f }{ \partial q_{1} } \hat{\mathbf{q}}_{1} + \frac{1}{h_{2}}\frac{ \partial f }{ \partial q _{2}}\hat{\mathbf{q}}_{2}+\frac{1}{h_{3}}\frac{ \partial f }{ \partial q_{3} } \hat{\mathbf{q}}_{3}=\sum \limits _{i=1} ^{3}\frac{1}{h_{i}}\frac{ \partial f}{ \partial q_{i}}\hat{\mathbf{q}}_{i}

hih_{i}はスケール因子だ。

公式

  • 直交座標系:

    h1=h2=h3=1 h_{1}=h_{2}=h_{3}=1

    ∇f=∂f∂xx^+∂f∂yy^+∂f∂zz^ \nabla f= \frac{\partial f}{\partial x}\mathbf{\hat{\mathbf{x}} }+ \frac{\partial f}{\partial y}\mathbf{\hat{\mathbf{y}}} + \frac{\partial f}{\partial z}\mathbf{\hat{\mathbf{z}}}

  • 円筒座標系:

    h1=1,h2=ρ,h3=1 h_{1}=1,\quad h_{2}=\rho,\quad h_{3}=1

    ∇f=∂f∂ρρ^+1ρ∂f∂ϕϕ^+∂f∂zz^ \nabla f = \frac{\partial f}{\partial \rho}\boldsymbol{\hat \rho} + \frac{1}{\rho}\frac{\partial f}{\partial \phi}\boldsymbol{\hat \phi} + \frac{\partial f}{\partial z}\mathbf{\hat{\mathbf{z}}}

  • 球座標系:

    h1=1,h2=r,h3=rsin⁡θ h_{1}=1,\quad h_{2}=r\quad, h_{3}=r\sin\theta

    ∇f=∂f∂rr^+1r∂f∂θθ^+1rsin⁡θ∂f∂ϕϕ^ \nabla f= \frac{\partial f}{\partial r} \mathbf{\hat r} + \frac{1}{r}\frac{\partial f}{\partial \theta} \boldsymbol{\hat \theta} + \frac{1}{r\sin\theta}\frac{\partial f}{\partial \phi}\boldsymbol{\hat \phi}

導出

三次元デカルト座標系で、以下の式を満たすa\mathbf{a}をffの勾配と名付け、∇f\nabla fと記されると定義した。

df=a⋅dr d f =\mathbf{a} \cdot d\mathbf{r}

任意の曲線座標系でもこのように定義する。ffの全微分は下の通りだ。

df=∂f∂q1dq1+∂f∂q2dq2+∂f∂q3dq3 d f = \frac{ \partial f}{ \partial q_{1} }dq_{1}+\frac{ \partial f}{ \partial q_{2}}dq_{2}+\frac{ \partial f}{ \partial q_{3}}dq_{3}

曲線座標系での位置ベクトルr\mathbf{r}の微小変化量は次の通りだ。

dr=h1dq1q^1+h2dq2q^2+h3dq3q^3 d\mathbf{r}=h_{1}dq_{1}\hat{\mathbf{q}}_{1}+h_{2}dq_{2}\hat{\mathbf{q}}_{2}+h_{3}dq_{3}\hat{\mathbf{q}}_{3}

これから、以下の式を満たすa\mathbf{a}を探求する。

df=a⋅dr \begin{equation} df=\mathbf{a} \cdot d\mathbf{r} \end{equation}

a=a1q^1+a2q^2+a3q^3\mathbf{a}=a_{1}\hat{\mathbf{q}}_{1}+a_{2}\hat{\mathbf{q}}_{2}+a_{3}\hat{\mathbf{q}}_{3}とすると、(1)(1)は下の通りだ。

∂f∂q1dq1+∂f∂q2dq2+∂f∂q3dq3=a1h1dq1+a2h2dq2+a3h3dq3 \frac{ \partial f}{ \partial q_{1} }dq_{1}+\frac{ \partial f}{ \partial q_{2}}dq_{2}+\frac{ \partial f}{ \partial q_{3}}dq_{3} = a_{1}h_{1}dq_{1}+a_{2}h_{2}dq_{2}+a_{3}h_{3}dq_{3}

ゆえにai=1hi∂f∂qia_{i}=\dfrac{1 }{h_{i}}\dfrac{ \partial f}{ \partial q_{i} }であり、次が成立する。

a=1h1∂f∂q1q^1+1h2∂f∂q2q^2+1h3∂f∂q3q^3 \quad \mathbf{a}=\frac{1 }{h_{1}}\frac{ \partial f}{ \partial q_{1} }\hat{\mathbf{q}}_{1}+\frac{1 }{h_{2}}\frac{ \partial f}{ \partial q_{2} }\hat{\mathbf{q}}_{2}+\frac{1 }{h_{3}}\frac{ \partial f}{ \partial q_{3} }\hat{\mathbf{q}}_{3}

これで上記のベクトルa\mathbf{a}をffの勾配と定義し、∇f\nabla fと表記する。

■

関連項目