%表題 磁気流体力学における磁場に関する保存則
%
%履歴 1990/07/03 吉田茂生      セミナーのレジュメ
%     2002/03/10 佐々木洋平    TeX 化開始
%     2002/07/05 佐々木洋平    残りはヘリシティの保存
%     2002/07/12 佐々木洋平    暫定版完成. あとは絵を追加すること. 
%

\documentclass[a4j,12pt,notitlepage]{jarticle}
\usepackage{times} % pdf ファイルを作成する場合にはコメントを外す.
\usepackage{amssymb}
\usepackage{Dennou6}

\Dtitle[磁場に関する保存則]{MHD における磁場に関する保存則}
\Dauthor[吉田・竹広・佐々木・林]
        {吉田茂生・竹広真一・佐々木洋平・林祥介}
\Ddate{2002 年 7 月 12 日}
\Dpath{/riron/mhd/jisoku/src/}
\Dnoparindent

\begin{document}

\maketitle
\tableofcontents
 
\newpage
\markright{\abstractname}
\begin{abstract}
 \noindent
 流体が完全導体とみなせる場合,
 磁場の時間変化は完全導体の誘導方程式
 \begin{equation}
  \DP{\Dvect{B}}{t} = \Drot(\Dvect{u}\times\Dvect{B}) \nonumber
 \end{equation}
 によって記述される. 
 これは完全流体の渦度方程式
 \begin{equation}
  \DP{\Dvect{\omega}}{t} = \Drot(\Dvect{u}\times\Dvect{\omega}) \nonumber
 \end{equation}
 と同じ形をしているので,
 完全流体の渦度に関する保存則は,
 完全導体の磁場についても成立する.

\vspace{1em}
\noindent 
 ここでは
 完全導体の磁場に関する保存則である,
 \begin{itemize}
  \item Alf\'ven の定理(磁束の保存, 磁力線凍結)
  \item  磁気ヘリシティー,クロスヘリシティーの保存則 
 \end{itemize} 
 を導出する.
 
\end{abstract}
%----------------------------------------------------------------------
\Dparskip

%======================================================================
\newpage

 \section{完全導体の誘導方程式}

  磁気流体力学(以下, MHD) における磁場の時間変化は誘導方程式
  \begin{equation}
   \Deqlab{誘導方程式}
    \DP{\Dvect{B}}{t} = \Drot(\Dvect{v}\times\Dvect{B}) + \eta\Dlapla\Dvect{B}
  \end{equation}
  及び磁場のソレノイダル条件
  \begin{equation}
   \Deqlab{磁場のソレノイダル条件}
    \Ddiv\Dvect{B} = 0
  \end{equation}
  によって記述される.
  ここで
  $\Dvect{B}$ は磁束密度,    
  $\Dvect{v}$ は速度ベクトルである.
  $\eta$      は磁気拡散率であり,
  透磁率 $\mu$, 電気伝導度 $\sigma$ を用いて
  \begin{equation}
   \eta \equiv \frac{1}{\mu \sigma}
  \end{equation}
  と定義される
  \footnote{
  磁気流体力学の定式化については
  シリーズ「磁気流体力学の定式化」
  を参照されたい.
  }.
  ここでは \Deqref{誘導方程式} の右辺第 2 項が無視できる条件を考える.
  
  \Deqref{誘導方程式} に対して
  時間のスケールを $\tau$, 
  長さのスケールを $l$,
  速度のスケールを $v\sim l/\tau$
  としてスケーリングすると
  \begin{equation}
   \DP{\Dvect{B}}{t_{*}}
    =   \Drot(\Dvect{v}_{*}\times\Dvect{B}) 
      + \frac{1}{R_{m}}\Dlapla\Dvect{B}
  \end{equation}
  となる.
  ここで $t_{*}$は無次元化された時間を, 
  $\Dvect{v}_{*}$ は無次元化された速度を表すとする. 
  $R_{m}$ は磁気レイノルズ数であり, 次式で定義される:
  \begin{equation}
   R_{m} \equiv \frac{l v}{\eta} = \frac{\tau_{M}}{\tau}.
  \end{equation}
  $\tau_{M}\equiv l^{2}/\eta$ は磁気拡散時間と呼ばれる.
  すなわち磁気レイノルズ数は
  磁気拡散時間と系の時間スケールとの比を表している.
  これが十分大きい場合には,
  誘導方程式 \Deqref{誘導方程式} の
  右辺第 2 項を無視することができる.

  磁気レイノルズ数が大きい, すなわち
  磁気拡散時間が十分長い
  \footnote{
  通常, 
  磁気レイノルズ数が十分大きい現象は
  磁気拡散時間が十分長い現象として解釈される.  
  なぜならば,
  MHD で記述される現象の時間スケールは
  「光速が伝播する時間スケールよりも十分大きくなければならない」
  という条件が, 時間スケールに対する下限として存在する.
  そのため磁気レイノルズ数が大きい現象を
  時間スケールが十分小さい現象とは考えない.
  }
  ための条件は
  \begin{itemize}
   \vspace{-1em}
   \item 磁気拡散率が十分小さい, 
	 すなわち電気伝導度 $\sigma$ が十分大きい.
   \item 系の空間スケール $l$ が十分大きい.
  \end{itemize}
  のいずれか, もしくは両方である.
  この時, 磁場の時間変化は
  \begin{equation}
   \Deqlab{完全誘導}
   \DP{\Dvect{B}}{t} = \Drot(\Dvect{v}\times\Dvect{B})
  \end{equation}
  によって記述される.
  特に
  電気伝導度が無限大 ($\sigma \to \infty$) の導体は完全導体と
  呼ばれるため, \Deqref{完全誘導} は「完全導体の誘導方程式」と呼ばれる.

  \Deqref{完全誘導} は完全流体
  \footnote{
  本文書では粘性が存在しない流体を完全流体と呼ぶ.
  しかしながらプラズマ物理学の本の幾つかでは,
  電気伝導度が無限大の場合の磁気流体を完全流体と呼び,
  粘性が存在しない流体はそのまま非粘性流体と呼ぶこともあるので
  注意されたい.
  }
  に関する渦度方程式
  \begin{equation}
   \Deqlab{渦度}
    \DP{\Dvect{\omega}}{t} = \Drot(\Dvect{v}\times\Dvect{\omega})
  \end{equation}
  と同じ形をしている 
  (ここで $\Dvect{\omega}\equiv\Drot\Dvect{v}$ は流体の渦度である).
  よって完全流体の渦度について成立する法則は
  \Deqref{完全誘導}で記述される磁場についてもそのまま成立する事になる.
  以下では,
  完全導体の誘導方程式 \Deqref{完全誘導}が成立する場合の
  磁場に関する保存則を導出する.
  
