一类TVD型组合差分方法及其在磁流体数值计算中的应用

根据太阳风数值模拟的特点，考虑到算法的质量（收敛速度、稳定性、精度等），结合磁流体数值计算的特性，对三维球坐标磁流体动力学（MHD）方程组中的流体部分采用一种修正Lax-Friedrichs差分法而对磁场部分采用MacComack格式，发展了一类快捷的具有TVD特性的组合数值新方法，作为格式的检验，在一维情况下，将其与PPM格式进行了比较，对一维快慢磁流体激波问题得到了与PPM格式精度相同的结果，然后将其诮到定态太阳风的数值模拟上，在不同等离子体β情形下，可得到理想的太阳风定态结构，为今后将此数值模式应用到具有复杂磁场位型或三维直实太阳风暴的数值模拟研究奠定了基础。     

近空间高超声速流场通量分裂型DSMC-IP方法

针对近空间高超声速流场的特点,采用Van Leer格式对直接模拟蒙特卡洛-信息保存(DSMC-IP)方法质量守恒方程中的计算格式进行改进.将局部马赫数作为分裂通量的标准,并重构单元分界面两侧的左右输运通量,使得计算格式具有通量分裂的特点,解决了DSMC-IP方法在高超声速流场计算中的应用问题.将改进后的通量分裂型DSMC-IP方法引入非结构网格中,对二维近空间高超声速流场进行数值模拟.计算结果表明:通量分裂型DSMC-IP方法所得出的数值结果与实验值及参考值符合较好,明显降低了直接模拟蒙特卡洛(DSMC)方法所带来的统计耗散.当来流气体的稀薄程度增加时,其非平衡效应也更加明显,而通量分裂型DSMC-IP方法的计算结果与参考值相差均在10%以内,良好地反映了非平衡条件下的流场特征,验证了通量分裂型DSMC-IP方法的可行性和有效性.      

USB测距测速数据误差建模和评估

统一S波段（UnifiedS-band,USB）系统是目前中国载人航天和月球探测任务的主要测控网。文章对中国USB系统的测距测速数据的误差构成进行了详细分析,给出了传播媒介、相位中心偏差、时标偏差、通道时延以及卫星星历误差的修正模型,并通过对实测数据的残差分析,对目前中国USB测量数据的精度进行了评估,给出了USB测距和测速数据的实际精度。最后,针对目前USB数据的处理和应用情况,提出了数据处理上需注意的一些问题。文章的结论对今后改善中国USB系统的测定轨精度具有一定的参考价值。     

重气体介质中超临界翼型跨声速流动特性

采用Peng-Robinson非理想气体状态方程模拟重气体介质的热力学特性,并与雷诺平均Navier-Stokes方程结合,形成封闭的重气体介质流动模型。针对超临界翼型流动问题,利用LU-SGS（lower-upper symmetric Gauss-Seidel）隐式时间推进格式和有限体积法,分别求解空气介质和重气体介质下的流动特性。数值模拟结果表明:在跨声速条件下重气体介质中超临界翼型的升阻力增大、超声速区域表面负压增加、边界层位移厚度减小、激波后移、表面摩擦阻力明显增大、后缘流动分离推迟。该研究为后续重气体介质中飞行器颤振特性研究及修正方法的发展提供了基础支持。      

高精度多维限制器的性能分析

目前常用的限制器大都是基于一维构造,无法在多维情况下保证物理量的单调特性进而导致非物理振荡.为弥补传统方法的这一构造缺陷,多维限制器(MLP)通过多维修正使单元通量值介于周围相邻单元通量的最大值和最小值之间,在保证求解精度的情况下有效避免了多维振荡.基于一维激波管、无黏涡及激波边界层干扰等算例,对高精度MLP的特性进行了研究分析.结果显示:3阶MLP在连续和间断区域均可有效地避免多维振荡;与高阶WENO(Weighted Essentially Non-Oscillatory)方法相比,3阶MLP不仅算法简单、易于实现,还可显著提高求解的精度、保单调性及收敛性.因此可用于工程及科学研究的复杂流动,具有较好的应用前景.      

基于CO2浓度监测的飞机火警探测方法

