平面二维溃堤水流泥沙数值模拟

基于一阶迎风格式平面二维河道溃口水流模拟,建立了非均匀沙模型并用于模拟平面二维溃堤水流泥沙及河床变形。给出了实际算例,证明该模型计算速度快,精度高,符合水流泥沙运动基本规律。     

明渠交汇水流阻力对平面二维数值模拟的影响

与顺直明渠水流相比,明渠交汇水流由于存在分离区和自由剪切面,其紊动特性引起的阻力较为复杂,交汇水流除床面阻力外还受到较强的阻力。基于数值模拟与实验比较分析表明,分离流与自由剪切流所产生的阻力(紊动粘性阻力)对于交汇水流数值模拟具有重要影响。同时,由于平面二维模型对二次流影响的忽略,断面环流较强时数值模拟会产生较大的阻力(环流阻力),而断面环流较弱时,交汇水流的三维特性相对较弱,采用平面二维模型和k-ε紊流模型能达到较好的模拟效果。因此,紊动及环流阻力的计算对交汇水流数值模拟的精度至关重要。     

山区河流平面二维流场的数值模拟

针对山区河流边界地形复杂、水位及流速变化剧烈的特点,用"混合五对角法"建立平面二维数学模型并求解,进行水位和流速的数值模拟.利用嘉陵江重庆金刚碑-朝阳桥河段地形及5个横断面的水位和水深平均流速资料对此模型进行了验证,模型快速、准确地再现了该河段的水位和水深平均流速情况.     

淹没丁坝群二维水流数值模拟新方法

为模拟淹没丁坝群平面二维水流运动,提出了淹没丁坝群二维水流数值模拟新方法并建立了数学模型。新方法的主要实施方案:① 将丁坝视为无厚度坝,用网格线概化丁坝;② 采用新的网格节点布置形式,即水深、流速节点布置于网格界面上,水位节点布置于网格中心,有别于一般交错网格节点布置。模型采用基于结构网格下的有限体积法对方程组进行离散,同时将淹没丁坝坝顶水深代入离散方程中进行求解。采用已有的水槽试验资料,进行了初步验证,模拟了长江下游东流水道已建丁坝群工程实施后河道的流场和水位场,结果表明计算和实测符合较好。     

通用河网二维水流模拟模式研究

通过对流域河网二维水流模拟问题的分析,提出了一个通用化的河网二维模拟概化模式,该模式将河网二维概化为"树状"河道计算单元、"环状"河道计算单元、"十字型"河道计算单元等河网二维基本单元,对这些河道计算单元的求解,构建流域型河网一二维耦合模式的求解模型,并由此建立了通用化的河网二维水流数学模型软件系统。通过"树状"河道二维计算单元求解模型的推导与实现,验证了概化模式的可行性,表明该思路可以应用于大型流域河网一、二维水流模拟,具有很好的应用推广价值。     

浅水湖泊平面二维水流-水质-底泥污染模型研究

从三维对流-扩散方程出发,推导包含底泥污染的水质控制方程,并与浅水流动方程耦合,建立相应的浅水湖泊平面二维水流-水质-底泥污染的数学模型。运用伽辽金加权余量法推导出数学模型方程组的有限元公式,开发了实时二维水动力、水质模型。利用模型计算玄武湖混合流、水质动态变化过程,结果与实测值较吻合,从而为浅水湖泊混合流、包含底泥污染的水质模拟提供了一种可操作的新途径。     

多沙河流平面二维泥沙数学模型研究

针对冲积河流河床演变速度快、边界调整幅度大的特点,采用积分方法数值求解水流泥沙基本方程,建立了适合于多沙河流泥沙及河床演变规律的二维数学模型,在数值方法方面考虑了物质守恒性及对多变边界的适应性,采用了比较简单的动边界处理模式。通过工程实例的计算验证,证实了模型所采用的基本方法是合理的,模型有较强的河床边界适应能力,可用以解决有关多沙河流河床演变预测及浑水洪水预报的有关问题。     

平面二维泥石流数值模拟方法的探讨

本文提出了泥石流二维平面流动方程组，基于研究泥石流堆积泛滥的流速与泥深分布图像的构思，笔者应用有限差法进行了数值解，其计算结果与实测结果基本吻合。     

水流紊动对泥沙絮凝的影响

主要研究水流紊动对泥沙絮凝的影响。认为水处理方面的絮凝动力学理论可用于研究泥沙絮凝问题。根据试验实测资料及国外有关试验结果研究表明,在紊动状态下,泥沙絮凝存在两个阶段,每个阶段中紊动对泥沙絮凝的影响也有所不同。     

冰盖下水流速度分布的二维数值模拟分析

基于标准的κ-ε模型,利用有限体积法,对冰盖下水流速度分布进行了二维数值模拟和分析.结果表明:当冰盖和河床糙率比值变化时,水流速度剖面发生变化;冰盖下的最大流速值比速度对称分布时要大,且偏于光滑的一边;在相对光滑壁处近壁流速值大于相应对称分布时的流速值,该结果与实验室实测速度剖面情况吻合.数值试验表明,最大流速点位置的偏移量仅在有限的范围内,冰盖区、床面区及冰盖下水流平均速度相差较小,表明了工程应用中假定各区流速相等是可行的.     

