|
/ ~; U: X! U `2 y7 K% ~( \ 本文意在介绍发生在海洋中的动力过程的方程组,阅读本文需要基本的牛顿力学知识即可
* ^7 q7 R+ F; p+ z 动量方程E1-E3
4 M/ [. }4 M& Z4 T! j 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 % i4 r* X- o4 ?" A$ _4 r
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
1 |# E9 A, i r, s 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 0 X8 P. d- c& p- F+ R1 _# X
上述三个方程分别是动量方程的x、y、z分量形式 " U/ Y$ l+ F9 L6 ?) S( S
也可以写成矢量形式: ( I9 g* w6 v. X+ ~
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
$ y/ c0 p: W" v. S0 \ 以下我将逐个解释各项含义 ! s; k4 }$ [9 q2 `2 Z
等式左边为速度对时间的全导数,以E1为例,u为速度的x方向分量,u是(x,y,z,t)的函数
4 e0 H2 i, t: i0 M9 ] 等式右边包括重力、压强梯度力、科氏力、黏性力、湍应力、天体引潮力 * R B! O& t8 Z: a4 b+ I3 C' s
重力不用过多分析,仅存在于z方向 7 J; z6 E8 z% K1 D* D
压强梯度力:x方向为例,
9 @6 O7 i2 f/ X. b, @! 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
5 z6 z1 u( r6 V+ o- ^# y$ Y9 D% Y 科氏力: F=−2Ω×VF=-2\Omega\times V 5 @5 q- X7 L- `2 V* L+ P
Ω=2π/day=7.27÷105m/s\Omega=2\pi/day =7.27\div10^5 m/s
7 z' {% m, g! s* L& U Ω(0,Ωcosφ,Ωsinφ)\Omega (0,\Omega cos\varphi,\Omega sin\varphi) 9 t/ A6 l; \/ f) y; ^$ Z; Z
φ=latitude\varphi=latitude # F" V6 s, R( N, L; \& u7 Y. P
近似计算 2 }" {1 g6 ?; s! j* ]# A) `
Fx=fvF_x=fv
, Z0 T, X+ }8 W" i3 T8 h! [ Fy=−fuF_y=-fu & a# j: \/ T4 b% M) P
ff 为科氏系数 f=2Ωsinφf=2\Omega sin\varphi
, L1 K7 @1 `8 T8 f 黏性力为黏合系数与梯度的乘积,湍应力由湍流的脉冲造成的,天体引潮力过于复杂(与日月等天体有关,暂不介绍)
+ r1 h+ e" l0 x) X' J5 V% ~& r0 X E4 连续性方程
C8 X3 L6 u8 N% [( b' o: i- K ∂u/∂x+∂v/∂y+∂w/∂z=0\partial u/\partial x+\partial v/\partial y+\partial w/\partial z=0 % `" e' E& j/ Y& ?
Eularian观点:定点处观察经过的流体质量变化 ) z! d% ]9 W1 c; [6 S) d9 N
∂ρ/∂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 , A0 G. D! x$ b1 J8 \) }+ H, L
转化为Lagrange观点:跟踪流体微团
/ j% k: d: o v' h6 x [" o 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
4 M3 R, p' \0 C% x, i% Y E5-E6盐守恒、热守恒 $ [5 l6 H" ^4 f7 \; r# Y
E7 状态方程
- M& Y& J& j9 P2 a% Y* i" l ∂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
( S! e2 N6 p! U' k( h5 ]# k# ?4 R6 R1 E7 G3 M; A7 o9 B B
. h9 f' ~' |# V8 U- X6 k
6 a5 Q" | F% o6 W
9 J2 T$ B% W& x0 p+ u |