目前飞机普遍采用感温、感烟火警探测方法,其存在误报率高的缺陷.分析了发生火灾时气态燃烧产物的生成规律,提出了以CO_2气体为火警探测参量,通过准确监测环境中CO_2气体浓度的变化,进行飞机火警探测的方法;根据飞机飞行所处环境条件的特殊性,研究了基于非色散红外(NDIR)吸收原理的CO_2气体浓度监测的飞机火警探测方法应用于飞机上所面临的温度补偿难题;根据飞机火警探测器测量电压与环境温度、参比电压之间的关系,采用偏最小二乘(PLS)法建立飞机火警探测器温度补偿模型;根据飞机火警探测器测量电压与CO_2气体浓度之间的关系,采用PLS法建立CO_2气体浓度计算模型;研究了基于CO_2气体浓度监测的报警算法,为飞机火警探测提供了一种新的方法.     

水下超声速气体射流线性稳定性研究

水下超音速气体射流稳定性是水下推进、水下焊接技术研究的焦点领域。通过建立线化小扰动气液混合流体控制方程,开展水下超音速气体射流线性稳定性研究。控制方程的形式表明射流基本流分布(速度与浓度分布)、雷诺数、无量纲密度、无量纲粘度能够影响水下超音速气体射流稳定性。使用Chebyshev配置点法对控制方程进行求解。计算结果表明射流基本流分布以及无量纲密度是影响射流稳定性的主要因素。射流基本流分布越平缓,射流最大扰动增长率越小;无量纲密度越小,即气液密度差越小,射流最大扰动增长率越小。随雷诺数增加,射流最大扰动增长率先减小后增加,当雷诺数大于105时,射流最大扰动增长率趋于常数。对不同无量纲粘度,射流扰动增长率-特征波数曲线高度重合。同时分析了射流不同截面的稳定性特征,计算结果表明射流最大扰动增长率随距喷口距离增加而减小。这些结论有助于理解水下超音速射流物理过程,为实际应用提供理论指导。     

间断Galerkin有限元隐式算法GPU并行化研究

为了提高间断伽辽金（discontinuous Galerkin, DG）有限元方法的计算效率,围绕求解Euler方程,构建了基于图形处理器（graphics processing unit, GPU）并行加速的隐式DG算法。算法结合Roe格式进行空间离散,采用人工黏性法处理激波等间断问题,时间推进选用下上对称高斯-赛德尔（lower-upper symmetric Gauss-Seidel, LU-SGS）隐式格式。为了克服传统隐式格式固有的数据关联依赖问题,借助于本文提出的面向任意网格的单元着色分组技术,先给出了LUSGS隐式格式的并行化改造,使得隐式时间推进能按颜色组别依次并行,由于同一颜色组内算法已不存在数据关联,可以据此实现并行化。在此基础上,再结合DG算法局部紧致等特点,基于统一计算设备架构（compute unified device architecture, CUDA）编程模型,设计了依据单元的核函数,并构建了对应的线程与数据结构,给出了DG有限元隐式GPU并行算法。最后,发展的算法通过了多个二维和三维典型流动算例考核与性能测试,展示出隐式算法GPU加速的效果,且获得...     

轴对称Boltzmann模型方程统一算法与空天飞行环境喷管流动模拟

拦截机动飞行器周围大范围区域存在主喷/侧喷流/羽流影响,而传统的纳维-斯托克斯（Navier-Stokes, N-S）方程不能很好模拟发动机喷管扩张段出口附近流动情况,需要一种新的方法来处理这种全流域流动问题。为解决该问题,针对特定轴对称喷管内流动,本文通过数学推导确立描述不同克努森数稀薄环境条件下的轴对称喷管内流动Boltzmann模型方程,初步建立适于该模型方程的数值格式与气体动理论统一算法。通过开展同轴圆筒间的定常/非定常旋转流动以及轴对称喷管内流动数值计算研究,发现统一算法计算流场与其他途径得到的结果吻合较好,验证了统一算法在全局克努森数喷管流动模拟的适应性和可靠性。通过与低密度风洞实验对比,喷管出口核心区羽流结构一致,羽流轴线压力分布一致,表明统一算法可以有效解决喷管入口压缩段到扩张段多流域混合,尤其是出口附近稀薄气体真空低压环境流动问题。