\documentclass[a4j,12pt,dvipdfmx]{jarticle}
\usepackage[dvipdfmx]{graphicx}
\usepackage{D6math,D6misc}
\usepackage{amsmath}

\Dnoparindent
\Dparskip

\begin{document}

\begin{center}
  \Large\bf
  移流テストの定式化
\end{center}

\begin{flushright}
  2026 年 6 月 4 日 竹広 真一
\end{flushright}

ここでは DCPAM のセミラグランジアンスキームのテストのための
移流問題の定式化を行う.

\section{非発散流れ場}

テストを行う流れ場は経度緯度 $\sigma$ 座標系で非発散であることを
要請する.
%
\begin{align}
  \frac{1}{a\cos\varphi}\DP{u}{\lambda}
  + \frac{1}{a\cos\varphi}\DP{}{\varphi}( v\cos\varphi )
  + \DP{w}{\sigma} = 0,
\end{align}
%
ここで $\lambda$, $\varphi$, $\sigma$ は
経度, 緯度, 鉛直 $\sigma$ 座標であり,
$(u,v,w)$ はそれぞれの速度成分である. 
速度場の満たすべき自然な境界条件として
%
\begin{align}
  u|_{\lambda=0} = u|_{\lambda=2\pi}, \quad
  v=0 \ \mbox{at} \ \varphi=\pm \frac{\pi}{2},
  \quad
  w=0 \ \mbox{at} \ \sigma=0,1.
\end{align}

一方, パッシブスカラー $q(\lambda,\varphi,\sigma,t)$ の移流は
%
\begin{align}
  \DP{q}{t}
  + \frac{1}{a\cos\varphi}\DP{}{\lambda}(qu)
  + \frac{1}{a\cos\varphi}\DP{}{\varphi}(qv\cos\varphi )
  + \DP{}{\sigma}(qw) = 0
\end{align}
%
に従う. 非発散条件を用いて変形すると, 
%
\begin{align}
  \DP{q}{t}
  + \frac{u}{a\cos\varphi}\DP{q}{\lambda}
  + \frac{v}{a} \DP{q}{\varphi} + w\DP{q}{\sigma}=0.
\end{align}

% ----------------------------------------------------------------------
%
%
\section{帯状方向の速度場}

$v=0$ の場合の帯状流れ場の例を考える.
2 次元流れなので非発散条件を満たすよう流線関数を次のように導入できる.
%
\begin{align}
  u = -\DP{\Phi}{\sigma}, \quad w = \frac{1}{a\cos\varphi}\DP{\Phi}{\lambda}.
\end{align}
%
このとき, パッシブスカラーの移流を表す式は,
%
\begin{align}
  \DP{q}{t} -\frac{1}{a\cos\varphi}\DP{\Phi}{\sigma}\DP{q}{\lambda}
            + \frac{1}{a\cos\varphi}\DP{\Phi}{\lambda}\DP{q}{\sigma}
  = \DP{q}{t} + \frac{1}{a\cos\varphi}J_{\lambda,\sigma}(\Phi,q) = 0.
\end{align}
%
ここで $\Ddsty J_{\lambda,\sigma}(f,g) =
\DP{f}{\lambda}\DP{g}{\sigma}-\DP{f}{\sigma}\DP{g}{\lambda}$ は
$(\lambda,\sigma)$ 系でのヤコビアンである. 
$q$ が流線関数 $\Phi$ と同じ分布の場合,
ヤコビアンが 0 となり分布が時間変化しないことに注意されたい. 


境界条件を流線関数として
%
\begin{align}
  \Phi(\lambda,\sigma)
  &= \Phi_0 \sin(n\pi\sigma)\cos(m\lambda)\cos\varphi,
  \\
  u(\lambda,\sigma)
  &= - n\pi\Phi_0 \cos(n\pi\sigma)\cos(m\lambda)\cos\varphi,
  \\
  w(\lambda,\sigma)
  &= -\frac{m\Phi_0}{a} \sin(n\pi\sigma)\sin(m\lambda).
\end{align}

% ----------------------------------------------------------------------
%
%
\section{子午面内の流れ場}

子午面内の流れは, 帯状平均した質量保存則を満たさねばならない.
%
\begin{align}
  \frac{1}{a\cos\varphi}\DP{}{\varphi}( v\cos\varphi )
  + \DP{w}{\sigma} = 0,
