エクマン境界層のソースを表示

{{参照方法|date=2024年3月}}
'''エクマン境界層'''（エクマンきょうかいそう、{{Lang-en-short|Ekman layer}}）は、[[コリオリの力]]と{{仮リンク|渦粘性|en|eddy viscosity}}がつりあう[[流体]]に見られる[[境界層]]である。[[ヴァン・ヴァルフリート・エクマン]]によって初めて記述された。

==エクマン境界層の一般解==
二次元回転平面流体において、圧力傾度力、コリオリの力、渦粘性が釣り合っている状況を考える。すなわち、平面流れ（''u'' , ''v''）が次の微分方程式系に従うと考える。

:<math>
\begin{align}
  -fv &= -\frac{1}{\rho_o} \frac{\partial p}{\partial x}+K_z \frac{\partial^2 u}{\partial z^2}, \\
   fu &= -\frac{1}{\rho_o} \frac{\partial p}{\partial y}+K_z \frac{\partial^2 v}{\partial z^2}
\end{align}
</math>

ここで[[地衡流]]（''u<sub>g</sub>'' , ''v<sub>g</sub>''）を考えると、

:<math>
\begin{align}
  -fv_g &= -\frac{1}{\rho_o} \frac{\partial p}{\partial x}\\
   fu_g &= -\frac{1}{\rho_o} \frac{\partial p}{\partial y}
\end{align}
</math>

非地衡流成分 ''u' ''= ''u'' - ''u<sub>g</sub>'' , ''v' ''= ''v'' - ''v<sub>g</sub>''が従う方程式は

:<math>
\begin{align}
  -fv^\prime &= K_z \frac{\partial^2 u^\prime}{\partial z^2}, \\
   fu^\prime &= K_z \frac{\partial^2 v^\prime}{\partial z^2}
\end{align}
</math>

以下では、非地衡流成分 ''u' '', ''v' ''を省略して ''u'' , ''v''と記す。

ここで複素速度 ''V''=''u''+''iv''（''i''は[[虚数単位]]）を用いると、次の二階連立偏微分方程式にまとめることができる。

:<math>
 K_z\frac{\partial^2 V}{\partial z^2}=ifV
</math>

この微分方程式の一般解は

:<math>
V =A\exp\left[\left(\frac{if}{K_z}\right)^{\frac{1}{2}}z\right]
          +B\exp\left[-\left(\frac{if}{K_z}\right)^{\frac{1}{2}}z\right]
</math>

ここで<math>h_E\equiv\sqrt{2K_z/f}</math>と置くと、<math>\sqrt{i}=(1+i)/\sqrt{2}</math>であるから

:<math>
V =A\exp\left[(1+i)\frac{z}{h_E}\right]+B\exp\left[-(1+i)\frac{z}{h_E}\right]
</math>

この''h<sub>E</sub>''は、エクマン境界層の厚さの目安のひとつになっている。
===海上を吹く風によって形成されるエクマン境界層===
海上を吹く風が海面に及ぼす風応力によって海面に運動がおこる。この運動が渦粘性によって次第に下層まで伝わり、エクマン境界層が生じる。風が長期間一定であると仮定したときのエクマン境界層について記述する。これは上記の一般解に2つの境界条件を課すと解くことができる。

境界条件のひとつは、海洋表層における力が風応力による摩擦によって生じるという仮定から導かれる。

:<math>
\rho K_z\frac{\partial V}{\partial z}\Big|_{z=0}=\tau,\quad \tau=\tau_x+i\tau_y
</math>

もうひとつは、いま考えている層が海面に沿った境界層であることから導かれる。海はエクマン層に比べて無限に深いとし、海底ではエクマン層に由来する流れが消えていると仮定する。

:<math>
V^\prime(-\infty)=0
</math>

2つ目の境界条件から''B''=0がわかる。次に1つ目の境界条件から、

:<math>
\begin{align}
    \rho K_z\frac{\partial V}{\partial z}\Big|_{z=0}=A\rho K_z\frac{1+i}{h_E}=\tau,\\
    A=\frac{\tau h_E}{\rho K_z(1+i)}=\frac{\tau h_E(1-i)}{2\rho K_z},\\
\end{align}
</math>

これを一般解に代入して整理し、実部と虚部に分けると、

:<math>
\begin{align}
u=\frac{h_E}{2\rho K_z}\left[(\tau_x+\tau_y)\cos\frac{z}{h_E}+(\tau_x-\tau_y)\sin\frac{z}{h_E}\right],\\
v=\frac{h_E}{2\rho K_z}\left[(-\tau_x+\tau_y)\cos\frac{z}{h_E}+(\tau_x+\tau_y)\sin\frac{z}{h_E}\right]
\end{align}
</math>

特に''z''=0において、地衡流を除いた風応力のみによる吹送流をみると、

:<math>
u_0=\frac{h_E}{2\rho K_z}(\tau_x+\tau_y)\qquad,v_0=\frac{h_E}{2\rho K_z}(-\tau_x+\tau_y)
</math>

すなわち海表面における流れの速度ベクトルは、風応力ベクトルの右45度の角をなす（北半球）。

==参考文献==
*{{Cite|和書|title=気象力学通論 |author=小倉義光 |publisher= |year= |isbn= |page= }}

==関連項目==
*{{仮リンク|エクマン螺旋|en|Ekman spiral}}
*[[エクマン輸送]]

{{Normdaten}}
{{DEFAULTSORT:えくまんきようかいそう}}
[[Category:気象]]
[[Category:水環境]]
[[Category:大気力学]]
[[Category:境界層]]
[[Category:エポニム]]
[[Category:ヴァン・ヴァルフリート・エクマン]]