皮托管式轴对称进气道跨音速内,外流场的AF—2方法计算

本文构造了轴向大扰动,径向小扰动轴对称速势方程的AF-2分解格式,计算了IC型皮托管式进气道的跨音速内外流场,并对这种分解格式进行了线化的稳定性分析。分析表明:AF-2方法的稳定性优于SLOR方法,其收敛速度至少比SLOR方法快10倍。     

旋成体轴对称跨声速全位势流的高效差分算法

本文用守恒型全位势方程,贴体坐标网格,对旋成体轴对称跨声速绕流的差分数值计算方法进行了研究;根据最佳收敛准则,提出了轴对称情形的AF2迭代算法,并将此算法应用于半球头柱体、弹体等各种外形的旋成体。与一般方法仅适用于亚声速自由流不同,本计算可从亚声速、跨声速自由流一直到低超声速自由流。计算结果表明,本文方法收敛快,与实验及其它方法的结果符合较好。     

CFD程序验证的虚构解方法及其边界精度匹配问题

采用网格收敛分析和虚构解理论对程序验证进行了理论探索与数值实验研究,并发展完善了该方法.给出了虚构解程序验证中一种比较简单的源项处理方案;给出了虚构解程序验证中气动方程组一种通用的虚构解系数选取方法,最后用虚构解方法对二维结构网格Euler流动计算程序ARC2D进行了验证,并发现了该程序边界条件实现精度与流场内部精度不相匹配问题,并给予了修正,数值实践证明虚构解方法是一种严格稳健有效的程序验证方法,用它完成的程序验证具有较高的可信度.     

亚音速机翼和机身大迎角下非线性气动载荷的两种计算方法

本文计算无粘亚音速流中任意平面形状薄翼和旋转体机身在大迎角下计入气流分离时的定常非线性气动载荷。采用两种方法:(1)定常流的松弛法;(2)陡然起动流的时间相依法。方法(1)对旋涡强度和自由涡的位置进行松弛迭代,取得了收敛的结果。方法(2)计算机翼陡然起动后一系列时刻的旋涡强度和自由涡的位置,直到结果趋近于常数解。亚音速流的压缩性效应是假设小扰动计算的。两种方法计算的结果很接近,和可以找到的实验数据也很接近。     

时间推进法求解叶轮机械内不可压流动

发展了一种用时间推进方法求解叶轮机械内部不可压流动的计算程序。在这个程序中, 使用分布体力方法简化湍流粘性的求解过程, 根据Ｂａｌｄｗ ｉｎ－Ｌｏａｍ ｘ 两层代数湍流模型确定粘性应力的大小。程序中还使用局部时间步长和多重网格数值技术加速收敛。火箭燃料主泵前诱导轮内部流场的实例计算表明, 采用这一计算程序能够正确预测叶轮机械内部不可压流场。     

出口宽高比对S形二元收敛喷管雷达散射截面的影响

在S形二元收敛喷管进出口面积、偏心距、面积变化规律及中心线变化规律不变的条件下,采用自适应超椭圆方法设计不同出口宽高比的S形二元收敛喷管.基于迭代物理光学法与等效边缘电磁流法自主开发计算腔体部件雷达散射截面（RCS）的程序;然后通过文献中的实验数据验证了计算程序的准确性和可靠性.通过计算程序分别分析了水平、垂直两种极化方式下不同S形二元收敛喷管的边缘绕射场与总散射场的电磁散射特性.结果表明:在水平、垂直两种极化方式下,喷管出口宽高比的变化对S形二元收敛喷管边缘绕射场的RCS影响较小,不同出口宽高比的S形二元收敛喷管边缘绕射场的RCS相差不超过4dB.喷管出口宽高比的变化对总散射场的RCS影响较大;正探测角时,宽高比为1.5时,在大部分探测角范围内总散射场有较低的RCS;负探测角时,宽高比为3.5时有较低的RCS.      

喷管型面对固体火箭发动机性能的影响

