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

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

光纤环及其骨架材料膨胀系数测量方法研究与实现

介绍了一种光纤环及其骨架膨胀系数的简易测量方法,并组建了测量系统。利用铝材和不锈钢,对系统的可行性进行了验证,同时用该系统测量了光纤环及不同骨架的膨胀系数,并对实验数据进行了分析研究,实验结果表明：该系统能够为找出与光纤环相匹配的骨架材料提供技术支撑,具有一定的实用和推广应用价值。     

一种膨胀映射的不动点定理

王尚志给出了若干种膨胀映射的概念，并证明了它们的不动点存在定理。对其中第一型膨胀映射的不动点定理做了进一步推广。     

正态分位点的矩阵算法浅析

假设标准正态的t分位点记为z(t),本文将根据z(t)的Taylor展开式中各个系数之间的矩阵递推关系,给出一种计算z(t)的近似值的方法,并给出具体的程序运算来验证本文的结论.     

可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展。但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题。界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因。为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案,本文从数值框架的选择、非守恒方程相容离散、高精度有界格式构造、界面压缩、界面-单介质分区计算方法等多个维度出发,综述近几年我们在可压缩多介质流动问题高精度数值模拟方法方面的研究进展。通过上述多个维度的工作,我们建立了一套适用于可压缩多介质流动问题的低耗散、基本无数值振荡的高精度欧拉数值方法,并成功应用于可压缩多介质大变形流动和界面不稳定性诱导湍流混合等问题的数值模拟。相关数值方法研究成果已集成至武器物理等领域工程数值模拟软件中,为相关工程任务提供了重要技术支撑。     

Analysis of the data processing strategies of spherical harmonic expansion model on global ionosphere mapping for moderate solar activity

Spherical harmonic (SH) expansion is widely used to model the global ionosphere map (GIM) of vertical total electron content (VTEC). According to the impact of different data processing methods of the SH expansion model on the VTEC maps, we specifically performed comprehensive analysis in terms of the data sampling rate, the time resolution, the spherical harmonic degree, and the relative constraint. One month of GPS data (January in 2016) from the International GNSS (Global Navigation Satellite System) Service (IGS) network in a moderate ionospheric activity period at the descending phase of Solar Cycle 24 was processed. To improve the computational efficiency of the daily GIM generation, the data sampling rate of 5?min was recommended allowing the GIM precision loss within 0.10 TECU (total electron content unit). The global VTEC map could be better represented in temporal and spatial domains with higher time resolution and higher spherical harmonic degree, especially at low latitude bands and in the southern hemisphere. The GIM precision improvement was about 10.91% for 1-h and about 15.15% for 0.5-h compared with the commonly used 2-h time resolution. The use of spherical harmonic degree 17 or 20 instead of 15 could improve the precision by 3.19% or 6.06%. We also found that an optimal relative constraint had to be found experimentally considering both the GIM precision and the GIM root mean square (RMS) map.     

失效概率矩独立全局灵敏度分析的高效算法

失效概率矩独立全局灵敏度分析对指导可靠性优化至关重要，本文从乘法降维及Edgeworth级数展开的角度提出一种失效概率矩独立全局灵敏度指标分析的高效算法。所提算法通过Edgeworth级数展开将失效概率矩独立全局灵敏度指标的求解转化为输出无条件及条件前四阶整数矩的求解。对于输出的无条件及条件四阶矩的计算，基于乘法降维的思想，本文推导了重复利用计算输出的无条件矩的输入输出样本来计算条件矩以及外层关于输入变量一维积分的计算公式，使得仅需重复利用输出无条件前四阶矩求解产生的积分网格内的数据即可同时求得输出的前四阶条件矩以及外层关于输入变量的一维积分。所提算法大大提高了失效概率矩独立全局灵敏度的分析效率。航空发动机涡轮盘以及汽车前轴的分析结果验证了所提算法的高效性及准确性。     

用粘性涡方法计算伞状物绕流

文章用粘性涡方法计算了二维伞状物绕流问题。首先阐述了涡方法的基本理论 ,包括涡元的对流、粘性扩散、物面涡元的脱落、物体所受气动力计算等 ,然后介绍了使用涡方法的一般步骤。最后将完全充满的伞衣简化为二维圆弧壳 ,计算了伞衣内外压力分布、阻力系数及升力系数     

基于数值模拟的孔构件挤压强化疲劳寿命预测

为研究飞机机身7050铝合金孔构件挤压强化后的疲劳增益,采用数值模拟与试验研究相结合的方法,对孔构件的挤压强化过程、疲劳加载过程、疲劳裂纹萌生和裂纹扩展过程进行研究。通过数值模拟探究了孔构件在不同状态下危险截面的应力分布和对应的疲劳行为,分析了残余应力场对疲劳性能的影响,探讨了残余应力与疲劳裂纹萌生和裂纹扩展的内在联系,建立了孔构件挤压强化疲劳寿命数值预测模型。结果表明：孔挤压强化引入的残余压应力可以减小孔构件在受载时孔壁最大拉应力,改变疲劳裂纹的萌生位置,抑制疲劳裂纹萌生和裂纹扩展,提高7050铝合金孔构件疲劳寿命近2倍,疲劳寿命数值预测模型误差在12%以内。     

间隙高度对超声速膨胀器流场及性能的影响

采用三维雷诺平均N-S方程和标准k-ε湍流模型,对不同间隙高度超声速膨胀器的流场和性能进行了数值研究,结果表明:间隙高度显著影响三维流道内的局部流动特性和超声速膨胀器的整体性能,随间隙高度增加,气流最高相对马赫数降低,高速区范围逐步缩小;泄漏涡增强,尺度变大,横向和径向运动明显,泄漏损失增加,但激波及激波附面层相互作用的损失降低;超声速膨胀器的膨胀比先增大后减小,等熵绝热效率持续降低.下端壁、吸力面附近低能流体之间以及与壁面的摩擦损失和间隙泄漏损失是有间隙超声速膨胀器三维流道内损失的主要来源,超声速膨胀器的间隙高度宜在0.9%h₀~1.5%h₀之间选取.