收藏本站 劰载中...网站公告 | 吾爱海洋论坛交流QQ群:835383472

海洋动力学 -什么叫海洋动力舰

[复制链接]
# {$ f0 X+ K, @6 K) a; b, q2 e

本文意在介绍发生在海洋中的动力过程的方程组,阅读本文需要基本的牛顿力学知识即可

' H3 H6 O6 k+ t q

动量方程E1-E3

. H, X3 u% Y0 \) L! N

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

, A8 z3 _* e* 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

8 {- L$ Q8 n% c, ~0 Q7 `

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

5 Z- h: j5 F- S$ l- z( d% b

上述三个方程分别是动量方程的x、y、z分量形式

4 U3 G3 t1 ]# A

也可以写成矢量形式:

1 _. z3 I1 x* r& p; `4 R

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

' u1 t( A# _) x

以下我将逐个解释各项含义

/ w2 k& X0 g3 o4 A2 Y- ~( {( S

等式左边为速度对时间的全导数,以E1为例,u为速度的x方向分量,u是(x,y,z,t)的函数

$ d- X" w+ p6 G; l; X! m# b

等式右边包括重力、压强梯度力、科氏力、黏性力、湍应力、天体引潮力

# `/ F4 x+ R( F2 ~2 Q2 Y

重力不用过多分析,仅存在于z方向

) Q7 E. h9 X8 S) t* g3 |, z) w7 Z

压强梯度力:x方向为例,

. d% ?$ S U4 T+ n

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

/ Q1 h9 ?1 |4 d; v; m

科氏力: F=−2Ω×VF=-2\Omega\times V

, e9 {5 f, _; k2 @( ~" D

Ω=2π/day=7.27÷105m/s\Omega=2\pi/day =7.27\div10^5 m/s

$ z1 b/ Q6 `- [

Ω(0,Ωcosφ,Ωsinφ)\Omega (0,\Omega cos\varphi,\Omega sin\varphi)

8 p7 _- D8 G& ]% e2 z

φ=latitude\varphi=latitude

: Y" @" h+ [; Y& n2 y

近似计算

; E$ m3 d& U7 M5 j# M

Fx=fvF_x=fv

( {; G- E3 `0 h. `! E. B* L) g- d

Fy=−fuF_y=-fu

0 N* Z/ A9 O: h/ S9 z& U3 K

ff 为科氏系数 f=2Ωsinφf=2\Omega sin\varphi

9 r" Z5 ]2 d' l

黏性力为黏合系数与梯度的乘积,湍应力由湍流的脉冲造成的,天体引潮力过于复杂(与日月等天体有关,暂不介绍)

a. M l/ o, t+ U) E% x& X

E4 连续性方程

, F. m% x: {. W, S

∂u/∂x+∂v/∂y+∂w/∂z=0\partial u/\partial x+\partial v/\partial y+\partial w/\partial z=0

! U! ?+ M: B4 [

Eularian观点:定点处观察经过的流体质量变化

$ l& x- Y' z& d

∂ρ/∂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

5 H. G; G* C! X9 _ V' E$ q3 F/ Y

转化为Lagrange观点:跟踪流体微团

- u7 f1 O6 H, S# @: H7 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

& u( |$ N/ A# c3 W

E5-E6盐守恒、热守恒

* a$ k4 d. s N

E7 状态方程

8 g. p6 S+ B4 n& C: Z: k

∂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

[& \7 t6 o. B. |6 g0 `, J) D 7 R8 _1 s# l# z$ i6 P5 a. f0 g' k9 r% |2 | " [7 _ B* _' _6 K, j$ {* {0 E N3 O% w8 Q" m6 i' w
回复

举报 使用道具

相关帖子

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
快乐幸福
活跃在前天 15:23
快速回复 返回顶部 返回列表