|
9 y- I5 v; n- G* ~7 S& u
本文意在介绍发生在海洋中的动力过程的方程组,阅读本文需要基本的牛顿力学知识即可 6 m' j. E X, l! u' a y
动量方程E1-E3 ; Y- Q, G o1 Q) t
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
9 i- H+ \0 b: a& y# [ 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 ) ~% s: q2 A4 ^1 E9 a+ J: k
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 : M/ X1 z+ F8 O& w. Q$ _' m
上述三个方程分别是动量方程的x、y、z分量形式
+ K4 a3 P/ n! \3 f+ |; Z 也可以写成矢量形式:
/ g, u0 l. R. c/ c+ s. 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
& x7 U% [% G/ z" _ 以下我将逐个解释各项含义 8 F% |# p/ X. [, V% z5 X
等式左边为速度对时间的全导数,以E1为例,u为速度的x方向分量,u是(x,y,z,t)的函数
' I5 B3 m9 [) s1 @" z 等式右边包括重力、压强梯度力、科氏力、黏性力、湍应力、天体引潮力
: Q" q x& c: C3 N3 a 重力不用过多分析,仅存在于z方向 7 C- E, W6 D- `; y" e
压强梯度力:x方向为例, - L9 d* m- [- I1 _5 Z& `( v2 l5 O
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 0 T: [, j; D; X5 [3 L! ?
科氏力: F=−2Ω×VF=-2\Omega\times V ; K) W8 h2 F9 U
Ω=2π/day=7.27÷105m/s\Omega=2\pi/day =7.27\div10^5 m/s
- i5 I6 v* _% g: W" d6 o Ω(0,Ωcosφ,Ωsinφ)\Omega (0,\Omega cos\varphi,\Omega sin\varphi) * |7 N0 H+ C. ]8 t, ?
φ=latitude\varphi=latitude - [4 I& B* ^% {( `0 H% t. @8 F# y
近似计算
8 [" i( q% L& m0 B! E Fx=fvF_x=fv
7 t8 G) T2 x3 i* `* j Fy=−fuF_y=-fu 8 Y! i% M/ v" q1 |0 u
ff 为科氏系数 f=2Ωsinφf=2\Omega sin\varphi ' t% g& m4 R' U" z
黏性力为黏合系数与梯度的乘积,湍应力由湍流的脉冲造成的,天体引潮力过于复杂(与日月等天体有关,暂不介绍)
( X8 [$ f" k8 h9 ^5 r& o/ f E4 连续性方程 ( ]1 } p9 K3 `- Y- e
∂u/∂x+∂v/∂y+∂w/∂z=0\partial u/\partial x+\partial v/\partial y+\partial w/\partial z=0
2 r/ Z* C; W# d( S& c Eularian观点:定点处观察经过的流体质量变化
. e# d5 F% o1 N$ j8 t6 @# w ∂ρ/∂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
) ]1 Y+ P4 m0 {+ D1 ~% A 转化为Lagrange观点:跟踪流体微团 : G5 _9 v. Y [/ U5 }1 ?
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
3 L: k5 s3 }' v8 F& y6 X0 F) `: ~ E5-E6盐守恒、热守恒 2 O% h3 ?2 u. u( T: N
E7 状态方程
. T, m: R' P1 z' N5 ^3 ?5 [ ∂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 ( z X1 r! O8 k1 k" v3 o# I
0 d( k# Q" Z' a
) ?& i7 q( ~6 B$ o2 [
& t2 w* D6 F* L
3 P9 r. v& W' s* Z+ v4 b |