多层快速多极子算法的改进措施

 为精确求解散射问题,采用混合场积分方程（CFIE）、多层快速多极子算法（MLFMA）和共轭梯度算法（CG）的收敛技术。基于传统多层快速多级子算法,详细研究了二维拉格朗日插值节点数对计算精度的影响,并改进了插值方法,在不同的层采用不同的插值节点数;提出了在不同的层采用不同的精度控制来计算多级子模式数;分析了稀疏矩阵的对称性对内存使用的影响以及磁场积分方程对迭代初始值的选择。数值计算结果表明以上改进可较大幅度地提高计算精度和计算效率,同时降低内存使用,可满足复杂目标电磁散射计算要求。     

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

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

HB与ISO标准中锥齿轮轮齿弯曲疲劳强度计算标准比较

为正确评估不同标准中的轮齿强度计算结果，分析和比较了航空工业标准（HB）与国际标准化组织（ISO）标准计算锥齿轮轮齿弯曲疲劳承载能力。通过计算方法比较和实例计算对比2种途径，找出了2种标准计算公式的差异、修正系数种类和取值上的差异，及其对计算结果的影响。结果表明：2种标准考虑的影响因素不同，ISO标准相较HB考虑的因素更为全面，特别是在动载系数、应力修正、尺寸系数等参数计算上存在较大差异，使得采用ISO标准计算齿根应力基本值更小，齿根应力与许用齿根应力更大；ISO标准中提供的动载系数根据转速区域的不同采用的计算公式不同，提供的尺寸系数考虑了不同组织和热处理工艺对系数取值的影响，对计算结果更有指导意义。     

湍流模型对叶片表面换热计算的影响

为了比较湍流模型对涡轮叶片外换热计算结果的影响,采用五种湍流模型及两种壁面函数计算了NASA-MarkII导向叶片表面换热并与实验数据进行了对比,结果表明:在相同边界条件下不同湍流模型的计算结果有很明显的差别,即使是同一种湍流模型,如果采用不同的壁面处理函数其计算结果也是有很大差别的;某些湍流模型的计算值只是在某个区域较为理想,还不能找到在整个叶片表面计算结果与实验数据较为接近的湍流模型.在尚无普适性较好的湍流模型的情况下,研究在不同计算域采用不同湍流模型和壁面处理函数的计算技术,是一种较好的可行的方案.      

直升机附件舱温度场的数值模拟

为考查某型直升机附件舱的热控制能力,采用有限元软件分析了其在强制冷却和自然冷却情况下的温度场分布.应用非结构化网格和有限体积法进行计算区域和控制方程的离散,采用标准k-ε模型封闭湍流控制方程.针对3种冷却系统设计方案,分析了不同冷却气流进口尺寸、相对位置以及不同附件散热量对附件舱温度分布的影响,计算了自然冷却情况下附件舱密闭时的温度场.计算和试验结果均表明,所采用的冷却方案可以满足附件舱冷却的需要.      

基于旋转爆震三维流场结构分析的计算模型对比研究

为了研究不同计算模型对三维旋转爆震数值模拟的影响,分别基于Euler方程、N-S方程、RANS方法和IDDES方法并结合滑移、无滑移壁面边界条件,耦合氢气/空气的有限化学反应速率模型（7组分8基元反应）,采用高分辨率的五阶有限差分格式WENO-PPM5离散对流项,对圆环筒型燃烧室内的三维旋转爆震波进行数值模拟,得到了不同模型下旋转爆震波的流场结构、传播特性和推力性能。对比了不同计算模型对流场结构的影响。当采用滑移壁面边界条件时,Euler方程和N-S方程的计算结果较为一致。当使用无滑移壁面边界条件时,边界层的存在会导致可燃混气与燃烧产物之间接触面上的爆燃燃烧区域沿壁面向上游渗透,增大爆燃区域的范围;且爆震波锋面非受限侧变窄、严重变形,不同计算方法计算的变形程度有所不同。     

关于混合气体输运系数计算方案的影响研究

