[논문 리뷰] Computational study of shock waves propagating through air-plastic-water interfaces

Mauricio J. del Razo, Randall J. LeVeque|arXiv (Cornell University)|2015. 03. 31.

Combustion and Detonation Processes참고 문헌 11인용 수 4

한 줄 요약

이 연구는 충격파가 공기-플라스틱-물 인터페이스를 통과할 때를 모의하기 위해 1차원 유한부피법과 Tammann 상태방정식, 하이브리드 HLLC-정확한 리만 해법을 사용한다. 뇌 손상의 원인으로 간주되는 폭발파에 의한 뇌 손상(TBI)을 배경으로 하며, 주요 결과로는 얇은 플라스틱 인터페이스는 모의에서 무시해도 되며, 이는 충격파 전파에 미치는 영향이 무시할 만큼 작기 때문이다. 이는 TBI 연구에서 복잡한 생물학적 인터페이스를 단순화하는 데 유리하다.

ABSTRACT

The following study is motivated by experimental studies in traumatic brain injury (TBI). Recent research has demonstrated that low intensity non-impact blast wave exposure frequently leads to mild traumatic brain injury (mTBI); however, the mechanisms connecting the blast waves and the mTBI remain unclear. Collaborators at the Seattle VA Hospital are doing experiments to understand how blast waves can produce mTBI. In order to gain insight that is hard to obtain by experimental means, we have developed conservative finite volume methods for interface-shock wave interaction to simulate these experiments. A 1D model of their experimental setup has been implemented using Euler equations for compressible fluids. These equations are coupled with a Tammann equation of state (EOS) that allows us to model compressible gas along with almost incompressible fluids or elastic solids. A hybrid HLLC-exact Eulerian-Lagrangian Riemann solver for Tammann EOS with a jump in the parameters has been developed. The model has shown that if the plastic interface is very thin, it can be neglected. This result might be very helpful to model more complicated setups in higher dimensions.

연구 동기 및 목표

• 저강도 폭발파가 경미한 뇌 손상(mTBI)을 유발하는 메커니즘을 이해하기 위해, 기존 영상 검사로는 손상 여부를 확인할 수 없는 경우를 대비한다.
• 실험 데이터만으로는 복잡한 인터페이스에서의 동역학을 해석할 수 없기에, 1차원 모델을 사용해 공기-플라스틱-물 인터페이스를 통과하는 충격파 전파를 시뮬레이션한다.
• 실험 장치(예: Transwell 용기)에서 사용되는 플라스틱 인터페이스가 충격파 전파에 미치는 영향을 평가하고, 계산 모델에서 이를 단순화하거나 생략할 수 있는지 판단한다.
• 특히 공기에서 물로의 전이에서 압력 증폭과 파동 왜곡을 정량화하여 뇌 내피세포 손상 위험 평가에 기여한다.
• 폭발에 의한 mTBI 연구를 위한 향후 3차원 모의 및 두상유동 모델의 기초를 마련한다.

제안 방법

• 압축성 유체의 에이우러 방정식을 해결하기 위해 보존형 유한부피법을 사용하며, 압축성 기체와 거의 불압축성 유체/고체를 모의하기 위해 Tammann 상태방정식(EOS)을 적용한다.
• 특히 Tammann 상태방정식을 고려할 때 물질 특성의 불연속성과 점프를 처리하기 위해 하이브리드 HLLC-정확한 에이우러-라그랑주 리만 해법을 개발하였다.
• 충격관 실험 설정을 재현하기 위해 1차원 모델을 구성하였으며, 공기, 얇은 플라스틱 층, 물로 구성된 시스템에서 충격파가 인터페이스에 상호작용하는 영역을 중심으로 분석하였다.
• 플라스틱 인터페이스 두께를 다양하게 변화시켜 전파된 충격파의 진폭과 형상에 미치는 영향을 평가하였다.
• 플라스틱 두께가 감소함에 따라 분석적 해에 수렴하는 것으로 나타나, 선형 음향학 이론과의 비교를 통해 모델의 타당성을 검증하였다.
• 전체 공기-플라스틱-물 모델과 비교하여, 단순화된 공기-물 인터페이스 모델을 테스트하여 플라스틱 층 생략의 정확도를 평가하였다.

실험 결과

연구 질문

• RQ1 플라스틱 인터페이스 두께가 공기-플라스틱-물 인터페이스를 통과하는 충격파 전파에 미치는 영향은 어떠한가?
• RQ2 직접적인 공기-물 인터페이스와 비교할 때 플라스틱 층은 전파된 충격파의 진폭과 형상에 얼마나 큰 영향을 미치는가?
• RQ3 특히 물 내 충격파 증폭을 고려할 때, 계산 모델에서 플라스틱 인터페이스를 생략해도 정확도에 큰 영향을 주지 않는가?
• RQ4 공기에서 물로 전이될 때 압력 증폭의 크기는 얼마이며, 이는 뇌 조직 손상 메커니즘과 어떻게 관련되는가?
• RQ5 얇은 플라스틱 층 내에서의 다중 반사 및 간섭 무늬는 하류의 시간에 따른 압력 프로파일에 어떤 영향을 미치는가?

주요 결과

• 공기-물 인터페이스의 경우, 감도계 4(물 내)에서 최대 압력 진폭이 284.26 KPa에 도달하였으며, 이는 초기 공기 충격파 184.06 KPa 대비 54% 증가한 것이다.
• 플라스틱 인터페이스 두께를 0.1 m로 줄였을 때, 물 내 전파된 충격파 프로파일은 공기-물 인터페이스의 경우와 거의 구분되지 않으며, 최대 진폭 오차는 0.38%에 불과하다.
• 플라스틱 두께가 감소함에 따라 플라스틱 층 내 다중 반사 간섭이 매우 빠르게 발생하여 압력 프로파일이 단일 충격파 형태로 수렴한다.
• 가장 얇은 플라스틱 케이스(0.1 m)와 플라스틱 없음 케이스(공기-물)의 압력 감도계 프로파일은 눈에 띄게 구분되지 않으며, 플라스틱 층의 영향이 거의 없음을 시사한다.
• 모의 결과는 선형 음향학 이론의 예측과 일치한다: 얇은 플라스틱 층을 통과한 전파 압력은 직접적인 공기-물 인터페이스의 예상 값에 수렴한다.
• 모델링 목적에서 공기-플라스틱-물 인터페이스는 매우 얇은 플라스틱 층을 고려할 경우, 공기-물 인터페이스로 단순화해도 정확도에 큰 영향을 주지 않음을 입증하였다.