类车体三维绕流场的数值模拟

轿车的空气动力学特性对其行驶的动力性、操纵稳定性有着重要影响，采用CFD技术对其性能进行预测，可以提高新产品技术性能，缩短开发周期。本文对类车体进行物理、数学建模、设定边界条件，计算其相关性能，并进行了定性定量分析，结果与参考文献相比较具有合理的一致性。     

高超声速双锥体绕流的数值计算与流场分析

 双锥几何外型的高超声速绕流问题涉及激波-激波相互干扰以及激波-边界层相互作用等复杂流动现象,是计算流体力学领域的热点问题,常被作为此类复杂流动的典型算例。本文从守恒型Navier-Stokes方程组出发,考虑了非平衡态气体的振动激发过程,采用高分辨率TVD有限差分格式研究了5种来流马赫数（15.56,10.34,9.59,8.20和8.06）,2种气体环境（纯N2和纯O2）,以及2种双锥外型（尖头体和钝头体）的气动热力学特性,较详细地给出了上述各种情况下壁面热流密度分布以及气动力系数等重要参数,所得结果与国外相关实验数据吻合较好,初步显示了本文所编程序的有效性。对差分格式精度和典型限制器（minmod、van Leer限制器）进行比较,发现数值计算所得分离区大小受格式耗散性影响很大,其中使用van Leer限制器所得到的分离区大小与实验结果吻合得较好。     

二维管道内交变流动的格子-Boltzmann方法模拟研究

采用格子-Boltzmann方法（LBM）对二维管道内不可压缩交变流动——Womersley流的流动与换热进行了数值模拟研究．结果表明：管道内压力梯度呈周期性变化，在周期较小时，管道的绝大部分区域速度分布相当平坦，在靠近壁面附近出现了速度的峰值，且靠近壁面处速度变化剧烈；在周期较大时，速度分布近似于抛物线，管道中心速度总大于壁面附近速度，且不同周期下的模拟结果与解析解都吻合得很好．另外，当速度交变时，温度表现为波动变化．速度与温度的相位差、温度的分布特性、温度的波动特性等都与压力梯度交变周期及幅值有关．研究结果表明，LBM可以用于交变流动的模拟．鉴于LBM的众多优点，该方法有望成为求解复杂流动与换热的一种有效的数值模拟手段．     

后置导叶型喷水推进泵总特性及三维流场

研究了带后置导叶的喷水推进轴流泵的三维流场,采用k-ε湍流模型和SIMPLEC算法对流场进行了数值模拟,分析了叶片表面的压力分布,模拟得到的总性能曲线与实验曲线比较吻合,模拟结果表明,导叶出口面上轴向速度与圆周速度的分布不均匀,与设计目的是一致的。     

超声速圆盘空化器超空泡流动数值计算方法

以超空泡射弹为研究背景,为分析在超声速条件下水的压缩性对圆盘空化器超空泡流场的影响,基于理想可压缩流体势流理论,提出了一套有限体积数值计算方法.采用可压缩流体满足的连续性方程和Tait状态方程,结合Riabouchinsky超空泡闭合模式,提出了针对超空泡流场反问题的一种求解方法.根据超空泡表面不可穿透条件,设计了一种新的超空泡外形迭代方式.在解决超声速圆盘空化器超空泡流场计算问题的基础上,分析了压缩性对超空泡形态和阻力系数的影响.在超声速条件下（马赫数为1.0~1.2）,流体压缩性将导致超空泡前后略微不对称,前端比尾端截面更窄,空泡最大截面略微向后移;在相同空泡数下,随着马赫数的增加,空泡长细比逐渐增加,压差阻力系数也不断增加.计算结果与经验公式及有关文献结果吻合较好.     

高雷诺数下绕钝体分离流动的一种新数值算法

本文提出一种新的高雷诺数下统钝体分离流动的数值计算方法.文中将流场分成邻近物面的内区和稍远离物面的外区两部分.内区用有限差分法求解,外区用涡法求解,并对内区向外区输运涡量的计算和内外区交界面边界条件的确定作了特殊的处理.计算结果表明本方法是正确有效的.     

交变流动蓄冷器阻力特性

交变流动的阻力特性不同于稳态流动,其阻力特性表现为动态参数的波动压降和相位差．本文给出了５０Ｈｚ交变流动蓄冷器的阻力特性实验结果．以雷诺数和无因次距离为准则数拟合了交变流动蓄冷器最大和平均压降因子的准则关系式．与稳态经验公式得出的压力降的值进行比较,发现交变流动蓄冷器的压力降为稳态流动的２至３倍．同时,给出了不同丝网目数对应的极值雷诺数和无因次相位差的准则关系式．     

一类振荡函数的数值积分方法

目的研究型如∫π－πｆ（ｘ）ｓｉｎ（ωｘ）ｄｘ（ω为正整数）的振荡函数的数值积分问题．方法用Ｐｅａｎｏ核和Ｔａｙｌｏｒ公式．结果给出了这种类型数值积分的两个计算公式和误差估计．结论数值计算表明,该计算公式与常用的Ｌｏｂａｔｏ法和Ｆｉｌｏｎ法相比具有计算量小和求积精度高的特点．     

