波瓣混合结构三维流场数值计算

在贴体曲线坐标系上，采用有限控制体的非交错网格设置，对波瓣混合结构的三维混合流场进行了数值研究。经计算得到了波瓣混合结构气流混合的不同轴向截面上的速度分布，可以看出，由于波瓣尾缘处轴向旋涡的作用，主流和次流形成高效率的对流型混合，其轴向速度分布型面发生了明显的变化，气流间梯度较大的参数分布趋于均匀，这些特征与实验结果是一致的。     

槽道湍流相干结构的研究

采用高阶差分方法求解非定常不可压Navier-Stokes方程，进行槽道流相干结构的数值研究，与常用的谱方法相比，不需限制边界条件，可用于更一般的粘性流体的流动问题，文中所用时间分裂法和三维耦合高阶差分格式适用于包括邻近边界点在内的整个计算区域，克服了三维各自用高阶中心差分格式不适用邻近边界点的困难。通过算例模拟了槽道流动中单个相干结构的演化发展规律，结果与实验一致。     

三维曲线坐标系N—S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错风格设置，该算法双ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的了出现，在压力修正方程中对压力梯度的离散附加一项１－δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

一种适合全流速的流场数值方法

发展一种适合于所有流速的压力修正方法,在考虑了密度的可压缩性后,压力修正方程夸为对流扩散模型的,从而与般标量一样使用二阶中心差分离散。整个算法具有二阶精度,进一步使用多步压力修正来增强网格非正交时的收敛性态,在对流项离散中,设计了独特的人工耗散,也同样适用于所有流速,利用本算法求解全三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,计算采用了非交错网格方案,使用Ｒｈｉｅ－Ｃｈｏｗ方法消除了速度－压力的不关联,     

二维突扩管流的有限分析解

运用有限分析法求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，获得了二维突扩管流的数值解。对于雷诺数为２５，５０，１００，１５０，２００的来流，预报的流动图谱具有相似性，回流区域及再附距离的计算也获得了比较满意的结果     

分区对接网格下多段翼型N-S方程计算

对多段翼型的复杂粘性流场，各翼段物面及尾迹区生成C型网格．外场采用O型网格，中间用H型附加区过渡．自动生成高效的分区对接C-O-H型网格。采用时间精确有限体积高精度格式解可压缩N—S方程。层流、紊流的N—S方程求解结果，说明了本方法的有效性。     

超声速和高超声速粘性绕流QPNS方程的GLS有限元解

从完全非定常Ｎ－Ｓ方程（ＦＮＳ）中略去流向粘性导数导出拟抛物Ｎ－Ｓ方程（ＱＰＮＳ），ＱＰＮＳ技术结合经典抛物推进和拟时松弛，拟时松弛能有效地抑制解的发散。本文导出的熵变量形式ＱＰＮＳ方程具有对称性和自动满足热力第二定律，这将提高解的稳定性，伽辽金／最小二乘法（ＧＬＳ）被用来构造ＱＰＮＳ的弱解式。对于钝前缘物体，ＱＰＮＳ在前缘区并不适用，本文代之以用ＦＮＳ求解该区，由离散解得到的非对称线性方程组，对     

双重互易边界元法在求解N—S方程中的应用

提出了用双互易边界元法求解Ｎ－Ｓ方程的新思路。用Ｐｅａｃｅｍａｎ－Ｒａｃｈｆｏｒｄ算子分理解将时间相依的Ｎ－Ｓ方程分解为线性和非线性子问题，线性问题用共轭梯度法解除压强－速度的耦合；非线性问题进行局部线化。对所得的Ｐａｓｓｉｏｎ方程及相关类型方程采用双重互易的边界元法求解，消除了传统Ｎ－Ｓ方程边界元解法的区域积分问题。     

热完全气体的热力学特性及其N-S方程的求解

首先用五次多项式拟合给出了温度在50～3000K范围内的热完全空气的焓值与温度之间的函数关系式，导出了其它热力参数e，cp，cv和γ的表达式。接着提出了热完全空气总温、总压的计算方法，并将其计算结果与量热完全空气的结果进行了比较。最后，将该热完全空气模型用于N-S方程求解，对NAPA软件进行了改进，并用该软件时高马赫数钝体绕流流场和乘波体流场进行了计算，分析了气体真实效应时流场结构及参数的影响。结果表明，本文提出的总温、总压计算方法及时NAPA软件的改进是成功的，可较准确地模拟高超声速流动的主要特征。     

汽车风洞试验段流场的试验研究

以上海地面交通工具风洞中心全尺寸气动声学风洞1:15模型风洞为研究对象,采用三维热线风速仪,测量了该模型风洞内不同工况下试验段内流场的分布情况.通过所测流场内平均速度的分析,了解该模型风洞试验段内气流的定常流动形态,对比典型射流结构进行分析,得到了该开口式低速汽车风洞试验段内喷口处射流剪切层所特有的结构形态;通过对风洞试验段内不同截面处流场内部湍流强度的分析比较,了解了喷口射流剪切层内以及收集口处的湍流强度大小与分布情况.     

燃烧室冷态流场的数值研究

在任意曲线坐标系下对环形燃烧室三维冷态流场进行数值模拟。采用偏微分方程法和区域法生成三维贴体网格,紊流模型采用标准k-ε模型和重整化群k-(εRNGk-ε)模型,在非交错网格坐标系中采用SIMPLE算法求解。数值分析两种紊流模型对火焰筒内气流流动的影响,并与实验结果相比较。结果表明,与标准k-ε模型相比,RNGk-ε模型更适用于模拟紊流流动。     

