직접 시뮬레이션 몬테카를로

전산유체역학 시뮬레이션 유체역학 마지막 수정: 2026-07-28 04:38:52

1. 개요[편집]

직접 시뮬레이션 몬테카를로(Direct Simulation Monte Carlo, DSMC)는 기체를 분자 집단으로 직접 표현하되 분자 간 충돌을 확률적으로 표본추출해 볼츠만 방정식의 해를 얻는 입자 기반 수치기법이다. 1963년 그레이엄 버드(G. A. Bird)가 제안했고, 오늘날 크누센수 전이 영역( $0.1 \lesssim \mathrm{Kn} \lesssim 10$ )에서는 사실상 유일한 실용 해법으로 통한다.

핵심 아이디어는 두 가지다. 첫째, 시뮬레이션 입자 하나가 실제 분자 $F_N$ 개를 대표한다. $F_N$ 은 보통 $10^{12}\sim10^{18}$ 이라, 실제로는 아보가드로 수 규모인 계를 입자 수백만 개로 흉내 낼 수 있다.¹ 둘째, 이동과 충돌을 분리한다. 한 시간 스텝 동안 모든 입자를 충돌 없이 직선으로 날린 뒤, 셀 단위로 충돌 쌍을 확률적으로 뽑아 속도를 바꾼다. 분자 간 힘을 적분하지 않으니 분자동역학보다 몇 자릿수 싸고, 그 대가로 희박 기체(이항 충돌 지배) 가정에 묶인다.

이 분리가 정당한 이유는 희박 기체에서 평균 충돌 시간이 분자 충돌 자체의 지속 시간보다 압도적으로 길기 때문이다. 그래서 시간 스텝은 평균 충돌 시간보다 작아야 하고, 셀 크기는 평균자유행로보다 작아야 한다 — 이 두 줄이 DSMC 격자 설계의 전부라고 해도 과언이 아니다.²

2. 알고리즘 한 사이클[편집]

이동: 모든 입자를 $\mathbf{x} \leftarrow \mathbf{x} + \mathbf{v}\,\Delta t$ 로 옮긴다. 도중에 벽을 만나면 그 자리에서 반사시킨다.
재색인: 옮겨진 입자를 다시 셀에 배정한다(정렬·해시).
충돌: 각 셀 안에서 충돌 쌍을 표본추출하고, 뽑힌 쌍의 속도를 충돌 후 값으로 갱신한다. 공간 정보는 쓰지 않는다 — 같은 셀이면 어디 있든 충돌 후보다.
표본화: 셀별로 $\rho$ , $\mathbf{u}$ , $T$ , 응력, 열유속을 입자 모멘트로 누적한다.

정상 유동이면 과도기가 지난 뒤 수만~수십만 스텝을 계속 누적해 잡음을 줄이고, 비정상 유동이면 같은 조건을 여러 번 반복(앙상블)해 평균한다.

3. NTC 충돌 표본추출[편집]

충돌 단계의 표준은 버드의 NTC(no-time-counter) 방식이다. 셀 부피 $V_c$ , 현재 입자 수 $N$ , 시간 평균 입자 수 $\bar{N}$ 일 때 후보 쌍 개수를

N_{\text{cand}} = \frac{1}{2}\,\frac{N \bar{N} F_N (\sigma_T c_r)_{\max} \Delta t}{V_c}

로 계산하고, 각 후보 쌍을 무작위로 뽑아 확률

P_{\text{acc}} = \frac{\sigma_T c_r}{(\sigma_T c_r)_{\max}}

로 수락한다. 여기서 $c_r$ 은 상대속도, $\sigma_T$ 는 총 충돌 단면적. 이건 전형적인 기각 표본추출(rejection sampling)이며, 앞선 시대의 시간 카운터(TC) 방식보다 훨씬 균일하고 통계적으로 건전해서 완전히 대체했다. 수락된 쌍의 충돌 후 속도는 질량중심 속도를 보존한 채 상대속도 벡터를 산란 모형에 따라 회전시켜 얻는다.

4. 분자 모형[편집]

