非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

星载大型网状天线非线性结构系统有限元分析

论述了大型网状可展开天线的非线性理论模型及其非线性模拟有限单元模型,并以一副4．2m天线为例进行了有限元建模分析,获得了这类天线的动力学振动特性。研究结果表明：拉索几何非线性和间隙接触非线性对大型网状天线系统的动力学特性具有较大影响。     

基于能量等效的直齿锥齿轮时变啮合刚度计算

针对直齿圆锥齿轮啮合刚度的计算问题,从微元思想出发,将变截面齿廓划分为若干等截面微段齿段,基于能量等效建立微段齿段啮合刚度计算模型,并利用积分方法获得单齿啮合刚度。此外,基于力平衡和变形协调条件进一步提出了齿轮时变啮合刚度计算模型,同时根据几何关系导出了相应传动误差计算式。利用有限元分析对解析计算模型进行了验证,并分析了误差来源。结果表明:利用该模型不仅能够将直齿锥齿轮啮合刚度计算精度保证在2%以内,还达到了快速求解的目的。      

动网格技术在增升装置优化中的应用

将基于Denaulay背景图的动网格方法应用于增升装置二维缝道参数优化设计过程中。针对二维增升装置的构型特点和翼型运动的特点,采用了添加辅助点的双层Denaulay背景图动网格方法,并详细论述了辅助点添加、运动技术。为了进一步提高动态网格质量,提出了针对前后缘运动位置中比较极端的两种位置生成两个初始网格,并采用两种初始网格进行动态网格的生成。结果表明新动态网格方法可以为二维增升装置提供高质量的动态网格。     

W型地下井发射环境数值模拟与分析

为研究导弹地下井发射环境,采用计算流体力学方法对W型井导弹发射过程的流场进行了数值模拟。通过将仿真数据与理论计算及试验结果进行对比,验证了计算模型和数值方法的准确性。采用动网格技术处理流场计算区域随导弹运动的变化过程。仿真结果表明,在整个发射过程中燃气流具有反射、引射和充分发展3个阶段,造成发射环境参数有较大变化。分析了压力和温度参数的变化规律。该数值算法和仿真结果可推广应用于发射井的设计、优化及改进。     

三维热化学非平衡流动非结构网格DSMC方法及其应用

研究了过渡流域三维热化学非平衡流动DSMC方法实现的过程。以四面体非结构网格为基本单元，提出一种新型的高效搜索算法，该算法不仅可以快速跟踪模拟分子在网格之间的迁移，而且可以准确判别分子与物面是否相互作用，避免了原有算法中分子表面反射的非确定论判据。设计了适合三维DSMC方法的动态局部时间步长技术，节约了计算时间。将碰撞距离的思想引人到非结构网格上来，有效地消除了网格尺度小于三分之一平均自由程的限制。利用Fortran90的动态分配内存技术编制了适用于任意外形的通用计算程序。最后对全尺寸航天飞机高超声速绕流进行了数值试验，计算结果验证了该算法的可行性及高效性。     

基于CFD动网格的内装式空射分离研究

研究了使用飞机发射运载火箭时的机箭分离过程.提出了一种可以解决内装式空射分离问题的整体解决方案,主要基于CAD三维实体建模方法、多体动力学方法和结合六自由度运动的计算流体力学(CFD)动网格方法.以CFD动网格计算为核心,进行了空射分离过程关键问题的研究,得到了火箭在载机气动干扰环境下的运动特性,为轨道卫星和弹道导弹的快速、廉价发射打下了研究基础.     

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.     

基于非结构自适应网格和嵌入边界法的残差格式研究

嵌入边界法由于在求解NS方程时能够简化网格生成问题而在计算流体领域受到越来越广泛的关注.简言之,嵌入边界法能够简化大变形和运动条件下多物理流动模拟、流固相互作用耦合问题,然而壁面边界条件的精确处理仍旧是该方法需要解决的问题.在本文工作中,为考虑壁面边界条件而在NS方程中增加了补偿项,同时采用非结构网格自适应技术保持了壁面边界条件的精度.