弯曲槽道边壁振动情况下湍流特性的大涡模拟

为了更好地理解透平机械流道内能量耗散的机理和湍流结构变化,利用大涡数值模拟方法,数值试验了时、空做正弦振动的壁面所围成的弯曲槽道内的湍流流动。研究了6种不同的振动参数下,槽道内压力分布、近壁涡量以及其它相关湍流统计量的分布特性。计算结果表明,受展向壁面振动的影响,近壁区速度梯度增加,粘性底层变薄,湍流强度增加,能量耗散加大。该研究还发现壁面展向振动而引起的能量耗散主要依赖于近壁区粘性底层的展向涡量大小。     

应用基于RNG方法的湍流模型数值模拟钝体绕流的湍流流动

首次采用基于RNG方法的湍流模型对绕流钝体后尾流流场进行了数值模拟。计算在任意非正交曲线坐标系下,采用非交错网格的有限体积法求解二维不可压N-S方程,计算结果与实验数据以及与采用标准的κ-ε两方程湍流模型计算的结果进行比较,结果表明基于RNG方法的湍流模型对于与时间相关的大尺度运动——旋涡脱落的尾流的详细结构,能够给出真实的模拟。     

四种湍流模型对空化流动模拟的比较

采用四种湍流模型对NACA66水翼空化流动进行模拟,分析了不同湍流模型和壁面函数对空化流动模拟结果的影响.分析结果表明:ASM模型和RNG k-ε模型对空化数不敏感,且计算精度较高;Realizable k-ε模型对空化数的变化很敏感,对不同的空化数计算精度相差很大;标准k-ε模型对空化数不敏感,但精度不高.增强壁面函数法对空化数不敏感,计算精度较高;标准壁面函数法和非平衡壁面函数法对空化数有一定的敏感性,但计算结果差异不大.     

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

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

采用RNG紊流模型计算静止环境中圆形负浮力射流

基于RNG方法的κ-ε湍流模型,运用混合有限分析算法对静止环境中圆形负浮力射流进行数值模拟。从射流的最大入侵高度、流速矢量图、温度等值线、湍动能等值线、横断面上的流速和温度分布以及轴线上温度变化等方面与有效的试验资料进行了较为全面的对比验证。两者的良好吻合表明,基于RNG方法的紊流模型和混合有限分析法的结合能较好地模拟变密度流动。     

含闸墩溢流坝三维过坝水流数值模拟

利用有限体积方法对含闸墩溢流坝过坝水流三维流场进行数值模拟,采用κ-ε两方程模型模拟湍流,利用流体体积(VOF)法确定自由水面线,分别对两种墩型过坝水流进行计算,给出了3种水头下的水面线及坝面压力,通过与实验结果的比较表明,建立的模型能够准确地模拟带闸墩溢流坝三维流场,从而为溢流坝的水力设计提供可靠的设计依据。     

油气藏水力压裂计算模拟技术研究现状与展望

水力压裂技术是油气藏尤其是页岩气开发中的核心技术,利用数值模拟方法进行压裂优化和产能预测又是水力压裂成功的关键。本文首先介绍了水力压裂技术的发展历程。然后从计算模型(二维模型、拟三维模型和全三维模型)和数值模拟方法(基于连续介质和基于非连续介质)两方面对油气藏开发领域的水力压裂计算模拟技术进行较全面的总结。最后,从以下3个方面指出现今研究的不足并提出了进一步的研究建议:(1)全三维模型的完善-全三维模型应当与真实的工程参数和监测数据结合,用于校正模型本身,而校正后的全三维模型又可预测和优化新的现场水力压裂作业; (2)数值模拟方法的选用-已有的水力压裂数值模拟方法种类繁多,需要针对各种方法的适用范围、计算效率和模拟效果等,进行全面的比较和优化; (3)页岩储层中天然裂缝网络的数值模拟-天然裂隙网络加剧了页岩储层力学性质的各向异性,同时水力裂缝沟通天然裂缝活化扩展是有利于储层的增渗增产,对压裂缝网的形态、尺寸和连通率等起着至关重要的作用。因此,数值计算过程中综合考虑页岩储层中天然裂缝与水力裂缝的相互作用,将是未来水力压裂模拟的热点。     

裂隙水流运动的入渗试验及数值模拟研究

