跨音压气机叶栅流场的附面层修正

本文发展了二维压气机叶栅亚音、跨音流中粘性——无粘相互作用程序。无粘计算在弦向大扰动、弦法向小扰动的假设条件下用跨音松弛法求解;粘性附面层用Nash方法计算。考虑了叶栅尾迹的影响。最后,连同通道区和尾迹区一起,求得了包括冲波影响的叶栅的总损失系数。计算结果表明,在跨音情况下,粘性对冲波强度、位置以及全流场的速度分布都有不可忽视的影响。     

跨音转子旋转失速现象研究

本文研究了跨音转子的旋转失速现象。用旋涡理论建立了旋转失速工况下的涡系流动模型及流场算方法,给出了低速转子和跨音砖子的流场计算实例。用热线风速仪和脉动压力测试仪相结合测定了跨音转子的速度场及特征参数。分析和研究了跨音转子中旋转失速现象发生、发展和消失过程及其流动特点,其中发现了一些新的,有研究价值的流动现象。     

一种跨声速定常流场求解加速方法

跨声速定常流场的隐式求解相当于使用牛顿迭代法求解一个非线性方程组。为满足牛顿迭代收敛性的要求,通常需要对所求解问题进行全局化处理。在同伦延拓的框架内,提出了一种基于拉普拉斯算子的方程延拓方法,提高了定常流场隐式求解收敛速度。针对定常流场通常初始化为均匀来流的特点,一方面利用拉普拉斯算子的椭圆性加快边界条件信息向流场内部的传播,另一方面利用拉普拉斯算子的线性和正定性改善延拓问题的正则性,综合两者增加拟牛顿算法的稳定性,提高可用CFL数,最终达到提高流场求解效率的目的。由于流场问题的复杂性和非线性,难以通过理论分析得出先验的最优非线性求解策略。因此,通过无黏NACA0012翼型、湍流RAE2822翼型和三维ONERA M6机翼等算例的数值实验,研究了拉普拉斯项参数对收敛效率的影响,给出了效率较优的参数组合,验证了本文方法在跨声速情况下相对于经典伪时间推进法可以节约20%以上的CPU计算时间。      

对旋轴流风机的稳态流场数值模拟

在任意曲线转动坐标系及重整化群(RNG)方法改进的k-ε两方程湍流模型基础上,提出了"通用进出口边界条件"的给法,并利用改进的强隐式(MSIP)方法进行计算迭代,提高了对多排叶轮机数值模拟的稳定性和准确性.通过对某对旋轴流风机定常全三维流场的计算,并与设计参数和实验性能参数进行对比分析后,验证了在风扇/压气机的设计中,不宜在叶尖区分配最大负荷,这种设计对对旋轴流压缩系统尤其不利.     

多块网格计算跨音压气机带叶尖间隙流场

开发了三维数值模拟程序研究轴流跨音转子叶尖间隙流动,应用高雷诺数k-ε湍流模型加壁面函数的方法,计算了轴流跨音转子NASA Rotor37在设计转速下的流场.叶尖间隙采用分区的H型网格和主流区连续对接耦合计算,没有用间隙模型,也没有考虑Vena收缩效应而减小间隙量.在用有限体积法对Navier-Stokes方程和湍流方程进行空间离散的过程中采用了交错网格的方法将N-S方程与湍流方程紧密地耦合在一起,从而提高了计算精度.计算结果和实验数据进行了详细的比较和分析.结果表明,中部叶展具有与实验结果非常一致的流场特征,根、尖区流场则因涡粘假设和激波问题的存在而使流动细节与实验结果略有偏差.     

跨音叶栅流场的相仿律

跨音叶栅绕流参数繁多,为了便于系统研究参数影响,也为了扩展某个气动性能良好的叶栅的使用范围,研究不同流场之间的相仿可能是有益的。多年前VonKarman曾给出了二维孤立翼型的跨音相仿律。然而,由于远前后方边界条件的不对称性,以及叶轮机械流场的周期性 所以叶栅的相仿律和孤立翼型的尚应有所不同。本文指出,为了维持两个叶栅流场的相仿,五个对应跨音相仿参数应分别相等。文中举例表明,相仿性有助于由一个巳知给定叶栅流场迅速得到很多在不同条件下运转的相仿跨音叶栅流场。     