用时间相关法数值计算跨音速喷管流动时,其控制方程分别采用守恒型和非守恒型方程.通过对JPL喷管进行计算发现:在相同的精度下,非守恒型方程的收敛速度比守恒型方程的收敛速度快一倍多.采用修正阻尼法(DASHCO法)进行数值计算,其收敛速度均提高六倍以上.最后,采用直接优化法研究了喷管收敛段型面及喉部型面对发动机性能的影响.     

翼型跨音速小扰动势流-边界层计算

翼型绕流的外部无粘流,采用跨音速小扰动速势方程和改进小扰动速势方程对三种差分格式进行对比计算:(1)Murman-Cole非守恒式;(2)守恒式;(3)Engquist-Osher式。在相同的初值条件下,它们的计算稳定性、收敛性和收敛解基本相同,但激波位置(1)在前,(2)在后,(3)在中间。对人工时间阻尼项εφ_(xt)的作用进行了线化分析,所得结论与数值实验结果相符。二级精度格式计算表明:采用这种方案能在几乎不增加计算机时的前提下,大大提高差分计算精度。本文最后采用边界层积分方程法,计算翼型边界层,并与势流进行迭代计算。在弱激波条件下,迭代收敛解与实验值接近。     

球面收敛调节片喷管RCS数值模拟研究

结合物理光学迭代法(IPO)和等效棱边电磁流法(EEC),自主开发了电大尺寸腔体电磁散射特性计算程序,并通过对某试验模型雷达散射截面积(RCS)进行计算,完成了程序的可靠性验证。利用该程序对带有电磁信号遮挡板且采用3种不同尾缘修形技术的球面收敛调节片矢量喷管(SCFN)进行了RCS数值模拟研究。结果表明:随着入射方位角的增大,喷管的RCS整体上有减小的趋势;喷管尾缘的修形只有在雷达波以较大方位角入射时对减小RCS效果明显;电磁信号遮挡板的存在很大程度上减弱了喷管的电磁信号特征。     

轴对称喷管内流场的加权余数有限元计算

本文采用非定常技术,用加权余数有限元方法给出喷管内流场的数值计算方法,对收敛一扩散喷管的实例计算表明:本文计算的壁面和中心线M数分布曲线与实验结果是一致的。     

一种CFD/CSD耦合计算方法

针对柔性大展弦比机翼气动弹性分析和主动弹性机翼(AAW)设计发展了一种计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法。其主要思想是采用在同一物理时间弱耦合求解CFD/CSD技术。气动力采用非定常N-S方程的双时间有限体积求解技术,结构响应则采用有限元数值求解技术。CFD和CSD耦合计算的边界信息(气动力和网格)由所设计的界面程序传输。网格信息传输采用守恒体积转换(CVT)方法将CSD计算结构响应位移插值到CFD网格点上。变形已有的CFD网格技术用以确定CFD的变形网格。以位移或载荷的迭代误差为判断耦合计算的收敛标准。最后得到了机翼在Ma=0.8395,α=5.06°时CFD/CSD耦合计算的收敛值。针对计算结果分析了机翼受静气动弹性过程中结构响应和气动特性随时间变化的效应。初步研究结果表明:这种弱耦合方法求解非线性气动弹性问题是可行的。     

二元喷管三维粘性流场的数值模拟

本文用数值方法对二元喷管的三维粘性流场进行模拟,将 Mac Cormack有限差分格式拓广到全三维 N-S方程的求解。在数值迭代中使用加快收敛速度的当地时间步长,并根据格式压缩计算机存储量和计算量。进行了两个算例的数值模拟,计算结果与实验相符。      

不同喷口修形的二元收敛喷管RCS数值模拟

用等效棱边电磁流(EEC)方法和物理光学迭代(IPO)法自主开发了计算程序,用来分析不同尾缘修形喷管的雷达散射特性;在计算中考虑喷管出口边缘的绕射场对雷达散射截面(RCS)的影响;提出了一种高效的判断面元间遮挡关系的算法和喷管出口轮廓线的辨别方法;在此基础上,研究了四种不同喷口修形的二元收敛喷管的RCS特性.结果表明:该开发程序具有较好的计算精度;适当的喷口修形可降低重点姿态角下喷管的RCS.      

二维流热耦合数值模拟程序及方法研究

