利用谱方法数值求解三维球坐标系下的Navier-Stokes方程

从三维球坐标系下的Navier-Stokes方程出发,采用谱配置方法建立了一个三维球坐标系下中层大气波动传播的全非线性动力学数值模式;构造了合理的纬向边界,并利用初始给定的重力波波包对该模式进行了数值检验．计算结果展示了重力波波包在大气中传播的三维图像,波包整体向上传播,波相关扰动速度随高度的增加而增长．     

三维球坐标系下重力波波包非线性传播过程听数值研究

采用全隐欧拉（FICE）格式，对三维球坐标系下重力波波包在中层大气的非线性传播过程进行了数值模拟，模拟结果完整地反映了重力波波包在中层大气的非线性传播过程，并与线性重力波理论进行了定量比较。     

三维球坐标系下重力波波包非线性传播过程的数值研究

采用全隐欧拉(FICE)格式，对三维球坐标系下重力波波包在中层大气的非线性传播过程进行了数值模拟.模拟结果完整地反映了重力波波包在中层大气的非线性传播过程，并与线性重力波理论进行了定量比较.     

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

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

大气波动传播问题的一种无反射数值边界格式

根据无反射边界条件的基本原则,考虑到高阶的数值边界格式可以减小边界误差和虚假反射,利用曲线拟合中最小二乘法的思想,提出了一种高阶的光滑拟合外推边界格式(SFEBS).因为大气运动的控制方程简化后可以归结为一个对流方程,所以作为边界格式之间的比较和检验,在一维情况下用对流方程和数值模拟中具有代表性的波包和激波作为算例,将其与传统的基于Taylor展开思想构造的边界条件(TEBS)进行了比较.计算结果表明,在高阶情况下,用SFEBS计算波动传播问题的虚假反射约为用同阶TEBS计算的1/6,说明高阶的SFEBS可以作为一种很好的无返射数值边界格式.为今后模拟大气波动提供了一种很好的数值边界格式.      

民用通风机气动设计中弯掠长叶片的数值模拟

采用数值模拟的方法,分别对弯掠叶片造型和直叶片造型两种情况下高性能民用通风机叶片内三维紊流流场进行了计算;计算程序中采用有限体积法,在任意非正交曲线坐标系下,用SIMPLE算法求解三维N-S方程和K-ε湍流方程.计算结果表明,弯掠叶片造型的民用通风机的流动参数和总性能参数更接近设计指标,而且采用弯掠叶片可以更有效地控制气流的流动,对减少二次流损失起了很大作用.      

三维重力波非线性传播数值模型的一种求解格式

采用交替方向隐格式（ADI）求解三维可压大气中重力波非线性传播的数值模型，经过几个波周期后，重力波波包括估然沿着线性理论射线方 稳定地传播，波包保持着很好的形状，内点和边界保持稳定；水平风与垂直运动的振幅随着高度的增加而增加。这说明ADI格式是模拟三维可压大气中重力波传播的一种可靠的格式，模拟结果表明，在非线性情况下，重力波仍可保持性传播的某些特征，非线性效应有现在对波包括传播速度的修正。     

基于HITRAN数据库的大气激光雷达信号仿真

信号仿真研究是大气激光雷达研究的重要环节，结果可为激光雷达系统的设计和研制提供基础.本文仿真系统可用于激光雷达回波信号模拟.该仿真系统利用HITRAN数据库中的吸收谱线计算大气分子对激光的吸收光谱，并采用大气辐射传输模型中的气溶胶模式模拟气溶胶对激光的衰减.利用激光雷达方程，数值模拟了355nm，532nm以及1064nm的回波信号，并将532nm的数值模拟结果与实测结果进行对比，评估了激光雷达系统的光学效率.      

模拟三维重力波传播过程的并行数值模式

从三维大气运动的基本控制方程出发, 基于MPI消息传递接口和区域分解的思想, 建立了模拟三维可压缩大气中重力波传播过程的并行数值模式. 在对垂直方向进行区域分解的基础上, 针对跳点网格的特点, 以垂直速度主格点线为子区域的下边界, 状态变量主格点线为子区域的上边界进行区域分解. 利用MPI消息传递接口来传递计算各子区域中变量时所要用到的相邻子区域上相应变量, 从而顺利地完成对整个计算域的并行计算. 根据线性重力波理论, 通过模拟小振幅重力波的传播过程对所建模式进行了验证. 结果表明, 模式可以很好地模拟小振幅重力波在三维空间中的传播过程, 模拟的重力波振幅随着传播高度的增加以指数形式增长, 对能量传播路径以及有效扰动位能和扰动动能的模拟结果均与理论预测的结果吻合很好. 另外, 随着进程个数的增加, 完成相同的计算量所用的计算时间也显著减少. 这表明本文建立的并行数值模式不但能够很好的模拟重力波的传播过程, 而且能够有效地节约计算时间.      

修正的动力盘模型与三维模拟螺旋桨滑流比较

