격자 기체 오토마타

전산유체역학 계산물리 시뮬레이션 마지막 수정: 2026-07-25 04:13:02

1. 개요[편집]

격자 기체 오토마타(Lattice Gas Cellular Automata, LGCA 또는 lattice gas automaton, LGA)는 유체를 연속적인 장(field)이 아니라 격자 위를 정해진 속도로 뛰어다니는 불리언 입자들의 집합으로 보고, 극도로 단순한 이동(streaming)·충돌(collision) 규칙만으로 거시적인 나비에-스토크스 방정식 거동을 창발시키는 유체 해석 기법이다.

핵심 주장이 상당히 도발적이다. “미시 규칙의 세부는 중요하지 않다. 보존 법칙(질량·운동량)과 격자의 대칭성만 맞으면, 거시 극한에서는 어차피 나비에-스토크스가 나온다.” 실제로 그렇게 된다는 것이 1986년에 증명되면서 한동안 학계가 뒤집혔다. 오늘날 LGCA를 실무에서 직접 쓰는 사람은 드물지만, 이 방법의 직계 후손인 격자 볼츠만 방법이 산업 CFD의 한 축이 되었기 때문에 LGCA는 역사적 선조로서 반드시 짚고 넘어가야 하는 자리에 있다. 이 문서는 그 셀룰러 오토마타 관점과 계보에 집중하고, 분포함수 기반의 현대적 정식화는 격자 볼츠만 방법에 위임한다.

128×48 육각 격자에서 FHP-I 격자 기체를 비트 연산으로 돌린다. 노드마다 6방향 점유를 6비트로 들고, 정면 2체 충돌(±60도 무작위 회전)과 120도 간격 3체 충돌을 64칸 룩업 테이블로 적용한 뒤 각 비트를 자기 방향 이웃으로 한 칸 옮긴다. 위 패널은 노드당 입자 수를 그대로 찍은 불리언 세계이고, 아래 패널은 같은 순간을 B×B 블록으로 공간 평균한 뒤 시간 평균까지 건 속도장이다 — 왼쪽 흰 곡선이 벽 사이 ux 프로파일이다. 위아래 벽은 bounce-back, 좌우는 주기경계이며, 일정 유량을 유지하도록 매 스텝 -x 입자 몇 개를 +x 로 뒤집어 구동한다. 상단의 총 입자 수는 매 스텝 다시 센 값으로, 충돌도 스트리밍도 개수를 정확히 보존하는 한 초기값에서 변하지 않는다.

2. HPP 모형 — 최초의 시도와 그 한계[편집]

시작은 1973년 아디 하디, 이브 포모, 올리비에 드 파치스의 HPP 모형이다. 정사각 격자 위에서 각 노드는 상하좌우 네 방향 링크를 갖고, 링크당 입자는 0개 또는 1개만 존재한다(배타 원리). 규칙은 두 줄로 끝난다.

이동: 모든 입자가 자기 방향의 이웃 노드로 한 칸 이동한다.
충돌: 한 노드에서 정면으로 마주 보는 두 입자만 있을 경우, 이 쌍을 90° 회전시킨다. 나머지는 그대로 통과.

이 충돌은 입자 수와 총 운동량을 정확히 보존한다. 상태가 전부 비트라서 부동소수점 반올림 오차가 원리적으로 존재하지 않고(부동소수점 연산 걱정이 없다), 계산이 비트 연산 몇 개로 끝나며, 완벽하게 국소적이라 병렬 컴퓨팅에 이상적이다. 1970년대 하드웨어 사정을 생각하면 매력적인 조합이었다.

문제는 물리였다. HPP는 두 가지 치명상을 안고 있었다.

