迎风型TVD格式在高速进气道粘性流动模拟中的应用

迎风型ＴＶＤ格式在高速进气道粘性流动模拟中的应用刘晶昌，徐建中（中国科学院工程热物理研究所北京１０００８０）邢君波（航天部３１所北京１０００８０）关键词超音速进气道，数值模拟，Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程１前言高速进气道流动做为一种非常复杂的流动，...     

二维卡门涡街的格子Boltzmann仿真

格子气自动机和格子Boltzmann方法的迅速发展提供了一类求解流体力学问题的新的方法。本文中，我们介绍了Boltzmann方法，解决了格子气方法中的缺点，通过选择适当平衡分布及其参数，导出了Navier－Stokers方程，并得到了声速和粘性系数。最后在微机上模拟了在无限长平板流动问题及绕单分离板的流动问题，得到了卡门涡街。结果表该模型有格子气方法及其它的数值方法所没有的优点，计算更精确、更直观、更有效。     

翼型大攻角低速分离流动的数值模拟

应用高收敛率、高精度和高分辨率的数值计算方法,通过求解非定常、可压缩雷诺平均的Ｎａｖｉｅｒ—Ｓｔｏｋｅｓ方程和ｑ-ω低雷诺数双方程湍流模型,数值预测了Ｓ８０９翼型在０－７０°攻角范围内大尺度分离与失速流场的流动与损失特性。     

气固两相三维圆柱绕流的直接数值模拟

本文采用直接数值模拟(DNS)方法研究了气周两相三维圆柱绕流的涡量场和颗粒扩散,并着重讨论了圆柱绕流中卡门涡街的形成和涡结构的转捩过程。同时分析了圆柱绕流中不同Stokes数的颗粒在涡结构作用下的横向扩散。结果显示Stokes数为1的颗粒主要分布在流场大尺度涡结构的外边界,而Stokes数为0．01的颗粒在涡结构的作用下,在流场中充分混合。     

绕圆柱体自由表面磁流体流动和传热的研究

本文对在不同雷诺数下,绕圆柱体的磁流体自由表面流动及传热进行了模拟,分析了磁场对绕流圆柱尾迹和涡分离的影响,获得了两种雷诺数下的电磁力密度、流场和温度场分布。结果表明,磁场不仅影响了流动的形态,而且对湍流有抑制作用,降低了自由表面的更新机制,从而减少了传热能力;在相同的Hartmann数下,相比低雷诺数下的流动换热情况,高雷诺数下的湍流不能被完全抑制,自由表面与尾迹的相互作用也较强,因而自由表面换热也较强。     

粘性涡方法在多圆柱绕流数值模拟中的应用

用粘性涡方法求解了涡量速度形式的二维不可压缩Navier—Stokes方程．利用分步方法分别求解无粘的对流运动和粘性扩散，粘性扩散用随机走步方法来处理，利用涡量交换技术来满足无滑移边界条件，同时为进一步减少计算量，在求解对流速度时，采用了多极子展开法加速技术使得计算量由N2的量级减少到NlogN的量级．数值模拟了单圆柱和双圆柱纵向排列的绕流流场．     

数值模拟圆柱绕流旋涡运动及尾流不稳定性分析

1引言流体绕过圆柱所产生的非定常旋涡运动以及由此引起的流动不稳定性在理论和实践上都具有重要的意义。数值模拟圆柱绕流旋涡产生及演化过程，探讨圆柱尾流涡街产生的机制，控制尾迹不同速度型以抑制涡街的产生，避免涡激振动在工程上造成破坏作用具有重大实际意义。为使问题简化，本文以二维圆柱绕流作为研究对象。该流动涉及到非定常分离，旋涡的形成、运动及发展，流动不稳定性质改变等许多未完全解决的问题。B。had等山对圆柱突然起动问题作了一系列实验研究。文献门对二维圆柱绕流问题作了系统数值研究。本文采用文献门提出的差分格…     

部分流泵整机非定常流动数值模拟

本文在设计工况下采用滑移网格技术对部分流泵进行了整机非定常流动数值模拟,分别分析了叶轮、蜗壳内的非定常流动规律。计算表明,部分流泵内流的非定常特性极其强烈,只有采用非定常计算方法才能反映其内流的实质。本文的计算为进一步研究部分流泵的内流现象、提高效率、减少水利损失提供了一定的理论依据。     

三维十五点格子Boltzmann模型仿真

格子气和格子Ｂｏｌｔｚｍａｎｎ方法的迅速发展提供了一类求解流体力学问题的新方法。格子Ｂｏｌｔｚｍａｎｎ方法在保留了格子气模型优点的同时,克服了它的不足之处。本文讨论了一种三维十五点格子Ｂｏｌｔｚｍａｎｎ模型,通过选择适当的平衡分布及参数,并用Ｃｈａｐｍａｎ－Ｅｎｓｋｏｇ展开和多尺度技术导出了Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程．在微机上模拟了工程中比较常见的管排绕流问题,并与实验观察到的结果进行了比较,结果表明该模型能较好的模拟复杂流动现象,并具有较好的工程应用背景。     

