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

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

纵向大扰动方程的变参数型变分裂法

 本文用变参数型变分裂法,解定常跨音速纵向大扰动方程。计算翼型气动力。对变参数型变分裂法以及它的参数选取方法作了必要的介绍。对NACA0012翼型,M_∞=0.8,α_A=0°;M_∞=0.95,α_A=0°;M_∞=0.75,α_A=1°,进行了计算。所得结果与其它结果进行了比较,表明本方法是快速收敛的、可靠的。     

高超声速无粘流场的推进计算

本文将文[1]中用于时间相关法计算的NND格式推广到定常超声速流动的空间推进计算,采用二步的预测、校正方法保证了推进方向的二阶精度,可以证明,这种二阶精度的NND格式具有TVD性质,是MacCormack二步显式格式的推广。本文首先将格式应用于二维平板上斜激波反射流场的推进计算,以检验格式捕捉激波的能力,同时研究了不同的通量分裂方法对格式捕捉激波能力的影响,得到了相当满意的结果。在此基础上,计算了航天飞机简化外形的身部超声速流场,给出了M_∞=10,α=0°,和M_∞=5,α=5°两种状态的部分结果,计算结果清楚地描绘了由于气流在机翼附近受到强烈压缩而产生的内嵌激波与外激波相交的复杂流场结构,与文[7]相比,流场结构更为清晰。     

新型高超音速静压探针的设计和实验

参考了近40年来国外有关超音速和高超音速静压探针的研究成果,设计了多种并发展了二种高超音速静压探针。它们可精确地用于持续工作流场仅为30ms的流场静压测量。实测静压与自由流静压比p_(s,m)/p_(s,∞)=1.009;马赫数比M_(∞m)/M_(∞,)=0.993。实验在加拿大多伦多大学宇航研究院的高超音速炮风洞里完成。炮风洞喷管出口气流特性是:(空气)自由流总压比_(t,∞)=25.5MPa,总温T_(t,∞)=1000K;马赫数M_∞=8.30;自由流静压p_(s,∞)=2.13kPa;雷诺数R_e=3.2×10~7。以探针直径D为基准的R_(e,D)=34,000～52,800。这种探针己被用于“管内高超音速流场发展实验研究”中了。     

小钝头锥无粘平衡、非平衡绕流流场的数值分析

本文给出钝头细长锥无粘、平衡和非平衡高超声速绕流流场的数值分析。对存在化学反应的空气,考虑了13种气体组分和22个化学反应,包括N_2、O_2和NO的振动松弛与离介过程之间的耦合反应。头部亚-跨声速区和后身高超声速区分别用直线法和有限差分法进行计算。在来流条件M_∞=20～26,Re_∞=(p_∞u_∞α)/μ_∞=1.5×10~3～5.5×10~6(α为头部半径)的参数范围内,给出了绕流流场和等离子体鞘层诸特性,如激波形状、电子密度和电子碰撞频率等。     

高超音速后掠激波与边界层干扰流场特性

应用表面油流和液晶温度显示、表面热流率和压力测量四种测试技术,在Ma_∞=7.8、单位长度雷诺数尺Re_∞=3.5×10~7/m条件下,研究了30°迎角尖前缘翼面诱导激波和湍流边界层干扰流场特性。结果表明:表面流动是准锥型,在强相互作用下,流动发生了二次分离和再附,导致表面油流线的流向角、温度、热流率和压力分布均出现明显的凹坑。     

低温烧蚀材料的烧蚀实验研究

本文给出了低温烧蚀材料的高超声速风洞烧蚀实验研究的结果。模型外形是半锥角为9°的球锥体,头部半径为15毫米和25毫米,材料为蜂蜡和樟脑,实验条件是M_∞=6,α=0°、2°、6°和8°,Re_∞为(1.90～4.42)×10~7米~(-1)。实验结果表明:随着Re_∞的增加,可以获得层流、转捩和湍流烧蚀表面形态和端头烧蚀外形;证实了在实验过程的后期出现平衡外形,在平衡外形中存在角点;融熔型材料和升华型材料的烧蚀表面形态差别较大;有攻角和零攻角的烧蚀外形有较大差异。     

高超声速稀薄气流中球阻力系数的一种简单公式

本文提出了以粘性参数K=M_∞/√Re_D为变量的球阻系数的简单关系式。它适用于高超声速稀薄气流的情况,在Re_D=5×10~2～10~5;M_∞=8～25的范围中与实验数据吻合。此外,按此式计算的球阻系数而得到的标定球的动压值与由皮托压力测定所得结果一致。     

