基于FDTD的时域混合方法及其在天线前门耦合数值模拟中的应用

强电磁脉冲环境下的平台-机载天线一体化耦合计算属于典型多尺度时域电磁计算问题,采用传统的FDTD方法数值模拟时,由于精细结构的存在导致网格量巨大,计算效率低下。介绍了一种将非均匀FDTD方法与细导线FDTD方法以及多网格集总元件FDTD方法相结合的时域混合方法,能够有效降低计算开销,结合并行计算技术,快速计算得到天线端口上耦合产生的瞬态电压和电流响应,并将该方法成功应用于无人机平台-天线一体化前门耦合数值模拟中。     

反应堆蒙特卡罗临界模拟中均匀裂变源算法的改进

在反应堆pin-by-pin精细建模及蒙特卡罗模拟计算研究中, 由于不同栅元的功率密度差异较大, 导致蒙特卡罗方法临界计算的样本在不同栅元之间的分配不均衡, 由此引起栅元内的各种计数的统计误差差异较大. 为使大部分栅元内计数的统计误差降至一个合理的水平, 单纯增加总样本已不是一个高效的解决方法. 通过在特定临界计算迭代算法的基础上改进并实现均匀裂变源算法的思想, 对大亚湾压水堆pin-by-pin模型取得了具有较高效率的数值结果. 本工作为具有自主知识产权的蒙特卡罗粒子输运模拟软件JMCT最终达到反应堆pin-by-pin模型(包括一系列国际基准模型)的模拟性能要求提供了一个有效的工具.     

二维球坐标系中子输运方程的一种并行SN算法

针对二维球坐标系下中子输运方程的SN算法, 提出基于(单元, 方向)二元组的有向图模型, 在已有的基于有向图的并行流水线算法基础上, 设计粒度可控多级并行SN算法。其中, 采用区域分解和并行流水线相结合的方式挖掘空间-角度方向的并行度, 提出能群流水并行方法, 并通过设置合适的流水线粒度来平衡有向图调度、通信和空闲等待开销。实验结果表明: 该算法可以有效地求解二维球坐标系下的中子输运方程。在某国产并行机1920核上, 对于96万网格、60个方向、24能群、数十亿自由度的典型中子输运问题, 获得了71%的并行效率。     

JMCT优化查找策略的实现

分析连续能量点截面中子输运程序JMCT中区间查找的时间复杂度.针对主能量网格,裂变释放中子数和不可分辨共振网格进行区间映射加速,并以等差和等比序列网格为例,对比分析不同散列网格的质量.最后选择寿期初与寿期中燃料棒栅元的k_∞问题进行对比测试.结果显示,对于核素数目较多的输运问题,优化查找策略具有很好的加速效果.     

动载下裂纹应力强度因子计算的改进型扩展有限元法

相较于常规扩展有限元法(extended finite element method,XFEM),改进型扩展有限元法(improved XFEM)解决了现有方法线性相关与总体刚度矩阵高度病态问题,在数量级上提升了总体方程的求解效率,克服了现有方法在动力学问题中的能量正确传递、动态应力强度因子数值震荡、精度低下问题.本文基于改进型XFEM,采用Newmark隐式时间积分算法,重点研究了动载荷作用下扩展裂纹尖端应力强度因子的求解方法,与静力学方法相比,增加了裂纹扩展速度项与惯性项的贡献.通过数值算例研究了网格单元尺寸、质量矩阵、时间步长、裂尖加强区域、惯性项、扩展速度项及相互作用积分区域J-domain的网格与单元尺寸对动态应力强度因子求解精度的影响,验证了改进型XFEM计算动态裂纹应力强度因子方法的有效性.针对文献中具有挑战性的"I型半无限长裂纹先稳定后扩展"问题,改进型XFEM给出目前为止精度最好的动态应力强度因子数值解.     

计算电磁学及其在复杂电磁环境数值模拟中的应用和发展趋势

概述计算电磁学的发展历程及其数值方法,并从学术研究、商业软件以及专用软件三个层次论述国内外研究的最新进展.结合复杂电磁环境的应用介绍复杂电磁环境数值模拟平台的研究现状,最后提出未来发展趋势的几点设想.     

强爆炸火球辐射流体自适应网格高精度数值模拟

提出一种自适应网格方法,应用于基于Euler方法的强爆炸辐射流体高精度数值求解。通过与Zinn数值结果对比,验证该方法的正确性。研究自适应网格对冲击波和光辐射输出模拟精度的影响,对比不同网格尺度下的计算耗时。在相同的条件下,使用自适应网格与均匀网格加密3倍得到的冲击波超压分布、光辐射输出演化接近,计算效率提升约8.5倍。由此可知,提出的自适应网格方法可用于强爆炸问题的高精度数值模拟,显著地提升模拟精度和模拟效率。     

一种面向舰船结构毁伤的大变形流固耦合数值计算方法

水下爆炸导致舰船结构毁伤是一个复杂的非线性大变形流固耦合过程,高精度的流固耦合计算是获得高置信模拟结果的关键。基于浸没边界思想,本文提出一种面向大变形壳理论的流固耦合数值方法,可精确刻画流固耦合界面并高效求解流固界面约束方程。基于该方法,本文提出了完整的适用于水下爆炸舰船结构毁伤的大变形流固耦合数值计算方案,并基于大规模并行编程框架,研发形成适用于舰船结构毁伤的流固耦合大规模并行计算软件。与泰勒平板理论解和水下爆炸结构冲击响应实验数据等进行对比表明,本文方法可有效模拟大变形流固耦合工程问题,具备较高数值求解精度。在此基础上,完成了水下爆炸整船结构毁伤过程大规模数值模拟。该方法可有效应用于舰船毁伤等级评估,应用前景广阔。     

过渡金属硫族化合物柔性基底体系的模型与应用

柔性基底体系是晶体外延生长领域于20世纪90年代提出的概念.其核心思想是利用超薄的基底,使其在外延生长时能同时与外延晶膜发生应变,以抵消二者之间的晶格失配,从而减少外延晶膜中的位错,提高晶膜的质量.但是人工制备性能优良的超薄基底往往需要较为复杂的工艺.另一方面,过渡金属硫族化合物由于其层状结构特性和层间较弱的范德瓦耳斯相互作用,是天然的柔性基底.本文介绍近几年来新发展的过渡金属硫族化合物柔性基底体系的模型及应用.以Au-MoS₂作为柔性基底外延生长的原型,结合密度泛函理论、线性弹性理论以及位错理论构建模型,并根据计算结果解释了早先利用透射电子显微镜观测到的Au薄膜在MoS₂上外延生长的相关实验现象.此外,本文还介绍了受到该理论模型启发的相关实验工作,特别是利用Au薄膜分离大面积、单层、高质量MoS₂的技术.最后,讨论了在该领域内值得关注和进一步探索的理论问题.     

非结构多面体二阶局部保界全局重映算法

提出一种三维非结构多面体二阶保界全局重映算法.在旧网格上选取模板利用最小二乘构造插值多项式,采用凸包算法计算多面体相交部分,最后使用局部保界修正技术修补重映后的越界量.多项数值实验表明这种格式同时具有高精度、高分辨率和高效率的特点.