基于结构网格点的Euler方程无网格自适应求解算法

提出了一种自适应无网格Euler方程求解算法。构造了基于无网格点云结构的布点加密技术,实现了借助压力梯度变化准则的流场局部加密;借助流场结构化剖分技术,给出初始无网格节点分布,并采用四步显式Runge—Kutta时间格式推进计算,求解了绕翼型的典型流动问题。数值算例表明,本文提出的方法能有效地提高流场的分辨率,如捕捉激波等流动特征清晰明了。     

Modelling of low-energy galactic electrons in the heliosheath

The modulation of cosmic ray electrons in the heliosphere plays an important role in improving our understanding and assessment of the processes applicable to low-energy galactic electrons. A full three-dimensional numerical model based on Parker’s transport equation is used to study the modulation of 10 MeV galactic electrons, in particular inside the heliosheath. The emphasis is placed on the role that perpendicular diffusion plays in causing the extraordinary large increase in the observed intensities of these electrons in the heliosheath. The modelling is compared with observations of 6–14 MeV electrons from the Voyager 1 mission. Results are shown for the radial intensity profiles of these electrons, as well as the modulation effects of varying the extent of the heliosheath by changing the location of the termination shock and the heliopause and the value of the local interstellar spectrum. We confirm that the heliosheath acts as a modulation ‘barrier’ for low-energy galactic electrons. The significance of this result depends on how wide the inner heliosheath is; on how high the very local interstellar spectrum is at these low energies (E < 100 MeV) and on how small perpendicular diffusion is inside the inner heliosheath.     

非定常气动力对飞行动力学特性影响分析

在模糊逻辑非定常气动力建模的基础上,研究了非定常气动力对飞行动力学特性的影响.首先分析了非定常气动模型与定常气动模型的差别.其次采用分支分析方法,对带定常与非定常气动模型的飞行动力学系统的陡尾旋状态平衡分支图进行计算、比较与分析.并通过仿真对结果进行验证.研究表明,在大迎角状态下,非定常气动模型与定常气动模型存在较大的差别,对飞行动力学特性将产生一定的影响,需要在大迎角控制律设计中考虑非定常气动力的影响.     

考虑局部应力效应的轮盘子午面破裂转速预测

平均应力法用于航空发动机轮盘破裂转速预测时结果一般偏大,对于带螺栓孔的轮盘破裂转速预测误差更大。考虑局部应力效应对平均应力法进行修正,引入修正因子来表征轮盘子午面平均周向应力和危险部位局部应力对轮盘破裂转速的加权影响。通过开展典型轮盘材料GH4169的光滑试样和缺口试样拉伸试验,并开展了两件模拟盘破裂转速测量试验,得到修正因子与缺口强度比的关系,然后同时采用这两种方法对某涡轮盘的破裂转速进行了预测,并开展了验证试验。试验结果表明:平均应力法预测结果偏大,误差达到了9.4%,而考虑局部应力修正后,预测值仅略低于试验值,误差仅2.67%。考虑局部应力效应的轮盘子午面破裂转速预测方法简单可靠,具有较高的工程实用价值。      

小波变换在大地回波噪声处理中的应用

基于小波变换的多尺度分辨特性,作者对连续信号、脉冲信号和宽平稳白噪声在小波域中不同分解尺度上的模极大值变化规律进行分析,根据它们的变化规律用Mallat快速算法设计出一套用于抑制低信噪比雷达回波中的干扰噪声、最大限度提纯真实回波信号的算法.仿真实验证明,该算法对抑制噪声作用明显,能较好地提纯雷达回波信号,这为低信噪比下雷达回波信号检测提供了一条有效途径.     

轮胎刚度特性对大型民机前起落架摆振影响研究

轮胎刚度作为轮胎动力学模型的重要参数之一,研究其对大型民机前起落架摆振的影响规律,可以从前起落架防摆设计的角度为轮胎刚度设计提供参考依据。基于前起落架摆振非线性数学模型,使用Matlab/Matcont 软件计算不同轮胎扭转刚度、轮胎侧弯刚度下的摆振区域图;研究轮胎扭转刚度、轮胎侧弯刚度对前起落架摆振的影响规律,并对比二者对前起落架摆振影响的敏感度。结果表明：对前起落架摆振影响的敏感性,轮胎扭转刚度大于轮胎侧弯刚度;轮胎回正力矩系数每减小1%,扭转摆振最大临界阻尼减小0.88%,侧向摆振中速区最小临界减摆阻尼增大33.87%;减小轮胎扭转刚度,增大轮胎侧弯刚度有利于抑制大型民机前起落架摆振。     

