二维二维非恒流数值模拟在河道冲淤变化中的应用.docxVIP

二维二维非恒流数值模拟在河道冲淤变化中的应用 防洪影响评估是河流非防洪工程建设中的一个重要技术问题。防洪评价计算可选用规范推荐的经验公式、数学模型和物理模型来进行计算模拟。对于简单的由于建筑物阻水影响洪水下泄的情况,一般可采用经验公式计算壅水,而壅水高度和范围对河段防洪影响较大的应采用数学模型或物理模型试验。由于物理模型投资大、周期长等原因,使其适用性受到一定限制;数学模型可细致模拟河道内工程位置附近水流泥沙的运动规律,且具有投资少、周期短等优势,可同时进行多方案分析计算,具备极好的灵活性,目前在此类问题研究中应用较为广泛。 二维浅水方程数值模拟在防洪影响评价中应用较多。从Leendertse开始,各种不同的算法被运用到深度平均的浅水方程数值模拟中。新的算法研究包括多方面的内容,如:空间差分方法,时间差分算法,方程求解技术和边界条件处理技术。开发新的数值算法必需在算法复杂性和计算效率两者直接找到平衡。开发的算法必需满足连续性、稳定性和收敛性。算法必需满足具有物理现实意义和总量守恒。 有限体积法正是基于守恒方程积分形式的一种算法。应用有限体积法求解N-S方程和浅水方程已经被广泛应用,对有限体积法的优势也做了充分的说明。有限体积法求解的关键问题是通过各单元边界法相数值通量的估算。常用的估算方法有:Osher格式,通量向量分裂格式(FVS),通量差分裂格式(FDS)。 本文采用有限体积法,运用FDS算法,建立了平面二维水流数学模型,对浙江省49省道温溪至鹤城段改建工程建设前后的河道水流流速、水位、流态、河床的变化情况进行准确模拟和计算,并由此分析工程建设对防洪的影响,为该工程提供技术支持。 1 数学模型的建立 1.1 控制方程 对于平面大范围的自由表面流动、水深尺度远小于平面尺度、无明显垂直环流、垂向流速小的浅水流动,可以引入静水压力假设,并沿水深方向进行积分来简化方程,简化后的方程即为二维浅水方程。其守恒形式为: ?q?t+?f(q)?x+?g(q)?y=b(q)(1)?q?t+?f(q)?x+?g(q)?y=b(q)(1) 式中:q=[h,hu,hv]T为守恒物理向量;f(q)=[hu,hu2+gh2/2,huv]T为x向的通量向量;g(q)=[hv,huv,hv2+gh2/2]T为y向的通量向量;h为水深;u和v分别为x向 和y向的垂线平均匀流速分量;b(q)=[0,gh(s0x-sfx),gh(s0y-sfy)],s0x和sfx分别为x向的河底坡度及摩阻坡度,s0y和sfy分别为y向的河底坡度及摩阻坡度。 1.2 积分离散 应用散度定理对方程(1)在任意单元Ω上进行积分离散,求得FVM的基本方程: ?Ω?q?tdΩ=-∫?ΩF(q)?ndl+?Ωb(q)dΩ(2)?Ω?q?tdΩ=?∫?ΩF(q)?ndl+?Ωb(q)dΩ(2) 式中:n为?Ω单元边外法向单位向量;dΩ和dL为面积分和线积分微元;F(q)·n为法向数值通量,F(q)=[f(q),g(q)]T。 假定物理量q和函数b(q)在单元Ω内均匀分布,线积分用沿单元边界的通量来近似,则式(2)可以写成如下的离散形式, Adqdt=-m∑j=1Fjn(q)Lj+Ab(q)(3)Adqdt=?∑j=1mFjn(q)Lj+Ab(q)(3) 式中:A为单元Ω的面积;Fjnjn(q)为通过单元边法相通量;Lj为单元边j的长度;m为单元Ω的边的数目。 单元边法向通量Fjnjn(q) 简记为Fn(q),其定义如下: Fn(q)=cosΦ?f(q)+sinΦ?g(q)(4)Fn(q)=cosΦ?f(q)+sinΦ?g(q)(4) 式中:Φ为法向单位矢量n与x轴的夹角。 根据通量向量f(q)与g(q)的旋转不变性,即满足下面的关系: Τ(Φ)Fn(q)=f[Τ(Φ)q]=f(ˉq)(5)T(Φ)Fn(q)=f[T(Φ)q]=f(qˉ)(5) 即 Fn(q)=Τ(Φ)-1f(ˉq)(6)Fn(q)=T(Φ)?1f(qˉ)(6) T(Φ) 和T(Φ)-1分别为坐标旋转变换矩阵及其逆阵,表达式为: Τ(Φ)=[10000cosΦsinΦ00-sinΦcosΦ00001]Τ(Φ)-1=[10000cosΦ-sinΦ00sinΦcosΦ00001](7)T(Φ)=??????10000cosΦ?sinΦ00sinΦcosΦ00001??????T(Φ)?1=??????10000cosΦsinΦ00?sinΦcosΦ00001??????(7) 将式(6)代入方程(3),可得: Adqdt=-m∑j=1Τ(Φ)-1f(ˉq)Lj+Ab(q)(8)Adqdt=?∑j=1mT(Φ)?1f(qˉ)Lj+Ab(q)(8) 式中:ˉqqˉ由向量q变换而来,其相应的流速分量分别为法向和切向。 由公式(8)可