无结构网格上平面二维水沙模拟的有限体积法

基于无结构网格有限体积法的算法框架,通过引入跨单元界面法向水沙数值通量的逆风分解,将悬沙与床沙交换以及分组挟沙力计算模式自然地嵌入二维水沙运动方程组的数值格式中,形成高精度、守恒性好的二维水沙有限体积算法。最后,利用该算法对谭江樟州河段的水沙输运和河床变形进行了数值模拟。结果表明,该算法能够较好地模拟复杂条件下河道水沙输运的往复特征和河床变形的动态过程,其精度满足河道工程后效分析的要求。     

水沙运动一、二维分区耦合虚拟点算法

建立了一种计算具有多空间尺度水沙运动的一、二维分区耦合虚拟点算法。该算法通过在一、二维区域的交界面定义虚拟点,并利用有限元的非结构化网格与单元叠加的特点,将所有一、二维区域连接点的影响系数集成到总体阵中虚拟点的元素上,实现不同分区之间的强耦合作用。同传统的使用简单匹配条件的一、二维分区弱耦合算法相比,该方法可以提高分区耦合计算的稳定性。对算法的正确性用一河口算例进行了检验。     

淤泥质河口水沙运动数学模型相关问题

为研究淤泥质河口的水沙运动规律,建立了用于模拟淤泥质河口水沙运动的二维数学模型。该模型采用基于无结构三角网格下的有限体积法对方程组进行离散,结合Roe-MUSCL方法及时间方向的预测-校正格式,使模型在时空方向具有二阶计算精度。模型中分别采用不同方法计算粘性和非粘性泥沙的输移源项,并引入粘性泥沙的起动流速和冲刷率计算公式。采用已有的概化水槽试验数据对模型进行了初步验证。然后模拟了1995年10月小潮及大潮期间海河口的潮流运动与泥沙输移过程,计算得到的潮位、潮流速及含沙量过程与实测过程符合较好,结果表明模型能够用来模拟淤泥质河口粘性和非粘性泥沙的不平衡输移过程。同时还比较了泥沙输移源项的不同处理方式对计算结果的影响,计算表明在淤泥质河口水沙运动数学模型中必须同时考虑粘性和非粘性泥沙的输移。     

二维明渠非恒定流的格子Boltzmann模拟

根据二维浅水波方程的特点,应用由格子Boltzmann方程展开的多尺度方程,建立了一个能模拟二维明渠非恒定流的格子Boltzmann模型。通过对一维溃坝波和二维圆形溃坝波的模拟分析,验证了该模型的有效性。以某水电站尾水渠为例,模拟分析了由机组甩负荷和增负荷引起的水位波动规律,为尾水渠的设计提供了依据。     

三门峡库区一维非恒定非均匀泥沙输移数学模型

黄河是高含沙河流,含沙量季节差异显著。因此研究黄河干流的非恒定水沙输移规律以及库区泥沙的淤积问题极为重要。应用圣维南方程组以及非恒定泥沙连续方程建立了非恒定非均匀泥沙含沙量计算公式,并根据沙量平衡方程推求出三门峡库区河底高程的变化规律。使用已有的资料对模型进行的验证表明:模型计算与实测资料符合良好,该模型具有较好的应用前景。     

方腔回流区水流运动特性三维数值分析

建立了一个三维水动力学模型,对方腔回流区水流运动特性进行了模拟分析。数学模型以三维浅水环流方程为基础,在垂向空间引入σ坐标变换,采用高分辨率半隐有限元离散格式。计算结果与解析解及水槽试验数据吻合良好。应用该模型对方腔回流区的水流运动特性进行了计算分析,进一步揭示了方腔回流运动的非定常非对称不封闭特性、流动结构表现为竖轴环流与立面环流相叠加、流速沿垂线分布相对均匀等流动规律。     

泥石流的二维数学模型

泥石流是在重力作用下,由砂粒石块和水等组成的固液混合物,是一种发生于山区的复杂的地质灾害现象。泥石流主要是由暴雨诱发引起的,它沿着复杂的三维地形高速流动,具有流体流动的特性。为了模拟泥石流的运动规律,预测降雨诱发的泥石流的到达距离和泛滥范围,减少和避免泥石流引起的灾害,把泥石和雨水组成的固液混合物假定为遵循均匀、连续、不可压缩的、非定常的牛顿流体运动规律。基于质量守恒方程和Naiver-stokes方程,采用深度积分方法,推导出了一个模拟泥石流运动的二维数学模型。所有方程式可用有限差分法来求解。结合GIS,该模型可用于预测泥石流的流动距离和泛滥范围,以及泛滥范围内的危险房屋和路段,也可以用于泥石流灾害的风险性分析。