非正交曲线同位网格三维水动力数值模型

建立一个自由表面三维水动力学计算模型，平面采用贴体坐标变换，竖直方向σ坐标变换，基于有限体积法，同位网格布置即流速，压力布置于同一控制点上，自由表面由二维水深积分方程中的水位方程求解，运用SIMPLEC法求解速度－压力耦合问题，动水压力校正方程校正压力和流速场，标准k-ε紊流模计算紊动粘性系统以封闭Reynolds应力项，运用该数值模糊模拟S形弯道水流三维流动。     

等宽明渠交汇水流数值计算

针对明渠交汇口缓流流动特性,分别应用水深平均H-L模型及k-ε模型建立了二维数值模型．水深方向采用静压假定,可模拟自由面和河床地形变化;应用交错网格上的有限体积法对方程组进行离散并采用水深-速度校正算法求解;每个代数方程求解采用TDMA逐行扫描法;针对交汇口几何特征,分区进行网格划分、应用实测资料对两种模型进行了验证比较,结果表明两种模型均可模拟明渠交汇口基本流动特性,但H-L模型所得的分离区计算结果比k-ε模型更为准确．     

用HH-SIMPLE方法求解明渠中三维浅水底孔泄流问题

本文针对明渠中三维底孔喷射流问题,提出了一个新的求解 N-S 方程的有限差分方法——HH-SIMPLE 方法.用该方法求解三维水流问题,可以保持 SIMPLE 方法的优点,即可将动量方程中的速度与压力解锅;求解离散化动量方程时可用隐式迭代格式,而无需作“刚盖近似”或“静水压力假设”,因而可很好地处理自由表面.本文中用该法求解的算例在国际计算流体力学领域中目前尚为空白,其三维计算结果可精确地描述在上游有来流情况下的浅水底孔喷流的流场、自由水面和压力场问题.计算结果和试验资料吻合较好.     

求解不可压N—S方程的完全压力校正方法

为解决在非交错网格上求解不可压N-S方程遇到的压力波动问题,将一种完全压力校正方法推广到三维非正交曲线坐标系中.该方法的校正压力由质量守恒项和压力平滑项组成.对带有移动顶盖的三维方腔层流和90°方截面弯管层流流动问题的数值模拟结果表明,该方法准确可靠,程序编制简易,计算效率明显提高.     

澳门周围水域潮汐流动的数值模拟

采用有限体积法求解二维浅水流动方程,建立了一种在非规则结构化网格上数值求解平面二维浅水流动的数值方法,其中应用Roe格式计算数值对流通量,并通过MUSCL途径插值界面状态,以对水位插值代替水深插值来消除地形变化对格式的影响;对源项采用分步法及半隐式求解,防止了源项引起的计算不稳定;处理动边界采用人工小水深的方法．应用于澳门水域潮流的模拟,计算结果与实测吻合较好．     

蓄滞洪区平面二维干河床洪水演进数值模拟

在地形较复杂的情况下,采用常用的处理干河床问题的方法有时会出现水流虚假流动的现象。本在洪水干河床演进模拟中提出了“水位平铺”假设进行处理,取得了比较好的模拟效果。模型建立过程中采用有限差分法的逆风格式离散控制方程,并用ADI方法求解代数方程组,模型计算速度快,节约计算内存,基本可以满足洪水实时预报要求。     

交汇水流三维数值模型

应用显式代数应力模型建立了交汇水流三维数值模型。方程离散采用交错网格上的有限体积法,并采用SIMPLE算法求解,代数方程求解采用TDMA逐行扫描法。针对交汇口几何特征分区进行网格划分,自由水面求解采用VOF方法。应用实测资料对计算结果进行了验证,同时与应用k-ε紊流模型的三维数值计算结果进行了比较。研究结果表明,文中所提出的数值模型可以更好地模拟交汇水流独特的三维流动特性。     

振动作用下单螺杆挤出机内聚合物熔体的流动

提出了一种求解二维耦合流场的有限差分数值模拟新方法，避免了求解压力Poisson方程。对振动力场作用下聚合物熔体在单螺杆挤出机螺槽内的流场进行了数值计算，得到了不同时刻螺槽内聚合物熔体的流动速度与剪切应力分布。结果表明，振动将导致螺槽内聚合物熔体的流动速度与剪切应力周期性变化，当振动扰动系数较大时，螺槽展开方向甚至会出现瞬时压力反流．最后，通过剪切应力的分布曲线证明本研究所得结果的正确性。     

三维弹塑性有限元法模拟薄带轧制过程(Ⅰ)——关于接触摩擦模型及其应用

为解决压力加工中的接触摩擦问题,根据金属流动反作用原理,本文提出了一个新的接触摩擦模型,经三维弹塑性有限元程序计算以及对凹辊面薄带轧制的金属三维变形分析表明,在不用对摩擦力大小、方向及中性点位置作任何假定的情况下,新模型能够有效地模拟中间物质剪切及压力传递。     

非结构网格下两步压力校正算法的潮流模拟

针对非结构网格下的潮流数值计算,采用有限体积法离散基本方程,对流-扩散项离散采用幂率格式,交叉扩散项引入超松弛校正方法,水位校正方程应用Rhie-Chow动量插值思想和SIMPLEC类算法导出.采用两步压力校正法将交叉扩散项单独考虑,以提高模型对高度畸变网格的适应能力.计算结果表明,该处理方法可提高压力欠松弛系数,并能有效改善复杂区域潮流模拟的健壮性.通过对黄浦江感潮河段潮流的模拟,表明该模型的预测结果与实测资料吻合良好,模型对潮流的模拟结果是可靠的.     

