涡方法分析并列圆柱的旋涡脱落现象

利用快速离散涡方法及流场可视化技术分析二维粘性不可压缩均匀流中并列圆柱体的旋涡脱落现象，识别作用在圆柱上的气动力及旋涡脱落频率。使用边界元方法及引入一个新颖的涡量边界条件计算表面涡量，利用系数影响矩阵技术获得无网格方法，采用自适应快速多极子算法快速计算涡元速度。数值算例再现了不同间距并列圆柱体的同步同相旋涡脱落、同步反相旋涡脱落、尾流的不对称性及间隙流的半稳定性等现象，数值结果与早期实验结果吻合得很好。     

升力面尾涡卷起的地面效应

本文就二维平板从理论上证明了改进涡格法的超收敛性。用该方法研究了近地升力面尾涡卷起的地面效应。尾涡面上的自由线涡方向应与当地流线方向一致,通过迭代满足这一非线性边界条件的确定尾涡的形状。地面效应以镜象涡系模拟。一系列计算结果表明迭代迅速收敛。本文讨论了升力面距地高度对尾涡位置和形状的影响,计算表明当该高度小于一定距离后,尾涡呈现不稳定。     

对流场自适应的确定性涡方法

本文系统地研究了随机走步涡团法所面临的主要问题，以统一的途径构造了自适应涡团模型，粘性扩散的确定性算法及确定边界涡量的新方法，数值结果是令人满意的。     

关于旋涡定义的思考

针对旋涡定义这个长期未解决的难题,在分析历史上提出的几种思路的基础上,本文提出:定义旋涡的总目标,应当体现其管状运动形态,并与其作为流体运动肌腱的动力学功能有机结合;该定义应能以统一的方式覆盖经典涡动力学业已成熟的结果,并且指引对湍流中复杂涡状结构及其相互作用的识别。据此,本文根据涡量场演化的已知规律,提出了旋涡的一个动力学定义,作为继续深化讨论的参考。同时提出,人们针对湍流结构研究构造的各种涡判据,虽然对湍流涡状结构的可视化起了重要作用,但因其难以遵循管式涡量场的因果演化而无法取代完备的涡定义。相比之下,近十年来形成的涡量面理论,有望在未来发展中使旋涡的定义和检验问题回归涡量动力学的自然规律。     

燃气轮机压气机涡量动力学理论及分析方法

阐述了压气机涡量动力学理论及分析方法,建立了压气机总性能参数与两个重要的涡量参数(边界涡量流(BVF)和周向涡量)的直接数学物理关系,研究了基于BVF和周向涡量的分析优化方法.将该方法应用于燃气轮机多级轴流压气机中进行涡量动力学分析,从涡量角度指出了改进的方向.结果表明:周向涡量能反映近壁面的高损失区,使周向涡量峰值束缚在近壁区有利于降低端区损失,在通流设计中可通过优化环量分布控制周向涡量分布,算例中基于周向涡量优化可使跨声压气机转子效率提高1.13%;边界涡量流BVF能反映旋涡的壁面根源,通过优化BVF的分布可控制涡量壁面根源,有利于抑制旋涡和流动分离,基于BVF优化可使转子效率提高1.12%.      

可压缩湍流边界层近壁区马蹄涡的演化

采用大涡模拟方法对可压缩槽道湍流进行了数值模拟.通过给定扰动波方法确定进口边界条件,很快得到了完全发展的湍流.利用动态亚格子模型、预处理技术以及高阶格式和方法,得到了精确的数值结果.计算得到的速度剖面和脉动分量与直接数值模拟的解进行了对比,验证了计算结果的可靠性.结果清楚地展示了湍流场中马蹄涡的生成和演化过程,包括对称涡腿、单侧涡腿的马蹄涡结构,以及流场中存在由亚谐波引起的拟序结构的交错现象.分析了湍流边界层近壁区可能出现的复杂的马蹄涡结构,及其涡结构形状也不严格对称的现象.     

尾涡中涡强与轴向速度的相互关系

大型飞机的翼尖涡会对飞入它尾流中的其它飞机造成巨大的危害。文章希望找到尾涡核中的轴向速度超过自由流速度的条件，从而深入了解这种尾涡的特性。在风洞实验中，我们采用的机翼模型的型号是NACA0015，它的半翼展展弦比为0．80。我们利用三传感器热线探针来测量尾涡核中的轴向速度，并将之与测得的涡的环量进行了比较，结果发现，轴向速度与一个无量纲环量参数呈线性关系。当这一环量参数较小时，轴向速度小于自由流速度，而当环量增大时，轴向速度也将逐渐增大，甚至能够超出自由流速度70％。我们还发现，轴向速度的大小对实验中所采用的两种端帽结构——平型和圆型十分敏感。另外，如果利用涡的半径来修正这一环量参数，那么两种端帽结构下的数据将会重合为一条直线。     

