격자 볼츠만 방법

전산유체역학 수치해석 마지막 수정: 2026-07-03 16:52:41

격자 볼츠만 방법 Lattice Boltzmann Method
약칭	LBM
이론적 기반	볼츠만 수송방정식
계산 스케일	중간자적 (mesoscopic)
대표 격자 모델	D2Q9, D3Q19, D3Q27
핵심 절차	충돌 (collision) · 스트리밍 (streaming)

1. 개요[편집]

유체를 통째로 풀지 말고, 아주 많은 가상 입자가 격자 위에서 춤추게 하라. 그러면 유체는 저절로 나타난다.

격자 볼츠만 방법(Lattice Boltzmann Method, LBM)은 유체를 연속체로 보고 나비에-스토크스 방정식을 직접 이산화하는 대신, 입자 분포함수의 볼츠만 수송방정식을 격자 위에서 푸는 중간자적(mesoscopic) 수치해석 기법이다. 기체를 이루는 개별 분자를 하나하나 추적하는 미시적(microscopic) 관점과, 밀도·속도장만 다루는 거시적(macroscopic) 관점의 딱 중간에 서서, 격자점마다 정의된 입자 분포함수 $f_i$ 의 시간 발전을 계산한다.

전통 CFD가 미분방정식을 이산화해 거대한 연립 대수방정식을 반복법으로 푸는 것과 달리, LBM은 매 시간 단계마다 “충돌”과 “스트리밍”이라는 국소적이고 단순한 두 연산만 반복한다. 압력-속도 연성을 위한 반복 계산(SIMPLE 알고리즘 같은)이 아예 필요 없고, 계산이 이웃 격자점끼리만 통신하는 국소 연산이라 병렬화가 극도로 쉽다. 이 단순함과 병렬성이 LBM을 GPU 시대의 총아로 만든 핵심 무기다.

2. 격자 속도 모델 (DnQm)[편집]

LBM의 무대는 균일한 정렬 격자다. 각 격자점에는 미리 정해진 유한 개의 이산 속도 방향 $\mathbf{e}_i$ 가 있고, 각 방향마다 그 방향으로 움직이는 입자의 양을 나타내는 분포함수 $f_i$ 가 얹힌다. 이 격자 구조를 DnQm 표기로 부르는데, $n$ 은 공간 차원, $m$ 은 이산 속도(방향) 개수다.

모델	차원	속도 방향 수	용도
D2Q9	2D	9	2차원 표준
D3Q15	3D	15	3차원 경량
D3Q19	3D	19	3차원 표준
D3Q27	3D	27	3차원 고정밀

가장 유명한 것이 2차원 표준인 D2Q9다. 정지 상태(1개), 상하좌우(4개), 대각선(4개)을 합쳐 9개 방향의 분포함수를 가진다. 각 방향에는 격자 대칭성과 등방성을 보장하도록 정해진 가중치 $w_i$ 가 부여되며, 이 가중치와 이산 속도 집합이 나중에 거시 방정식을 올바르게 복원하는 데 결정적 역할을 한다. 속도 방향이 많을수록(D3Q27 등) 정밀도와 안정성은 올라가지만 메모리와 연산이 함께 늘어나는, 익숙한 트레이드오프가 여기서도 성립한다.

3. 충돌-스트리밍[편집]

LBM의 심장은 놀랍도록 단순한 두 단계의 반복이다. 격자 볼츠만 방정식은 다음과 같이 쓴다.

f_i(\mathbf{x} + \mathbf{e}_i \Delta t,\; t + \Delta t) = f_i(\mathbf{x}, t) + \Omega_i

이 한 줄이 스트리밍(streaming)과 충돌(collision)로 분해된다.

충돌 단계: 각 격자점에서 분포함수 $f_i$ 가 국소적으로 재분배된다. 입자들이 부딪혀 평형을 향해 이완하는 과정으로, 위 식의 충돌항 $\Omega_i$ 가 담당한다. 순전히 국소 연산이라 이웃과 통신이 없다.
스트리밍 단계: 충돌 후의 분포함수가 각자의 이산 속도 방향 $\mathbf{e}_i$ 를 따라 이웃 격자점으로 그대로 이동한다. 그냥 값을 옆 칸으로 옮기는 것뿐이라 부동소수점 연산조차 거의 없다.

이 “부딪히고(collision) → 흘려보내고(streaming)“의 반복이 LBM 계산의 전부다. 무거운 행렬 역산도, 전역 압력 방정식도 없다. 이 간결함이 코드를 짜기 쉽게 만들고, 무엇보다 병렬화를 거의 공짜로 만든다.¹

4. BGK 근사[편집]

문제는 충돌항 $\Omega_i$ 다. 원래 볼츠만 방정식의 충돌항은 입자 간 충돌을 다루는 복잡한 적분항이라 그대로는 감당이 안 된다. 그래서 실무 LBM은 이를 극도로 단순화한 BGK 근사(Bhatnagar-Gross-Krook, 1954)를 쓴다. 핵심 아이디어는 “충돌이란 결국 분포함수를 국소 평형 분포 $f_i^{eq}$ 쪽으로 끌어당기는 이완 과정일 뿐”이라는 것이다.

\Omega_i = -\frac{1}{\tau}\left( f_i - f_i^{eq} \right)

