隧道中高速列车空气阻力的数值计算

对高速列车在隧道中行驶时所引起的湍流场进行了数值模拟，计算了列车所受到的空气阻力。数值计算采用变密度、轴对称的ｋ～ε模型和ＰＩＳＯ算法，并将粗糙面速度分布公式和尼古拉兹公式引入壁函数法中，以反映壁面粗糙度的影响，提出了对不规则边界的速度修正法。针对不同的车速、阻塞比和车头细长比，分别进行了模拟，获得了全场的压力分布，以及空气阻力与基本参数的关系，并与传统的经典方法阻力计算结果进行了对比分析。     

在1.4米航空风洞中模拟大气边界层

在风洞中正确模拟大气边界层的流动特性是风工程风洞试验结果可信的必要条件,本次试验研究的目的是在短试验段的航空风洞中建立大比例的大气边界层模拟流场。通过适当的方式延长1.4m×1.4m航空风洞的试验段长度,并利用尖塔、粗糙元等边界层发生装置,在该风洞中建立了边界层流场,测量了流场的平均风速剖面、湍流强度剖面、脉动风速的自相关系数、风谱等参数,讨论了湍流积分尺度的处理和大气边界层几何模拟比例的确定,用谱拟合法和自相关系数积分法求出了湍流积分尺度。结果分析表明:试验所得流场是合理的大气边界层模拟流场,其平均风速剖面幂指数α=0.3,大气边界层模拟比例为1∶500,为后续的建筑物模型动态风荷载试验提供了前提条件。     

一类平衡大气边界层边界条件研究

基于标准k-ε湍流模型,首先利用湍流粘度方程和剪切应力在整个边界层内恒定的假设,推导出一类耗散率表达式,并根据常用的湍动能入口剖面方程以及平均风速剖面方程,计算获得相应的耗散率方程;然后在输运方程中添加自定义源项,通过已经确定的平均速度方程、湍动能方程、耗散率方程计算得到相应输运方程的自定义源项表达式,并进行空风洞数值模拟,从而得到了一类满足平衡大气边界层的来流边界条件．通过将这种边界条件与由湍流平衡条件得到的边界条件进行比较,表明本方法获得的边界条件更适用．并且,本方法无需考虑修正壁面函数和修正湍流模型常数,因而计算更为简单,可为平衡大气边界层的研究提供一种新的思路．     

边长为6m的立方体其周围流场的数值研究

使用CFD方法研究立方体周围流场.Silsoe研究院建造了一个边长为6 m的立方体作为研究钝体周围的压力场和速度场的实验对象.相应于这个立方体的全尺度实测结果、风洞实验结果和数值计算结果在很多文献中可以查到.本文使用标准k-ε模型和RNGk-ε模型,考察当来流速度为对数率速度剖面时,不同的来流湍动能设置对计算结果的影响,所得计算结果与已有的全尺度实测结果作比较,得到了有意义的结论.     

地面效应影响下的高速火箭橇气动特性数值与风洞试验研究

精确预示地面效应下高速火箭橇的气动特性及流场规律对高速火箭橇的设计和评估具有重要意义。本文应用有限体积方法,研究了湍流模型对火箭橇气动特性计算精度的影响,建立了基于realizable k-ε湍流模型的火箭橇气动特性的高精度数值计算方法;结合风洞试验方法,研究了雷诺数和地面效应对高速火箭橇流场流动规律的影响,分析了火箭橇气动特性。结果表明,火箭橇阻力系数随雷诺数增大而减小,升力系数和俯仰力矩系数随雷诺数增大而增大,但雷诺数对高速火箭橇气动特性的影响较小;地面效应会使火箭橇流场发生激波-激波干扰、激波-边界层干扰和激波反复反射等复杂气动现象,大幅提升了火箭橇的升力系数和俯仰力矩系数,但对阻力系数的影响较小。研究为高速火箭橇气动外形的设计及运动稳定性的评估提供依据。     

空化、超空化流动的数值模拟方法研究

基于结构化网格,运用可压缩流N—S方程及k-ε湍流模型对流场进行求解,在低压区域引入一种基于混合密度函数的空化模型对轴对称体的空化、超空化流动进行了数值模拟．通过将半球圆柱的计算结果与实验数据和前人的计算结果进行对比,验证了所发展的数值方法的可靠性．最后,采用非定常的数值方法,研究了钝头体射弹的空化、超空化流动特性,并模拟了其超空泡的发展过程．     

三维湍流流动的涡量-向量势法数值模拟