二维凹槽超声速湍流流动数值模拟

本文采用时间、空间均为二阶精度的NND差分格式,以及Menterk ωSST二方程湍流模型,数值模拟了二维开式凹槽超声速粘性流动。用组分扩散方程计算了凹槽内粒子驻留时间。本文给出了M∞=3,长深比L/D=3,后壁倾斜角为90°和30°两种几何形状凹槽的计算结果,与相应的实验符合甚好。计算亦表明,该模型较高地估计了涡粘性,从而使流动更稳定。     

超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

驱动方腔内涡流的数值模拟

本文介绍一种求解涡度流函数形式的定常不可压缩Navier-Stokes方程的半隐式指数型差分格式。涡度方程用上述格式求解,而流函数方程用多层网格法求解。这种组合求解方式具有很好的稳定性及较快的收敛速度。本文对Re=100,400,1000,3000,5000的驱动方腔及Re=400,1000的驱动长方腔进行了数值模拟,结果与现有的计算相吻合。     

钝体化学非平衡粘性流数值模拟

用有限差分法求解Ｎ Ｓ方程,数值模拟了二维钝楔头部和近尾迹非平衡化学反应流场。模拟条件为Ｍａ∞＝２５,飞行高度ｈ＝７０ｋｍ,半楔角θ＝２０°,钝度比为０６０;气体组分为７个组元,共１１个化学反应式。头部流场结果与已有计算结果符合较好。在Ｍａ∞＝２５条件下,底部近尾迹较大区域内具有较高的温度和较大的电子密度,且下降趋势缓慢。     

DCS5在计算气动声学中的应用研究

对流项采用五阶线性耗散紧致格式(DCS5)离散,粘性项采用6阶中心紧致格式离散,时间项采用5级5阶龙格-库塔时间推进格式求解流体动力学方程组,并由此建立近声场的直接模拟方法。为了应用结构网格在复杂几何形状条件下对流场进行计算,发展了高精度的多块特征对接边界技术。完成了低速无粘振荡圆柱和粘性圆柱绕流等标准计算声学算例,获得了令人满意的结果。     

绕三维物体亚音速振荡位流及其在飞行器动导数计算上的应用

 本文介绍处理不同外形三维物体亚音速振荡绕流一种统一的方法。本方法的主要特点是采用有限元素法直接解由格林定理导出的物面速势积分-微分方程以求得物面的速势分布,然后再用有限差分法对速势进行微分求物面的压力分布。 由于本方法理论上比较严格,适用于复杂外形物体绕流的计算,所得结果又比较准确,因此近几年来在国外得到越来越广泛的应用。在本文中,采用了与有关文献相同的基本方程,但在气动影响系数的计算上略有不同,本文并将这一方法应用于飞行器动导数的计算,所得的结果与实验结果符合。     

三维高超声速喷流干扰流场的数值模拟

采用ＬＵ隐式方法求解了有限体积法离散的完全Ｎ－Ｓ方程,数值模拟了三维高超声速喷流干扰流场。计算中采用了Ｊａｍｅｓｏｎ的当地极值递减（ＬＥＤ）格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,与实验结果相比较,数值方法对附体流动计算较准确,对分离区附近流场计算有一定误差。     

紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。     

有侧壁干扰的机翼半模型风洞实验的N-S方程数值模拟

应用三维可压、雷诺平均Navier-Stokes方程数值模拟了机翼半模实验风洞侧壁干扰和三维机翼半模与安装侧壁结合部流场。计算采用中心有限体积多步Runge-Kutta时间步长格式,将Baldwin-Lomax两层代数紊流模型扩展应用于三维拐角区流动。用本文方法计算了3个算例 ,并与国外的有关实验数据进行比较,计算和实验吻合良好。     

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

本文在机翼钝前缘处用精确速势方程和精确的边界条件,其他地方用纵向大扰动而横向小扰动的速势方程和相应的边界条件,联立求解。数值算例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°。计算与实验接近。 本文建立跨音速定常小扰动速势差分方程的线松弛改进迭代在局部线化假设下的稳定性条件和松弛解收敛到原来的微分方程解的条件。这些条件大多数与数值实验相符。     

TVD格式数值模拟斜腔内旋涡运动

 采用ＴＶＤ有限体积格式,通过求解非定常可压缩Ｎ－Ｓ方程,对恒定出入口压力比驱动下的斜腔内部流动进行了数值模拟。计算后得到清晰的流场涡系结构及其发展过程；显示了空腔流动中的自维持振荡现象；并分析了剪切层不稳定对周期性振荡的作用机制。