边条翼前缘涡非定常涡场特性研究中PIV技术的应用

描述了应用PIV技术在水槽中对边条机翼上旋涡及破裂旋涡流场进行的测量和分析。实验是在北航水槽中进行的。通过PIV技术的测量，揭示了旋涡及破裂旋涡中的非定常特性，这种非定常特性同飞机上机翼、尾翼的抖振密切相关。实验结果表明，对于未破裂的边条涡，存在着两种非定常特性，其一是剪切层中不断地有小涡沿剪切层输运和合并。其二是由一次涡诱导的二次涡与剪切层中的小涡互相诱导引起的非定常现象。对于破裂涡，则发现与未破裂的涡相比，截面上涡量分布的区域突然扩大很多，最大涡量的绝对值也比上游未破裂区截面上的涡量最大值小。此外还发现在涡量分布区域出现反涡量，这同涡破裂后出现涡核螺旋变形有关。对于同一截面处涡量分布是非定常的。     

侧风条件下地面涡来源和气动特性的数值研究

为了研究侧风条件下地面涡的涡量源和地面涡的气动特性，选取缩比进气道并对侧风条件下的流场进行了数值模拟。结果表明:90°侧风条件下，尾涡和环境涡量都是地面涡的重要来源，即使不存在环境涡量，尾涡仍然可以单独形成地面涡，尾涡-地面涡流动模型能够很好地解释这一现象。地面涡涡量随来流速度的增加先增大后减小，临界速度比随离地间隙的增大而增大。离地间隙增大一倍，临界速度比增大65%，地面涡涡量峰值减小30%。受到地面涡本身以及流动分离的影响，进气道的总压损失随着来流速度比的减小而增大。      

关于涡方法的讨论

本文讨论了如何提高二维涡方法的计算精度和减小涡方法计算量的问题。文中指出了改进的途径,同时列举了当前这些方面的研究进展情况。     

高超声速边界层的转捩及预测

首先描述了边界层转捩的基本过程及研究内容。在此基础上,指出了高超声速边界层不同于不可压缩边界层的流动不稳定性特性,并介绍了边界层的转捩机理与感受性特征;给出了高超声速三维边界层中预测转捩的常用方法,并着重介绍了多用于工程实际的e N方法以及对e N方法的理性改进,同时列举了在高超声速三维边界层中应用e N方法实现转捩预测的多个实例。最后,分析并总结了高超声速边界层转捩预测所存在的困难及需要解决的问题。     

时间依赖边界条件在亚声速自由射流数值模拟中的应用

对亚声速自由射流进行数值模拟时，无法预先给定远场边界和无反射入流边界上的速度值。以往采用的黎曼不变量在自由边界上给定无反射条件的处理方法，只给出了这些边界法线方向上的一维特征关系。时间依赖边界条件是一种高精度的边界处理方法，可以很方便地处理计算中遇到的各种边界，广泛地应用于DNS中。我们将这种边界条件用于定常流动的计算中，计算了出口马赫数为0．4和0．7的轴对称自由射流，获得了与实验结果吻合的数值模拟结果。     

Marangoni boundary layers

This paper deals with the theory of Marangoni boundary layers developing near the free surface of a nonuniformly heated liquid.

In the Introduction the notion of a thermocapillary effect is considered, and the conditions, under which thermocapillary flow is dominant, are elucidated. The similarity criteria of the problems are revealed, and the effects resulting from melt convection are considered.

In the second section the mathematical aspects of the theory of Marangoni boundary layers are considered. The boundary layer effect on the form of a free surface is investigated, and some examples are presented.

In Section 3 a qualitative flow pattern in a Marangoni boundary layer is analyzed, similarities and differences in the problems of Prandtl and Marangoni boundary layer theory are discussed. The possibility of Marangoni boundary layer degeneration at the outer boundary, as well as the initiation of flow containing a counter-stream zone are emphasized.

The fourth section outlines the calculation scheme on the basis of boundary-layer theory for convective flows with directed crystallization and in the absence of mass forces. Some examples are presented.

