%表題   流体力学: 境界条件
%
%履歴   1989-04-21 竹広真一       流体力学の基礎 / 流体力学の基礎方程式
%       1990-04-23 保坂征宏
%       1996-04-23 林 祥介        「基礎法則」へ
%       2000-05-24 竹広真一       「境界条件」を分離
%                                 

%\documentstyle[a4j,12pt,ascmac,twoside,dennou,Depspic]{jarticle}
\documentclass[a4j,12pt]{jarticle}
\usepackage{Dennou6}

\Dtitle{流体力学: 境界条件}
\Dauthor{林 祥介, 竹広 真一}
\Ddate{2000 年 05 月 24 日}
\Dpath{/riron/renzoku/housoku/src/}
\Dnoparindent
\Dparskip

\begin{document}

%\pagenumbering{roman}
\maketitle
\tableofcontents
%\clearpage
%\pagenumbering{arabic}

\begin{abstract}

ここでは, 流体力学において通常よく用いられる境界条件の例を示す. 

流体の運動を完全に記述をするためには, 
質量保存則・運動量保存則・エネルギー保存則に加えて, 
境界条件が必要である. 

境界で与えるべき情報は, 
質量フラックス, 運動量フラックス, エネルギーフラックスである. 
しかし普通, 境界条件としては, 
境界におけるフラックスを直接与えるのではなく, 
物理的意味付けのはっきりした変数, 
例えば, 圧力, 速度の境界における値を与えることが多い. 
その場合, 
与えられた条件が境界における各フラックスの値を決めているかどうかは, 
問題に応じて考察, 確認しなければならない\footnotemark

\footnotetext{
実戦的には, 1次の差分スキームをつくって検討するとよい.}

\end{abstract}

%--------------------------------------------------------------------
\Dparskip
%--------------------------------------------------------------------
\newpage

\section{固体境界}

静止した固体表面において, 流体が満たすべき条件を挙げる. 

質量保存則より, 固体表面を通しての質量流束は $ 0 $ であるから, 

\begin{equation}
   \Dvect{v} \cdot \Dvect{n} = 0 \ \ \ \ \ \ \mbox{at surface}.
\end{equation}

$ \Dvect{n} $ は固体表面の法線ベクトルである. 
\vspace{5mm}

粘性流体の場合は, 
これだけでは条件が不足して解くことができないので, 
境界条件を付け加える必要がある. 
固体表面では流体は密着していると考えて（粘着条件）, 


\begin{equation}
    \Dvect{v} = 0 \ \ \ \ \ \ \mbox{at surface}.
\end{equation}

一般に, 動いている物体の表面での条件を, 
物体が静止している場合をもとに, 考える. 

境界を通しての質量フラックスが 0 であるから, 
境界面の法線方向の速度成分が一致する. すなわち, 

\begin{equation}
   \Dvect{v} \cdot \Dvect{n} = \Dvect{v}^{(s)} \cdot \Dvect{n}
   \ \ \ \ \ \ \mbox{at surface}.
\end{equation}

ただし $\Dvect{v}^{(s)}$ は物体の表面速度である. 

粘性流体の場合の対応する粘着条件は, 
流体の速度と物体表面の速度が一致することである. 

\begin{equation}
  \Dvect{v} = \Dvect{v}^{(s)} \ \ \ \ \ \ \mbox{at surface}.
\end{equation}

さらにこの他に, 熱に関する境界条件が必要である. 
これは, 問題設定により異なるので, ここでは述べない\footnotemark[2]. 

\footnotetext[2]{
  例えばよく用いられるのは, 
  断熱条件 $ \  q_k =0. \ \ \  \mbox{at surface}$ である. 
  あるいは, Benard 問題のように, $T$ を fix にすることである. }
 
%----------------------------------------------------------------------
\section{2種類の流体の境界面}

2種類の流体の境界面を考える
\footnote{この当たりの詳細な議論は「多成分多相系の流体力学」を参照すべし}. 
境界面が静止している系で考えたとき
境界面を通しての質量流束が同じでなければならない. 
\begin{equation}
   \rho^{(1)}\Dvect{v^{(1)}} \cdot \Dvect{n} 
 = \rho^{(2)}\Dvect{v^{(2)}} \cdot \Dvect{n} \ \ \ \ \ \ \mbox{at surface}.
\end{equation}
ただし, 
添え字 $(1), (2)$ は, 
それぞれ, 流体 $(1), (2)$ の物理量であることを示す. 

境界面が動いている系からだと
\begin{equation}
 = \rho^{(1)}(\Dvect{v^{(1)}}-\Dvect{v^{(s)}}) \cdot \Dvect{n} 
 = \rho^{(2)}(\Dvect{v^{(2)}}-\Dvect{v^{(s)}}) \cdot \Dvect{n} 
   \ \ \ \ \ \ \mbox{at surface}.
\end{equation}

となる. 

特に混じりあわない流体の場合には, 面を通しての質量流束が
 $ 0 $ である. 
