用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。     

亚、超音速小侧滑角机翼升力面理论和数值计算方法

本文推导出侧缘与自由流不平行的常值面涡半无限角翼基本解,其特征是从面涡小块背风侧缘拖出的涡不是顺着自由流,而是沿着侧缘拖出,从而克服了用Woodward方法计算侧滑翼绕流中使用对称分划小块的困难,提出了适用于从亚音速直到超音速范围的推广的Woodward方法。本方法将问题分解为若干个对称与反对称问题,大大地节省了计算机内存量或计算时间。本文还从基本方程,或者直接从数值计算方法出发,推导出准相似律。计算结果与精确解以及和其他方法的结果符合得很好。     

基于Euler方程的翼型优化设计高效基因算法研究

把基因算法 ( GAs)用于翼型优化设计 ,构造了用于优化设计的适应度函数 ,发展了一种具有全局优化性能的翼型优化设计方法。为提高算法的搜索效率 ,提出了改进型的实数编码方式 ,设计了新型的 Bézier样条翼型基因表达方式 ,并发展了高效隐式 Euler方程解算器。用此方法进行了基于升阻比极大化的跨音速和基于阻力极小化的超音速翼型优化设计数值模拟 ,取得了令人满意的结果。     

轴对称进气道外罩跨音速非守恒型全位势流AF3算法

本文针对轴对称进气道外罩跨音速绕流问题的特点,采用了贴体坐标网格;并对这类网格拓扑构造了三种Baker的隐式近似因式分解AF3有限差分迭代算法。通过理论分析与数值实验,找到了既能计算亚临界情形,又能计算超临界情形的AF3格式;并就此格式讨论了适合于进气道绕流特性的AF3迭代格式中间变量的提法,具体提出了格式实施的相应方法,研制了可供分析任意轴对称进气道外罩跨音速特性的计算程序。由于采用了任意曲线坐标系上的轴对称非守恒型全位势方程适合于用不同方法生成的计算网格,因此,本文的工作还有可能进一步推广到计算更复杂的轴对称进气道跨音速绕流问题。计算表明,本方法收敛快,结果好。     

翼-身组合体跨音速守恒型全位势流AF2高效算法的应用

 本文用Joukowski交换将圆截面机身的每一横截面逐一转绘为垂直切缝,因而将翼-身组合体转绘为某单独机翼。我们用作者采用过的数值保角转绘法对转换的单独机翼生成O型贴体坐标网格,然后用逆转绘构造翼-身组合体绕流的O型贴体坐标网格。跨音速计算采用守恒型金位势方程,精确边界条件和AF2高效有限差分迭代算法。本文指出,只需少量修改,就可将单独翼分析程序扩展为可同时适用于翼-身组合体。本文考虑的组合体其机身为有限或无限长,圆形截面,机翼前后机身轴线允许弯曲;机翼有任意平面形状,但翼尖弦长不能太短。     

跨音速翼型非守恒全位势流AF3高效算法的应用

本文利用Baker的AF3高效迭代算法,改进了Jameson提出的计算跨音速翼型非守恒全位势流的有限差分方法,研制了任意翼型跨音速分析的高效算法程序。许多试验算例表明,本算法,收敛速度快,且计算结果与Jameson的原方法非常一致。     

匹配数据的平滑滤波技术

讨论了匹配测量数据的平滑，滤波和外推问题，所谓匹配测量数据是指位置参数和速度参数呈求导的匹配关系，提出了利用数据的匹配特性这一先验信息来改善多项式最小二乘滤波器的设计，对测量数据中的位置和测速数据进行统一处理，在最小方差意义下给出了平滑滤波公式及其效率方差比，由于充分利用了先验信息，减少了待估计参数的个数，显著提高了参数的估计精度，仿真计算结果表明该方法能显著改善目标运动参数的估计精度。     

用全位势方程计算机翼的亚声速,跨声速和超声速绕流

对大后掠小展弦比细长机翼,本文对机翼纵轴垂直的每一横流截面生成O型网格,形成对机翼流场的H-O型网格,用守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕机翼的可压缩位流。自由流可从亚声速直到低超声速的全部跨声速范围。本算法要求机翼前缘有大后掠角,后缘可稍许后掠或前掠。本文算例表明,所研制的计算程序已可提供工程实用。     

旋成体大迎角分离流的Euler方程数值模拟

用三维Euler方程计算旋成体的大迎角绕流,我们发现:亚声速自由流时,计算得到的是无分离流动,与实际的分离流不相符合,因而计算的法向力大大小于有分离的实验值;而当自由流为超声速,且横向自由流M数较高时,由于背风表面出现激波,无粘Euler方程的数值计算能自动捕捉到背风面上的分离涡,计算法向力与实验值接近。为了对亚声速自由流用Euler方程数值模拟大迎角旋成体背风面上分离涡,根据问题的物理图像,在假设的分离线上,提出强加Kutta条件的方法,成功地计算出有分离流动,使计算的法向力与实验值接近。