In the fifth section some problems of the theory of Marangoni boundary layers that need further investigation are discussed.     

碳／碳材料烧蚀对电离边界层的影响

本文给出了碳/碳复合材料热化学烧蚀对层流和湍流电离边界层影响的计算方法。边界层中考虑了20个气体组元,并假设所有化学反应平衡。应用了边界层方程的新解法和简化的输运性质处理。算例表明,烧蚀对边界层的电离特性有很大影响,烧蚀产物中,碱金属的电离并不总占优势。     

超临界层流机翼边界层及气动特性分析

高空长航时无人机设计巡航状态的雷诺数较小,黏性边界层对气动特性的影响较大。详细分析了雷诺数对机翼边界层和气动力的影响,用数值方法对超临界层流机翼三维层流-转捩-湍流混合边界层特性进行了研究,分析比较了高空小雷诺数和中空大雷诺数情况下机翼三维边界层的特性,尤其是边界层转捩点位置、表面摩阻和气动特性的雷诺数效应。研究表明雷诺数对于高空无人机机翼边界层厚度、摩擦阻力和升阻比影响较大;对层流机翼的转捩点位置和升力系数影响较小;自然层流机翼技术可以应用于高空无人机设计。     

径流/斜流叶轮全三维反问题方法中的边界条件和级数截断

研究了脉动流场不同类型的进出口边界条件,比较了端壁边界条件的不同给法,并用基于定常雷诺平均Navier-Stokes(RANS)的商用计算流体动力学(CFD)软件数值研究了进出口、端壁的边界条件和级数截断误差对气动设计的影响.研究表明:进出口边界条件的影响可忽略不计,端壁边界条件的影响较大.取傅里叶级数约前15项后,截断误差的影响可忽略不计.      

实现边界元并行分布计算的一种新方法

 给出了实现边界元并行分布计算的一种新方法。它将边界元离散节点分为与计算网络所拥有的处理机台数相对应的若干组, 以此实现任务分配, 独立地计算形成每组节点对应的离散边界元积分方程, 最后利用迭代方法行求解分布的边界元方程组。以弹性问题为例, 在Inmos T800 Transputer 系统上实现了边界元计算。验证了本文方法的可行性与有效性。     

高超声速锥柱裙模型边界层转捩的弹道靶实验

为研究高超声速边界层转捩现象、给边界层计算提供可靠的对比数据,在中国空气动力研究与发展中心趟高速弹道靶上开展了锥柱裙模型高超声速边界层转捩的自由飞实验。所采用的锥柱裙模型全长105mm,飞行速度1．94km／s（Ma=5．65）,单位雷诺数4．32×10^7～1．20×10^8m-1。使用激光阴影成像技术,获得了锥柱裙模型边界层转捩和湍流边界层发展的图像,测得的湍流边界层厚度在0．6～2．2mm之间,湍流涡的流向尺寸与边界层厚度的比值介于0．3～0．8之间且沿流向呈下降趋势。实验结果表明：弹道靶实验能够获得给定飞行环境下的高超声速边界层转捩图像,从图像中可以清晰判断转捩位置或区域、测量边界层厚度和分析湍流涡的尺寸。     

小攻角高超声速钝锥边界层内扰动演化的数值模拟

首先从有限差分格式出发,给出了基本无振荡的高阶激波捕捉格式,然后,采用数值模拟方法对马赫数为6的2°攻角高超声速钝锥边界层的稳定性进行了研究。计算发现,由于攻角的存在,钝锥的稳定性特征与零攻角时有本质的差别,比如背风面的扰动比迎风面增长更快,但扰动增长最慢的地方并不是迎风面,而是侧面的某个位置;又比如背风面主要是长波起作用,迎风面和侧面主要是短波起作用;斜模式不稳定在整个钝锥边界层中起最主要作用。     

Transverse curvature effects in turbulent boundary layer

The effect of transverse surface curvature on the turbulent boundary layer is reviewed by recourse to experiments on axial flow along a circular cylinder. Three flow regimes are identified depending on values of the two controlling parameters, namely, the Reynolds number and the ratio of the boundary layer thickness to cylinder radius. The boundary layer flow resembles a wake when both parameters are large. As expected, the effect of curvature is small when the Reynolds number is large and the boundary layer is thin. When the boundary layer is thick and the Reynolds number is small, which is typical of laboratory investigations, the effect of transverse curvature is felt throughout the boundary layer with evidence for relaminarization at the low Reynolds numbers. This review describes the experimental evidence and points out gaps that remain.