用离散涡方法模拟喷嘴出口射流流场的数值模型

根据流体力学原理,用离散涡方法建立了在复平面上模拟喷嘴出口流场的一套完整模型。在此基础上,用ＦＯＲＴＲＡＮ语言编写了计算机程序,模拟了喷嘴出口流场的速度分布,利用涡元矢量图对计算结果进行了分析,得出了射流在喷嘴出口分离流动及流场中旋涡形成及发展的规律。射流在有锐缘的地方会发生流动分离,形成的剪切边界层逐渐卷起,形成旋涡结构。在喷嘴出口与哨嘴出口处形成的旋涡结构,在势流流场和其它涡元的作用下逐渐向下游运动,并向两边扩散。由于粘性扩散及计算误差的影响,在远喷距处会出现涡旋结构变形     

基于小波分析后台阶流动中的湍动能能谱分布

利用小波分析良好的时频双局域性特点,对准周期性非定常流场的速度信号进行了多辨析研究.使用大涡模拟方法得到了后台阶流动的流场,借助正交离散小波变换将流场中各处的速度信号在不同频带上逐层分解.提出了以特征频带名义频率和特征频带速度幅值作为描述流场各点处湍动能在不同频带上能谱分布的量化指标的新方法.以分析后台阶流动的脉动速度场为例,结果表明,所提出的方法可以为改进大涡模拟方法中的网格分布提供理论依据.     

建筑结构非定常绕流风场的数值模拟方法

给出了一种适用于计算建筑结构非定常绕流风场的大涡数值模拟算法.该算法基于有限差分法,采用曲线坐标结构网格,能精确描述形状复杂的物面边界,为下一步准确模拟含有因结构受风变形所致运动边界的绕流场奠定了基础.该算法采用投影法解耦纳维-斯托克斯方程中的压力和速度,对非定常流场的时间步进采用二阶Adams-Bashforth方法.采用同位网格以减少计算所需内存,为了平抑同位网格下中心差分格式导致的固有压力波动现象,计算对流速度时采用Rhie-Chow动量插值方法,利用编制的曲线坐标系下大涡模拟数值计算程序,对德州理工大学(TTU)建筑足尺模型绕流风场进行了模拟,所得结果与现场实测和风洞试验结果进行了比较.结果表明,本文算法是建筑结构非定常绕流风场数值模拟的有效方法.     

四阶紧致格式有限体积法湍流大涡模拟

为准确预测不可压复杂湍流,提出了一种应用于大涡模拟的高精度有限体积法。该方法在非交错网格上数值求解大涡模拟方程,空间离散采用有限体四阶紧致格式,时间推进用四阶Runge-Kutta法,压力速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动和振荡平板上方的流动,证实了该方法具有近四阶精度;并在此基础上,采用动力Smagorinsky亚格子应力模式,成功地完成了充分发展槽道湍流的大涡模拟计算,所得结果与直接数值模拟结果吻合良好。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

双拉格朗日模型模拟气固两相双圆柱绕流

以湍流剪切层中大尺度旋涡结构与颗粒群的相互作用及运动规律为研究对象,提出了一种基于离散涡方法与颗粒碰撞模型相耦合的双拉格朗日气固两相流模型.用次循环调整颗粒运动计算的时间步长,将两相计算耦合在一起,从而解决了离散涡方法求解单相流场与颗粒运动方程的时间相容性问题.首先用颗粒在兰金涡中的运动验证了所提的算法,随后将该方法应用于高雷诺数湍流旋涡结构中的颗粒水平输运过程研究,模拟了圆柱尾迹中颗粒的沉降过程及其被旋涡卷吸、扬起和团聚的过程,说明在有大尺度旋涡存在的流场中颗粒更容易被流体向下游榆运.当来流速度与重力方向一致时,颗粒分布特性受颗粒Stokes数(St)影响.当St小于1时,颗粒响应流场变化的时间短,随St增加颗粒对流场变化的响应变得迟钝,重力所起的作用更加明显;当St远远大于1时,颗粒的分布几乎不受旋涡的影响,而与自由沉降的分布十分接近.     

管束间复杂区域流场的离散涡数值模拟