三维抛物化N-S方程全隐格式矢通量分裂求解

采用全隐格式的矢通量分裂算法,在通用曲线坐标系内求解三维无因次抛物化强守恒型N-S(Navier-Stokes)方程.所建立的计算程序,可计算较为复杂边界的三维流场,如求解亚、跨、超音速共存的流场,分区模拟内外流共存流动及其混合剪切尾流.给出了一轴对称亚跨超音速喷管内流场的计算结果,包括计算网格,流矢分布,压力、温度和马赫数等值线分布.也给出了一喷管内外流及其混合尾流的计算结果,包括计算网格,流矢分布,无因次压力、密度和温度等值线分布.     

气流掺混的数值模拟

采用显式的TVD MacCormack格式的时间分裂形式求解了平板上横流流场的二维N-S方程和组分输运方程,得到了流场中气体超音速流动和混合的结构图谱.为了获得对近壁附面层分离较好的模拟效果,采用了Lam和Bremhorst的低雷诺数ｋ-ε模型来模拟流动中的湍流粘性.计算结果表明,该方法可以用来求解横流流场中的气流掺混问题.     

基于湍流尺度的混合RANS/LES模型

比较BL(Baldwin-Lomax)湍流模型的内层模型与Smagorinsky亚格子模型,两者的主要差别在于湍流尺度的计算方法不同,而RANS(Reynolds Average Numerical Simulation)方程组与LES(Large Eddy Simulation)方程组在形式上又是相同的.基于这两点,发展了RANS/LES混合模型.该模型是在网格滤波尺度与到壁面的距离相等处将流场分区,在近壁区直接求解RANS方程组,在远离壁面区域,求解LES方程组.为了保证湍流粘性系数在内外层分界处的光滑过渡,内外层湍流粘性系数通过混合函数连接.混合方法减轻了对壁面附近网格分辨率的要求和时间步长的限制,因此比LES所需的计算成本少.用该混合模型模拟了绕双椭球的高超音复杂流场,结果表明采用混合模型可以准确地模拟可压缩湍流流场.      

一种战斗机外形涡流场Euler方程数值模拟

 用三步显式格式时间推进求解物理空间曲线网格上有限体积离散的Euler方程,数值模拟一种战斗机外形的涡流场.计算的机翼表面压力分布与实验符合较好,用Euler方程捕获到机翼前缘分离涡(主涡)及翼面上的二次分离涡.Euler方程解中出现的二次涡可能是逆压梯度、尖锐边缘和人工粘性共同作用的结果,它的流谱与实验定性符合.     

二维微槽道电渗流的增强混合

微尺度流动中的混合一直是微流系统研制中的一个大问题,对壁面zeta势非均匀分布情况下的二维微槽道电渗流进行了数值分析,重点分析了这种情况下的微流动混合增强.基于电渗流的Helmholtz-Smoluchowski滑移速度模型,应用SIMPLE方法计算了4种滑移速度在整个壁面上的分布不均匀一致时的微槽道电渗流动的流场及被动标量在其中的混合.计算结果表明,滑移速度的非均匀一致导致微槽道电渗流中的流线发生剧烈的转折,从而使得微槽道流动的混合效率有很大的提高.流场中流线的转折越剧烈、转折的次数越多,微槽道流动混合效率的提高就越大.     

谐波平衡法在低速非定常流模拟中的应用

谐波平衡法是一种有效的周期性非定常流的计算方法.采用基于可压缩流的谐波平衡方程在计算低速不可压流动时,会由于对流通量计算格式中的数值粘性污染,降低解的精度和收敛性.采用预处理技术,使得基于可压缩流的谐波平衡方程可以直接用于低速周期性非定常流的计算中.选取典型的不可压方腔驱动流和低雷诺数圆柱绕流为例进行了时间推进法和谐波平衡法的计算对比.计算结果表明预处理后的谐波平衡方程适合于低速流的计算,在谐波平衡法中采用较少阶数的谐波计算就可以还原出几乎准确的非定常流场.     

低速流预处理Roe格式中的数值粘性

