斜爆轰波并行数值模拟的化学加速算法性能研究

doi:10.13675/j.cnki.tjjs.200724

推进技术 ›› 2021, Vol. 42 ›› Issue (4): 765-775.DOI: 10.13675/j.cnki.tjjs.200724

斜爆轰波并行数值模拟的化学加速算法性能研究

姜长磊¹，董刚¹，吴锦涛²

1.南京理工大学瞬态物理重点实验室，江苏南京 210094;2.北京宇航系统工程研究所，北京 100076

出版日期:2021-04-15 发布日期:2021-04-15
作者简介:姜长磊，硕士生，研究领域为可压缩反应流的数值模拟。E-mail：137868777@qq.com
基金资助:
航天一院CALT基金项目（CALT2020-ZD01）；南京理工大学瞬态物理重点实验室开放基金（6142604180205）。

Performance Study of Chemical Acceleration Algorithm for Parallel Numerical Simulation of Oblique Detonation Wave

1.Key Laboratory of Transient Physics，Nanjing University of Science and Technology，Nanjing 210094，China;2.Beijing Institute of Astronautical System Engineering，Beijing 100076，China

Online:2021-04-15 Published:2021-04-15

摘要/Abstract

摘要： 考虑基元反应条件下的爆轰波精细结构的数值模拟计算量巨大，发展高精度和高效率的计算方法十分必要。以Ma=7的H₂/O₂/N₂预混气来流形成的斜爆轰波为数值模拟对象，研究了一种基于并行计算架构的用于加速化学反应计算的建表算法的计算性能，考察了不同建表策略，即TP（Transposed processing）策略和PLP（Purely local processing）策略，以及不同并行分区数量对算法性能的影响。研究结果显示，本文采用的建表算法能够很好地再现斜爆轰结构，其计算精度不受建表策略和并行计算分区数量的影响；而算法的计算效率则取决于不同分区对应的数据表之间操作的同步性，其中，数据表中节点数据的取回率和数据表设定的尺寸上限都会影响数据表操作的同步性。本文采用的两种建表策略的计算结果表明，TP策略数据表中节点取回率高于PLP策略，故计算效率更高；而计算分区数量越少，则分区对应的数据表尺寸上限越大，数据表的同步性就越好，计算效率也越高。

关键词: 斜爆轰波;建表算法;并行计算;计算性能;数值模拟

Abstract: The numerical simulations of fine detonation wave structures with detailed reaction mechanism are computationally expensive， so it is necessary to develop an algorithm with high-accuracy and high-efficiency. In present study， the numerical simulations of oblique detonation wave generated by a H₂/O₂/N₂ incoming flow with Ma=7 were used to examine computational performance of a tabulation algorithm， which is used to accelerate the computations of chemical reactions on the parallel frame. The effects of the tabulation strategy of algorithm， that is， TP （Transposed processing） and PLP （Puely local processing） strategies， and the parallel zone number on the algorithm performance were focused on. The results show that algorithm in present study can well reproduce the structures of oblique detonation wave， and that the computational accuracy is independent of tabulation strategy and parallel zone number of the algorithm. The computational efficiency of the algorithm is determined by the operation synchronization among the data tables that are mapped into the corresponding computational zones. Both the retrieval ratio of nodes in the tables and the upper limit of table size have the influence on the operation synchronization among the tables. For TP and PLP strategies of the algorithm in present study， the retrieval ratio of nodes in TP strategy is larger than that in PLP strategy and therefore gives rise to the higher computational efficiency. The simulation for case of fewer zone number which has larger upper limit of table size leads to the better synchronization and subsequently the higher computational efficiency.

Key words: Oblique detonation wave;Tabulation algorithm;Parallel computation;Computational performance;Numerical simulation

姜长磊，董刚，吴锦涛. 斜爆轰波并行数值模拟的化学加速算法性能研究[J]. 推进技术, 2021, 42(4): 765-775.

JIANG Chang-lei¹, DONG Gang¹, WU Jin-tao². Performance Study of Chemical Acceleration Algorithm for Parallel Numerical Simulation of Oblique Detonation Wave[J]. Journal of Propulsion Technology, 2021, 42(4): 765-775.

参考文献

