|
! S6 r3 ?- V3 L1 z5 `& ~& E 本文意在介绍发生在海洋中的动力过程的方程组,阅读本文需要基本的牛顿力学知识即可 % T" i2 R6 U' D8 R5 ^. q
动量方程E1-E3
" x4 R1 ^0 s) T5 P* ?: u8 p* a0 W 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 2 R7 B0 b9 L" d, Q7 b
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 ! E, ~. T% T* 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 % p* e2 D7 b5 g$ Q/ e
上述三个方程分别是动量方程的x、y、z分量形式 ' w7 E9 X6 c: h0 c( Y8 [* U+ ]- N
也可以写成矢量形式:
$ B G6 o q; p! B+ k 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 # Z+ e7 d3 F% _8 P
以下我将逐个解释各项含义
. Y6 [$ o; I5 k! _/ j1 I/ ` 等式左边为速度对时间的全导数,以E1为例,u为速度的x方向分量,u是(x,y,z,t)的函数
' b+ w% X2 Q3 U; i9 S u" @) t 等式右边包括重力、压强梯度力、科氏力、黏性力、湍应力、天体引潮力 % { o- G, y. x* y l0 D k
重力不用过多分析,仅存在于z方向
3 y. Y- L* Z* z# k1 {) `$ J 压强梯度力:x方向为例,
$ p/ `3 J' R5 N% ^! s2 Y0 B' T 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 * p E0 U7 k2 r7 s; y' s" t1 a8 Z
科氏力: F=−2Ω×VF=-2\Omega\times V & `* p% g* b4 c+ W6 C
Ω=2π/day=7.27÷105m/s\Omega=2\pi/day =7.27\div10^5 m/s # U k8 |( x/ J
Ω(0,Ωcosφ,Ωsinφ)\Omega (0,\Omega cos\varphi,\Omega sin\varphi) $ R$ ^; {" M& f
φ=latitude\varphi=latitude
( _: |- L% y( p5 T, ~ 近似计算 / z G2 y9 [' `
Fx=fvF_x=fv
( R/ m% l3 p% s$ L2 H" [ Fy=−fuF_y=-fu
- Z i g) i3 X$ b ff 为科氏系数 f=2Ωsinφf=2\Omega sin\varphi / \. n, ?- Y$ i, @6 x4 o
黏性力为黏合系数与梯度的乘积,湍应力由湍流的脉冲造成的,天体引潮力过于复杂(与日月等天体有关,暂不介绍) 4 X7 d" _3 l" b( [
E4 连续性方程 $ n4 X) o$ C. b" c5 K9 s) `4 p
∂u/∂x+∂v/∂y+∂w/∂z=0\partial u/\partial x+\partial v/\partial y+\partial w/\partial z=0 5 U. J; E/ n: d. b" O1 @, N* B
Eularian观点:定点处观察经过的流体质量变化
3 P5 ^) `3 v$ S6 }' z' l3 } ∂ρ/∂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
( p8 S* s7 J# k* G 转化为Lagrange观点:跟踪流体微团
: C {$ z- A( \* ]9 e, d ^! w# Z 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
7 f$ _) o a$ j5 Q1 G E5-E6盐守恒、热守恒
' I0 T4 m. q2 Q; p3 D, J E7 状态方程
+ Z5 [2 E p6 y' G+ ?; U ∂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 1 E, C3 i0 r9 h! C5 v; J( J" @& G
. i) P( y3 @8 a% S5 p6 V
2 h; o) O, n0 g' ] s/ n9 ^+ u: Q9 @0 E
5 R* _ ]# M/ s& @' w$ Y: i |