|
# S& ]3 n2 N2 l- P
本文意在介绍发生在海洋中的动力过程的方程组,阅读本文需要基本的牛顿力学知识即可
8 X+ a, X$ c7 L- S# O+ \# V- Z* }7 r 动量方程E1-E3
: o# L4 p, j" F. k- ~9 u, q 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
6 f J8 V+ A, q/ A 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 + n X2 G, d6 ` l
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
) G$ ~- u+ Q6 \1 ? 上述三个方程分别是动量方程的x、y、z分量形式
e# T1 E0 F; G 也可以写成矢量形式: 5 A& v y$ F& q: t
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 " O4 F9 V4 ], F$ d+ t
以下我将逐个解释各项含义 ( e: q" u1 K, h; _: u$ Q1 y
等式左边为速度对时间的全导数,以E1为例,u为速度的x方向分量,u是(x,y,z,t)的函数 - J8 l; ~; m3 g
等式右边包括重力、压强梯度力、科氏力、黏性力、湍应力、天体引潮力 ! T; {- J( i# a
重力不用过多分析,仅存在于z方向
; N# p( m6 k( K7 m3 I 压强梯度力:x方向为例, 2 @$ b( w }! {5 F* F# ]3 X0 l
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
! v, O! d$ d' s& p 科氏力: F=−2Ω×VF=-2\Omega\times V
9 t* G/ t6 t7 H! G Ω=2π/day=7.27÷105m/s\Omega=2\pi/day =7.27\div10^5 m/s 3 B2 s/ O, k7 z+ s
Ω(0,Ωcosφ,Ωsinφ)\Omega (0,\Omega cos\varphi,\Omega sin\varphi)
/ B! ~' D! }8 N/ o* b- }; Y$ } φ=latitude\varphi=latitude
# C: z) t& s0 h; z( j 近似计算 , H( ]1 ?+ o; ^6 {& o" J& g
Fx=fvF_x=fv
' @, ?' Y1 u2 m- } Fy=−fuF_y=-fu
2 N8 _' e( ?8 ^9 B. A: W9 N% c( [ ff 为科氏系数 f=2Ωsinφf=2\Omega sin\varphi
' F2 [! V/ t# X, } 黏性力为黏合系数与梯度的乘积,湍应力由湍流的脉冲造成的,天体引潮力过于复杂(与日月等天体有关,暂不介绍) 2 g5 _' l; _9 w: ^, f) u5 V5 [
E4 连续性方程 / w/ i3 T8 C* {$ N" c
∂u/∂x+∂v/∂y+∂w/∂z=0\partial u/\partial x+\partial v/\partial y+\partial w/\partial z=0 4 m* `6 D2 e# ^5 X: F# V1 e( t9 |
Eularian观点:定点处观察经过的流体质量变化 - `* W5 y& |/ U/ I ]
∂ρ/∂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 y, y k" S6 }
转化为Lagrange观点:跟踪流体微团
$ E% d. m$ ]# m 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
6 Y4 y* a$ M7 ]$ V/ X$ d( i E5-E6盐守恒、热守恒
$ z* @# _$ M$ L% V) h( e2 [- S& c' T E7 状态方程 - r9 ]: ~, ~2 | H
∂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
5 G. k* q3 v/ p/ R) u3 B' n. E# {$ ~' i1 P Y0 {, Q3 t
) N4 U. ? Y; b' a* s! y% R; |+ X8 I6 j" f; P4 \8 a0 t
! j/ c! x7 k6 J; w; f% v9 i |