[论文解读] Simulations of detonation waves with smoothed dissipative particle dynamics

Gérôme Faure, Jean‐Bernard Maillet|arXiv (Cornell University)|Sep 11, 2017

Gas Dynamics and Kinetic Theory被引用 3

一句话总结

该论文通过为每个介观粒子引入一个随化学动力学演化的反应进度变量，将平滑耗散粒子动力学（SDPD）方法扩展至模拟化学反应体系。该方法成功模拟了硝基甲烷中的爆轰波与冲击波转爆轰过程，准确恢复了爆速和反应区宽度等关键物理特性，在分辨率足够时，结果对反应速率指前因子的变化具有鲁棒性。

ABSTRACT

Smoothed Dissipative Particle Dynamics (SDPD) is a mesoscopic method which allows to select the level of resolution at which a fluid is simulated. The aim of this work is to extend SDPD to chemically reactive systems.To this end, an additional progress variable is attached to each mesoparticle and evolves according to chemical kinetics. This reactive SDPD model is illustrated with numerical studies of the shock-to-detonation transition in nitromethane as well as the stationary behavior of the reactive wave.

研究动机与目标

• 将SDPD方法扩展至介观尺度下模拟化学反应体系。
• 实现对类似硝基甲烷等反应性材料中爆轰波及冲击波转爆轰过程（STD）的模拟。
• 评估反应速率指前因子对爆轰波结构与动力学的影响。
• 评估粒子分辨率对反应SDPD模拟精度的影响。
• 证明SDPD可在低于传统MD或DPDE方法所需分辨率下捕捉爆轰波的关键物理特性。

提出的方法

• 为每个介观粒子引入一个反应进度变量λ，以追踪反应体系中的化学转化程度。
• 将DPDE方法中的反应机理（Maillet et al., 2007）适配至SDPD框架。
• 采用阿 Arrhenius 速率定律，结合温度依赖的反应速率，控制λ的演化。
• 应用移动窗口技术，模拟爆轰波的长时间动力学行为。
• 采用热力学一致的SDPD公式，确保能量守恒与流体动力学力的正确性。
• 通过有限数量的介观粒子结合光滑粒子插值方法，求解Navier-Stokes方程。

实验结果

研究问题

• RQ1 SDPD能否有效扩展至介观尺度下基于化学动力学的反应体系建模？
• RQ2 Arrhenius指前因子如何影响硝基甲烷中爆轰波的结构与动力学？
• RQ3 粒子分辨率在多大程度上影响快速反应下反应SDPD模拟的精度？
• RQ4 反应SDPD模型能否正确再现冲击波转爆轰过程及稳定爆轰波的特性？
• RQ5 在分辨率足够时，该方法是否能在反应速率变化下实现尺度不变性？

主要发现

• 当分辨率足够时，所有测试指前因子下的爆速均与理论预测值相差小于1%。
• 反应区宽度近似与指前因子成反比，表明时间与长度尺度可实现重标度。
• 当指前因子为4×10¹⁵ s⁻¹时，若采用粗分辨率，爆速偏差达15%，表明对快速动力学而言分辨率不足。
• 在相同高指前因子下，采用更高分辨率模拟（K=10）得到的爆速为6777 m·s⁻¹，接近理论预期值。
• 经重标度后，不同指前因子下的压力分布曲线完全重合，证实系统的尺度不变性。
• 当介观粒子过大时，反应机理无法准确模拟快速反应，凸显需根据分辨率进行参数调优。