非结构重叠网格显式装配算法

针对重叠网格中洞映射法占用过多物理内存的问题,发展了一种改进型洞映射法;基于相邻单元搜索法,发展了一种基于相邻阵面的贡献单元搜索法;通过将割补法与隐式切割技术相结合提出了一种非结构重叠网格显式装配算法。该算法首先生成一套包围物面的笛卡儿网格,其次存储所有与物面边界相交的笛卡儿网格信息,最后根据所存储笛卡儿网格与所需判断的网格单元的相对位置来判断其是否为洞内单元。在成功判断出所有洞内单元后,以当前洞边界为初始阵面推进,同时以各个网格单元的物面距离为判别标准对重叠区域进行优化,生成最终插值边界。所提算法优化了传统非结构重叠网格装配过程,具有物理内存占用低,贡献单元搜索次数少以及计算效率高等特点。通过2个典型复杂流动算例验证了所提算法的准确性与适用性。      

洞映射方法的研究和改进

针对传统洞映射方法存储大的缺点,对"广义封闭"的概念进行扩展,提出了最小洞映射方法,该方法允许挖洞曲面结束于网格截断面,有效缩小了洞映射区域,节省了存储开销.对适用于广义封闭的洞映射单元识别方法进行了分析,指出了用物面信息判断映射单元属性可靠性较差.发展了一种新的特别适合广义封闭问题的识别方法:"Inverse mark",使用计算网格结点信息自动识别洞外单元,再作为初始点在网格内部推进.研究表明:"Inverse mark"方法计算效率高,可靠性好,自动化程度高,有效提高了重叠网格方法对缝隙等局部复杂结构的适应性     

混合重叠网格插值方法的改进及应用

混合重叠网格间合理插值是保证流场计算正确的基础之一。本文针对重叠区网格尺寸匹配度较差时插值误差较大问题,发展了一种新型混合重叠网格插值方法,通过使用二阶精度插值、按单元尺度区分和扩充模板,改善了插值区网格尺寸匹配度较差时的插值精度。本文方法适用于任意单元类型的混合网格重叠,对各类单元处理透明,实现简单。计算结果表明,采用本文方法在网格重叠区流场变量传递正确有效,插值区网格尺寸匹配度较差时,相比原始方法,等值线过渡更为光滑,变量经过插值区耗散更小,计算与试验值符合更好。     

基于虚网格的格心ADT搜索法

针对重叠网格寻点问题,提出了一种直接在格心网格下操作,基于虚网格的ADT(Alternating Digital Tree)搜索方法.通过对计算网格进行扩展,建立虚边界网格,解决了格心网格因覆盖区域不完整而无法直接建立ADT数据结构的问题.搜索结果即为包含解的格心单元集合,可直接对结果列表遍历以得到合理贡献单元,故完全摒弃了可靠性差的Stencil Walk方法.由虚网格的定义,使寻点在边界附近的处理更为灵活,可以准确给出边界附近贡献单元的有效信息,同时简化了虚网格体系的构建.扩展了重叠边界类型,构建的搜索空间能完整覆盖网格范围,解决了对称面重叠问题.算例研究表明:该方法可靠性好,边界处理能力强,有效提高了重叠网格方法对大型复杂网格的解算能力.     

B样条应用于超限插值的网格生成方法

采用B样条曲线曲面生成技术与传统网格生成方法中超限插值技术,并将两者结合起,形成一种计算网格的生成方法.它发挥了两者的优点和长处:能充分反映复杂外形的完全信息,提高对外形的适应程度;可控制切矢的边界条件,提高生成网格的质量.算例表明了上述特点.      

关于非结构网格生成中的背景函数

背景信息对于用阵面推进法生成非结构网格是一个关键的环节.文章基于众多文献中提到的背景信息提供方法,提出一种更为实用的背景信息提供技术.此方法不仅避免了不必要的机时,内存的消耗,又使网格具有光滑性,网格的生成具有灵活性,方便性的特点.基于生成的网格采用有限体积法研究无粘高超音速化学非平衡流的流场分布,生成的网格能够成功的运用于高超音速化学非平衡流的数值计算。     

二维非结构网格的可压缩Euler方程流场解

讨论了非结构网格的生成方法,用Delaunay三角化方法和推进阵面法相结合的方法生成了非结构网格,提高了生成非结构网格的自动化程度和网格的质量.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼形和两段翼型的算例,并与实验结果进行了比较和分析.     

燃烧室三维贴体网格的生成方法研究

研究了燃烧室三维贴体网格人机交互生成的方法.使用该方法可以生成包括扩压器、内外环通道和火焰筒在内的燃烧室流场计算所需的整个网格.可以对网格的疏密和正交性进行调整.可以处理火焰筒上的主燃孔,掺混孔和气膜孔.并能够输出边界曲面上网格点的类型,用于流场计算.网格边界线的输入和边界线上的网格点分布采用人机交互方式生成.使用VISUAL BASIC语言编写程序界面.把整个网格生成及图形显示功能集成在一个软件中.     

二维有限元网格全自动生成

利用ＰＥＤＳ（ＰａｒａｍｅｔｒｉｃＥｎｇｉｎｅｅｒｉｎｇＤｅｓｉｇｎＳｙｓｅｍ）造型系统的几何，拓扑信息实现了二维有限元网格自动生成，在实现过程中，改进了Ｗａｓｔｏｎ三角化算法，把此方法应用到有限元剖分中，还利用了迭代的Ｌａｐｌａｃｉａｎ算法对网格进行了均匀化处理，另外，还提供了可靠，快速对网格编辑，局部细化的方法，最终使有限元网格生成完全自动化，而且速度快，质量好，密度易于控制。     

