基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

基于开路电压法的电池荷电状态估计误差校正

本文提出分别研究电池开路电压在充电结束后停置过程和放电结束后停置过程的不同变化规律,并考虑了不同停置时间在两个过程中的作用和影响,从而达到提高开路电压法估算荷电状态精度的效果。实验验证表明,此方法的精度比传统开路电压法更加准确。     

Multiscale topology optimization using feature-driven method

This paper presents a multiscale design method for simultaneous topology optimization of both macrostructures and microstructures. Geometric features are extended as design primitives at both macro and micro scales and represented by Level Set Functions (LSFs). Parameters related to the locations, sizes, and orientations of macro and micro features are considered as design variables and optimized simultaneously. In the overlapping areas of different macro features, embedded microstructures are optimally figured out as the solution of the corresponding sub-optimization problem. In this study, the eXtended Finite Element Method (XFEM) is implemented for structural and sensitivity analyses with respect to design variables. This method has the advantage of using a fixed grid independent of the topology optimization process. The homogenization procedure is applied to calculate the effective properties of considered microstructures in each macro feature. Numerical examples are presented to illustrate the effectiveness of the proposed method. Results depict that the multiscale design cannot obviously improve structural stiffness compared with a solid-material design under the linear elastic condition.     

机匣处理作用下高负荷轴流压气机转/静子匹配设计研究

为了探索机匣处理作用下转/静子的轴向匹配方法以进一步提高压气机级的失速裕度,研究了静子的叶型安装角及"弯"、"掠"规律对压气机性能的影响,针对机匣处理与优化静子的组合结构进行了非定常数值模拟,阐述了该结构的扩稳机理以及压气机新的失速机制。研究结果表明,在机匣处理作用下,静子成为压气机失速的触发因素,通过对静子叶型安装角及"弯"、"掠"规律的优化均可进一步提高压气机级的失速裕度,其中改变静子"弯"型对压气机级失速裕度的改善最大。组合应用机匣处理与尖部反弯根部正弯静子后,压气机效率基本不变,失速裕度提升了80.2%,较单独使用机匣处理提升30.9%。在该组合结构作用下,压气机的失速由静子触发,静子叶根吸力面在激波作用下发生附面层分离,且与轮毂表面附面层相互作用形成角区涡,接近失速边界时,静子叶根形成"前缘溢流,尾缘反流"现象,造成静子通道的大范围堵塞,诱发压气机失速。压气机级的扩稳应充分考虑机匣处理的影响,对静子进行优化设计。     

加筋壁板大开口结构稳定性分析及试验验证

为了研究加筋壁板大开口结构的压缩承载极限问题,规划了两组具有不同设计参数的试验件进行了静力试验,并同时采用了线性特征值法以及非线性的弧长法对其压缩极限进行了对比分析。在采用弧长法分析时,引入了模态扰动技术,而且考虑了材料的弹塑性性质。计算结果与试验的对比表明:针对复杂结构稳定性问题,屈曲特征值法具有较大局限性,而结合模态扰动的弧长法对这种复杂的加筋壁板类结构压缩极限的分析精度较好。通过分析给出了加筋壁板大开口结构设计的一些建议。     

大气环境对直升机旋翼桨-涡干扰噪声辐射特性的影响

建立了一个适用于旋翼桨-涡干扰气动载荷计算的计算流体力学(Computational fluid dynamics,CFD)/自由尾迹耦合模型,为提高计算效率,提出了一种高效的耦合策略进行不同计算域间的信息交换策略。在此基础上,结合基于声类比法的FW-H方程构建了旋翼桨-涡干扰噪声的计算方法。应用所建立的方法,以OLS(Operational loads survey)旋翼为研究对象,深入分析了大气环境对旋翼噪声辐射特性的影响。研究发现:随着飞行高度的增加,旋翼噪声辐射特性发生了明显的改变,逐渐由桨盘前行侧转变为指向桨盘前方,噪声幅值先增大后减小。文中从桨-涡干扰距离、干扰位置变化角度计算分析了大气压力、音速及空气密度等环境参数对旋翼桨-涡干扰噪声辐射特性的影响,并得出了一些有实际意义的影响规律。     

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构，层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟，研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型，基于结构化拼接网格，对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性，雷诺数对翼型气动特性具有显著影响。随雷诺数降低，时均分离泡外形增大、位置后移，平均阻力系数增大，特别是在较低雷诺数下，翼型升阻力系数随时间出现振荡现象。进一步研究表明，造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低，流动黏性增大，导致分离剪切层速度梯度减小，流动发生转捩及再附位置后移，直至翼型表面不再发生转捩和再附。     

压电阻抗技术监测固体推进剂老化研究

为了进一步研究表面粘贴式的主动传感技术——压电阻抗技术(EMI, Electro-Mechanical Impedance)监测固体推进剂的老化,构建了基于线粘弹杆结构的一维机电耦合模型,根据粘弹性波的理论与边界条件进行分析,建立一种线粘弹结构动态模量与其机械阻抗的表征关系,并对构建的机电耦合阻抗模型进行了数值计算及试验验证。针对HTPB固体推进剂开展高温热加速老化试验及压电主动激励试验,根据监测所得导纳频谱并通过提取结构机械阻抗进行分析,结果表明:压电陶瓷片与推进剂结构在高频机电耦合运动时(300kHz,280kHz)有明显的共振现象,结构固有频率高于机电耦合共振频率,并且在200kHz～400kHz频段与700kHz～900kHz频段内,结构机械阻抗频谱峰值会因固体推进剂的热老化时间增长而降低,且与热老化时间之间满足线性关系。由此可见,压电阻抗法能够与固体推进剂在力学性能上建立表征关系,通过得到固体推进剂结构机械阻抗能够监测其老化损伤。     

插入式机翼下壁板对接附加弯矩研究

插入式机翼下壁板对接具有双剪传力稳定、疲劳性能好的优点，但其结构中心线在对接区变化明显，会带来附加弯矩。为尽量减小对接区的附加弯矩，提出了在建立飞机骨架模型时即优化中央翼下翼面外形面相对外翼下翼面的位置方法。基于插入式机翼下壁板对接结构的特点，阐述了对接结构偏心的来源和附加弯矩的形成；针对某A型飞机的对接结构计算了偏心值，并利用力法对附加弯矩在对接区的分布进行了计算分析。以某A型飞机的对接结构为基础，建立了4组插入式下壁板对接结构的模型，每组模型的中央翼下翼面位置相对外翼的不同；分别用力法和有限元法对附加弯矩进行了计算。结果表明：可以通过优化中央翼下翼面外形面的相对位置达到减小对接区附加弯矩的目的。描述了另外两种下壁板对接形式的附加弯矩情况，并和插入式的进行了简单比较。最后，总结了为减小区域附加弯矩及其不利影响在对接结构设计上需要注意的点。     

三维叶轮机叶片黏性反问题设计改进方法

提出了一种改进的适用于三维黏性流场的叶轮机叶片反问题设计方法。该方法假设叶片的中弧线具有虚拟移动速度,其位移量由目标载荷与实际载荷的差值计算,并利用黏性底层厚度进行限制。采用三次B样条曲线插值方法对叶片中弧线进行光顺,新叶型通过更新后的中弧线和给定的叶片厚度得到。算例验算结果表明:该方法能够根据设计意图对叶型进行修改,鲁棒性强,收敛速度快,叶片的可变自由度高,不依赖于特定的网格和求解器,具备一定的通用性。