对管束间复杂区域流场的数值模拟方法进行了分析,针对管束绕流的具体困难,对现有基于Lagrange框架的离散涡方法进行了改进并提出相应的数值算法.改进后的算法只需用较少的涡元数量,使管束绕流的快速计算成为可能.将数值模拟结果与实验进行对比,发现无论是瞬时流场还是平均速度和脉动速度分布的结果,两者都非常吻合.证明提出的离散涡方法模拟管束间的单相流动是完全可行的.     

180°方形弯管中湍流流动的大涡模拟

把大涡模拟应用于复杂几何形状的流域是一个重要研究课题。对１８０°方形弯管中的湍流流动作大涡模拟,以网格节点为中心的交错网格作为离散网格,用有限体积方法和中心差分进行空间离散,用Ａｄａｍｓ－Ｂａｓｈｆｏｒｔｈ格式进行时间离散,用简单的Ｓｍａｇｏｒｉｎｓｋｙ模型来作为次网格湍流模型。用多重网格方法求解压力方程以提高计算效率。数值模拟结果与实验结果取得了很好的一致。     

时间步长对格子玻尔兹曼法模拟室内气流精度的影响

基于格子玻尔兹曼法的大涡模拟（LBM-LES）是湍流模拟的新方法,但不恰当的时间步长δt可能会影响其计算精度。首先理论总结了δt可能对LBM-LES湍流模拟造成的影响,阐明过大的δt会导致速度场产生压缩性误差,而过小的δt会导致超松弛碰撞产生速度场的数值振荡。其次,通过对等温室内气流案例进行LBM-LES模拟,定量讨论了δt引起的压缩性误差和数值振荡问题。结果表明,δt较大时流场密度变化剧烈,且格子玻尔兹曼单位的马赫数（M）超过0.3的区域中速度场产生了明显的压缩性误差。同时,过小的δt导致平均及脉动风速均产生了数值振荡,这在网格分辨率较高时尤为明显。建议模拟时在确保δt足够小以满足最大风速区域的M<0.3的基础上,尽量增大δt以防止产生数值振荡。     

用离散涡方法模拟喷嘴出口射液流场的数值模型

根据流体力学原理,用离散涡方法建立了在复平面上模拟嘴出口流场的一套完整模型。在此基础上,用ＦＯＲＴＲＡＮ语言编写了计算机程序,模拟了喷嘴出口流场的速度分布,利用涡失矢量图对计算结果进行了分析,得出 射流嘴出口分离动及流场中旗旋涡形成发展的规律。     

湍流分离流中颗粒的扩散机制

后台阶流动包含分离流重要的流动特性, 采取欧拉-拉格朗日耦合算法对后台阶分离流动中颗粒扩散运动进行数值研究. 气相场采取大涡模拟方法, 亚格子模式基于标准的Smagorinsky 模式, 颗粒相运动采取轨道法模拟. 计算所得气相的流向平均速度和平均脉动速度与实验结果吻合较好, 验证了模型和方法的正确性. 基于此, 数值分析后台阶两相流动的特性以及流场涡结构的发展和演化过程. 结果表明: 两相流中颗粒的扩散特性既受到颗粒粒径的影响, 又与颗粒和涡结构的相互作用时间有关. 后台阶流场中增加结构物时, 流场涡结构发生变化, 即与扰动源保持一定距离后, 涡数量增多, 流场中颗粒分布不均匀, 较多颗粒聚集在涡的外缘.     

三维水翼初始梢涡空泡数的尺度效应

针对梢涡流场和初始梢涡空泡数的尺度效应问题,利用大涡模拟(LES)湍流模型对三维水翼的梢涡流场流速进行模拟计算;为减少误差,对梢涡流域网格进行了局部加密处理,对未发生空化时梢涡内的轴向速度和切向速度进行计算.结果表明,LES湍流模型的流场流速计算结果与实验值吻合较好.同时,介绍了经典初始梢涡空泡数尺度效应公式的推导过程,并利用数值计算的速度环量和涡核半径修正尺度效应公式.     

圆柱空腔内旋转流动涡破裂随Re变化的实验研究