A curved no-slip boundary condition for the lattice Boltzmann method

We present a generalization of the no-slip boundary condition by Lätt et al. [J. Lätt, B. Chopard, O. Malaspinas, M. Deville, A. Michler, Straight velocity boundaries in the lattice Boltzmann method, Physical Review E 77 (5) (2008) 056703] from straight to curved geometries for the lattice Boltzmann Bhatnager–Gross–Krook method (LBGK). The boundary condition is based on a reconstruction of the populations from the density, velocity and rate of strain. For curved boundaries, the reconstruction reduces the question of accuracy to a technical issue of interpolation. We present a method of interpolation allowing a very accurate representation of the curved boundary. The resulting boundary condition is verified for three different test cases: Taylor–Couette flow in-between rotating cylinders, laminar flow around a cylinder and flow past an impulsively started cylinder, demonstrating its second order accuracy and low error constant. The present boundary is stable for relaxation frequencies close to two.     

复合管道内液-固两相流浆体水击压强计算

对于单一管道的单相液体水击现象计算已经比较成熟,但是很少有对于复合管道的液固两相流系统的浆体水击 计算方面的报道。根据实际工程设计需要,提出了具有复合管道的液-固两相流浆体水击压强的计算方法。同时对某电厂 液一固两相流输灰系统浆体水击压强进行了计算与分析,提出了防止浆体水击的措施。     

用SPIV技术测量压气机转子尖区复杂流动

本文将数字式SPIV技术应用到低速大尺寸压气机实验台上,并在设计状态和近失速状态下,对转子槽道内叶尖区域多个截面内的三维瞬态速度场进行了成功测量。测量结果表明SPIV透过机匣视窗及即可直接测量旋涡和二次流的瞬态结构,而且能够解决多级压气机内部流场的测量.在压气机SPIV测量中,激光反光的控制和示踪粒子的均匀稳定散播是决定成败的关键因素。     

气液两相流体绕方柱流动的数值模拟

本文详细给出了气液两相流动双流体模型下的各控制方程组和相间作用力,通过将气泡引起的紊流与剪切引起的紊流线性叠加给出了气液两相流动的紊流模型。通过给出的边界条件,采用ＩＴＡ－ＩＰＳＡ算法对垂直上升的气液两相流横向冲刷方柱引起漩涡证替脱落的过程进行了数值模拟,得到了流场和局部含气率变化的示意图,数值模拟与试验的结果基本一致。     

弯曲时片压气机时栅内二次流的数值研究

本文应用Ｂｅａｍ－Ｗａｒｍｉｎｇ近似隐式因子分解格式以及ＭＭＬ代数湍流模型,采用拟压缩性方法求解雷诺平均拟压缩Ｎ－Ｓ方程组,对弯曲叶片压气机叶栅内的三维粘性流场进行了数值研究．结果表明,正弯曲叶栅内通道涡较直叶栅强,其诱导产生的壁面涡较弱且位于吸力面与壁端的角区内,与壁角涡接近．在正弯曲叶栅出口处,通道涡开始处于破裂前期,叶栅总损失增加。反弯曲叶栅通道涡较弱,其诱导的壁面涡很强,位于通道涡左上方,壁面涡和通道涡的有利作用,使通道涡更稳定,叶栅总损失比直叶栅和正弯曲叶栅要小．     

纳米尺度圆柱绕流现象的分子动力学模拟研究

本文采用分子动力学模拟方法,利用Lennard-Jones(L-J)势能模型,模拟了Re为28时Ar流体流过Pt纳米圆柱的绕流现象.采用小步长时均值作为涡的瞬时值的方法,在纳米尺度、纳秒量级下得到了涡的周期性产生、组合、发展和脱落现象;在大步长时均条件下,得到了稳定的对称涡,表现出了绕流现象在不同时间范围内的不同特征.绕流过程还体现了流体的密度变化,圆柱上游密度大于下游密度、对称的两侧离轴线越远密度越大.结果表明,分子动力学模拟方法可以有效地细节刻画圆柱绕流现象.     

加力室隔热屏流场计算

加力室隔热屏流场计算赵坚行，刘全忠（南京航空航天大学动力工程系南京２１００１６）关键词隔热屏，改值计算，紊流反应流１引言本文采用数值计算的方法模拟带有隔热屏、外冷却气流、尾喷口的加力室热态流场。计算中采用修正卜。紊流模型来预估粘性系数。燃烧模型采用卜...     

LB方法的一个简单多速模型

研究一个简单的二维三速格子Boltzmann模型。将微观的离散模型与宏观的连续统模型相结合,以宏观的流体力学方程作为约束条件确定局部平衡分布函数。运用Chapmann-Enskog展开法和多尺度技术导出Euler方程和Navier-Stokes方程以及模型的粘性系数表达式。