悬浮物扩散和地形冲淤变化的数值模拟方法.docxVIP

悬浮物扩散和地形冲淤变化的数值模拟方法

悬浮物扩散和地形冲淤变化数学模型建议采用较成熟的二维模式。

D.1基本方程组

D.1.1连续方程：

（D.1）

D.1.2动量方程：

（D.2）

D.1.3悬沙扩散方程：

（D.3）

D.1.4床面冲淤变化方程：

（D.4）

式（D.1）、（D.2）、（D.3）、（D.4）中：

—总水深(m)，；

-水位(m)；

—水深(m)；

—代表x方向的流速分量（m/s）；

—代表y方向的流速分量（m/s）；

f—柯氏力系数（s-1）；

—风摩擦系数；

—代表海平面上10m处，x方向的风速分量（m/s）；

—代表海平面上10m处，y方向的风速分量（m/s）；

—谢才系数，；

—曼宁系数；

—代表x方向的海水紊动粘滞系数（m2/s）；

—代表y方向的海水紊动粘滞系数（m2/s）；

ω--沉速（m/s）；

—含沙量（kg/m3）；

—挟沙力（kg/m3）；

—代表x方向的泥沙紊动扩散系数（m2/s）；

—代表y方向的泥沙紊动扩散系数（m2/s）；

—泥沙干容重（kg/m3）；

—冲淤厚度（m）。

D.2基本参数确定

D.2.1水流紊动粘滞系数Kx、Ky

由试验确定，或取经验值50～500m2/s。

D.2.2紊动扩散系数Dx、Dy

由试验确定，或采用Elder公式

（D.5）

（D.6）

D.2.3挟沙力S*

由于经验公式较多，这里推荐以下几个常用公式：

（1）张瑞瑾公式

（D.7）

式中：

—挟沙力（kg/m3）；

-系数；

—流速（m/s）；

-水力半径；

--沉速（m/s）；

-经验系数；

g-重力加速度（m/s2）

（2）刘家驹公式

（D.8）

式中：

—挟沙力（kg/m3）；

—系数；

—泥沙干容重（kg/m3）；

—流速（m/s）；

-海平面上10处风速（m/s）；

—水深(m)；

g-重力加速度（m/s2）

（3）窦国仁公式

（D.9）

式中：

—挟沙力（kg/m3）；

-代表泥沙密度（kg/m3）；

-代表海水的密度（kg/m3）；

-曼宁系数；

—流速（m/s）；

-海平面上10处风速（m/s）；

—波高（m）；

—水深(m)；

-沉速（m/s）；

-水温（℃）

（4）曹祖德公式

（D.10）

式中：

—挟沙力（kg/m3）；

—系数；

—流速（m/s）；

—波速（m/s）；

-粒径（mm）；

--沉速（m/s）

D.2.4时间步长确定

（D.11）

式中：-系数，可取1～3；

-网格最小边长；

-计算域最大水深

D.3初始条件和边界条件

D.3.1初始条件

（D.12）

（D.13）

（D.14）

（D.15）

D.3.2固边界条件

流速：（D.16）

泥沙通量：（D.17）

D.3.3开边界条件

潮位：（D.18）

或流速：（D.19）

悬沙流入：