[논문 리뷰] Poiseuille flow of liquid methane in nanoscopic graphite channels by molecular dynamics simulation

Martin Horsch, Jadran Vrabec외 2명|ArXiv.org|2009. 11. 29.

NMR spectroscopy and applications참고 문헌 9인용 수 3

한 줄 요약

이 연구는 대규모 분자 동역학 시뮬레이션을 통해 LS1/Mardyn 코드를 사용하여 나노스코픽한 극성 탄소 채널 내 액체 메탄의 포아줄 유동을 연구하였으며, 나노스케일의 경계 지배 유동에서 연속체 행동으로의 전이 기간 동안 Darcy의 법칙이 유지됨을 입증하였다. 시뮬레이션 결과, 20에서 50 nm 사이의 채널 폭에서 슬립 지배 유동과 연속체 유역 사이에 급격한 전이가 발생하는 것으로 나타났으며, 초고성능 컴퓨터를 통해 최대 480만 개의 상호작용 중심을 포함한 시스템에서 정확한 스케일링이 달성되었다.

ABSTRACT

MD simulations of methane confined between graphite walls with up to 4,800,000 interaction sites, i.e., carbon atoms and methane molecules, are conducted, where the channel width is varied to include both the boundary-dominated regime and the transition to the continuum regime. This proves that MD can be used today to cover the entire range of characteristic lengths for which continuum methods fail.

연구 동기 및 목표

• 분자 동역학를 사용하여 나노스코픽한 극성 탄소 채널 내 액체 메탄의 유체역학적 거동을 연구하기 위해.
• 대규모 유체-고체 시스템에 대한 LS1/Mardyn MD 시뮬레이터의 확장성을 평가하기 위해.
• 제한된 기하구조 내에서 나노스케일의 슬립 지배 유동과 연속체 유사 유동 사이의 전이 영역을 규명하기 위해.
• 탄소 나노소재에 봉인된 메탄에서 나노스케일에서 연속체로의 전이 기간 동안 Darcy의 법칙의 적용 가능성을 검증하기 위해.
• 최적화된 병렬 계산 기법을 사용하여 탄소 나노소재 내 유체 이동의 대규모 MD 시뮬레이션을 가능하게 하기 위해.

제안 방법

• 메탄-메탄 및 메탄-탄소 상호작용을 모델링하기 위해 잘라내고 이동한 레나르드존(비례) 잠재 에너지(LJ/TS)를 사용하였으며, 메탄의 열역학적 성질을 검증한 파arameter를 사용하였다.
• 메탄과 탄소 기반 탄소 원자 간의 이질 상호작용 파라미터를 유도하기 위해 Lorentz-Berthelot 혼합 규칙를 사용하였다.
• 실험 값(1.421 Å 대비 1.461 Å)과의 결합 길이 편차를 보정하기 위해, 탄소 기반 벽면에 재조정된 Tersoff 잠재 에너지를 적용하였다.
• 비균일한 시스템에서 부하 균형 및 통신 효율성을 최적화하기 위해, LS1/Mardyn 내에서 등방성 및 채널 기하학 기반의 공간 분할을 구현하였다.
• 외부 가속도를 조절하기 위해 PI 제어기를 사용하여 z축 방향의 목표 유속을 유지하였으며, 포아줄 및 코우에트 유동을 시뮬레이션하였다.
• 초고성능 컴퓨터(cacau 및 HP XC6000)를 사용하여 최대 480만 개의 상호작용 중심을 포함한 시스템에서 등엔트로피 상태에서 시뮬레이션을 수행하였다.

실험 결과

연구 질문

• RQ1 액체 메탄이 고체 탄소 나노채널 내에서 슬립 지배 나노스케일 흐름에서 연속체 유사 행동으로의 흐름 거동 전이 과정은 어떻게 이루어지는가?
• RQ2 다양한 채널 폭에서 봉인된 메탄 흐름에서 압력 강하-유속 관계를 Darcy의 법칙이 어느 정도 정확하게 기술하는가?
• RQ3 탄소 나노소재를 포함한 대규모 유체-고체 시스템에 대해 LS1/Mardyn MD 시뮬레이터의 확장성은 어떠한가?
• RQ4 다양한 공간 분할 전략(등방성 대비 채널 기반)이 비균일한 MD 시뮬레이션에서 성능에 어떤 영향을 미치는가?
• RQ5 메탄이 고체 탄소 채널 내에서 나노스케일과 연속체 영역 사이의 전이가 일어나는 데 필요한 임계 채널 폭 범위는 무엇인가?

주요 결과

• 액체 메탄이 고체 탄소 채널 내에서 슬립 지배 유동에서 연속체 유사 유동으로의 전이가 20에서 50 nm 사이의 좁은 폭 범위에서 발생한다.
• Darcy의 법칙은 고체 탄소 채널 내 메탄 흐름에서 나노스케일 및 연속체 영역 모두에서 압력 강하-유속 관계를 정확하게 기술한다.
• 2 nm 이하의 채널에서는 경계 슬립 속도가 평균 유속의 99% 이상을 초과하여, 가장 작은 척도에서 강한 비연속체 효과가 확인되었다.
• LS1/Mardyn MD 시뮬레이터는 초고성능 컴퓨터에서 뛰어난 확장성을 보였으며, 최대 480만 개의 상호작용 중심을 포함한 시스템을 처리하였다.
• 채널 기하학 기반의 도메인 분할 전략이 비균일한 유체-고체 시스템에서 등방성 분할 전략보다 우수했으며, 부하 균형 및 통신 효율성을 향상시켰다.
• 재조정된 Tersoff 잠재 에너지를 사용한 탄소 벽면은 결합 길이를 1.461 Å에서 1.421 Å로 보정하여 구조적 정확도를 향상시켰다.