[1] Wolański P. Detonative Propulsion［J］. Proceedings of the Combustion Institute， 2013， 34： 125-158.
[2] 范宝春. 极度燃烧［M］. 北京：国防工业出版社， 2018.
[3] Choi J Y， Kim D W， Jeung I S， et al. Cell-like Structure of Unstable Oblique Detonation Wave from High-Resolution Numerical Simulation［J］. Proceedings of the Combustion Institute， 2007， 31： 2473-2480.
[4] Teng H H， Jiang Z L， Ng H D. Numerical Study on Unstable Surfaces of Oblique Detonations［J］. Journal of Fluid Mechanics， 2014， 744： 111-128.
[5] Teng H， Ng H D， Li K， et al. Evolution of Cellular Structures on Oblique Detonation Surfaces［J］. Combustion and Flame， 2015， 162： 470-477.
[6] Teng H， Ng H D， Jiang Z. Initiation Characteristics of Wedge-Induced Oblique Detonation Waves in a Stoichiometric Hydrogen-Air Mixture［J］. Proceedings of the Combustion Institute， 2017， 36： 2735-2742.
[7] Fang Y， Zhang Y， Deng X， et al. Structure of Wedge-Induced Oblique Detonation in Acetylene-Oxygen-Argon Mixtures［J］. Physics of Fluids， 2019， 31（2）.
[8] 陈楠， Esfehani S A， Bhattrai S，等. 当量比对斜爆轰波诱导区特性影响的数值模拟研究［J］. 推进技术， 2018， 39（12）： 2798-2805.
[9] Pope S B. Computationally Efficient Implementation of Combustion Chemistry Using in Situ Adaptive Tabulation［J］. Combustion Theory and Modelling， 1997， 1（1）： 41-63.
[10] Dong G， Fan B C， Chen Y L. Acceleration of Chemistry Computations in Two-Dimensional Detonation Induced by Shock Focusing Using Reduced ISAT［J］. Combustion Theory and Modeling， 2007， 11（5）： 823-837.
[11] Dong G， Fan B C. Chemistry Acceleration Modeling of Detonation Based on the Dynamical Storage/Deletion Algorithm［J］. Combustion Science and Technology， 2009， 181（9）： 1207-1216.
[12] Lu L， Lantz S R， Ren Z， et al. Computationally Efficient Implementation of Combustion Chemistry in Parallel PDF Calculations［J］. Journal of Computational Physics， 2009， 228（15）： 5490-5525.
[13] Wu Jintao， Dong Gang， Li Baoming. Parallel Chemistry Acceleration Algorithms Based on ISAT Method in Gaseous Detonation Computations［J］. Computers and Fluids， 2018， 167： 256-284.
[14] Wu Jintao， Dong Gang， Li Yi. Parallel Chemistry Acceleration Algorithm with ISAT Table Size Control in the Application of Gaseous Detonation［J］. Shock Waves， 2019， 29： 523-535.
[15] Jiang G S， Shu C W. Efficient Implementation of Weighted ENO Schemes［J］. Journal of Computational Physics， 1995， 126（1）： 202-228.
[16] Balsara D S， Shu C W. Monotonicity Preserving Weighted Essentially Non-Oscillatory Schemes with Increasingly High Order of Accuracy［J］. Journal of Computational Physics， 2000， 160（2）： 405-452.
[17] Teng H， Zhang Y， Jiang Z. Numerical Investigation on the Induction Zone Structure of the Oblique Detonation Waves［J］. Computers and Fluids， 2014， 95（3）： 127-131.
[18] Burke M P， Chaos M， Ju Y， et al. Comprehensive H₂/O₂ Kinetic Model for High-Pressure Combustion［J］. International Journal of Chemical Kinetics， 2012， 44（7）： 444-474.
[19] Brown P N， Byrne G D， Hindmarsh A C. VODE： A Variable-Coefficient ODE Solver［J］. SIAM Journal on Scientific and Statistical Computing， 1989， 10（5）： 1038-1051.
[20] 董刚，范宝春，朱旻明，等. 动态存储方法在气相爆轰波数值模拟中的应用［J］. 爆炸与冲击， 2008， 28（1）.

斜爆轰波并行数值模拟的化学加速算法性能研究

Performance Study of Chemical Acceleration Algorithm for Parallel Numerical Simulation of Oblique Detonation Wave

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 0

编辑推荐

Metrics

本文评价