したがって境界面の速度を $\Dvect{v^{(s)}}$ とすると

\begin{equation}
   \Dvect{v^{(s)}} \cdot \Dvect{n} 
 = \Dvect{v^{(1)}} \cdot \Dvect{n} 
 = \Dvect{v^{(2)}} \cdot \Dvect{n} \ \ \ \ \ \ \mbox{at surface}.
\end{equation}

すなわち, 境界面の法線方向の速度が等しい. 

粘性流体の場合には粘着条件がよく用いられる. 

\begin{equation}
   \Dvect{v^{(1)}} = \Dvect{v^{(2)}} \ \ \ \ \ \ \ \mbox{at surface}.
\end{equation}

流体同士の境界面では, 
境界面に垂直な方向の運動量流束が連続でなければならない. 
境界面 $df$ を通しての運動量流束は 
$ \Pi_{ik} n_k df =(\rho v_i v_k - \sigma_{ik} ) n_k df$ 
である. 境界面が静止している系で考えると
\begin{eqnarray*}
\Pi_{ik}^{(1)} n_k df & = & \Pi_{ik}^{(2)} n_k df ,  \\
(\rho^{(1)} v_i^{(1)} v_k^{(1)} - \sigma_{ik}^{(1)} ) n_k df
 & = & (\rho^{(2)} v_i^{(2)} v_k^{(2)} - \sigma_{ik}^{(2)} ) n_k df , 
\end{eqnarray*}
である\footnotemark[1]. 
\footnotetext[1]{
これが互いに相手の流体に働く力になる. 
この境界条件は, 
相手の流体に働く力同士は, 大きさが等しく, 
向きが反対になること---作用反作用の法則---と同値である. }

まず, 理想流体の場合には, 
\begin{eqnarray*}
(\rho v_i^{(1)} v_k^{(1)} + p \delta_{ik}^{(1)} ) n_k df
 & = & (\rho v_i^{(2)} v_k^{(2)} + p \delta_{ik}^{(2)} ) n_k df \\ 
(\rho^{(1)}  v_i^{(1)}  v_k^{(1)}  )n_k + p^{(1)} n_i df
 & = & (\rho^{(2)} v_i^{(2)} v_k^{(2)} ) n_k + p^{(2)} n_i df .  
\end{eqnarray*}

混じりあわない場合は (30) より,  
$v_k^{(1)} n_k^{(1)} =v_k^{(2)} n_k^{(2)} = 0$ だから, 
境界条件は, 
\begin{equation}
   p^{(1)} = p^{(2)} \ \ \ \ \ \ \mbox{at surface}.
\end{equation}
すなわち, 圧力が連続であることが境界条件となる. 

粘性流体の場合には, 
粘着条件より $ \Dvect{v}^{(1)} = \Dvect{v}^{(2)} = 0$ だから, 
\begin{equation}
    \sigma_{ik}^{(1)} n_{k} = \sigma_{ik}^{(2)} n_{k}
     \ \ \ \ \ \ \mbox{at surface}. 
\end{equation}
が境界条件である. 

%======================================================================
\newpage
\section{参考文献}

\begin{description}
  \item Batchelor,G.K., 橋本英典 他 訳 : 入門流体力学, 
        東京電機大学出版局, 614pp.

  \item Landau,L.D., Lifshitz,E.M., 竹内 均 訳, 1970 : 流体力学1, 
        東京図書, 280pp.

  \item 今井　功, 1973 : 流体力学(前編), 裳華房, 428pp.

  \item Glansdorff,P.,Prigogine,I., 松本 元 ,竹山 脇三 訳, 1977 : 
        構造・安定性・ゆらぎ. 
        みすず書房, 297pp.
\end{description}
\vspace{2em}

%======================================================================
\newpage
\section{謝辞}

本稿は 1989 年から 1993 年に東京大学地球惑星物理学科で行われていた, 
流体理論セミナーでのセミナーノートがもとになっている. 
原作版は竹広真一による「流体力学の基礎方程式」 (1989-04-21) であり, 
保坂征宏による改定 (1990-04-23) を経て, 
林祥介/竹広真一によって「連続体力学: 基礎法則」として書き直された 
(1996-04-23). 
さらに竹広真一によって「流体力学: 境界条件」に分割・編集された
(2000-05-24). 
構成とデバッグに協力してくれたセミナー参加者のすべてに
感謝するものである. 

本ドキュメントは
\begin{quote}
  http://www.gfd-dennou.org/library/rironn/renzoku/housoku/pub/
\end{quote}
において, 無保証無責任を原則として公開している. 
原著作者ならびにその他の資源提供者(図等の版元等を含む)
の諸権利に抵触しない(不利益を与えない)限り, 
資源は自由に利用していただいて構わない. 
\copyright 林祥介・竹広真一 (Y.-Y. Hayashi and S. Takehiro) 1989. 



\end{document}

% Local Variables:
% TeX-command-default: "pLaTeX"
% End: