轴对称进气道外罩跨音速非守恒型全位势流AF3算法

本文针对轴对称进气道外罩跨音速绕流问题的特点,采用了贴体坐标网格;并对这类网格拓扑构造了三种Baker的隐式近似因式分解AF3有限差分迭代算法。通过理论分析与数值实验,找到了既能计算亚临界情形,又能计算超临界情形的AF3格式;并就此格式讨论了适合于进气道绕流特性的AF3迭代格式中间变量的提法,具体提出了格式实施的相应方法,研制了可供分析任意轴对称进气道外罩跨音速特性的计算程序。由于采用了任意曲线坐标系上的轴对称非守恒型全位势方程适合于用不同方法生成的计算网格,因此,本文的工作还有可能进一步推广到计算更复杂的轴对称进气道跨音速绕流问题。计算表明,本方法收敛快,结果好。     

双轴旋转式激光捷联惯导系统的转位方案研究

惯性导航系统的误差随时间累积,旋转调制技术可以有效地提高惯导系统的长航时精度,旋转调制方案是决定旋转式捷联惯导系统导航精度的一个重要因素.针对双轴旋转惯导系统,相较于16次序转位方案,提出了一种新的32次序双轴旋转调制方案.根据捷联惯导系统的误差方程,推导出旋转捷联惯导的误差方程,分析了误差补偿的机理,研究了惯性器件常值偏置误差、标度因数误差和安装角误差的传播特性.仿真结果表明,32次序双轴旋转调制方案相对于16次序转位方案有明显的优势,可以有效地降低姿态角误差和经纬度误差.     

欧洲计算流体力学——个人观点

介绍了70年代初到20世纪末欧洲科学家对CFD发展的主要贡献,并重点关注了对国际流体力学数值方法会议(ICNMFD)和国际计算力学学术会议(ISCFD)的成功举办具有重大贡献的科学家.ICNMFD始于1969年(新西伯利亚),由参与了美国和苏联太空竞赛的两国科学家联合创办,而ISCFD则作为ICNMFD在亚太地区的同类会议初创于1995年(东京).由Oshima教授和本人提议,得益于上述两个会议全体科学委员会成员的共同努力,ICNMFD和ISCFD自2000年开始,合并成为每两年举办一次的国际计算流体力学会议(ICCFD).绝大部分本文讨论所涉及的学者通常都在受邀进行大会报告后,成为了前述三大国际会议中的一个或多个科学委员会的成员.作者对文中讨论所涉及的研究人员的圈定负有全权责任,上述人员的选择更多地基于应用而非理论计算流体力学,因此尚存在不完善之处.此外由于文中所讨论的学者大都是作者在40多年计算流体力学研究生涯中的私交好友,因此并没有将年轻一代的科学家考虑进来,但这并非由于他们的贡献不重要.另外值得注意的是欧洲的计算流体力学同美国的计算流体力学有非常紧密的联系,本文讨论所涉及的研究学者有些最早从欧洲开始涉足计算流体力学而现在则生活和工作于美国(例如Glowinski,Roe,Van Leer等人),而有些则是在美国取得了博士学位(例如Launder,Napolitano),他们都同美国的科学家进行过合作,因此很难清晰地对欧洲和美国的计算流体力学进行界定.对那些现已移居美国的欧洲科学家们,本文优先考虑了他们移居之前的贡献.最后但并非是最不重要的,作者通过电子邮件一一联系了文中所涉及的科学家,希望他们本人能够给出其自身贡献的评价草稿,从而帮助作者以最权威的视角在文中相应章节对其工作进行评述.本文向提供自述评价的诸位致以谢意,同时也希望其他同仁能够对本人不恰当的擅自评述给予谅解.     

一种新的实时参数估计方法

研究导弹控制系统仪器测试时动态模型的实时参数估计 ,利用递归傅里叶变换作实时数据分析 ,提出了一种基于线性状态空间模型的实时参数估计方法。实际应用表明该方法可以得到准确的模型参数估计和适当的误差边界 ,且计算量少、动态模型估计参数收敛性非常好。     

两类激波捕捉格式的性能分析

考虑了两类典型的激波捕捉格式:特征形式的MUSCL(Monotone Upstream-centred Schemes for Conservation Laws)格式和WENO(Weighted Essentially Non-Oscillatory)格式.MUSCL格式在作特征变换时使用了局部线性化的思想,并且针对波系的性质施加相应的限制器;通过逐维重构实现有限体积法的WENO格式.针对一维、二维和三维Euler系统进行数值实验.在一维和二维的情况下,特征形式的MUSCL格式在接触间断的捕捉上具有较明显的优势,而对于激波的捕捉则差别不大.对于三维问题则是WENO格式对流场的分辨更精细.最后对上述结果给出解释,并且提出了可能的改进方法.     

Design for aircraft engine multi-objective controllers with switching characteristics

The aircraft engine multi-loop control system is described and the switching control theory is introduced to solve the regulating and protecting control problems in this paper. The aircraft engine multi-loop control system is firstly described and the control problems are formulated. Secondly, the theory of the smooth switching control is devoted and a new extended scheme for the smooth switching of a switched control system is introduced. Then, for the key technologies of aero-engines switching control, a design algorithm is presented which can determine which candidate controller should be put in feedback with the plant to achieve a desired performance and the procedure to design the aircraft engine multi-loop control system is detailed. The switching performance objectives and the switching scheme are given and a family of PID controllers and compensators is designed. The simulation shows that using the switching control design method can not only improve the dynamic performance of the aircraft engine control system and reduce the switching times, but also guarantee the stability in some peculiar occasions.     

执行器饱和的多智能体一致性控制

针对执行器饱和约束下的一类离散多智能体系统，设计了一种基于线性矩阵不等式的一致性控制律，并对线性矩阵不等式解的存在性进行了分析。由于饱和环节的影响，导致多智能体系统的一致性与个体的初始状态有关，不同的控制律所能允许的个体之间的初始偏差的范围不同。为了定量分析控制律所对应的初始偏差范围，引入吸引域的概念对所设计的控制律进行评估，并对控制律参数进行优化以扩大吸引域的范围。最后，通过数值仿真说明了所设计控制律的有效性。     

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.     

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

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

基于流场/声爆耦合伴随方程的超声速公务机声爆优化

基于自主研发的大规模并行结构化网格CFD求解器PMB3D以及并行化伴随方程求解器PADJ3D，开展了流场/声爆伴随方程的求解研究。首先采用标准算例，对内部CFD代码PMB3D软件和声爆预测代码进行了声爆计算可信度验证，以及声爆强度对近场声压梯度的校核。针对并行环境下多块对接网格的近场声压提取操作的复杂性，提出了"包围盒"的方法实现并行环境下近场声压装配单元编号、网格块编号以及对应的进程编号确定，基于声爆计算坐标将并行传递的数据进行一维排序，为声爆预测、伴随方程以及梯度求解提供输入条件。通过线性插值雅克比矩阵实现均匀坐标系梯度信息向非均匀坐标转换，并进一步根据结构化网格特征提出了插值原则，简化了近场声压转换雅克比矩阵的变分。通过装配单元记录，实现声爆强度对流场守恒变量的变分结果向各个进程装配，将装配结果作为流场伴随方程的右端项实现流场声爆耦合伴随方程的求解。此外，对小型超声速公务机开展了声爆优化，对比分析了设计前后的声压及其频谱特性。