基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

基于Fibonacci搜索方法的航空发动机性能寻优

介绍了Fibonacci搜索方法应用于发动机性能优化中的原理,建立带两个出口面积、气流马赫数可变的引射器的涡喷发动机的总体性能计算模型,应用Fibonacci搜索方法在优化参数区间内对该发动机进行了最大推力、最低油耗、最低涡轮前温度模式的寻优.举例说明了Fibonacci搜索方法应用的具体过程.给出了不同高度、马赫数下的寻优结果.结果表明:最优目标性能满足的约束边界与飞行高度、马赫数密切相关.不同的飞行高度、马赫数对应了不同的最优控制策略.但不论飞行任务如何,根据该方法进行的优化控制所获得的目标性能是最佳的.该优化方法不仅适用于该类涡喷发动机,也为其他类型航空发动机的优化控制打下一定的基础.      

Dynamic modeling and analysis of vortex filament motion using a novel curve-fitting method

Applications of a novel curve-fitting technique are presented to efficiently predict the motion of the vortex filament, which is trailed from a rigid body such as wings and rotors. The governing equations of the motion, when a Lagrangian approach with the present curve-fitting method is applied, can be transformed into an easily solvable form of the system of nonlinear ordinary differential equations. The applicability of Be′zier curves, B-spline, and Lagrange interpolating polynomials is investigated. Local Lagrange interpolating polynomials with a shift operator are proposed as the best selection for applications, since it provides superior system characteristics with minimum computing time, compared to other methods. In addition, the Gauss quadrature formula with local refinement strategy has been developed for an accurate prediction of the induced velocity computed with the line integration of the Biot–Savart law. Rotary-wing problems including a vortex ring problem are analyzed to show the efficiency, accuracy, and flexibility in the applications of the proposed method.     

基于格子Boltzmann方法表面形貌对微通道对流换热的影响

为推动微通道冷却技术在航空发动机中的应用,进一步提高航空发动机性能,本文通过扫描电子显微镜、粗糙度扫描仪等多方式扫描得到实际加工物体表面形貌,结合分形理论进行处理,得到真实、光滑、分形插值、康托尔集四类圆形截面管模型,管径100μm~400μm;采用格子Boltzmann方法进行数值模拟,首先通过与文献中实验结果的对比,验证了此方法的正确性,然后对空气在这四类圆形截面微通道中的对流换热特性进行了数值模拟,计算雷诺数Re范围80~640;结果发现：表面形貌是微通道中对流换热不可忽略的因素,相对粗糙度越大越有利于换热,真实扫描管的换热性能要比相同条件下光滑管高2.22%;分形理论可用于微通道表面形貌的构建,在相同迭代次数下,分形插值模型与真实扫描形貌换热性能相差0.03%,而康托尔集模型与真实扫描形貌换热性能相差1.88%,分形插值模型较康托尔集模型更能真实地反映实际形貌;微通道中,随着雷诺数的增加Nu将不再是常数,而是有增加的趋势。     

10cm氙离子推力器变推力特性研究

非重力阻尼的连续、快速、高精度补偿是实现重力梯度测量卫星精细重力场测量的关键技术之一,直接影响到整星工程任务的成败。针对重力梯度测量卫星在轨飞行期间对电推进系统宽范围连续变推力能力的应用需求,分析了10cm氙离子推力器推力调节响应特性。在此基础上,通过对阳极电流、励磁电流和阳极流率等推力高敏感响应参量的组合调节,开展了推力调节试验研究,验证了10cm氙离子推力器宽范围连续变推力调节能力,获得了1～20mN范围内的推力调节性能及其变化规律。试验结果表明：在采用地面供电、供气设备条件下,10cm氙离子推力器能够在100～597W的功率范围内实现0.98～20.29mN的推力宽范围调节,比冲175～3500s,推力分辨率优于50μN。研究为建立10cm氙离子电推进系统的推力控制数学模型及调节控制算法奠定基础。     

基于模型的装备维修性设计控制技术