本文发展了一套二维流热耦合数值模拟程序,流体域中采用结构化网格,固体域中采用非结构化网格,流体域和固体域耦合面上的网格为非匹配连接。为了使该程序能够同时模拟低速不可压流动,程序中还加入了预处理方法。此外,本文还发展了一种能较好保持通量守恒的界面插值技术,来处理流体域和固体域耦合面上的非匹配网格。数值计算结果与实验值或商业软件结果的对比表明:本程序中的预处理方法有效地改善了低速不可压流动下的计算收敛性和收敛精度,使本程序能够同时模拟高速和低速不可压流动;扩压器和叶栅的流热耦合计算表明本文程序在流热耦合计算方面具有较好的模拟能力,同时也表明本文发展的界面插值技术的可靠性。     

亚超音速全机纵向气动系数数值计算方法的改善与发展

本报告是在亚超音速位流解基础上改进发展的数值计算方法。用面分布有限基本解计算作用在物面上的压力,用数值积分方法计算作用在物体上的力与力矩。本文采用联合流场、有效剖面吸力概念以及顺流向片条作可压缩紊流附面层的积分解,来改进阻力特性的数值计算。 数值方法已编制成计算程序。 本文提供了用这个计算程序计算的若干算例,并与实验结果作了对比。比较表明,计算结果与实验值有良好的一致性。     

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

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

气动热数值模拟中的网格相关性及收敛

针对气动热数值模拟中的网格相关性和收敛问题,采用计算流体力学(CFD)计算方法,以典型的钝头为算例,进行了数值模拟研究。研究结果表明:网格是气动热数值模拟中的关键因素,壁面附近法向网格间距最为敏感,热流结果随网格间距不同会出现数倍乃至数量级的差异,在得出正确的壁面压力结果的网格上不能保证得到正确的热流结果。分析认为壁面附近客观存在的巨大温度梯度是网格要求敏感的原因。计算表明:气动热计算的收敛过程比压力和流场收敛慢得多;在压力和流场完全收敛时,热流与最终的迭代收敛值偏差达20%以上;以方程残值的下降量级和压力收敛作为热流收敛的判别标准是不合适的。在气动热计算中,应直接观察热流数据的收敛,确保得到真正的收敛解。     

紊流射流掺混流场的有限差分计算

本文给出了二维紊流射流掺混流场的有限差分数值计算方法。用这种方法可以计算等截面和变截面二维(平面的和轴对称的)可压缩紊流射流掺混的流场分布。本方法已经用FORTRAN算法语言编制成计算程序,在FELIX C-256电子计算机上进行运算。通过算例和实验结果表明,本文提供的方法是可行的。     

GMRES算法在悬停旋翼数值模拟中的应用

发展了一种基于广义极小残差(GMRES)算法的悬停旋翼数值模拟方法,并对GMRES算法中矩阵与向量乘积的两种计算方法进行了分析和讨论。应用该方法在旋转坐标系中采用非结构混合网格和格点格式有限体积法对以绝对速度为变量的欧拉方程进行了直接求解,其中对流项的离散应用了基于Roe的Riemann近似解的迎风格式。对Caradonna-Tung旋翼跨声速悬停流场进行了数值模拟,计算结果与相关实验数据吻合较好,并与LU-SGS方法进行了对比,表明GMRES算法可以有效地加速流场的收敛,提高计算效率。     

关于跨音速零升力翼型绕流的松弛和两种自修正风洞的收敛性

本文采用横向小扰动而纵向大扰动速势方程,计算了跨音速零升力翼型的绕流。在线松弛的数值实验中,φ_γ的差分式用简单迭代和φ_(xx)的差分式用改进迭代时,稳定性较好。此结论与文献的线化理论分析相符。 本文用混合差分法数值模拟,证明了基于两个控制面上的静压和基于一个控制面及翼面上的静压的跨音速零升力翼型自修正风洞的收敛性。对前一种方案,NACA0012翼型,M∞=0.9,RAE104翼型,M∞=0.8,对后一种方案,NACA0012翼型,M∞=0.72,0.8,在迎角为零和风洞高度与翼弦之比为3时,均能收敛到无洞壁干扰的自由流。