二维流热耦合数值模拟程序及方法研究

本文发展了一套二维流热耦合数值模拟程序,流体域中采用结构化网格,固体域中采用非结构化网格,流体域和固体域耦合面上的网格为非匹配连接。为了使该程序能够同时模拟低速不可压流动,程序中还加入了预处理方法。此外,本文还发展了一种能较好保持通量守恒的界面插值技术,来处理流体域和固体域耦合面上的非匹配网格。数值计算结果与实验值或商业软件结果的对比表明:本程序中的预处理方法有效地改善了低速不可压流动下的计算收敛性和收敛精度,使本程序能够同时模拟高速和低速不可压流动;扩压器和叶栅的流热耦合计算表明本文程序在流热耦合计算方面具有较好的模拟能力,同时也表明本文发展的界面插值技术的可靠性。     

基于非结构动网格的跨声速非定常流场求解

对用于非结构动网格生成的弹簧近似方法进行了研究.通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入防扭转、防挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进.采用改进后的方法并通过耦合求解了基于ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程,模拟了做谐和振动的刚性无后掠机翼及某翼身组合体的跨声速绕流.结果表明,改进后的方法可以大大提高网格变形能力和网格质量,且计算结果与实验结果吻合良好.     

N-S方程数值模拟方法研究及边界层中相干结构的非线性演化

建立了一套三维不可压流动的直接数值模拟方法。该算法在x及y向构造了基于非等间距网格的紧致有限差分格式和非线性项的迎风紧致型格式,在z向采用Fourier谱方法。并研究、提出了三维、非定常流体运动下游边界的无反射出流条件。另外,本文还构建了多涡结构的初始扰动理论模型,利用直接数值模拟的方法研究了边界层中多涡结构的非线性演化特点。结果表明多涡结构非线性演化中的许多特点,如雷诺应力分布、高剪切层的形成、喷射和扫掠等与壁湍流中相干结构的发展规律和现象十分相似。同时,通过对边界层中多涡结构的非线性演化问题的直接数值模拟,验证了该算法具有计算精度高,稳定性好,收敛速度快,出流边界影响小等优点。     

主动对等计算系统的路由机制

针对无结构对等计算系统路由机制中用户需求粒度过细,忽略需求模糊性和多元性的问题,提出了基于用户需求的路由机制.借鉴模糊理论,提出了从文件类型、文件质量、下载速度多方面综合评价节点能力的路由模型.给出了节点能力因素及节点能力的定义,采用隶属度的方式来描述节点提供资源的能力与用户需求的符合程度.分析和实验结果表明,该机制在满足用户需求和减少网络开销方面能够达到很好的均衡.      

Jameson有限体积法对非结构网格推广的改进

本文将Ｄｅｌａｕｎａｙ法用于复杂物形的非线性结构网格生成,并利用黄明格提出的布点原理避免边界穿透。绕流的计算采用Ｅｕｌｅｒ方程和Ｊａｍｅｓｏｎ的有限体积法。本文对三维问题提出新的激波敏感器,对方法做了改进,并用于ＡＧＡＲＤ Ｍ６机冀、导弹模型和飞机模型的绕流计算。     

非结构网格上Euler方程的区域分裂算法及并行计算

改进了“波阵面”区域分裂算法,并应用于流场区域的划分;对于子区域边界的不光滑现象,为尽量减少通讯消耗,提出了一种边界并行优化策略。利用PVM并行环境,探讨了非结构网格上求解Euler方程的分区并行算法。根据改进的区域分裂算法及优化策略,运用Jameson有限体积法,对二维翼型流场进行了分区并行求解,多区计算的结果与单区计算的结果作了比较,表明了本文研究方法的有效性。     

分级型进化算法及其在翼型跨声速气动优化中的应用

提出了变精度模型的分级型进化算法并对初始翼型为NACA0012的二维翼型进行了跨声速流动下的形状增升优化。借鉴自然进化中不同进化阶段个体适应度评估环境不同的机理，构造了分级精度模型以加速优化过程。对翼型进行了给定气动条件下的形状增升优化，给出了优化结果，并与传统基因算法及单精度模型算法结果进行了对比。     

基于非结构重叠网格的螺旋桨滑流非定常数值模拟

基于各向异性非结构混合网格及多套网格重叠技术,通过求解非定常雷诺平均Navier-Stokes(URANS)方程,分别研究了螺旋桨动力特性和螺旋桨滑流对涡桨飞机气动特性影响。运用物面相交准则实现多套网格间的挖洞,采用距离权和三线性插值技术传递重叠区变量信息,并针对定轴旋转优化了网格重叠边界。采用双时间步求解控制方程,内迭代计算采用LU-SGS(lower-upper symmetric Gauss-Seidel)方法。模拟了单独螺旋桨旋转的非定常流动,计算的拉力系数和扭矩系数与试验结果一致,表明采用的数值方法和网格技术能有效模拟螺旋桨滑流效应。数值模拟了某涡桨飞机,对比了有无滑流下飞机的气动特性,给出了滑流对飞机影响。结果表明:滑流导致其后部气流加速和旋转,改变飞机气动特性,增大飞机的阻力和俯仰力矩,导致气流下洗,影响机翼及平尾压力分布。      

飞机外挂物离机过渡过程的有限元算法.

研究飞机外挂物在按给定轨迹下沉旋转的离机过渡过程中，流场非结构网格的调整和再生成方法以及过渡过程欧拉方程有限元数值解法。中指出外挂物运动时流场网格调整或再生成仅需在围绕在外挂物的一小区内进行以减少计算量，而且固定网格的欧拉方程的有限元数值解法也可适用于运动网格情况，只要将固定网格的守恒通量项Fk用运动网格的Ek来代替，并建立了Fk和Ek之间的关系。末给出了带外挂物飞机的实例计算结果。     

Numerical Simulation of Unsteady Flow Around Forward Flight Helicopter with Coaxial Rotors

Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of the two coaxial rotors both become worse because of the aerodynamic interaction between them, and the influence of the top rotor on the bottom rotor is greater than that of the bottom rotor on the top rotor. The downwash velocity at the bottom rotor plane is much larger than that at the top rotor plane, and the downwash velocity at the top rotor plane is a little larger than that at an individual rotor plane. The downwash velocity and thrust coefficient both become larger when the collective angle of blades is added. When the spacing between the two coaxial rotors increases, the thrust coefficient of the top rotor increases, but the total thrust coefficient reduces a little, because the decrease of the bottom rotor thrust coefficient is larger than the increase of the top rotor thrust coefficient.