三维可压缩非定常流的壁面分离判据及其分离线附近的流动形态

通过引入依赖于密度的物面法向速度变换wr=-1/ρh1h2∫ozh1h2(e)ρ/(e)tdz,描述物理速度空间(ua,va,wa=u,v,w+wr)具有无滑移壁面条件的三维可压缩非定常连续方程可转换成变换速度空间(u,v,w)内具有无滑移条件定常连续方程.因此,采用定常壁面分离的分析方法和结论,再通过变换和研究wr的贡献,给出了三维可压缩非定常壁面分离的判则以及分离线附近的流动形态.研究指出,二维和三维情况下,都出现伴有壁外附着的壁面分离情况.数值模拟证实了理论和结论.     

适用于参数可调结构的非定常气动力降阶建模方法

基于计算流体力学(CFD)的非定常气动力降阶模型(ROM)可以极大提高气动弹性分析效率,然而现有的ROM只能针对固定参数结构,即只适合于固定模态振型,这使得现有ROM在气动弹性优化设计和不确定性分析等结构变参问题中应用受限。针对该问题,在文献[20]基础上提出了一种新的适用于任意模态振型的非定常气动力建模方法。首先将待设计/分析的结构进行参数化抽样和模态分析,之后通过主成分分析(PCA)得到若干基振型,再将这些基振型线性叠加即可拟合抽样空间内任何参数下结构的前若干阶振型。当结构参数改动时,仅仅是叠加系数发生变化。算例表明,仅用很少的基振型就能达到理想的拟合精度。经典的气动力降阶方法可用于基振型坐标下的气动力降阶,进一步变换可得到适用于不同结构的ROM,这意味着,结构参数可以在抽样空间内任意调节改动,而ROM却是通用的。该方法能广泛用于气动弹性优化设计和不确定性分析工作,可提高颤振分析精度和效率。     

隐式Runge-Kutta方法在非定常化学非平衡流动中的应用

推导了求解NS方程组的隐式龙格-库塔方法,并从稳定性、精度和效率三方面展开讨论。通过定常和非定常算例对该方法进行验算,分析表明该方法计算结果正确、稳定性好、精度高、效率高,适合求解大型复杂非定常化学非平衡流动。     

基于线性弹性体的网格变形方法研究

网格变形方法在计算流体力学的各个方面都发挥着极其重要的作用,是各种应用问题解决方案的基本组成部分.研究和实现了基于线性弹性体的网格变形方法,对基于弹性体法的网格变形方法进行了综述,从弹性方程理论、高效求解方法等方面进行了介绍,给出了方法的最优参数选择,并从基本变形能力、多段翼型设计、非定常动态特性模拟等方面,展示基于线性弹性体网格变形方法的研究和应用进展.结果表明:实现的线弹性体网格变形方法特别善于处理几何大变形,尤其是网格扭转比较大的应用,并且变形后网格品质高、鲁棒性强.     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

Contribution of geometry of interaction between interplanetary and terrestrial magnetic fields into global magnetospheric state and geomagnetic activity

The paper presents results of our study of dependence of geomagnetic activity from geoeffective parameters taking into account mutual orientation of the interplanetary magnetic field, electric field of the solar wind and geomagnetic moment. We attract a reconnection model elaborated by us made allowance for changes of geometry of the solar wind–magnetosphere interaction during annual and diurnal motions of the Earth. We take as our data base the interplanetary magnetic field and solar wind velocity measured at 1 a.u. at ecliptic plane for the period of 1963–2005 and K_p, D_st, a_m indices. Taken as a whole a geoeffective parameter suggested by us explains 95% of observed variations of the indices. Changes of the geometric factor determined by mutual orientation of the solar wind electric field and geomagnetic moment explain larger than 75% of observed statistical variations of D_st and a_m indices. Based on our results we suggest a new explanation of semi-annual and UT variation of geomagnetic activity.     

对旋式通风机压力脉动特性的数值研究

对矿用对旋式通风机内部的三维非定常湍流流动进行了数值模拟,采用FFT方法得到了包括通风机进口、集流器、前级电机、两级叶轮、后级电机、扩散器和通风机出口等在内的共10个监测面上监测点处压力脉动的时域和频域分布。分析结果表明:对旋式通风机两级叶轮及其附近区域是流场中压力脉动较强的区域;前后两级叶轮之间的动-动干涉、两级叶轮与其相邻部件之间的动-静干涉以及流动的周期非定常特性是通风机内部压力脉动产生的主要根源;设计良好的集流器和扩散器可以显著降低通风机进出口处的压力脉动。研究结果可为通风机的噪声预估和控制提供参考。     

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。     

仿鸟柔性扑翼气动特性与能耗的数值研究

建立了适当的三维仿鸟柔性扑翼模型,并以配平重力和平衡阻力为条件,数值计算了它的低雷诺数非定常流场.研究揭示了翼面初始扭转角度、动态俯仰幅度等重要设计参数与飞行性能的关系,表明扑翼平面的初始扭转程度、扑翼柔性材料的选择以及两者之间的合理搭配对扑翼机的成功飞行至关重要.研究分析了仿鸟扑翼的流场涡结构、升力推力产生原理,下扑过程附着上翼面的前缘涡是升力产生的重要机制.对扑翼气动功率的比较分析也发现,人造扑翼机需要的气动功率明显高出同等大小的鸟类,在效率方面尚不及扑翼飞行生物.