本文采用涡量-向量势法,结合常规的 k-ε双方程湍流模型,并运用重新推导的固体壁面涡量边值计算法,对扁平射流煤粉燃烧器内的冷态三维流场进行了数值模拟。计算结果与实验数据相吻合。本文工作证明,用涡量-向量势法求解三维湍流问题是可行的,该法具有不应忽视的优越性,值得进一步开发。     

一种用于分析封闭腔内湍流自然对流换热的k-ε新模型

由于目前用于求解湍流自然对流流动与传热的k-ε模型在应用过程中存在不足,结合高雷诺数k-ε模型需要借助壁面函数法来确定壁面上相关参数值和低雷诺数k-ε模型在近壁区布置更多节点以便获得粘性底层详细信息的特点,重新定义了湍流普朗特数σt的计算式,提出了一种修正的k-ε新模型;利用该模型对封闭方腔内的湍流自然对流流动与传热进行了数值分析。结果表明:与文献中数值模拟结果相比,当108≤Ra≤1014时本文模型所得壁面平均努塞尔特数更接近文献中的实验值,与实验值之间的相对误差在8%以内;壁面的局部努塞尔特数与文献中的实验值吻合得较好。这说明本文模型用于求解封闭腔内湍流流动与传热问题是合适的,比其它湍流模型更能准确地描述封闭腔内湍流自然对流换热中边界层发展与壁面传热特性之间的内在联系。     

物面引射对滑移边界条件的影响

讨论壁面有引射的滑移边界条件的提法．利用Ｃｈａｐｍａｎ Ｅｎｓｋｏｇ速度分布函数,通过分析Ｋｎｕｄｓｅｎ层外缘和壁面处的质量、动量和能量等通量守恒,得到了有壁面引射的多组元气体有催化反应的壁面滑移边界条件方程组．利用粘性激波层方程和所得的边界条件对钝头体绕流驻点区的流场进行了计算,讨论了物面引射对边界上诸量及流场的影响     

k-ε双方程湍流模型对制退机内流场计算的适用性分析

为研究不同双方程湍流模型对制退机内复杂流场计算的适用性,以某火炮制退机为研究对象,建立了实际结构下的三维计算模型,利用动网格与滑移网格技术,实现了火炮实际后坐速度下的制退机内部三维运动流场的数值计算。分别采用标准k-ε模型、RNGk-ε模型和Realizable k-ε模型计算制退机内部各腔室压力,与实验曲线对比,结果表明,应用标准k-ε模型对后坐冲击过程的制退机内部压力计算的误差最小,与实验结果吻合最好。     

基于滑移速度壁模型的复杂边界湍流大涡模拟

采用滑移速度壁模型实现了浸入边界方法与壁模型相结合的大涡模拟.本文首先分别采用平衡层模型和非平衡壁模型对周期山状流进行数值模拟,以考查在壁模型中考虑切向压力梯度的作用.数值结果表明,流场的压力对本文所采用的壁模型形式并不敏感,但是考虑切向压力梯度可以显著改进壁面摩擦力的计算结果,并且能够准确的预测强压力梯度区以及分离区内的流动平均统计特性.不考虑压力梯度效应的平衡层模型显著低估了壁面摩擦力的分布,同时无法准确预测分离区内的平均速度剖面.非平衡模型的修正项正比于切向压力梯度和壁面法向距离,因此在强压力梯度区或者网格较粗时,计算得到的平均压力和摩擦力分布以及流动的低阶统计量均与参考的实验和计算结果吻合.在此基础上,通过回转体绕流的大涡模拟考查了该方法用于模拟高雷诺数壁湍流的适用性,非平衡壁模型可以准确地捕捉流动的物理结构并较准确地预测其水动力学特性.结果表明,将浸入边界方法与非平衡滑移速度壁模型相结合的大涡模拟,有望成为数值模拟复杂边界高雷诺数壁湍流的工具.      

床面上直立圆柱的三维湍流数值模拟

从数值预报桥墩等结构物床面局部冲刷的角度发展绕直立圆柱的三维湍流的数值模拟技术. 基于Wilcox的k-ω两方程湍流模式，采用基于有限体积法的压力修正SIMPLE算法, 计算了绕床面直立圆柱的三维湍流流场，分析了光滑和粗糙床面两种情况下的流动情况. 通过系列的验证计算，表明该计算模型能够比较准确地反映不同外来流条件下绕直立圆柱的流场. 计算结果揭示了床面粗糙度对绕圆柱的湍流流动的影响.     

多孔介质壁面剪切湍流速度时空关联的研究

