北京大学学报自然科学版 ›› 2024, Vol. 60 ›› Issue (1): 89-108.DOI: 10.13209/j.0479-8023.2023.022

上一篇    下一篇

高焓流动中的可压缩颗粒两相流并行求解器: 数值方法及其验证

李青1,2, 余钊圣3, 刘朋欣1, 李婷婷4, 陈坚强1, 袁先旭1,†   

  1. 1. 空天飞行空气动力科学与技术全国重点实验室, 中国空气动力研究与发展中心, 绵阳 621000 2. 天目山实验室, 杭州 310000 3. 浙江大学航空航天学院, 杭州 310027 4. 西安交通大学化学工程与技术学院, 西安 710049
  • 收稿日期:2022-12-12 修回日期:2023-01-29 出版日期:2024-01-20 发布日期:2024-01-20
  • 通讯作者: 袁先旭, E-mail: yuanxianxu(at)cardc.cn
  • 基金资助:
    国家自然科学基金(12272398)、中国博士后科学基金(2021MD703970)和国家重点研发计划(2019YFA0405200)资助

MPI Solver of Particle-Laden Compressible High Enthalpy Flow: Numerical Method and Validation

LI Qing1,2, YU Zhaosheng3, LIU Pengxin1, LI Tingting4, CHEN Jianqiang1, YUAN Xianxu1,†   

  1. 1. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000 2. Tianmushan Laboratory, Hangzhou 310000 3. School of Aeronautics and Astronautics, Zhejiang University, Hangzhou 310027 4. School of Chemical Engineering and Technology, Xi’an Jiaotong University, Xi’an 710049
  • Received:2022-12-12 Revised:2023-01-29 Online:2024-01-20 Published:2024-01-20
  • Contact: YUAN Xianxu, E-mail: yuanxianxu(at)cardc.cn

摘要:

在极端力学环境下可压缩点力颗粒两相流理论方程的基础上, 提出一种基于动态链表数组的并行颗粒求解器, 使其与可压缩流体求解器耦合。与基于欧拉坐标系的携带流体相不同, 基于拉格朗日坐标系的求解器采用动态链表数组对弥散颗粒相进行内存分配, 可以解决因使用全局数组解决拉格朗日/欧拉坐标系转换带来的弥散颗粒群内存利用率低和计算效率低下问题。最后对多物理效应的可压缩颗粒两相流求解器进行验证, 并在马赫数渐进趋于零的情况下, 对两个不可压缩槽道颗粒湍流标模进行验证。

关键词: 动态链表数组, 拉格朗日/欧拉坐标系转换, 可压缩点力颗粒两相流直接数值模拟求解器

Abstract:

Based on the theoretical equations of compressible particle-laden flow in multi-physical conditions, a MPI-Particle Solver based on dynamic linked list array is developed, which is coupled with a compressible flow solver of carried phase. The dispersed particle phase is different from the carried fluid phase, the previous one is based on the Lagrange coordinate system, while the later one is based on the Euler coordinate system. The traditional way is to use the global array to assign the memory to the dispersed particle phase to achieve the transfer between two-phase, but at expense of low computational efficiency. This paper uses dynamic linked list array to allocate memory to dispersed particle phase, achieve both problems. A series of DNS validations with references case have been done, in terms of multi-physical effects and two canonical cases at incompressible limit.

Key words: dynamic linked list array, Lagrange/Euler coordinate system conversion, particle solver of particle-laden compressible flow