部分平均Navier-Stokes方法模拟超声速斜面空腔流动

采用部分平均Navier-Stokes (PANS)方法对超声速斜面空腔流动进行了数值模拟研究,评估其对于超声速非定常湍流模拟的性能,并与雷诺平均N-S(RANS),脱体涡模拟(DES)的计算结果及实验数据进行对比.研究表明:模化湍动能比例全场可变的PANS方法预测的速度、壁面压力和摩擦力系数分布与DES的结果非常接近,都与实验值吻合得较好,优于RANS的计算结果;非定常雷诺平均N-S(URANS)计算的自由剪切层近乎为二维稳态的,而DES和PANS方法可以求解出更为丰富的流动结构;模化湍动能比例全场统一的PANS方法虽然可以求解出相比RANS更多尺度的流动结构,但在近壁区不能回归到RANS模型,预测的湍流边界层的速度型偏离对数律,后续的流动计算也偏离实验数据.      

高纬极区离子速度分布函数多项式解及非相干散射谱计算

为更准确描述高纬极区电离层离子分布函数, 分别采用弛豫碰撞模型和麦克斯韦分子碰撞模型描述玻耳兹曼方程的碰撞项, 通过求解两种模型下的输运方程, 分别得到两种模型基于麦克斯韦分布下离子分布函数的13矩近似和基于双麦克斯韦分布下离子分布函数的16矩近似. 进一步根据Sheffield理论, 利用两种模型下离子分布函数的13矩和16矩近似, 计算了非相干散射谱, 并对结果进行对比分析. 结果表明, 相对于弛豫碰撞模型, 麦克斯韦分子碰撞模型能更好地描述电离层E层中离子与中性成分的相互作用. 相对于离子分布函数的13矩近似, 16矩近似更适合描述由于电场增加导致的离子温度各向异性特征.      

抛物化稳定性方程在曲面边界层中的应用

为了研究抛物化稳定性方程在曲面边界层中的应用,选取了等曲率圆环管道、等曲率板和NACA0012翼型3种典型的曲面边界层.通过计算小幅值扰动波的演化,抛物化稳定性方程的计算结果和线性稳定性理论、数值模拟扰动方程结果相符,说明可以使用抛物化稳定性方程研究曲面边界层的小幅值扰动波的演化和稳定性分析;通过计算有限幅值扰动波的演化,线性阶段抛物化稳定性方程计算结果和扰动方程相符,在非线性很强的阶段,抛物化稳定性方程计算发散,发散的位置可作为转捩位置.      

基于共形几何代数的空间并联机构位置正解

将共形几何代数(CGA)引入空间并联机构位置正解中,提出了一种空间3-RPS并联机构位置正解新算法。以任意一条支链轴线与静平台平面的夹角为待求变量,基于点的CGA表达方法建立了该支链与动平台连接的铰接点关于待求变量的数学表达式;通过2次构造2个空间球和1个平面的外积,分别获得动平台其余2个铰接点的点对;利用距离公式,只需简单的平方运算可直接推导出该问题关于待求变量的一元16次输入输出方程,进而获得了该机构的全部16组解析解,无增无漏。该方法没有繁琐的坐标变换和矩阵计算,以及复杂的多元高次非线性方程组消元求解。通过数字实例计算表明,求解过程较清晰地揭示出机构运动的几何特点,几何直观性好。     

CMGs驱动空间机械臂的自适应终端滑模控制

摘要： 针对V构型控制力矩陀螺（CMGs）驱动的空间机械臂的轨迹跟踪控制问题,研究一种自适应非奇异终端滑模（ANTSM）控制方法.利用基于Kane方程的递推组集算法建立了系统的动力学模型.以跟踪误差为变量,构造非奇异滑动面,以保证跟踪误差在滑动面上有限时间收敛.针对系统质量特性参数与关节处干扰力矩的不确定性,设计自适应控制器用以调节控制增益.该控制方法无需不确定性的上界,且闭环系统具有最终一致有界性.仿真结果表明,该控制器可使系统准确跟踪期望轨迹,并对质量特性参数不确定性和关节干扰力矩具有良好的鲁棒性.     

基于修正N-S方程本构关系的气体动理论耦合方法

随着稀薄程度的增加,Navier-Stokes方程的线性本构关系难以正确描述稀薄气体输运特性,高阶非线性本构关系往往数学形式极为复杂,对数值求解造成稳定性差等问题。为了发展适宜于近空间飞行器气动特性分析的高超声速稀薄流动模拟方法,本文利用求解Boltzmann模型方程的气体动理论统一算法(Gas Kineitc Unified Algorithm,GKUA)对应力张量、热流等宏观量数值积分求解的优势,提出了一种基于数值修正N-S方程本构关系的气体动理论耦合方法。通过将GKUA获得的应力张量及热流用于修正N-S方程的本构关系,实现了存在局部稀薄效应的流动模拟,并且通过可压缩平板边界层、圆柱绕流问题的数值模拟,验证了方法的有效性。     

共轴刚性旋翼的悬停气动性能和流场干扰

建立了一种基于RANS (雷诺平均Navier Stokes)方程的共轴刚性旋翼悬停流场数值模拟方法。使用旋翼亚声速和跨声速悬停实验结果,验证了该方法的准确性。对刚性旋翼XH 59A流场模拟表明:气动特性与飞行实验结果比较吻合,共轴旋翼特性优于同实度单旋翼,具备更高的悬停效率。对比半实度单旋翼气动性能,受下旋翼对流动的抽吸影响,上旋翼性能略有下降;下旋翼性能下降更甚,主要是因为处于上旋翼下洗流中,其有效迎角减小。数值方法获得该旋翼最高悬停效率为67%,总距角为14°,14°。对比研究表明,共轴刚性旋翼较常规共轴旋翼极间距小,悬停性能更高。      

叶轮式通风器分离效率计算方法

为研究叶轮式通风器分离特性并建立满足工程设计要求的分离效率计算方法,在对油滴进行了受力分析后,推导了油滴运动轨迹模型,结合油滴捕获条件,给出了不同工况条件下油滴的最小分离直径,进而实现了叶轮式通风器分离效率的理论计算。研究表明:在不同转速、通风量和通风器入口温度时,分离效率和最小分离直径计算值与实验值趋势一致;分离效率最大误差分别为0.95%、1.01%和1.01%;最小分离直径最大误差分别为1.33%、2.29%和6.20%。计算结果与实验结果吻合较好,该方法可为叶轮式通风器设计及分离性能优化提供参考。