一种中心型间断有限元MMALE方法

针对多介质可压缩流体动力学问题,提出了一种单元中心型二维MMALE(Multi-Material Arbitrary Lagrangian-Eulerian)方法.在拉氏步,流体力学方程组采用中心型间断有限元方法求解.对于混合网格,采用Tipton压力松弛模型更新物理量,用等参坐标法更新物质中心点坐标.界面重构采用一种...     

某型无人机发射段飞行轨迹研究

用火箭助推起飞是广泛应用的无人机起飞的一种方式,它和用起飞车助跑起飞方式相比有一些优越之处。本文总结了在某型无人机采用火箭助推起飞时的飞行轨道研究中探讨过的几个问题,包括发射参数的选择原则和发射起飞过程中的一些特殊问题.如力矩平衡问题,无人机和滑轨之间相对运动问题等。     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究,将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元,采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元,更好地表述了真实物面特征;物面边界相邻体网格单元相应构造为高阶体网格单元,其余体网格单元采用八节点六面体单元,以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解,验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外,由于采用了隐式计算方法,仅需较少的时间步就能迭代收敛。     

Nutation instability of spinning solid rocket motor spacecraft

The variation of mass,and moment of inertia of a spin-stabilized spacecraft leads to concern about the nutation instability.Here a careful analysis on the nutation instability is performed on a spacecraft propelled by solid rocket booster (SRB).The influences of specific solid propellant designs on transversal angular velocity are discussed.The results show that the typical SRB of End Burn suppresses the non-principal axial angular velocity.On the contrary,the frequently used SRB of Radial Burn could amplify the transversal angular velocity.The nutation instability caused by a design of Radial Burn could be remedied by the addition of End Burn at the same time based on the study of the combination design of both End Burn and Radial Burn.The analysis of the results proposes the design conception of how to control the nutation motion.The method is suitable to resolve the nutation instability of solid rocket motor with complex propellant patterns.     

某风扇/增压级非设计点性能计算及进气畸变影响分析

在原流线曲率法程序的基础上, 增加了适用于高速风扇/压气机的双激波损失模型.真实地反映了高速风扇/压气机的实际工作状况, 提高了预测高速风扇/压气机非设计点性能的准确性, 完善了分析进气畸变对风扇/压气机性能及稳定性影响的准三维计算模型.针对一台外涵转子叶尖马赫数大于1.4的大涵道比风扇/增压级, 计算分析了其非设计点性能以及进气畸变的影响.计算结果与试验结果吻合较好, 精度可靠, 可应用于工程实际.      

大涵道比涡扇发动机风扇/ 增压级试验件结构 设计及验证

某大涵道比风扇/增压级试验件为1级风扇加3级增压级结构,为实现风扇/增压级部件特性录取,结构设计在满足气动设计要求的基础上,重点考虑了合理性和可靠性问题,采用理论结合工程经验的方法完成了转子支撑、静子连接和轴承润滑及密封等难点设计,兼顾大涵道比风扇/增压级结构及气动特点,采用宽弦风扇叶片设计、榫头/榫槽圆弧榫连及装配式风扇转子流路板结构。顺利完成试验件研制工作,经试验验证,压比、效率等性能参数达到或超过设计指标。     

大涵道比风扇/增压级进气畸变数值模拟研究

为探究进口总压畸变下大涵道比风扇/增压级内部流场的主要特征,基于三维彻体力模型的思想,开发了一套能够实现风扇/增压级内外涵联算功能的三维数值计算程序。利用该程序模拟了某大涵道比风扇/增压级在周向总压畸变进气下的三维流场。计算结果表明：大涵道比风扇单转子不同叶高处的畸变传递特征差异较大,转子出口总压畸变强度由叶根到叶尖逐渐降低,在叶尖处衰减为最小值1.5%;在转子出口相应诱导出的总温畸变强度由叶根到叶尖逐渐升高,在叶尖处达到最大值1.4%;进口周向总压畸变导致风扇转子总压比下降0.5%,而风扇转子出口形成的总压总温复合畸变导致增压级总压比下降2%;总压畸变在增压级中呈逐级衰减趋势,而高温畸变区的周向范围在逐级增加。     

长拉杆-止口连接弯曲刚度损失及对转子系统振动响应影响

航空燃气轮机为了实现高负荷、轻质化的追求,在转子结构设计中,趋向于提高转速和加大长径比。这使得转子系统弯曲模态临界转速降低,转子在工作转速范围内不可避免会产生一定的弯曲变形。转子弯曲变形会影响连接结构界面接触特性的变化,使其连接结构局部弯曲刚度产生损失。因此,对于工作转速靠近弯曲临界的高速转子系统,需要考虑连接结构界面接触状态变化对转子系统振动特性的影响。以高负荷的长拉杆-止口连接转子系统为对象,分析连接界面接触应力分布特性,提出连接结构弯曲刚度损失修正方法,以此为基础建立界面连接转子动力学模型。通过对止口连接三级轴流压气机转子结构动力学特性的仿真和试验研究表明,在靠近弯曲振型临界转速下,转子连接界面接触状态的变化会产生弯曲刚度损失,对转子动力学特性具有显著影响。     

Maxwell方程的高阶间断有限元数值解法

采用高精度方法求解时域Maxwell方程,方程的空间离散采用基于计算流体力学(Computational fluid dynamics,CFD)领域的高阶间断有限元格式,非定常时间迭代采用四步龙格-库塔格式。为了提高计算效率,本文采用了Quadrature-free implementation和网格分区并行技术。数值结果表明,采用高阶格式的情况下,采用稀疏网格便可以得到高精度数值解。另外由于本文的方法基于非结构网格,因此非常适合计算复杂外形的情况。