任意多面体的快速非结构网格生成算法

对于包含复杂气动外形的CFD数值模拟,网格生成是其中的关键,快速稳定的网格生成技术是其中的重要组成部分。本文首先建立了基于可视面的适用于任意多面体的快速初始化算法;然后改进了Delaunay生成算法后使用Delaunay改进生成算法细化网格,讨论了网格质量判定依据对网格生成的影响,通过开发Delaunay面交换技术优化网格生成过程;最后建立了基于顶点弹簧理论的网格优化方法,以提高网格生成的质量。结果表明,本文建立的算法效率较高。     

PVM平台下三维超声速钝头体底部流场的并行数值模拟及其拓扑分析

用拼接法划分网格，并运用基于消息传递的并行软件平台ＰＶＭ分配计算任务，实现了较大规模的三维超声速钝头体底部流场的并行数值模拟。结果表明，底部压强特性与实验符合较好，流场的拓扑结构合理。     

混合尺寸目标电磁计算的时域伪谱/有限体积混合方法

针对单一算法对混合尺寸目标进行时域电磁分析的困难,提出一种时域伪谱(PSTD)同时域有限体积(FVTD)混合方法。FVTD可方便地分析复杂的几何结构和材料,但是难以计算电大尺寸的目标,PSTD则特别适合计算电大尺寸的规则结构,但在模拟复杂的几何结构尤其是带有曲边结构以及电大、电小共存结构时存在困难。混合方法克服了单独算法的缺点,融合各自的优势,提高了算法的求解能力和应用范围。为了减小两种算法连接边界带来的反射,采用了FVTD计算面均值的二次函数重构方法,给出了交叠网格和非交叠网格两种混合方案。数值试验表明,混合方法有较高的精度,具有时域分析混合尺寸目标电磁问题的能力。     

数值模拟金属靶抗侵彻中的网格密度影响研究

在数值模拟中,有限元网格的数量对计算结果和计算成本都有较大的影响。采用MSC.Patran建立了弹体侵彻金属靶板的计算模型,并利用MSC.Dytran对所建模型进行了计算;结合相关的理论,对不同网格密度划分的有限元模型的计算结果进行了分析,并和实验进行了对比,得到了数值模拟弹体侵彻金属装甲中靶板的最佳网格数量。当网格无因次量μ=0.10时,所得到的数值模拟结果与实际吻合得较好:当厚度方向的网格数取为12个时,文中所采用的模型能较好地模拟弹体对靶板的侵彻效果。     

弹塑性流体计算中的一种实现方法

针对弹塑性流体计算问题，在ALE方法基础上提出了一种数值模拟的实现途径，主要包括含弹塑性的流体力学数学模型的建立，滑移面处理的方法，网格生成等，最后给出计算结果。     

基于新型非结构有限体积解耦算法的激波诱导燃烧数值模拟

把基于结构网格有限差分方法建立的化学非平衡流动模拟的新型解耦算法推广到基于非结构网格的有限体积法,从而可用于模拟复杂几何构型下的化学反应流动.对H2/ Air预混气体中激波诱导振荡燃烧的Lehr实验进行了数值模拟,计算得到的振荡频率与实验结果符合很好,表明计算方法具有时间和空间二阶精度.通过对不同几何外形发射体的计算发...     

一种新的基于方向导数的二维自适应网格生成算法

结合形变原理及网格迭代思想,利用方向导数计算控制函数,提出一种新的二维空间自适应网格生成算法。数值实验表明,该算法能较好地适应解函数的空间剧烈变化。与其他自适应算法比较,其主要优点是该算法逻辑简单,避免了解网格偏微分方程,节约了网格计算时间。     

基于STAR-CCM+的全附体潜艇尾流场数值分析

为充分了解潜艇周围的粘性流场特征,基于STAR-CCM+仿真软件对SUBOFF全附体模型的尾流场进行了数值模拟。通过将数值模拟结果与试验结果进行比较,系统分析了棱柱层网格的划分与湍流模型的选取对计算结果的影响。结果表明:SSTk-ω模型较适合应用在潜艇模型尾流场计算中,并且棱柱层的总厚度不能小于15mm,第一层计算网格节点到表面的无因次距离y+应在25到170之间。在此基础上,应用滑移网格技术对潜艇尾流场进行数值模拟,通过加密网格和减小时间步长使得轴向伴流分数与整体定常计算值吻合较好。     

平面运动激波绕障碍物流动的非结构Euler解

用三角形非结构网格及有限体积法求解Euler方程 ,对运动激波跨越不同截面形状的二维障碍物绕流进行了数值模拟。计算结果给出了不同时刻各流场的复杂结构 ,该结构与实验结果一致。计算表明用非结构算法模拟复杂外形绕流具有简单、准确的特点 ,并能避免尖角处非物理结果的产生     

