高频不稳定燃烧的声学数值仿真

高频不稳定燃烧一直是液体火箭发动机研制过程中所要面临的重大难题之一.采用具有低色散低耗散特点的计算气动声学方法,对自燃推进剂变轨发动机(OME)的高频不稳定燃烧进行时域下的数值仿真.由Crocco的压力时滞模型对燃烧热释放和声波之间的耦合进行模拟,并对不同的时滞模型参数对稳定性结果的影响进行了分析,给出发动机的稳定性极限图,确定一阶切向及一阶径向振型为主的不稳定振型,与地面试车实验捕捉到的不稳定振型相一致.结果表明:采用计算气动声学方法对带有Crocco压力时滞模型的声波扰动方程进行时域下的数值求解,可以对发动机的高频不稳定燃烧进行成功地预测.      

基于非湍流脉动动能方程的高超声速转捩预测

通过求解基于稳定性理论的非湍流脉动动能输运方程,预测了来流Ma为5.91的裙锥绕流流动的层流-湍流转捩点,并采用SST(Shear-Stress Transport)湍流模型结合代数间歇因子模型对转捩流场进行了数值模拟.结果表明非湍流脉动动能输运方程模型可以很好地捕捉逆压梯度对转捩前期边界层内的不稳定波动频率的影响.等温壁面条件下计算得到的转捩位置与静风洞实验结果基本一致.壁面冷却促进第2模态不稳定波动的增长,使转捩提前发生.代数间歇因子模型模拟高马赫数流动的转捩区长度较短,温度峰值偏低.     

飞轮控制的欠驱动刚体航天器姿态控制器设计

针对"三正交加斜装"反作用轮系统中某两个本体轴上的飞轮失效的欠驱动情况,研究了航天器的姿态控制问题.在系统初始角动量为零的条件下,设计分段解耦控制律,实现了姿态稳定.采用欧拉角描述法建立了欠驱动航天器的姿态动力学方程和运动学方程.在系统初始角动量为零的条件下,通过分析方程的解耦特性,设计了分段解耦控制律.该方法经过6次机动控制,可实现姿态稳定.数值仿真验证了方法的有效性.     

卫星定位和精度因子的改进方法

在卫星导航定位系统中,在精度因子计算和采用最小二乘法进行定位求解时,传统上采用测量矩阵直接求逆方法来进行.为了克服矩阵求逆带来的计算量大和数值稳定性差的不足,利用测量矩阵的对称正定性,提出了一种基于矩阵 UTDU 分解的定位解算和精度因子计算方法.改进方法具有严格的数学理论基础,保证了方法的正确性和有效性.数值分析结果表明,相对直接求逆的传统方法而言,在定位解算时,该方法能降低约 60%的运算量,而在精度因子计算中,约能降低36%的运算量.且改进方法能大大降低求解矩阵的条件数,提高了求解的数值稳定性.      

基于天线驱动组件的多体卫星MIMO控制及稳定裕度应用研究

针对两输入两输出的多体卫星控制系统,采用工程中的天线驱动组件为天线转动提供驱动力矩,在天线小角度机动条件下,对此线性模型采用多变量频域设计法中的序列回差设计技术,设计卫星本体俯仰姿态和天线转动的多输入多输出(MIMO)控制器,利用Nyquist稳定判据证明系统的稳定性.在此基础上分别利用回差阵奇异值法和逆回差阵奇异值法计算MIMO控制系统的稳定裕度,将两种方法所得结果相结合以减小稳定裕度计算的保守性,最后通过数值仿真验证了控制算法和控制系统稳定裕度结果的有效性.     

近壁区湍能耗散率输运方程的理论模型

提出在湍流边界层的近壁区采用三维不稳定波来描述湍流扰动速度,然后根据理论模型对湍能耗散率的生成、粘性破坏以及扩散等特性进行了系统定量的计算和分析.理论计算结果与直接数值模拟结果符合较好,特别是在y<10的区域中明显优于传统湍流模型的计算结果.这不仅在理论上有益于对湍流物理机制的了解,而且可能为湍流模型计算开辟一条新的途径.     

气动谐振管加热机理数值研究

