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

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

海陆交界非平稳场景星载ScanSAR扇贝效应抑制方法

扫描式合成孔径雷达(ScanSAR)因不受光线、天气影响,能够提供宽测绘带图像而被广泛应用在海洋监测领域。但ScanSAR图像的扇贝效应(Scalloping)和子带拼接效应(Inter-Scan Banding,ISB)严重地降低图像质量,尤其在海陆相接一类的非平稳场景下,干扰更为严重,影响图像的应用效能。针对上述问题,文章给出了一种基于图像分层和卡尔曼滤波的ScanSAR图像扇贝效应和子带拼接效应的抑制方法。针对非平稳场景,先进行图像分层处理使其满足线性卡尔曼滤波需求;然后,利用加性二维周期函数对子图层的扇贝效应和子带拼接效应进行建模;再将卡尔曼滤波器应用于条纹抑制;最后,对真实ScanSAR图像数据的有效处理结果验证了算法的有效性。     

转子结构系统界面失效分析及稳健设计方法

针对结构非连续转子中连接界面接触状态变化引起界面力学特征改变,使转子系统动力特性不稳健的问题进行分析,并提出了一种界面接触状态的稳健设计方法。结果表明:该方法通过分析结构非连续转子系统中界面损伤失效的力学过程,对结构特征参数进行优化设计,降低界面力学特性对载荷环境的敏感度。以涡轮转子套齿连接为例进行优化设计得到,在加工装配误差所导致的应力极限情况下,最优解的接触应力水平更加远离应力约束边界,使带有连接界面的结构非连续转子系统的力学特性更加稳健,具有较好的工程应用价值。      

典型进气函数下转静盘腔瞬态响应特性研究

为研究航空发动机转静盘腔非稳态流动的瞬态响应特性,建立了转静盘腔CFD计算模型,采用基于Fluent用户自定义函数(UDF)编程的非稳态数值计算方法,分别研究了进口压力以阶跃函数、斜坡函数和正弦函数变化时转静盘腔的瞬态响应特性,并提出了相应的瞬态响应特性评价指标。结果表明:不同进气函数下的转静盘腔瞬态响应特性存在明显差异,同一进气函数下盘腔中各流动参数的瞬态响应存在不同程度的滞后。阶跃函数的无量纲跃升幅值由1.025变化到1.3时,盘腔平均总压的响应时间增加了3369.2%;斜坡函数的斜坡时间由0.01s增加到0.16s时,盘腔平均总压的响应时间增加了163.9%;而正弦函数的角速度由22rad/s变化到55rad/s时,盘腔平均总压的响应时间缩短了45.6%。不同进气函数下转静盘腔平均总温随时间的变化曲线均呈现出明显的过冲现象,其峰值时间和超调量变化规律与进气函数类型密切相关。     

锥导乘波体构型的气动特性不确定度分析

为研究锥导乘波体偏离设计条件下气动特性变化情况,采用稀疏的非嵌入式混沌多项式方法,对乘波体气动特性进行了不确定性量化及全局非线性灵敏度分析。首先,采用CATIA二次开发技术对锥导乘波体进行参数化建模;其次,在来流速度、温度、密度和迎角满足特定扰动的条件下,通过拉丁超立方试验设计生成样本,并采用CFD进行计算;最后,根据试验设计样本建立响应面,通过混沌多项式分析得到了乘波体气动力系数的不确定度。灵敏度分析结果表明,迎角在锥导乘波体的气动特性变化中起主导作用。对马赫数和压强的流场不确定性分析结果表明,气动特性变化主因是乘波体前缘处的压力泄漏,影响了上表面压力分布,导致了气动性能的改变。     

微型离心压气机热效应评估模型及应用

基于微型离心叶轮非绝热边界假设,通过对其内部实际换热过程的详细分析,合理将叶盘、叶片换热等效为环肋、直肋换热,重构了微型离心叶轮的换热过程;根据导热微分方程和肋片散热方程,对叶轮内部热传导及对流换热过程控制方程进行了适应性修正;采用全三维数值模拟与模型预估结果进行对比,结果表明:较之零维热网络模型,该模型能够将预估精度至少提高4%;原有只在三维仿真阶段考虑热边界影响的非绝热压气机设计方法相比,结合了该模型的设计方法,可将各自最佳效率点压比和效率分别提高11%和30%,同时,设计周期降为原有方法的14.3%。      

凸约束非凸二次规划问题的分枝定界方法

针对凸约束非凸二次规划问题,给出了一个分枝定界方法。首先,我们构造一个多胞体包含可行域,然后根据凸集上非凸二次规划问题的整体最优解在可行域边界达到的性质,对锥所包含的可行域的边界构造一个包含它的超矩形体,并对这个超矩形体构造一个外接球。我们通过求解球约束非凸二次规划问题的整体最优解来确定下界,并把锥的棱与可行域的边界交点的目标函数值的最小值作为上界,把锥剖分技术与外逼近方法结合起来寻找原问题的整体最优解。最后,我们对这个方法进行收敛性分析。     

非均匀栅距对压气机叶片非定常气动力的影响

研究了轴流压气机叶片采用非均匀栅距对叶片非定常力影响的问题,利用高分辨率格式求解准三维雷诺平均N-S方程。非定常流场计算采用了双时间步法与LU-SGS隐式解法耦合的方法。研究了带进口导叶的单级三排叶栅非定常流动,其中动叶栅距不变,而导叶和静叶采用非均匀栅距。对不同结构的计算分析表明,在导叶和静叶采用一定的非均匀栅距组合型式后,在保证压气机级的性能变化不大的同时,动叶表面所受非定常气动力的状况得到改善,这将有助于提高动叶片的疲劳寿命。     

动态负载平衡的二维非结构网格DSMC并行算法研究

基于PC-CLUSTER群机并行体系结构与消息传递库MPI并行环境,研究了二维非结构网格DSMC并行算法。提出一类非结构网格动态分区策略,保证各子区域的分子数量大致相等,实现计算进程间的动态负载平衡。利用MPI库函数构造了符合DSMC并行原理的单步通讯法。采用单控制多数据流(SPMD)以及Master/Slave并行模式,设计了二维非结构网格DSMC整体并行算法。在程序的编制过程中,充分展现了Fortran90高级语言的主要特性,引人动态数组、指针、链表及派生类型数据。最后对过渡流域高超声速绕流进行了并行计算,数值试验的结果表明本文设计的并行算法可以取得较为理想的加速比。     

基于非结构嵌套网格方法的旋翼地面效应数值模拟

 建立了一个基于非结构嵌套网格的流场求解器,用来精确模拟复杂的旋翼近地流场,为更好地分析地面效应(IGE)对旋翼气动特性的影响提供一套计算方法。在该求解器中,控制方程采用惯性坐标系下的非定常N-S方程,空间方向上采用二阶迎风格式,并用背景网格的一个面模拟地面作用,以方便地面边界条件的处理。应用所建立的模型,首先针对有实验结果可供对比的旋翼无地面效应（OGE）流场进行了数值模拟,以验证计算方法。然后着重计算了IGE下流场中的桨尖涡空间位置和旋翼拉力增益,并对旋翼在小速度前飞状态下的近地流场及地面涡形成过程进行了模拟和分析。在此基础上,得出了一些有意义的结论。