[논문 리뷰] Molecular Dynamics Simulation on Stability of Converging Shocks

В. В. Жаховский, K. Nishihara외 1명|ArXiv.org|2009. 04. 22.

Laser-Plasma Interactions and Diagnostics인용 수 3

한 줄 요약

이 연구는 밀도가 높은 레너드존 유체에서 수렴하는 원통형 충격파에 대한 최초의 분자역학(MD) 시뮬레이션을 제시하며, MD가 원자 척도에서 유체역학적 불안정성을 높은 정밀도로 포착할 수 있음을 보여준다. 이 방법은 머스크럼 형성, 초음속 제트, 모드에 따라 달라지는 리히트마이어-메쉬코프 불안정성 성장 등을 밝혀내며, 후기 단계의 역행으로 인해 늦게 발생하는 모드 3의 리파일은 일반적인 유체역학 예측에 도전한다.

ABSTRACT

Molecular Dynamic (MD) approach is applied to study the converging cylindrical shock waves in a dense Lennard-Jones (LJ) fluid. MD method is based on tracking of the atom motions and hence it has an fundamental advantages over hydrodynamic methods which assumes shocks as a structureless discontinuity and requires an equation of state. Due to the small thickness of shock fronts in liquid the two million particles is enough to simulate propagation of a cylindrical shocks in close detail. We investigate stability of converging shocks with different perturbation modes and its mixture. It was shown that in a case of relatively large initial ripples the Mach stems are formed. Supersonic jets generated by interaction of reflected shocks in downstream flow are observed. We also study the Richtmyer-Meshkov (RM) instability of an interface between two Lennard-Jones liquids of different mass densities. Surprisingly, mode 3 ripples grow very slow in comparison with higher mode numbers and growth rate of a higher mode decay slower.

연구 동기 및 목표

• 밀도가 높은 유체에서 수렴하는 원통형 충격파의 안정성을 원자 척도 분자역학 시뮬레이션을 통해 조사하기 위해.
• 충격파를 불연속성으로 간주하고 상태방정식에 의존하는 유체역학적 방법의 한계를 극복하기 위해.
• 특히 다양한 모드 수를 가진 표면 불안정성의 성장 역학을 고려하여, 원자 척도에서 리히트마이어-메쉬코프 불안정성을 검토하기 위해.
• 충격 압축 및 반사 과정에서 초기 불안정성 모드와 그 비선형 진화의 역할을 평가하기 위해.
• 원자 해상도로 극한 조건(최대 50 GPa, 12,000 K)을 시뮬레이션하여 관성형 핵융합 과학의 기준 자료를 제공하기 위해.

제안 방법

• 시뮬레이션은 NEC SX-5 슈퍼컴퓨터에서 실행된 완전한 벡터화된 MD 코드를 사용하며, 축 방향으로는 주기적 경계 조건을 갖는 원통형 영역 내에서 2,097,148개의 레너드존 원자를 처리한다.
• 레너드존 포텐셜은 r_cut = 2.5σ에서 잘라내고, 물리적으로 비합리적인 힘을 방지하기 위해 스무딩 항목을 추가하여 수정하였으며, 매개변수는 아르곤에 캘리브레이션된 상태(σ = 3.405 Å, ε/kB = 119.8 K)를 사용한다.
• 각도에 따라 위치가 변하는 피스톤 경계 R(φ) = R₀[1 + δ sin(mφ)]를 통해 초기 불안정성 모드 수 m과 상대적 진폭 δ = 0.05를 갖는 불안정성을 생성한다.
• 계산은 T = 0.72 mdu (≈83 K) 및 n = 0.79 mdu⁻³로 초기화되며, 두 가지 다른 물질(경량 m_a = 48, 중량 m_b = 16m_a)을 모의하기 위해 표면에서 질량과 속도를 급격히 변화시킨다.
• 준비 단계 동안 랭지에인 열역학기구가 열적 평형을 유지하며, 충격파 발생은 t = 0에서 피스톤을 내측으로 이동시켜 시작된다.
• 시간에 따른 불안정성 진폭과 혼합 영역 두께를 추적하여 성장률을 계산하며, 비선형 진화는 표면 위치 및 Δ(t)의 시계열 분석을 통해 분석된다.

실험 결과

연구 질문

• RQ1 다양한 초기 불안정성 모드(m)가 밀도가 높은 유체에서 수렴하는 원통형 충격파의 안정성과 진화에 어떤 영향을 미치는가?
• RQ2 충격파 수렴 과정에서 어떤 원자 척도의 구조—예를 들어 머스크럼 및 반사된 충격파 패tern—이 나타나는가?
• RQ3 두 물질로 구성된 레너드존 시스템에서 리히트마이어-메쉬코프 불안정성은 어떻게 진화하는가? 성장률은 모드 수에 따라 어떻게 달라지는가?
• RQ4 왜 모드 3는 더 높은 모드들에 비해 비정상적으로 느리게 성장하는가? 그리고 위상 역행은 어떤 역할을 하는가?
• RQ5 특히 극한의 압축 조건에서, 분자역학이 기존의 유체역학 코드보다 유체역학적 불안정성을 더 높은 정밀도로 해석할 수 있는가?

주요 결과

• 초기 불안정성 모드 수 m ≥ 5일 경우 충격파 전면의 불안정성 진폭이著しく 증가하지만, m < 5일 경우 성장이 억제되어 모드에 따라 달라지는 안정성 임계값이 존재함을 시사한다.
• 초기 리파일이 크면 충격파 수렴의 매우 초기 단계에서 머스크럼이 형성되며, 위상 역행 기간 동안 m 또는 2m개의 변을 가진 다각형 충격파 전면이 나타난다.
• 후행 유동에서 반사된 충격파 간의 상호작용으로 인해 초음속 제트가 생성되며, 이는 관성형 핵융합 타깃의 박리면 불안정성의 증폭 가능성을 암시한다.
• 혼합 영역 두께 Δ(t)는 m = 8일 경우 Δ₈ ∼ t^0.45, m = 5일 경우 Δ₅ ∼ t^0.3, m = 3일 경우 Δ₃ ∼ ln(t)로 증가하여 시간이 지남에 따라 성장률이 감소함을 나타낸다.
• 모드 3의 성장률은 t^−1로 감소하며, 이는 더 높은 모드들보다 훨씬 느리게 감소함을 보이며, 이는 매우 늦은 시점에 위상 역행이 발생하기 때문이다.
• 수렴의 최점에서 유체는 약 2.2배의 압축, 약 12,000 K의 온도, 약 50 GPa의 압력을 경험하며, 충격파 수렴은 약 20 ps 내에 발생한다.