黄河口泥沙异重流基本控制参数的数值试验

《黄河口泥沙异重流基本控制参数的数值试验》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、黄河口泥沙异重流基本控制参数的数值试验摘要：为了定量地评价泥沙异重流基本控制参数对泥沙异重流流动的影响，本文设计了一系列数值试验来展示各种参数与泥沙异重流流动特征之间的定量关系。采用ADI法求解三维方程经垂向积分后的平面二维泥沙异重流方程。所选取的基本控制参数包括水下斜坡坡角，总阻力系数，水力挟带系数，泥沙饱和浓度和泥沙沉速。通过数值试验得出水下斜坡坡角θ、泥沙饱和浓度Es越大，总阻力系数cd、泥沙沉速ω越小，相应地泥沙异重流流速、流动厚度以及悬沙浓度越大。其中斜坡坡角对泥沙异重流影响最大。当泥沙异重流流动处于超临界流状态，水力挟带系数Ealstyle="TEXT-A

2、LIGN:left"align=left>　　泥沙异重流及形成的地质体在沉积学和石油地质中有重要意义。在黄河口，泥沙异重流是泥沙搬运的重要方式，是河口区沉积动力的重要过程，是认识河口三角洲发育演化的重要途径。为使研究深入一步，我们进行了泥沙异重流的数值模拟。近年来，国内对异重流数模研究较少，仅见方春明等(1997)关于泥沙异重流潜入时的立面二维数值模拟以及邱晨霞(1995)对盐水异重流的二维两层数值计算。1平面二维泥沙异重流控制方程及数值模拟　　目前，包含三个基本控制方程(流体质量、动量守恒方程和泥沙质量守恒方程)的泥沙异重流数值模拟被广泛应用，本文把一维方程推广到平

3、面二维的情形，成为四方程模型，如下所示[1](1)(2)(3)(4)其中u，v分别为x，y方向速度，θx，θy分别为x，y方向坡度，Ewx与Ewy分别为x、y方向的水力挟带系数，τx、τy为x方向与y方向的剪切摩擦力，kx678910111213、ky分别为x方向与y方向的悬沙扩散系数，t为时间，h为异重流厚度，s为深度平均的悬沙含量，ρs，ρalstyle="TEXT-ALIGN:left"align=left>在方程(1)中，水力挟带系数Ealstyle="TEXT-ALIGN:left"align=left>　　Ashidaetal.(1975):Ealstyl

4、e="TEXT-ALIGN:left"align=left>　　Parkeretal.(1987):Ealstyle="TEXT-ALIGN:left"align=left>　　Fukushimetal.(1985):Ealstyle="TEXT-ALIGN:left"align=left>　　以及俞维升(1991)：Ealstyle="TEXT-ALIGN:left"align=left>在方程(2)、(3)中，剪切摩擦力τx、τy是作用在异重流上下界面摩擦力的总和，它们由下式所示。其中cd=cf(1+α′),被称为总阻力系数，范家骅等(1980)的cd值约0.00

5、3，Parkeretal.(1987)的cd值在0.001～0.058之间，俞维升(1991)的cd值在0.003～0.004之间，可见cd值有很大的不确定性。　　在方程(4)中，泥沙饱和浓度Es是描述泥沙异重流流动特征参数的函数，对黄河口可用下式表示　中国水科院(1997)公式[2，3]s≤15kg/m3,Es=123(V2/gh)0.36(V/ω)-0.33(1-Δh/h)0.2s>15kg/m3,Es=9.7(V2/gh)0.01(V/ω)0.16(1-Δh/h)0.22V为合流速,h为水深,ω为泥沙沉速，Δh为潮差。　　张青玉经验式Es=9.83(V2

6、/h)0.23定解条件：　　边界条件：陆地边界：Vn＝0(n为岸线的法线方向),s′=0(′表示陆地)　　水边界：h*(t)=h(t),s*=s(*表示水边界)；河口边界：u=1.5(m/s)s=50(kg/m3)　　初始条件：u=0v=0s=0h=0.5(m)　　采用有限差分法中的ADI法求解方程组(1)、(2)、(3)、(4)。事实上，我们应当同时求解描述泥沙异重流流动的方程，因为它们是一组相互耦合的方程组。然而同时求解相当复杂，我们便采取一种叫冻结系数法[4]的近似耦合法来求解。　　计算黄河口泥沙异重流时，泥沙异重流起点在黄河入海处。数值计算区域如图1所示，x、

7、y方向空间步长均取1/3千米，时间步长取67891011121324秒；计算区海底地形如图2所示；柯氏参数f＝2Ωsinφ，Ω取7.2722×10-5s-1，φ取36.7°；泥沙密度ρs取2.65g/cm3，水体密度ρ3，泥沙异重流密度ρf取1.040g/cm3。图1泥沙异重流数值计算区域Areaofnumericalsimulationofsedimentdensitycurrent图2计算区海底地形图Mapofareaofnumericalsimulation2平面二维泥沙异重流的参数试验　　为了定量地评价基本控制参数对泥沙异重流流动的影响，本文