沙运动规律数学模型究及清淤方案优选.docVIP

1 基本理论和方法 1.1 模型的基本方程 1.1.1 水流模型方程 水流连续方程： (1)水流运动方程：(2)(3)其中：x、y为空间坐标，t为时间,h为水深，z为水位z=z0+h,z0为床面高程，u,v为流速在x、y方向上的分量，f为柯氏力系数f=2ωsinφ,ω为地球自转速度，φ为纬度，Cs为谢才系数，曼宁公式Cs=h1/6/n,n为糙率，g为重力加速度，λ为涡动粘滞系数，，拉普拉斯算子。1.1.2 泥沙模型方程 ??? 悬沙连续方程： (4)??? 河床冲淤演变过程：(5)其中：C为悬移质含沙量，Kx,Ky分别为x、y方向的泥沙扩散系数，P为水流作用下泥沙的起悬量，D为水流作用下泥沙的沉降量，Cbx、Cby分别为x、y方向的河口浮泥层推移输沙率，z0为床面高程，Cm为床面极限含沙量。 1.2 泥沙模型方程中物理量的确定 1.2.1 泥沙沉降量D的确定 ??? 泥沙沉降量D与悬移含量C和泥沙沉速ω*有关，可表示为： D=Φ1Cω*(6)其中：Φ1为校正值。 1.2.2 泥沙起悬量P的确定 ??? 泥沙起悬量P与饱和挟沙状态下的含沙量C*和泥沙沉速ω*有关，可表示为： P=Φ2C*ω*(7)??? C*通过Lane—Kalinske公式确定为： C*=5.55I0[0.5(u*T/ω)exp(-(ω*/u*T)2)]1.61(8)式中u*T为潮流与波浪合成摩阻流速， u2*T=(u2+v2)/C2s+fwu20/4(9)u0=πH/[Tsh(2πh/L)]，波浪底部最大轨迹流速H—波高，T—波周期，L—波长，fw—波浪阻力系数。 ??? ab=H/2[sh(2πh/L)]-1，近底波浪质点振幅，r*=120d50床面粗度，d50床面泥沙中径；??? ,泥沙的平均沉速，ω1，ω2分别为泥沙最大(d1)与最小(d2)颗粒的沉速；I0为床面活动泥沙含量。??? 这样，将床面层附近泥沙起悬量和沉降量区分开来处理的方法，具有明确的物理意义，可以清楚地反映出悬移泥沙与床面层内运动泥沙及床面活动泥沙发生相互交换过程中矛盾的两个方面。在波浪掀沙和潮流输沙过程中，平衡只是相对的、暂时的，而不平衡则是绝对的、长期的。泥沙起悬量作为一个起悬函数(Pickup Function)的形式出现，有助于反映泥沙在床面层附近的交换机理，同时采用Lane-Kalinske公式，其优越之处在该公式在建立时，充分考虑了底部水流湍流脉动特性，反映了一定的力学机理，从目前流体力学研究来看，所不完善的是尚未区分大尺度拟序结构和小尺度随机脉动的作用，这方面我们正进行深入研究，但作为工程应用，该公式已经过了许多实验室和现场天然实测资料的验证，证明是可靠的[1,2]。总之，这样选取，有助于逐步弄清泥沙在床面层附近的交换机理，也有助于从力学机理上反映掀沙和输沙过程的不平衡特性。 1.3 模型方程的数值求解方法 1.3.1 水流方程的数值解 ??? (1)杂交有限元方法 ??? 杂交有限元方法是三角形单元和等参四边形单元相交错结合的有限元方法，这种方法对于具有复杂边界和长宽尺寸悬殊的水域是非常合适的[3,4]。 ??? 设单元中的近似函数取： (10)其中Φi(x,y)为基函数，Z、U、V取相同的基函数；如果单元为三角形，Φi(x,y)为线性三角元基函数；如果单元是等参四边形，则Φi(x,y)为四结点等参基函数。 ??? 将单元插值近似函数(10)代入方程(1)、(2)、(3)的Galerkin积分式中，可得杂交有限元方程： (11)(12)(13)??? 式中A、B1、C1、D1、G、B2、C2、D2、B4、C4分别为系数矩阵，其元素分别为：Aij=(ΦiΦj),，，，，，，，，，B3的元素与B1相同；C3的元素与C2相同。()代表在三角形单元或四边形单元内积分。对于不同的区域——三角形单元或等参四边形单元；i,j,k分别为1，2，3或1，2，3，4. ??? (2)质量集中方法 ??? 在过行长系列时间过程量模拟时，为节省计算时间和计算机内存，提高计算效率，按以下三种单元组合情况(图1)，分别采用集中质量法对单元系数矩阵进行化简，具体方法可参见文献[3]。 (a)三角形单元(b)等参四边形单元(c)三角形-等参四边形单元图1??? (3)时间方向的差分离散和半隐格式 ??? 为计算稳定，采用半隐差分格式，先利用Un、Vn和Zn通过离散式(12)和(13)，计算出Un+1和Vn+1值： AUn+1=(A-ΔtB1-ΔtG)Un-(C1·Δt-ΔtfA)Vn-ΔtD1Zn(14)AVn+1=(-ΔtfA-ΔtB2)Un-(C2·Δt-A+