CN101770526A - 在pc机上的并行化实现分子污染直接模拟蒙特卡罗方法 - Google Patents
在pc机上的并行化实现分子污染直接模拟蒙特卡罗方法 Download PDFInfo
- Publication number
- CN101770526A CN101770526A CN200810188075A CN200810188075A CN101770526A CN 101770526 A CN101770526 A CN 101770526A CN 200810188075 A CN200810188075 A CN 200810188075A CN 200810188075 A CN200810188075 A CN 200810188075A CN 101770526 A CN101770526 A CN 101770526A
- Authority
- CN
- China
- Prior art keywords
- node
- molecule
- collision
- molecules
- zone
- 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
Links
Images
Abstract
一种在PC机上并行实现分子污染直接模拟蒙特卡罗的方法,具体包括以下步骤:(1)在各节点上进行初始化工作,每个节点上的区域信息互不相关;(2)各个节点计算分子的自由运动及与边界面的碰撞,其中与壁面碰撞计算仅在包含壁面的节点上进行;(3)各节点将所有分子重新排序;(4)各节点计算分子间的碰撞过程;(5)各节点进行宏观量的统计,之后进行本节点下一步取样,即返回(2);(6)当所有节点取样次数达到预设置,各个节点将计算结果回传到同一节点上,再由该节点进行输出,结束计算。在这种情况下,并行计算显示了容量大、时间短的优越性,不仅解决单个计算机内存不足的难题,而且大大地提高了运算速度。
Description
技术领域
本发明涉及一种应用于航天器上的分子污染直接模拟蒙特卡罗方法,特别涉及在PC机上的并行化实现分子污染直接模拟蒙特卡罗方法。
背景技术
用于卫星整星污染模拟的直接模拟蒙特卡罗(Direct Simulation MonteCarlo,DSMC)方法是一种昂贵的计算方法,所占内存大、运算时间长,并且,工程应用中,航天器自身的几何结构决定了模型巨大的计算量,传统的台式计算机串行算法很难满足运算要求。在这种情况下,并行计算显示了容量大、时间短的优越性,不仅解决单个计算机内存不足的难题,而且大大地提高了运算速度。这是因为:1它可以加快速度,即在更短的时间内解决相同的问题或在相同的时间内解决更多更复杂的问题,特别是对一些新出现的巨大挑战问题,不使用并行计算时根本无法解决的;2节省投入,并行计算可以以较低的投入完成串行计算才能够完成的任务;3物理极限的约束,光速是不可逾越的速度极限,设备和材料也不可能做得无限小,只有通过并行才能够不断提高速度。
发明内容
一种在PC机上并行实现分子污染直接模拟蒙特卡罗的方法,基于DSMC方法,在整个计算区域采取区域分裂法,将区域按照Z向划分为n个节点,各个节点只处理本节点包含的区域,计算内容完全相同,包括分子的自由运动和与壁面的碰撞、分子由来流进入、分子间的碰撞以及宏观参数的统计;各节点上的网格号保持不变,同时确定分子进入的区域,直接将速度、转动能量、振动能量及分子种类标识传入该节点,节点间的信息交换只存在于相邻节点之间;具体包括以下步骤:(1)在各节点上进行初始化工作,每个节点上的区域信息互不相关;(2)各个节点计算分子的自由运动及与边界面的碰撞,其中与壁面碰撞计算仅在包含壁面的节点上进行,当某分子的网格超出该节点的计算边界,则说明该分子已经离开原来的计算区域而进入新的计算区域,根据分子的网格编号确定出分子进入的新的计算区域;在分子原来所在区域内,将该分子作真空边界处理,而在新的区域内,将该分子作为该区域的N+1,该区域内的分子数目为N个分子处理;此处N为正整数;同时,将该分子离开原节点的分子信息发送至新节点;(3)各节点将所有分子重新排序;(4)各节点计算分子间的碰撞过程;(5)各节点进行宏观量的统计,之后进行本节点下一步取样,即返回(2);(6)当所有节点,即从节点取样次数达到预设置,各个节点将计算结果回传到主节点上,再由主节点进行输出,结束计算。
附图说明
根据下面提供的示例实施例的详细说明和附图,本发明可被更加充分地理解,所述详细说明和附图仅仅是为了说明的目的,因而并非是限制本发明。
图1DSMC方法外部并行消息流程图;
图2DSMC的并行方案流程图。
具体实施方式
用于卫星整星污染的DSMC方法算例计算量很大,单台PC机已不能满足工程化要求。因此考虑采用串行程序并行化的一个最基本原则就是要有优良的计算通讯比,即计算量大于通讯量,这样才能有利可图;否则,虽然计算量分配在多个计算机上,但计算机之间的通讯却占据了大部分的时间,结果反而并行计算所耗费的时间大于串行计算,这样的并行是没有任何意义的。另外,并行计算不能够影响到计算结果的准确性,其结果必须与串行程序的计算结果相一致,否则,并行效率再高,也不具任何意义。DSMC方法外部并行消息流程图如图1所示。
本发明对单个程序进行并行化处理。方案的流程图如图2所示。该方案结合DSMC方法的特点,在整个计算区域采取区域分裂法。考虑到网格是按Z向一层一层划分的,因此,区域按照Z向划分为n个节点。各个节点只处理本节点包含的区域,计算内容完全相同,包括分子的自由运动和与壁面的碰撞、分子由来流进入、分子间的碰撞以及宏观参数的统计等。各节点上的网格号保持不变,所以该方案除了各节点之间存在信息传递外,在区域信息及分子编号上无任何关系。这样该方案只用判断分子是否离开本区域(节点)而进入其它区域(节点),同时确定分子进入的区域,直接将速度、转动能量、振动能量及分子种类标识等传入该节点。由于DSMC方法中时间步长是按照网格尺寸选取的,即要保证分子在一个时间步长内走过的路程最多不能超过网格尺寸,一般为网格尺寸的1/3,所以分子离开本节点计算区域后,只可能进入其相邻的节点而不会跨区,这样,节点间的信息交换只存在于相邻节点之间。
具体步骤如下:
(1)在各节点上进行初始化工作,每个节点上的区域信息(包括边界信息,网格信息)互不相关。
(2)各个节点计算分子的自由运动及与边界面的碰撞(与壁面碰撞计算仅在包含壁面的节点上进行),当某分子的网格超出该节点的计算边界,则说明该分子已经离开原来的计算区域而进入新的计算区域,根据分子的网格编号确定出分子进入的新的计算区域。在分子原来所在区域内,将该分子作真空边界处理,而在新的区域内,将该分子作为该区域的N+1(该区域内的分子数目为N)个分子处理。同时,将该分子离开原节点的分子信息(运动后的位置坐标、速度、转动能量、平动能量和分子所处的网格号等)发送至新节点。
(3)各节点将所有分子重新排序。
(4)各节点计算分子间的碰撞过程。
(5)各节点进行宏观量的统计,之后进行本节点下一步取样,即返回(2)。
(6)当所有节点,即从节点取样次数达到预设置,各个节点将计算结果回传到主节点上,再由主节点进行输出,结束计算。
Claims (3)
1.一种在PC机上并行实现分子污染直接模拟蒙特卡罗的方法,基于DSMC方法,在整个计算区域采取区域分裂法,将区域按照Z向划分为n个节点,各个节点只处理本节点包含的区域,计算内容完全相同,包括分子的自由运动和与壁面的碰撞、分子由来流进入、分子间的碰撞以及宏观参数的统计;各节点上的网格号保持不变,同时确定分子进入的区域,直接将速度、转动能量、振动能量及分子种类标识传入该节点,节点间的信息交换只存在于相邻节点之间;具体包括以下步骤:
(1)在各节点上进行初始化工作,每个节点上的区域信息互不相关;
(2)各个节点计算分子的自由运动及与边界面的碰撞,其中与壁面碰撞计算仅在包含壁面的节点上进行,当某分子的网格超出该节点的计算边界,则说明该分子已经离开原来的计算区域而进入新的计算区域,根据分子的网格编号确定出分子进入的新的计算区域;在分子原来所在区域内,将该分子作真空边界处理,而在新的区域内,将该分子作为该区域的N+1,该区域内的分子数目为N个分子处理;此处N为正整数;同时,将该分子离开原节点的分子信息发送至新节点;
(3)各节点将所有分子重新排序;
(4)各节点计算分子间的碰撞过程;
(5)各节点进行宏观量的统计,之后进行本节点下一步取样,即返回(2);
(6)当所有节点,即从节点取样次数达到预设置,各个节点将计算结果回传到主节点上,再由主节点进行输出,结束计算。
2.如权利要求1所述的方法,其中的区域信息包括边界信息和网格信息。
3.如权利要求1所述的方法,其中分子信息具体包括运动后的位置坐标、速度、转动能量、平动能量和分子所处的网格号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200810188075 CN101770526B (zh) | 2008-12-29 | 2008-12-29 | 在pc机上的并行化实现分子污染直接模拟蒙特卡罗方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200810188075 CN101770526B (zh) | 2008-12-29 | 2008-12-29 | 在pc机上的并行化实现分子污染直接模拟蒙特卡罗方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101770526A true CN101770526A (zh) | 2010-07-07 |
CN101770526B CN101770526B (zh) | 2013-03-20 |
Family
ID=42503384
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200810188075 Expired - Fee Related CN101770526B (zh) | 2008-12-29 | 2008-12-29 | 在pc机上的并行化实现分子污染直接模拟蒙特卡罗方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101770526B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268120A (zh) * | 2014-09-26 | 2015-01-07 | 安徽大学 | 一种聚合反应动力学过程蒙托卡罗模拟的并行计算方法 |
CN109509100A (zh) * | 2018-10-22 | 2019-03-22 | 广东工业大学 | 一种基于集群的风险价值Monte Carlo并行计算框架设计方法 |
CN111354086A (zh) * | 2018-12-24 | 2020-06-30 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种适用于dsmc方法网格位置属性判断的双向三维扫描方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020041566A1 (en) * | 2000-08-11 | 2002-04-11 | Jen-Shun Yang | Dynamic scheduling for packet data network |
CN101268935A (zh) * | 2008-04-25 | 2008-09-24 | 清华大学 | 可调用多种成像算法的在体荧光分子成像建模方法 |
-
2008
- 2008-12-29 CN CN 200810188075 patent/CN101770526B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020041566A1 (en) * | 2000-08-11 | 2002-04-11 | Jen-Shun Yang | Dynamic scheduling for packet data network |
CN101268935A (zh) * | 2008-04-25 | 2008-09-24 | 清华大学 | 可调用多种成像算法的在体荧光分子成像建模方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268120A (zh) * | 2014-09-26 | 2015-01-07 | 安徽大学 | 一种聚合反应动力学过程蒙托卡罗模拟的并行计算方法 |
CN104268120B (zh) * | 2014-09-26 | 2017-04-12 | 安徽大学 | 一种聚合反应动力学过程蒙特卡罗模拟的并行计算方法 |
CN109509100A (zh) * | 2018-10-22 | 2019-03-22 | 广东工业大学 | 一种基于集群的风险价值Monte Carlo并行计算框架设计方法 |
CN111354086A (zh) * | 2018-12-24 | 2020-06-30 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种适用于dsmc方法网格位置属性判断的双向三维扫描方法 |
CN111354086B (zh) * | 2018-12-24 | 2023-04-14 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种适用于dsmc方法网格位置属性判断的双向三维扫描方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101770526B (zh) | 2013-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2016169286A1 (zh) | 一种离散制造系统的车间布局方法 | |
CN103324780B (zh) | 颗粒流动仿真系统和方法 | |
CN103345580B (zh) | 基于格子Boltzmann方法的并行CFD方法 | |
CN106339351B (zh) | 一种sgd算法优化系统及方法 | |
CN102521854B (zh) | 一种适用于二维流场的并行流线放置方法 | |
Dubey et al. | Challenges of Extreme Computing using the FLASH code | |
CN101770526B (zh) | 在pc机上的并行化实现分子污染直接模拟蒙特卡罗方法 | |
Ho et al. | A design methodology for efficient application-specific on-chip interconnects | |
CN102945198A (zh) | 一种表征高性能计算应用特征的方法 | |
CN102508640A (zh) | 基于任务分解的分布式rfid复杂事件检测方法 | |
CN106202383A (zh) | 一种应用于网络爬虫的网络带宽占比动态预测方法与系统 | |
CN108595255A (zh) | 地理分布式云中基于最短路径算法的工作流任务调度方法 | |
Cui et al. | An efficient algorithm for online soft real-time task placement on reconfigurable hardware devices | |
CN101436204A (zh) | 一种基于并行元胞自动机的城市演变模拟实现方法 | |
Minkenberg et al. | End-to-end modeling and simulation of high-performance computing systems | |
Wright | Solving optimization problems on computational grids | |
Nikitin et al. | Analytical performance modeling of hierarchical interconnect fabrics | |
Liu et al. | Business process modeling and analysis based on queue logical Petri nets | |
CN1870545A (zh) | 自动进行方案对比的可信网络仿真系统 | |
Xie et al. | Modeling traffic of big data platform for large scale datacenter networks | |
Soumya et al. | Area-performance trade-off in floorplan generation of Application-Specific Network-on-Chip with soft cores | |
Li et al. | Analytical performance modeling for null message-based parallel discrete event simulation | |
CN102890642B (zh) | 基于匹配矩阵的异构重构计算的性能分析方法 | |
Ozen et al. | Genetic algorithm based NoC design with voltage/frequency islands | |
CN102253861A (zh) | 一种分步运算插件的执行方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130320 Termination date: 20201229 |