可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展。但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题。界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因。为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案,本文从数值框架的选择、非守恒方程相容离散、高精度有界格式构造、界面压缩、界面-单介质分区计算方法等多个维度出发,综述近几年我们在可压缩多介质流动问题高精度数值模拟方法方面的研究进展。通过上述多个维度的工作,我们建立了一套适用于可压缩多介质流动问题的低耗散、基本无数值振荡的高精度欧拉数值方法,并成功应用于可压缩多介质大变形流动和界面不稳定性诱导湍流混合等问题的数值模拟。相关数值方法研究成果已集成至武器物理等领域工程数值模拟软件中,为相关工程任务提供了重要技术支撑。     

速度脉动下分层旋流火焰动态传播结构提取

为了深刻理解分层旋流火焰的动态结构，开展了分层旋流火焰动态的模态分析和参数化研究。实验在BASIS （北京航空航天大学发展的具有台阶隔离段的独立分层旋流燃烧器）上开展，以甲烷为燃料预混燃烧，使用高速摄像捕捉了速度脉动下的CH*化学发光信号动态图像。利用传统的本征正交分解（POD）和重定向POD分析了火焰动态特性，并利用重定向POD提取了火焰动态中沿着轴向和径向的周期性传播结构，传播结构的频率与速度脉动频率保持一致，且轴向模态的传播结构最为明显。同时，研究了总当量比对重定向POD模态的影响，分析表明：总当量比增大时，重定向POD模态的拓扑结构会发生变化；重定向POD模态中出现了错位现象，这种现象与相平均图像中的火焰传播行为保持一致。重定向POD方法能够有效提取旋流火焰动态中的传播结构，为旋流火焰燃烧不稳定性的机理研究提供了一种分析工具。     

FlowStar:国家数值风洞(NNW)工程非结构通用CFD软件

计算流体力学（CFD）仿真软件是流体相关的数学物理知识和工程实践经验的数字化表达，是工业数字化转型的重要助推。然而，大型工业CFD软件研发难度极高，需要同时兼顾功能多样、系统稳定、性能优越、交互友好等特征。依托国家数值风洞（NNW）工程，研发出一款通用流场模拟软件NNW-FlowStar，并在航空、航天等工业部门大力推广使用。软件基于非结构有限体积求解方法和大规模并行计算技术开发，结合现代化软件工程思维设计，具备先进的数值方法、高效的计算效率和友好的用户操作界面，可满足各类复杂外形的高效气动模拟。独特的重叠网格技术配合六自由度运动模块，可帮助实现武器分离、舱门定轴转动等各类气动-运动协同仿真需求。多类标模案例和复杂工程应用表明，FlowStar软件算法鲁棒、精度可靠，是一款高精度、高效率、高可靠性的通用CFD仿真软件。通过对软件的架构设计和功能应用进行介绍，使相关从业人员能更好地了解FlowStar软件，最终促进国产自主CFD软件生态的良性发展。     

Weibull分布参数估计值对细节疲劳额定强度的影响

将新材料、新工艺应用于民机结构设计时，细节疲劳额定强度（DFR）法仍是重要的评估方法。根据新型铝合金疲劳试验数据，全面分析了Weibull分布形状参数对DFR计算参数的影响；分析标准S-N曲线斜度参数对DFR的影响时，应考虑可靠性寿命的变化，采用低可靠度寿命计算曲线斜度参数可降低Weibull分布形状参数的影响；比较了采用不同分布参数估计值与传统上采用波音公司给定值进行DFR计算的差异。结果表明，对于通过符合性检验的试验数据，使用不同Weibull分布的参数估计值计算得到的DFR都大于按传统方法（用波音公司给定值）的计算结果，其中，三参数Weibull分布计算的DFR最大，至少增加10%以上。     

Multi-scale design and optimization for solid-lattice hybrid structures and their application to aerospace vehicle components

By integrating topology optimization and lattice-based optimization, a novel multi-scale design method is proposed to create solid-lattice hybrid structures and thus to improve the mechanical performance as well as reduce the structural weight. To achieve this purpose, a two-step procedure is developed to design and optimize the innovative structures. Initially, the classical topology optimization is utilized to find the optimal material layout and primary load carrying paths. Afterwards, the solid-lattice hybrid structures are reconstructed using the finite element mesh based modeling method. And lattice-based optimization is performed to obtain the optimal cross-section area of the lattice structures. Finally, two typical aerospace structures are optimized to demonstrate the effectiveness of the proposed optimization framework. The numerical results are quite encouraging since the solid-lattice hybrid structures obtained by the presented approach show remarkably improved performance when compared with traditional designs.     

Effects of wing flexibility on aerodynamic performance of an aircraft model

