跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

跨音速非等熵全位势方法及粘流/无粘流相互作用计算

跨音速流动是现代航空军机和民机设计中必须计算的一种流动。目前采用求解N-S方程的方法尚不能提供很经济而又非常精确的解,也就不能成为工程设计使用的常规手段。在实际飞行中,高雷诺数的事实表明,粘流/无粘流相互作用计算的方法可以成为工程设计计算的有力工具。 欧拉方程是描述无粘流的准确方程,但与全位势方程相比其求解过程要复杂得多,计算费用也高。为取得最优外形工程设计,要作大量的流动模拟计算,故推广全位势方法并提高其准确度是很有意义及具吸引力的。全位势方法假定流动是等熵和无旋的,故只适用于无激波或仅有较弱激波的流场计算。为改进传统全位势方法,应对其作两方面的修正。分析表明,旋度修正比熵修正更为高阶,作为改     

跨声速非等熵势函数方程的数值求解

基于动量方程的特点,本文将只与温度有关的“等熵密度”与熵的变化分离,而由散度型动量方程直接计算出气体经过激波的熵增。与经典势函数方程迭代,可以方便地求解非等熵跨声速流动。从对任意回转面叶栅的跨声速流动计算来看,这种非等熵势函数得到的激波位置可比等熵势函数前移1～2个网格,强度有所减弱,和试验结果较接近。由于各流线的熵不相等,对Kutta条件和出口边界条件作了修正。本文建议的方法也可用于计算流函数的熵增。     

跨声速非匀熵全速位方程的有限元数值解(英文)

本文改进了传统全速位方程的有限元数值解,既保留了全速位方程简单的优点又考虑了出现激波后流场不是匀熵的特点,对密度和后缘库塔条件进行修正,改善了计算结果。文中给出翼型和机翼的算例,其结果与传统全速位方程的结果相比,激波的位置偏往上游,强度较弱,这个趋势使结果较接近于欧拉方程的解。     

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

位势场与涡轮端区二次流的相互作用机理研究

采用数值模拟方法,研究了某亚声速涡轮内部下游导叶位势场与动叶端区二次流和叶尖泄漏流的非定常相互作用,详细分析了不同时刻的流动特性。结果表明:受下游导叶的非定常位势作用影响,动叶出口瞬时时刻二次涡结构呈现明显的周期性变化,并导致动叶出口熵增和效率均随时间改变,文中利用旋涡运动学理论对产生这一现象的流动机理进行了深入探讨。     

绕侧滑锥体,三维机翼与机身的超音速非线性全位势流计算

本文在Grossman的非守恒全位势方程的方法中,提出以固定网格为基础的头激波安装法,节省计算机时。计算结果与实验和欧拉方程解符合很好,与线化理论结果的比较,显示了线化理论的不足。     

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

宽速域翼尖涡及其与斜激波相互作用

针对宽速域飞行器产生的翼尖涡及其与斜激波相互作用问题,采用数值模拟方法,在来流马赫数0.2～6.0范围内,研究了翼尖涡的演化特性,并探讨了斜激波诱导翼尖涡破碎的机理。结果表明,在宽速域范围,翼尖涡的切向速度和环量分布仍然符合低速流动中的自相似关系式;在超声速至高超声速流动中,翼尖涡强度沿流向更快地减弱;当翼尖涡与斜激波相互作用时,涡核处压力越低的翼尖涡更容易破碎,而经典的涡破碎理论难以在宽速域范围内准确预测涡破碎现象;通过引入涡核处的压力亏损效应修正经典理论,大幅提升了预测斜激波诱导翼尖涡破碎的能力。     

旋翼VPM/CFD耦合模型的建立及其在小前进比气动分析中的应用

为了提高旋翼计算流体力学（CFD）流场计算效率及克服其在尾迹捕捉上的不足,将旋翼黏性涡粒子方法（VPM）与CFD分析相结合,建立了一个新的旋翼VPM/CFD耦合气动分析模型。在该模型中,采用VPM分析以实现对旋翼尾迹中黏性涡的高效捕捉而不引入数值耗散,而采用CFD分析用于精确地模拟旋翼桨叶近体区域内复杂的流动现象,同时也为VPM分析提供一个较高精度的涡源模型。至于两者信息交换,则使用集中涡源法将CFD信息传递至VPM分析,而VPM计算得到的尾迹信息则通过边界修正施加至CFD域的远场边界上,从而可以鲁棒地实现CFD域与VPM域的耦合计算。在此基础上,对“Helishape 7A旋翼”小前进比前飞桨-涡干扰（BVI）状态进行了较为深入的气动分析,计算结果表明：与全CFD计算比较,建立的VPM/CFD耦合分析模型可以有效地避免旋翼尾迹区桨尖涡的数值耗散,从而更加可靠地捕捉桨-涡干扰状态下的桨叶非定常气动载荷脉动,同时对于本文算例,计算效率可以提高30%以上。     