作为一个基础统计量,时空关联函数在湍流问题的研究中有着广泛的应用,是研究湍流噪声、湍流中物质扩散和大涡模拟亚格子模型等问题的重要参考.本文通过建立三维多孔结构壁面剪切湍流模型,采用含Darcy-Brinkman-Forchheimer作用力项的格子Boltzmann方程对无穷大多孔介质平行板之间壁湍流进行了数值模拟,进而研究其速度脉动时空关联函数的统计特性.一方面,根据计算得到的流场数据,对比分析了常规槽道湍流与多孔介质壁面槽道湍流的时间关联函数.另一方面,计算并讨论了不同孔隙率和渗透率的多孔介质壁面对速度脉动时空关联性的影响.通过研究表明:多孔结构壁面剪切湍流的时空关联函数等值线与椭圆理论相符;在研究参数范围内,多孔介质壁面的速度时空关联系数随着孔隙率增大而增大,随着渗透率增大而减小.同时发现在槽道壁面的近壁区、过渡区、对数律区和中心区等不同位置处,速度时空关联呈现较大差异性:越远离壁面位置(对数律区和中心区),其时空关联函数所呈现的关联等值线椭圆越细长,高值相关等值线越集中.多孔介质主要改变速度时空关联椭圆图像的椭圆率,说明多孔介质壁面主要影响湍流横扫速度.     

基于滑移速度壁模型的复杂边界湍流大涡模拟

采用滑移速度壁模型实现了浸入边界方法与壁模型相结合的大涡模拟.本文首先分别采用平衡层模型和非平衡壁模型对周期山状流进行数值模拟,以考查在壁模型中考虑切向压力梯度的作用.数值结果表明,流场的压力对本文所采用的壁模型形式并不敏感,但是考虑切向压力梯度可以显著改进壁面摩擦力的计算结果,并且能够准确的预测强压力梯度区以及分离区内的流动平均统计特性.不考虑压力梯度效应的平衡层模型显著低估了壁面摩擦力的分布,同时无法准确预测分离区内的平均速度剖面.非平衡模型的修正项正比于切向压力梯度和壁面法向距离,因此在强压力梯度区或者网格较粗时,计算得到的平均压力和摩擦力分布以及流动的低阶统计量均与参考的实验和计算结果吻合.在此基础上,通过回转体绕流的大涡模拟考查了该方法用于模拟高雷诺数壁湍流的适用性,非平衡壁模型可以准确地捕捉流动的物理结构并较准确地预测其水动力学特性.结果表明,将浸入边界方法与非平衡滑移速度壁模型相结合的大涡模拟,有望成为数值模拟复杂边界高雷诺数壁湍流的工具.     

射流元件控制通道流动的三维数值模拟

针对一种射流元件控制通道的复杂结构 ,采用分块对接技术和网格“融合”技术生成计算网格 ,并运用五步显式格式的 Runge-Kutta法和多重网格法对含全 N-S方程、RNG k-ε湍流模型和两层分区壁面模型的流动模型进行数值求解。通过对控制通道内部流动的数值模拟和流场特性分析 ,提出了改进方案     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

基于人工神经网络的湍流大涡模拟方法

大涡模拟方法(LES)是研究复杂湍流问题的重要工具,在航空航天、湍流燃烧、气动声学、大气边界层等众多工程领域中具有广泛的应用前景.大涡模拟方法采用粗网格计算大尺度上的湍流结构,并用亚格子(SGS)模型近似表达滤波尺度以下的流动结构对大尺度流场的作用.传统的亚格子模型由于只利用了单点流场信息和简单的函数关系,在先验验证中...     

k-ε双方程湍流模型对制退机内流场计算的适用性分析

为研究不同双方程湍流模型对制退机内复杂流场计算的适用性,以某火炮制堪机为研究对象,建立了实际结构下的三维计算模型,利用动网格与滑移网格技术.实现了火炮实际后坐速度下的制退机内部三维运动流场的数值计算.分别采用标准k-ε模型、RNGk-ε模型和Realizablek-ε模型计算制退机内部各腔室压力,与实验曲线对比,结果表...     

物面引射对滑移边界条件的影响1)

讨论壁面有引射的滑移边界条件的提法.利用Chapman-Enskog速度分布函数,通过分析Knudsen层外缘和壁面处的质量、动量和能量等通量守恒,得到了有壁面引射的多组元气体有催化反应的 壁面滑移边界条件方程组.利用粘性激波层方程和所得的边界条件对钝头体绕流驻点区的流场进行了计 算,讨论了物面引射对边界上诸量及流场的影响.