任意旋成面跨音速叶栅流动的一种高效率解法

提出了任意旋成面叶栅跨音速绕流的一种高效率解法。导出了以Ｖｏｎ Ｍｉｓｓｅｓ 坐标为自变量的任意旋成面叶栅中流动的流线控制方程；简单而有效地处理了速度的双 值问题。用数值方法研究了流面半径和流片厚度随轴向变化对于流动的影响。     

交汇水流三维数值模型

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

基于大涡模拟的开孔潜体流噪声数值模拟

基于大涡模拟（LES）方法和Lighthill声学类比方法,对一倒T型开孔进行流场及流激噪声进行数值模拟,计算所得结果与公开发表文献上的数值及实验结果吻合良好,验证了基于LES的流噪声数值模拟方法的可行性.在此基础上,对某三维开孔潜体（原型）及其改型进行了流噪声数值模拟,对潜体内腔流场及其噪声辐射结果进行分析,从数值模拟的角度验证了改型潜体在阻力和噪声性能上都要优于原型潜体.
     

摇荡物体周围粘性流的数值模拟

提出了一种新的求解非定常不可压Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的数值方法，在每个时间步长上采用人工可压性技术，通过对动量方程取散度，从而获得离散化的Ｐｏｉｓｓｏｎ方程，当关于虚拟时间的迭代过程达到收敛时，便得到了无散速率场，采用该方法模拟了一个纵摇椭圆柱体周围的非定常不可压性流动。     

小幅振动旋成体跨音速绕流的数值分析

假设绕小幅振动旋成体所导致的物理量的变化为平均定常扰动量的高阶小量,那么将物理量(速度、密度及压力)摄动展开后代入Euler方程,便导出了前两阶摄动量的控制方程和相应的边界条件,用TVD格式进行数值求解,得到了平均定常跨音速流(零阶)和基频(一阶)非定常分量的数值解。     

行波磁场作用下金属液体运动行为的数值模拟

将行波磁场理论和流体运动N-S方程相结合,建立了行波磁场作用下材料制备用容器内三维流场的数学模型,并引入了新的电磁力计算模型.通过与文献中实验结果的对比,验证了本模型及计算程序的可靠性.研究还发现,行波磁场作用下装置内纵截面形成了两个对称的大环流;当径高比降低时,容器内的流场结构没有改变,仍然为两个对称的旋流,形状由椭圆形变成长方形.当电磁力增大时,轴向速度曲线由平顶变为尖顶的抛物线,并且电磁力越大,波顶离底部壁面的距离越近.     

液流通过节流孔的流场结构的数值解析

为了研究液压控制阀的动态特性,对从初始稳定流开始的液流通过节流孔的过渡过程进行了数值解析计算。对压力差分别为成倍增加、减少至零、反方向阶跃变化3种情况下的管内液流的过渡过程进行仿真,作出了流场随时间变化的流线图。用数值计算结果来研究流场结构及变量的变化。     

渗流问题数值模拟计算的新方法

采用边界拟合坐标系，取渗流区域的边界为坐标线（二维问题）或坐标面（三维问题），用 Ｔｈｏｍｐｓｏｎ法将不规则渗流域转变为规则域，从而在规则域上求解转变后的渗流方程．避免了在不规则域上划分单元及迭代运算中生成新网格的繁琐过程，尤其对三维有自由面的无压渗流问题，本方法具有精度高而工作量小等优点。     

二维圆柱层流绕流及其控制数值模拟

采用有限体积法求解二维N-S方程,通过改变时间、空间计算参数,对雷诺数100的二维圆柱非定常流场进行了数值模拟,将所得结果与实验数据以及已有的计算结果进行了对比.并对在尾部多个位置加入不同长度阻隔板控制流动的圆柱流场进行了考察.结果发现:在合适的位置放置阻隔板,可以有效抑制尾迹中涡脱落的形成.     

微型环形燃烧室数值模拟及优化设计

为了优化微型环形燃烧室的设计方案,利用流体分析软件Fluent,采用RNGk-ε双方程湍流模型分别对两个微型环形燃烧室的内部流场进行了三维的冷态数值模拟,得出了它们的内部流场分布。两个模型分别是采用了航空发动机大型燃烧室设计理论及公式设计而成的燃烧室模型,和已经在实际中成功应用的国外某型微型涡喷发动机燃烧室模型。该文分析两个模型各燃区内计算流场的差异,确定了满足设计性能要求的目标流场,并对原有设计参数进行局部调整,使其流场组织向目标流场靠近。通过此方法,确定了燃烧室优化后的设计参数。     

天然河道岸边排污的三维数值模拟

对某磷肥厂废水排入某江干流－河段形成的流场和浓度场的进行了三维数值模拟，并采用原观测资料对数学模型和计算方法进行了验证，三维计算中采用了实际的地形资料和标准k－ε模型方程组，成功地模拟了天然河流岸边水面排放形成的流场与浓度场，计算结果与原型观测数据具有良好的一致性。