여기서 $\tau$ 는 단일 이완시간(relaxation time)으로, 분포함수가 평형에 얼마나 빨리 수렴하는지를 조절한다. 놀랍게도 이 이완시간 $\tau$ 가 유체의 동점성계수 $\nu$ 와 직접 연결된다.

\nu = c_s^2 \left( \tau - \frac{1}{2} \right) \Delta t

즉 점성을 넣고 싶으면 $\tau$ 하나만 조절하면 된다. 단일 이완시간이라 SRT(Single Relaxation Time) 모델이라고도 부르며, 단순하지만 고레이놀즈수에서 안정성이 떨어지는 약점이 있어 이를 개선한 MRT(다중 이완시간), 정규화 LBM 등의 개량판이 다수 존재한다.

5. Chapman-Enskog와 나비에-스토크스 복원[편집]

여기서 자연스러운 의문이 든다. “격자 위에서 가상 입자 놀이를 하는 게 진짜 유체 방정식을 푸는 것과 무슨 상관인가?” 이 연결고리를 증명하는 것이 채프먼-엔스코그 전개(Chapman-Enskog expansion)다.

이 다중 스케일 점근 해석은 격자 볼츠만 방정식을 크누센수(Knudsen number)에 대해 전개하여, 거시 극한에서 분포함수의 모멘트가 정확히 질량 보존(연속 방정식)과 운동량 보존(나비에-스토크스 방정식)을 만족함을 보인다. 분포함수의 0차·1차 모멘트가 각각 밀도와 운동량에 대응한다.

\rho = \sum_i f_i, \qquad \rho \mathbf{u} = \sum_i \mathbf{e}_i f_i

다시 말해 LBM은 나비에-스토크스를 “직접” 풀지 않지만, 그 해를 “복원”한다. 미시적 규칙(충돌-스트리밍)에서 거시적 유체 거동이 창발하는 셈인데, 이 점근 등가성이 성립하려면 앞서 언급한 격자 등방성(D2Q9의 가중치 설계)과 저마하수(low Mach number) 조건이 필요하다. LBM이 본질적으로 약압축성·저속 유동에 맞춰진 방법이라는 사실이 바로 여기서 나온다.²

6. 강점과 약점[편집]

LBM은 만능이 아니다. 잘하는 것과 못하는 것이 뚜렷이 갈린다.

6.1. 강점[편집]

복잡 경계·다공성 매질: 계단식 경계로 복잡한 형상을 손쉽게 표현하고, 벽면 처리는 분포함수를 되튕기는 바운스백(bounce-back)이라는 단순 규칙으로 끝난다. 다공성 매질처럼 형상이 미쳐 날뛰는 문제에서 격자 생성 부담이 전통 CFD보다 훨씬 가볍다.
병렬성(GPU): 충돌은 완전 국소, 스트리밍은 최근접 이웃 통신뿐이라 GPU와 궁합이 환상적이다. 수백만 격자점을 그래픽 카드에 갈아 넣어 실시간에 가깝게 돌리는 사례가 많다.
다상·다성분: Shan-Chen 모델 등으로 상 분리를 자연스럽게 다룰 수 있어 다상유동 연구에서도 활발히 쓰인다.

6.2. 약점[편집]

압축성·고마하수: BGK-LBM은 태생이 약압축성 저속 유동용이다. 초음속이나 강한 압축성이 얽힌 문제는 표준 LBM의 사정권 밖이다.
고레이놀즈수 안정성: 레이놀즈수가 높아지면 이완시간이 $\tau \to 1/2$ 로 붙으면서 수치 불안정이 고개를 든다. MRT나 엔트로피 LBM 같은 처방으로 버티지만, 근본적으로 까다로운 영역이다.
균일 격자 제약: 전통적 LBM은 균일 정렬 격자를 전제해, 경계층처럼 국소적으로만 촘촘한 격자가 필요한 곳에서 국소 격자 세밀화(refinement)를 넣기가 유한체적법보다 번거롭다.

정리하면, LBM과 전통 CFD는 대체재라기보다 서로 다른 특기를 가진 도구다. 유한체적법 기반 상용 코드가 여전히 산업 표준의 왕좌에 앉아 있지만, 복잡 형상·다공성·GPU 대량 병렬이라는 특정 전장에서는 LBM이 확실히 한 방을 가지고 있다.³

7. 관련 문서[편집]

“행렬 역산이 없다”는 말은 CFD를 오래 한 사람에게는 거의 종교적 해방으로 들린다. 압력 포아송 방정식과 씨름해 본 적이 있다면 이 단순함의 가치를 안다. ↩
LBM이 저마하수 방법이라는 건 종종 오해를 부른다. “볼츠만 방정식 기반이니 희박 기체나 초음속에 강하겠지”라고 기대하고 왔다가, 정작 표준 LBM이 비압축성 저속 유동에 최적화되어 있다는 걸 알고 당황하는 것이 국룰. ↩
학회에서 LBM 진영과 전통 FVM 진영이 만나면 “그거 장난감 아니냐” vs “너희는 격자 짜다 인생 끝나잖아”의 신경전이 벌어지곤 한다. 물론 둘 다 각자의 문제에서 훌륭하게 작동한다. ↩