一种改进的网格信任模型

提出了针对网格计算环境改进的信任模型,给出信任计算模型的信任值更新函数,根据网格实体的交易行为特征适当增加或减少其信任值。该模型对域内信任关系和域间信任关系分别采取不同的方法进行处理,能够评估实体与实体间的信任关系;在信任的计算中引入风险机制,防止协同作弊和诋毁的安全隐患。改进后的模型精化了信任算法的粒度,提高了信任评价的准确度。该可信度模型具有较好的可扩展性,可较好地解决网格环境中存在的某些安全问题。     

复杂平面区域的三角网格生成算法

为生成复杂平面区域的有限元网格,提出了基于网格细化的三角网格生成算法.该算法首先采用耳尖移除法对区域边界做三角划分,得到粗略的初始网格.提出Delaunay优化平分方法,根据网格密度细化初始网格,该网格细化方法结合最长边平分技术与Delaunay边交换技术,可有效提高内点生成与单元细分的质量.实验表明,基于Delaunay 优化平分的三角网格生成算法可对任意平面域进行网格剖分,生成符合有限元计算要求的高质量三角网格.     

树网结构上一种新的矩阵迭代求逆并行算法

运用树网结构可以完成矩阵的并行快速求逆,其中迭代法是一种非常重要的方法。本文给出了一种新的迭代格式,对任意非奇异矩阵Ａ,运用新的迭代格式对Ａ求逆相对于经典牛顿迭代法,在同样精度要求下,时间可减少一半。     

含气膜孔涡轮叶片寿命可靠性设计优化中的网格参数化方法

为了自动实现含气膜孔涡轮叶片寿命的可靠性设计优化,减少其有限元计算时间成本,提出了含气膜孔涡轮叶片的非结构局部网格参数化方法.该网格参数化方法将叶片划分为参数化区域和非参数化区域,分别设计非参数化区域的网格以及参数化区域的网格控制点,更新气膜孔的形状和位置后,保持非参数化区域的网格不变,重构参数化区域的网格.相比于结构...     

中尺度数值天气预报模式MM5分布式并行计算

中尺度数值预报模式是进行中尺度天气预报的有效手段。中尺度模式MM5是国际上应用最广泛的中尺度预报模式之一。数值天气预报的巨大计算量和实效性要求必须通过高性能分布式并行计算来实现。分析了MM5串行算法的特点,研究了其并行算法的实现,讨论了算法的一些改进,给出了MM5模式在分布式并行巨型计算机上的测试结果。     

非结构动网格分区并行计算方法

针对包含运动边界的非定常流动问题,基于弹簧近似法设计了一套非结构动网格分区并行计算流程,引入了METIS软件进行非结构网格二次剖分,研究了对接边界通信类型及其处理方法,基于消息传递并行编程模式实现了非结构动网格求解器的并行化。测试结果表明:该求解器具有较高的并行效率,网格变形模块对非结构动网格并行计算性能具有重要影响。     

三角域上带多个形状参数的二次Bézier曲面片

对于三角域上二次Bézier曲面的形状调整问题,给出了带6个形状参数的拟二次Bernstein基函数。分析了该组基函数的性质,并由此组基函数定义了三角域上带形状参数的二次Bézier曲面片,该曲面片不仅具有三角域上二次Bézier曲面片的性质,而且在控制网格保持不变的条件下,具有形状可调性。形状参数具有明确的几何意义。实例表明：构造的曲面为几何造型中的曲面设计提供了有效的方法。     

任意多边形的非结构网格快速生成算法研究

对于包含大尺度运动边界的CFD数值模拟,网格重构是其中的关键,快速稳定的网格生成技术是其中的重要组成部分。建立了基于有向边的适用于任意多边形的快速三角初始化算法;证明了最长边剖分网格细化算法在一定条件下发散,并结合Delaunay边交换技术使细化算法封闭;建立了基于顶点弹簧理论的网格优化方法,以提高网格生成的质量。结果表明,算法具有较好的鲁棒性和高效性。     

基于蛇怪蜥蜴踩水机理的仿生推进装置数值计算方法研究

为突破传统排水型两栖车辆“阻力墙”现象,提高两栖车辆水上航速,从蛇怪蜥蜴高速踏水机理出发,设计一种新型轮一叶复合式水上仿生推进装置,采用滑移网格和动网格方法对不同转速下的推进装置三维模型的动力学特性进行数值计算,并对仿真结果进行对比和分析．结果表明：滑移网格和动网格方法在计算结果上获得比较一致的结果,动网格方法更能反映模型的真实运动情况,而滑移网格方法在计算过程中表现出较好的稳定性,同时显著减少计算时间