Global aerodynamic design optimization based on data dimensionality reduction

In aerodynamic optimization, global optimization methods such as genetic algorithms are preferred in many cases because of their advantage on reaching global optimum. However, for complex problems in which large number of design variables are needed, the computational cost becomes prohibitive, and thus original global optimization strategies are required. To address this need, data dimensionality reduction method is combined with global optimization methods, thus forming a new global optimization system, aiming to improve the efficiency of conventional global optimization. The new optimization system involves applying Proper Orthogonal Decomposition (POD) in dimensionality reduction of design space while maintaining the generality of original design space. Besides, an acceleration approach for samples calculation in surrogate modeling is applied to reduce the computational time while providing sufficient accuracy. The optimizations of a transonic airfoil RAE2822 and the transonic wing ONERA M6 are performed to demonstrate the effectiveness of the proposed new optimization system. In both cases, we manage to reduce the number of design variables from 20 to 10 and from 42 to 20 respectively. The new design optimization system converges faster and it takes 1/3 of the total time of traditional optimization to converge to a better design, thus significantly reducing the overall optimization time and improving the efficiency of conventional global design optimization method.     

Viscoelastic behaviour of pumpkin balloons

The lobes of the NASA ULDB pumpkin-shaped super-pressure balloons are made of a thin polymeric film that shows considerable time-dependent behaviour. A nonlinear viscoelastic model based on experimental measurements has been recently established for this film. This paper presents a simulation of the viscoelastic behaviour of ULDB balloons with the finite element software ABAQUS. First, the standard viscoelastic modelling capabilities available in ABAQUS are examined, but are found of limited accuracy even for the case of simple uniaxial creep tests on ULDB films. Then, a nonlinear viscoelastic constitutive model is implemented by means of a user-defined subroutine. This approach is verified by means of biaxial creep experiments on pressurized cylinders and is found to be accurate provided that the film anisotropy is also included in the model. A preliminary set of predictions for a single lobe of a ULDB is presented at the end of the paper. It indicates that time-dependent effects in a balloon structure can lead to significant stress redistribution and large increases in the transverse strains in the lobes.     

水平表面气流剪切作用下的水膜厚度

飞机结冰表面上的液态水受气流吹拂作用会发生向后溢流,从而影响结冰区域范围及防冰系统设计;为了获得水膜流动规律,对水平平板表面上气流剪切驱动的水膜流动进行了实验测量和建模分析。通过水膜流动风洞试验台产生高速气流驱动水膜的流动,使用色散共焦位移计测量同一位置的水膜在不同时刻的厚度变化,结果表明气-液界面由底层薄水膜和多种尺度的波动组成,具有变化速度快随机性强的特点。通过水膜厚度随气流速度及水膜雷诺数的变化规律,发现平均水膜厚度与两者均呈现出单调非线性的依赖关系。基于薄水膜流动理论和平均水膜厚度实验结果,提出了高速气流剪切作用下的气-液波动界面剪切因子计算式,适用于风速17.8~52.2m/s,水膜雷诺数26~128之间的平板水膜流动计算。     

Optimization of entry-vehicle shapes during conceptual design

During the conceptual design of a re-entry vehicle, the vehicle shape and geometry can be varied and its impact on performance can be evaluated. In this study, the shape optimization of two classes of vehicles has been studied: a capsule and a winged vehicle. Their aerodynamic characteristics were analyzed using local-inclination methods, automatically selected per vehicle segment. Entry trajectories down to Mach 3 were calculated assuming trimmed conditions. For the winged vehicle, which has both a body flap and elevons, a guidance algorithm to track a reference heat-rate was used. Multi-objective particle swarm optimization was used to optimize the shape using objectives related to mass, volume and range. The optimizations show a large variation in vehicle performance over the explored parameter space. Areas of very strong non-linearity are observed in the direct neighborhood of the two-dimensional Pareto fronts. This indicates the need for robust exploration of the influence of vehicle shapes on system performance during engineering trade-offs, which are performed during conceptual design. A number of important aspects of the influence of vehicle behavior on the Pareto fronts are observed and discussed. There is a nearly complete convergence to narrow-wing solutions for the winged vehicle. Also, it is found that imposing pitch-stability for the winged vehicle at all angles of attack results in vehicle shapes which require upward control surface deflections during the majority of the entry.     

基于Kriging模型插值的孔位修正策略

