[논문 리뷰] Simulations of detonation waves with smoothed dissipative particle dynamics

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

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 시뮬레이션 정확도에 미치는 영향 평가하기.
• 전통적인 MD나 DPDE에 비해 더 굵은 해상도로도 폭발파의 핵심 물리적 성질을 잘 포착할 수 있음을 보여주기.

제안 방법

• 메조입자당 진행 변수 λ를 도입하여 반응성 시스템에서 화학 전환을 추적한다.
• Maillet 등(2007)의 DPDE에서 제안한 반응 메커니즘을 SDPD 프레임워크에 적응시킨다.
• 온도에 따라 변화하는 반응 속도를 가지는 아레니우스 속도 법칙을 사용하여 λ의 진화를 제어한다.
• 폭발파 장시간 동역학을 시뮬레이션하기 위해 이동 창 기법을 적용한다.
• 에너지 보존 및 유체역학적 힘을 고려한 열역학적으로 일관된 SDPD 수식을 구현한다.
• 유한한 수의 메조입자를 사용하고 스무딩 입자 보간법을 적용하여 나비에-스토크스 방정식을 해결한다.

실험 결과

연구 질문

• RQ1 SDPD는 메조스케일에서 화학 동역학을 고려한 반응성 시스템을 효과적으로 확장하여 모델링할 수 있는가?
• RQ2 아레니우스 계수는 메탄올의 폭발파의 구조 및 동역학에 어떤 영향을 미치는가?
• RQ3 메소입자 해상도가 빠른 반응에 대한 반응성 SDPD 시뮬레이션 정확도에 어느 정도의 영향을 미치는가?
• RQ4 반응성 SDPD 모델은 충격파에서 폭발파로의 전이 및 정적 폭발파 성질을 정확히 재현하는가?
• RQ5 해상도가 충분할 경우 반응 속도 변화에 관계없이 척도 불변성(스케일 인variant) 행동을 달성할 수 있는가?

주요 결과

• 해상도가 충분할 경우 모든 테스트된 계수에서 폭발 속도가 이론적 예측값과 1% 이내로 유지된다.
• 반응 영역 너비는 계수의 역수 비례로 증가하며, 시간 및 길이 척도의 재스케일링을 나타낸다.
• 계수 4×10¹⁵ s⁻¹일 경우 굵은 해상도를 사용할 경우 폭발 속도가 15% 이상 틀어지며, 이는 빠른 동역학에 대해 해상도가 부족함을 시사한다.
• 동일한 고계수를 사용하되 더 높은 해상도 시뮬레이션(K=10)에서는 폭발 속도가 6777 m·s⁻¹로 이론적 기대치에 가깝게 도출된다.
• 다양한 계수에 대한 압력 프로파일은 스케일링 후 동일한 곡선으로 수렴하며, 시스템의 척도 불변성을 확인한다.
• 메조입자가 너무 클 경우 반응 메커니즘이 빠른 반응을 정확히 모의하지 못하며, 이는 해상도에 따라 조정이 필요함을 시사한다.