基于有序BDL树集的船舶电网潮流算法及执行时间分析

针对嵌入式船舶电力模拟训练系统在频繁模拟操作指令下对潮流计算的快速性需求以及嵌入式模拟训练系统极其有限的运算资源之间的矛盾,提出一种基于有序BDL树集的船舶电网潮流算法。该算法在常规回推前推法基础上对节点编号进行优化,将多电站船舶电力系统在实际开环运行模式下等效为由多个辐射状电网构成的有序BDL树集,树集中的每个有序BDL树包含单个电站中的主配电板层、配电层和负载层,并把发电机潮流整合入主配电板层,避免了常规潮流算法中节点编号突变问题。对该算法的时间复杂度分析结果及实际运行结果进行比较分析,证明该算法比常规船舶电网潮流算法更优。     

A self-similar solution of a curved shock wave and its time-dependent force variation for a starting flat plate airfoil in supersonic flow

The problem of aeroelasticity and maneuvering of command surface and gust wing interaction involves a starting flow period which can be seen as the flow of an airfoil attaining suddenly an angle of attack. In the linear or nonlinear case, compressive Mach or shock waves are generated on the windward side and expansive Mach or rarefaction waves are generated on the leeward side. On each side, these waves are composed of an oblique steady state wave, a vertically-moving one-dimensional unsteady wave, and a secondary wave resulting from the interaction between the steady and unsteady ones. An analytical solution in the secondary wave has been obtained by Heaslet and Lomax in the linear case, and this linear solution has been borrowed to give an approximate solution by Bai and Wu for the nonlinear case. The structure of the secondary shock wave and the appearance of various force stages are two issues not yet considered in previous studies and has been studied in the present paper. A self-similar solution is obtained for the secondary shock wave, and the reason to have an initial force plateau as observed numerically is identified. Moreover, six theoretical characteristic time scales for pressure load variation are determined which explain the slope changes of the time-dependent force curve.     

涡扇发动机低压部件通流耦合计算

为了研究各部件匹配与过渡段几何参数对发动机总体性能的影响,针对某涡扇发动机风扇、低压涡轮、内外涵道、掺混段和尾喷管等部件,基于周向平均Navier-St okes通流控制方程,采用时间推进有限体积法进行了部件耦合计算求解,将沿程各部件性能计算参数与试验做了对比,结果表明:在设计状态下发动机沿程各截面计算特征参数与试验值最大误差不超过5.0%,说明该通流耦合计算方法能够快速预估发动机总体性能,完善各部件匹配关系,在发动机总体设计初期有一定工程应用价值。     

跨声速层流翼型的混合反设计/优化设计方法

跨声速层流翼型设计须兼顾优良的超临界特性和自然层流特性,因而对设计方法提出了更高的要求。针对现有反设计方法和直接优化设计方法的不足,发展了一种适用于跨声速层流翼型的混合反设计/优化设计方法。该方法引入了基于经验的局部流场特征作为反设计目标,翼型性能指标作为直接优化设计目标,然后加权形成了混合反设计/优化设计总目标,并同时考虑了气动和几何约束。优化算法采用基于自适应并行加点技术的代理优化,流动数值模拟采用耦合基于线性稳定性理论的e^N转捩自动判定的雷诺平均Navier-Stokes（RANS）方程求解器。针对现代中短程民用客机需求,以NPU-LSC-72613翼型为基准,开展了层流翼型减阻的混合反设计/优化设计。分别将局部目标压力分布、总阻力作为反设计和直接优化设计目标,得到了较好的优化结果,验证了方法的有效性。经过2轮优化结果显示混合反设计/优化设计总目标显著下降。所设计翼型吸力面局部压力分布与目标压力分布基本一致,总阻力下降15.5%;吸力面和压力面层流范围均大于55%倍弦长,激波强度显著减弱,说明所设计翼型同时具有优良的超临界和层流特性。将所设计翼型配置到机翼上,通过三维数值模拟进行校验,结果显示所设计跨声速层流机翼升阻比提高了6.64%,在一定升力系数范围内,气动性能均有显著提高,验证了所设计跨声速层流翼型在机翼设计中的适用性。     

高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。     

地面燃气轮机单管燃烧室流量分配试验

对设计的100kW地面燃气轮机单管燃烧室空气流量分配进行了试验研究。在常压环境下采用堵孔法分别得到了旋流器、主燃孔和掺混孔的流量特性曲线,分析得到燃烧室不同结构的流量系数;试验研究了不同进口流量条件下燃烧室的流量分配,测量得到了燃烧室总压损失。研究发现:随着进口空气流量的增加,燃烧室的流量分配比例基本保持稳定,并且与燃烧室的设计空气比例基本吻合;随着燃烧室进口流量（雷诺数）的增加,旋流器、主燃孔和掺混孔的流量系数呈线性降低;随着进口流量（雷诺数）的增加,燃烧室总压损失逐渐增大;对主燃孔和掺混孔的流量特性测量中,两种测量方法得到的试验结果稍有差别,总体上看两种测量方法的试验结果较为接近。通过对比分析证明两种试验测量方法真实可靠,该研究结果可为100kW地面燃气轮机燃烧室的设计与优化提供依据。      

宽高比及表面粗糙度对矩形微尺度通道流动特性的影响

实验研究了不同宽高比矩形金属微尺度通道的流动特性,并与硅通道对比来探究表面粗糙度的影响。金属微尺度通道宽为0.4mm,通道宽高比分别为0.50、0.67、和1.00,相对表面粗糙度范围为0.1609%~0.2145%。硅通道宽度及高度都为0.4mm,相对表面粗糙度为0.00325%。实验工质为空气,实验的雷诺数范围是250~3000。实验结果表明:在粗糙微尺度通道会发生转捩提前,并且转捩雷诺数随着宽高比增大而减小。在层流区,微尺度矩形通道摩擦因数随着雷诺数增大而减小,在转捩区出现明显增大后再减小。表面粗糙度很小时,未发现转捩提前现象。