面向领域应用的空间域和频域分解模式并行计算

发布时间：2018-04-24 06:30

本文选题：并行计算 + 区域分解　；参考：《吉林大学》2016年博士论文

【摘要】：数值计算在现代科学研究中发挥着重要作用,是帮助研究人员分析越来越复杂过程模型的重要手段。并行计算技术可以用来解决科学计算应用中耗时长和模拟模型规模受阻问题。这方面的研究包括基础计算平台和并行应用两个方面,并行应用是结合具体领域问题发挥并行计算效用的落脚点。开发一个并行数值计算软件,需要将原有串行计算模式转化为并行化计算模式。区域分解是最常见的并行化方法,区域在实际问题中分为空间域和频率域等。设计并实现相应的区域分解方法,是开发并行计算模式的基础。本文关注这两种并行模式的具体设计和并行计算的实际应用。分别面向有限体积法地下复杂环境耦合过程和电磁场谐波平衡有限元计算问题,设计了相应的并行计算方法并完成了实际领域的并行计算应用。本文在地下复杂环境并行数值计算的研究工作受到中国地质调查局计划项目“全国二氧化碳地质储存潜力评价与示范工程”的支持。以美国劳伦斯伯克利实验室(LBNL)TOUGHREACT三场耦合模拟模型为基础,与吉林大学环境与资源学院课题组合作开发分布式内存平台上的高性能计算软件THC-MP。该软件同LBNL内部开发的基于共享内存并行模式的TOUGHREACT-OMP相比,在并行扩展度和计算效率方面有很多优势,在公开文献中首次实现了分布式内存平台上的耦合过程并行计算。该软件成功实现了鄂尔多斯盆地储层大型三维模型的模拟计算,为二氧化碳地质封存示范性项目数值模拟工作提供了支持。本文在电磁场分析计算的研究工作受到博士研究生联合培养项目支持,与澳大利亚格里菲斯大学工学院合作解决了谐波平衡有限元方法的高效计算问题。通过在频率域上对多次谐波进行分解,设计了频域并行求解方法,并将该算法应用到变压器直流偏磁问题计算。为稳态电磁场复合频率问题和工业产品设计的高效模拟计算提供了一种快速解法。本文工作的主要贡献有:1.针对环境科学领域多相流体、溶质运移和化学反应耦合模拟过程,基于空间区域分解方法,完整实现了整个耦合过程的并行计算。将流体、运移、化学三部分的数据按照子空间网格集合进行重组,设计了基于统一空间区域分解的耦合过程并行计算框架。在核心方程求解阶段实现了矩阵运算和网格处理计算的混合并行模式,优化了耦合过程并行计算的负载均衡。2.将耦合过程的并行计算框架实用化,开发了热-流体-化学三场耦合的模拟器THC-MP。和二氧化碳封存模拟计算模块ECO2N结合,在松辽盆地储层模型上验证了并行计算的数值可靠性。并将该模拟器应用到了鄂尔多斯盆地石千峰组储层的二氧化碳地质封存过程模拟,在分布式内存并行平台上完成了场地大规模三维模型的高效并行模拟。3.对电磁场计算中的谐波平衡有限元分析问题,在频率域上对数值模拟计算进行分解,设计了基于谐波分解的块矩阵类Jacobi迭代并行求解方法。并针对一个含有多次谐波的高频变压器直流偏磁问题,使用等同于谐波次数并行度的计算对该方法进行验证。实验结果显示,频率域并行计算方法减小了含有多次谐波的模型矩阵规模同时有效缩短了模拟计算时间。
[Abstract]:Numerical calculation plays an important role in modern scientific research. It is an important means to help researchers to analyze more and more complex process models. Parallel computing technology can be used to solve the problem of time consuming and the scale of simulation model in the application of scientific computing. This research includes two aspects of basic computing platform and parallel application, Parallel applications are the foothold of playing the utility of parallel computing in combination with specific domain problems. Developing a parallel numerical computing software needs to transform the original serial computing model into a parallel computing model. Regional decomposition is the most common parallelization method, and the region is divided into space domain and frequency domain in practical problems. The method of regional decomposition is the basis of developing parallel computing model. This paper focuses on the concrete design of the two parallel modes and the practical application of parallel computing. This paper has designed the parallel computing method of the finite volume method for the coupled process of underground complex environment and the finite element calculation of electromagnetic harmonic balance, and completed the practical field. The research work on parallel numerical calculation in the underground complex environment is supported by the China Geological Survey Project "national carbon dioxide geological storage potential evaluation and demonstration project", based on the three field coupling simulation model of the Laurence Berkeley Laboratory (LBNL) of the United States and the Jilin University ring. A high performance computing software THC-MP. on distributed memory platform is developed by the combination of frontier and resource colleges. Compared with TOUGHREACT-OMP based on shared memory parallel mode developed in LBNL, the software has many advantages in parallel scalability and computational efficiency. In the public literature, the first implementation of the distributed memory platform coupling is for the first time. The software has successfully implemented the simulation calculation of the large three-dimensional model of the Ordos Basin, which provides support for the numerical simulation of the carbon dioxide geological sequestration demonstration project. This paper is supported by the joint training program of the doctoral graduate students in the electromagnetic field analysis and calculation, with Griffith, Australia. The University Institute of technology has worked together to solve the problem of efficient calculation of the harmonic balance finite element method. By decomposing the multiple harmonics on the frequency domain, a parallel solution method in frequency domain is designed. The algorithm is applied to the calculation of the DC bias of the transformer. The high efficiency simulation calculation for the complex frequency problem of the steady-state electromagnetic field and the design of the industrial product is carried out. A fast solution is provided. The main contributions of this work are as follows: 1. the simulation process of multiphase fluid, solute transport and chemical reaction in the field of environmental science, and the complete realization of the parallel computation of the whole coupling process based on the spatial domain decomposition method. The data of the three parts of the fluid, transport and chemical are set up according to the subspace grid. The parallel computing framework of coupling process based on the unified spatial domain decomposition is designed. The hybrid parallel mode of matrix operation and grid processing is realized in the phase of core equation solving. The parallel calculation framework of the parallel computation of the coupling process is optimized. The parallel calculation framework of the coupling process is practiced, and the thermal fluid chemical three is developed. The field coupled simulator THC-MP. and the carbon dioxide sequestration simulation module ECO2N are combined to verify the numerical reliability of the parallel computation on the Songliao basin reservoir model, and the simulator is applied to the simulation of the carbon dioxide geological sequestration process of the stone thousand Peak Reservoir in the Ordos Basin and completed on the distributed memory parallel platform. The high efficiency parallel simulation of the large scale three-dimensional model of the field.3. is used to solve the problem of harmonic balance finite element analysis in electromagnetic field calculation. The numerical simulation calculation is decomposed on the frequency domain. A block matrix class Jacobi iterative parallel solution method based on harmonic decomposition is designed. The method is verified by the calculation of the parallel degree of harmonic frequency. The experimental results show that the parallel computation method of frequency domain reduces the scale of the model matrix with multiple harmonics and reduces the simulation time effectively.

【学位授予单位】：吉林大学
【学位级别】：博士
【学位授予年份】：2016
【分类号】：TP338.6

【相似文献】