昆虫悬停飞行的动稳定性:理论分析与数值模拟

首先,从昆虫(身体和拍动翅)的动力学方程和N-S方程出发,在一定假设下,将运动方程简化为6自由度刚体的方程,并用线化理论给出了反映昆虫悬停飞行纵向动稳定性的理论解.然后,用完全的动力学方程和N-S方程偶合的数值解对上述理论的简化假设进行了检验.     

三维控制翼周围粘性分离流动的数值模拟

本文首先分析了中心型隐式格式和迎风型隐式格式的优缺点,在此基础上发展了一种自由参数可确定的反扩散隐式格式。用此格式计算了绕三维控制翼的层流和湍流分离流动。出发方程是三维简化N-S方程,用Baldwin-Lomax湍流模型。结果表明,这种隐式格式效率高,能满意地刻划流场。     

ALE方法求解圆柱的涡致振动

数值求解基于拉格朗日-欧拉(ALE)描述的不可压缩流体的N-S方程,计算了较小雷诺数下圆柱的涡致振动现象.N-S方程的对流项和扩散项分别采用三阶迎风紧致格式和四阶中心紧致格式离散,计算网格采用H-O非交错网格系统,并结合分块耦合方法.柱体运动简化为弹簧-阻尼-质量系统,柱体运动方程采用经典龙格-库塔方法求解.通过模拟柱体和流体之间的非线性耦合作用,成功地捕捉到了"锁定"、"拍"和"相位开关"等现象,并与试验数据相吻合.另外,本文详细分析了Re＝200时圆柱尾涡形态、升力、阻力以及圆柱位移等随圆柱自振频率变化的过程,捕捉到尾涡结构变化的频率转折点.     

三维层流分离流动Navier-Stokes方程的数值解法

本文研究了配平翼和直角板的三维层流分离流动。利用反扩散的显隐式差分格式和显式混合差分格式对N-S方程和简化N-S方程作了计算,求得了流场的物理量分布。与已有数值结果比较,本文的结果是满意的。     

超声速二维湍流分离流动的数值解法

本文采用层流中引导简化N-S方程的方法,建立了简化的湍流运动方程组,分别用零方程、一方程模型计算了绕二维压缩拐角的湍流分离流动。数值求解用显隐式的反扩散格式。     

单转子风扇的三维反问题气动设计

采用准三维流函数反问题设计方法与三维 N-S求解方法的相互迭代对单转子风扇进行气动设计。在叶片初步计算中得到叶栅进出口气流角沿径向分布 ,并将它作为本文的目标函数。采用准三维反问题求解方法 ,依次构造出各个 S1流面上的叶片几何形状和气流角分布。然后再采用 N-S方程的求解方法 ,对叶片进行全三维流场的数值计算。通过 N-S方程计算结果与目标函数的对比 ,重新修正叶片出口气流角分布 ,并作为下一次反问题设计的目标函数。经过反问题与 N-S方程求解的反复迭代 ,最终得到满足设计要求的叶型。     

二维舵面绕流的N-S方程数值计算研究

本文基于结构化网格,应用Baldwin-Lomax湍流模型进行绕舵面流动的二维N-S方程数值模拟,将部分计算结果与风洞试验结果进行了对比分析研究,分析了计算方法的可靠性和计算结果的合理性.为开展三维舵面绕流的N-S方程数值模拟和铰链力矩数值计算研究奠定了基础.     

二维可压粘流N—S方程的边界元解法

本文给出了一种求解二维可压Navier-Stokes(简称N-S)方程的边界元法。使用线化技术把控制方程变成变系数的线性偏微分方程,其基本解由组合代数的方法构造,从而给出了N-S方程的积分表达式。本文计算了几例低R_e数的跨声速翼型绕流,结果表明边界元法可以应用到可压粘流中。     

超声速飞行器侧向喷流干扰流场传统数值模拟方法的误差分析