粘性非结构网格生成方法及绕流计算

运用非结构网格进行计算时,需要在物体附面层区域内生成高展弦比的三角形网格。本文给出了一种基于物理学中弹簧原理的粘性非结构网格生成新方法－－挤压法,它能有效地实现粘性区域内位伸网格到附面层外一般非拉伸网格的光滑过渡。然后在非结构粘性网格上对Ｎ－Ｓ方程运用有限体积法进行了数值计算,ＮＡＣＡ－００１２翼型以及平板层流的数值试验结果较为满意。     

非正交曲线坐标湍流族流的数值模拟

对锥形燃烧室的冷态流场进行了数值模拟，采用Ｔｈｏｍｐｓｏｎ的非正交贴你坐标系统，用坐标变换的方法处理锥形燃烧室曲线壁面边界，用ＳＩＭＰＬＥ算法求解连续方程和Ｎ－Ｓ方程。亲流模型采用工程上常用的ｋ－ε模型。还运用贴体坐标系统数值研究了不同旋流数和不同进口紊流动能对湍流旋流流动的影响，并与实验结果相比较获得较满意的结果，对新型燃烧室的设计提供了理论依据。     

低速高湍流度90°弯管流动数值模拟

分析了曲率对弯曲管道流动的影响 ,给出了曲率修正的双层 k-ε湍流模型 ,并数值模拟了非均匀来流低速高湍流度 90°弯管内流动。文中比较了曲率修正双层 k- ε湍流模型的计算结果与实验结果 ,以及双层 k- ε湍流模型和逆压力梯度修正双层 k- ε湍流模型的计算结果。通过比较 ,发现经过曲率修正后的双层 k- ε湍流模型更好地模拟了非均匀来流低速高湍流度 90°弯管内流动。同时比较了壁面参数的插值方法 ,发现在弯曲管道流动模拟中对壁面参数进行一阶插值比零阶插值能更好地模拟真实流场     

风雨耦合下大型冷却塔流场特性与表面气动力

现行冷却塔结构抗风设计均忽略了降雨带来的影响,但在强风暴雨极端气候条件下,暴雨亦会直接影响塔筒内、外表面气动力并改变脉动风的湍流效应,而传统研究大多仅关注风驱动雨对于结构表面的冲击效应。为解决该问题,以国内已建世界最高220 m大型冷却塔为例,以风-雨双向耦合算法为核心,首先采用计算流体动力学(Computational fluid dynamics,CFD)技术对冷却塔周围风场进行数值模拟,并将表面风压分布与规范及实测曲线进行对比验证了风场模拟的有效性;再添加离散相模型(Discrete phase model,DPM)并进行雨滴和风场的同步迭代计算。在此基础上,系统研究了塔筒内外表面风驱雨量、雨滴附加作用力和雨致压力系数等影响规律,揭示了风雨场中塔筒表面速度流线、湍动能强度、雨滴运行速度和轨迹的作用机理。最终提出了基于风雨双向耦合算法的风-雨致等效压力系数新模型及其分布特性。研究结论可为此类冷却塔在极端气候下的表面荷载取值提供参考。     

串列叶栅粘性紊流流场有限分析数值研究

本文将近年发展起来的有限分析方法(Finite Analytic Method)用于曲线坐标系上紊流N-S方程的数值计算。计算了来流雷诺数为2.0×10~5,气流攻角分别为0°,10°,-10°三种情况的二维粘性紊流串列叶栅流场。文中用k—ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。数值计算结果与实验结果相比较,吻合程度令人满意。     

圆柱体尾流结构的控制

本文根据水洞观察的结果,分析了圆柱体尾流中,卡门涡街形成的原因是由于流动的不稳定性所产生。根据这一原则,找到了两种控制尾流结构的方法。当雷诺数小于大约200的情况下,在圆柱体“死水区”一侧的边沿,放置一个直径很小的控制圆柱体,可以抑制卡门涡街的产生。另一种方法是在尾流的平均轴线上,放置一个宽度与圆柱体直径相等的平板,同样能够抑制涡街的产生。风洞试验的结果表明,即使在高雷诺数下(亚临界范围内),由于尾流中放置了平板,使得非定常旋涡受到抑制,从而显著地改变了圆柱体上的气动力。     

多普勒激光测速仪在轴流风扇流场测量中的应用

介绍了多普勒激光测速仪在风扇转子流场测量中的一些基本方法，说明采用频谱分析法要比目前使用较多的计数法多普勒激光测速仪有更多的优越性，利用ＢＳＡＤＡＮＴＥＣ频谱测量系统实测了转子叶尖附近的流场，并获得了一些重要的流场信息。     

反流及径向流动相互作用下刚性旋翼气动特性研究

采用雷诺平均(Reynolds-averaged Navier-Stokes,RANS)方程针对直升机前飞不同前进比状态下四叶片刚性旋翼开展数值模拟研究,对比前进比0.1和0.6时的旋翼气动特性差异。计算结果表明,前飞时桨盘后行侧根部附近出现反流流动区域,翼型截面压强系数呈现非常规分布,该区域桨叶几乎不提供升力,且反流区面积随前进比的增大而增加。以静态前掠反流翼段为研究对象,采用脱体涡(Detached eddy simulation,DES)方法研究其非定常空气动力学特性,发现反流翼段表面出现特殊复杂的附着涡结构,在展向流动的影响下,翼段根部与尖部的涡结构发生耦合作用;反流翼段的升力系数随桨距角的增加而增大,且在失速迎角后并未下降。