裂隙对喀斯特水流起着重要控制作用。选取具有垂向与水平向交叉裂隙的岩石剖面,做单环注水入渗试验,基于Navier-Stokes(N-S)方程组建立入渗水流数值模型,运用试错法,以实测单环水位变化为目标推求裂隙水力隙宽,计算入渗水量在不同裂隙中的比例;利用立方定律,估计出裂隙饱和渗透系数,基于地下水流连续性方程构建裂隙水流数值模型;比较两种数值模型计算结果。模拟结果表明:N-S方法推求的裂隙等效水力宽度远小于实测裂隙宽度;裂隙水流大部分沿垂向裂隙下渗,但也有部分水流沿横向裂隙渗流;立方定律与N-S方法数值模拟结果差异不大,前者估计的裂隙下渗率稍小于后者,相对而言,N-S方法能更精确地描述裂隙水流运动过程。     

基于Godunov格式的溃坝水流数学模型

为了更好地把握溃坝洪水风险,减小因溃坝洪水而造成的人员生命和财产损失,建立了基于Godunov格式的一维、二维溃坝水流耦合数学模型。一维溃坝水流模型采用HLL格式的有限体积法求解,二维溃坝水流模型采用基于非结构网格的Roe格式离散求解,在一维、二维模型的链接处采用重叠计算区域的方法实现一维模型和二维模型之间的水力要素信息交换。经弯道溃坝算例和断面突变溃坝算例验证,该耦合模型具有良好的可靠性和适用性,验证后的耦合模型为大尺度的溃坝水流数值模拟打下了基础。     

某应急水源地地下水资源评价

采用地下水模拟系统软件GMS6.0,建立评价区地下水数值模拟模型。通过地下水流模拟,直观地反映并预测出开采条件下地下水降落漏斗的形态。     

静止均匀环境中气泡射流特性的研究——形成区及形成后区特性的数值研究(Ⅱ)

采用气-液两相流两方程湍流模型,结合混合有限分析法,对静止环境中气泡射流进行数值模拟和分析。相对各种积分模型,这种模拟气泡射流的微分方程方法具有较强的预报性和普遍的适应性。数值计算结果与试验资料的较好吻合,从而验证了湍流模型和计算方法的有效性和正确性。     

基于结构网格的溃坝水流数值模拟

针对溃坝水流数值模拟面临的复杂地形和不规则边界等问题,基于结构网格建立了适应复杂地形和不规则边界的溃坝水流数值模拟有限体积模型（HydroM2D）。模型基于具有守恒特性的二维浅水方程,利用HLLC格式的近似Riemann解计算网格界面通量,利用MUSCL-Hancock法不断向前积分,使模型在时空上具有二阶精度;对源项进行离散处理确保模型的稳定性;模型引入有效干湿边界和不规则地形边界处理方法,准确模拟了干湿单元的动态交替和复杂边界上的水流特性。最后分别利用水槽试验、物理模型和实际算例对模型进行验证。结果表明,该模型对不同情景下的溃坝洪水模拟结果和实测资料以及现有模型模拟结果具有较高的一致性,模拟精度较高,稳定性较好,具有推广应用价值。     

静止均匀环境中气泡射流特性的研究--形成区及形成后区特性的数值研究(Ⅱ)

采用气-液两相流两方程湍流模型,结合混合有限分析法,对静止环境中气泡射流进行数值模拟和分析.相对各种积分模型,这种模拟气泡射流的微分方程方法具有较强的预报性和普遍的适应性.数值计算结果与试验资料的较好吻合,从而验证了湍流模型和计算方法的有效性和正确性.     

高拱坝宽尾墩三维流场数值模拟

高拱坝表孔宽尾墩收缩射流水舌的扩散特性与其出口流速分布密切相关,为了求得收缩射流的水力特性,应用RNGk ε双方程湍流模型对高拱坝表孔不同体型及过流量条件下可能出现的降水曲线流态、急流冲击波流态、水跃壅水流态、缓流流态4种流态进行了数值模拟计算,水流自由表面的模拟采用VOF法,得到了溢流水面及流速分布。结果表明,无论在何种流态下,计算所得的水面曲线均与物理模型试验值吻合良好,所得的出口流速分布结构也验证了高拱坝宽尾墩收缩射流水舌不同流态下分散的机理。     

非结构网格上的三维浅水流动数值模型