비등방성: 정사각 격자의 회전 대칭군은 4회 대칭뿐이다. 나비에-스토크스의 이류항을 복원하려면 4계 텐서 $\sum_i c_{i\alpha}c_{i\beta}c_{i\gamma}c_{i\delta}$ 가 등방이어야 하는데, 4회 대칭 격자는 이 조건을 만족하지 못한다. 결과적으로 거시 유동이 격자 축 방향을 기억한다 — 원기둥 뒤 후류가 격자에 맞춰 네모나게 나오는 식.
가짜 보존량(spurious invariants): HPP에서는 각 격자 행·열을 따라가는 운동량 성분이 별도로 보존된다. 진짜 유체에는 없는 보존량이라, 이것이 비물리적인 장기 상관을 만든다.

3. FHP 모형 — 육각격자가 답이었다[편집]

1986년 위리엘 프리슈, 브로스 하슬라허, 이브 포모가 낸 짧은 논문 하나가 판을 바꿨다. 격자를 정삼각형 격자(각 노드가 6개 이웃을 갖는 육각 대칭)로 바꾸기만 하면, 필요한 4계 텐서가 등방이 되어 거시 극한에서 제대로 된 나비에-스토크스가 나온다는 것이었다.

FHP-I은 노드당 6비트(6개 속도)를 쓴다. 충돌 규칙은 여전히 소박하다.

2체 정면충돌: 마주 보는 한 쌍을 ±60° 회전시킨다. 두 결과가 대칭이므로 무작위로(또는 격자 위치·시각에 따라 결정론적으로) 하나를 고른다. 이 무작위화가 없으면 또 다른 가짜 보존량이 생긴다.
3체 대칭충돌: 120° 간격의 세 입자가 만나면 나머지 세 방향으로 반사한다.

여기에 정지 입자를 추가하면 FHP-II(7비트), 가능한 충돌을 전부 넣으면 FHP-III이 된다. 충돌 규칙이 많을수록 점성이 낮아져(= 유효 레이놀즈수가 올라가) 더 흥미로운 유동을 볼 수 있다. 구현은 노드 상태 7비트를 인덱스로 하는 룩업 테이블 한 방이면 끝나므로, 곱셈 없이 유체를 푸는 진기한 코드가 된다.¹

3차원은 훨씬 골치 아팠다. 3D 정규 격자 중에는 필요한 등방성을 주는 것이 없어서, 데뮈에르·랄망·프리슈는 4차원 면심입방(FCHC) 격자를 쓰고 이를 3D로 사영하는 우회로를 택했다. 노드당 속도 24개, 상태 $2^{24}$ 개 — 룩업 테이블이 곧바로 감당 불가 크기가 되는 지점이다.

4. 거시 극한 — 왜 나비에-스토크스가 나오는가[편집]

LGCA의 미시 동역학을 앙상블 평균하면 각 방향 점유 확률 $f_i = \langle n_i \rangle$ 에 대한 격자 볼츠만 방정식이 되고, 여기에 채프먼-엔스코그 전개를 적용하면 크누센수 1차에서 오일러 방정식, 2차에서 나비에-스토크스가 복원된다. 이때 점성은 미시 규칙(충돌 집합)이 결정하는 유도량이지 사용자가 넣는 물성치가 아니다.

다만 복원된 방정식은 진짜 나비에-스토크스와 미묘하게 다르다. 이류항 앞에 밀도 의존 계수 $g(d)$ 가 붙는다.

\partial_t (\rho \mathbf{u}) + g(d)\,\nabla\!\cdot(\rho \mathbf{u}\mathbf{u}) = -\nabla p + \nu\nabla^2(\rho\mathbf{u}), \qquad g(d) = \frac{D}{D+2}\cdot\frac{1-2d}{1-d}

여기서 $d = \rho/b$ 는 링크당 평균 점유율(페르미-디랙형 평형분포에서 온다), $b$ 는 노드당 속도 수. 이 $g \neq 1$ 때문에 갈릴레이 불변성이 깨진다. 밀도를 고정하고 저마하수 영역에 머물면 시간·점성 재규격화로 넘어갈 수 있지만, 밀도 변화가 큰 문제나 자유표면 문제에서는 구조적 약점이다. 게다가 $d=1/2$ 에서 $g=0$ , 그 위에서는 부호가 뒤집힌다(입자-구멍 대칭).

