CN108052738B - 色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 - Google Patents
色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 Download PDFInfo
- Publication number
- CN108052738B CN108052738B CN201711324508.7A CN201711324508A CN108052738B CN 108052738 B CN108052738 B CN 108052738B CN 201711324508 A CN201711324508 A CN 201711324508A CN 108052738 B CN108052738 B CN 108052738B
- Authority
- CN
- China
- Prior art keywords
- local
- time
- vector
- format
- global
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 22
- 239000002612 dispersion medium Substances 0.000 title claims abstract description 13
- 239000013598 vector Substances 0.000 claims abstract description 57
- 238000000034 method Methods 0.000 claims abstract description 31
- 239000006185 dispersion Substances 0.000 claims abstract description 17
- 239000003245 coal Substances 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 32
- 230000005684 electric field Effects 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000010287 polarization Effects 0.000 claims description 9
- 239000002609 medium Substances 0.000 claims description 8
- 230000001419 dependent effect Effects 0.000 claims description 7
- 230000035699 permeability Effects 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 7
- 230000004907 flux Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 230000005672 electromagnetic field Effects 0.000 claims description 5
- 238000010521 absorption reaction Methods 0.000 claims description 4
- 239000002184 metal Substances 0.000 claims description 4
- 229910052751 metal Inorganic materials 0.000 claims description 4
- 238000012805 post-processing Methods 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000001737 promoting effect Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 13
- 230000000875 corresponding effect Effects 0.000 abstract description 9
- 230000002596 correlated effect Effects 0.000 abstract description 2
- 239000000463 material Substances 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 5
- 238000000576 coating method Methods 0.000 description 4
- 239000011248 coating agent Substances 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000005288 electromagnetic effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 239000010970 precious metal Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 229910052723 transition metal Inorganic materials 0.000 description 1
- 150000003624 transition metals Chemical class 0.000 description 1
- 230000035899 viability Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Discrete Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法,应用于三维电磁分析数值领域;本发明对待分析的包含复杂色散煤质的多尺度目标进行三维建模,建立对应的几何结构模型;并采用四面体网格对建立的几何结构模型进行剖分;通过引入由精细网格相关的新的指数形式未知向量;从而获得全局显式局部无条件稳定的高阶时间格式;本发明的全局显式局部无条件稳定的高阶时间格式中,全局时间步长不再受到精细网格区域严苛的稳定性条件限制,显著增大的全局时间步长可大大减少时间迭代步数,从而有效缩短仿真分析时间。
Description
技术领域
本发明属于三维电磁分析领域,特别涉及一种三维电磁分析数值求解技术。
背景技术
随着隐身技术、物联网电磁兼容及电磁屏蔽等技术的发展,各类复杂色散材料越来越得到重视。从军事隐身技术中针对雷达、红外甚至可见光等不同应用环境的各类隐身涂层材料,到民用物联网设备的抗电磁干扰和电磁屏蔽设计、太阳能薄膜电池等领域中都有大量采用。随着电磁环境的日趋复杂,对这些材料的多频谱性能需求也越来越高。特别对于复杂电磁环境下的多尺度目标而言,其本身在几何结构上具有典型的多尺度特性,其中通常又包含贵金属、过渡金属或复合材料涂层等具有色散或各向异性性质的材料,这些涂层材料不仅会继承目标本身的几何上的多尺度特性,还不得不面临不同材料之间、材料与金属之间材料上的多尺度性(其精确仿真需要在材料交界面上进行网格加密,进而产生所谓的材料上的多尺度性)。这些特性使得包含复杂色散媒质的多尺度目标在复杂电磁环境下将产生复杂的电磁效应,对目标本身的电磁特性、稳定性和可靠性甚至战场生存能力提出严峻挑战。因此,精确地获得这类目标的复杂电磁环境电磁响应特征具有极其重要的意义。
要真实、精确地模拟这类多尺度目标的复杂电磁环境电磁响应特征,有必要在仿真中考虑材料的色散和各向异性特性,加之这些材料具有多频谱特性,其精确时域宽带电磁仿真分析对于传统计算电磁学方法而言非常具有挑战性。近年来兴起的时域间断伽辽金法,在空间离散和时间离散方面非常灵活,并具有高度的并行性,特别适合用于这类复杂问题的三维仿真分析。但是目前时域间断伽辽金方法理论和技术离真正实现包含复杂色散媒质的多尺度目标的高精度、高性能电磁仿真分析仍有较大的差距。对于这类具有典型多尺度特征的大型问题而言,自适应网格加密是保证计算精度和减少计算开销的最重要的途径之一。然而局部网格加密使得显式时间格式的稳定性问题变得更加严峻。如果采用全局显式时间格式,均一化的全局时间步长将带来巨大的计算资源消耗,因为最大时间步长必须满足最小网格的稳定性条件,导致全局时间步长非常小,时间迭代步数大大增加。而这个由最小网格确定的时间步长,对于大尺寸网格并不必要。因此在局部网格加密问题中,全局时间步长将在大尺寸网格区域带来大量的不必要计算花销。全隐式时间格式似乎能够解决这个问题,其无条件稳定的性质使得时域间断伽辽金法能够采用比显式时间格式更大甚至大得多的时间步长来获得同样精度的结果。但是,其代价是需要在每次时间迭代求解全局线性方程组或者对其系数矩阵求逆,巨大的计算资源消耗使得全隐式时域间断伽辽金法很难用于多尺度复杂电磁环境问题这类实际大型问题的三维电磁仿真分析。因此,迫切需要针对多尺度目标的复杂电磁环境电磁响应特征的时域间断伽辽金分析研究更加高效时间格式。
发明内容
为解决传统计算电磁学方法对包含复杂色散媒质的多尺度目标进行高效高精度三维电磁仿真分析困难的问题,本发明提出一种色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法,在几何细节部分具有无条件稳定特性,能够显著增大几何细节部分的时间步长,减少时间迭代次数,具有高效的全局显式特性,大大缩短仿真时长。
本发明采用的技术方案为:色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法,包括:
S1、对待分析的包含复杂色散煤质的多尺度目标进行三维建模,建立对应的几何结构模型;
S2、采用四面体网格对建立的几何结构模型进行剖分,得到若干四面体网格;
S3、构造边值问题,根据伽辽金过程,推导出各四面体网格中局部边值问题的间断伽辽金弱形式;
S4、采用基函数对间断伽辽金弱形式进行离散,得到间断伽辽金半离散格式;
S5、将所有场分量整合为一个未知向量,根据得到的间断伽辽金半离散格式构造一个常微分方程;
S6、根据网格大小,将步骤S2得到的若干四面体网格划分为精细网格和粗糙网格,并根据精细网格和粗糙网格对常微分方程进行未知量分离;
S7、引入由精细网格相关的新的指数形式未知向量替换步骤S6中的未知量,经处理,得到全局显式局部无条件稳定的高阶时间格式;
S8、对得到全局显式局部无条件稳定的高阶时间格式进行推进,获得时域电磁场解,并进行相应的后处理获得所需时域或频域电磁响应特征。
进一步地,步骤S3所述构造边值问题,具体为:根据时域Maxwell方程,通过引入极化电流向量得到混合的Maxwell-Drude方程;采用理想电壁描述金属边界条件;并采用Silver-Muller吸收边界条件截断计算区域;
更进一步地,所述各四面体网格中局部边值问题的间断伽辽金弱形式,具体为:
其中,t为时间变量,ε∞为频率无穷大时的相对介电常数,ωd为等离子体频率,γd为碰撞频率,为第i个网格单元给定面上的外法向单位矢量,为第i个网格单元内部的磁场矢量,为第i个网格单元内部的电场矢量,为第i个网格单元内部的极化电流向量,为第j个网格单元给定面上的外法向单位矢量,为第j个网格单元内部的电场,Fh为所有的面集合,为第i个网格单元对应的基函数矢量,ε0为真空介电常数,μ为计算区域中介质的磁导率。
进一步地,步骤S4所述间断伽辽金半离散格式,具体为:
其中,第i个网格单元内部的未知电场矢量,第i个网格单元内部的未知磁场矢量,第i个网格单元内部的未知极化电流矢量,ε∞为频率无穷大时的相对介电常数,ωd为等离子体频率,γd为碰撞频率,μ为计算区域中介质的相对磁导率,ε0为真空介电常数,t为时间变量,Vi表示四面体Ki中交界面的编号集合,表示电场相关质量矩阵,表示磁场相关质量矩阵,Sii表示交界面Γij的自作用通量矩阵,Sij表示第i个网格单元与第j个网格单元的互作用通量矩阵,其中Γij表示第i个网格单元与第j个网格单元的交界面。
进一步地,步骤S7具体为:
S71、引入由精细网格相关的新的指数形式未知向量替换步骤S6中的未知量,获得局部无条件稳定的指数时间积分格式;
S72、采用高阶显式时间格式对步骤S71中获得的局部无条件稳定的指数时间积分格式进行时间离散,得到全局显式局部无条件稳定的高阶时间格式。
本发明的有益效果:本发明的色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法,引入由精细网格相关的新的指数形式未知向量替换常微分方程中的未知向量,得到局部无条件稳定的指数时间积分格式,采用诸如低存储龙格库塔格式的高阶显式时间格式对局部无条件稳定的指数时间积分格式进行时间离散,获得全局显式局部无条件稳定的高阶时间格式;该全局显式局部无条件稳定的高阶时间格式中,全局时间步长不再受到精细网格区域严苛的稳定性条件限制,显著增大的全局时间步长可大大减少时间迭代步数,从而有效缩短仿真分析时间;从而构造出的高阶局部隐式时间格式有效提升了包含复杂色散媒质的多尺度目标的电磁仿真分析效率。
附图说明
图1为本发明方案的流程图。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
如图1所示为本发明的方案流程图,本发明的技术方案为:色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法,包括:
S1、对待分析的包含复杂色散煤质的多尺度目标进行三维建模,建立对应的几何结构模型。选取待分析的包含复杂色散媒质的多尺度目标,例如军车、战斗机和无人机。根据目标几何特征和材料属性,建立三维模型,同时加入截断吸收边界条件,形成计算区域的几何结构模型。
S2、采用四面体网格对建立的几何结构模型进行剖分,得到若干四面体网格;具体的:采用四面体网格剖分步骤S1所建几何结构模型计算电磁学方法中的一种公知过程,因此本步骤不再详细描述。剖分后的计算区域被分割为若干个四面体网格,从而连续的几何结构空间被转化为离散的网格空间。
S3、构造边值问题,根据伽辽金过程,推导出各四面体网格中局部边值问题的间断伽辽金弱形式;
所述构造边值问题,具体为:首先,从如下时域Maxwell方程出发:
对于复杂色散媒质,本发明实施例采用Drude色散模型进行描述。该色散模型中,复电位移可表示为
其中
这里ε∞为频率无穷大时的相对介电常数,ωd为等离子体频率,γd为碰撞频率,三者皆为Drude模型参数,与材料特性直接相关。将式(2)进行傅里叶变换,并带入到时域Maxwell方程,可得
又根据其定义式(3)可求得其傅里叶变换
本发明实施例采用理想电壁描述金属边界条件:
采用Silver-Muller吸收边界条件截断计算区域:
至此,式(6)、式(7)和式(8)共同构成了本方法的边值问题。
其中,t为时间变量,ε∞为频率无穷大时的相对介电常数,ωd为等离子体频率,γd为碰撞频率,为第i个网格单元给定面上的外法向单位矢量,为第i个网格单元内部的磁场矢量,为第i个网格单元内部的电场矢量,为第i个网格单元内部的极化电流向量,为第j个网格单元给定面上的外法向单位矢量,为第j个网格单元内部的电场,Fh为所有的面集合,为第i个网格单元对应的基函数矢量,ε0为真空介电常数,μ为计算区域中介质的磁导率。伽辽金过程为一种公知过程,在有限元法和时域间断伽辽金法中广泛应用,这里不再详细描述。
S4、采用基函数对间断伽辽金弱形式进行离散,得到间断伽辽金半离散格式;具体为:间断伽辽金法支持多种类型的基函数,这里以高阶叠层矢量基函数为例。在四面体网格单元内部,局部电磁场可由Φil基函数的线性组合表示为
其中,di为局部未知量的个数,由基函数阶数和类型确定,eil和hil为基函数系数。那么,根据上式,对步骤D中得到的间断伽辽金弱形式进行离散,可得到四面体网格Ki中的局部间断伽辽金半离散格式
其中,第i个网格单元内部的未知电场矢量,第i个网格单元内部的未知磁场矢量,第i个网格单元内部的未知极化电流矢量,ε∞为频率无穷大时的相对介电常数,ωd为等离子体频率,γd为碰撞频率,μ为计算区域中介质的相对磁导率,ε0为真空介电常数,t为时间变量,Vi表示四面体Ki中交界面的编号集合,表示电场相关质量矩阵,表示磁场相关质量矩阵,Sii表示交界面Γij的自作用通量矩阵,Sij表示第i个网格单元与第j个网格单元的互作用通量矩阵。此外
Γij表示第i个网格单元与第j个网格单元的交界面。
其中,每个对角块分别对应于一个局部质量矩阵(电场相关质量矩阵)或(磁场相关质量矩阵);矩阵K为块矩阵,其对角块为非对角块为将场分量和整合到一个未知向量u,则全局半离散格式(12)可转化为如下常微分方程形式
其中
S6、根据网格大小,将步骤S2得到的若干四面体网格划分为精细网格和粗糙网格,并根据精细网格和粗糙网格对常微分方程进行未知量分离;具体为:
根据网格尺寸大小,将步骤C中得到的网格划分为精细网格和粗糙网格两类。那么根据此分类情况,步骤F中常微分方程的未知量可分离为
其中矩阵P是对角矩阵,对角元素为0或1,这里1用来标记精细网格部分相关未知量,I为单位对角矩阵。则上述常微分方程(13)可分离为
其中Cf=CP和Cc=C(I-P)分别由精细和粗糙网格相关未知量构成。
S7、引入由精细网格相关的新的指数形式未知向量替换步骤S6中的未知量,经处理,得到全局显式局部无条件稳定的高阶时间格式;包括以下分步骤:
S71、引入由精细网格相关的新的指数形式未知向量替换步骤S6中的未知向量,去除精细网格部分对该常微分方程时间格式的显式依赖,获得局部无条件稳定的指数时间积分格式;
引入由精细网格相关的新的指数形式未知向量
而右端项亦可转化为
则化简后,得到局部无条件稳定的指数时间积分格式
至此,通过引入式(17)的新变量,去除了时间格式对精细网格部分相关矩阵Cf的显式依赖,Cf不需要单独进行时间迭代。并且,和的乘积为单位矩阵,因此式(20)中矩阵是粗糙网格部分相关矩阵Cc的相似矩阵,具有相同的特征谱。换而言之,精细网格部分不再影响式(20)的稳定性。那么式(20)进一步时间离散之后,对精细网格部分而言,所获得时间格式是无条件稳定的。
S72、采用高阶显式时间格式对步骤S71中获得的指数时间积分格式进行时间离散,得到全局显式局部无条件稳定的高阶时间格式;
采用诸如低存储龙格库塔格式的高阶显式时间格式对式(20)进行时间离散,可获得全局显式局部无条件稳定的高阶时间格式。该时间格式中,全局时间步长不再受到精细网格区域严苛的稳定性条件限制,显著增大的全局时间步长可大大减少时间迭代步数,从而有效缩短仿真分析时间。因此,通过该方式构造出的高阶局部隐式时间格式将有效提升包含复杂色散媒质的多尺度目标的电磁仿真分析效率。
S8、对得到全局显式局部无条件稳定的高阶时间格式进行推进,获得时域电磁场解,并进行相应的后处理获得所需时域或频域电磁响应特征。具体的:对步骤S7中得到全局显式局部无条件稳定的高阶时间格式进行时间推进。每次时间迭代均可获得相应时刻的各类场分量,通过后处理计算可得所需瞬态或频域电磁响应特征。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (1)
1.色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法,其特征在于,包括:
S1、对待分析的包含复杂色散煤质的多尺度目标进行三维建模,建立对应的几何结构模型;
S2、采用四面体网格对建立的几何结构模型进行剖分,得到若干四面体网格;
S3、构造边值问题,根据伽辽金过程,推导出各四面体网格中局部边值问题的间断伽辽金弱形式;步骤S3所述构造边值问题,具体为:
S31:采用理想电壁描述金属边界条件;并采用Silver-Muller吸收边界条件截断计算区域;根据时域Maxwell方程,通过引入极化电流向量得到混合的Maxwell-Drude方程;其中,为辅助电位移矢量,t为时间变量;
S32:根据伽辽金过程,推导出各四面体网格中局部边值问题的间断伽辽金弱形式,具体为:
其中,t为时间变量,ε∞为频率无穷大时的相对介电常数,ωd为等离子体频率,γd为碰撞频率,为第i个网格单元给定面上的外法向单位矢量,为第i个网格单元内部的磁场矢量,为第i个网格单元内部的电场矢量,为第i个网格单元内部的极化电流向量,为第j个网格单元给定面上的外法向单位矢量,为第j个网格单元内部的电场矢量,Fh为所有的面集合,为第i个网格单元对应的基函数矢量,ε0为真空介电常数,μ为计算区域中介质的磁导率;
其中,di为局部未知量的个数,eil和hil为基函数系数;采用该基函数对间断伽辽金形式进行离散,所得间断伽辽金半离散格式具体为:
其中, 第i个网格单元内部的未知电场矢量,第i个网格单元内部的未知磁场矢量,第i个网格单元内部的未知极化电流矢量,ε∞为频率无穷大时的相对介电常数,ωd为等离子体频率,γd为碰撞频率,μ为计算区域中介质的相对磁导率,ε0为真空介电常数,t为时间变量,Vi表示四面体Ti中交界面的编号集合,表示电场相关质量矩阵,表示磁场相关质量矩阵,Sii表示交界面Γij的自作用通量矩阵,Sij表示第i个网格单元与第j个网格单元的互作用通量矩阵,其中Γij表示第i个网格单元与第j个网格单元的交界面;
S6、根据网格大小,将步骤S2得到的若干四面体网格划分为精细网格和粗糙网格,并根据精细网格和粗糙网格对常微分方程进行未知量分离;
S7、引入由精细网格相关的新的指数形式未知向量替换步骤S6中的未知量,经处理,得到全局显式局部无条件稳定的高阶时间格式;步骤S7具体为:
S71、引入由精细网格相关的新的指数形式未知向量替换步骤S6中的未知量,获得局部无条件稳定的指数时间积分格式;具体为:
所述新的指数形式未知向量表达式为:
全局显式局部无条件稳定的高阶时间格式表达式为:
其中,Cf表示精细网格部分相关矩阵,Cc表示粗糙网格部分相关矩阵;
S72、采用高阶显式时间格式对步骤S71中获得的局部无条件稳定的指数时间积分格式进行时间离散,得到全局显式局部无条件稳定的高阶时间格式;
S8、对得到全局显式局部无条件稳定的高阶时间格式进行推进,获得时域电磁场解,并进行相应的后处理获得所需时域或频域电磁响应特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711324508.7A CN108052738B (zh) | 2017-12-13 | 2017-12-13 | 色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711324508.7A CN108052738B (zh) | 2017-12-13 | 2017-12-13 | 色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108052738A CN108052738A (zh) | 2018-05-18 |
CN108052738B true CN108052738B (zh) | 2021-10-15 |
Family
ID=62132228
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711324508.7A Active CN108052738B (zh) | 2017-12-13 | 2017-12-13 | 色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108052738B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109190169B (zh) * | 2018-08-02 | 2022-07-29 | 电子科技大学 | 一种三维时域电磁学杂交时域间断伽辽金数值方法 |
CN109657288B (zh) * | 2018-11-28 | 2022-07-26 | 电子科技大学 | 一种三维显隐时域电磁学数值方法 |
CN109726433B (zh) * | 2018-11-30 | 2022-10-14 | 电子科技大学 | 基于曲面边界条件的三维无粘低速绕流的数值模拟方法 |
CN109492341B (zh) * | 2018-12-25 | 2023-03-24 | 南京邮电大学 | 表面等离激元波导的光热效应仿真方法 |
CN109684740B (zh) * | 2018-12-27 | 2022-10-14 | 电子科技大学 | 一种基于混合网格及时间步长的电磁学多尺度计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2784798A1 (en) * | 2013-03-27 | 2014-10-01 | Max-Planck-Gesellschaft zur Förderung der Wissenschaften e.V. | Determining an electromagnetic response of a sample |
CN104268322A (zh) * | 2014-06-27 | 2015-01-07 | 北京航空航天大学 | Weno差分方法的一种边界处理技术 |
CN105760597A (zh) * | 2016-02-03 | 2016-07-13 | 天津工业大学 | 基于DG算法的二维色散介质Crank-Nicolson完全匹配层实现算法 |
CN106295120A (zh) * | 2016-07-26 | 2017-01-04 | 上海无线电设备研究所 | 一种目标电磁散射特性快速计算方法 |
CN107305536A (zh) * | 2016-04-22 | 2017-10-31 | 南京理工大学 | 混合阶时域不连续伽略金方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2378444B1 (de) * | 2010-04-13 | 2015-07-29 | CST-Computer Simulation Technology AG | Verfahren, Vorrichtung und Computerprogrammprodukt zur Bestimmung eines elektromagnetischen Nahfeldes einer Feldanregungsquelle eines elektrischen Systems |
US9279314B2 (en) * | 2011-08-11 | 2016-03-08 | Conocophillips Company | Heat front capture in thermal recovery simulations of hydrocarbon reservoirs |
-
2017
- 2017-12-13 CN CN201711324508.7A patent/CN108052738B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2784798A1 (en) * | 2013-03-27 | 2014-10-01 | Max-Planck-Gesellschaft zur Förderung der Wissenschaften e.V. | Determining an electromagnetic response of a sample |
CN104268322A (zh) * | 2014-06-27 | 2015-01-07 | 北京航空航天大学 | Weno差分方法的一种边界处理技术 |
CN105760597A (zh) * | 2016-02-03 | 2016-07-13 | 天津工业大学 | 基于DG算法的二维色散介质Crank-Nicolson完全匹配层实现算法 |
CN107305536A (zh) * | 2016-04-22 | 2017-10-31 | 南京理工大学 | 混合阶时域不连续伽略金方法 |
CN106295120A (zh) * | 2016-07-26 | 2017-01-04 | 上海无线电设备研究所 | 一种目标电磁散射特性快速计算方法 |
Non-Patent Citations (3)
Title |
---|
A New Family of Exponential-Based High-Order DGTD Methods for Modeling 3-D Transient Multiscale Electromagnetic Problems;Hao Wang et al.;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20171130;第65卷(第11期);第5960-5974页 * |
High Order DGTD Solver for the Numerical Modeling of Nanoscale Light/Matter Interaction;Stéphane Lanteri et al.;《Spectral and High Order Methods for Partial Differential Equations ICOSAHOM 2016》;20170831;第243-255页 * |
Model order reduction based solver for discontinuous Galerkin element approximation of time-domain Maxwell’s equations in dispersive media;Kun Li et al.;《IMACS2016 - 20th IMACS WORLD CONGRESS》;20161231;第1-55页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108052738A (zh) | 2018-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108052738B (zh) | 色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 | |
CN107992696B (zh) | 一种复杂色散媒质中的改进的指数时间积分构造方法 | |
Musco et al. | Primordial black hole formation in the early universe: critical behaviour and self-similarity | |
Pétri | The pulsar force-free magnetosphere linked to its striped wind: time-dependent pseudo-spectral simulations | |
Abdelrahman et al. | A robust and accurate solver for some nonlinear partial differential equations and tow applications | |
CN109684740B (zh) | 一种基于混合网格及时间步长的电磁学多尺度计算方法 | |
Poole et al. | Control point-based aerodynamic shape optimization applied to AIAA ADODG test cases | |
Andrade et al. | On the stability of gravity with Dirichlet walls | |
CN113158527B (zh) | 一种基于隐式fvfd计算频域电磁场的方法 | |
CN111079278B (zh) | 三维时域杂交间断伽辽金方法外加电磁源项的处理方法 | |
CN107944141A (zh) | 基于杂交时域间断伽辽金法的时域计算电磁学数值方法 | |
CN106294894B (zh) | 快速分析非均匀目标电磁散射特性的有限元边界积分方法 | |
CN107908903B (zh) | 复杂色散媒质中指数时间积分法的矩阵指数时间去耦方法 | |
CN104915326B (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
El-Monier et al. | Dust ion acoustic shock and solitary waves in a magnetized multi-component plasma with superthermal electrons | |
Cheong et al. | General-relativistic Radiation Transport Scheme in Gmunu. I. Implementation of Two-moment-based Multifrequency Radiative Transfer and Code Tests | |
Chen et al. | A three-dimensional HIE-PSTD scheme for simulation of thin slots | |
CN107943755B (zh) | 复杂色散媒质中指数时间积分法的矩阵指数降维方法 | |
CN116776695A (zh) | 基于超高阶有限元技术的一维电磁计算方法、系统及设备 | |
Ding et al. | An accelerated surface discretization‐based BEM approach for non‐homogeneous linear problems in 3‐D complex domains | |
Loutsenko et al. | New exact results in spectra of stochastic Loewner evolution | |
Cheamsawat et al. | Free energy dependence on spatial geometry for (2+ 1)-dimensional QFTs | |
Ghosh et al. | Conservation laws and new exact solutions to the maccari’s modulation equations | |
Alhaidari et al. | The Dirac–Coulomb problem: a mathematical revisit | |
Song et al. | Localized solitary waves and their dynamical stabilities in magnetized dusty plasma |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |