|
# R0 F6 A/ z8 @
本文意在介绍发生在海洋中的动力过程的方程组,阅读本文需要基本的牛顿力学知识即可
# u# a% G0 @) r) S* `, Z) A' ^6 o2 c/ m 动量方程E1-E3 * U) V8 u4 D! Y# R3 f1 K6 Q+ S4 M; D
E1:∂u/∂t+u∂u/∂x+v∂u/∂y+w∂u/∂z=−1/ρ⋅∂p/∂x+fv+υΔu+∂(AH∂u/∂x)/∂x+∂(AH∂u/∂y)/∂y+∂(Az∂u/∂z)/∂z+FxE1:\partial u/\partial t+u\partial u/\partial x+v\partial u/\partial y+w\partial u/\partial z=-1/\rho\cdot\partial p/\partial x+fv+\upsilon\Delta u+\partial (A_H \partial u/\partial x)/\partial x+\partial (A_H \partial u/\partial y)/\partial y+\partial (A_z \partial u/\partial z)/\partial z+F_x
* x* F( J5 n5 s) G9 h3 b: {7 ` E2:∂v/∂t+u∂v/∂x+v∂v/∂y+w∂v/∂z=−1/ρ⋅∂p/∂y−fu+υΔv+∂(AH∂v/∂x)/∂x+∂(AH∂v/∂y)/∂y+∂(Az∂v/∂z)/∂z+FyE2:\partial v/\partial t+u\partial v/\partial x+v\partial v/\partial y+w\partial v/\partial z=-1/\rho\cdot\partial p/\partial y-fu+\upsilon\Delta v+\partial (A_H \partial v/\partial x)/\partial x+\partial (A_H \partial v/\partial y)/\partial y+\partial (A_z \partial v/\partial z)/\partial z+F_y
: H) w. {% j# f* l% L7 Y E3:∂w/∂t+u∂w/∂x+v∂w/∂y+w∂w/∂z=g−1/ρ⋅∂p/∂z+υΔw+∂(AH∂w/∂x)/∂x+∂(AH∂w/∂y)/∂y+∂(Az∂w/∂z)/∂z+FzE3:\partial w/\partial t+u\partial w/\partial x+v\partial w/\partial y+w\partial w/\partial z=g-1/\rho\cdot\partial p/\partial z+\upsilon\Delta w+\partial (A_H \partial w/\partial x)/\partial x+\partial (A_H \partial w/\partial y)/\partial y+\partial (A_z \partial w/\partial z)/\partial z+F_z - j3 }9 G/ @+ R- w t7 x8 j; H
上述三个方程分别是动量方程的x、y、z分量形式 3 C, l" Y# ?. K2 l/ J0 k/ K" z
也可以写成矢量形式: 9 C4 x, u: \# p q- C
dV¯/dt=g−1/ρ⋅(hamilton)P+Ω×V¯+υΔ(hamilton)barV+Ft+Frd\bar{V}/dt=g-1/\rho\cdot(hamilton)P+\Omega \times \bar{V}+\upsilon\Delta(hamilton)bar{V}+F_t+F_r 0 k& n' F9 g2 |0 {
以下我将逐个解释各项含义 ( [% T( y( e/ r+ e
等式左边为速度对时间的全导数,以E1为例,u为速度的x方向分量,u是(x,y,z,t)的函数 ' ]8 T$ K$ M; d# V+ u2 k2 R# Q
等式右边包括重力、压强梯度力、科氏力、黏性力、湍应力、天体引潮力
4 G/ S# i6 p% w1 i$ B# O 重力不用过多分析,仅存在于z方向
% N" E/ j5 d* v# \* L 压强梯度力:x方向为例, # f2 p) q$ Z0 \
a=F/m=(p−(p+δp))⋅δyδz/ρ⋅δxδyδz=−1/ρ⋅∂p/∂xa=F/m=(p-(p+\delta p))\cdot\delta y\delta z/\rho\cdot \delta x\delta y\delta z=-1/\rho\cdot \partial p/\partial x / U7 k8 A$ |% x. p% p2 ?7 R
科氏力: F=−2Ω×VF=-2\Omega\times V
* c# I7 h, P! c5 p3 D% u* W Ω=2π/day=7.27÷105m/s\Omega=2\pi/day =7.27\div10^5 m/s
% E4 E G2 h% H. \3 W I Ω(0,Ωcosφ,Ωsinφ)\Omega (0,\Omega cos\varphi,\Omega sin\varphi)
( Y5 b j3 Z+ f" V φ=latitude\varphi=latitude B, h8 h$ Z9 C7 }. _8 S
近似计算 & J" W& i. f; C8 L! J% {
Fx=fvF_x=fv
8 Q C: e' h9 u. i+ m2 m k Fy=−fuF_y=-fu " V A6 k& e! q0 |
ff 为科氏系数 f=2Ωsinφf=2\Omega sin\varphi
4 ?4 U. f; @* [ 黏性力为黏合系数与梯度的乘积,湍应力由湍流的脉冲造成的,天体引潮力过于复杂(与日月等天体有关,暂不介绍)
2 P* x! I' ]2 r1 x. E E4 连续性方程 ! P" D% m- P0 s* @% \ ~
∂u/∂x+∂v/∂y+∂w/∂z=0\partial u/\partial x+\partial v/\partial y+\partial w/\partial z=0
% I' K2 D: m+ E3 e. S; p8 n5 Z Eularian观点:定点处观察经过的流体质量变化 6 \- I/ v1 C4 G2 Q
∂ρ/∂t+(∂(ρu)∂x+∂(ρv)/∂y+∂(ρw)/∂z=0\partial \rho/\partial t+(\partial(\rho u)\partial x+\partial(\rho v)/\partial y+\partial (\rho w)/\partial z=0
& |- L" u! [! G7 \* t 转化为Lagrange观点:跟踪流体微团
) \7 D; j2 N6 t. _0 d% q 1/ρDρ/Dt+(∂u/∂x+∂v/∂y+∂w/∂z)=01/\rho D\rho /Dt +(\partial u/\partial x+\partial v/\partial y+\partial w/\partial z)=0 8 @0 D/ Z4 L$ ^ ]
E5-E6盐守恒、热守恒 ! w0 s$ \0 o& x/ L9 j
E7 状态方程
2 y+ s. S+ n0 X( C8 v ∂s/∂t+u∂s/∂x+v∂s/∂y+w∂s/∂z=kDΔs+∂(kH∂s/∂x)/∂x+∂(kH∂s/∂y)/∂y+∂(kH∂s/∂z)/∂z\partial s/\partial t+u\partial s/\partial x+v\partial s/\partial y+w\partial s/\partial z=k_D\Delta s+\partial(k_H \partial s/ \partial x)/\partial x+\partial(k_H \partial s/ \partial y)/\partial y+\partial(k_H \partial s/ \partial z)/\partial z 3 {- R% l; ^$ g0 D
* A- a; E6 v, H2 Y! [; \
; n9 q, ?# y' G0 G- K, K5 {% }
& ]2 [( w/ L) d; ?8 ~5 D9 Z4 o2 a X
|