论定常不可压粘流的积分方程解法

本文对定常不可压粘性绕流积分解法中的涡量迭代过程、涡量边界条件的提法以及压强系数的计算等方面都作了改进,使数值解更易于稳定收敛并提高精度。为了验证本文的方法,对低雷诺数下的圆柱绕流进行了计算,结果表明计算效率是高的,迭代过程是稳定的。将本文方法推广到三维流动没有原则上的困难。     

多守恒性的差分格式的分析和比较

球面大气浅水波方程(在忽略摩擦和外源强迫的条件下)具有总能量、总质量、总位涡、总位涡拟能和总角动量守恒等五个重要的物理守恒性.然而,将大气浅水波方程差分离散化后,常常不能保持这些守恒性.本文对新发展的两种较好的差分格式进行比较;一种可称之为修正总能量守恒格式,它可以较好地保持上述五个守恒性中的四个;而另一种可称之为准辛格式,它可以使上述五种守恒性均近似保持.文中对这两类格式作了具体分析,并作相应的数值试验和比较,结果表明这两类格式均值得推广应用.     

Two-dimensional viscous vortex flow around a circular cylinder

A code based on a vortex method for simulating viscous vortex flow around a circular cylinder has been developed. The original method to obtain the no-through and no-slip condition introduced here is described in detail. The diffusive part of the vorticity transport equation is treated using a deterministic method based on the solution of the heat equation. Results have been obtained and are presented concerning the vorticity field, the velocity field, the evolution of drag and lift coefficients and the Strouhal number.     

光照处理中改进的Gamma矫正方法

在图像处理中，光照始终是影响图像质量的重要因素。为了有效弱化光照对图像质量的影响，利用非线性函数叠加的方法构造了一个合理的Gamma值变化曲线，使Gamma矫正方法充分适应了图像中光照变化的实际情况，并有效防止了Gamma矫正可能产生的图像失真。实验证明该方法在改善图像光照条件、提高图像质量方面具有良好的性能。     

弯曲流道内二次流动的运动学分析和 数值可视化方法

通过运动学分析,得到了涡量与二次流之间的关系方程,介绍了弯曲流道内二次流结构的可视化方法.以截面法向涡量为基础,可求解得到与涡量对应的二次流速度场.该方法具有以下几个特点:从运动学求解得到的涡量与二次流流函数之间的方程,直接体现了二次流作为涡的特性,且作为分析流动结构的后处理方法,能够量化地求解出二次流速度以及分布情况.该方法消除了位势流的影响,使二次流显示更加真实清晰,这种方法比速度矢量图法更加准确完善.      

一种基于三消隐点的摄像机自标定算法

提出了一种对摄像机进行自标定的方法.该方法是利用自然场景中三组两两正交的平行线所获得的消隐点之间具备的约束条件来实现标定的,并能够对获取的消隐点进行优化校正.利用该方法开展了实验研究,并与采用传统标定方法得到的结果进行了比较,结果表明本文所提出的算法具有较好的应用前景.     

小流量下周向弯曲叶片叶顶泄漏流动的PIV研究

对带有周向前弯和周向后弯叶片的低压轴流风机,采用PIV技术测量小流量工况下叶顶区域瞬态速度场。基于小流量压力峰值工况下的瞬态测量结果,详细阐述了周向弯曲叶片叶顶泄漏涡的演化过程。结果表明：小流量压力峰值工况下,周向弯曲叶片叶顶泄漏涡存在瞬态特性,泄漏涡起源向前缘方向迁移,涡核在向下游发展的过程中不断破碎,沿端壁向通道中部迁移,并发生径向迁移;周向后弯叶片的泄漏涡较周向前弯叶片迁移明显,涡核破碎过程剧烈;周向前弯叶片有利于叶顶泄漏涡的控制。