% % 表題 SPMODEL サンプルプログラムドキュメント % igwave.tex % % 履歴 2004/05/27 小高 正嗣 % 2004/06/08 小高 正嗣 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \documentclass[a4j,12pt]{jarticle} % % プリアンブル % \usepackage{Dennou6} % dennou-sty-6 を使用 \Dtitle[ {\small 2 次元ブシネスクモデル:内部重力波の計算}] % ヘッダ部タイトル, モデル名と同じ { {\large SPMODEL サンプルプログラム} \\ % ここは変更しない 水路領域 \\ 2 次元ブシネスク方程式モデル:\\ % モデル名 内部重力波の計算 \\ {\large igwave1.f90} % プログラム名 } % \Dauthor[小高 正嗣]{小高 正嗣} % 著者 \Ddate{2004 年 6 月 8 日} % 日付 \Dnoparindent % 段落の字下げをしない % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % 本文開始 % % 基本的な章節構成は以下の通り. % % 1. はじめに % 2. 支配方程式系 % 3. 離散化 % 4. 利用モジュールとその他の設定 % 5. 数値実験 % 6. 参考文献 % -. 謝辞 % % 必要に応じて構成を変更すること. % \begin{document} \maketitle % タイトルの作成 \tableofcontents % 目次の作成 \pagebreak %---------------------------------------------------------------------- \Dparskip %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \newpage % 改ページ \section{概要} SPMODEL サンプルプログラム『igwave1.f90』に用いられている基礎方程式と境界 条件, および, このプログラムを用いた数値実験の方法について解説する. 基礎 方程式は 2 次元のブシネスク方程式系である. 計算はスペクトル法を用いて行 い, 展開関数は水平方向にはフーリエ級数, 鉛直方法には境界条件にあわせてフー リエ正弦級数またはフーリエ余弦級数を用いる. 波数切断は三角切断であ る. スペクトル変換と逆変換および微分演算には, SPMODEL ライブラリ(spml) の esc\_module を用いている. 数値実験では水平方向に双極子型に並んだ正負 の渦度強制により励起された内部重力波の計算を行う. {\bf プログラム名} \\ % プログラム名 {\footnotesize igwave1.f90 } {\bf プログラム取得元}\\ % プログラム取得先 {\footnotesize http:\slash \slash www.gfd-dennou.org\slash arch\slash spmodel\slash 2d-channel-esc\slash boussinesq\slash igwave\slash SIGEN.htm } {\bf SPMODEL サンプルプログラム目次}\\ % サンプルプログラム目次(変更しない) {\footnotesize http:\slash \slash www.gfd-dennou.org\slash arch\slash spmodel\slash sample.htm } {\bf SPMODEL の使い方}\\ % SPMODEL の使い方(変更しない) {\footnotesize http:\slash \slash www.gfd-dennou.org\slash arch\slash spmodel\/ } %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \newpage % 改ページ \section{支配方程式系} ここでは支配方程式系と境界条件を記す. \subsection{支配方程式系} 支配方程式系は 2 次元のブシネスク方程式系である. \begin{eqnarray} && \DP{u}{t}+u\DP{u}{x}+v\DP{u}{y} = -\frac{1}{\rho _{0}}\DP{p}{x} + \nu \Dlapla u + F_{u}, \\ && \DP{v}{t}+u\DP{v}{x}+v\DP{v}{y} = -\frac{1}{\rho _{0}}\DP{p}{y} - \frac{\rho}{\rho _{0}}g + \nu \Dlapla v + F_{v}, \\ && \DP{u}{x}+\DP{v}{y} = 0,\\ && \DP{\rho}{t}+u\DP{\rho}{x}+v\DP{(\rho+\rho _{0})}{y} = \kappa \Dlapla \rho. \end{eqnarray} 各記号の定義は\Dtabref{変数, 物理定数の定義}に表す. 以下では基本場の温度と密度の鉛直勾配は一定とする. 渦度 $\zeta$ と流線関数 $\psi$ \begin{eqnarray} &&\zeta =\DP{w}{x}-\DP{u}{y}, \\ &&v=\DP{\psi}{x}, \quad u=-\DP{\psi}{y}, \end{eqnarray} を導入し, 渦度と密度の予報方程式として支配方程式系を表すと以下のようになる. \begin{eqnarray} && \DP{\zeta }{t} + J(\psi, \zeta) + \frac{g}{\rho _{0}}\DP{\rho }{x} = \nu \Dlapla \zeta + F_{\zeta}, \Deqlab{渦度方程式}\\ && \DP{\rho}{t} + J(\psi, \rho) + \DP{\psi}{x}\DP{\rho _{0}}{y} = \kappa \Dlapla \rho. \Deqlab{密度の予報方程式} \end{eqnarray} ここで $J(A, B)$ はヤコビアン \[ J(A,B) = \DP{A}{x}\DP{B}{y} - \DP{A}{y}\DP{B}{x} \] であり, \[ F_{\zeta} = \DP{F_{v}}{x}-\DP{F_{u}}{y} \] である. \begin{table}[h] \begin{center} \begin{tabular}{cll} \hline 記号 &\qquad\qquad& 変数/物理定数 \\ \hline $x$ && 水平座標 \\ $y$ && 鉛直座標 \\ $t$ && 時間 \\ $u$ && $x$ 方向速度 \\ $v$ && $y$ 方向速度 \\ $\rho_{0}$ && 密度(基本場) \\ $\rho $ && 密度偏差 \\ $\psi$ && 流線関数 \\ $\zeta$ && 渦度 \\ $F_{u}, F_{v}$ && 強制項 \\ $g$ && 重力加速度 \\ $\nu$ && 動粘性係数 \\ $\kappa$ && 拡散係数 \\ \hline \end{tabular} \caption{変数, 物理定数の定義} \Dtablab{変数, 物理定数の定義} \end{center} \end{table} \subsection{境界条件} 境界条件は水平に周期境界条件, 鉛直方向には $y=0, y_{m}$ に置いた剛体壁 面で $v=0$, 応力なし, 熱フラックスなしとする. 水平方向の境界条件は計算領域をそれぞれ $x_{m}$ とすると \begin{equation} \zeta(x+x_{m},z) = \zeta(x,z) \end{equation} などと表される. 鉛直方向の境界条件は $y=0, y_{m}$ において, \begin{eqnarray} &&\psi = \mbox{Const.}, \\ &&\zeta = 0, \\ &&\DP{\rho}{y} = 0 \end{eqnarray} である. 境界で与える $\psi $ の値は簡単のため $0$ とした. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \newpage % 改ページ \section{離散化} この節では方程式系の空間離散化および使用した時間積分法について説明し, プロ グラム内で実際に用いられている方程式を記述する. なお, 強制項 $F_{\zeta}=0$ の与え方については「5. 実験設定」において示す. \subsection{空間離散化} 支配方程式\Deqref{渦度方程式}, \Deqref{密度の予報方程式}で, 強制項 $F_{\zeta}=0$ とした場合の離散表現は以下のようになる. \begin{eqnarray} && \DP{\zeta _{i,j}}{t} + \DP{\psi_{i,j}}{x}\DP{\zeta _{i,j}}{y} - \DP{\psi_{i,j}}{y}\DP{\zeta _{i,j}}{x} + \frac{g}{\rho _{0}}\DP{\rho _{i,j}}{x} = \nu \Dlapla \zeta_{i,j} , \Deqlab{離散化された渦度方程式}\\ && \DP{\rho_{i,j}}{t} + \DP{\psi_{i,j}}{x}\DP{\rho _{i,j}}{y} - \DP{\psi_{i,j}}{y}\DP{\rho _{i,j}}{x} + \DP{\psi _{i,j}}{y}\DP{\rho _{0}}{y} = \kappa \Dlapla \rho_{i,j}. \Deqlab{離散化された密度の予報方程式} \end{eqnarray} ここで添字 $i,j$ は格子点 $(x_{i}, y_{j})$ 上の値であることを表す. \subsection{空間方向のスペクトル表現} 空間離散化した支配方程式\Deqref{離散化された渦度方程式}, \Deqref{離散 化された密度の予報方程式}をスペクトル法を用いて表現する. スペクトル 展開は水平方向にフーリエ級数, 鉛直方法には境界条件にあわせてフーリエ正 弦級数またはフーリエ余弦級数を用いて行う. 非線形項を扱う場合は, 先に格 子点上での非線形項の値を計算し, その値のスペクトルを求める方法(変換法) を用いる. 浮力項についても同様に扱う. 以下では $k,l$ をそれぞれ$x,z$ 方向波数, $K,L$ を切断波数, $I,J$ を格子点数とする. $\zeta_{i,j}, \psi _{i,j}, \rho _{i,j}$ はスペクトル逆変換によって以下 のように展開される. \begin{eqnarray} \zeta_{i,j} &=& \sum _{k=-K}^{K}\sum _{l=0}^{L} \exp\left(\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right)\hat \zeta_{k,l}, \\ \psi_{i,j} &=& \sum _{k=-K}^{K}\sum _{l=0}^{L} \exp\left(\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right)\hat \psi_{k,l}, \\ \rho_{i,j} &=& \sum _{k=-K}^{K}\sum _{l=0}^{L} \exp\left(\frac{2\pi i k x_{i}}{x_{m}}\right) \cos\left(\frac{\pi l y_{j}}{y_{m}}\right)\hat \rho_{k,l}. \end{eqnarray} スペクトル係数 $\hat \zeta_{i,j}, \hat \psi _{i,j}, \hat \rho _{i,j}$ は スペクトル変換によって以下のように与えられる. \begin{eqnarray} \hat \zeta_{k,l} &=& \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right)\zeta_{i,j}, \\ \hat \psi_{k,l} &=& \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right)\psi_{i,j}, \\ \hat \rho_{k,l} &=& \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \cos\left(\frac{\pi l y_{j}}{y_{m}}\right)\rho_{i,j}. \end{eqnarray} 以上を用いると支配方程式\Deqref{離散化された渦度方程式}, \Deqref{離散 化された密度の予報方程式}のスペクトル表現は以下のようになる. \begin{eqnarray} \DP{\hat \zeta _{k,l}(t)}{t} &=& - \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right) \left(\DP{\psi}{x}\right)_{i,j} \left(\DP{\zeta}{y}\right)_{i,j} \nonumber \\ && + \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right) \left(\DP{\psi}{y}\right)_{i,j} \left(\DP{\zeta}{x}\right)_{i,j} \nonumber \\ && - \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right) \frac{g}{\rho _{0}}\left(\DP{\rho}{x}\right)_{i,j} \nonumber \\ && - \nu \left[ \left(\frac{2\pi k}{x_{m}}\right)^{2}+ \left(\frac{\pi l}{y_{m}}\right)^{2}\right]\hat \zeta _{k,l}, \Deqlab{渦度方程式のスペクトル表現}\\ \DP{\hat \rho _{k,l}(t)}{t} &=& - \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \cos\left(\frac{\pi l y_{j}}{y_{m}}\right) \left(\DP{\psi}{x}\right)_{i,j} \left(\DP{\rho}{y}\right)_{i,j} \nonumber \\ && + \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \cos\left(\frac{\pi l y_{j}}{y_{m}}\right) \left(\DP{\psi}{y}\right)_{i,j} \left(\DP{\rho}{x}\right)_{i,j} \nonumber \\ && - \frac{\pi l}{y_{m}}\DP{\rho _{0}}{y}\hat \psi _{k,l} - \kappa \left[ \left(\frac{2\pi k}{x_{m}}\right)^{2}+ \left(\frac{\pi l}{y_{m}}\right)^{2}\right]\hat \rho _{k,l}. \Deqlab{密度の予報方程式のスペクトル表現} \end{eqnarray} ここで \begin{eqnarray} \left(\DP{\zeta}{x}\right)_{i,j} &=& \sum _{k=-K}^{K}\sum _{l=0}^{L} \exp\left(\frac{2\pi i k }{x_{m}}\right) \sin\left(\frac{\pi l }{y_{m}}\right) \frac{2\pi i k}{x_{m}}\hat \zeta _{k,l}, \\ \left(\DP{\zeta}{y}\right)_{i,j} &=& \sum _{k=-K}^{K}\sum _{l=0}^{L} \exp\left(\frac{2\pi i k }{x_{m}}\right) \cos\left(\frac{\pi l }{y_{m}}\right) \frac{\pi l}{y_{m}}\hat \zeta _{k,l}, \\ \left(\DP{\psi}{x}\right)_{i,j} &=& \sum _{k=-K}^{K}\sum _{l=0}^{L} \exp\left(\frac{2\pi i k }{x_{m}}\right) \sin\left(\frac{\pi l }{y_{m}}\right) \frac{2\pi i k}{x_{m}}\hat \psi _{k,l}, \\ \left(\DP{\psi}{y}\right)_{i,j} &=& \sum _{k=-K}^{K}\sum _{l=0}^{L} \exp\left(\frac{2\pi i k }{x_{m}}\right) \cos\left(\frac{\pi l }{y_{m}}\right) \frac{\pi l}{y_{m}}\hat \psi _{k,l}, \\ \left(\DP{\rho}{x}\right)_{i,j} &=& \sum _{k=-K}^{K}\sum _{l=0}^{L} \exp\left(\frac{2\pi i k }{x_{m}}\right) \cos\left(\frac{\pi l }{y_{m}}\right) \frac{2\pi i k}{x_{m}}\hat \rho _{k,l}, \\ \left(\DP{\rho}{y}\right)_{i,j} &=& - \sum _{k=-K}^{K}\sum _{l=0}^{L} \exp\left(\frac{2\pi i k }{x_{m}}\right) \sin\left(\frac{\pi l }{y_{m}}\right) \frac{\pi l}{y_{m}}\hat \rho _{k,l} \end{eqnarray} である. \subsection{時間積分} ここでは時間積分法について記述し, プログラム内で実際に用いられている方 程式を記述する. 以下では $\Delta t$ を時間格子間隔, 時刻 $\tau \Delta t$ における $\hat \zeta_{k,l}$ の値を $\hat \zeta_{k,l}^{\tau}$ 等と表す. 時間方向の離散化は Euler スキームを用いて行う. 時空間方向に離散化され た方程式は以下のように表される. \begin{eqnarray} &&\hat \zeta _{k,l}^{\tau +1} = \hat \zeta _{k,l}^{\tau} + \Delta t \hat F_{k,l}^{\tau}, \\ &&\hat \rho _{k,l}^{\tau +1} = \hat \rho _{k,l}^{\tau} + \Delta t \hat G_{k,l}^{\tau}, \end{eqnarray} \begin{eqnarray} \hat F_{k,l}^{\tau} &=& - \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right) \left(\DP{\psi}{x}\right)_{i,j}^{\tau} \left(\DP{\zeta}{y}\right)_{i,j}^{\tau} \nonumber \\ && + \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right) \left(\DP{\psi}{y}\right)_{i,j}^{\tau} \left(\DP{\zeta}{x}\right)_{i,j}^{\tau} \nonumber \\ && - \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right) \frac{g}{\rho _{0}}\left(\DP{\rho}{x}\right)_{i,j}^{\tau} \nonumber \\ && - \nu \left[ \left(\frac{2\pi k}{x_{m}}\right)^{2}+ \left(\frac{\pi l}{y_{m}}\right)^{2}\right]\hat \zeta _{k,l}^{\tau}, \Deqlab{渦度の増分} \\ \hat G_{i,j}^{\tau} &=& - \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \cos\left(\frac{\pi l y_{j}}{y_{m}}\right) \left(\DP{\psi}{x}\right)_{i,j}^{\tau} \left(\DP{\rho}{y}\right)_{i,j}^{\tau} \nonumber \\ && + \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \cos\left(\frac{\pi l y_{j}}{y_{m}}\right) \left(\DP{\psi}{y}\right)_{i,j}^{\tau} \left(\DP{\rho}{x}\right)_{i,j}^{\tau} \nonumber \\ && - \frac{\pi l}{y_{m}}\DP{\rho _{0}}{y}\hat \psi _{k,l}^{\tau} - \kappa \left[ \left(\frac{2\pi k}{x_{m}}\right)^{2}+ \left(\frac{\pi l}{y_{m}}\right)^{2}\right]\hat \rho _{k,l}^{\tau}. \end{eqnarray} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \newpage % 改ページ \section{使用モジュールとその他の設定} スペクトル変換と逆変換, 微分演算は SPMODEL ライブラリ (spml) の esc\_module に含まれる関数を用いて行う. フーリエ正弦および余弦変換, そ れらの逆変換の際の数値積分は台形公式を用いて行う. spml が下位で使用す る ISPACK の仕様から, 格子点数 $I,J$ は偶数で, かつ $I/2, J/2=2^{a}3^{b}5^{c}$ ($a, b, c$ は 0 または整数) でなければならな い. 非線形項の計算によって生じるエリアジングを防ぐため, 格子点数 $I,J$ と切断波数 $K,L$ は $I>3K,J>3K/2$ を満たすように与える. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \newpage % 改ページ \section{数値実験} 数値実験では水平方向に双極子型に並んだ正負の渦度強制により励起された内部 重力波の計算を行う. 渦度強制は以下のように与える. \begin{equation} F_{\zeta} = \DP{F_{v}}{x} = \DP{}{x}\left[A_{0}\exp \left(-\frac{(x-x_{0})^{2}+(y-y_{0})^{2}}{\sigma^{2}} \right) \sin \omega t \right] \end{equation} このとき強制項 $F_{\zeta}$ は水平方向にフーリエ展開, 鉛直方向にフーリエ 正弦展開される. \Deqref{渦度の増分}の右辺には \[ \frac{1}{I}\frac{1}{J} \sum _{i=0}^{I-1}\sum _{j=0}^{J-1} \exp\left(-\frac{2\pi i k x_{i}}{x_{m}}\right) \sin\left(\frac{\pi l y_{j}}{y_{m}}\right) F_{\zeta, i,j}^{\tau} \] が追加される. パラメータは\Dtabref{使用したパラメータの値}にまとめた値を 用いる. 格子点数 $I, J$ と切断波数 $K, L$ はそれぞれ $I=128, J=64, K=L=42$ と する. 時間格子間隔 $\Delta t$ は $0.1$ sec, 計算ステップ数は 4,000 ステップである. %\Dfigref{密度分布の時間変化}に計算された密度分布の時間変化を示した. 室 %内実験との対比がしやすいように $0\leq x\leq 60$ cm までの範囲を示してあ %る. \begin{table}[h] \begin{center} \begin{tabular}{cll} \hline パラメータ &\qquad\qquad& 数値 \\ \hline $\rho_{0}$ && 1.0 kgm${}^{-3}$ \\ $\DP{\rho _{0}}{y}$ && $-1.0$ kgm${}^{-2}$ \\ $g$ && 9.8 msec${}^{-2}$ \\ $\nu$ && 10${}^{-5}$ m${}^{2}$sec${}^{-1}$ \\ $\kappa$ && 10${}^{-9}$ m${}^{2}$sec${}^{-1}$ \\ $x_{m}$ && 2.0 m \\ $y_{m}$ && 1.0 m \\ $A_{0}$ && 1.0$\times 10^{-5}$ kgm${}^{-3}$ \\ $x_{0}$ && $x_{m}/2$ \\ $y_{0}$ && $y_{m}/2$ \\ $\sigma$ && 2.0$\times 10^{-2}$ m \\ \hline \end{tabular} \caption{使用したパラメータの値} \Dtablab{使用したパラメータの値} \end{center} \end{table} 渦度強制の振動数 $\omega$ は, ブラント--バイサラ振動数 \[ N \equiv \sqrt{\left|\frac{g}{\rho _{0}}\DP{\rho _{0}}{y}\right|} \] に対する比として与える(\Dtabref{渦度強制の振動数}). \begin{table}[h] \begin{center} \begin{tabular}{cll} \hline 実験番号 &\qquad\qquad& $\omega $ の値 \\ \hline case 1 && $0.3 N$ \\ case 2 && $0.6 N$ \\ case 3 && $0.9 N$ \\ case 4 && $1.2 N$ \\ \hline \end{tabular} \caption{渦度強制の振動数} \Dtablab{渦度強制の振動数} \end{center} \end{table} それぞれの $\omega $ を与えた場合の $t=400$ sec における密度偏差の分布を \Dfigref{case 1 の密度偏差分布}$\sim$\Dfigref{case 4 の密度偏差分布}に示 す. 内部重力波の伝播にともない, 密度偏差の分布は×印状のパターンを示す (\Dfigref{case 1 の密度偏差分布}). 強制の振動数をブラント--バイサラ振動 数に近づけるにつれて, ×印のパターンは閉じていく(\Dfigref{case 2 の密度 偏差分布}, \Dfigref{case 3 の密度偏差分布}). 強制の振動数がブラント--バ イサラ振動数を越えると波が発生しなくなる(\Dfigref{case 4 の密度偏差分布}). \begin{figure}[p] \begin{center} \Depsf[15cm][]{./figs/igwave1-omega0.3.ps} \vspace*{-1.0cm} \caption{$t=400$ sec における密度偏差の分布. case 1 の $\omega = 0.3 N$ の場合.} \Dfiglab{case 1 の密度偏差分布} \Depsf[15cm][]{./figs/igwave1-omega0.6.ps} \vspace*{-1.0cm} \caption{$t=400$ sec における密度偏差の分布. case 2 の $\omega = 0.6 N$ の場合.} \Dfiglab{case 2 の密度偏差分布} \end{center} \end{figure} \begin{figure}[p] \begin{center} \Depsf[15cm][]{./figs/igwave1-omega0.9.ps} \vspace*{-1.0cm} \caption{$t=400$ sec における密度偏差の分布. case 3 の $\omega = 0.9 N$ の場合.} \Dfiglab{case 3 の密度偏差分布} \Depsf[15cm][]{./figs/igwave1-omega1.2.ps} \vspace*{-1.0cm} \caption{$t=400$ sec における密度偏差の分布. case 4 の $\omega = 1.2 N$ の場合.} \Dfiglab{case 4 の密度偏差分布} \end{center} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \newpage % 改ページ \newpage \section{参考文献} \begin{description} \item GFD-online (酒井 敏, 飯澤 功, 荒巻 英治), 1997: 実験室の中の空と海, 内部重力波, http:\slash \slash www.gfd-dennou.org\slash arch\slash gfd-exp\slash gfd\_exp\slash exp\_j\slash doc\slash iw\slash guide01.htm. \item 木村竜治, 1983: 地球流体力学入門, 気象学のプロムナード 13, 東京堂出版, pp.247. \item Mowbray, D. E., Rarity, B. S. H., 1967: A theoretical and experimental investigation of the phase configuration ofinternal waves of small amplitude in a density stratified liquid. {\it J. Fluid Mech.}, {\bf 28}, 1--16. \item 竹広真一, 石岡圭一, 森川靖大, 小高正嗣, 石渡正樹, 林祥介, SPMODEL 開発グループ, 2004: 階層的地球流体力学スペクトルモデル集 (SPMODEL), http:\slash \slash www.gfd-dennou.org\slash arch\slash spmodel\slash , 地球流体電脳倶楽部. \end{description} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \newpage % 改ページ \section*{謝辞} % 謝辞. このまま利用する. \markright{謝辞} % ヘッダに表示 \addcontentsline{toc}{section}{謝辞} % 目次に表示 本資源は, 地球流体電脳倶楽部のインターネット上での学術知識の集積と活用 の実験の一環として \begin{center} http:\slash \slash www.gfd-dennou.org\slash arch\slash spmodel\slash \end{center} において公開されているものである (\copyright 地球流体電脳倶楽部スペクトルモデルプロジェクト spmodel@gfd-dennou.org 2002. ). 本資源は, 著作者の諸権利に抵触しない(迷惑をかけない)限りにおいて自由に 利用していただいて構わない. なお, 利用する際には今一度自ら内容を確かめ ることをお願いする(無保証無責任原則). 本資源に含まれる元資源提供者(図等の版元等を含む)からは, 直接的な形で の WEB 上での著作権または使用許諾を得ていない場合があるが, 勝手ながら, 「未来の教育」のための実験という学術目的であることをご理解いただけるも のと信じ, 学術標準の引用手順を守ることで諸手続きを略させていただいて いる. 本資源の利用者には, この点を理解の上, 注意して扱っていただける ようお願いする. 万一, 不都合のある場合には \begin{center} spmodel@gfd-dennou.org \end{center} まで連絡していただければ幸いである. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \end{document}