TRIP2.0软件的确认：DPWⅡ复杂组合体的数值模拟

 采用“亚跨超声速计算流体力学软件平台”（TRIP2.0）数值模拟了阻力预测小组（AIAA CFD Drag Prediction Workshop Ⅱ,DPWⅡ）翼/身/架/舱复杂组合体运输机构型,数值模拟采用的多块对接网格、测压和测力的试验结果均来自DPWⅡ,对比计算采用了CFL3D的结果。重点针对DLR-F6翼/身/架/舱复杂组合体构型,详细研究了网格密度和湍流模型对总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致。采用SA一方程和SST两方程模型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布有一定的影响。     

耦合运动非定常气动模型对飞机飞行特性仿真的影响

 在大振幅偏航滚转单自由度强迫振荡和耦合强迫振荡运动风洞试验的基础上,分别运用准定常气动力建模方法和非定常气动力建模方法,获取飞机的气动力模型;运用两种不同的气动力模型,对飞机的纵向正弦振荡机动飞行进行了仿真。结果表明,准定常模型在小迎角下的仿真结果与非定常模型的结果基本一致;在较大迎角时,准定常模型仿真结果趋于稳定,而非定常模型仿真结果使飞机进入横航向振荡飞行状态,这一现象与F-16XL的试飞结果类似,其主要原因是耦合运动的非定常模型与准定常模型所产生的阻尼特性不一致。这也表明,分析大迎角飞行特性时应考虑耦合运动下的非定常气动特性。     

模型加力室大涡模拟的并行计算

采用基于MPI(消息传递库)的并行算法,在贴体网格下对带V形槽稳定器模型加力燃烧室紊流化学反应流场进行数值模拟,湍流模型采用k方程亚网格尺度模型,燃烧模型采用亚网格EBU模型,采用热通量辐射模型估算辐射通量。在程序设计中,采用动态内存分配、分区算法和多点重合交错网格系统,并行计算的结果与单机运行结果的对比表明计算结果是正确的,可以明显的提高运算效率,是解决复杂燃烧流动大规模数值模拟的有效手段。      

低泄漏高寿命刷式封严的数值计算

通过空气动力学计算方法探索新型低泄漏高寿命三层片式刷式封严结构。利用商用计算流体力学软件STAR-CD,用定常可压缩N-S方程对刷封进行三维流动模拟,计算不同压差下的泄漏情况。以泄漏量作为评价指标,确定新型刷式封严结构的密封特性并对其进行改进。计算结果表明:随着进出口压差的增大,泄漏量增大;改进的封严刷环结构泄漏性能优于简单的片式刷封结构;"70度模型"的泄漏性能略好于"75度模型",说明增加层与层间的夹角可以达到减少泄漏量的效果。      

波瓣喷管结构参数对引射混合器性能影响的数值研究

通过三维CFD数值计算,研究了波瓣喷管几何结构参数对波瓣喷管引射-混合器的影响规律。计算结果表明:瓣宽增加导致引射流量和波瓣出口处速度环量减小,热混合效率降低,但可以减小混合流动损失,提高总压恢复系数。波瓣扩张角的改变,对引射流量的改变因不同混合管尺寸而异,在混合管直径较小时,随着波瓣扩张角在20~°90°范围内增大,引射流量呈先增加后减小的趋势;在混合管直径较大时,引射量持续提高。扩张角的增加可提高速度环量,但是流动混合损失增加,总压恢复系数减小。      

Improvement of Baldwin-Lomax turbulence model for supersonic complex flows

Entropy represents the dissipation rate of energy. Through direct numerical simulation (DNS) of supersonic compression ramp flow, we find the value of entropy is monotonously decreasing along the wall-normal direction no matter in the attached or the separated region. Based on this feature, a new version of Baldwin-Lomax turbulence model (BL-entropy) is proposed in this paper. The supersonic compression ramp and cavity-ramp flows in which the original Baldwin-Lomax model fails to get convergent solutions are chosen to evaluate the performance of this model. Results from one-equation Spalart-Allmaras model (SA) and two-equation Wilcox k-x model are also included to compare with available experimental and DNS data. It is shown that BLentropy could conquer the essential deficiency of the original version by providing a more physically meaningful length scale in the complex flows. Moreover, this method is simple, computationally efficient and general, making it applicable to other models related with the supersonic boundary layer.     

收敛无穷计算及其应用

经典计算不能很好地刻画无穷计算的行为.基于形式系统序列及其极限,讨论一类称为收敛无穷计算的问题,旨在建立刻画无穷计算在变化的环境中如何交互与演化以及演化的极限状态的逻辑理论基础.提出了收敛无穷计算的一个逻辑和推理系统,其表达能力超过一阶逻辑.还基于经典计算模型图灵(Turing)机和形式系统序列及其极限,提出了收敛无穷计算的模型,称为过程模式.在极限计算的意义下,其计算能力超过了Turing机和实数机器.讨论了上述研究在数据挖掘中的应用,用收敛无穷计算研究了数据挖掘的极限行为.     

三维粘弹性大变形动力响应分析

将大变形增量问题的Ｔｏｔａｌ Ｌａｇｒａｎｇｉａｎ法应用于粘弹性结构的动力响应分析，利用Ｐｒｏｎｇ级数导的Ｋｉｒｃｈｈｏｆｆ应力和Ｇｒｅｅｎ应变表示的粘弹性积分型增量本构方程及其递推公式，从而给出一种计算三维粘弹性大变形动力响应的有限元增量迭加方法，经算例验证，所述方法及其相应程序具有机器内存少、计算量小和效率高的优点。     

超长航时太阳能无人机关键技术综述

超长航时太阳能无人机（UAV）以其高效节能、原理上可实现无限巡航的特点受到广泛关注，而其独特的设计指标与任务特性也对各项关键技术提出了较高要求。多设计要素的高度耦合意味着不同于常规飞行器的总体设计方法，低密度、低速度的飞行条件使其具有明显的低雷诺数气动特性，柔性超大展弦比机翼带来了复杂的气动弹性问题，低翼载荷特性与较大的风场扰动增加了控制难度，极端的飞行环境与苛刻的任务指标对能源、动力系统带来了新挑战，飞行性能对能源系统的高度依赖开辟了飞行轨迹优化的研究方向。本文梳理了超长航时太阳能无人机关键技术的研究现状，在此基础上对各项技术中的难点问题进行了阐释，并对超长航时太阳能无人机未来发展趋势进行了展望。     

Airfoil Aeroelastic Flutter Analysis Based on Modified Leishman-Beddoes Model at Low Mach Number

Based on modified Leishman-Beddoes (L-B) state space model at low Mach number (lower than 0.3), the airfoil aeroelastic system is presented in this paper. The main modifications for L-B model include a new dynamic stall criterion and revisions of normal force and pitching moment coefficient. The bifurcation diagrams, the limit cycle oscillation (LCO) phase plane plots and the time domain response figures are applied to investigating the stall flutter bifurcation behavior of airfoil aeroelastic systems with symmetry or asymmetry. It is shown that the symmetric periodical oscillation happens after subcritical bifurcation caused by dynamic stall, and the asymmetric periodical oscillation, which is caused by the interaction of dynamic stall and static divergence, only happens in the airfoil aeroelastic system with asymmetry. Validations of the modified L-B model and the airfoil aeroelastic system are presented with the experimental airload data of NACA0012 and OA207 and experimental stall flutter data of NACA0012 respectively. Results demonstrate that the airfoil aeroelastic system presented in this paper is effective and accurate, which can be applied to the investigation of airfoil stall flutter at low Mach number.