\end{align}
%
ここで $v,w$ はそれぞれ速度の緯度および鉛直($\sigma$)成分である. 
これを満たすべく, ここでの流線関数 $\Psi$ を次のように定義する.
%
\begin{align}
  v = - \frac{1}{\cos\varphi}\DP{\Psi}{\sigma}, \quad
  w = \frac{1}{a\cos\varphi}\DP{\Psi}{\varphi}.
\end{align}
%
パッシブスカラー $q(\varphi,\sigma,t)$ の移流は, 
%
\begin{align}
  \DP{q}{t} -\frac{1}{a\cos\varphi}\DP{\Psi}{\sigma}\DP{q}{\varphi}
            + \frac{1}{a\cos\varphi}\DP{\Psi}{\varphi}\DP{q}{\sigma}
  = \DP{q}{t} + \frac{1}{a\cos\varphi}J_{\varphi,\sigma}(\Psi,q) = 0.
\end{align}
%
ここで $\Ddsty J_{\varphi,\sigma}(f,g)=\DP{f}{\varphi}\DP{g}{\sigma}$ は
$(\varphi,\sigma$ 系でのヤコビアンである.
$q$ が流線関数と同じ分布の場合, ヤコビアンが 0 となり
分布が時間変化しないことに注意されたい. 

境界条件を満たす流れ場として $\Ddsty \varphi=\pm\frac{\pi}{2}$ で
特異でない解を選択すると, 例えば, 
%
\begin{align}
  \Psi &= \sin(n\pi\sigma) \cos^2\varphi P_l'(\sin\varphi),
  \\
  v &= - n \pi \cos(n\pi\sigma) \cos^2\varphi P_l'(\sin\varphi),
  \\
  w &= - \frac{l(l+1)}{a}\sin(n\pi\sigma) P_l(\sin\varphi).
\end{align}
%
ただし $P_l$ および $P_l'$ は $l$ 次のルジャンドル多項式およびその微分
$\Ddsty P_l'=\DD{P_l(x)}{x}$ である.

この流れでパッシブスカラーが移流されると元の分布を保てないが,
流れが時間的に振動すると, その周期でもとに戻されるだろう.
すなわち,
%
\begin{align}
  \Psi &= \sin(n\pi\sigma) \cos^2\varphi P_l'(\sin\varphi)\cos\left(\frac{2\pi t}{T}\right),
  \\
  v &= - n \pi \cos(n\pi\sigma) \cos^2\varphi P_l'(\sin\varphi)\cos\left(\frac{2\pi t}{T}\right),
  \\
  w &= - \frac{l(l+1)}{a}\sin(n\pi\sigma) P_l(\sin\varphi)\cos\left(\frac{2\pi t}{T}\right).
\end{align}
%
ここで $T$ は流れの振動周期である.


\begin{itemize}
\item $n=1$, $l=1$ の場合 : $P_1=x$, $P_1'=1$ なので\\
  \begin{align}
    \Psi &= \sin(\pi\sigma) \cos^2\varphi\cos\left(\frac{2\pi t}{T}\right),
    \\
    v &= - \pi \cos(\pi\sigma) \cos\varphi\cos\left(\frac{2\pi t}{T}\right),
    \\
    w &= - \frac{2}{a}\sin(\pi\sigma) \sin\varphi\cos\left(\frac{2\pi t}{T}\right).
  \end{align}

\item $n=1$, $l=2$ の場合 : $\Ddsty P_2=\frac{1}{2}(3x^2-1)$, 
  $P_2'= 3x$ なので
  \begin{align}
    \Psi &= \sin(\pi\sigma) 3\sin\varphi\cos^2\varphi \cos\left(\frac{2\pi t}{T}\right),
    \\
    v &= - \pi \cos(\pi\sigma) 3\sin\varphi\cos\varphi\cos\left(\frac{2\pi t}{T}\right),
    \\
    w &= - \frac{3}{a}\sin(\pi\sigma) (3\sin^2\varphi-1)\cos\left(\frac{2\pi t}{T}\right).
  \end{align}
  
\end{itemize}

\begin{figure}
  \centering
  \includegraphics[width=0.6\textwidth]{./figures/meridflow.png}
  \caption{
    子午面の流れの例. 上が $n=1$, $l=1$ の場合, 下が $n=1$, $l=2$ の場合. 
  }
\end{figure}

\end{document}

%%% Local Variables:
%%% mode: japanese-latex
%%% TeX-master: t
%%% End:
