|
$ b. }7 y: S( a) h, |3 C 本文意在介绍发生在海洋中的动力过程的方程组,阅读本文需要基本的牛顿力学知识即可
* _% L2 `* A5 E2 _% p5 c! N1 w4 s 动量方程E1-E3 + E2 |( ~* t, l; r, H
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
7 o! t9 M( O0 ?" `9 f, N 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 * { A! F- I# q' h5 X$ Q& D' F7 r
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
7 L2 A6 s8 [ {* j9 L 上述三个方程分别是动量方程的x、y、z分量形式 & ]) z* ^( \, c
也可以写成矢量形式: 5 w( O; ^2 @0 z0 ?
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
+ q0 q" @& f" m# z' [. N8 ~9 @: ~8 } 以下我将逐个解释各项含义 % C% E# k0 Z J H; }* e+ O) y
等式左边为速度对时间的全导数,以E1为例,u为速度的x方向分量,u是(x,y,z,t)的函数
+ U! D: d5 D. _, g6 h+ K3 l 等式右边包括重力、压强梯度力、科氏力、黏性力、湍应力、天体引潮力 3 n, v2 y- T% y' f6 l( S8 }0 y
重力不用过多分析,仅存在于z方向
( U& \: {( k7 o- C) c" f3 t 压强梯度力:x方向为例, * V* O9 c& C' l* t- i& ?
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
6 H/ n4 T' D! K+ a 科氏力: F=−2Ω×VF=-2\Omega\times V * T+ X6 {9 l3 h8 p! a; B
Ω=2π/day=7.27÷105m/s\Omega=2\pi/day =7.27\div10^5 m/s 7 X# b# f/ e4 b4 X( @2 G
Ω(0,Ωcosφ,Ωsinφ)\Omega (0,\Omega cos\varphi,\Omega sin\varphi) 0 m( R& U% c9 n/ o+ r
φ=latitude\varphi=latitude + Z) A N/ z- A0 G
近似计算 & b2 _; ^/ H2 Y& }- [
Fx=fvF_x=fv / S* j$ ^; r- `8 B4 w# O
Fy=−fuF_y=-fu 4 V8 A' H1 v8 R
ff 为科氏系数 f=2Ωsinφf=2\Omega sin\varphi ' O' c* e6 C9 u8 X! M0 `! l) H
黏性力为黏合系数与梯度的乘积,湍应力由湍流的脉冲造成的,天体引潮力过于复杂(与日月等天体有关,暂不介绍)
4 \* i, U; i5 _2 \; ? ^ E4 连续性方程 / Z: C" Z6 y+ u
∂u/∂x+∂v/∂y+∂w/∂z=0\partial u/\partial x+\partial v/\partial y+\partial w/\partial z=0 # u/ e% F; i2 l# A: o
Eularian观点:定点处观察经过的流体质量变化 / G$ g% Q" U" _* N/ k* G5 |
∂ρ/∂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
8 F5 Q" }. `1 x- V$ u6 d 转化为Lagrange观点:跟踪流体微团 ! C ~( U( ]0 L! L" r
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
5 q1 a; O$ G5 X: i) q; z E5-E6盐守恒、热守恒 / o8 E4 D( G- Q1 |
E7 状态方程 * |* I* F. B% T3 Z+ w
∂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 - N0 J* p3 E/ f$ F4 j7 m5 Z
8 ~3 T" [. B7 W
2 x) j' x. j1 u. O0 |& I
- v8 p+ D/ [" N# e. H
) J" G) a- s" K( j9 i |