大气中激光等离子体通道的传播特性分析

对大气中激光等离子体通道的传播终端特性进行了研究。建立了大气中激光等离子体通道的近似电磁模型,基于圆柱坐标系中纵向分量所满足的波动方程,给出了纵向场与横向场的关系,由边界条件导出了大气中激光等离子体通道的特征方程。特征方程的数值解获得了相应色散曲线。     

探针式光纤探头的有效检测深度

研究了探针式光纤探头的有效检测深度。采用扩散方程描述光在强散射生物组织中的传播规律,运用纽曼边界条件,使用灵敏度矩阵对有效检测深度进行了研究。灵敏度矩阵和光源与检测器的空间位置分布密切相关。光源与检测器对的灵敏度矩阵定义为：光源位置和检测器位置分别作为光源的光场分布的乘积。设计了6组针对不同参数的模型,参数主要包括光源、检测器光纤芯径、光源、检测器光纤中心距、光纤探头所处的深度、生物组织的吸收系数和散射系数。对探针式光纤探头在各组模型中的有效探测深度进行了仿真,重点研究了六种参数对检测深度的影响。所得结论可以推广到分析光纤检测器的空间分辨率和多光纤检测器的研究中。     

基于带波动算子非线性Schr(o)dinger方程数值分析的守恒差分算法

对一类带波动算子的非线性Schr(o)dinger方程进行了数值分析,提出了一个含参数的二阶守恒差分格式,根据参数选取的差异,该格式既可隐式计算也可显式计算.对初值条件进行了中心差分离散,使其具有二阶精度,从而与守恒格式的精度一致.利用矩阵理论证明了差分解的存在惟一性,并利用一个重要的不等式在先验估计的基础上,运用能量估计的方法证明了该格式按无穷范数以二阶精度收敛到真实解.数值实验表明该格式具有较高的计算效率.     

离散速度方法及其在跨流域问题中的应用研究进展

往返大气层跨流域飞行器的气动力和气动热计算一直是当前流体力学研究的难点和热点之一。由于该类飞行器面临的流动场景不再是单一的连续流或稀薄流,采用Navier-Stokes方程求解器或DSMC方法均不能获得全流域的准确结果。近年来,以不依赖于连续性假设的Boltzmann模型方程为基础,通过在位置空间和速度空间同时离散求解该方程,实现了跨流域问题的统一求解。本文对该类算法的研究进展进行回顾和分析,着重介绍气体动理论统一算法（Gas kinetic unified algorithm, GKUA）、统一气体动理学格式（Unified gas kinetic scheme, UGKS）和改进离散速度方法（Improved discrete velocity method, IDVM）3种数值途径,分析它们的基本假设和实现方式,关注它们目前取得的进展和应用情况。同时,本文还将IDVM进一步扩展到非定常情形,以便用于非定常跨流域问题求解。最后,本文对该类算法存在的一些问题进行讨论,期望在后续的研究中能予以解决。     

亚声速大迎角模型试验洞壁干扰修正方法研究

通过测量洞壁附近的压力分布来模拟透气壁试验段的洞壁边界条件,采用数值求解Euler方程的方法模拟模型在风洞中的绕流场,然后将洞壁边界条件用远场边界条件替换,用同样方法计算出模型在自由流中的绕流场,从而计算出洞壁干扰对模型的绕流场和气动力的影响。针对大迎角模型的洞壁干扰问题的特殊性,采用重叠网格技术,对飞机大迎角标模和SB-03模型进行了亚声速大迎角的洞壁干扰计算与修正。     

用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法,前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率,而后者保证了网格间边界附近流场求解精度;本文非定常Euler方程的求解采用有限体积法,通量计算采用Osher格式,时间积分采用双时间步长推进。为验证本文提出的方法,分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率,保证了网格边界处的精度,具有较高的工程应用价值。     

航空发动机高速球轴承动态特性分析

 高速滚动轴承的动态特性对航空发动机转子系统的性能有重要的影响,为尽可能真实模拟轴承实际运转情况,建立轴承零件相互作用模型得到各零件之间力和力矩,使用较精确的弹流模型计算拖动力,采用拟动力学法建立任意受载的球轴承动力学分析模型,提出了工程中非线性方程组的高效解法。编制调试了计算程序READ ( Rolling Element bearing Advanced Dynamics tool),以276927NK1W1(H)航空轴承为算例分析了工况参数和结构参数对滚动轴承动态特性的影响规律,最后通过Gupta的动力学程序算例验证了程序的准确性。轴承动态特性的研究为研究航空发动机转子轴承集成系统动态性能分析奠定了基础。     

从闭合式格林函数推导微带表面波本征值方程

从闭合式格林函数推导出微带介质结构的表面波本征值方程,给出了一种新的表面波本征值方程的推导方法。在闭合式格林函数中,通过推导出平面介质结构的表面波本征值方程,可以方便的从中提取出表面波的贡献,从而使其用Prony’s方法或GPOF方法近似时,可使曲线更平滑。     

Influence of shock attenuation on tailored operation in free piston shock tunnels

The free piston shock tunnel is a type of shock tunnel with high performance. For this type of tunnel, the influence mechanism of shock wave attenuation on tailored operation is explored by numerical simulation and theoretical analysis. By introducing the normalized velocity, the simple constraint equation for shock wave under the tailored operation is deduced. Moreover, the real gas effect is also taken into account in this equation. Based on the equation, the tailored operation of shock tunnels can be predicted with very few calculations. The present study shows that the change rate of the thermodynamic state of the gas behind the shock wave is inconsistent with the attenuation rate of the shock wave, which is the fundamental reason why the wind tunnel achieves tailored operation at a lower Mach number of shock waves. This lower Mach number of shock waves differs from the corresponding ideal value by a factor, which is about the square root of shock attenuation rate.     

熵变量欧拉方程组跨超声速流SUPG有限元数值解

通过变量代换将拟线性守恒变量欧拉方程组转换为具有对称性的熵变量欧拉方程组，对该方程组构造SUPG的弱解式，其中权函数摄动项中包含流向算子和间断捕获算子。对翼型的跨超声速流动算例表明，本法具有较好的稳定性和能有效地抑制激波两侧振荡。应用带预处理的GMRES迭代法求解非对称方程组能显著地提高收敛速度。