振荡流底边界层运动的数值研究

将大涡模拟法简化为二维形式,利用SGS格子涡模型封闭二维Navier-Stokes方程水流运动方程组,得到平底振荡流动边界层立面二维水流数值模拟,控制方程采用SMAC法求解,计算结果与水槽实验的实测资料比较较表明,该模型能够模拟振荡流边界流动特性,并进一步讨论了振荡流边界的紊动特性沿垂线分布和随相位变化的情况。     

取水池内三维流场的数值计算

本文应用不可压缩流体三维Navier—Stokes方程组求解明渠内的三维水流运动．利用k-ε双方程湍流模式计算涡粘性系数．控制方程组采用隐式差分离散，SIMPLEC程式求解．文中对矩形明渠开挖槽沟后挖槽区的流场进行了计算．所得结果与有关文献给出的测量值比较一致．作为应用，对某电厂循环供水系统中取水池内的三维流场进行了计算．计算结果可供工程界参考．     

拱坝坝面动水压力的加权余量解法

本文假定库水为无粘性的不可压缩流体,将坝-水系统作为三维问题.利用加权余量法给出了一种求解拱坝坝面动水压力的计算方法.文中研究了矩形河谷上拱坝坝面动水压力的分布规律,并讨论了河谷宽高比和拱坝中心角等因素对坝面动水压力的影响.算例结果表明,本文所建议的方案简便易行,计算效率高,是求解拱坝坝面动水压力的一种较好方法.     

基于非结构网格的三维浅水模型

采用水位积分平衡法并引入边壁滑移系数,解决σ坐标系下有限体积离散时产生的静水压力项与底坡项的平衡问题和边壁阻力的模拟问题,建立非结构网格上三维浅水模型.风生流试验计算结果表明,模拟的纵向流速的垂向分布与理论分布吻合较好;Holtz丁坝试验计算结果表明,引入边壁滑移系数能够模拟坝体附近上游小回流区和下游小回流区;淹没丁坝...     

内燃机缸内冷态湍流流动及传热的二维数值计算

阐述内燃机缸内冷态湍流流动和传热的一种二维数学模型，即用数值方法求解控制缸内气体流动和传热的一组偏微分方程组，得到缸内流场及温度场，进而求出壁面处温度的变化规律．计算是基于控制容积法，采用可在轴向伸缩变化的轴对称动坐标系．在压力较正方程和Ｋ－ε湍流模型中均考虑了流体的压缩性效应，旋流也作为因变量之一纳入了计算．对平顶及带凹坑的燃烧室的模拟计算结果与实验数据相吻合．     

应用滑移界面方法求解动／静相干 叶排内非定常流动

基于滑移界面技术，发展了求解二维动／静相干叶排内非定常粘性流动的计算方法和软件，该方法克服了混合平面方法所引入的方法误差．压力与速度耦合采用非迭代PISO算法；计算区域的离散采用非交错斜交网格．计算了二维不可压缩Couette流动及无限长相对运动动／静叶排内的非定常流动，验证了用该方法和软件计算动／静相干叶排内非定常流动的有效性．     

连铸结晶器中钢液三维紊流流动的数值模拟

针对连铸结晶器中钢液的紊流流动特点，建立了有限元求解这一限定空间射流的三维紊流时均场的方法，用水模拟试验结果验证了文中所用模型和计算方法的可靠性，且对实际工况条件下的板坯连铸结晶器中的三维紊流流动进行了模拟计算．其方法可用于分析各种结晶器液池中的流动特征，并为凝固壳厚度的模拟计算和板坯质量控制奠定了基础．此外，指出了前人计算中在边界条件上的不合理性     

应用滑移界面方法求解动/静相干叶排内非定常流动

基于滑移界面技术,发展了求解二维动／静相干叶排内非定常粘性流动的计算方法和软件,该方法克服了混合平面方法所引入的方法误差．压力与速度耦合采用非迭代ＰＩＳＯ算法;计算区域的离散采用非交错斜交网格．计算了二维不可压缩Ｃｏｕｅｔｅ流动及无限长相对运动动／静叶排内的非定常流动,验证了用该方法和软件计算动／静相干叶排内非定常流动的有效性     

喷动流化床流动特性的数值模拟

本文在实验研究的基础上,建立了喷动流化床流动特性的三维湍流数学模型。该模型可揭示气固两相喷动流化床的流动特性,弥补实验手段的不足。模型的气相流场采用欧拉方法,固相颗粒场采用离散单元直接数值模拟的方法,并以SIMPLE算法为基拙,对模型进行了求解。     

澳门周围水域潮汐流动的数值模拟

采用有限体积法求解二维浅水流动方程,建立了一种在非规则结构化网格上数值求解平面二维浅水流动的数值方法,其中应用Roe格式计算数值对流通量,并通过MUSCL途径插值界面状态,以对水位插值代替水深插值来消除地形变化对格式的影响;对源项采用分步法及半隐式求解,防止了源项引起的计算不稳定;处理动边界采用人工小水深的方法.应用于澳门水域潮流的模拟,计算结果与实测吻合较好.