高精度微分求积曲梁单元的建立与应用

首先由能量原理导出曲梁弯曲问题的控制微分方程，在此基础上应用微分求积法原理分别给出了曲梁内点和端点的微分求积方程，由此形成曲梁单元的刚度方程，从而建立了微分求积曲梁单元，并给出了曲梁结构刚度方程的边界条件。通过算例分析，得到了微分求积单元法结构离散时应使单元数量少的原则和求解精度与单元长度比基本无关的性质。与有限元方法的结果比较表明，本文导出的曲梁单元是一种具有很高求解精度的单元。》     

一种新的涡扇发动机动力特性计算模型

以等熵流动原理为出发点,结合发动机短舱入口、出口及远场特征边界条件,建立了一种新的涡扇发动机动力特性计算模型.在基于非结构网格的并行计算软件中实现了模型,对轴对称超高涵道比涡扇模型(CRUF)、轴 、对称涡轮动力模型(TPS)以及某型民用飞机整机带发动机模型进排气流场进行了数值模拟,计算结果与传统的发动机动力特性模型计算结果以及实验数据进行比较,三者吻合较好,从而验证了构造的涡扇发动机动力特性计算模型的正确性、可靠性以及工程实用性.     

含摩擦和加热的跨声速流动数值模拟研究

喷管中跨声速流动的反应流体动力学是实验和理论研究的一个重要领域.本文主要由两部分构成,在第一部分中对文献[2]中的Murman格式[1]的一般形式进行了简化,得到了守恒形式的欧拉方程,并对其求解方法进行了讨论,同时针对Chattot等的工作展开了分析,并给出了带激波的流动计算结果;而文章的第二部分则研究了包含摩擦、加热以及耗散的跨声速流动,并重点关注了一步Arrhenius化学反应模型下的爆震现象.本文还探究了摩擦力和带燃烧的跨声速流动之间的相互作用,同时详细给出了新的欧拉方程数值求解算法并同已有文献算例进行了对比验证.     

Numerical simulation of aerodynamic interaction for a tilt rotor aircraft in helicopter mode

A rotor CFD solver is developed for simulating the aerodynamic interaction phenomenon among rotor, wing and fuselage of a tilt rotor aircraft in its helicopter mode. The unsteady Navier–Stokes equations are discretized in inertial frame and embedded grid system is adopted for describing the relative motion among blades and nacelle/wing/fuselage. A combination of multi-layer embedded grid and ‘‘extended hole fringe" technique is complemented in original grid system to tackle grid assembly difficulties arising from the narrow space among different aerodynamic components, and to improve the interpolation precision by decreasing the cell volume discrepancy among different grid blocks. An overall donor cell searching and automatic hole cutting technique is used for grid assembly, and the solution processes are speeded up by introduction of Open MP parallel method. Based on this solver, flow fields and aerodynamics of a tilt rotor aircraft in hover are simulated with several rotor collective angles, and the corresponding states of an isolated rotor and rotor/wing/fuselage model are also computed to obtain reference solution.Aerodynamic interference influences among the rotor and wing/fuselage/nacelle are analyzed,and some meaningful conclusions are drawn.     

用非结构直角网格和欧拉方程计算运载火箭绕流

采用八叉树结构，物面附近用投影法，生成运载火箭绕流计算的非结构直角网格。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。算例表明，网格生成的时机化费很少，质量好；欧拉方程解算的收敛质量也好。     

Riccati方程最大解的一个重要定理

Riccati 方程在最优控制和滤波理论中起着非常重要的作用。最近,Kawasaki 和 Shimereura 引用了 Riccati 方程最大解的一个性质,但他们的证明是不正确的。本文对这个问题做了进一步研究,并提出和证明了 Riccati 方程最大解的一个重要定理。     

物体穿越斜激波干扰特性的数值研究

当物体以超声速穿越同向运动的斜激波时，它的飞行Mach数和攻角都会发生较大的变化。穿越过程中斜激波与物体脱体激波发生相互作用，形成复杂的激波干扰波系结构，基本的干扰形式为复合Mach反射结构。不同的穿越位置时，由于不同的来流条件，使物体上尾部凹坑内的流场和波系发生很大的变化，随着物体的穿越过程，当头部和凹坑下游壁面与翼面穿越激波时，阻力系数有两次剧烈的变化，法向力基本呈线性下降，压力中心先大幅度后移，物体逐渐穿过激波后又快速前移。因此穿越过程使物体加速，给物体的姿态带来很大的扰动。翼面穿过激波也造成阻力、法向力下降，压心后移。如果对有控制的飞行器，穿越激波时气动力特性和稳定性的急剧变化，将给控制系统带来很大的困难。     

无碰撞激波的理想MHD模拟

采用了WENO格式数值求解一维理想MHD方程组，模拟了行星际无碰撞激波，研究了垂直无碰撞激波与行星际反向磁场结构和高密度等离子体团的相互作用过程，并与粒子模拟的结果进行比对，两者的结果非常类似．模拟结果表明，对涉及无碰撞激波的大部分现象中理想MHD模拟是准确且可行的，同时相对于粒子模拟又有很好的计算效率，便于扩展至二维或三维的情形．     

Implicit discontinuous Galerkin method on agglomerated high-order grids for 3D simulations

High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is inves-tigated in combination with discontinuous Galerkin (DG) method for solving the 3D compressible Euler and Navier-Stokes equations. In this method, a fine linear mesh is first generated by standard commercial mesh generation tools. By taking advantage of an agglomeration method, a quadratic high-order mesh is quickly obtained, which is coarse but provides a high-quality geometry represen-tation, thus very suitable for high-order computations. High-order discretizations are performed on the obtained grids with DG method and the discretized system is treated fully implicitly to obtain steady state solutions. Numerical experiments on several flow problems indicate that the agglomer-ated high-order mesh works well with DG method in dealing with flow problems of curved geome-tries. It is also found that with a fully implicit discretized system and a p-sequencing method, the DG method can achieve convergence state within several time steps which shows significant effi-ciency improvements compared to its explicit counterparts.     

多喷口环量控制翼型流动的研究

提出了多喷口环量控制翼型的想法，并进行了数值研究。结果表明，多喷口环量控制方法可以在中、小Ｃμ吹气时提高翼型的升力和减小能耗，并能改善环量控制翼型在较高来流马赫数下的气动性能，避免“压缩性失速”。