一种求解涡轮平面叶栅的跨音松驰法

一、引　　言跨音松驰法是七十年代初发展的一种新方法。文献 [1 ]是一篇关于跨音松驰法的奠基性文章 ,文章提出了求解小扰动势方程的混合差分格式。文献 [2 ]发展了求解全位势方程的旋转混合差分格式。近年来 ,求解跨音全位势方程的方法有了更全面的发展 [4、5]。作者于 80年发展了一种求解大弯度二维叶栅的跨音松驰法[3] 。本文将该法推广到求解涡轮平面叶栅跨音带激波流场 ,直接求解全位势方程 ,进一步提高边界区差分格式的精度 ,并注意初场的给法。本法具有节约机时和内存的优点。数值算例表明 ,对于头部不大的涡轮平面叶栅 ,应用本方法…     

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

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

绕10°、20°锥柱体跨音速流动的混合有限差分计算

对于10°、20°锥柱体,用跨肯速小扰动位势流的混合有限差分格式,用跨音速完全位势流的旋转差分格式与简易差分格式,做了数值计算。完全位势方程旋转差分格式与简易差分格式得到的数值结果是一致的,但前者的计算马赫数范围为宽。20°锥柱体完全位势方程的数值结果与实验结果,基本一致,只在顶部与肩部有些差异。小扰动位势方程的数值计算对于10°锥桂体是适用的,对于20°锥柱体足不适用的。     

绕钝前缘翼型定常跨音速小扰动势流的一种计算方法

提高钝前缘翼型的跨音速压力分布计算方法的精度与效率对翼型设计十分重要。国内对跨音速小扰动势流方法进行了研究,但对M_∞>0.8的超临界情况尚未计算。国外文献[3]指出当M_∞>0.9时完全速势方程的一种方法失败了。 本文试图改进经典的小扰动势流方程,探索比较稳定的迭代方法和超松弛方法,以克服超临界流计算常常不易收敛的困难,使小扰动势流计算的应用范围扩大到更高的M     

ISOMS——计算粘性流场的一种隐式差分格式

用稳定化方法求解可压缩完全Navier-Stokes方程(CCNS方程)是计算复杂流场的一种有效的手段。我们在文献[1]、[2]中利用新的开关函数,构造了一种新的混合格式,称为SOMS格式,用于数值求解CCNS方程。本文是将SOMS格式推广到隐式混合格式,称为ISOMS格式(The implicit second order monotone scheme)。它具有如下几个特点:ⅰ)它是隐式格式,稳定性良好;ⅱ)穿过激波是单调的;ⅲ)对于     

带操纵面机翼的非定常跨声速流数值解法

本文介绍一种带操纵面机翼的非定常跨声速流的有限差分计算方法。采用的方程是三元非定常跨声速修正小扰动位势方程,使用时间积分法,求解格式是近似因式分解交替方向隐式(ADI)格式。使用这种方法计算了F-5机翼在来流马赫数为0.9和0.925时的定常气动力和操纵面振荡的非定常气动力。计算结果与国外的NLR试验结果和XTRAN3S方法的计算结果进行了比较,表明计算是成功的。     

跨音速小扰动方程的多重解

本文利用Engquist-Osher格式,离散跨音速小扰动方程(TSD)和边界条件,选择不同初场计算NACA0012翼型二维定常位势流场,得到了与试验值较符合的结果,同时用Murman-Cole非守恒和守恒两格式计算并进行了比较。着重用数值实验进一步探讨多重解问题,发现并总结出了几点很有参考价值的规律和现象。     

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

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

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

超声速主流中逆向喷流流场的数值模拟

本文在差分格式NND-2基础上^[1]，提出了一个修正格式，形式简单，应用方便。从薄层近似的NS方程出发，用该格式计算了超声速主流中存在逆向喷流的钝体绕流的粘性流场，成功地捕捉到流场内的各种波系和涡系结构。计算的结果与文献[2]中的实验结果比较表明二者符合很好。     

旋成体零攻角纵向大扰动势流的AF-2迭代及其并行算法

 <正> 1.引言 为提高跨音速差分计算效率,人们在计算格式的设计、改进方面做了大量的工作,并取得了可喜的成就。并行计算机的出现和发展,使我们有可能进一步提高计算效率。第9期阂赛金等:旋成体零攻角纵向大扰动势流的人F一2迭代及其并行算法A弓n因为并行处理机具有处理数据能力强、计算效益高的特点,但这种机器必须结合具体问题和机器特点加以考虑才能发挥其优势。目前,利用我国设计的“YH一1”(“银河一1”)亿次并行处理机,开展跨音速差分计算的并行算法研究具有理论和运用上的重要意义     