比较用于研究螺旋桨滑流对机翼影响的4种数值模拟方法,即未进行修正的动力盘数值模拟、修正的片条理论应用于动力盘模型的数值模拟、有旋转的修正的片条理论应用于动力盘模型的数值模拟及定常三维实体螺旋桨数值模拟求解RANS(Reynolds Averaged Navier-Stokes)方程.以某半径为1.008 m三叶螺旋桨在转速为2 575 r/min,飞行速度为 66.889 m/s 工况下螺旋桨后的轴向和环向诱导系数作为评估依据进行比较,结果表明:修正的动力盘模型和三维实体模型数值模拟能预测螺旋滑流区内的流动情况,而三维实体模型能更多地反映流动细节,但是会增加计算成本,表明动力盘模型可以替代三维实体螺旋桨进行数值模拟.     

下行重力波波包在可压大气中的传播

应用ADI格式模拟下行重力波在二维可压大气中非线性传播。结果显示，重力波能量向下传输，波振幅明显衰减，与线性方程的数值模拟结果的比较表明，在非线性情况下，重力波仍可保持线性传播的某些特征，蜚 线性效应随波振幅的衰减对波包传播速度的修正作用很快变弱。     

重力波波包在三维非等温大气和水平风场中非线性传播的数值研究

采用全隐欧拉格式（FICE）对重力波波包在三维非等温大气、均匀和剪切风场中的非线性传播进行了数值模拟，给出了重力波波包三维非线性传播的全过程，分析了重力波的传播特性及背景温度场、风场对重力波传播的影响。结果表明：波包扰动速度振幅的增长比在WKB条件下振幅的增长要慢；波包非线性传播的路径、能量传输速度不同于WKB近似下的结果，非线性效应导致了重力波的传播特性的改变；温度场的非均匀性会改变重力波传播的路径和速度；剪切风场使扰动速度振幅的增加变得缓慢，使垂直波长减小。     

重力波波包在真实大气中传播特性的数值研究

采用二维全隐欧拉（FICE）格式，对重力波波包在真实大气中的非线性传播和演变过程进行了数值模拟，模拟结果表明，在中层大气下部激发的向上传播的重力波波包在传播到中层顶之前，波相关能量沿着射线路径传播，非线性效应和背景温场对波能量传播路径的影响很小，当波包传播到低热层大气后，波包饱和波相关能量几乎完全沿着水平方向传播，垂直方向的能量传播受到抑制，这与在无耗大气下，WKB近似条件下的线性重力波理论的预言相差很大，深入的分析表明抑制重力波波包向上传输能量的关键因素是大气分子粘性的垂直非均匀性。非线性和背景温场的影响不足以完全抑制波能量的向上传播，此外，在波包的整个传播过程中，由于非线性，背景温度和背景耗散的共同作用，重力波波包的垂直波长随时间明显减小，这些结果说明大气的分子粘性特别是分子粘性的垂直不均匀性对重力波波包在中、高层大气的非线性传播过程起着重要作用。     

重力波波包相互作用时空演变的数值研究

本文设计了一种研究二维空间中重力波波包共振相互作用演变的数值模式.利用此模式获得二个大振幅重力波波包通过碰撞而发生的完整的相互作用过程.数值结果表明,能量上行重力波可以通过波-波相互作用激发能量下行重力波;在相互作用的时空演变过程中,3个波包的能量收支情况和相互作用的强度具有局地性;由于考虑了空间传播效应,相互作用过程不再具有周期性,而是出现一个具有不同物理含义的相互作用特征时间,这个新的相互作用特征时间和相互作用最终的发展程度不仅由相互作用系数和初始波包振幅决定,还受到波包空间尺度和波包群速度的相对大小控制.      

重力波波包在非等温和可压大气中的非线性传播

本文采用二维全隐欧拉格式对具有高斯分布的重力波波包在非等温,可压大气中的非线性传播过程进行数值模拟。     

驾驶员－飞机闭环系统特性综合的一种方法

为获得驾驶员－飞机良好的匹配特性，根据Ｎｅａｌ－Ｓｍｉｔｈ闭环系统评价准则，建立了相应的边界条件方程．采用牛顿迭代和最速下降法相结合，迭代求得驾驶员模型参数和飞机模型参数．按良好飞行器品质要求，对闭环特性进行综合．驾驶员模型选用俯仰跟踪时的形式，飞机模型采用等效系统模型．     

Alfvén脉动的控制方程

本文由磁流体力学方程组导出了行星际空间中不可压缩小尺度脉动的控制方程组。对这一组方程的讨论表明,在小振幅极限下,脉动幅度的径向变化可由通常文献中引用的Alfvén波在缓变磁流体介质中传播的WKB解来描述。当脉动幅度与平均磁场强度之比为有限值时,在一般情况下,控制方程中的非线性项不能略去,因而不能用WKB解来描述有限振幅脉动幅度的径向变化。这一结论可以解释为什么在0.3—0.9AU实测Alfvén脉动振幅和谱的变化与WKB解不一致。在这组控制方程的基础上提出了一个定性的模式。将Alfvén脉动看作主要由许多具有不同波矢k的向外传播的有限振幅的Alfvén模式组成。它们在传播过程中将产生一级小量的波动。零级量与一级量的非线性相互作用使波能向高频区串级。零级量波动振幅的变化不仅受到太阳风慢变化的影响,而且受到这一非线性相互作用的影响。这一模式可以定性解释主要的Alfvén脉动的观测事实。