图象识别法在磁暴预报中的应用

本文将所研究的全部太阳和地球物理资料时间划分为442个时间段,其中前195段为“下降期”,后247段为“低年期”.按一定的选暴标准将其分为有磁暴时间段“D”和无磁暴时间段“N”.以表的形式提出与磁暴的发生有一定相关的12个物理先兆.利用计算机进行“D”和“N”的识别.结果表明:1)在“下降期”中,综合耀斑指数CFI≥6、持续时间大于10分钟的IV型射电爆发和行星际磁场的扇形边界过地球等现象,是引起磁暴发生的主要原因;2)在“低年期”,当太阳上的耀斑及其相关活动现象不发生时,产生磁暴的主要条件是反映27天重现性的太阳活动源(冕洞发出的高速流)和行星际磁场的扇形边界过地球等先兆;3)在“试验时段”内,“D”和“N”类识别的识别率为73—82%,在“预报时段”内用此图象识别法识别率为73—80%,误识率均不大于30%.用多组阈值(K)进行了正交设计试验,并经过“显著性”检验,本方法成功率的置信度大于97.5%.      

The space-developed dynamic vertical cutoff rigidity model and its applicability to aircraft radiation dose.

We have developed a dynamic geomagnetic vertical cutoff rigidity model that predicts the energetic charged particle transmission through the magnetosphere. Initially developed for space applications, we demonstrate the applicability of this library of cutoff rigidity models for computing aircraft radiation dose. The world grids of vertical cutoff rigidities were obtained by particle trajectory tracing in a magnetospheric model. This reference set of world grids of vertical cutoff rigidities calculated for satellite altitudes covers all magnetic activity levels from super quiet to extremely disturbed (i.e., Kp indices ranging from 0 to 9+) for every three hours in universal time. We utilize the McIlwain "L" parameter as the basis of the interpolation technique to reduce these initial satellite altitude vertical cutoff rigidities to cutoff rigidity values at aircraft altitudes.     

基于改进关键帧选择的ORB-SLAM2算法

针对同时定位与建图（SLAM）算法精度不高且跟踪易失败的问题，提出了一种改进关键帧选择的ORB-SLAM2算法。通过ORB-SLAM2算法计算帧间相对位姿；在原有算法的基础上，增加旋转与平移量作为判定依据，决定是否创建新关键帧；针对移动机器人所安装的相机与机器人产生相对运动引发误拍摄，导致劣质关键帧生成的问题，设计了劣质关键帧剔除算法；基于RGB-D数据集与自主研发的移动机器人进行了实验验证。实验结果表明:改进的关键帧选择算法能够准确及时地选择关键帧，最优情况下定位误差约为原误差的51.9%，有效消除了相机与机器人之间相对运动产生的影响，直线误差仅为原误差的82.1%。改进算法能够有效提高定位精度，减少跟踪失败。      

基于矢量场旋度的流带可视化研究

流带是计算流体力学(CFD,Computational Fluid Dynamics)矢量场可视化的重要技术之一.在对现有的几类流带生成算法进行比较分析的基础上,重点研究了基于矢量场旋度的流带生成算法,并且根据复合函数微分法的原理,针对结构化网格的定常三维矢量场,提出了算法具体实现中求解旋度的方法.流带可视化的结果有助于研究流场中涡的结构.      

重叠网格中隐式装配策略的改进

作为重叠网格技术的重要组成部分，隐式装配没有显式的“挖洞”过程，仅需在搜索贡献单元的同时通过对比不同单元的品质来分类单元。通过对传统隐式装配过程进行改进，提出了一种更高效的隐式装配策略，减少了贡献单元的搜索次数；同时，提出了一种基于笛卡儿网格映射的局部贡献单元搜索法，提高了网格装配效率。首先，针对每一个子网格，在计算其到自身物面的最短物面距离的同时也计算到其他所有物面的最短物面距离；然后，通过比较同一单元到不同物面的最短物面距离来控制洞边界的位置；最后，仅对插值单元进行局部的贡献单元搜索，避免了针对所有单元进行全局贡献单元搜索的过程。通过3个典型复杂流动算例验证了所提方法的准确性与高效性。      

Adjustment of systematic errors in ALS data through surface matching

Surface matching is a well researched topic in both Computer Vision (CV) and terrestrial laser scanning (TLS) or ground based light detection and ranging (LiDAR), but the extent of the range images derived from these technologies is typically orders of magnitude smaller than those derived from airborne laser scanning (ALS), also known as airborne LiDAR. Iterative closest point (ICP) and its variants have been successfully used to align and register multiple overlapping views of the range images for CV and TLS applications. However, many challenges are encountered in applying the ICP approach to ALS data sets. In this paper, we address these issues, explore the possibility of automating the algorithm, and present a technique to adjust systematic discrepancies in overlapping strips, using geometrical attributes in a given terrain. In this method, the ALS point samples used in the algorithm are selected depending on their ability to constrain the relative movement between the overlapping laser strips. The points from overlapping strips are matched through modified point to plane based on the ICP method.     

高超音速热流数值计算中的误差匹配原则

通过对粘性边界层的数值模拟和误差分析,提出了差分格式的误差匹配原则.在边界层内的计算网格一般都具有很大的长宽比,而计算误差与网格长宽比密切相关.为了使不同坐标方向上的误差保持在同一量级上,需根据误差匹配原则确定计算网格的尺寸.对绕钝锥的高超音速流动进行了数值模拟,进而计算了钝锥表面的热流分布.数值计算表明,参照误差匹配原则进行网格划分是有益的,可在保持流场数值计算精度和提高热流计算准确性的同时,使网格的数量和分布更为合理.