非线性双曲型守恒律的一类应用通量限制器的二步二阶格式

本文提出了计算非线性双曲型守恒律弱解的一类新的二步二阶精度有限差分格式,这类高精度非线性差分格式的构造应用了三点单调差分格式的数值通量和通量限制器。文章证明,在拟CFL条件限制下,格式为总变差衰减差分格式(TVD格式),因此不发生常用的诸如二步二阶Lax-Wendroff格式、Maccomark格式和其他古典二阶格式计算不连续解时在其周围出现的剧烈振荡现象。 为了减少Roe推广技术引起的计算复杂性,文中我们基于Davis,S.F.的思想,使用了Rusanov数值通量,把标量差分格式推广到方程组,构造了求解双曲型守恒律方程组的新的预测-修正格式,较之其他一般的预测-修正格式而言,仅增加了很少的计算复杂性,初步的数值试验表明,本文提出的格式无振荡和具有极好的激波分辨率。     

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。     

基于体积分数形式的多介质流动数值模拟

把非守恒形式的体积分数方程与Euler方程组耦合，用交错中心型格式求解流体力学方程组。对一维情况下两种气体的相互作用和二维气泡的变形及塌陷问题进行了数值模拟，并把模拟结果与γ-model方法和NND（Non—oscillatory and non—free—parameter dissipation difference scheme，无振荡无自由参数耗散）格式的模拟结果进行比较，达到了抑制界面两侧非物理振荡和保持体积分数正性的目的。     

Euler方程的动力学通矢量分裂格式

文中基于气体分子动力学理论，研究和发展了一类新型逆风差分格式——动力学通矢量分裂格式（ＫＦＶＳ），并论证了格式的逆风性质和Ｌ１稳定性。最后给出了高精度格式的构造，并就一维激波管问题作了数值试验，还与精确解作了比较，结果令人满意。     

导弹垂直发射系统三维超声速射流的数值模拟

对导弹垂直发射系统(VLS,Vertical launching system)内的三维湍流流场进行了数值模拟.采用了单元中心有限体积法和三阶迎风TVD差分格式离散Navier-Stokes方程,采用双时间步和预处理过程提高计算收敛速度,引入对角化的双时间ADI算法求解控制方程.仿真结果准确地模拟了发射系统内的激波形成过程,真实地反映出VLS中的复杂流动.研究表明,压力室底部射流中心的波后压强以大约1 200Hz的频率振荡.      

格子玻尔兹曼和气体动理学通量算法及其应用进展

采用守恒律方程求解流体流动问题时,单元界面通量的计算尤为关键,该过程也被称为通量重构。由于离散控制方程的物理量定义在解点上,如何利用解点上的值来计算单元界面的通量,是计算流体力学最为关心的问题之一。针对该问题的研究已发展了各式各样的计算格式,例如完全基于数学重构的差分近似、基于部分物理重构的黎曼通量求解器以及近年发展起来的基于完全物理重构的气体动理学格式、格子Boltzmann通量算法和气体动理学通量算法。本文首先对几种典型的通量重构算法进展进行回顾和分析;然后着重介绍格子玻尔兹曼通量算法和气体动理学通量算法的研究进展及其相关应用;最后就该类算法存在的挑战和可能的研究方向进行展望。     

SAR极坐标格式成像算法对运动目标响应特性

运动目标的合成孔径雷达（SAR）成像特征是SAR地面运动目标指示（SAR/GMTI）系统进行运动目标检测、成像和重定位的基础。为此,从信号二维解耦合校正距离徙动角度推导和分析了极坐标格式算法（PFA）对运动目标的响应特性,揭示了极坐标格式转换过程中的距离插值和方位插值对运动目标距离徙动的校正原理。理论分析表明：PFA在完成对静止目标成像的同时,还能够自动校正所有运动目标的线性距离徒动,且校正过程无需目标运动信息,因此对于径向运动目标,在忽略波前弯曲条件下PFA也能够对其进行完全聚焦。最后,通过仿真数据处理对理论分析结果进行了验证,表明PFA作为一种运动目标成像预处理方法具有很大的应用前景。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格 库塔方法。采用四阶龙格 库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。     

预处理AUSM+格式在全速域反应流模拟中的应用

对预处理修正的AUSM+(advection upwind splitting method+)格式进行了详细讨论, 并将其应用于全速域范围内的化学反应流场数值模拟.基于多组分Navier-Stokes方程, 结合控制方程时间导数项的预处理, 采用预处理修正的AUSM+格式发展了一种模拟全速域范围内化学反应流场的有效算法.验证算例包括超声速、亚声速以及极低速情形的化学反应流场, 计算结果表明基于预处理AUSM+格式发展的算法能够有效地用于全速域范围内化学反应流动问题的模拟.