为了研究混合气体分子输运系数的不同近似计算方案对边界层及射流剪切层流场计算的影响,以含双方程k-ωSST湍流模式的质量平均Navier-Stokes方程为控制方程,分别采用两种计算方案:公式计算混合气体导热系数、扩散系数的直接计算方案和由普朗特数和施密特数等关联混合气体导热系数、扩散系数的关联计算方案,对Burrow和Kurkov实验流场进行了数值模拟。结果表明:与采用直接计算方案相比,采用关联计算方案,流场数值模拟在获得相当的计算结果的同时,大大降低了计算量;抛开与控制方程的相容性问题,混合气体扩散系数各直接计算方案对流场参数计算影响不明显;在湍流计算时,边界层/剪切层区域的混合气体湍流输运系数大于分子输运系数,由此弱化了混合气体分子输运系数不同的近似计算方案对流场数值模拟的影响;相对输运系数近似计算方案,湍流可压缩修正模型对组分剖面的影响更为明显。计算表明,在边界层近壁流场、剪切层流场以及层流流场的精细化数值模拟中,混合气体分子输运系数的计算方案有重要的影响。     

翼型和机翼阵风响应的数值模拟

采用有限体积法来数值求解非定常Euler和NS方程,并在非定常计算中通过引入"网格速度"来计入阵风的影响,实现了对刚性翼型和机翼阵风响应的有效数值模拟。首先采用该方法来对翼型攻角突然增大的阵风、陡边型阵风这两种典型阵风情况下的翼型气动力响应进行了计算,结果表明:计算结果与现有的一些理论结果都符合较好,与文献中计算结果也吻合,NS方程的计算结果与Euler方程计算结果在变化趋势上是一致的,但具体的数值有差别。最后针对攻角突然增大的阵风,对不同马赫数下不同展弦比平直机翼的阵风响应进行了数值模拟,得出的计算结果与文献中的计算结果一致。     

TRIP软件的静气动弹性计算模块开发及精度验证

在结构网格流场求解软件TRIP基础上,发展了一套静气动弹性计算模块。首先简要介绍了静气动弹性计算模块总体架构、构成单元和耦合策略,然后详细介绍了构成单元中采用的一些数值计算方法。最后利用大展弦比机翼、DLR-F6翼身组合体模型和HIRENASD机翼模型对静气动弹性计算模块进行了测试和精度验证。大展弦比机翼计算结果表明柔度矩阵、模态叠加和有限元数值求解在模态选取合适的情况下能够得到一致的变形计算结果。DLR-F6和HIRENASD模型不同计算软件的结果一致性表明本文静气动弹性计算模块采用了正确的算法流程,而计算与试验的结果一致性表明静气动弹性计算模块具有很好的预测精度。     

自吸气脉冲爆震发动机超声速进气道气动性能的数值研究

采用有限体积法计算了某轴对称超声速进气道,在不同幅度和频率的出口扰动压强下的进气道内结尾正激波的运动情况,得出了进气道内结尾正激波运动特性和扰动压强的频率和幅度的关系。在计算中,本文采用了多块结构化网格,控制体积的界面无粘通量采用三阶迎风格式插值获得,同时采用了minmod通量限制器,以确保在激波处的解的物理特性;扩散通量采用二阶中心差分格式插值获得。定常计算采用当地时间步法,非定常计算采用双时间步法。离散的代数方程采用交替方向迭代法求解。     

二维激波斜壁问题的数值研究

1.前言 在数值模拟激波构造时,通常将流体看作理想流体,即对Euler方程用有限差分法求解。有限差分法的优点是通用性强,但将其用于计算激波问题时,会在不连续波面附近产生数值衰减、过冲或欠冲,难以进行精确计算。FCT(Flux-Corrected Transport Method)法是Boris等开发的一种新的计算方法,它能充分发挥有限差分法的长处,可有效地计算激渡问题。但是,能对二维复杂形状的绕流进行数值研究还是在Thompson等     

Non-Fourier heat conduction analysis of a 2-D plate with inner cracks at arbitrary direction angles

This paper investigates the problem of a 2-D plate containing internal cracks at arbitrary angles to the heating surface under thermal shock. The finite difference method is applied to solve the temperature distribution and a difference scheme of the oblique crack is developed. For the first time, a non-Fourier heat transfer analysis method of a 2-D plate with inner cracks at arbitrary direction angles is established. Specifically, for a strip with an inner crack paralleling to the surface,numer...     

单边膨胀球面2 元喷管雷达隐身修形研究