绕钝前缘机翼跨音速势流的松弛计算及其稳定性与收敛性讨论

本文在机翼钝前缘处用精确速势方程和精确的边界条件,其他地方用纵向大扰动而横向小扰动的速势方程和相应的边界条件,联立求解。数值算例1为矩形机翼,展弦比λ=12,翼剖面为NACA0012,自由流的马赫数M_∞=0.63,迎角α=2°,翼根剖面压力分布的计算结果与二元亚音速精确数值解（Sells,1968）接近。算例2为NACA RM A51G31实验的机翼,垂直于1/4弦线的翼剖面为NACA64A010,其后退角χ_1/4=45°,λ=3,根梢比η=2,M_∞=0.4,0.8,0.9,α=2°。计算与实验接近。 本文建立跨音速定常小扰动速势差分方程的线松弛改进迭代在局部线化假设下的稳定性条件和松弛解收敛到原来的微分方程解的条件。这些条件大多数与数值实验相符。     

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

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

跨音速翼型风洞缝壁开度选择初探

对跨音速风洞洞壁干扰问题的研究由来已久。经典的洞壁干扰理论是以考虑压缩性影响的线化亚音速理论为基础的。将此线性理论用于跨音速风洞,存在着根本缺陷。另一方面,由于人们缺乏对透气壁流体流动及粘性损失特性的足够了解,过去的干扰修正理论都只能基于半经验的均匀线化边界条件。随着计算流体力学的迅速发展,出现直接求解跨音速非线性方程确定洞壁干扰的方法,但这类方法,由于使用均匀线化边界条件,其适用性仍十分有限。     

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

可控扩散叶型的扩稳优化

采用正问题方法,将叶型几何参数化、叶型性能分析程序与遗传算法相结合,对某传统可控扩散叶型（CDA）进行叶型损失和攻角范围的综合性能优化设计。结果表明:优化叶型与原叶型具有相近的设计点损失,而其攻角范围由原来的11°增大至17.5°,其中负攻角范围增大了近4.5°。另外,优化叶型的损失随攻角变化也更为平缓,意味着可以在更为宽广的攻角范围内保持稳定的性能。分析表明:吸力面速度峰值位置由原叶型40%弦长处前移至20%弦长处,增加了减速区的长度,使减速更为平缓,是正攻角裕度增大的主要原因。负攻角裕度增加有两方面原因,优化叶型喉道面积增大且喉道位置与设计点吸力峰值位置错开,具有较大的堵塞裕度;压力面前部区域速度较为平缓,甚至略微加速,直至30%弦长后才开始减速扩压,避免了因压力面前缘处的较大速度尖峰以及随后的持续扩压导致附面层的过早分离。      

计算绕薄翼型跨音速非定常流的积分方程法

 本文从小扰动方程出发推导出绕薄翼型的跨音速非定常流动的积分方程,为使此方程适用于具有激波的流场引入人工粘性,并对其数值求解。通过算例讨论了人工粘性和计算域大小对计算结果的影响,与实验及其它数值计算结果比较,表明本方法的准确度令人满意。计算量小,收敛性好是此法的特点。     

翼型近尾迹流动的PIV研究—运动学特性

利用在线式PIV系统(ParticleImageVelocimetry),在低速风洞中对NACA0012翼型在雷诺数2.39×105,0°和4°攻角下的近尾迹流动进行了实验研究。实验结果表明,在较高的雷诺数下翼型近尾迹流动是一种以旋涡的运动学和动力学特性为主导的湍流剪切流。在测量范围内,翼型的尾缘处是近尾迹涡街的形成区;尾缘后0.5倍弦长的区域存在类似于卡门涡街的有序结构,是旋涡发展区域,旋涡具有较好的稳定性;距翼型尾缘0.5倍弦长至1倍弦长的区域,是翼型近尾迹流动由有序走向无序区域,旋涡开始破裂。翼型表面边界层对翼型近尾迹湍流剪切流的演化有重要影响。实验结果还给出了近尾迹流动的平均速度、湍流强度和剪切应变变化率,以及速度脉动量的二阶关联量u'u',u'v'和v'v' 的分布。