二维柱几何中子输运方程的并行区域分解方法

分析不同的区域分解方法及优先级插入算法对二维柱几何下中子输运方程S_n间断有限元方程并行效率的影响,给出基于最小面体比的正方形区域分解方法及沿径向的优先级插入算法,并通过将正方形区域分解方法与径向优先级插入算法进行组合,形成新的算法.新算法更适应于二维柱几何下输运方程S_n间断有限元方法的并行计算.数值试验表明,在通信延迟较高的大型国产并行机上,新算法用数百个CPU还可以取得较好的并行效果,比已有方法具有更良好的可扩展性.     

高次方程异步并行迭代算法

本文给出一个基于共享存储MIMD计算机模型的求解任意实系数高次方程全部实数根的异步并行迭代算法，并分析算法的复杂度。     

中子输运方程源项反演的瑞利商迭代算法

中子输运方程源项反演问题中,有效倍增因子的求解经过球谐函数展开和差分离散后,转化为求解大型矩阵的特征值问题.针对矩阵的特点,利用Gerschgorin圆盘定理,给出反幂法迭代初值的选取方法,并进行相应的扰动分析.针对有效倍增因子是矩阵最大特征值的特点,采用瑞利商迭代算法进行求解.理论和算例表明,该方法不需选取特定初值,达到三阶收敛速度.     

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

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

中子输运方程基于界面修正的源迭代并行计算方法

中子输运方程的计算量非常大。在现有的计算机条件下，进行精密物理的数值模拟所需要的中子计算仍是非常的费时间和费内存的，不采用并行计算是难以承受的。并且，由于中子输运隐式离散纵标方法引起的数据强相关，以及计算过程必须严格沿中子运动方向进行(否则会出现计算不稳定)，因此会出现相当严重的算法同步的问题，使得隐式格式在大型并行计算机上实施时所能得到的并行度十分有限，严格限制了其实现具有高并行度的迭代计算的可能性。因此，对中子输运方程隐式差分格式进行并行改造是十分必要的。     

扩散方程区域分解的多步算法

利用分数步法进行内边界值的多步计算,改进二维扩散方程的区域分解算法,形成新的并行算法,放宽稳定性条件.其中采用分数步空间大步长离散格式计算内边界点值.算法精度与隐格式相当.与改进前相比,稳定性条件放宽了q倍(g为两个相邻时间步之间执行分数步内边界值计算的次数).利用离散极值原理,严格证明了算法的收敛性.在并行机上进行数...     

波动方程基于自然边界归化的区域分解算法

提出了无界区域波动方程的区域分解算法,基于自然边界归化,分别研究了重叠型与非重叠型区域分解算法,首先将控制方程对时间进行离散化,得到关于时间步长离散化格式,对每一时间步长给出了Dirichlet-Neumann和Schwartz交替算法,对Schwartz交替算法,给出了算法的收敛性,对圆外区域研究了压缩因子,并给出了数值例子。     

高次方程的一个分布式并行迭代算法

本文给出一个基于分布式存储MIMD树机模型求解实系数高次方程全部实数根的并行迭代算法并分析算法的复杂度。     

二维柱几何非定常中子输运方程基于格式的界面预估校正并行算法

针对二维柱几何非定常中子输运方程的Sn-间断有限元方法,提出基于格式的界面预估校正并行算法.数值算例表明,该并行算法在精度与并行度等诸方面均具有良好的性质,与已有的基于隐式格式的并行扫描算法相比,对于二维中子输运大规模计算问题,并行计算效率较高,并行加速比可增加-倍以上,且可保持原隐式格式的计算精度.     

粒子输运离散纵标方程基于界面修正的并行计算方法

为了改造粒子输运方程求解的隐式格式,研究设计适应大型并行计算机的并行计算方法,介绍一类求解粒子输运方程离散纵标方程组的基于界面修正的源迭代并行计算方法.应用空间区域分解,在子区域内界面处首先采用迎风显式差分格式进行预估,构造子区域的入射边界条件,然后,在各个子区域内部进行源迭代求解隐式离散纵标方程组.在源迭代过程中,在内界面入射边界处采用隐式格式进行界面修正.数值算例表明该并行计算方法在精度、并行度、简单性诸方面均具有良好的性质.     

块三对角线性方程组的一类二维区域分解并行不完全分解预条件

基于二维重叠区域分解,对每个子区域上局部不完全LU分解所得到的上、下三角因子分别进行组合,给出一类全局并行不完全分解型预条件.所给出的并行化方法适用于任何不完全LU分解型预条件.对采用二维区域分解与一维区域分解时所得并行预条件的并行计算性能进行分析比较.实验结果表明,提出的并行化方法普遍优于加性Schwarz并行化方法,且当处理器个数相对较多时采用二维区域分解优于一维区域分解.     

非定常Navier-Stokes方程基于完全重叠型区域分解的有限元并行算法

基于完全重叠型区域分解技巧,提出三种求解非定常Navier-Stokes方程的有限元并行算法.其基本思想是首先对空间施行完全重叠区域分解,然后各个处理器使用向后Euler格式独立并行求解关于时间t的常微分方程;对于非线性的对流项,分别采用半隐格式和全隐格式进行处理.算法中每个处理器所负责的子问题是一个全局问题,它定义在整个求解区域上,但绝大部分自由度来自其所负责的子区域,从而使得算法实现简单,通信需求少.数值算例验证了算法的有效性及其良好的并行性能.     

三维粗糙面电磁双站散射的直接型区域分解计算

提出三维粗糙面双站电磁散射的直接型有限元-区域分解方法.首先建立含有迭代Robin边界条件(IRBC)的区域分解法耦合模型,再用内视法导出高度稀疏分块的分区耦合矩阵,之后给出缩减耦合矩阵带宽的子区域排序方法和IRBC的FFT加速算法.用有限元-完全匹配层和未分区的有限元-IRBC方法验证数值结果.     

抛物方程显隐修正迎风区域分裂差分法

给出抛物方程一种有效的区域分裂差分格式,提高了计算效率.对一阶项采用二阶迎风差分格式,内边界点和各子区域分别采用显隐差分格式.在较弱的稳定性条件下,得到离散l2模误差估计结果.最后给出具体的数值算例,以验证方法的实用性.     

二维抛物型问题的Strang型交替分段区域分裂格式

给出求解二维抛物型方程的Strang型的交替分段区域分裂格式。交替分段思想可以将区域分为一些不重叠的子区域,Strang型算子分裂技巧通过将高维问题的求解分解为几个低维问题的求解来降低其求解的复杂度。方法是无条件稳定的,理论分析了截断误差。数值算例说明格式的有效性及时空的二阶精度.