针对旋转封闭圆柱空腔内的轴对称涡破裂现象,设计了可变旋转Ｒｅ实验系统。实验中空腔高径比Ｈ／Ｒ＝１．７８。采用激光相位多普勒测量技术测量了圆柱空腔轴线上轴向速度分量分布随Ｒｅ变化的情况,定量观测和分析了涡破裂现象的发生、发展过程,测量到涡破裂现象出现和消失这两种临界状态２。     

平面射流拟序结构的离散涡模拟

采用离散涡方法对平面射流流场进行了数值模拟。显示了平面射流中的大尺度涡结构的卷起、旋转和合并过程,模拟结果揭示了平面射流中的大尺度拟序结构的时空演化特征和发展变化规律。得到了定点的瞬时速度和流场的时均速度分布,时均速度场的计算结果和实验值吻合,瞬时速度分布揭示了平面射流拟序结构的非定常、不稳定的非线性特征。     

考虑旋转和曲率修正的螺旋桨梢涡流场模拟

应用显式代数雷诺应力湍流模型对螺旋桨尾流中梢涡流场分布进行了数值研究,为了避免过高地预报梢涡涡核内湍流黏性耗散,对湍流模型进行了旋转和曲率修正．应用全六面体网格对螺旋桨计算域进行网格划分,为了避免数值离散误差,对梢涡区域进行了网格加密处理．计算结果表明：提出的尾流中梢涡流场分布数值模拟方法能够准确预报螺旋桨梢涡流场的分布及涡核位置,并准确反映了梢涡形成和发展过程中梢涡内主涡和次涡的关系,与实验测量结果基本一致．     

混合层流动的数值模拟

对混合层流动进行了数值模拟。采用了离散涡方法。计算中采用扩散速度法处理粘性作用项。计算结果表明：与通常采用随机涡方法或核膨胀法处理粘性作用项相比,采用扩散速度法的计算结果与实验更吻合,同时更准确地反映了混合层流动的自相似性。在混合层流动中,邻近大尺度涡结构之间存在着连续不断的合并;采用二维离散涡方法预报的横向脉动速度均方根大于实验值是由二维离散涡方法忽略了真实流动的三维特性引起的。     

圆柱绕流远场涡对的数值模拟

基于离散涡方法求得非定常、不稳定流场,数值模拟了三种不同时刻高雷诺数下圆柱绕流结构的发展,从流谱图、等涡量线图和涡谱图可以清晰地看出从近场的初生卡门涡街,过渡到远场的二次涡街的过程,计算结果发现：远场离散涡有形成二个涡的涡对及三个涡的涡对的趋势,计算结果说明了流体运动中涡对结构的本质：由于来流是均匀的,没有加入任何拔动,当流体流过钝体时产生具有剧烈分离的不稳定流动,因此在远场形成的二次涡对及卡门涡     

离散涡法模拟建筑物绕流流场

应用离散涡数值方法对高雷诺数下建筑物绕流流场进行了数值模拟，得到了不同相对位置、相对距离的建筑物在不同来流方向的分离流场和涡量分布，总结出受建筑物影响的流函数变化规律，取得了较为合理和可靠的计算结果，说明离散涡方法是研究建筑物的绕流问题的有效手段     

轴流通风机叶顶区域流动的实验研究

将粒子图像测速(PIV)技术应用到低速轴流通风机实验台上,在设计工况下对轴流通风机转子叶顶区域的瞬态速度场进行了实验测量,对3个不同叶顶间隙高度,分别测量得到了在3个不同周向平面上的速度．同粒子多普勒测速仪(PDA)时均测量结果进行比较,认为PIV可以得到与PDA相近的测量精度．根据PIV实验测量结果,重点研究了叶顶泄漏流动的涡旋结构,叶顶泄漏涡位置的不稳定性和锁相平均后的叶顶泄漏涡涡心的轨迹,并将实验结果与原有模型进行了比较．结果表明,叶顶泄漏涡涡心的运行轨迹与原有模型得出的估算公式的计算结果接近。     

超临界雷诺数下圆柱分离旋涡运动的数值研究

本文在文献[1]的基础上,用离散涡模型与湍流边界层理论相结合的方法,研究了超临界雷诺数下圆柱突然起动后的分离旋涡运动。文中考虑了二次涡的影响和涡旋的粘性扩散。计算结果与实验结果相符。