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

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

Chebyshev谱方法在声散射数值计算中的应用

基于声学边界元基本理论,并结合Chebyshev谱方法进行了声散射数值计算研究.采用2阶边界单元,谱点离散选用不包含边界的CG(Chebyshev-Gauss)配点法,克服了传统边界元在表面单元节点处出现的法向及法向导数不连续的现象;应用CHIEF(Combined Helmholtz Integral Equation Fomulation)方法进行了数值非唯一性处理,提高了计算的精度.与传统边界元方法进行了精度和效率的对比分析,证明该方法具有计算快速、精度高的特点.     

模拟弹塑性碰撞的修正SPH方法

从基本的无网格光滑粒子法SPH(Smoothed Particle Hydrodynamics)近似出发,修正了模拟固体力学中大变形弹塑性碰撞的SPH方法.在边界处采用修正的边界条件,弹塑性分析过程中采用增量理论计算应力,迭代过程中用守恒光滑法进行滤波修正消除拉力不稳定.对SPH方法进行了程序实现,给出了杆弹塑性碰撞的算例.计算分析表明,SPH方法节点的影响域较大、精度较相同节点间距有限元法的结果有一定差距,但是通过增加粒子数量可以提高SPH的精度,保持了其简单性和计算大变形的特性.      

NET-RAM的设计与实现

对SNOW机群系统的分布式共享存储管理机制进行适当的简化,就可实现网络虚拟存储环境.它以本地节点的主存作为cache,以网络上其它节点机的主存作为本地节点机主存的后援存储,取代本地节点机的硬盘或其它介质的存储设备.充分利用网络上节点的物理存储空间,使一些在单机上由于内存限制而无法完成的计算成为可能.这对于时间敏感性不高但需使用大内存的计算程序尤其适合.本文在介绍了网络虚存的概念之后,详细阐述了一个网络虚拟环境NET-RAM总体结构,软cache的实现方法,并对该网络虚存环境进行了测试,分析了测试结果.最后对今后的工作进行了探讨.      

基于DSI插值的三角网格质量优化

通过对三角网格的单元顶点进行几何位置调整,提高了网格的质量,实现了网格的质量优化.几何位置调整是使用离散点光滑插值(DSI,Discrete Smooth Interpolation)实现的,针对在计算时影响质量优化的邻接边界的单元顶点,采用了在边界处补偿三角形的方法,消除了单元收缩,提高了网格的质量.与加权拉普拉斯算法进行了比较和分析,优于拉普拉斯算法;为了使三角网格在位置调整时保持原始网格的几何细节特征,在插值算法中施加了控制点约束.最后使用算例对算法进行了验证.      

粘性流体换热问题用直角坐标网格的贴体解法

采用直角坐标网格但通过在边界附近保留不规则控制体的非结构化直角坐标网格方法求解粘性不可压流体流动和换热问题.在边界附近切割单元的求解采用常规的方法处理,而不采用插值方法处理,并在同位网格上对控制方程组进行离散.因而该算法是贴体的,并可适用于解决复杂边界流动换热问题.通过对同轴环行空间和二维倾斜空腔内的对流换热问题的计算对该数值方法在传热问题中的应用进行验证.对比结果表明,计算结果与精确解和STAR-CD的结果在一定Re数和网格数时是很接近的,可以满足一定的精度要求,说明该数值计算方法可用于传热问题.      

DSMC方法的边界采样与统计技术

研究了适合在三维直接模拟蒙特卡罗(DSMC,Direct Simulation Monte Carlo)方法计算中根据入口边界条件合理地布置进入计算区域内粒子参数的边界采样技术.提出了一种根据穿越入口边界、出口边界以及撞击壁面边界的粒子统计数据获取边界位置处宏观流场参数的新型边界统计技术,弥补了DSMC计算中传统统计技术仅能够获取流场内部参数,不能够获取界面处流场参数的不足.这种界面统计技术具有普遍的适用性,不仅限于DSMC方法.以无穷远来流圆球绕流为例进行计算,验证了DSMC方法的边界采样与统计技术的有效性.     

等效原理算法在不同等效面下的误差分析

等效原理算法(EPA)借鉴了区域分解的思想,并基于惠更斯等效原理实现了对大规模散射问题的计算.本文对EPA中的等效原理算子(EPO)进行了误差分析,选取了立方体、球体和进行了光滑处理的立方体等作为等效面寻找误差来源,并提出了误差改进的建议.根据计算结果,等效磁流在等效面上描述的误差一般集中在等效面表面的几何不连续处,且主要受等效面上的法向向量不连续性和基函数的选取的影响.故在选取等效面时应尽量考虑使用平滑的等效面,而不是直接使用立方体等效面.      

基于神经网络的撞球机器人控制器设计

对撞球机器人的母球控制问题展开研究,设计了一种基于神经网络(NN)的控制器,使机器人能够控制母球在击打目标球后按照预定的模式运动至目标点——即完成走位。针对该问题非线性且非光滑的特点,对坐标系进行阐述并给出机器人击球的模型;在光滑的假设下使用理论分析的方法建立母球的运动学模型与边库反弹的理想镜像模型;进而使用神经网络方法对理想模型进行修正,并对不同的轨迹模式进行分析与分类。测试结果表明:经过训练的机器人能够掌握各种模式的走位,统计结果与模型分析结果相吻合;相比于单一使用神经网络方法,本文使用理论分析与神经网络相结合的方法能够有效地提升网络的品质,降低训练的误差。     

