重力波波包在时变背景场中传播的数值研究

采用二阶算子分裂格式对高斯型重力波波包在时变背景温场和时变背景风场中的非线性传播过程分别进行了数值模拟和分析.研究表明,时变背景场对波包空间结构演变影响不明显,各物理量的偏振关系与线性理论偏振关系较一致;在时变背景温场中传播,水平群速度和垂直群速度的变化趋势与线性理论结果相反,波频率的变化也与线性理论色散方程计算结果有一定的差别,这表明了线性重力波理论的局限性;在时变背景风场中传播,波频率的变化比时变背景温场中更显著.传播路径与在稳恒背景场中比较,更接近线性理论计算出的路径.      

液体中心式同轴离心喷嘴液膜破碎特性仿真研究（液体火箭推进）

为了实现对不同工况下液体中心式同轴离心喷嘴液膜破碎特性的数值模拟研究,采用网格自适应加密技术、耦合的Level-set和Volume of Fluid（CLSVOF）方法对气液界面进行捕捉,利用改进延迟分离涡模拟（IDDES）方法模拟湍流。分析了液膜的破碎模式、喷雾锥角、破碎长度以及流场特性。通过观察分析得到：随着气液相互作用的增强,液膜破碎模式依次经历主导表面波发展导致的液膜破碎、Rayleigh-Taylor（R-T）和Kelvin-Helmholtz（K-H）不稳定性引起的液膜破碎,以及气动破碎模式。随着气液动量通量比（Momentum Flux Ratio,MFR）的增大,喷雾锥角和破碎长度逐渐减小且呈渐进趋势,发现无量纲喷雾锥角和破碎长度均与MFR_^-A成正比例关系。相同液膜破碎模式而不同工况时,主要流场特征一致。     

使用高阶逆风通量差分裂格式的悬停旋翼流场数值模拟

为减少Jameson二阶中心差分有限体积法导致的旋翼尾迹的数值耗散,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,建立了一个三维雷诺平均N-S方程数值模拟悬停状态旋翼流场的方法。为了充分考虑旋翼尾迹对流场的影响,采用周期性边界条件和由动量理论导出的远场边界条件。为进一步减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。然后,进行了算例计算,给出了桨叶表面的压强分布,与可得到的试验数据及二阶中心差分方法的计算结果进行了对比,并针对几种不同桨尖形状的旋翼悬停流场进行计算,数值结果显示:后掠桨尖可减弱超临界流动。此外,还计算和分析了旋翼下方不同位置上的涡量分布。计算结果表明,本文方法能够有效地减少旋翼尾迹的数值耗散。      

金属振动陀螺调平方案设计与比较

金属振动陀螺是一种低成本、 轻小型的新型固体波动陀螺,在战术级应用领域具有广大的应用前景.金属振动陀螺谐振子的频率分裂直接反映陀螺的性能指标,频率分裂可以通过机械调平的方式进行修正.对金属振动陀螺的调平方法进行了梳理和比较,提出了质量修正、刚度修正2种调平思路和增加、 去除材料等5种修正方法.通过对各种调平修正方法的比较,选择激光去重法对谐振子顶面进行质量调平修正,并进行实验验证.实验结果表明,该方法修正后,在保持品质因数基本不变的情况下,谐振子的频率分裂由5.3Hz降低到0.9Hz,陀螺的零偏稳定性由55(°)/h降低到6(°)/h,性能指标提高了1个数量级.     

基于WENO-分段线性格式的直升机流场数值模拟