飞机装配时的孔位精度直接影响其最终装配质量，而待制孔零件易发生变形和整体偏移，因此一般通过基准孔信息对待制孔位进行补偿。提出一种基于Kriging插值的孔位修正方法，通过基准孔建立孔位理论坐标与实际偏差值的Kriging模型，进而预测待制孔的孔位偏差，并给出每个位置的预测标准误差，以指导基准孔的增添和布置。最后通过有限元和数值实验分别验证了零件在极限变形和整体平移旋转情况下该方法的有效性，实现了孔位估计误差小于0.3 mm。     

金属多层膜超晶格结构研究

用高分辨的三晶体x—射线衍射仪测量了Zr、Cu金属交替的多层膜的x—射线小角反射谱。用扩展了的Born—wolf光学模型对Zr、Cu多层膜的x—射线小角反射数据进行了定量分析。结果表明，该金属多层膜在其生长方向具有超晶格结构。在两层金属膜之间的界面上，由于内扩散，形成了界面交混层，它使x—射线小角反射所产生的高阶Bragg峰的强度减弱。由于两种金属膜的晶格不匹配，使高阶Bragg峰加宽。界面的粗糙度可用来模拟x—射线小角反射曲线的阻尼振荡效应。     

一种大面积比喷管的分段式设计与数值分析

为了满足大推力上面级发动机大面积比喷管的设计需求,采用了排放冷却前段和单壁辐射冷却尾段的分段式设计方案。在排放冷却前段传热计算的基础上,通过对内流场进行数值模拟,重点研究了单壁尾段在引入上游排放冷却气氢情况下的冷却特性和喷管效率。结果表明:对于大面积比喷管,采用带二次流的单壁金属喷管延伸段是现实可行的,有望达到较好的冷却保护效果并提高喷管效率。     

气膜孔内局部堵塞对气膜冷却特性的影响

在吹风比为0.3~1.5,堵塞比为0.1~0.5范围内,对平板上单排倾斜气膜孔内局部堵塞所引起的冷气射流流动和冷却特性变化进行了数值研究,分析了堵塞比、堵塞位置和吹风比对绝热气膜冷却效率的影响.在研究参数范围内的研究结果表明:气膜孔出口端前缘的局部堵塞有利于抑制肾形涡,对气膜孔下游的绝热气膜冷却效率有改善作用,并且随着堵塞比和吹风比的增加,对冷却效果的改善越为明显.而在气膜孔出口端尾缘或侧边的堵塞则在较大的堵塞比下削弱气膜冷却效果.相对于气膜孔出气端的局部堵塞,在气膜孔进气端和中部的堵塞对绝热气膜冷却效率的影响要微弱得多.     

超燃冲压发动机燃烧室主动冷却设计研究

超燃冲压发动机燃烧室主动冷却模型制造和热考核试验的费用高昂、周期长,为了降低试验风险,采用三维计算和经济的验证试验相结合的方法开展了超燃冲压发动机燃烧室主动冷却设计。冷却结构基本参数设计和侧壁冷却流动设计是确定设计方案的两大基础,前者采用三维传热计算结合冷却面板传热验证试验完成,后者采用超临界燃料流动三维并行计算结合水流动验证试验完成。在此基础上,经过多轮的结构设计与三维传热及强度计算评估迭代,确定了最终的燃烧室主动冷却结构。设计的主动冷却燃烧室在来流马赫数2.5,总温1700K条件下成功通过200s热考核试验,表明所采用的设计方法、验证试验和计算工具是有效和可信的。     

肋角度对气膜冷却特性的影响

在光滑二次流通道的基础上,分析对比了两种带肋通道(135°肋和45°肋)对气膜冷却特性的影响.采用瞬态液晶测试技术获得了气膜孔下游表面传热系数比与气膜冷却效率分布.使用Fluent软件RANS数值方法对相应结构进行了数值模拟,并使用了realizable k-ε湍流模型.光滑二次流通道模型中,气膜孔内流线呈螺旋状分布,导致较大的孔内速度分离与流动损失.冷气射流分成两部分,其中一股形成一对偏斜的对转涡.135°肋结构中,二次流通道上部分的旋转涡为顺时针方向,使得气流易于流入气膜孔,气膜孔内流线呈直线分布.45°肋结构中,二次流通道上部分旋转涡为逆时针方向,增强了气膜孔内旋转涡.45°肋结构中冷气流入气膜孔之后的流动结构与光滑二次流通道结构相似.135°肋结构气膜冷却效率最大而表面传热系数比最低.