경질구(HS): $\sigma_T$ 가 상대속도와 무관. 구현은 제일 쉽지만 점성의 온도 의존성이 $\mu \propto T^{0.5}$ 로 고정돼 실제 기체( $\omega \approx 0.7\sim0.8$ )와 안 맞는다.
VHS(variable hard sphere, Bird 1981): 유효 지름을 $d \propto c_r^{-\nu}$ 로 두어 $\sigma_T \propto c_r^{-2\nu}$ . 산란은 여전히 등방이지만 $\omega = \nu + 1/2$ 로 점성의 온도 지수를 실측에 맞출 수 있다. DSMC의 사실상 기본값.
VSS(variable soft sphere): 산란각을 $\cos\chi = 2b_{\text{rel}}^{1/\alpha}-1$ 로 편향시켜, 점성뿐 아니라 확산계수(슈미트수)까지 동시에 맞춘다. 혼합기체에서 필수.

내부 자유도는 보르그나케-라슨(Borgnakke–Larsen) 현상론 모형으로 처리한다. 충돌의 일정 비율( $1/Z_{\text{rot}}$ , $1/Z_{\text{vib}}$ )만 비탄성으로 지정하고, 그 충돌에서는 병진+내부 에너지 총합을 보존한 채 평형 분포에서 재분배한다. 완화 수 $Z$ 가 곧 회전·진동 완화 시간을 결정하는 튜닝 손잡이고, 고온 재진입에서는 여기에 TCE(total collision energy) 화학반응 모형과 전리 모형까지 얹는다.

5. 경계조건과 통계 잡음[편집]

벽 조건은 경면 반사(법선 속도만 뒤집음, 마찰 없음)와 확산 반사(벽 온도의 맥스웰 분포에서 새 속도를 뽑음)를 수용계수 $\sigma$ 로 섞어 쓴다. 실제 공학 표면은 $\sigma \approx 0.8\sim1.0$ 로 거의 완전 확산에 가깝다. 더 정교하게는 CLL(Cercignani–Lampis–Lord) 모형으로 법선·접선 에너지 수용을 따로 준다. 유입 경계에서는 벽으로 향하는 반쪽 맥스웰 분포에서 유입 플럭스를 표본추출한다.

DSMC의 진짜 아킬레스건은 통계 잡음이다. 셀의 평균 속도는 열속도의 표본 평균이라, 표본 $M$ 개에서 상대오차는 대략

\frac{\epsilon_u}{u} \sim \frac{1}{\mathrm{Ma}\sqrt{M}}

로 스케일한다. 상온 공기의 열속도는 수백 m/s인데 MEMS 유동 속도는 mm/s 수준이라 $\mathrm{Ma} \sim 10^{-5}$ , 즉 필요한 표본 수가 $1/\mathrm{Ma}^2$ 로 폭증한다. 초음속 재진입에서는 아무 문제 없던 방법이 저속 미세유동에서는 못 쓸 물건이 되는 이유다.³ 이를 우회하려고 정보보존(IP) 방법, 그리고 평형 맥스웰 분포로부터의 편차만 입자로 표현하는 저분산 편차 DSMC(LVDSMC) 같은 변형이 나왔다.

계통 오차도 잊으면 안 된다. 셀 크기와 시간 스텝은 각각 수송계수에 2차 오차를 만들고, 셀당 입자 수가 너무 적으면 같은 입자 쌍이 반복 충돌해 상관이 생긴다. 관례는 셀 크기 $\lesssim \lambda/3$ , 시간 스텝 $\lesssim \tau_c/3$ , 셀당 입자 20~30개 이며, 이 세 방향 모두에 대해 수렴 확인을 해야 한다는 점에서 메시 독립성 연구의 정신이 그대로 적용된다.

6. 병렬화, 하이브리드, 그리고 응용[편집]