发展了一套适合于格心格式求解器的基于加权本质无振荡(WENO)-分段线性格式的旋翼/机身气动干扰高精度CFD计算方法。应用该方法对多种旋翼/机身组合模型前飞算例进行了数值模拟,计算得到的机身表面压力系数分布与实验结果吻合良好,说明了该方法对旋翼/机身气动干扰研究的有效性。之后进一步将该方法应用到悬停状态的X3构型复合式高速直升机旋翼/机翼/螺旋桨组合模型的复杂流场模拟中,并与悬停状态孤立旋翼和旋翼/机翼组合模型的流场进行了对比。研究发现:机翼对旋翼下洗流起阻滞作用,引起机翼下方不规则流动,且螺旋桨滑流与旋翼下洗流会相互干扰产生一定的偏折,旋翼下洗流速度更大,螺旋桨滑流会产生明显的向下偏折。      

高超声速球锥粘性激波层压力分裂法数值求解

本文将压力分裂法引入细长钝体高超声速粘性激波层数值求解中,较好地解决了这种流场远后身部分的推进求解问题。针对非解析球锥体绕流数值求解的困难,作者提出了实用型初始激波形状计算方法,达到了简化数值求解步骤,提高数值计算效率的目的,数值求解结果与实验结果相符合。     

高超声速层流气动热预测混合算法研究

提出了兼顾激波和边界层模拟的混合算法--在激波区逼近NND格式,在流动光滑区逼近三阶偏置格式,采用了Steger-Warming、KFVS和Vanleer三种典型的矢通量分裂方法,对来流马赫数为10、基于头部半径的雷诺数为1×105和1×106的球头绕流的壁面热流进行了计算.计算结果表明:混合算法明显降低了热流对网格的依赖性.在此基础上完成了钝双锥和双椭球壁面热流的预测,得到了与实验测量一致性比较好的结果.     

航天器表面污染物质沉积变化和控制因子评估

航天器进入空间环境以后,空间环境分子污染和颗粒污染形成了航天器表面污染层,从而对航天器的各技术分系统产生不同的负面影响.介绍了中外中轨道航天器表面污染物质沉积变化在轨探测结果,并对污染物质沉积量变化和控制因子做初步评估.结果表明,污染物质沉积量在航天器入轨初期的1~2年内受航天器自身出气物质量、放气速率、表面温度及所处的气流方向等因子所控制.初期沉积量大,正是受到航天器入轨后自身出气量大、放气速率较高的控制,同时迎风面比背风面沉积量大.入轨后期表面沉积量长期变化呈现出明显的降变或缓慢涨落,而且具有全向性特征,因此探讨了具有全向性影响能力的控制因子相关特性,其中高能粒子通量和太阳紫外辐射通量变化可能是主要控制因子.      

基于SAMPEX卫星观测的南大西洋异常区高能质子动态分布特征

利用SAMPEX卫星1992年7月至2004年6月19～27MeV高能质子数据对南大西洋异常区的分布特征进行研究, 发现南大西洋异常区高能质子分布随高度及F_10.7的变化十分显著. 在540±25km高度上, 地磁较为平静时期南大西洋异常区高能质子微分通量随着F_10.7的增大而减小, 同时在F_10.7≥115sfu时减小趋势较为平缓. 对中等及以上磁暴进行统计分析发现, 磁暴期间南大西洋异常区高能质子微分通量和SYM-H指数的绝对值存在明显的反相关关系, 且地磁暴对南大西洋异常区高能质子微分通量存在明显的持续影响效应. 磁暴发生期间高能质子微分通量明显减少. 磁暴恢复相及其之后高能质子微分通量呈现较为显著的恢复过程.      

基于粘接元技术的区域分裂解法及其应用

针对线性椭圆型问题的，发展了一种无重叠区域分裂方法，首先引入了Lagrange乘子，以使子域交界处解的连续性约束条件获得弱满足。接着重点讨论了Lagrange乘子的近似空间，即粘接元（mortar elements)空间的建立，以及所引起的离散鞍点问题的共轭梯度迭代解法。这种区域分裂方法非常适用于网格不匹配情形。然后，对基于二阶鼓包（bump)函数近似的当地事后误差估算方法提出改进，以使之适应于区分裂情形，我们还将上述方法应用于广义Stokes问题。最后，给出了以误差估算值为准确则的自适应网格上的数值结果。