多维三次样条函数及其应用

本文将一维三次样条函数推广到多维问题，它保留了三对角矩阵方程便于求解的特性，并能得到满意的插值、一阶和二阶偏导数。关于边界条件问题，若能给出边界节点处的一阶偏导数边界条件，则精度最高，但在大多数应用问题中，往往给不出此条件。为此，本文提出一种改进方法，即用拉格朗日三点插值法由域边附近的节点数据计算出边界节点处的一阶偏导数。算例表明，此法可改善精度。     

激波与物质界面相互作用的数值模拟

计算激波与物质界面相互作用问题的关键是精确捕获运动物质界面.采用level set捕获物质界面,用改进虚拟流(modified ghost fluid method)方法定义虚拟流节点参数,从而将多介质流问题分解成多个单一介质流问题求解.应用该方法对激波与物质界面相互作用系列问题进行了数值模拟,并给出了物质界面的演化过程.计算结果发现激波强度、物质界面对涡的形式有很大影响,界面两边的密度比决定了界面是否反转.     

第二代低阶面元法及在动力短舱上的应用

在气动力设计和计算中,面元法是主要实用计算方法之一。新近发展的第二代低阶面元法实用性强和计算精度较高。本文介绍了该面元法的理论,通过对比早期面元法,着重探讨了第二代低阶面元法的特点及其优越性。运用自行研制的计算程序,对低速动力短舱进行了计算和分析研究,结果表明与国外资料所提供的计算结果吻合得较好。     

瞬态传热问题的微分求积和精细积分求解方法

给出了瞬态传热问题的高效高精度求解方法,该方法分别用微分求积法(DQM)和精细积分法(PIM)离散空间域和时间域.微分求积方法除了精度高、效率高之外,处理复杂边界条件的灵活性也优于有限元法(FEM).用精细积分法处理一阶瞬态传热微分控制方程,不需要增加额外自由度,还可以达到计算机精度.给出了隔热结构4种边界条件下的数值结果.并就上表面恒温、其他面绝热边界条件计算结果与有限元分析结果进行了对比,算例分析表明,采用微分求积和精细积分法布置少量的节点就可以达到较高的精度.      

含局部缝隙与重叠特征的Mindlin板等几何分析

等几何分析(IGA)中非均匀有理B样条(NURBS)被同时用作计算机辅助设计(CAD)中的建模工具以及有限元分析(FEA)中的逼近函数。NURBS模型中常见的缝隙和重叠问题使得分析变得困难。基于Mindlin板理论,对含有缝隙与重叠部分的NURBS模型进行等几何分析,采用Nitsche方法处理模型交界面上的非协调问题,并通过标准数值仿真算例的计算结果与解析解对比验证方法的可行性。研究结果表明:基于Nitsche的等几何方法可以用来对含局部缝隙与重叠特征的非协调Mindlin板模型进行分析;NURBS次数越高,等几何分析计算结果越精确,并且收敛速度越快。     

三维发动机进气道水滴撞击特性分析

用带数值扩散项的欧拉法对某型三维发动机进气道水滴撞击特性进行了研究.实现了可压流场中三维复杂形体的局部水滴收集系数的计算.数值扩散项的引入解决了欧拉法计算产生奇异结果发散的问题.通过对不同进气道流量的算例结果的比较,发现进气道流量对进气道唇口附近的水滴撞击区和撞击量都有影响.进气量大,最大局部水滴收集系数也大.通过对有无整流罩的计算结果进行比较,发现此型发动机整流罩对水滴撞击特性的影响很小.      

椭圆轨道卫星编队燃料最优机动实时算法

文章针对椭圆轨道编队卫星最优燃料机动的实时计算问题进行研究。利用LGL伪谱法对问题进行快速求解。由于最优燃料控制为非连续控制,如直接采用伪谱法则只有当离散点数目很大时,才能实现较高精度的编队机动控制,而这会影响算法的实时计算效率。为了克服这方面的不足,文章提出了一种时间尺度变换策略,并结合模型预测的滚动优化方法,设计了一种编队实时机动控制算法。该算法在兼顾解的全局最优性的同时改善了解的局部精确性,最后通过数值仿真验证了所设计算法的有效性。     

基于OWRR-PF的机载图像中导弹跟踪方法

针对利用图像信息实现飞行器的导弹预警跟踪中导弹运动2D建模、运动模型非线性和目标干扰非高斯等问题,研究了比例导引下导弹在三维空间中的运动在成像面上的运动模式,建立了比例制导导弹的2D投影运动状态模型.由于模型中速度和加速度等主要物理量的非线性,并考虑导弹运动过程中受到的风向、风力、气旋和气流等随机干扰的非高斯性,采用粒子滤波方法实现导弹跟踪;并针对粒子滤波在跟踪过程中存在的粒子退化问题,提出有序权值残差重采样粒子滤波(OWRR-PF,Ordered Weight Residual Resampling Particle Filter)方法,该方法缓解了粒子退化问题,提高了跟踪的准确度.利用所建立的导弹运动模型进行连续视频试验,与标准粒子滤波相比跟踪精度提高了70%左右,与残差重采样粒子滤波方法相比跟踪精度提高了15%左右.     

某光电检测仪典型问题分析及解决

某光电检测仪器用于实时精确检测设备的动态安装方位准确度。就该仪器使用过程中出现的光能量分布不均、二次像点现象严重及测角零位漂移等几个典型问题进行分析,提出并采取了相应的解决措施,为以后解决类似的问题提供了技术途径。