CN107239586A - 对无条件稳定时域有限差分法有效的区域分解并行方法 - Google Patents

对无条件稳定时域有限差分法有效的区域分解并行方法 Download PDF

Info

Publication number
CN107239586A
CN107239586A CN201610188141.XA CN201610188141A CN107239586A CN 107239586 A CN107239586 A CN 107239586A CN 201610188141 A CN201610188141 A CN 201610188141A CN 107239586 A CN107239586 A CN 107239586A
Authority
CN
China
Prior art keywords
mrow
electric field
field value
msubsup
region
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.)
Granted
Application number
CN201610188141.XA
Other languages
English (en)
Other versions
CN107239586B (zh
Inventor
陈如山
陈星光
丁大志
樊振宏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201610188141.XA priority Critical patent/CN107239586B/zh
Publication of CN107239586A publication Critical patent/CN107239586A/zh
Application granted granted Critical
Publication of CN107239586B publication Critical patent/CN107239586B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06F9/00—Arrangements for program control, e.g. control units
    • G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/46—Multiprogramming arrangements
    • G06F9/54—Interprogram communication
    • G06F9/546—Message passing systems or structures, e.g. queues

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种对无条件稳定时域有限差分法有效的区域分解并行方法。本发明利用因果区域分解方法实现了蛙跳交替方向隐格式时域有限差分法的高度并行，保持了无条件稳定性的同时又可高度并行计算。本发明可以有效地节省时域有限差分的仿真计算时间，编程简单具有很强的实际工程应用价值。

Description

对无条件稳定时域有限差分法有效的区域分解并行方法
技术领域
本发明属于隐格式时域有限差分算法数值计算技术，属于大规模并行计算领域，具体是一种对无条件稳定时域有限差分有效的区域分解并行算法。
背景技术
随着计算机硬件技术的飞速发展，单核计算已不能满足目前的需求，多核大规模并行计算成为目前的主流方向。传统的显式时域有限差分法虽然满足天然的并行性可以实现大规模并行计算，但显式时域有限差分法受稳定性条件的限制，无法放大时间步长，导致了仿真时间的加长。隐格式时域有限差分法例如蛙跳交替方向隐格式时域有限差分法不受稳定性条件限制可以放大时间步长，但在求解过程中需求解三对角矩阵，限制了大规模高效并行计算。追求无条件稳定时域有限差分法的大规模高效并行计算一直是研究的热点方向。
发明内容
本发明的目的在于提供一种对无条件稳定时域有限差分法有效的区域分解并行方法。
实现本发明目的的技术解决方案为：一种对无条件稳定时域有限差分法有效的区域分解并行方法，步骤如下：
第一步，消息传递接口MPI初始化，确定进程的总数目以及每个进程的编号；
第二步，执行程序前处理，设置仿真参数，设置吸收边界参数，根据进程的总数目，划分区域，使进程号与子区域一一对应，通过时间步长大小确定缓冲区厚度，建立相邻子区域之间的通信索引；
第三步，根据蛙跳交替方向隐格式时域有限差分的迭代公式进行时间迭代，一次时间步迭代里有两次通信，分别通信缓冲区内正确的电场值和磁场值，该通信只发生在相邻区域之间；迭代结束得到空间中的电场和磁场值；
第四步，数据后处理，根据计算出的场值提取相关的物理参数；
第五步，MPI并行结束。
本发明与现有技术相比，其显著优点为：(1)可以实现隐格式时域有限差分法的大规模并行计算：(2)可以获得与传统显式时域有限差分相比拟的并行效率，同时可以放大时间步长，大大减小了计算时间。
附图说明
图1为计算区域划分示意图。
图2为F15模型示意图。
图3为F15模型雷达散射截面积比较图。
具体实施方式
下面结合附图对本发明作进一步详细描述。
本发明是一种对无条件稳定时域有限差分法有效的区域分解并行方法，步骤如下：
第一步，MPI(消息传递接口)初始化，将三维计算空间按照每个方向进行区域划分，确定进程的总数目，对每个进程进行编号；
第二步，执行程序前处理，对仿真参数进行设置包括吸收边界条件的设置，平面波的加入等。附图1中的实线为步骤1中的计算区域划分的边界，每个计算区域向外进行区域扩展使之与相邻的计算区域有相交的计算区域形成缓冲区域，例如区域5的实际计算区域为位于计算区域4、9、10中虚线所共同围成的计算区域。对每个区域实线与虚线所围成的缓冲区域进行编码，建立相邻计算区域之间的通信索引；
第三步，开始时间迭代
采用公式(1)更新整个计算区域的电场辅助变量采用公式(2)、(3)以及公式(1)计算出的辅助变量的值对整个计算区域的电场值进行更新。整个计算区域的边界采用卷积完全匹配层吸收边界条件，每个子区域的扩展边界采用理想电壁截断。每个子区域的扩展区域为缓冲区，如附图1区域7中周围黑色区域为区域7的缓冲区，缓冲区内正确的电场值分别由相邻子区域计算得到，区域7将内部黑色虚线到缓冲区内边界的电场值分别发送给相邻8个区域，来更新周围八个区域缓冲区内的电场值。区域7缓冲区内的正确电场值从相邻八个区域获得，来更新区域7缓冲区内的场值。子区域所处的空间位置不同通信的子程序个数也不同，区域10与周围3个子区域进行信息传递，区域14的通信对象为周围5个子区域。整个计算区域内所有的子区域在每个时间步均按此方法进行场值更新，空间中每点均可得到正确电场值，至此一个时间步中电场的计算全部完成。
采用公式(4)更新整个计算区域的磁场辅助变量采用公式(5)(6)以及公式(4)计算出的辅助变量值计算出整个空间的磁场值，磁场缓冲区场值的更新与上述电场更新方式相同。重复上述步骤直至时间迭代结束；
第四步，时间迭代结束数据后处理，根据计算出的场值提取相关的物理参数；
第五步，MPI并行结束；
为了验证本发明的正确性与有效性，下面分析了F15战斗机的电磁散射特性。
算例：F15战斗机模型，几何尺寸4.78m×3.35m×1.06m。激励源设置为：调制高斯脉冲，中心频率300MHz，带宽600MHz，入射波的方向θ＝180°,(飞机头部方向照射)，观察角度0°≤θ≤180°,飞机剖分尺寸0.01m，三个方向网格数为578×435×207，时间步长为CFLN＝5，即5倍于传统显示时域有限差分法时间步长。从附图3可以看出本发明方法在300M时的双站RCS与仿真软件CST吻合较好。表1给出了本发明提出方法的并行效率，从表中可以看出本发明实现了无条件稳定时域有限差分法的高度并行。
表1.并行效率

Claims (3)

1.一种对无条件稳定时域有限差分法有效的区域分解并行方法，其特征在于步骤如下：
第一步，消息传递接口MPI初始化，确定进程的总数目以及每个进程的编号；
第二步，执行程序前处理，设置仿真参数，设置吸收边界参数，根据进程的总数目，划分区域，使进程号与子区域一一对应，通过时间步长大小确定缓冲区厚度，建立相邻子区域之间的通信索引；
第三步，根据蛙跳交替方向隐格式时域有限差分的迭代公式进行时间迭代，一次时间步迭代里有两次通信，分别通信缓冲区内正确的电场值和磁场值，该通信只发生在相邻区域之间；迭代结束得到空间中的电场和磁场值；
第四步，数据后处理，根据计算出的场值提取相关的物理参数；
第五步，MPI并行结束。
2.根据权利要求1所述的对无条件稳定时域有限差分法有效的区域分解并行方法，其特征在于：所述步骤二中，执行前处理程序对仿真参数进行设置，确定时间迭代步数和步长，根据进程总数目将时域有限差分FDTD计算区域划分为与进程数相同个数的子区域，对每个子区域根据时间步长的大小进行区域扩展形成缓冲区，对每个子区域的缓冲区的电场和磁场进行整合编码，建立相邻子区域之间的通信索引。
3.根据权利要求1所述的对无条件稳定时域有限差分法有效的区域分解并行方法，其特征在于：所述步骤三中，根据公式(1)(2)(3)计算x、y、z方向电场值
<mrow> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>e</mi> <mi>i</mi> <mi>j</mi> </mrow> <mi>n</mi> </msubsup> <mo>=</mo> <msub> <mi>b</mi> <mi>j</mi> </msub> <msubsup> <mi>&amp;psi;</mi> <mrow> <mi>e</mi> <mi>i</mi> <mi>j</mi> </mrow> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>+</mo> <msub> <mi>a</mi> <mi>j</mi> </msub> <msub> <mi>&amp;delta;</mi> <mi>j</mi> </msub> <msubsup> <mi>H</mi> <mi>k</mi> <mi>n</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>E</mi> <mi>&amp;xi;</mi> <mrow> <mi>n</mi> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>e</mi> <mi>&amp;xi;</mi> <mrow> <mi>n</mi> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>E</mi> <mi>&amp;xi;</mi> <mrow> <mi>n</mi> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
其中，采用公式(1)更新整个计算区域想x、y、z方向的电场辅助变量采用公式(2)、(3)以及公式(1)计算出的辅助变量的值对整个计算区域的电场值进行更新，Δt为时间步长，κη为完全匹配层参数，δη为微分算子，为辅助电场值，为磁场场值；整个计算区域的边界采用卷积完全匹配层吸收边界条件，每个子区域的扩展边界采用理想电壁截断；每个子区域的扩展区域为缓冲区，区域7周围黑色区域为区域7的缓冲区，缓冲区内正确的电场值分别由相邻子区域计算得到，区域7将内部黑色虚线到缓冲区内边界的电场值分别发送给相邻8个区域，来更新周围八个区域缓冲区内的电场值；区域7缓冲区内的正确电场值从相邻八个区域获得，来更新区域7缓冲区内的场值；子区域所处的空间位置不同通信的子程序个数也不同，区域10与周围3个子区域进行信息传递，区域14的通信对象为周围5个子区域；整个计算区域内所有的子区域在每个时间步均按此方法进行场值更新，空间中每点均可得到正确电场值，至此一个时间步中电场的计算全部完成；磁场的处理方式与电场相同，重复上述步骤直至时间迭代结束。
CN201610188141.XA 2016-03-29 2016-03-29 对无条件稳定时域有限差分法有效的区域分解并行方法 Active CN107239586B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610188141.XA CN107239586B (zh) 2016-03-29 2016-03-29 对无条件稳定时域有限差分法有效的区域分解并行方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610188141.XA CN107239586B (zh) 2016-03-29 2016-03-29 对无条件稳定时域有限差分法有效的区域分解并行方法

Publications (2)

Publication Number Publication Date
CN107239586A true CN107239586A (zh) 2017-10-10
CN107239586B CN107239586B (zh) 2020-12-04

Family

ID=59982832

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610188141.XA Active CN107239586B (zh) 2016-03-29 2016-03-29 对无条件稳定时域有限差分法有效的区域分解并行方法

Country Status (1)

Country Link
CN (1) CN107239586B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108228970A (zh) * 2017-12-11 2018-06-29 上海交通大学 结构动力学分析显式异步长并行计算方法
WO2022237218A1 (zh) * 2021-04-22 2022-11-17 东南大学 一种用于ota仿真的复包络蛙跳隐式时域有限差分法
CN117573375A (zh) * 2024-01-15 2024-02-20 上海交通大学四川研究院 一种面向自适应解耦方程的动态负载平衡并行计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006139723A (ja) * 2004-11-15 2006-06-01 Ntt Docomo Inc Ｆｄｔｄ演算装置、ｆｄｔｄ演算方法
CN103969627A (zh) * 2014-05-26 2014-08-06 苏州市数字城市工程研究中心有限公司 基于fdtd的探地雷达大规模三维正演模拟方法
CN104951580A (zh) * 2014-03-28 2015-09-30 南京理工大学 无条件稳定和有条件稳定混合的时域谱元电磁分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006139723A (ja) * 2004-11-15 2006-06-01 Ntt Docomo Inc Ｆｄｔｄ演算装置、ｆｄｔｄ演算方法
CN104951580A (zh) * 2014-03-28 2015-09-30 南京理工大学 无条件稳定和有条件稳定混合的时域谱元电磁分析方法
CN103969627A (zh) * 2014-05-26 2014-08-06 苏州市数字城市工程研究中心有限公司 基于fdtd的探地雷达大规模三维正演模拟方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
T.STEFANSKI 等: "Parallel Implementation of ADI-FDTD on Shared and Distributed Memory Computers", 《THE SECOND EUROPEAN CONFERENCE ON ANTENNAS AND PROPAGATION》 *
吴建斌 等: "一种基于PDD算法的ADI-FDTD算法研究", 《计算机工程与应用》 *
段鑫 等: "基于WinSock和多线程技术的高性能并行FDTD", 《信息与电子工程》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108228970A (zh) * 2017-12-11 2018-06-29 上海交通大学 结构动力学分析显式异步长并行计算方法
CN108228970B (zh) * 2017-12-11 2020-07-14 上海交通大学 结构动力学分析显式异步长并行计算方法
WO2022237218A1 (zh) * 2021-04-22 2022-11-17 东南大学 一种用于ota仿真的复包络蛙跳隐式时域有限差分法
CN117573375A (zh) * 2024-01-15 2024-02-20 上海交通大学四川研究院 一种面向自适应解耦方程的动态负载平衡并行计算方法
CN117573375B (zh) * 2024-01-15 2024-04-02 上海交通大学四川研究院 一种面向自适应解耦方程的动态负载平衡并行计算方法

Also Published As

Publication number Publication date
CN107239586B (zh) 2020-12-04

Similar Documents

Publication Publication Date Title
CN107239586A (zh) 对无条件稳定时域有限差分法有效的区域分解并行方法
CN105654422A (zh) 点云配准方法和系统
Yue et al. A multi‐grid method of high accuracy surface modeling and its validation
CN102930071A (zh) 基于非匹配网格的周期结构的三维电磁场仿真模拟方法
CN107688680A (zh) 一种高效的时域有限元区域分解并行方法
CN102156764A (zh) 一种分析天线辐射和电磁散射的多分辨预条件方法
CN108566208A (zh) 一种层次格网的编码方法及装置
CN116738891B (zh) 一种增强飞行器流场模拟稳定性的lu-sgs改进方法
CN110287450A (zh) 基于积分变分节块法的响应矩阵并行处理实现方法
Elling A possible counterexample to well posedness of entropy solutions and to Godunov scheme convergence
CN105631094A (zh) 基于分段线性循环卷积的一维左手材料Crank-Nicolson完全匹配层实现算法
Wang et al. Research on a surrogate model updating-based efficient multi-objective optimization framework for supertall buildings
CN105205299B (zh) 电大目标电磁散射特性快速降维分析方法
CN105760597A (zh) 基于DG算法的二维色散介质Crank-Nicolson完全匹配层实现算法
CN103279612B (zh) 复杂目标雷达回波快速获取的多重网格预条件方法
CN111339688B (zh) 基于大数据并行算法求解火箭仿真模型时域方程的方法
CN106156475A (zh) 电大尺寸目标的瞬态电磁特性快速提取方法
CN107368652A (zh) 一种基于cndg算法截断等离子体的完全匹配层实现算法
CN106156394B (zh) 基于显式差分格式的电磁特性提取方法
Huang et al. Automatic frame field guided hexahedral mesh generation
CN115329975A (zh) 仿真方法、装置、设备及存储介质
CN106951683A (zh) 一种用于核电厂厂房屏蔽计算的高效并行扫描方法
CN102880590A (zh) 二阶波动方程的非分裂完全匹配层的构造方法
CN107609274B (zh) 基于传输线与级别调度法的二维静磁场并行有限元方法
Cui et al. Development of a parallel multi-dimensional grid-based Vlasov solver for plasma plume simulation

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