%======================================================================
\newpage
\section{種々の量の定義}

 前節で述べたように, 
 完全流体の渦度方程式について成立する事は
 完全導体の誘導方程式についても成立する.
 よってその類似性を比較しやすいように,
 流体における渦線, 渦管, 循環にあたる量として,
 磁力線, 磁力管, 磁束を定義する.

%--------------------------------------------------------------------- 
 \subsection{磁力線}
 
 磁力線は次式で定義される:
 \begin{equation}
  \Deqlab{磁力線の定義}
  \frac{\Dd x}{B_{x}} = \frac{\Dd y}{B_{y}} = \frac{\Dd z}{B_{z}}
 \end{equation}
 ここで $\Dd x, \Dd y, \Dd z$ は各々 $x,y,z$ 方向の線素,
 $B_{x},B_{y},B_{z}$ は各々磁場の $x,y,z$ 成分である.
 すなわち, 磁力線は
 曲線上の各点において磁場が曲線の接線になるように定義される.

%--------------------------------------------------------------------- 
 \subsection{磁力管}
 磁力管は任意の閉曲線 $C$ 上の各点を通る磁力線によって形成される曲面
 からなる. 
 磁力管の形成する面上の任意の点 $\Dvect{x}$ における法線ベクトルを 
 $\Dvect{n}(\Dvect{x},t)$ とするとき, \Deqref{磁力線の定義}より
 \begin{equation}
  \Deqlab{磁力管と直行}
  \Dvect{n}(\Dvect{x},t)\perp\Dvect{B}(\Dvect{x},t)
 \end{equation}
 が成立する.

%--------------------------------------------------------------------- 
 \subsection{磁束}
 
 磁束 $\Phi$ は任意の閉曲線 $C$ によって囲まれる面積を $S$ とするとき,
 \begin{equation}
  \Phi 
   \equiv \int_{S}\Dvect{B}\cdot\Dd \Dvect{S}
     =    \oint_{C} \Dvect{A}\cdot\Dd \Dvect{l}
 \end{equation}
 で定義される. 
 ここで $\Dvect{A}$ は磁気ポテンシャルであり,
 \begin{equation}
  \Dvect{B} = \Drot\Dvect{A}
 \end{equation}
 をみたす.

%--------------------------------------------------------------------- 
 \subsection{磁力管の強さ}
 
 流体における渦管の強さに対応する量として磁力管の強さを定義する.
 磁力管を輪切りにする断面 $S_{i}$ で磁束 $\Phi_{i}$ 作る:
 \begin{equation}
  \Phi_{i} 
   \equiv \int_{S_{i}} \Dvect{B}\cdot\Dd \Dvect{S}
   \quad (i=1,2,3,\cdots)
 \end{equation}
 ここで磁力管中を任意の二つの断面 $S_{i},S_{j}~(i\neq j)$
 によって切り取り, できた閉曲面 $S$ で磁場を面積分すると,
 \begin{eqnarray}
  \int_{S} \Dvect{B} \Dd \Dvect{S} 
   &=&   \int_{S_{i}}\Dvect{B}\cdot \Dd \Dvect{S}_{i}
       + \int_{S_{側面}}\Dvect{B}\cdot \Dd \Dvect{S}_{側面}
       - \int_{S_{j}}\Dvect{B}\cdot \Dd \Dvect{S}_{j} \nonumber \\
   &=& \int_{V} \Ddiv\Dvect{B} \Dd V \nonumber \\
   &=& 0
 \end{eqnarray}
 となる(途中の式変形ではガウスの定理を用いた).
 ここで \Deqref{磁力管と直行}より, 
 磁力管の側面において 
 \begin{equation}
  \int_{S_{側面}}\Dvect{B}\cdot \Dd \Dvect{S}_{側面} = 0
 \end{equation}
 となる. よって
 \begin{equation}
  \Phi_{i} = \int_{S_{i}}\Dvect{B}\cdot \Dd \Dvect{S}_{i}
   = \int_{S_{j}}\Dvect{B}\cdot \Dd \Dvect{S}_{j} = \Phi_{j} 
 \end{equation}
 となる. 
 すなわち磁束 $\Phi_{i}$ は断面 $S_{i}$ のとり方によらない
 磁力管固有の量である. これを磁力管の強さという.

