用隐式近似因子分解法计算喷管跨音速流场

本文用隐式近似因子分解法结合贴体曲线坐标计算了定常、无粘、跨音速喷管流场,选择了参数变化比较剧烈的小喉部曲率半径喷管,陡壁收敛段喷管和潜入喷管三种算例.计算结果表明,该方法收敛快,并且具有良好的精度.Courant数可以取至7.这种隐式法与显式法相比较,花费的机时少得多.本文的方法可望推广到求解N-S方程、PNS方程,和二相有粘喷管流动计算.     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

单级跨音速压气机三维非定常粘性流动计算

本文以数值求解雷诺平均非定常N-S方程方法完成了单级跨音速压气机动静叶相干三维非定常粘性流动计算,给出了高分辨率NND格式的时间推进有限差分解。对某跨音速单级压气机进行了定常和非定常流的详细数值计算,设计转速下定常粘性流场计算所得到的性能参数与实验结果吻合较好。在此基础上,对最高效率状态进行了动/静叶相互干扰的非定常流动计算,并对结果进行了分析。      

跨音速叶栅流场的无粘-边界层迭代计算

二维跨音速叶栅流场的计算是在无粘时间相关法计算程序的基础上采用无粘-边界层迭代而得。马赫数从亚音速到超音速范围的几种数值试验均与测试结果吻合一致。在边界层无强烈脱流时可获得很好的结果。     

机翼跨音速无粘绕流数值计算

本文介绍机翼跨音速定常无粘绕流的一种数值计算方法。选用精确速势方程作为问题的数学模型。在直角坐标系中,经过适当的坐标变换,先把后掠翼变成矩形翼,再把无限的物理空间变成有限的计算空间。在计算空间中,用有限差分格式(亚音速区用中心差分,超音速区用旋转差分)对精确速势方程离散化。差分方程形成的代数方程组用线松弛迭代求解。     

二元喷管跨音速流场的数值计算

本文采用有限体积法进行空间离散,Runge-Kutta法进行时间推进的办法,通过求解Euler方程来模拟二元收敛—扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡,在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解。计算结果与试验数据和其他数值方法吻合良好。      

亚跨超音速喷管内流场的整体矢通量分裂算法

从N－S方程出发，采用先进的矢通量分裂算法，对亚跨超音速喷管内流场进行了整体模拟计算。文中对几种不同计算格式的计算效率进行了比较，给出了超音速占优混合流动的直接分割求解方法。计算采用代数湍流模型，跨音速段的计算结果同其它实验和计算结果进行了比较。计算中发现对出口边界条件进行亚音速、超音速和有回流分别处理很有必要；边界层内的M数沿壁面法向变化非常快；不同壁温条件和不同的流动Re数均会改变边界层的厚度。     

火箭发动机塞式喷管流场的数值研究

用时间相关法和张涵信的无振荡、无自由参数ＮＮＤ显式格式,对火箭发动机塞式喷管二维轴对称定常流场进行了数值模拟。计算得到了有粘和无粘中两种情况下,喷管流场的马赫数压强、温度、流线以及速度矢量分布图,计算网络采用代数方法生成,并通过求解Ｌａｐａｃｅ方程对网格进行修正。计算结果表明：ＮＮＤ格式用于计算塞式喷管流场具有较高的效率和良好的分辨率;无粘和有粘情况得到的流场结构有较大的差别;基础喷管的出口马赫数     

固体火箭喷管两相粘性跨音速流场计算

本文进行了固体火箭喷管两相粘性湍流跨音速流场计算.粘性的气相控制方程用隐式近似因子分解法求解,粒子方程采用跟踪粒子轨迹的特征线法求解,粘性湍流选用代数模型,气相和凝相充分地偶合.CFL数可以取500左右,收敛速度快,使粘性两相跨音速喷管流场计算耗费的机时达到工程计算可以接受的程度.通过计算,获得了在粘性和粒子同时作用下的流场参数分布.这对固体火箭发动机流场需要同时考虑粒子和粘性作用的专题研究大有裨益.     

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

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

飞机／进气道系统内外流场计算

为了对飞机/进气道系统内外流场进行数值模拟,本文提出并研制了计算定常(或非定常)外压式飞机/进气道系统内外流场的计算方法和程序。它适用于任意几何外形的飞机/进气道在亚音速、跨音速或超音速来流条件下的内外流场和气动参数计算。求解的主管方程是在一般曲线坐标下的 Euler 方程、简化/雷诺平均的 Navier-Stokes方程。用隐式时间差分和近似因式分解方法离散和求解.针对飞机/进气道的复杂流场,本文采用了贴体网格生成技术和分区域方法.对于湍流流动,把 Baldwin-Lomax 湍流模型推广到飞机/进气道内外流场计算并使用了壁面律计算公式。应用上述方法对绕翼型的亚、跨音速、绕可压缩拐角的超音速湍流以及飞机进气道系统内外流场进行了计算,计算结果与实验数据和其它计算结果符合较好。     