使用Roe格式计算低马赫数流动,需采用预处理技术.在预处理技术中,需要选择适合的预处理参数,但有时很难找到,因而得不到合理的计算结果.通过数值实验,发现使用预处理方法时,格式的数值粘性对低速流动的模拟有显著的影响.提出了在预处理Roe格式中加入可调参数,用以控制数值粘性对计算结果的影响.带有可调节数值粘性参数的预处理Roe格式的数值粘性更接近低速流动的真实耗散,反映了低速流动的真实物理现象.对方腔流动和圆柱绕流的低雷诺数定常流动所做的数值实验,得到了比较满意的结果,表明这一改进是有效的.对同一低雷诺数不同低马赫数的定常圆柱绕流所做的数值实验,均得到了满意的结果,表明了方法的适用性.     

高超音速热流数值计算中的误差匹配原则

通过对粘性边界层的数值模拟和误差分析,提出了差分格式的误差匹配原则.在边界层内的计算网格一般都具有很大的长宽比,而计算误差与网格长宽比密切相关.为了使不同坐标方向上的误差保持在同一量级上,需根据误差匹配原则确定计算网格的尺寸.对绕钝锥的高超音速流动进行了数值模拟,进而计算了钝锥表面的热流分布.数值计算表明,参照误差匹配原则进行网格划分是有益的,可在保持流场数值计算精度和提高热流计算准确性的同时,使网格的数量和分布更为合理.     

数值模拟高超音速化学非平衡流的隐式格式

针对高超音速化学非平衡流提出了两种隐式格式.一种是近似因子分解算法对角化形式在隐式TVD(全变差递减)格式上的推广,另一种可视为隐式TVD格式的LU分解形式.数值试验表明,本文给出的算法保持了TVD格式的高分辨率特性,同时计算量大大减少,是高超音速化学非平衡流数值模拟的有效方法.     

Orbital error propagation considering atmospheric density uncertainty

Due to the influence of various errors, the orbital uncertainty propagation of artificial celestial objects while orbit prediction is required, especially in some applications such as conjunction analysis. In the orbital error propagation of artificial celestial objects in low Earth orbits (LEOs), atmospheric density uncertainty is one of the important factors that require special attention. In this paper, on the basis of considering the uncertainties of position and velocity, the atmospheric density uncertainty is also taken into account to further investigate the orbital error propagation of artificial celestial objects in LEOs. Artificial intelligence algorithms are introduced, the MC Dropout neural network and the heteroscedastic loss function are used to realize the correction of the empirical atmospheric density model, as well as to provide the quantification of model uncertainty and input uncertainty for the corrected atmospheric densities. It is shown that the neural network we built achieves good results in atmospheric density correction, and the uncertainty quantization obtained from the neural network is also reasonable. Moreover, using the Gaussian mixture model - unscented transform (GMM-UT) method, the atmospheric density uncertainty is taken into account in the orbital uncertainty propagation, by adding a sampled random term to the corrected atmospheric density when calculating atmospheric density. The feasibility of the GMM-UT method considering atmospheric density uncertainty is proved by the further comparison of abundant sampling points and GMM-UT results (with and without considering atmospheric density uncertainty).     

导弹“×”翼状态翼面流动分析

本文从流态结果的研究出发提出对导弹“×”翼状态翼面流动的分析方法,指出有迎角和侧滑角情况下四片翼面各不相同的有效迎角和前缘有效后掠角对于翼面流动以及气动力的影响。本文的分析方法对于分析“×”翼在有迎角和侧滑角情况下的操纵面效能及其导弹的制导有一定的帮助,也可以为有关计算模型的建立提供一些物理依据。     

一种利用正交流线坐标系求解跨音速管流的新方法

本文以正交流线坐标系为工具,采用位流模型,导出了正交流线坐标系下平面或轴对称流动的控制方程组,利用流函数与位函数的存在性,进而导出了两个首次积分。设计了一套求解匀熵管流的分凝迭代方案,模拟了管内无激波的平面或轴对称喷管的跨音速流动。数值实验表明该方法能够处理小喉道曲率半径平面或轴对称喷管跨音速流动,与已有的数值计算结果及实验结果符合较好。