5. 잡음, 그리고 격자 볼츠만으로의 계승[편집]

LGCA의 진짜 실용적 약점은 통계적 잡음이었다. 각 링크의 상태가 0 아니면 1이므로, 한 노드의 순간 속도는 사실상 표본 크기 $b$ 인 표본평균이다. 상대 요동은 $1/\sqrt{N}$ 로 줄어들 뿐이라, 매끈한 속도장 하나를 얻으려면 수백~수천 개 노드와 시간 스텝에 걸쳐 평균해야 했다. 유효 해상도가 그만큼 깎이는 셈이다.²

해법은 결국 “평균을 미리 취해 버리자”였다. 1988년 맥나마라와 자네티는 불리언 점유수 $n_i \in \{0,1\}$ 를 실수 분포함수 $f_i \in [0,1]$ 로 바꿔 버렸다. 잡음이 통째로 사라졌고, 대신 부동소수점 연산과 반올림 오차가 돌아왔다. 이어 이게라와 히메네스가 충돌 연산자를 국소 평형 주위로 선형화해 룩업 테이블을 없앴고, 1992년 첸·첸·매튜스와 첸·다위, 그리고 콴·데뮈에르·랄망이 BGK 단일 완화시간 충돌을 도입하면서 오늘날의 격자 볼츠만 방법(LBGK)이 완성됐다. 평형분포를 마하수 2차까지 다항식으로 잡으면 $g \to 1$ 이 되어 갈릴레이 불변성 문제도 해소된다.

정리하면 계보는 이렇다.

단계	상태 변수	잡음	대표 문제
HPP (1973)	4비트	큼	비등방, 가짜 보존량
FHP (1986)	6~7비트	큼	갈릴레이 불변성 깨짐
LBE (1988)	실수 $f_i$	없음	반올림 오차, 안정성
LBGK (1992)	실수 $f_i$	없음	현대 격자 볼츠만 방법

그렇다고 LGCA가 완전히 박제된 것은 아니다. 불리언 특성상 열역학적 요동이 자연스럽게 들어 있다는 점 때문에 요동 유체역학이나 미소 스케일 문제에서는 오히려 장점이 되고, 복잡한 다공질 구조(다공성 매질 유동) 내부 유동이나 반응-확산계, 교통류 모형에서는 지금도 CA 접근이 살아 있다. 그리고 “국소 규칙에서 연속체 방정식이 창발한다”는 이 발상 자체가 계산물리의 지적 자산으로 남았다.³

6. 관련 문서[편집]

그래서 1980년대 후반에는 아예 LGCA 전용 하드웨어를 만들자는 시도까지 있었다. MIT의 CAM(Cellular Automaton Machine)이 대표적. 부동소수점 유닛 없이 비트 연산만으로 유체를 푸는 기계라니, 지금 봐도 낭만적이다. ↩
노드 하나의 속도를 그대로 찍어 보면 그냥 노이즈다. 그래서 초기 LGCA 논문의 그림들은 죄다 ” $16\times16$ 블록 평균, 100 스텝 평균” 같은 단서가 붙어 있다. 격자를 100만 개 깔고 실효 해상도는 만 개인 상황. ↩
프리슈·하슬라허·포모 논문(1986)이 나온 뒤 “이제 CFD는 끝났다”는 분위기가 잠깐 돌았다. 결과적으로 유한체적법은 멀쩡히 살아남았고 LGCA는 LBM으로 진화해 옆자리에 앉았다. 혁명이 기존 체제를 갈아엎는 대신 새 분야 하나를 추가로 만들어 낸 흔한 사례. ↩