二元喷管跨音速流场的数值计算

本采用有限体积法进行空间离散，Runge-Kutta法进行时间推进的办法，通过求解Euler方程来模拟二元收敛－扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡，在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解，计算结果与试验数据和其他数值方法吻合良好。     

跨音速翼型分离绕流的粘性/无粘干扰迭代计算

本文用有粘/无粘干扰迭代的概念计算了跨音速任意翼型的绕流问题。位流的速位方程用AF2格式求解,而边界层微分方程用C-S盒式法求解,逆算法的引用可以克服边界层方程在分离点处的奇性问题,对分离区湍流代数模型的修正可以得到与实验更吻合的结果。计算结果表明有粘/无粘干扰迭代概念在小分离泡的情况下也是适用的。     

跨音速翼型多参数多约束优化设计方法

本文通过采用EXTREM数值优化方法并基于求解全速势方程的跨音速粘流翼型计算方法, 研究发展了一种应用形状函数修形的跨音速翼型多参数多约束数值优化设计方法, 应用该方法可以从普通低速翼型和超临界翼型出发通过数值优化后得到在跨音速区阻力系数最小化的翼型几何外形。设计表明, 该方法具有收敛快、调用目标函数次数少等优点。      

翼型跨音速绕流的数值计算

本文选用精确速势方程作为翼型跨音速无粘绕流的数学模型。在变换过的直角座标中,用有限差分(亚音速区用中心差分,超音速区用旋转差分)对精确速势方程离散化。差分方程形成的代数方程组用列松弛迭代法求解。 为了计入粘性效应,利用精确速势方程和附面层动量积分方程进行联合迭代求解。这一点对超临界翼型的计算显得特别重要。 算例与其它数值解及实验作了比较,对于一般翼型和超临界翼型,本文结果是良好的。     

包括粘性影响的跨音速纵向大扰动非定常翼型绕流计算

本文提出一种计算效率高、并改进小扰动理论的二维跨音速定常和非定常流的计算方法——非定常纵向大扰动流速势方程和边界条件的数值解。本方法还考虑了包括边界层位移厚度以及激波-边界层干扰的粘性影响。文中给出了NACA 0012翼型和NLR 7301超临界翼型绕流的算例,计算结果与实验作了比较。     

跨音速低展弦比轴流风扇转子内部流场的数值研究

采用三维雷诺平均Ｎ－Ｓ方程程序求解跨音速压气机风扇转子内部流场。这个程序采用有限体积显式时间推进方法求解控制方程,方程自变量选取在控制体顶点,采用简单的Ｈ型网格离散方程,离散后的代数方程用两步Ｒｕｎｇｅ－Ｋｕｔａ法进行时间推进求解,使用多重网格和当地时间步长方法加速计算收敛。程序中使用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型模拟湍流影响。用该方法计算了ＮＡＳＡ６７＃跨音速风扇转子的性能曲线,并重点计算最高效率点工况下内部流场。计算结果与实验结果的比较,证明该法可以用来模拟跨音速风扇转子内部流动。     

跨音速颤振计算方法研究

采用快捷的动态弹性变形技术生成三维块结构动态贴体网格。结构运动方程采用二级精度的龙格—库塔时间推进。气动力求解则用非定常欧拉方程的双时间有限体积推进,外时间为物理时间,与结构运动方程同步;对每个真实物理时间步采用五步Runge Kutta时间推进进行虚拟定常迭代。针对可压流颤振计算中存在的质量不相似问题,采用变质量、变刚度的方法计算出质量匹配点处的颤振速压值,并根据颤振速压随质量或刚度倍数的变化趋势,得出可用的飞机跨音速颤振速压。对一飞机简化外形的跨音速颤振特性计算得到了合理的计算结果。     

低超声速来流中机翼跨声速非定常气动力数值解法

本文给出一种低超声速来流中机翼跨声速非定常流的计算方法。该方法使用跨声速小扰动理论,正确地计入了沿流向平面流动特征线的压力波的运动边界条件。流动方程是使用时间精确交替方向隐式(ADI)有限差分解法求解。试算了M=1.1的矩形翼和M=1.1、1.34的F5机翼的定常和非定常气动力,并与国外的风洞试验和计算结果作了比较,符合很好。表明本方法是有效的。     

计算翼型跨音速有势绕流的压强极小积分有限元法

本文采用压强极小积分有限元法计算了NACA 0012翼型的跨音速流场。出发方程为全速位方程。用人工密度的方式引入人工粘性。计算网格是通过解析/数值变换的方法自动生成的。采用线松弛法和人工时间相关法求解有限元方程。给出了有关的数值结果。