跨声速机翼的设计计算方法

本文介绍一种三维跨声速机翼设计计算方法,它以跨声速小扰动理论为基础,以给定的压力分布与分析计算的压力分布之差作为目标函数,应用Green定理建立起目标函数与机翼剖面形状修改量的线性代数方程组,从而求解出翼剖面的修改量,加到原来机翼各剖面上,再用分析计算程序计算,经过多次的反复,就可以达到设计目标。本方法的特点是分析计算程序与设计程序相对独立,只作小量修改就可联合使用。本文介绍用此方法设计的两个算例,说明设计结果是好的,计算机时与其它方法相比是大大缩短了。     

跨声速欧拉方程的一种隐式有限元解法

本文用有限体积Ｇａｌｅｒｋｉｎ法求解了跨声速二维欧拉方程。求通量运用了Ｏｓｈｅｒ格式，通过外插把原来只有一阶精度的逆风格式改进成二阶精度。并采用基于一性化流量的隐式格式和非结构化网格，使程序具有较高的收敛效率。     

一种跨声速翼型设计方法及设计诸例

本文采用正反迭代、余量修正的设计思想,发展了一种能迅速达到指定目标压力分布的跨声速翼型设计方法。提出了推导各积分的解析表达式取代数值积分、加权光顺过程、直接含粘性迭代和进行矩阵状态分析等有效的改进措施,克服了现有设计方法中的主要缺陷。在设计软件中,正计算程序作为一个“插件”与反设计部分耦合,因此可以方便地与新的正计算方法组合以提高整个设计系统的先进性。对多种超临界、无激波、自然层流压力分布及低速翼型的设计实例表明,本方法只需少数几次设计迭代即可以设计出准确收敛于目标的新翼型。     

平面大头叶栅跨声速流动的一种有效解法

准确数值预估叶片头部表面的压力分布,一直是叶片气膜冷却设计中十分重要而又尚未圆满解决的问题。本文针对这个问题提出了一个高准确度、高效率的计算方法。 在文献[1]的基础上,本文从任意正交曲线坐标系中的基本方程出发,引入Von Misses变换,导出了相应的流线控制方程,并具体地在拟边界层坐标系中计算了RA叶栅流场,结果与实验符合甚好,据作者所知,本文方法是已见到的计算同类叶栅方法中准确度最高、花费机时最少的。     

联合使用激波捕获—拟合法求解激波间断的跨声速流场

吸取激波捕获和激波拟合两种方法各自的优点,发展了一种求解有激波的跨声速流场的方法,并对回转叶栅中定常跨声速流场做了计算。计算中,通过捕获法确定初始激波位置,然后经反复修正,拟合出确定的流场通道激波。计算表明,这种方法可自动获得清晰的激波,流场中气流各参数分布合理。这一方法可适用于复杂边界和不同进口M数,计算时间仅比势函数方法多一倍左右。     

FL—23风洞洞壁对跨声速颤振的干扰影响

本文提供了在FL-23风洞中进行的试验研究结果。试验采用动力相似的四种不同尺寸的60°三角机翼平板半模型。结果表明:洞壁干扰提高了跨声速颤振临界速压,当模型展长对试验段宽度的比值Lm/B为0.687、0.600和0.512时,使颤振临界速度的压缩性修正系数分别减小3.8％、2.5％和1.4％。改变开孔率的试验定性验证了上述结果的合理性。     

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

跨声速翼型多目标优化设计方法

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

低雷诺数翼型和跨声速翼型边界层计算

本文利用边界层动量积分方程和平均流动能积分方程兼容计算了翼面的层流和湍流边界层流动。文中采用一种e~9型转捩判别公式预测翼面存在层流分离气泡情况的转捩位置。并引入对剪应力系数C_r的滞后方程以体现湍流中雷诺应力的历程效应。为避免在翼面流动分离时边界层方程的奇异,采用反解法计算。算例表明,计算与测量结果吻合良好。     

机翼跨声速绕流的Galerkin有限元解法

本文给出了翼型和机翼跨声速绕流的Galerkin带权余数有限元解法。利用一系列适当形状的元素将计算区域离散,在元素内采用线性插值函数。对于超临界流动,文中采用了适合于有限元法的超声速上风技术,具有捕捉激波的能力。将这种技术应用于线松弛迭代解法中,在计算钝头翼型和机翼绕流时,都获得了成功。     

跨声速三维Euler方程的显式／隐式格式有限元解法

本文采用有限体积Ｇａｌｅｒｋｉｎ法和非结构网格数值求解了跨声速三维Ｅｕｌｅｒ方程,通量的求解采用ＶａｎＬｅｅｒ－阶逆风矢能量发裂方式,并通过外插使其上升为二阶精度,分别采用了显式及隐式两种格式求解。     

高速风洞洞壁干扰减小与修正研究

亚跨声速风洞试验的洞壁干扰问题是影响风洞试验结果准确度的一个重要因素。南京航空学院ＮＨ－１高速风洞首先使用了两种可变开闭比风洞壁来减小洞壁干扰，然后发展了一系列用壁压信息法对剩余洞壁干扰效应进行修正的方法，对国内外大量高速风洞实验数据进行了洞壁干扰修正，将修正结果和ＮＡＳＡ的非线性修正法结果。自修正风洞实验结果及无洞壁干扰的Ｎ－Ｓ方程计算结果进行了比较。结果表明，本文的修正方法结构正确，而计算量远     

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

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

数值分析跨声速主流与边界层迭代

本文首次提出新的用于分析跨声速无粘主流与边界层迭代过程稳定性的物理模型和数学模型,在局部冻结系数的意义下从理论上细致分析了主流采用跨声速流函数法,而边界层分别采用正模式、逆模式、半逆模式计算时主流与边界层迭代过程的稳定性和放大因子。指出上述这些迭代并不总收敛,并给出了使这一过程达到收敛应具备的条件,这对指导电子计算机上的数值试验有特别重要的意义。文中完成了4个典型算例,得到了与实验较为吻合的数值结果。     

跨声速翼型激波／边界层干扰被动控制的数值模拟计算

本文应用数值计算的方法着重研究了跨声速翼型开孔壁的减阻效果。计算采用了边界层与位流相互作用的模型,藉以了解开孔壁对激波强度及结构的影响和对边界层控制的效果。通过对NACA0012翼型的计算表明,本文采用的自然吹吸的开孔模型能够显著地削弱激波的强度,改善激波的结构,但会使粘性损失增加。在马赫数较小时,翼型开孔后总阻力会增加,而在大马赫数时,开孔翼型的减阻效果才表现出来。这种趋势是与实验结果相吻合的。     

喷管跨声速流场计算

本文用时间相关法求解定常跨声速喷管流场。对守恒型和非守恒型方程采用全场统一时间步长、当地时间步长和修正超阻尼法进行了计算。计算表明:采用当地时间步长可大大提高收敛速度,采用全场统一时间步长时,数值解非守恒型方程收敛较快;采用修正超阻尼法也能提高收敛速度,但与阻尼系数和时间间隔有关。