放置在高速喷流中的半开口气动谐振管在一定的条件下会出现加热现象.针对这一现象,利用数值计算方法计算了一维非定常非线性Euler方程组,采用MacCormark格式,对不同管口条件下谐振管内的振荡过程进行了数值研究.研究发现,当管口压力变化频率处在谐振管的线性谐振频率附近时,管内形成的激波以及激波的反射增强导致了谐振管端壁区域明显的加热效果.对不同的驱动气体和不同的谐振管型式(圆柱形和圆锥形)的数值研究也给出了有价值的结论.研究结果与实验事实或现象是相吻合的.      

捷联式天线稳定平台动力学建模与仿真分析

对于体积有限制的战术导弹应用,导引头天线伺服系统的体积受到严格限制,传统的速率陀螺稳定平台的应用受到挑战,为此提出采用捷联稳定方式来解决视线稳定问题.根据弹载捷联式天线平台的结构特点及稳定原理,应用坐标变换方法,完整地建立了两自由度捷联式天线平台的运动学与动力学模型.通过分析平台框架运动学耦合关系,表明捷联稳定方式与速率陀螺稳定方式的基本区别在于信息的获取和控制方式上,速率陀螺稳定方式是一种直接硬件稳定方法,而捷联稳定方式是一种软件补偿方法.在此基础上,对框架动力学耦合特性进行了数值仿真,结果表明:弹体对框架的耦合较强,需要进行弹体解耦,而框架之间的交叉耦合相对较弱,实际设计时可以忽略.所得结果为进一步研究捷联式天线稳定平台提供了理论基础,也为小型化导引头视线稳定平台系统的工程设计提供了一种解决方案.     

交替LU分裂算法及其在CFD中的应用

在CFD(Computational Fluid Dynamics)时间相关算法中,为了保证计算的稳定性,时间步长的取值通常会很小,这将导致计算过程收敛缓慢.针对这一问题,提出了一种新的迭代算法—交替LU分裂(ALUS,Alternating Lower-Upper Splitting)算法,可以有效加速收敛,提高计算效率.ALUS算法将系数矩阵分裂成上、下三角矩阵,因此仅需要利用追赶法求解两个三角矩阵,计算量较小,容易实现.给出了ALUS算法收敛的定理,并且通过线性问题以及CFD圆柱绕流的数值模拟对ALUS算法进行了检验.理论分析和数值实验的结果均表明:ALUS算法计算量小,大大节省了计算时间,而且该算法是鲁棒的.因此ALUS算法是高效的、稳定的算法,适用于CFD数值模拟.     

背景太阳风数值模拟的熵守恒格式

背景太阳风研究是根据行星际扰动传播情况预测空间天气状况的基础,磁流体（MHD）模拟是背景太阳风研究的重要手段.采用一种新的数值计算方式,利用Ideal GLM-MHD将计算过程中产生的磁场散度以c_h的速度向计算区域外传播,从而消去磁场散度;重构部分使用受约束的最小二乘法,将磁场散度作为约束条件添加到重构中,进一步对重构后的磁场梯度进行修正;通量计算采用满足热力学第二定律的熵守恒格式,该格式能够确保在计算过程中熵不增,保证数值稳定.研究结果表明,该方法应用于太阳风数值模拟的求解得到了更加稳定的结果.      

基于模糊理论的加速度计标度因数稳定性分析

运用模糊理论描述了影响加速度计标度因数稳定性诸多因素的不确定性问题,针对其中的磁钢磁感应强度和石英摆片热膨胀系数两类稳定性影响因素,给出了其三角模糊数区间,并通过该型加速度计的磁-结构耦合ANSYS分析求出了在模糊参数不同组合下的仿真输出,根据仿真输出结果获得了标度因数模糊稳定区间,给出了加速度计标度因数稳定性分析流程及实例.分析说明了基于模糊理论的稳定性分析方法可以为加速度计加速稳定试验截止条件提供依据,对提高加速度计标度因数长期稳定性水平具有重要意义.     

高温化学非平衡计算中源项处理方法

高超声速三维化学非平衡绕流流场的数值模拟中,在计算初期容易发生组元密度出现负值的非物理现象,另外源项的刚性是影响计算稳定性和收敛速度的主要原因。根据源项线性化理论,采用2种源项线性化方法处理化学反应源项.一种为根据流场内化学反应物理规律构造的线性化方法,该方法能抑制在计算过程中组元密度出现负值现象,提高了计算稳定性,加速了收敛速度.同时采用了时间预处理矩阵的线性化方法,较好地解决了非平衡化学反应与流场耦合的刚性问题。证明了所构造源项线性化方法相容性.数值实验表明该方法有效地避免了源项计算中密度出负的问题,加快了计算的收敛速度,从而提高了计算效率.      

