다중입자 충돌 동역학을 이용한 지구 대기 가스 정수압 시뮬레이션

초록

본 연구는 다중입자 충돌 동역학(MPCD) 기법을 활용해 2차원 박스 안에서 중력만을 작용시키는 이상 기체의 정수압 평형을 모사한다. 주기적, 반사, 바운스백 경계조건을 적용하고, 다양한 중력 가속도(g=1~3.5) 하에서 5000 타임스텝까지 시뮬레이션하였다. 결과는 낮은 g에서는 밀도와 속도 분포에 큰 요동이 나타나지만, g≈1.5–2.5 구간에서 이론적 기체법칙과 스케일 높이(H) 식에 부합하는 정수압 프로파일을 재현한다.

상세 분석

이 논문은 MPCD를 대기 정수압 문제에 최초 적용한 시도라 할 수 있다. MPCD는 입자 기반 메소드로, 입자를 격자 셀에 모아 스토캐스틱 회전 연산을 수행함으로써 충돌을 모델링한다. 저자는 2차원 박스(10a×20a, a=1) 안에 평균 5개의 입자가 셀당 존재하도록 설정하고, 입자 질량을 m=1, 온도에 해당하는 평균 속도 v≈1.5958, 음속 c_s≈1.1832 로 정규화하였다. 초기 속도는 Box‑Müller 변환을 이용해 정규분포(μ=0, σ=thermal speed)에서 샘플링하였다.

경계조건은 물리적 의미를 고려해 설계되었다. 좌·우는 주기적(Periodic)으로 설정해 수평 방향 무한성을 구현하고, 상단은 스펙큘러 반사(Specular)로 입자와 수직 방향 반사, 하단은 바운스백(Bounce‑back)으로 완전 흡수/반사를 가정하였다. 이는 실제 대기 상부와 지표면의 반사 특성을 단순화한 모델이다.

중력 가속도는 입자 질량에 비례하도록 설정했으며, g 값을 1, 1.5, 2, 2.5, 3, 3.5 로 변화시켜 실험하였다. 시뮬레이션은 5000 타임스텝까지 진행했으며, 2500~5000 단계의 평균값을 밀도 프로파일에 사용했다. 에너지와 선형운동량 보존이 확인되었으며, 평균 절대 오차(MAE)를 통해 정수압 평형 도달 시점을 2000 타임스텝 전후로 판단하였다.

결과 분석에서 g≤1.5 구간은 밀도와 속도 분포에 큰 요동이 남아 정수압 이론과 차이가 크지만, g≈2~2.5에서는 이론적 지수 감소식 ρ(z)=ρ₀exp