多体全自由飞过程的数值模拟

以三维非结构粘性直角网格重构为构架,数值求解了固定空间网格瞬间运动边界的N-S方程,逼近了流场的非定常过程,结合六自由度的弹道方程,对单体或多体运动过程进行了全自由飞数值模拟,数值模拟结果和级间分离的风洞自由飞试验进行了对比,初步的计算结果说明,全自由飞数值模拟在一定程度上可以描述风洞试验的飞行轨迹,并且适应能力更强,整个全自由飞数值模拟系统可以为更复杂的单体机动飞行或多体运动过程提供快速实用的数值模拟方法.     

Adaptive mesh refinement computation of acoustic radiation from an engine intake

A block-structured adaptive mesh refinement (AMR) method was applied to the computational problem of acoustic radiation from an aeroengine intake. The aim is to improve the computational and storage efficiency in aeroengine noise prediction through reduction of computational cells. A parallel implementation of the adaptive mesh refinement algorithm was achieved using message passing interface. It combined a range of 2nd- and 4th-order spatial stencils, a 4th-order low-dissipation and low-dispersion Runge–Kutta scheme for time integration and several different interpolation methods. Both the parallel AMR algorithms and numerical issues were introduced briefly in this work. To solve the problem of acoustic radiation from an aeroengine intake, the code was extended to support body-fitted grid structures. The problem of acoustic radiation was solved with linearised Euler equations. The AMR results were compared with the previous results computed on a uniformly fine mesh to demonstrate the accuracy and the efficiency of the current AMR strategy. As the computational load of the whole adaptively refined mesh has to be balanced between nodes on-line, the parallel performance of the existing code deteriorates along with the increase of processors due to the expensive inter-nodes memory communication costs. The potential solution was suggested in the end.     

Shape control of spacecraft formation using a virtual spring-damper mesh

This paper derives a distance-based formation control method to maintain the desired formation shape for spacecraft in a gravitational potential field. The method is an analogy of a vir-tual spring-damper mesh. Spacecraft are connected virtually by spring-damper pairs. Convergence analysis is performed using the energy method. Approximate expressions for the distance errors and control accelerations at steady state are derived by using algebraic graph representations and results of graph rigidity. Analytical results indicate that if the underlying graph of the mesh is rigid, the convergence to a static shape is assured, and higher formation control precision can be achieved by increasing the elastic coefficient without increasing the control accelerations. A numerical exam-ple of spacecraft formation in low Earth orbit confirms the theoretical analysis and shows that the desired formation shape can be well achieved using the presented method, whereas the orientation of the formation can be kept pointing to the center of the Earth by the gravity gradient. The method is decentralized, and uses only relative measurement information. Constructing a distributed virtual structure in space can be the general application area. The proposed method can serve as an active shape control law for the spacecraft formations using propellantless internal forces.     

航空发动机冰晶结冰研究进展

飞机在高空区域巡航时可能遭遇冰晶气象，导致发动机吸入冰晶颗粒发生结冰故障，欧美等国已将冰晶结冰纳入适航 取证要求，中国也考虑在后续发动机适航取证中增加冰晶结冰考核项目。为了准确掌握国内外冰晶结冰研究现状，通过调研文 献，对比了冰晶结冰与过冷水滴结冰在结冰位置、结冰条件和结冰机理方面的区别，阐述了试验研究和数值仿真方面的重要进展， 其中，试验研究分别介绍了冰晶在静态部件和动态部件内的撞击与结冰过程，数值仿真按照空气流场计算、冰晶运动轨迹和撞击 特性计算及结冰计算的顺序，总结了欧拉法与拉格朗日法的区别、冰晶运动相变与黏附过程及结冰热物理模型，并介绍了冰晶结 冰软件在数值模拟中的广泛应用，提出了目前在冰晶碰撞、结冰热力学模型、结冰探测及风洞试验等方面需要进一步研究的内容， 以期为中国冰晶结冰研究提供参考。     

离轴三反(TMA)相机在轨成像的偏流角计算与控制

推导了离轴三反(TMA)相机在轨成像的偏流角计算公式,并根据推导过程给出了偏流角的闭环和开环两种控制方式。分析了转序定义对控制的影响,获得了最佳转序。     

基于离散等价方程的非结构网格有限差分法

激波装配法把解析关系式嵌入流场避免了间断引起的理论问题，使全场一致高精度的实现成为可能，有望解决超声速流动转捩研究中的感受性模拟难题。但装配复杂激波结构以后，被分割的流场空间常出现不规则的几何形状，这给常规的基于结构网格的有限差分法应用带来困难。基于离散等价方程理论，提出了一种新的基于非结构网格的有限差分法，在空间二维离散点附近仅用3条网格线就可以构造出一阶迎风格式。数值算例表明收敛过程对网格质量不敏感，解决了激波装配法模拟激波相交出现小夹角后使用结构网格进行计算存在的难题。根据这种方法的特点展望了未来的应用前景。     

非结构混合网格自适应并行技术

计算流体力学（CFD）模拟实际工程问题所采用的网格规模可达千万量级，并行技术是减少计算时间的有效方法。耦合流场信息的网格自适应技术能有效动态优化计算网格，被NASA视为一项亟待发展的CFD关键技术。混合网格自适应系统包含网格分布优化、表面网格投影和空间网格匹配等关键技术。针对以上3项关键技术分别建立了高效的并行算法。首先，提出了"先唯一后同一"的两步法策略实现了网格单元分布优化过程的并行相容性；其次，基于局部曲面拟合思想，实现了曲面重构和新增物理网格点投影的完全并行；再次，提出了空间网格匹配技术的半并行算法，快速解决了网格单元交错问题。为了提高后续流场计算的并行效率，发展了基于并行重分区-网格数据迁移方法的动态负载平衡技术，并采用圆柱激波流场自适应模拟对动态负载平衡技术进行初步验证。最后，采用三角翼自适应加密测试了自适应系统的并行效率。结果表明，建立的混合网格自适应系统并行效率较高，且相比流场求解耗费总时间的比例低于1%。     

克服旋转飞行器螺旋运动的方法研究

旋转飞行器虽然具有内在稳定性，但在飞行过程中一些不确定因素会导致飞行器产生螺旋运动，而非理想的关于角动量矢量纯粹的旋转运动。结合旋转飞行器自身的特点，借鉴自旋卫生的被动章动阻尼的方法，提出了一种通过主动控制飞行器内部质量块的运动使得飞行器碑的惯量主轴发生偏移，惯量矩阵由原来的对称阵变为含有非零非对角元素的矩阵，从而改变飞行器的姿态角运动来克服螺旋运动的方法，这将有助于提高旋转飞行器的资态控制精度。以某型旋转飞行器为背景进行仿真研究，计算结果证实了该方法的有效性和可行性。     

机场道面动态响应分析

 本文首先研究了粘弹性地基上四边自由矩形板在运动载荷作用下的动态响应问题,得到了问题的近似解答。然后从该解答出发讨论了机场道面的动态响应问题,并以××机场道面板为计算模型,将理论计算结果和实验结果进行比较,二者符合较好。最后计算了道面的固有频率、临界速度,讨论了主要物理参数对道面动挠度的影响规律。     

非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。