三维球坐标系下波包非线性传播过程的一种数值求解方法

采用全隐欧拉(FICE)格式,将三维球坐标系下的大气基本运动方程作为控制方程,建立起了三维球坐标系下可压大气中的全非线性数值模式;对纬向边界也进行了合理的处理,并对该模式进行了数值检验.计算结果显示出重力波波包在大气中非线性传播的三维图像:波包整体向上传播,波相关扰动速度随高度增加而增长.     

固定剖分算法一致收敛的一个基于表征支集的判别方法

固定剖分（stationarysubdivision）方法是CAD／CAM、计算机图形学、图象处理等学科中一种新兴而有巨大发展潜力的几何造型方法。本文利用收缩性及样条函数给出了根据表征支集判别固定剖分方法一致收敛的一个充分条件。     

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

块结构的MUSCL算法求解三维进气道流场

本文将ＭＵＳＣＬ格式扩展为三维算法,并用于湍流流场的数值模拟。通过求解超声双斜楔压缩面流动,得到了含有激波、嵌入激波、剪切流和湍流边界层等流动特征的复杂流动结构,数值解与实验数据吻合得很好。文中将此算法与块结构网格技术相结合,使之能求解具有三维复杂几何形状的流场。并应用该算法模拟了Ｆ－１６战斗机进气道在５种典型飞行状态下的三维湍流流场,计算结果清楚地揭示了这些情况下的三维流动现象,得出了进气道的性     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

高阶精度方法下的湍流生成项对跨声速流动数值模拟的影响研究

利用五阶空间离散精度的WCNS格式和多块结构网格技术,通过求解雷诺平均NS方程,开展了SST两方程模型不同湍流生成项组合方式对跨声速流动数值模拟影响的计算分析。研究的主要目的是为高阶精度格式在复杂外形上的工程应用提供技术支撑。计算模型采用了RAE2822超临界翼型和DLR-F6翼身组合体构型。研究内容主要包括不同湍流生成项对残差收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布以及模型整体气动力特性的影响。不同湍流生成项组合方式的流场计算结果还与风洞试验数据进行了对比。研究结果表明:对于小迎角不存在明显分离的跨声速流动,不同湍流生成项对流场的高精度计算结果的影响很小,可以不用考虑。     

低耗散高可信度格式及其在气动力热计算中的应用

本文简要介绍了TVD类型格式中近似黎曼解的耗散机理,该耗散机制对于粘性问题的计算有较大的负面影响.通过引入低耗散的格式,所发展的数值方法对粘性流场的气动力热问题的计算具有较高的精度.以该数值方法为基础发展的软件系统,在亚跨超和高超声速领域中的气动力热问题上均得到了验证,计算结果与实验进行了对比验证,表明该计算软件系统具有较高的可信性.     

高超声速三维化学非平衡流动的数值模拟

高温非平衡化学反应对高超声速再入飞行器绕流流场的气动热特性有较大影响.通过采用有限体积NND(Non-oscillatory, containing No free parameters and Dissipative scheme)格式求解三维N-S方程,分别计算了小攻角条件下球锥的完全气体粘性流场和化学反应非平衡气体流场.比较并分析了化学非平衡作用对三维高超声速钝头体绕流流场的影响.所得结果与国外文献结果符合的较好,证明了计算方法与程序的有效性和实用性.      

一种基于方差的自适应火星图像阈值选取算法

在基于光学成像的深空探测自主导航中,图像边缘处理直接影响视线矢量的提取,从而对自主导航的精度产生较大影响。传统的Otsu算法更侧重同区域灰度的均匀性,适用于图像中目标区域和背景区域面积相差不大的情况,因此在自主导航初始阶段,导航精度较低。根据深空探测巡航段拍摄到的火星图像的特点,基于已有火星探测任务的实拍图像,在Otsu算法的基础上,重新设计准则函数,提出了一种基于方差的火星图像阈值自适应选取算法。该算法将灰度值高于待定阈值的区域的方差表示为待定阈值的函数,以函数一阶微分的最大值对应的灰度值作为最优阈值。该方法具有计算量小、图像分割精确的优点。仿真结果表明,相较于传统的Otsu算法,通过该算法得到的图像阈值能够实现更高的视线矢量提取精度及自主导航精度。     

动态地球磁层的时空剖分编码

动态地球磁层时空剖分模型借鉴了广泛应用于地学大数据的地球格网模型思想,实现了面向地球磁层的动态、非规则化物理空间的多层级剖分,且剖分格网形变稳定,可以实现一定范围内地球磁层时空特性的形式化。在此基础上,对剖分得到的时空格网进行编码表达,以便于计算机存储和处理,这是构建地球磁层大规模观测数据统一管理基础时空框架的另一个关键问题。由于地球磁层特殊的时空特性,经典、单一思路的地学编码方案较难完整地反映格网间的时空关系,因此难以支持格网间的基础时空关系计算。本文融合了整数坐标编码与Morton曲线编码的基本思路,实现了对漂移壳剖分格网的高效编码方案设计,完成了动态地球磁层的时空框架构建,从而为地磁数据的高效组织和处理奠定了基础。实验证明,本文提出的编码方案编码效率较高,且可以支持高效的相邻关系计算,为动态地球磁层多源、多层级、异构的大规模观测数据的组织与计算提供了一种解决方案。