针对当前复杂环境水流模拟的需求,建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动力模型。模型采用有限体积法离散sigma坐标下的三维浅水方程,运用Roe黎曼近似解评估水平界面通量。模型网格拟合边界能力强,可根据需要局部加密;格式数值性能优良,具有守恒性、单调迎风性、高数值分辨率等特性。同时,应用干湿判别法处理动边界,以适应浅滩地形漫/露过程模拟的需要。封闭水池内部风生环流、干河床上溃坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点,结果表明它能够准确地预测水流的三维流动结构,而且计算简单高效,具有良好的数值稳定性。     

基于三维渗流模拟的坝基渗漏和渗透稳定分析——以哇沿水库为例

水库坝基渗漏不但损失水库蓄水量,更重要的会引起大坝渗透变形失去稳定。应用三维地下水流数值模拟计算坝基渗漏量和分析渗透稳定,较好地解决了周边与底面边界的不规则问题和垂向上地层的非均质问题,提高了计算精度。以青海省哇沿水库为例,在充分分析水库坝基水文地质条件的基础上,建立了坝基渗漏三维地下水流数值模型,模拟了设置不同深度防渗墙时的渗漏量和坝后出溢段水力坡度。评价了大坝的渗透稳定性,为防止哇沿水库发生渗透变形,设置防渗墙深度应大于50m。当防渗墙深度为50m时,水库渗漏量为12 012m~3/d。     

各向异性密度流模型及其验证

针对密度流流动各向异性的特性,采用非线性k-ε模式求取雷诺应力,结合RNG模型及浮力修正建立了三维密度流模型.应用该三维模型计算了水库温差异重流,分析了温差异重流演进过程,并就流速分布与水温下泄过程与实测资料进行了比较.结果表明,所建模型由于引入非线性RNG模式,能有效反映密度流的各向异性特性,较好地模拟了流场和温度场强烈耦合的水库温差异重流,能准确预测温度分层过程与异重流潜行时间,具有较高精度.     

水生植物对水沙运动影响的三维湍流模型

水生植物的存在改变了明渠内水流的流动结构,也影响着泥沙的输运。建立三维湍流模型,在水流控制方程中加入植被阻力项和植被密度项来考虑刚性植物对水动力特性和泥沙输运特性的影响。应用该三维数值模型计算了矩形水槽内淹没植被对水流水平时均流速垂向分布的影响、复式明渠边滩栽种挺水植被对水流深度、平均流速分布以及植被岛周围泥沙床面的冲淤变形的影响。数值计算结果与实测结果吻合良好,表明本模型可以有效地描述刚性水生植物对水流泥沙运动的影响。     

非饱和土水-力本构模型及其隐式积分算法

在已有工作基础上建立了水力-力学耦合的非饱和土本构模型,在硬化方程中考虑饱和度的影响,同时在土水特征曲线中考虑了塑性体变的影响,从而使模型可以反映非饱和土中的毛细现象与土中弹塑性变形现象的耦合行为。采用隐式积分方法,建立了非饱和土耦合模型的数值模型,并推导了得到了水力-力学耦合的非饱和土的一致切线模量。利用该算法编制了本构模型计算的子程序,使其能向外输出切线刚度矩阵,用于有限元计算。为了验证该算法和程序的正确性,用所编制程序对不同路径下的土体行为进行了预测。通过预测结果与试验结果相对比,表明程序预测结果与试验数据相吻合,模型可以较好地模拟土体的水力-力学耦合行为特性。     

基于Godunov格式的排水管网水流数值模拟

为克服基于"LINK-NODE"计算模式的排水管网模型不能细致模拟管道内部水流运动过程的缺陷,采用Preissmann窄缝理论开发完成了一套新的基于Godunov有限体积格式的一维排水管网水动力模型。模型采用HLL近似Riemann解计算单元界面处的数值通量,细致区分了管网节点水位对与其相连接管道的不同影响,并采用严格的质量守恒方程和特征线理论来更新管网节点水位。通过3个经典算例验证了模型能够很好地模拟有压瞬变流和瞬变混合流等复杂管网流态,具有良好的精度和稳定性。将模型应用于陕西省西咸新区天福和园小区实测场次降雨排水过程模拟,取得了理想的计算结果。新研发模型可提供所有计算单元中心和管网节点处的水力要素变化过程,为精细化城市洪涝过程模拟奠定基础。