光源光谱对干涉式光纤陀螺零漂的影响

以采用Y波导集成光学调制器的保偏型干涉式光纤陀螺(IFOG,Interference Fiber Optic Gyroscopes)为研究对象,根据各光学元器件的参数,建立了各器件的琼斯矩阵以及光路传输模型,在此基础上进行了光路偏振误差的理论分析.通过推导,得到了保偏型干涉式光纤陀螺的偏振误差表达式.首次借助光源偏振度的改变,分析了不同形状的光源光谱对偏振模式耦合误差引起的陀螺零漂的影响.通过计算得出,当光源的偏振度在0~3%之间呈线性变化以及光路中其它参数不变时,谱宽较窄的光谱引起的零漂较大,谱宽较宽的光谱引起的零漂较小.仿真结果给光纤陀螺光源光谱的选择提供了理论依据.      

星敏感器光轴热稳定性仿真分析方法

热稳定性用来衡量星敏感器光轴指向抵抗温度变化影响的能力,常用试验及仿真的方法进行分析验证.仿真方法因其高效率、低成本、可贯穿产品设计全周期等优点,已成为热稳定性分析验证的重要手段.开发了一种热 力 光联合仿真分析的方法进行星敏感器热稳定性计算.该方法集合了热分析、变形分析、成像仿真及热漂移计算等多项仿真分析过程.以某高精度星敏感器为例,进行了热稳定性计算,结果表明了该仿真方法的有效性.     

三阶中心无振荡格式性能分析

三阶中心无振荡格式Cn3(Centered nonoscillatory scheme of third order)使用对称模板构造具有三阶精度的插值公式.利用单调区域与精确区域修正原始插值参数,该格式能够获得间断附近无振荡、光滑区域高精度的计算结果.通过一维和二维典型算例,将Cn3格式与三阶和五阶WENO格式(Weighted Essentially Non-Oscillatory schemes)进行比较,重点分析了3种格式的间断分辨率、计算稳定性和数值耗散性.分析可见,Cn3格式能够精确、稳定地捕捉激波和接触间断,同时对光滑流动区域的小尺度流动结构保持较低的耗散,值得进一步研究及推广应用.     

桨根柔性无铰旋翼桨叶气弹稳定性建模分析

建立了任意剖面无铰旋翼桨叶气弹稳定性分析模型,分析了复合材料柔性梁对无铰旋翼桨叶气弹稳定性的影响.将旋翼桨叶简化为二维截面特性线性分析和一维梁非线性分析,二维特性分析考虑了面内和面外翘曲,一维分析采用中等变形梁理论,考虑了剪切变形和与扭转相关的翘曲.采用准定常气动力,运用有限元理论,得出了无铰旋翼桨叶稳定性分析的有限元列式.研究了柔性梁不同铺层方式对气动弹性稳定性的影响.算例表明复合材料柔性梁铺层角对旋翼桨叶气动弹性稳定性影响明显,充分地利用这些影响变化规律,能够设计出更为先进的复合材料悬翼桨叶.      

人工释能法求解线性代数方程组

本文建立求解线性代数方程组的人工释能法迭代格式,并给出该迭代格式收敛的充分必要条件和最佳释能时机,同时给出实现人工释能法的数值计算方法及其数值稳定性条件。     

垂直下降状态下的旋翼三维流场数值模拟

应用激励盘模型,采用计算流体力学(CFD, Computational Fluid Dynamics)方法结合动量叶素理论,计算了旋翼在垂直下降时的流场特征与气动性能.计算过程中,旋翼载荷被描绘成沿桨盘径向分布、与当地流动参数相关的压力源项.在直角坐标系中,运用有限体积方法直接求解了定常不可压N-S(Navier-Stokes)方程,湍流模型为S-A(Spalart-Allmaras)一方程模型,重点使用诱导速度经验公式计算出了压力源项.旋翼流场模拟结果和性能预计结果同实验测量数据吻合较好,表明这种CFD方法可以有效地分析旋翼下降飞行时的空气流动特性,为进一步的涡环状态研究提供了参考.