简介:文章通过对EFM(effectivefieldmodeling)模型进行简化,消除了原模型的非守恒性项和非双曲性特性项,发展了一种基于密度的气液两相流模拟方法:ρ-VOF方法.利用体积分数信息对控制单元内的自由界面进行重构,得到了控制单元内流体的空间分布,并采用AUSM^+-up格式获得考虑气液流体接触间断信息的对流通量.新方法可统一处理激波间断和接触间断的相互作用,保持自由界面的尖锐性,并且其计算量与自由界面的空间复杂度无关.最后,数值模拟了液体激波管气液激波管和气体激波跨二维液滴传播等问题,并与文献结果进行对比,验证了本方法在气液两相流模拟中的准确性.
简介:由大粗糙元引起的髙超声速边界层强制转捩在航天技术中有实际应用,因而近年来受到人们的广泛关注.虽然目前导致该转捩过程的内在机理尚不完全清楚,但有一点是明确的,即粗糙元的尾迹流场中存在强对流不稳定性.文章的出发点是研究这种对流不稳定模态是如何触发转捩的.首先通过CFD方法,计算出髙超声速边界层中粗糙元的尾迹流场,并对其进行二维稳定性分析.结果发现,在传统不稳定Tollmien-Schlichting(T-S)模态出现的临界Reynolds数之前,存在髙增长率的无黏不稳定模态,表现为对称的余弦模态和反对称的正弦模态.然后对该不稳定模态在粗糙元尾迹流中的演化进行了模拟,验证了二维稳定性分析的结果,并考察了非平行性效应的影响.最后通过直接数值模拟,研究由这些不稳定模态触发转捩的全过程.结果表明,对流不稳定模态确实是导致边界层转捩的关键机制.该转捩过程的特点是,局部湍斑首先在不稳定模态特征函数的峰值附近出现,然后向全流场扩散.就文章研究的工况而言,余弦和正弦模态的相互作用对转捩的影响并不明显,且后者在转捩过程中起主导作用.
简介:文章给出了一种真正多维的HLLRiemann解算器.采用AUSM分裂将通量分解成为对流通量和压力通量,其中对流通量的计算采用迎风格式,压力通量的计算采用HLL格式,且将HLL格式的耗散项中的密度差用压力差代替,从而使得格式能够分辨接触间断.为了实现数值格式真正多维的特性,分别计算了网格界面中点和角点上的数值通量,并且采用Simpson公式加权组合中点和角点上的数值通量得到网格界面的数值通量.为了减少重构角点处状态时的模板宽度,计算中采用基于SDWLS梯度的线性重构获得2阶空间精度,而时间离散采用2阶保强稳Runge—Kutta方法.数值实验表明,相比于传统的一维HLL格式,文章的真正多维HLL格式具有能够分辨接触间断,以及更大的时间步长等优点.与其他能够分辨接触间断的格式(例如HLLC格式)不同,真正多维的HLL格式在计算二维问题时不会出现激波不稳定现象.
简介:使用高阶间断Galerkin(discontinuousGalerkin,DG)方法求解双曲守恒律方程组时,非物理效应常常导致计算过程的中断,这在很大程度上制约着该方法在计算流体力学中的应用.文章结合局部单元上原始流动变量的Taylor展开,设计了一种新型的限制器,通过对各阶空间导数的重构,有效地消除了非物理振荡的不利影响.对二维Euler方程的计算结果表明,该限制器不仅能够捕捉高质量的激波,而且能够保证残值的有效收敛.
简介:基于状态空间模型的许多传统滤波算法都基于Rn空间中的高斯分布模型,但当状态向量中包含角变量或方向变量时,难以达到理想的效果。针对J.T.Horwood等提出的nS?R流形上的GaussVonMises(GVM)多变量概率密度分布,扩展了狄拉克混合逼近方法,给出了联合分布的GVM逼近方法,推导了后验分布的GVM参数计算公式,设计了量测更新状态估计算法。将J.T.Horwood等的时间更新算法与所提出的量测更新算法相结合,可实现基于GVM分布的递推贝叶斯滤波器(GVMF)。仿真结果表明,当状态向量符合GVM概率分布模型时,GVMF对角变量的估计明显优于传统的扩展卡尔曼滤波器。
简介:开式凹腔作为超燃冲压发动机中增加掺混和稳焰的装置,其流动稳定性的研究对深入理解凹腔增加掺混和稳焰机理以及凹腔的设计有着重要的学术意义和工程应用价值.基于大涡模拟方法对超燃冲压发动机开式凹腔流动进行数值模拟,分别米用动力学模态分解(dynamicmodedecomposition,DMD)和本征正交分解方法(properorthogonaldecomposition,POD)对自激振荡流动进行稳定性分析.DMD方法可准确提取凹腔的振荡频率,与Rossitei'模型以及压力脉动FFT分析得到的频率吻合较好,且DMD中对应Roster前3阶频率的模态在流动中的主导作用顺序也与FFT分析结果一致,自激振荡中RossiterH模态占据主导作用,同时DMD方法对Rossiter3阶以上模态频率的预测能力明显强于FFT分析方法.在对低频的提取方面,DMD方法比Rossiter模型更具有优势.与前6阶Rossiter模态对应DMD模态均缓慢收敛,主要表现为剪切层中的分离涡结构和中部及下游区域中的涡结构.前3阶不稳定模态中的分离涡结构主要集中在中部剪切层以及后缘附近区域.POD方法中较少的模态包含流场绝大部分的能量.但是,通过POD方法提取的模态频率在分辨率上效果不佳,提取到最低频率为Rossiter3阶模态对应的频率,且模态中均存在次频,次频与主频之间的耦合导致模态的形态相差较大.另外,与DMD方法相比POD方法无法判断所提取的模态的稳定性.
简介:针对传统基于g信息的粗对准的捷联惯导系统中,受传感器噪声的影响,存在效视运动无法提取和双向量共线的缺点,提出了一种基于改良Kalman滤波的参数辨识粗对准方法。该方法通过构建视在重力在初始载体系中的映射模型,利用改良Kalman滤波进行模型参数辨识,然后通过识别参数重新构建视在重力在初始载体系中的映射,解决了由于传感器噪声导致有效视运动无法正常提取的缺点。利用识别参数具有随估计次数增多得到优化的特点,构造初始时刻和最终时刻向量,避免双向量共线问题。利用改良Kalman滤波算法的自适应特点,优化参数识别精度与速度。转台实验表明,采用改良Kalman滤波方法航向对准精度为-0.0414°,标准差为0.041°,而传统RLS方法得到的航向精度为-0.0738°,标准差为0.128°。由此可知,本文提出的方法性能更优。
简介:由条带和流向涡的循环再生构成的近壁自维持过程(self-sustainingprocess,SSP)是壁湍流产生和维持的重要机制.文章通过对最小槽道的直接数值模拟(directnumericalsimulation,DNS)获得近壁自维持过程的流场数据,采用正规正交分解法(properorthogonaldecomposition,POD)对该数据进行分析,获得了不同流向和展向尺度的特征模态,通过将Navier—Stokes方程在这些模态上进行投影,得到近壁自维持过程的降阶模型,并采用DNS数据对降阶模型的预测能力进行了评价.该模型被初步应用于大涡模拟近壁模型的构造.
简介:为提高光电平台的控制性能和稳定性,以平台反馈回路所用的光纤陀螺传感器为研究对象,对光纤陀螺角速率的历史输出、当前量测以及随机漂移进行融合补偿。采用双自回归模型确定了光纤陀螺时间序列输出的自回归多项式和光纤陀螺随机漂移的自回归关系。以陀螺当前输出为量测量,结合卡尔曼滤波算法将陀螺历史输出和历史随机漂移融合进状态方程,并进行随机漂移在线估计补偿。实验结果表明,光纤陀螺随机漂移的AR模型能达到90%拟合效果,经卡尔曼滤波补偿后随机漂移能降到1/10。该方法能很好地抑制光电平台三个框架轴光纤陀螺的随机漂移,补偿率为80%~90%。
简介:针对带有攻击角约束的多导弹同时攻击机动目标问题,提出了一种带有攻击角约束的协同制导律。首先基于平面内的导弹-目标相对运动方程,建立了带有攻击角约束的协同制导模型;其次,把协同制导律的设计过程分离为两个部分:一是基于图论的有关内容,运用有限时间一致性理论设计沿着视线方向上的加速度指令来保证所有导弹与目标的相对距离在有限时间内到达一致,进而保证所有的导弹同时击中机动目标;二是利用非齐次干扰观测器对机动目标的加速度进行估计,并运用滑模控制设计视线法向上的加速度指令来保证每枚导弹与目标间的视线角速率收敛到零和视线角收敛到期望的终端视线角,即每枚导弹以期望的终端视线角成功击中目标;最后,对三枚导弹同时打击同一机动目标的情况进行仿真,仿真结果表明本文设计的带有攻击角约束的协同制导律的有效性和正确性。
简介:在飞行器的气动外形优化设计中,参数化方法和优化算法具有十分重要的作用,对优化的计算时间、设计空间的数学特性有着深刻的影响.类别形状函数(classandshapetransformation,CST)方法是一种简洁高效的参数化方法,但对于复杂曲面很难使用统一的CST方法进行拟合.文章首先介绍了CST方法的三维实现,分析了其数学性质,提出了分块CST参数化方法,保留CST方法的特性,实现了分块曲面之间的光滑连接.针对气动外形优化设计的复杂情况,需要根据具体的飞行任务提出设计目标,并处理不同目标的矛盾问题.其次采用Pareto策略自动寻找最优方案集,并基于分块CST参数化方法、遗传算法和气动力快速计算方法,对类乘波翼身组合飞行器进行了优化设计,并改变原有问题的设定条件优化得到了全新外形.研究结果表明分块CST方法参数少,精度高,Pareto策略处理多目标准确有效,是气动外形优化设计中非常有用的工具.
简介:大视角图像匹配算法的鲁棒性与实时性直接影响飞行器对远距离目标定位的性能。针对目前仿射不变图像匹配算法实时性较差的问题,提出一种惯性信息辅助的快速大视角图像匹配方法。该方法对现有的快速图像匹配算法进行改进,避免了构建高斯金字塔,提高了算法效率。然后利用机载惯性导航信息求解实时图与参考图之间的单应性矩阵,并对实时图进行模拟视角变换以此减小图像间视角差异,克服了现有的大视角图像匹配算法盲目多次的匹配计算,实现了大视角图像的快速匹配。实验结果表明,惯性信息辅助的大视角图像匹配算法与现有的快速仿射不变性匹配算法相比,匹配效率提高了至少2倍。