多级局部时间步进技术的任意高阶混合网格时域不连续伽辽金方法技术

技术编号:25551070 阅读:30 留言:0更新日期:2020-09-08 18:50
本发明专利技术公开了一种多级局部时间步进技术的任意高阶混合网格时域不连续伽辽金方法,该方法选取麦克斯韦方程组为基础数值模型,结合任意高阶导数(ADER)时间步进方案,对计算域采用合理的四面体/六面体混合网格进行剖分,各剖分元胞根据稳定性条件分别地自动决定合适的时间迭代步长,可实现任意多个、任意比例的时间迭代步长大小,各元胞电磁场量按照自己的时间迭代步长进行迭代更新,直至所有元胞场量迭代到规定的时间点,对得到的时变电磁场量进行后处理,得到相应的S参数、雷达散射截面积和电磁场空间分布。本发明专利技术缓解了时域电磁分析方法的时间步长受限于最小离散网格尺寸所带来计算效率低的问题,不仅提高了计算精度,还减少了计算时间,特别适用于空间多尺度电磁问题的快速分析。

【技术实现步骤摘要】
多级局部时间步进技术的任意高阶混合网格时域不连续伽辽金方法
本专利技术属于电磁仿真
,具体是一种不连续伽辽金时域有限元算法数值计算技术,是一种模拟空间多尺度电磁问题的高效算法。
技术介绍
随着工程设计和电磁仿真环境的日益复杂,当多个时空尺度在感兴趣的问题上具有研究意义时,由于传统的数值计算技术缺乏足够的效率、精确性和鲁棒性来求解宏观多尺度Maxwell方程,因此对其电磁瞬态问题进行分析时将会面临着困难和挑战。
技术实现思路
本专利技术的目的在于提出一种多级局部时间步进技术的任意高阶混合网格时域不连续伽辽金方法。实现本专利技术目的的技术解决方案为:一种多级局部时间步进技术的任意高阶混合网格时域不连续伽辽金方法,步骤如下:第一步,选取含有局部精细或高电性结构的空间多尺度问题作为电磁仿真模型,采用合适尺寸的四面体和六面体单元对相应的区域进行空间离散,获取模型的结构信息;第二步,以麦克斯韦方程组作为基础控制方程,采用不连续伽辽金技术建立矩阵方程,时间上采用ADER迭代公式进行时间离散,得到显示求本文档来自技高网...

【技术保护点】
1.一种多级局部时间步进技术的任意高阶混合网格时域不连续伽辽金方法,其特征在于:/n第一步,选取含有局部精细或高电性结构的空间多尺度问题作为电磁仿真模型,采用四面体和六面体单元对相应的区域进行空间离散,获取模型的结构信息;/n第二步,以麦克斯韦方程组作为基础控制方程,采用不连续伽辽金技术建立矩阵方程,时间上采用任意高阶导数ADER迭代公式进行时间离散,得到显示求解矩阵迭代方案;/n第三步,根据ADER显式时间方案近似遵循的稳定性条件,依次将各个剖分单元按照任意整数比例划分至多个不同的计算域,自动的确定各自区域满足的迭代时间步长;/n第四步,计算域各自按时间步长从小到大依次进行时间迭代;/n第五...

【技术特征摘要】
1.一种多级局部时间步进技术的任意高阶混合网格时域不连续伽辽金方法,其特征在于:
第一步,选取含有局部精细或高电性结构的空间多尺度问题作为电磁仿真模型,采用四面体和六面体单元对相应的区域进行空间离散,获取模型的结构信息;
第二步,以麦克斯韦方程组作为基础控制方程,采用不连续伽辽金技术建立矩阵方程,时间上采用任意高阶导数ADER迭代公式进行时间离散,得到显示求解矩阵迭代方案;
第三步,根据ADER显式时间方案近似遵循的稳定性条件,依次将各个剖分单元按照任意整数比例划分至多个不同的计算域,自动的确定各自区域满足的迭代时间步长;
第四步,计算域各自按时间步长从小到大依次进行时间迭代;
第五步,按照规定时间步数完成电场和磁场的迭代求解,提取观察面上的电磁场信息,获取系统相应的参数和空间分布。


2.根据权利要求1所述的多级局部时间步进技术的任意高阶混合网格时域不连续伽辽金方法,其特征在于:所述步骤二中,遵循不连续伽辽金时域有限元方法的标准分析步骤,得到显式的半离散形式矩阵方程:



将上述方程记为:



其中,u=[e,h]T是所需求解的未知向量,是块对角特性的逆矩阵,是系统的右边矩阵。
各个矩阵元素具体的表达式:


















上式中ε、μ分别表示离散单元的介电常数和磁导率,Ni、Nj为有限元四面体或六面体矢量叠层测试基函数和场量的展开基函数,Nj+表示相邻体测试基函数,波阻抗和导纳的关系为
时间上采用ADER迭代格式完成离散,每个离散单元上场值的时间偏导项进行Taylor级数展开:

...

【专利技术属性】
技术研发人员:陈如山丁大志樊振宏包华广何姿王蒙蒙
申请(专利权)人:南京理工大学
类型:发明
国别省市:江苏;32

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1