周期激励Ueda电路中Hopf分岔与混沌的参数区域

周期激励Ｕｅｄａ电路是一种具有奇异混沌行为的强非线性电路。本文基于Ｈｏｐｆ分岔条件给出了一个确定混沌参数区域的方法。首先使用谐波平衡法获得基谐波幅度随时间变化的方程;然后由该方程平衡点的稳定性得出平衡点的Ｈｏｐｆ分岔出现条件;最后通过数值方法,计算出在Ｈｏｐｆ分岔曲线周围系统的Ｌａｙａｐｕｎｏｖ分量,确定Ｕｅｄａ电路出现奇异混沌吸引子的参数区域。结果表明,Ｕｅｄａ电路的Ｈｏｐｆ分岔曲线附近确实存在着奇异混沌吸收子,这为预测Ｕｅｄａ电路奇异混沌吸收子的出现提供了一种有效的方法。文末还讨论了混沌与Ｈｏｐｆ分岔之间的关系。     

三维凸域中候选卫星点集的确定算法

局部网格生成是无缝有限元并行计算的关键问题,然而,不恰当的局部网格生成算法会产生不一致网格.借助Delaunay三角剖分的唯一性,指出了三维凸域中Delaunay网格局部生成的不一致性,分析其产生的原因,建立避免局部网格不一致性的关键算法--搜索候选卫星点集的优化探索球算法,为网格并行生成做好数据准备.并从理论上证明了任意节点的优化探索球域包含其所有卫星点,从而,保证全域网格没有不一致性.通过算例分析了算法的正确性和局部性.     

不带稳定翼飞船返回舱俯仰动稳定性研究

联盟号飞船返回舱,如果不存在稳定翼,从高超声速到低Mach数时,力矩曲线有由一个平衡迎角(α1),经二个(a1、α'2)到三个(α1、α2、α3)平衡迎角的情况出现.当有一个平衡迎角时,可能有一临界Mach数Mcr,当M∞＞Mcr时,平衡迎角处是动稳定的;M∞＜Mcr,变为动不稳定的,同时出现极限环.也可能不存在这个Mcr,这时随Mach数降低到M∞=M'cr,即2个平衡迎角的情况,在力矩曲线与α轴相切的平衡点α'2,会出现鞍、结点分叉.M∞＜M'cr后,对于存在三个平衡点的情况,在α和α的相图上,它是结点-鞍点-结点的结构.当α1处是稳定结点形态时,只有在很小的扰动情况下,它才能返回α1的状态.如果α1为不稳定结点状态,在α1附近,微小的扰动,使迎角由α1迅速变成α3或者更大.这是混沌的先兆.     

自适应网格局部加密方法在可压缩流场模拟中的应用

本文基于M.J.Berger和J.Oliger发展的自适应网格局部加密算法,它采用多个组合网格的思想,用于求解双曲型方程组。它结合显式有限差分格式,采用Richardson外推技术自动进行局部截断误差估计,对精度低的区域产生新细网格局部加密,或去除无需再加密的旧细网格,以最小的运算量达到给定的精度要求。网格可一层层加密下去,按层次覆盖,每个网格是有任意方向的矩形均匀网格。这套算法独立于求解所用的差分格式,很容易和各种格式结合。 我们从M.J.Berger和J.Oliger程序出发,参考M.J.Berger和P.Colella,对守恒型方程组,实现了粗细网格交接面上数值通量守恒,从而能计算间断解。对网格边界处理、误差估计和网格重新生成等做了一些改进,并且可以灵活地在任何感兴趣区域指定加密。 我们用不定常Euler方程组计算了带前台阶的二维管道中Mach-3流动,采用MacCormack显式附加人工粘性差分格式。共用了四层网格结构:Δ=1/10、1/20、1/80和1/320。其中基本网格是Δ=1/10;到台阶上第一次激波反射区域从Δ=1/20加密至Δ=1/80,计算到定常所需运算量仅是同等均匀网格△=1/80的28％;我们还用Δ=1/320网格对台阶拐角和切向间断处指定加密,结果表明可明显改善台阶上的激波Mach反射,切向间断清晰可辨,所需的计算量是上述均匀网格Δ=1/