%======================================================================
\newpage
\section[Alf\'enの定理]{Alf\'ven の定理}

 Alf\'en の定理とは,
 \begin{itemize}
  \vspace{-1em}
  \item 磁束 $\Phi$ がラグランジュ保存量であること.
  \item 磁力線が流体に凍結されていること.
  \vspace{-1em}
 \end{itemize}
 の二つである. 
 これらは保存量を
 磁力管に注目して記述するか,
 磁力線に注目して記述するかの違いであり,
 両者は同じ事を言い替えただけである.

 \subsection{磁束の保存}

 流体とともに動く系からみると磁束は不変である:
 \begin{equation}
  \Deqlab{磁束の保存}
  \frac{D}{Dt}\Phi = 0.
 \end{equation}
 ここで $D/Dt$ はラグランジュ微分
 \begin{equation}
  \frac{D}{Dt} = \DP{}{t} + \Dvect{v}\cdot\Dgrad
 \end{equation}
 である.
 \Deqref{磁束の保存} は
 磁束はラグランジュ保存量である, とも言う.
 以下では\Deqref{磁束の保存}を示す.

 流体に固定された面 $S = S(t)$ における磁束 $\Phi$ の
 時間変化を考えると,
 \begin{eqnarray}
  \frac{D}{Dt}\Phi  
  &=&\frac{D}{Dt}\int_{S(t)}\Dvect{B} \cdot \Dd \Dvect{S}(t) \nonumber \\
  &=& \lim_{\delta t \to 0} \frac{1}{\delta t}
      \biggl[
        \int_{S(t+\delta t)}\Dvect{B}(t+\delta t) \cdot \Dd \Dvect{S}
      - \int_{S(t)}\Dvect{B}(t) \cdot \Dd \Dvect{S}
      \biggr] \nonumber \\
  &=& \lim_{\delta t \to 0}\frac{1}{\delta t}
      \biggl[
        \int_{S(t+\delta t)}\Dvect{B}(t+\delta t) \cdot \Dd \Dvect{S}
      - \int_{S(t+\delta t)}\Dvect{B}(t) \cdot \Dd \Dvect{S}
      \biggr] \nonumber \\
  && {}+\lim_{\delta t \to 0}\frac{1}{\delta t}
      \biggl[
      \int_{S(t+\delta t)}\Dvect{B}(t) \cdot \Dd \Dvect{S}
      - \int_{S(t)}\Dvect{B}(t) \cdot \Dd \Dvect{S}
      \biggr] \nonumber \\
  &=& \int_{S(t)}\DP{\Dvect{B}}{t} \cdot \Dd \Dvect{S}
    +\lim_{\delta t \to 0}\frac{1}{\delta t}
      \biggl[
      \int_{S(t+\delta t)}\Dvect{B}(t) \cdot \Dd \Dvect{S}
      - \int_{S(t)}\Dvect{B}(t) \cdot \Dd \Dvect{S}
      \biggr] 
 \end{eqnarray}
 となる. 
 ここで 
 $S(t +\delta t)$ と $S(t)$ とで囲まれる体積を $V$, その側面を $\tilde{S}$
 とすると,
 \begin{eqnarray}
  && 
     \int_{S(t+\delta t)}\Dvect{B}(t) \cdot \Dd \Dvect{S}
   - \int_{S(t)}\Dvect{B}(t) \cdot \Dd \Dvect{S} \nonumber \\
  && =  
   \biggl[
     \int_{S(t+\delta t)}\Dvect{B}(t) \cdot \Dd \Dvect{S}
   - \int_{S(t)}\Dvect{B}(t) \cdot \Dd \Dvect{S} 
   + \int_{\tilde{S}}\Dvect{B}(t) \cdot \Dd \Dvect{S}
   \biggr]
   - \int_{\tilde{S}}\Dvect{B}(t) \cdot \Dd \Dvect{S} \nonumber \\
  && =
     \int_{V}\Ddiv\Dvect{B}\Dd V
   - \oint_{C}\Dvect{B}(t) \cdot (\Dd \Dvect{l}\times\Dvect{v}\delta t)
   \quad 
   (\because \Dd \tilde{\Dvect{S}} = \Dd \Dvect{l} \times \Dvect{v}\delta t )
   \nonumber \\
  && =
   - \oint_{C}(\Dvect{v}\delta t\times \Dvect{B}(t))\cdot \Dd \Dvect{l}
   \nonumber \\
  && =
   - \int_{S}\Drot(\Dvect{v}\delta t\times \Dvect{B}(t))\cdot \Dd \Dvect{S}
 \end{eqnarray}
 となる. 
 ここで 磁場のソレノイダル条件 \Deqref{磁場のソレノイダル条件} を用いた.
 このことから
 \begin{equation}
  \frac{D}{Dt}\Phi 
   =\int_{S(t)}\biggl[
   \DP{\Dvect{B}}{t} 
   - \Drot(\Dvect{v}\times \Dvect{B}(t))
   \biggr]
      \cdot \Dd \Dvect{S}
 \end{equation}
 となる. 
 すなわち完全導体の誘導方程式 \Deqref{完全誘導} 
 が成立するなら磁束の保存 \Deqref{磁束の保存}が成立する.
 
 以上の導出過程は面に固定された座標 $(a,b,c)$ を用いる事で
 幾分簡単に行える. ある時刻に $c=const$ である面を選ぶ. 
 このとき面素片ベクトルは
 \begin{equation}
  \Dd S_{i} = \varepsilon_{ijk}\DP{x_{j}}{a}\DP{x_{k}}{b}\Dd a\Dd b
 \end{equation}
 となり, その時間発展は
 \begin{eqnarray*}
  \frac{D}{Dt}\Dd S_{i} 
   &=& 
     \varepsilon_{ijk}\DP{v_{j}}{a}\DP{x_{k}}{b}\Dd a\Dd b
   + \varepsilon_{ijk}\DP{x_{j}}{a}\DP{v_{k}}{b}\Dd a\Dd b \\
   &=& 
     \varepsilon_{ijk}
     \biggl(
         \DP{v_{j}}{x_{l}}\DP{x_{l}}{a}\DP{x_{k}}{b}
       + \DP{x_{j}}{a}\DP{v_{k}}{x_{l}}\DP{x_{l}}{b}
     \biggr)\Dd a\Dd b \\
   &=& 
     \varepsilon_{ijk}\DP{v_{j}}{x_{l}}
     \biggl(
          \DP{x_{l}}{a}\DP{x_{k}}{b}
	- \DP{x_{k}}{a}\DP{x_{l}}{b}
     \biggr)\Dd a\Dd b  \\
   &=& 
     \varepsilon_{ijk}\DP{v_{j}}{x_{l}}
     \biggl(
      \delta_{ln}\delta_{km}- \delta_{lm}\delta_{kn} 
     \biggr) 
     \DP{x_{n}}{a}\DP{x_{m}}{b}\Dd a\Dd b \\
   &=& 
     \varepsilon_{ijk}\varepsilon_{nmr}\varepsilon_{lkr}
     \DP{v_{j}}{x_{l}}\DP{x_{n}}{a}\DP{x_{m}}{b}\Dd a\Dd b \\
   &=& 
     \DP{v_{j}}{x_{j}}\Dd S_{i} - \DP{v_{j}}{x_{i}}\Dd S_{j}
 \end{eqnarray*}
 であるから,
 \begin{eqnarray}
  \Deqlab{磁束変化2}
  \frac{D}{Dt} \int_{S}\Dvect{B}\cdot\Dd \Dvect{S}
   &=&
   \int_{S}
   \biggl[
     \frac{D\Dvect{B}}{Dt}\cdot\Dd \Dvect{S} 
   + \Dvect{B}\cdot\biggl(\frac{D}{Dt}\Dd \Dvect{S}\biggr)
   \biggr] \nonumber \\
 &=&
   \int_{S}
   \biggl[
     \biggl\{
        \DP{\Dvect{B}}{t}+(\Dvect{v}\cdot\Dgrad\Dvect{B})
     \biggr\}\cdot\Dd \Dvect{S} 
     + B_{i}\DP{v_{j}}{x_{j}}\Dd S_{i}
     - B_{j}\DP{v_{j}}{x_{i}}\Dd S_{j}
     \biggr] \nonumber \\
 &=&
   \int_{S}
   \biggl[
        \DP{\Dvect{B}}{t}
	+ (\Dvect{v}\cdot\Dgrad)\Dvect{B}
        + (\Ddiv\Dvect{v})\Dvect{B}-(\Dvect{B}\cdot\Dgrad)\Dvect{v}
   \biggr]\cdot\Dd \Dvect{S}  
 \end{eqnarray}
 となる. 
 一方で \Deqref{完全誘導}の右辺を展開すると,
 \begin{eqnarray}
  \DP{\Dvect{B}}{t} 
   &=& \Drot(\Dvect{v}\times\Dvect{B}) \nonumber \\
   &=&   (\Dvect{B}\cdot\Dgrad)\Dvect{v} + (\Ddiv\Dvect{B})\Dvect{v}
       - (\Dvect{v}\cdot\Dgrad)\Dvect{B} - (\Ddiv\Dvect{v})\Dvect{B},
       \nonumber 
 \end{eqnarray} 
 ここで磁場のソレノイダル条件
 \Deqref{磁場のソレノイダル条件} を用いることで
 \begin{equation}
  \Deqlab{完全誘導2}
   \DP{\Dvect{B}}{t} 
      =  (\Dvect{B}\cdot\Dgrad)\Dvect{v} 
       - (\Dvect{v}\cdot\Dgrad)\Dvect{B} - (\Ddiv\Dvect{v})\Dvect{B}.  
 \end{equation}
 よって, \Deqref{磁束変化2} と \Deqref{完全誘導2}を比較する事で
 \Deqref{磁束の保存}が成立する事がわかる.

%----------------------------------------------------------------------
\newpage
 \subsection{磁力線の凍結}
 
 磁力線が流体に凍結されている事を示す.
 
 「磁力線が流体に凍結されている」とは, 次の事を意味している.
 ある瞬間に限りなく接近した $2$ つの流体粒子が
 $1$ つの磁力線の上にあるならば,
 それらの流体粒子はいつでもその磁力線の上に存在している.
 さらにその磁力線の強さ $\Dvect{B}/\rho$ 
 はそれら $2$ つの粒子間の距離に比例して変化する.
 直観的には, 
 磁力線が流体粒子に貼り付いて一緒に移動することを意味している.
 
 完全導体の誘導方程式 \Deqref{完全誘導2} をラグランジュ微分を
 用いて書き直すと
 \begin{equation}
  \Deqlab{完全誘導3}
  \frac{D\Dvect{B}}{Dt} 
   =  (\Dvect{B}\cdot\Dgrad)\Dvect{v} + (\Ddiv\Dvect{v})\Dvect{B}  
 \end{equation}
 となる.
 ここで流体の連続の式
 \begin{equation}
  \DP{\rho}{t} + \Ddiv(\rho\Dvect{v}) =0
 \end{equation}
 より
 \begin{eqnarray*}
  &&\DP{\rho}{t} + \Ddiv(\rho\Dvect{v}) =0  \\
  &&{} \to  
    {}\DP{\rho}{t} +(\Dvect{v}\cdot\Dgrad)\rho+ \rho\Ddiv\Dvect{v}=0 \\
  &&{} \to  
    {}\frac{D \rho}{Dt} + \rho\Ddiv\Dvect{v}=0
 \end{eqnarray*}
 となるから, 
 この式と \Deqref{完全誘導3}より $\Ddiv\Dvect{v}$ を消去する事で,
 \begin{equation}
  \Deqlab{完全誘導4}
  \frac{D}{Dt}\biggl(\frac{\Dvect{B}}{\rho}\biggr)
   = \biggl(\frac{\Dvect{B}}{\rho}\cdot\Dgrad\biggr)\Dvect{v}
 \end{equation}
 が得られる. 

 この式は流体中の微小線素 $\delta \Dvect{r}$ の時間発展 
 \begin{equation}
  \Deqlab{線素発展}
  \frac{D}{Dt}(\delta \Dvect{r}) = (\delta \Dvect{r}\cdot \Dgrad )\Dvect{v}
 \end{equation}
 と同じ形をしている.
 従って, 
 ある時刻 $t=t_{0}$ において, 
 $\Dvect{B}/\rho$ と $\delta \Dvect{r}$ 
 という二つのベクトルの向きが一致しているならば,
 その後も常に一致しており, 
 その長さは互いに比例して変化する事になる.
 
 これを数学的に表すならば,
 ある時刻 $t=t_{0}$ に 位置 $\Dvect{a} = (a,b,c)$ において 
 \begin{equation}
  \frac{\Dvect{B}}{\rho}(\Dvect{a},t_{0}) 
   = \alpha \delta \Dvect{r}(\Dvect{a},t_{0}), \quad \alpha=const
 \end{equation}
 ならば, 任意の時刻 $t$ と 位置 $\Dvect{x}$ において
 \begin{equation}
  \frac{\Dvect{B}}{\rho}(\Dvect{x},t) 
   = \alpha \delta \Dvect{r}(\Dvect{x},t)
 \end{equation}
 が成立することになる.

 よって磁場の時間発展の形式的な解は
 \begin{eqnarray}
  \frac{\displaystyle{\frac{B_{i}(\Dvect{x},t)}{\rho(\Dvect{x},t)}}}
   {\displaystyle{\frac{B_{j}(\Dvect{a},t_{0})}{\rho(\Dvect{a},t_{0})}}}
   = \frac{\delta r_{i}(\Dvect{x},t)}{\delta r_{j}(\Dvect{a},0)}  
 \end{eqnarray}
 もしくは
 \begin{equation}
  \frac{B_{i}(\Dvect{x},t)}{\rho(\Dvect{x},t)}
   =\frac{B_{i}(\Dvect{a},t_{0})}{\rho(\Dvect{a},t_{0})} \DP{x_{i}}{a_{j}},
   \quad
   \delta r_{i}(\Dvect{x},t) = \delta r_{j}(\Dvect{a},t_{0})\DP{x_{i}}{a_{j}}
 \end{equation}
 となる.
 
 次に磁力線をなしている線素 $\delta \Dvect{r}$ について
 $\displaystyle{ \frac{1}{\rho} \frac{B_{i}}{\delta r_{i}} }$
 の時間発展を考えると:
 \begin{eqnarray}
  \frac{D}{Dt}\biggl(\frac{1}{\rho}\frac{B_{i}}{\delta r_{i}}\biggr)
   &=&   \frac{1}{\delta r_{i}}\frac{D}{Dt}\biggl(\frac{B_{i}}{\rho}\biggr)
       - \frac{B_{i}}{\rho}\frac{1}{\delta r_{i}^{2}}\frac{D}{Dt}(\delta r_{i}).
 \end{eqnarray}
 ここで \Deqref{完全誘導4}, \Deqref{線素発展}より,
 \begin{eqnarray}
  \frac{1}{\delta r_{i}}\frac{D}{Dt}\biggl(\frac{B_{i}}{\rho}\biggr)
    - \frac{B_{i}}{\rho}\frac{1}{\delta r_{i}^{2}}\frac{D}{Dt}(\delta r_{i}) 
  &=& 
  \frac{1}{\delta r_{i}}\frac{B_{j}}{\rho}\DP{v_{i}}{x_{j}}
    - \frac{B_{i}}{\rho}\frac{1}{\delta r_{i}^{2}}\delta r_{j}\DP{v_{i}}{x_{j}}
    \nonumber \\ 
  &=& 
  \frac{1}{\rho \delta r_{i}}\DP{v_{i}}{x_{j}}
    \biggl(B_{j} - B_{i}\frac{\delta r_{j}}{\delta r_{i}}\biggr) 
    \nonumber \\ 
  &=& 0 \quad 
   \biggl(\because \frac{\delta r_{j}}{\delta r_{i}} = \delta_{ij}\biggr)
 \end{eqnarray}
 となる. よって
 \begin{equation}
  \Deqlab{磁力線伸縮}
  \frac{D}{Dt}\biggl(\frac{1}{\rho}\frac{B_{i}}{\delta r_{i}}\biggr) = 0.
 \end{equation}
 が成立する.
 
 \Deqref{磁力線伸縮} の直観的な意味は, 
 非圧縮性流体の場合を考えると理解しやすい. 
 流体が非圧縮性流体の場合には $\rho = const$ であるから,
 $B_{i}/\delta r_{i}$ が常に一定となる. 
 つまり磁力線が伸びるとそれに比例して磁場が強くなる.
 このことは「磁力線は伸びると強くなる」と言われる.
 
%=============================================================
\newpage
 \section{ヘリシティの保存}

\iffalse
 保存力場中の非粘性順圧流体において,
 境界で流れの法線成分が存在しないならばヘリシティ
 \begin{equation}
  \Deqlab{ヘリシティ}
  \int_{V} \Dvect{u}\cdot\Dvect{\omega}\Dd V
 \end{equation}
 が保存する事が知られている.
 本節では \Deqref{ヘリシティ} に対応する保存量として,  
 \fi
 境界で磁場の法線成分が存在しない場合に,
 磁気ヘリシティ
 \begin{equation}
  \Deqlab{磁気ヘリシティ}
  \int \Dvect{A}\cdot\Dvect{B} \Dd V
 \end{equation}
 及びクロスヘリシティ
 \begin{equation}
  \Deqlab{クロスヘリシティ}
  \int \Dvect{v}\cdot\Dvect{B} \Dd V =\int \Dvect{A}\cdot\Dvect{\omega} \Dd V
 \end{equation}
 がラグランジュ保存量であることを示す. 
 以下では, 流体は完全導体かつ非粘性順圧磁気流体であり,
 保存力場中にあるとする.

 先ず磁気ヘリシティについて考察する.
 完全導体の誘導方程式 \Deqref{完全誘導} を
 磁気ポテンシャル $\Dvect{A}$ と静電ポテンシャル $\phi$ を用いて書き直すと
 \begin{equation}
  \Deqlab{誘導A}
   \frac{D \Dvect{A}}{D t}
   = (\Dvect{v}\cdot\Dgrad)\Dvect{A} 
   + \Dvect{v}\times\Dvect{B}-\Dgrad \phi
 \end{equation}
 となる. 

 \Deqref{完全誘導4} 及び \Deqref{誘導A} より
 \begin{eqnarray}
  \frac{D}{Dt} \biggl(\frac{\Dvect{B}}{\rho} \cdot \Dvect{A} \biggr)
   &=&
 \Dvect{A}\cdot\frac{D}{Dt} \biggl(\frac{\Dvect{B}}{\rho}\biggr)
 +
 \biggl(\frac{\Dvect{B}}{\rho}\biggr)\cdot\frac{D\Dvect{A}}{Dt} 
 \nonumber \\
 &=& 
 \Dvect{A}\cdot\biggl(\frac{\Dvect{B}}{\rho}\cdot\Dgrad\biggr)\Dvect{v}
 +
 \biggl(\frac{\Dvect{B}}{\rho}\biggr)\cdot
 \{(\Dvect{v}\cdot\Dgrad)\Dvect{A} 
   + \Dvect{v}\times(\Dvect{B})+\Dgrad \phi\}
 \nonumber \\
 &=& 
  \biggl(\frac{\Dvect{B}}{\rho}\cdot\Dgrad\biggr)
    (\Dvect{A}\cdot\Dvect{v} - \phi).
 \end{eqnarray}
 これを全体積で積分すると
 \begin{eqnarray}
  \frac{D}{Dt}
   \int_{V} \Dvect{A}\cdot\Dvect{B} \Dd V
   &=&
   \int_{V} \frac{D}{Dt}\biggl(\frac{\Dvect{A}\cdot\Dvect{B}}{\rho}\biggr)\rho\Dd V
   \nonumber \\
   &=&
   \int_{V} (\Dvect{B}\cdot\Dgrad)(\Dvect{A}\cdot\Dvect{v} - \phi)\Dd V
   \nonumber \\
   &=&
   \int_{S} (\Dvect{A}\cdot\Dvect{v} - \phi)\Dvect{B}\cdot \Dd \Dvect{S}
   \nonumber \\
   &=&
   0. \quad (\because \Dvect{B}\cdot \Dd \Dvect{S} = 0~\mbox{on boundary}) 
 \end{eqnarray}
 よって 
 磁気ヘリシティ\Deqref{磁気ヘリシティ}はラグランジュ保存量であることが示された
 \footnote{
 磁気ヘリシティはゲージ変換に対して不変である.
 ゲージ変換によってポテンシャルが各々
 \begin{equation}
  \Dvect{A} \to \Dvect{A}' + \Dgrad \chi, \quad
   \phi \to \phi' + \DP{\chi}{t}
 \end{equation}
 へ変換されたとしても, \Deqref{誘導A} は
 \begin{eqnarray}
  &&\frac{D\Dvect{A}'}{Dt} + \frac{D}{Dt}\Dgrad \chi
   = (\Dvect{v}\cdot\Dgrad)\Dvect{A}'+(\Dvect{v}\cdot\Dgrad)\Dgrad\chi
   + \Dvect{v}\times(\Drot\Dvect{A}')
   +\Dgrad \phi' +\Dgrad\DP{\chi}{t} \nonumber \\
  &&{}\to
  \frac{D\Dvect{A}'}{Dt}
   = (\Dvect{v}\cdot\Dgrad)\Dvect{A}'
   + \Dvect{v}\times(\Drot\Dvect{A}')+\Dgrad \phi' 
 \end{eqnarray}
 となり, その形を変えない. よって
 磁気ヘリシティの保存はゲージの取り方によらず成立する.
 }.

 次にクロスヘリシティについて考察する.
 流体が非粘性順圧磁気流体であるとした場合の運動方程式は
 \begin{equation}
  \Deqlab{運動方程式}
  \frac{D\Dvect{v}}{Dt} = -\Dgrad P + \frac{1}{\rho}\Dvect{J} \times\Dvect{B}
   ,\quad P = \Omega + \int \frac{\Dd p}{\rho} 
 \end{equation}
 となる. ここで $\rho, p$ はそれぞれ流体の密度,圧力,
 $\Omega$ は保存力場のポテンシャルである.
 また, $\Dvect{J}$ は電流密度ベクトルである. 
 
 \Deqref{運動方程式} 及び \Deqref{完全誘導4} より
 \begin{eqnarray}
  \frac{D}{Dt} \biggl(\Dvect{u}\cdot\frac{\Dvect{B}}{\rho}  \biggr)
   =
   \biggl(\frac{\Dvect{B}}{\rho}\cdot\Dgrad\biggr)
    \biggl(-P + \frac{1}{2}\Dvect{v}^{2} \biggr).
 \end{eqnarray}
 これを全体積で積分すると
 \begin{eqnarray}
  \frac{D}{Dt}
   \int_{V} \Dvect{u}\cdot\Dvect{B} \Dd V
   &=&
   \int_{V} \frac{D}{Dt}\biggl(\frac{\Dvect{u}\cdot\Dvect{B}}{\rho}\biggr)\rho\Dd V
   \nonumber \\
   &=&
   \int_{V} 
    (\Dvect{B}\cdot\Dgrad)\biggl(-P + \frac{1}{2}\Dvect{v}^{2} \biggr)   
   \Dd V
   \nonumber \\
   &=&
   \int_{S} \biggl(-P + \frac{1}{2}\Dvect{v}^{2} \biggr)\Dvect{B}\cdot \Dd \Dvect{S}
   \nonumber \\
   &=&
   0. \quad (\because \Dvect{B}\cdot \Dd \Dvect{S} = 0~\mbox{on boundary}) 
 \end{eqnarray} 
 よってクロスヘリシティ
 \Deqref{クロスヘリシティ}はラグランジュ保存量であることが示された.

 クロスヘリシティの物理的な意味は, 
 磁力線 $C_{B}$ と 渦線 $C_{\omega}$ が一回絡んだ状況を考える
 と理解しやすい.
 $C_{B}, C_{\omega}$ 以外 では $\Dvect{B} = 0, \Dvect{\omega}=0$ とする.
 この時, クロスヘリシティは
 \begin{eqnarray}
  \int_{V} \Dvect{u}\cdot\Dvect{B} \Dd V 
   &=&
      \int_{V_{1}} \Dvect{u}\cdot\Dvect{B} \Dd V 
    + \int_{V_{2}} \Dvect{u}\cdot\Dvect{B} \Dd V \nonumber \\
   &=&
      \Phi\oint_{C_{\omega}} \Dvect{u} \cdot\Dd \Dvect{l} 
     + K\oint_{C_{B}} \Dvect{B}\cdot \Dd \Dvect{l} \nonumber \\
   &=& 2 \Phi K
 \end{eqnarray}
 となる. 
 ここで $\Phi$ は磁力線 $C_{B}$ の強さ, 
 $K$ は渦線 $C_{\omega}$の強さ(循環)である.
 よってクロスヘリシティが保存することより,
 ある時刻に一回絡んでいた磁力線と渦線は,
 流れによって変形しても一回絡んだままであることになる.
 すなわちクロスヘリシティの保存は, 
 磁力線と渦線の「絡み」の数が保存する, という幾何学的な意味を持つ.

 より一般に, 流体中に磁力線が $C_{i},i=1,2,...$ と無数に存在している
 場合を考えても,
 \begin{eqnarray}
  \int_{V} \Dvect{u}\cdot\Dvect{B} \Dd V
   = \sum_{i}\Phi_{i}\oint_{C_{i}}\Dvect{u}\cdot\Dd \Dvect{l} 
   = \sum_{i}\Phi_{i}K_{i}.
 \end{eqnarray}
 ここで $K_{i}$ は $C_{i}$ により作られる面を通過する
 渦度のフラックスであり 
 \Deqref{運動方程式}より定数である事がわかる:
 \begin{eqnarray}
  \frac{D}{Dt}\oint_{C_{i}}\Dvect{u}\cdot\Dd \Dvect{l} 
   &=& 
  \oint_{C_{i}}
  \biggl(
     \DP{\Dvect{u}}{t}-\Dvect{u}\times\Dvect{\omega}
  \biggr)\cdot\Dd \Dvect{l}
  \nonumber \\
   &=&
  \oint_{C_{i}}  
  \biggl\{
    \frac{\Dvect{J}\times\Dvect{B}}{\rho} 
     -\Dgrad\biggl(P+\frac{1}{2}\Dvect{u}^{2}\biggr) 
  \biggr\}\cdot\Dd \Dvect{l}
    \nonumber \\
   &=& 0 \quad \mbox(\because \Dvect{B}\parallel \Dd \Dvect{l})
 \end{eqnarray}
 よって一般の場合でも, 磁力線の「絡み」の数が保存することとなる.
 

%======================================================================
\newpage
\addcontentsline{toc}{section}{参考文献}
\section*{参考文献}
\markright{参考文献}
  \begin{description}
    \item Landau, L.D. and Lifshitz, E.M.(著), 
          井上健夫, 安河内昂, 佐々木健(共訳), 1962:
          『電磁気学 -- 連続媒質の電気力学 --』第 1 巻.
          ランダウリフシッツ理論物理学教程, 東京図書.

    \item 吉田茂生, 1990: 電磁流体力学入門. GFDノート『電磁流体力学』第 2 章.

    \item Moffatt,H.K, 1969: The degree of knottedness of tangled vortex lines.
	  {\em J.Fluid Mech.}, {\bf 35}, 117--129
  \end{description}

%======================================================================
\newpage
\addcontentsline{toc}{section}{謝辞}
\section*{謝辞}
\markright{謝辞}

本稿は 1989 年から 1993 年に東京大学地球惑星物理学科で行われていた, 
流体理論セミナーでのセミナーノートがもとになっているものである. 
原作版は吉田茂生による「電磁流体力学入門」 (1990-07-03) であり, 
竹広真一, 佐々木洋平, 吉田茂生, 林祥介によって
「MHD における磁場に関する保存則」として書き直された (2002-07-12). 
セミナー参加者および校正とデバッグに協力してくれたすべての方々に
感謝するものである. 

本資源は, 
地球流体電脳倶楽部の
インターネット上での学術知識の集積と活用の実験の一環として
\begin{center}
 http://www.gfd-dennou.org/arch/riron/mhd/teishiki/pub
\end{center}
において公開されているものである.
\copyright 吉田茂生・竹広真一・佐々木洋平・林祥介
 (S. Yoshida , S. Takehiro, Y. Sasaki and Y.-Y. Hayashi) 2001. 
本資源は, 
著作者の諸権利に抵触しない(迷惑をかけない)限りにおいて
自由に利用していただいて構わない. 
なお, 利用する際には今一度自ら内容を確かめることをお願いする
(無保証無責任原則). 

本資源に含まれる元資源提供者(図等の版元等を含む)からは, 
直接的な形での WEB 上での著作権または使用許諾を得ていない
場合があるが, 
勝手ながら, 「未来の教育」のための実験という
学術目的であることをご理解いただけるものと信じ, 
学術標準の引用手順を守ることで諸手続きを略させていただいている. 
本資源の利用者には, この点を理解の上, 
注意して扱っていただけるようお願いする. 
万一, 不都合のある場合には
\begin{center}
 riron@gfd-dennou.org
\end{center}
まで連絡していただければ幸いである. 

\end{document}

%%% Local Variables: 
%%% mode: yatex
%%% TeX-master: t
%%% End: 