DSMC는 입자 단위 작업이라 공간 분할 병렬화가 자연스럽지만, 입자가 도메인을 넘나들며 부하 불균형을 만들기 때문에 동적 재분할이 필수다. 격자도 국소 $\lambda$ 에 맞춰 적응적으로 세분하는 게 표준(밀도가 100배 차이 나는 노즈 캡과 후류에 같은 셀을 쓸 수는 없다). 비용은 $\mathrm{Kn}$ 이 작아질수록 잔인하게 오른다 — 셀 수는 $\lambda^{-3}$ , 시간 스텝은 $\lambda$ 에 비례해 줄어서, 연속체에 가까워질수록 계산량이 발산한다. 그래서 연속체 붕괴 파라미터( $\mathrm{Kn}_{GL} > 0.05$ 같은 기준)로 영역을 나눠 CFD와 DSMC를 붙이는 하이브리드 해석이 정석이 됐다.

이론적 뒷받침도 있다. 입자 수를 무한대로 보내는 극한에서 DSMC가 볼츠만 방정식의 해로 수렴한다는 것이 1992년 바그너(Wagner)에 의해 증명됐다. “몬테카를로 땜질”이 아니라 일관된 이산화라는 뜻이다.⁴

대표 응용은 희박기체역학이 필요한 곳 전부다. 고고도 재진입 공력가열과 항력 예측, 우주선 자세제어 스러스터의 플룸 충돌과 태양전지판 오염, 진공 펌프·증착 장비의 컨덕턴스 설계, 마이크로 노즐, 혜성 코마·행성 외기권 같은 우주 환경. 코드로는 버드 본인의 DS2V/DS3V, 샌디아의 오픈소스 SPARTA, OpenFOAM 기반 dsmcFoam+ 등이 널리 쓰인다.

7. 비슷해 보이는 것들과의 구분[편집]

몬테카를로 방법: DSMC는 적분값을 추정하는 게 아니라 시간 전진하는 물리 시뮬레이션이고, 확률은 충돌 쌍 선택과 산란에만 개입한다. “MC로 적분한다”와는 결이 다르다.
이산요소법: DEM 입자 하나는 실제 알갱이 하나이고 접촉력은 결정론적으로 적분된다. DSMC 입자 하나는 분자 $10^{15}$ 개의 대표이고, 입자끼리 힘을 주고받지 않는다.
분자동역학: MD는 모든 분자의 궤적을 퍼텐셜로 적분한다. 밀집 기체·액체에서는 MD가 맞고 DSMC는 틀리며, 희박 기체에서는 그 반대로 MD가 감당 불가능하다.
플라즈마 시뮬레이션의 PIC: 장을 격자에서 풀고 입자를 밀어주는 구조는 비슷하지만, PIC의 상호작용은 자기장·전기장 경유이고 DSMC는 이항 충돌이다. 실제로 중성 입자 충돌이 중요한 전기추진 해석에서는 PIC-DSMC로 둘을 합쳐 쓴다.

8. 관련 문서[편집]

여기서 초심자가 자주 하는 오해가 ” $F_N$ 을 키우면 계산이 싸진다”는 것. 맞긴 한데, 셀당 입자 수가 20개 아래로 떨어지는 순간 충돌 통계가 무너지고 수송계수가 조용히 틀어진다. 싸진 건 계산 시간이지 물리가 아니다. ↩
이 두 줄을 어기면 어떻게 되냐면, 분자가 셀 여러 개를 건너뛰며 날아가 “만난 적도 없는 분자와 충돌”하게 된다. 결과는 물리적으로 그럴싸해 보이는데 점성만 슬쩍 틀린 해다. 세상에서 제일 위험한 종류의 오답이다. ↩
그래서 이 바닥에는 “DSMC는 빠를수록 쉽다”는 역설적인 격언이 있다. 마하 20 재진입은 며칠이면 끝나는데, 초당 1 mm로 기어가는 마이크로채널은 슈퍼컴퓨터를 몇 주 돌려도 잡음에 파묻힌다. ↩
버드 본인은 1963년 논문을 낼 때 이 수렴 증명이 없었다. 30년 가까이 “이게 왜 되는지”를 모른 채 NASA와 유럽우주국이 이 방법으로 재진입체를 설계했다는 뜻인데, 공학이 이론을 앞지르는 흔한 광경이다. ↩