一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

CE／SE并行计算研究及在Euler方程中的应用

时空守恒元和解元方法(CE/SE)是一种高精度的数值求解方法.采用改进的CE/SE方法和并行分区技术,针对非结构网格,发展了一套求解二维Euler方程的并行程序.对NACA0012翼型的亚跨声速流动和多段翼型复杂绕流进行了分区并行求解,多区并行计算的结果与实验的结果吻合较好,表明了本文研究方法的正确性和有效性.同时文中的分区方法能实现各处理器之间的负载均衡,有效地节省计算时间.     

双框架控制力矩陀螺群的建模与分析

 研究采用双框架控制力矩陀螺群在卫星高精度姿态控制中的建模问题。对于中心刚体带有多个双框架控制力矩陀螺,构造了该系统的拉格朗日函数。采用拟拉格朗日方程建立了卫星姿态动力学模型。在建模中,陀螺的内、外框架未进行无惯量的假设,使得双框架运动带来的动力学特性得到了准确的反映;设计了适宜的标记,使得方程形式简洁并且统一。该模型为高精度的数学仿真和分析提供了基础。基于系统的拉格朗日函数,应用拉格朗日方程,分别建立了沿内、外框架轴和转子轴的轴向动力学方程,为执行机构伺服系统的动力学特性分析和控制器设计提供模型的依据。     

超燃冲压发动机内流性能的一维评估

在准一维Eu ler方程的基础上,通过增加反映面积变化、摩擦、添质和化学反应放热的源项,发展了超燃冲压发动机推进流道计算的一维软件,可以在很短的时间内提供发动机的推力和摩擦力,同时给出发动机参数沿轴向的变化和出口值。通过日本国家航空与航天实验室(NAL)的双模态模型氢反应流场计算验证了软件的可靠性后,针对煤油发动机,分析了隔离段和燃烧室扩张角、燃料喷射位置与方式、化学反应放热项的两种不同求解方法、隔离段高度对燃烧室性能的影响,为飞行器设计配备了能稳定激波且性能较好的发动机构型。     

非线性四阶椭圆边值问题解的存在性与唯一性

分别利用全局同胚理论和动力系统理论的一些结论,研究了非线性四阶椭圆边值问题解的存在性与唯一性。得到的结果改进和推广了非线性四阶椭圆边值问题的相应结果。     

改善航空发动机特性计算收敛性的方法

详细分析了航空涡轮发动机部件共同工作的非线性方程组的求解收敛性问题,统计出不收敛的类型并分析了不收敛的机理.针对常见的不收敛问题,以牛顿迭代法为基础,分别提出了以发动机部件特性图为基础的独立变量值限制法、独立变量初值拟合法、变步长牛顿法以及部件特性扩展法等解决措施.特性计算程序中采取这些改进的方法后,对某型定几何单轴涡喷发动机和某型尾喷管喉部面积可调的变几何双轴涡扇发动机特性计算结果表明,发动机特性计算的收敛性和收敛速度得到大幅度提高.      

T型尾翼民机深失速气动特性的计算研究

以非结构混合网格上的雷诺平均N-S方程(RANS)求解方法为基础,对采用T型尾翼的某民机外形在大迎角深失速情况下的粘性绕流进行了数值模拟,获得了该机在巡航和降落外形在深失速情况下的气动特性变化规律.通过对T型尾翼深失速情况下的流动特点进行分析,解释了T型尾翼布局飞机在深失速情况下俯仰力矩系数会随迎角出现不利于纵向稳定性的变化趋势的原因.     

直升机旋翼桨叶翼型积冰的数值模拟

本文对直升机旋翼桨叶翼型积冰进行了数值模拟.建立了空气、水滴两相流控制方程,通过求解气体欧拉方程得到空气流场,求解水滴控制方程得到翼型表面的水滴收集率.根据结冰表面的对流换热特性,使用积冰数值模拟的热力学模型确定翼型表面的积冰量,利用冰层时间推进法模拟积冰过程,采用积冰法向生长模型预测积冰形状;对NACA0012翼型在不同环境温度下的积冰进行了预测,最后对预测结果进行了分析.     

基于N-S方程串并行计算的机翼优化设计

对典型的大展弦比和中等展弦比机翼在0.7~0.9马赫数区间进行了串行和并行计算的数值校验,证实了以N-S (Navier-Stokes)方程为主控方程的串行和并行流场求解器的正确性,并讨论了并行效率和加速比.结合Powell算法,讨论了在确定的机翼平面形状和翼型的条件下,以升阻比最大为目标的三维机翼截面翼型最大厚度与扭角的优化设计.算例结果表明,厚度的非线性分布和负的扭角会改善机翼流场的流动状态,使机翼的升阻比得到提高,优化设计方法是可行的.