使用CFD方法,分别就真实喷管边界和简化喷口边界,计算超声速飞行器侧向喷流干扰流场,研究边界条件对干扰流场及气动力的影响.使用k-ε湍流模型封闭雷诺平均N-S方程,利用非结构网格对流场进行空间离散.通过对比,计算结果与实验值吻合良好,证明该方法具有一定可靠性.进一步研究表明喷流边界条件对喷流干扰流场具有一定影响:相对于简化喷口边界,真实喷管边界喷流出口的非均匀性导致喷口上游分离涡和激波位置较为靠前,从而引起附加气动力和力矩的变化;由于摩擦阻力的作用,真实喷管静推力存在损失;喷流压比为500时,总法向力和总俯仰力矩在两种边界条件之间的误差分别为8.21%和22.4%,误差较大.在进行侧向喷流干扰流场的精确计算时,需要考虑边界条件的影响.      

冲压发动机进气道压力振荡过程的数值研究

针对来流马赫数3.5,0°迎角飞行条件下冲压发动机的简化模型,数值模拟了进气道整流罩开启过程的非定常流动。计算采用AUSMPW有限体积格式,基于结构与非结构的混合网格技术,求解了非定常N-S方程。计算结果表明,在进气道整流罩打开以后,超声速气流会在燃烧室通道内形成振幅衰减较慢的低频剧烈振荡。在同样计算条件下采用Euler方程进行对比性研究,数值实验发现对这类大雷诺数、强非定常性的问题,Euler方程和N-S方程得到的结果相差不大。     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

非惯性非正交曲线坐标系中的边界层方程

应用张量分析方法将 N-S方程写于作非惯性运动的非正交曲线坐标系中 ,然后用匹配渐近展开法导出此坐标系中的边界层方程。消去压强项后 ,与惯性坐标系中的边界层方程相比 ,非惯性系中的边界层方程在 x1和 x2 方向的动量方程中分别只多了一哥氏力项 ,它们是由物体转动角速度的物面法向分量产生的 ,其它惯性力项均不进入边界层方程。方程中所含曲线坐标系的拉梅系数及其导数只需在物面上取值 ,物面曲率影响不进入边界层方程     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

跨声速非等熵势函数方程的数值求解

基于动量方程的特点,本文将只与温度有关的“等熵密度”与熵的变化分离,而由散度型动量方程直接计算出气体经过激波的熵增。与经典势函数方程迭代,可以方便地求解非等熵跨声速流动。从对任意回转面叶栅的跨声速流动计算来看,这种非等熵势函数得到的激波位置可比等熵势函数前移1～2个网格,强度有所减弱,和试验结果较接近。由于各流线的熵不相等,对Kutta条件和出口边界条件作了修正。本文建议的方法也可用于计算流函数的熵增。     

球柱三维分离及涡流的数值模拟

对三维可压非定常Ｎ－Ｓ方程进行有限体积中心离散,并用Ｒｕｎｇｅ－Ｋｕｔａ显式格式进行数值求解。数值模拟了球柱分离及涡流的形态结构,给出了具有激波、边界层、大分离区、涡流及其间相互干扰的复杂流型     

一种求解核函数方程的新方法—辅助核法

本文提出了求解核函数方程的辅助核法。此法克服了偶极子栅格法对网格分布的限制,具有更大的适用范围。数值计算表明,此法对求解绕翼型的跨声速激波流动是很有利的。     

发动机经验故障方程的建立与应用

本文提出了根据发动机故障样本建立故障方程的方法。发动机经验故障方程是不同于发动机小偏差方程的另一类故障方程,它为发动机故障诊断提供了一条简便易行的途径,并且使发动机故障诊断范围扩展到气路分析方法难以适用的场合。文中讨论了经验故障方程的建立方法及其在发动机故障诊断中的应用,给出了利用经验故障方程进行发动机故障诊断的实例。      

NUMERICAL STUDY OF TOPOLOGICAL STRUCTURE OF 3-D TRANSONIC VISCOUS FLOW FIELD (TVFF) INSIDE TURBINE CASCADE

ＮＵＭＥＲＩＣＡＬＳＴＵＤＹＯＦＴＯＰＯＬＯＧＩＣＡＬＳＴＲＵＣＴＵＲＥＯＦ３ＤＴＲＡＮＳＯＮＩＣＶＩＳＣＯＵＳＦＬＯＷＦＩＥＬＤ（ＴＶＦＦ）ＩＮＳＩＤＥＴＵＲＢＩＮＥＣＡＳＣＡＤＥＧｕｏＹａｎｈｕ，ＳｈｅｎＭｅｎｇｙｕ，ＷａｎｇＢａｏｇｕｏ（Ｄｅ...