针对单边膨胀球面2元喷管,采用斜切出口的方式进行雷达隐身修形设计。利用迭代物理光学(IPO)方法和等效边缘电磁流(EEC)方法进行电磁散射仿真研究,分析喷管修形设计中出口斜切角度对喷管后向散射场的影响,经仿真分析得出,水平极化下斜切角度为25°的喷管具有较低的散射特性,而垂直极化下斜切角度为15°的喷管后向隐身效果较好。     

旋转叶片上的三维附面层

在附壁的非正交曲线坐标系下,见图1。张量形式的流体运动方程见参考文献。 我们假设:(1)流体呈定常运动;(2)忽略彻体力;(3)TSL条件成立。在基本方程中,引入逆变物理分量代替逆变分量,并根据假设进行简化和量级分析,得到一组普遍适用的附面层方程:     

Large-eddy simulation of shock-wave/turbulent boundary layer interaction with and without SparkJet control

The efficiency and mechanism of an active control device ‘‘Spark Jet" and its application in shock-induced separation control are studied using large-eddy simulation in this paper.The base flow is the interaction of an oblique shock-wave generated by 8° wedge and a spatially-developing Ma = 2.3 turbulent boundary layer.The Reynolds number based on the incoming flow property and the boundary layer displacement thickness at the impinging point without shock-wave is20000.The detailed numerical approaches were presented.The inflow turbulence was generated using the digital filter method to avoid artificial temporal or streamwise periodicity.The numerical results including velocity profile,Reynolds stress profile,skin friction,and wall pressure were systematically validated against the available wind tunnel particle image velocimetry(PIV) measurements of the same flow condition.Further study on the control of flow separation due to the strong shock-viscous interaction using an active control actuator ‘‘Spark Jet" was conducted.The single-pulsed characteristic of the device was obtained and compared with the experiment.Both instantaneous and time-averaged flow fields have shown that the jet flow issuing from the actuator cavity enhances the flow mixing inside the boundary layer,making the boundary layer more resistant to flow separation.Skin friction coefficient distribution shows that the separation bubble length is reduced by about 35% with control exerted.     

轴对称跨音速进气道外罩绕流的高效差分算法

 本文考虑轴对称跨音速进气道外流场的有限差分计算方法。采用守恒型位流方程、贴体坐标网格和精确边界条件,根据最佳收敛准则对轴对称情形设计出新的近似因式分解迭代格式,并将所设计的格式应用到计算皮托式进气口的轴对称跨音速流场。对几种典型进气口的计算表明,本文格式收敛快,计算结果与实验符合很好。     

高阶精度非线性格式WCNS-E-5在二维流动中的应用研究

本文采用具有5阶精度的加权紧致非线性显式格式(WCNS-E-5)对定常与非定常二维流动进行数值模拟,研究表明该格式对各类间断有很好的分辨捕捉能力,而且对强间断如激波的计算,即使在高马赫数与高雷诺数条件下它仍具有很好的收敛性与可靠的计算结果.此外,WCNS-E-5在粗网格条件下也体现出优越性.类如WCNS-E-5的高精度激波捕捉方法将为以后开展湍流数值模拟工作提供坚实的技术保证.     

粉尘气体中平面激波在固壁反射的高分辨率数值模拟

本文在稀颗粒群双流体简化模型下，采用算子分裂技术，综合应用高分辨率ＧＲＰ格式，全二阶精度的ＦＣＴ格式及预估校正格式，取三种不同粒相固壁边界条件，数值模拟了粉尘气体激波在刚性固壁上的正反射，确证了由理想平衡混合气体分析所预示的三类激波反射结构，给出了非平衡反射时波后流场的松弛特性，并评价了不同的粒相固壁边界条件假定的适用性。     

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

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

非定常环境中叶栅总压恢复系数的数值研究

利用以Ni格式和隐式双时间步法求解二维非定常N-S方程为基础的数值模拟方法研究了流动的非定常性对涡轮叶栅总压恢复系数的影响,并与文中提出的定常模型和准定常模型进行了对比。结果发现,非定常流动能够提高叶栅的总压恢复系数,但这是建立在与准定常流相比较的层面上,而和定常流比较时,非定常流只会使总压恢复系数下降。