伴随数字化设计方法日趋成熟,飞机装备维修性设计面临如何与三维数字化设计环境下功能性能设计同步的难题。对比分析了数字化环境维修性设计与传统维修性设计的主要差别,提出了基于飞机装备数字化设计环境维修性设计流程,研究了维修性设计参数的表达模型以及维修性参数信息在三维设计数字环境中的转化、集成与利用方式。基于图论和网络节点框图,提出了一种以节点网络图、单节点维修时间模型、维修过程时间模型为表征的飞机装备维修性设计与控制方法。经实例验证,提出的这套基于数字模型的维修性设计技术,对装备在数字化设计环境下同步开展维修性设计具有一定的工程指导价值。     

基于断口分析的钛合金轮内部缺陷损伤容限

为进行某轮损伤容限设计,开展了裂纹扩展断口分析和仿真分析研究。由断口分析可知:疲劳源为一处内部自然缺陷;依据疲劳辉纹确定了裂纹扩展速率;在裂纹长度为2 mm附近,裂纹扩展速率明显增大,为第一、第二加载阶段转换区域;裂纹稳定扩展区裂纹长度与裂纹扩展速率呈双对数线性关系;应用列表梯度法和Paris公式法反推了第二加载阶段的疲劳寿命,与该阶段实际循环次数的最大误差是163%。裂纹稳定扩展阶段裂纹扩展仿真值与断口反推值吻合;非稳定扩展阶段仿真值与断口反推值的最大误差为-215%;基于以上研究,合理确定了某离心轮内部裂纹表面扩展停机检测周期。该类轮非稳定、失稳扩展阶段寿命占内部裂纹表面扩展阶段寿命的比例达248%~357%,因此准确计算具有重要意义。     

0.3 m低温风洞液氮供给系统研制

基于系统级的一维热流体模拟分析,建立了适用于研究分析0.3 m低温风洞液氮供给系统的数学模型,并开展了系统漏热、两相流及缓冲罐中液氮容积等流体动力学分析;在系统现有控制策略及试验数据的基础上,基于该数学模型开展了系统压力动态响应分析,获得了在阀门动态调节过程中管网压力的瞬态响应,计算结果与试验值的总体误差控制在10%以内。喷射压力一致化改造避免了阀间干扰,添加的回流管道消除了供给末端的两相流现象,使喷射压力控制精度达到1.1%,调节时间减少到23 s,实现了风洞总温快速安全调节和精确控制。     

退化量有缺失的无失效小样本轴承可靠性评估

针对有试验性能退化量参数记录、且(或)有部分数据缺失情况的小样本无失效轴承试验问题,通过由Taylor和Thompson提出的数据模拟法实现补全样本退化量,结合Bootstrap自助法扩大样本量,再根据基于性能退化轨迹的补充信息方法来进行其可靠性评估。选取7组受试轴承的振动退化量,对比在完整数据和带有缺失数据情况下的分析结果,发现可合理利用原舍弃不合规试验的部分有效信息,使之增加可靠性评估的样本数,从而得到较不用这些数据更为准确的结果,且所得结果较用完整数据结果绝对值相差在01以内。对比由极大似然估计法和加权E-Bayes法分析试验寿命数据的结果,发现该方法所得评估结果更优,与试验实际相比误差在10%以内,对于提高评估精度及降低试验成本有积极的实际意义。     

Fluid structure interaction simulation of supersonic parachute inflation by an interface tracking method

An Arbitrary Lagrangian–Eulerian (ALE) approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation. A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished. In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact, a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces. In addition, the breakpoint due to the fluid mesh with negative volume is losslessly restarted by the conservative interpolation method. Based on this method, fluid and structural dynamic behaviors of a highly folded disk-gap-band parachute are obtained. Numerical results such as maximum Root Mean Square (RMS) drag, general canopy shape and the smallest canopy projected areas in the terminal descent state are in accordance with the wind tunnel test results. This analysis reveals the inflation law of the disk-gap-band parachute and provides a new numerical method for supersonic parachute design.