The low-speed wind tunnel experiment is carried out on a simplified aircraft model to explore the influence of wing flexibility on the aircraft aerodynamic performance. The investigation involves the measurements of force, membrane deformation and velocity field at Reynolds number of 5.4 × 10⁴–1.1 × 10⁵. In the lift curves, two peaks are observed. The first peak, corresponding to the stall, is sensitive to the wing flexibility much more than the second peak, which nearly keeps constant. For the optimal case, in comparison with the rigid wing model, the delayed stall of nearly 5° is achieved, and the relative lift increment is about 90%. It is revealed that the lift enhanced region corresponds to the larger deformation and stronger vibration, which leads to stronger flow mixing near the flexible wing surface. Thereby, the leading-edge separation is suppressed, and the aerodynamic performance is improved significantly. Furthermore, the effects of sweep angle and Reynolds number on the aerodynamic characteristics of flexible wing are also presented.     

A review of topology optimization for additive manufacturing: Status and challenges

Topology optimization was developed as an advanced structural design methodology to generate innovative lightweight and high-performance configurations that are difficult to obtain with conventional ideas. Additive manufacturing is an advanced manufacturing technique building as-designed structures via layer-by-layer joining material, providing an alternative pattern for complex components. The integration of topology optimization and additive manufacturing can make the most of their advantages and potentials, and has wide application prospects in modern manufacturing. This article reviews the main content and applications of the research on the integration of topology optimization and additive manufacturing in recent years, including multi-scale or hierarchical structural optimization design and topology optimization considering additive manufacturing constraints. Meanwhile, some challenges of structural design approaches for additive manufacturing are discussed, such as the performance characterization and scale effects of additively manufactured lattice structures, the anisotropy and fatigue performance of additively manufactured material, and additively manufactured functionally graded material issues, etc. It is shown that in the research of topology optimization for additive manufacturing, the integration of material, structure, process and performance is important to pursue high-performance, multi-functional and lightweight production. This article provides a reference for further related research and aerospace applications.     

煤油燃料超燃冲压发动机燃烧特性实验研究

为研究空气节流时序对超燃冲压发动机点火和火焰稳定的影响，本文通过实验方法研究了13个状态的煤油燃料超燃冲压发动机的燃烧特性，煤油燃烧通过先锋氢气和节流空气增强稳定性。通过两个固定位置的压力传感器来监测火焰稳定状态，采用纹影和OH-PLIF相结合的测量手段，获得了流场结构和火焰发展信息。发动机入口来流条件为Ma = 2.0，总温950 K，总压0.82 MPa。在空气节流的作用下，煤油被先锋火焰引燃；在先锋氢撤除后，煤油仍然可以稳定燃烧。在扩张段中，空气节流和燃烧共同作用产生的激波串移动速度约为52 m/s，但在凹槽内其速度仅为3.7 m/s。通过监测点压力变化情况可以区分所研究状态的火焰稳定与否，通过对13个研究状态的考察，获得了火焰稳定临界曲线。当所研究状态点在临界曲线右上方区域时，火焰状态稳定；当所研究状态点在临界曲线左下方区域时，火焰将被吹熄；当所研究状态点在临界曲线上时，火焰不稳定，在空气节流撤除之前将被吹熄。     

Discrete time multiple spacecraft formation flying attitude optimal control in presence of relative state constraints

This paper addresses the challenge of synchronized multiple spacecraft attitude reorientation in presence of pointing and boundary constraints with limited inter-spacecraft communication link. Relative attitude pointing constraint among the fleet of spacecraft has also been modeled and considered during the attitude maneuvers toward the desired states. Formation fling control structure that consists of decentralized path planners based on virtual structure approach joint with discrete time optimal local controller is designed to achieve the mission’s goals. Due to digital computing of spacecraft’s onboard computer, local optimal controller based on discrete time prediction and correction algorithm has been utilized. The time step of local optimal algorithm execution is designed so that the spacecraft track their desired attitudes with appropriate error bound. The convergence of the proposed architecture and stability of local controller’s tracking error within appropriate upper bound are proved. Finally, a numerical simulation of a stereo imaging scenario is presented to verify the performance of the proposed architecture and the effectiveness of the algorithm.     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法，能够在一维计算域上准确捕捉到湍流统计规律，且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究，若要将其拓展用于模拟高速可压缩湍流，需对建模方法进行深度改进。相比于不可压ODT方法，本文基于欧拉参考框架，针对可压缩湍流的特性，将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量，并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进，开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程，在随机性模块中构造一维涡时，将三联映射的作用对象也相应地由原始变量改为守恒通量，并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对，验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合，精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM），且该方法的降维处理使其在降低计算成本方面具有显著优势。