갈릴레이 불변성

물리 수치해석 전산유체역학 마지막 수정: 2026-08-04 04:53:50

1. 개요[편집]

갈릴레이 불변성(Galilean invariance)은 등속으로 서로 스쳐 지나가는 모든 관성계에서 물리 법칙의 형태가 똑같다는 원리다. 좌표 변환

\mathbf{x}' = \mathbf{x} - \mathbf{V}t, \qquad t' = t, \qquad \mathbf{u}' = \mathbf{u} - \mathbf{V}

에 대해 방정식이 형태를 유지한다는 뜻이며, 뉴턴 역학과 나비에-스토크스 방정식은 이를 정확히 만족한다. 배 안에서 공을 위로 던지면 배가 멈춰 있든 등속으로 달리든 손으로 되돌아온다는, 갈릴레이가 1632년에 쓴 그 이야기다.

문제는 여기서부터다. 연속체 방정식이 갈릴레이 불변이라고 해서 그것을 이산화한 코드가 자동으로 불변인 것은 전혀 아니다. 격자를 까는 순간 우리는 “격자가 정지해 있는 계”라는 특권적 기준계를 하나 만들어 버린다. 그래서 똑같은 물리 문제를 배경 흐름 속도만 바꿔 두 번 풀면, 이론상 같아야 할 답이 다르게 나온다. 이 문서는 그 깨짐이 실제로 어디서 얼마나 일어나는지를 다룬다.

2. 연속체에서는 왜 성립하는가[편집]

핵심은 물질도함수가 변환에 대해 그대로라는 데 있다. $\mathbf{u}' = \mathbf{u}-\mathbf{V}$ 이고 $\nabla' = \nabla$ , $\partial_{t'} = \partial_t + \mathbf{V}\cdot\nabla$ 이므로

\frac{\partial}{\partial t'} + \mathbf{u}'\cdot\nabla' = \frac{\partial}{\partial t} + \mathbf{u}\cdot\nabla

가 된다. 시간 미분에서 튀어나온 여분의 $\mathbf{V}\cdot\nabla$ 를 이류항의 감소분이 정확히 상쇄한다. 압력 기울기와 점성 항은 속도의 차이와 미분에만 의존하므로 애초에 $\mathbf{V}$ 를 보지 못한다. 결국 나비에-스토크스는 형태를 그대로 유지한다.

거꾸로 말하면 갈릴레이 불변성이 깨지는 항은 언제나 속도의 절대값이 직접 등장하는 항이다. 이 관점을 쥐고 있으면 아래 사례들이 전부 같은 이야기의 변주라는 게 보인다.

주의할 점 하나. 이 불변성은 관성계끼리만 성립한다. 회전하는 계로 넘어가면 코리올리력·원심력이라는 겉보기 힘이 반드시 등장하며, 이는 근사의 문제가 아니라 정의상 그렇다(회전 좌표계 기법). 그리고 맥스웰 방정식은 갈릴레이 불변이 아니다 — 광속이 기준계에 무관하다는 사실과 충돌하기 때문이며, 이 모순을 해소한 것이 로렌츠 변환과 특수상대성이론이다.

3. 뇌터 정리가 붙여 주는 보존량[편집]

뇌터 정리에 따르면 연속 대칭 하나마다 보존량 하나가 대응한다. 시간 병진은 에너지, 공간 병진은 운동량, 회전은 각운동량이라는 표는 익숙한데, 갈릴레이 부스트에도 짝이 있다.

\mathbf{K} = \sum_i m_i\mathbf{x}_i - t\sum_i m_i\dot{\mathbf{x}}_i = M\mathbf{X}_{cm} - t\,\mathbf{P}

$d\mathbf{K}/dt = 0$ 은 곧 질량중심이 등속 직선운동한다는 진술이다. 이 보존량은 시간을 명시적으로 포함하는 다소 특이한 형태라 표에서 자주 빠지지만, 엄연히 뇌터 전하다.¹

실용적 함의는 기하 수치적분과 이어진다. 이산 라그랑지안이 부스트 대칭을 유지하도록 짜면 변분 적분기가 이 보존량도 이산 수준에서 지킨다. 분자동역학 코드가 총 운동량을 기계 정밀도로 붙잡고 있는지 확인하는 습관이 여기서 나온다 — 총 운동량이 새면 계 전체가 슬금슬금 드리프트하고, 그 순간부터 “정지한 상자”라는 전제가 무너진다.

4. 격자 볼츠만 — 마하수 제곱의 청구서[편집]

격자 볼츠만 방법의 표준 평형분포는 맥스웰-볼츠만 분포를 속도에 대해 2차까지만 자른 것이다.

f_i^{eq} = w_i\rho\left[1 + \frac{\mathbf{e}_i\cdot\mathbf{u}}{c_s^2} + \frac{(\mathbf{e}_i\cdot\mathbf{u})^2}{2c_s^4} - \frac{|\mathbf{u}|^2}{2c_s^2}\right]

여기 $\mathbf{u}$ 는 격자 기준계의 절대 속도다. 3차 이상의 항을 버렸으니 부스트에 대해 형태가 유지될 리 없고, 설령 3차 항을 되살려도 D2Q9·D3Q19 같은 표준 격자는 3차 모멘트 텐서를 연속체 값과 일치시킬 만큼의 속도 벡터를 갖고 있지 않다. 그 결과 운동량 플럭스에

\Delta\Pi_{\alpha\beta} \sim \partial_\gamma\!\left(\rho\, u_\alpha u_\beta u_\gamma\right)

형태의 입방 속도 결함(cubic velocity defect)이 남는다. 주도항 $\rho u_\alpha u_\beta$ 대비 상대 크기가 $O(u^2/c_s^2) = O(Ma^2)$ 이므로, LBM에 붙는 저마하수 제약 $Ma \lesssim 0.1$ 은 경험칙이 아니라 이 오차를 1% 수준으로 묶으려는 정량적 선택이다. 유도의 상세는 채프먼-엔스코그 전개 문서에 있고, 그 조상뻘인 격자 기체 오토마타에서는 훨씬 더 심각한 형태( $O(1)$ 계수 $g(\rho)\ne1$ )로 나타났다.

체감되는 증상은 이렇다. 정지한 유체 속을 속도 $U$ 로 지나가는 와류를 풀 때와, 같은 와류를 정지시켜 놓고 균일류 $U$ 를 불어 줄 때 — 물리적으로 완전히 같은 문제인데 — LBM은 서로 다른 답을 낸다. 처방으로는 3차 에르미트 항을 지원하는 고차 격자(D3Q27·D3Q39), 결함을 상쇄하는 보정항 추가, 그리고 중심 모멘트(cascaded/cumulant) 충돌 연산자가 있다. 마지막 것이 개념적으로 가장 깔끔한데, 모멘트를 격자 기준계가 아니라 국소 유체와 함께 움직이는 계에서 이완시켜 충돌 단계 자체를 갈릴레이 불변으로 만들기 때문이다.²

5. 이류 도식 — 절단오차는 계산틀을 안다[편집]

이건 LBM만의 병이 아니다. 가장 단순한 선형 이류 $u_t + a u_x = 0$ 에 1차 상류 차분을 쓰면 수정 방정식이 이렇게 나온다.

u_t + a u_x = \frac{a\,\Delta x}{2}\left(1 - \frac{a\,\Delta t}{\Delta x}\right)u_{xx} + O(\Delta x^2)

원래 PDE에는 없던 수치확산 계수가 $a\Delta x/2$ , 즉 속도에 비례한다. 원 방정식은 완벽히 갈릴레이 불변인데 그 이산화의 오차는 계산틀에 정면으로 의존하는 것이다. 고차 도식으로 올라가도 사정은 같아서, 2차 중심차분의 주도 분산 오차는 $\tfrac{a\Delta x^2}{6}u_{xxx}$ 로 역시 $a$ 에 비례한다. 애초에 상류 차분이라는 발상 자체가 “흐름의 방향”이라는 계산틀 의존적 정보로 스텐실을 고르는 것이라, 불변성이 남아 있을 여지가 없다.³

한 문장으로: 해는 불변이지만 오차는 불변이 아니다. 그리고 우리가 실제로 보는 것은 해가 아니라 해와 오차의 합이다.

천체물리 쪽의 유명한 시연이 있다. 켈빈-헬름홀츠 불안정 문제에 균일 속도를 얹어 전체를 부스트하면, 고정 격자 오일러 코드에서는 불안정 성장률이 눈에 띄게 달라지거나 아예 억눌린다. 반면 격자가 유체와 함께 움직이는 이동 격자 코드에서는 부스트 전후 결과가 거의 그대로다. 격자가 부여한 특권 기준계가 물리 결과를 바꾼 것이다.

6. 입자법은 어떤가 — SPH와 물질점법[편집]

SPH는 이 지점에서 태생적으로 유리하다. 커널이 상대 위치 $\mathbf{r}_i-\mathbf{r}_j$ 에만, 인공점성이 상대 속도 $\mathbf{v}_i-\mathbf{v}_j$ 에만 의존하고, 무엇보다 이류항을 이산화하지 않는다(입자가 유체와 함께 움직이므로 이류가 곧 입자 이동이다). 균일 부스트를 걸어도 모든 상호작용 항이 문자 그대로 동일하다. SPH 진영이 “우리는 갈릴레이 불변”이라고 자랑하는 근거가 이것이고, 실제로 정당한 자랑이다. 물론 SPH의 문제는 다른 데(인장 불안정, 접촉 불연속 처리)에 있다.

물질점법과 FLIP-PIC 유체 계열은 중간에 있다. 입자↔격자 전달에 쓰이는 형상함수는 단위분할( $\sum_p N_p = 1$ )을 만족하므로 균일 속도장은 전달 과정에서 정확히 재현된다. 즉 부스트 자체는 오차 없이 통과한다. 그러나 배경 격자는 여전히 정지해 있고, 입자가 셀 경계를 넘을 때 발생하는 셀 횡단 오차는 입자가 얼마나 빨리 셀을 가로지르느냐에 비례한다. 배경 속도를 키우면 단위 시간당 횡단 횟수가 늘어나므로 오차 축적률이 달라진다. 1차 형상함수의 기울기가 셀 경계에서 불연속인 것이 원인이라, GIMP·CPDI·B-스플라인 MPM 같은 매끄러운 형상함수 계열이 이 잡음을 줄이는 방향으로 발전했다.

7. 실무에서 이걸로 무엇을 하는가[편집]

첫째, 공짜 검증 시험이 된다. 문제 전체에 균일 속도를 더하고(초기조건과 경계조건 모두) 다시 푼 뒤, 함께 움직이는 계로 되돌려 원래 해와 비교한다. 차이가 곧 그 코드의 계산틀 의존성이며, 정답을 몰라도 되는 몇 안 되는 검증 및 확인 시험이다. LBM·저차 이류 도식·PIC 계열은 여기서 반드시 무언가를 뱉는다.

둘째, 기준계를 고를 수 있으면 골라야 한다. 오차가 $|u|$ 에 비례하니 평균 흐름을 0에 가깝게 두는 계에서 푸는 쪽이 유리하다. 회전하는 로터를 ALE 기법이나 이동 격자로 다루는 이유, 와류 이류 벤치마크를 함께 움직이는 계에서 돌리면 결과가 극적으로 좋아지는 이유가 모두 같다. 뒤집어 말하면, 논문의 눈부신 와류 보존 그림이 도식의 우수함이 아니라 기준계 선택의 산물일 수 있으니 배경 속도를 확인하고 볼 일이다.

셋째, 이동·곡선 격자에서는 자유류 보존이 최소 요건이 된다. 격자가 제멋대로 움직여도 균일한 흐름은 균일한 채로 남아 있어야 하고, 이를 이산 수준에서 보장하는 조건이 기하 보존 법칙(GCL)이다. 이것을 만족하지 못하는 이동 격자 코드는 아무것도 없는 균일류에서 가짜 파동을 만들어 낸다 — 갈릴레이 불변성 위반이 가장 노골적으로 드러나는 형태다.⁴

8. 관련 문서[편집]

시간이 명시적으로 들어간 보존량이라 처음 보면 반칙 같지만, 뇌터 정리의 요구 조건은 “라그랑지안이 전미분만큼만 변할 것”이지 “대칭 변환이 시간을 안 볼 것”이 아니다. 갈릴레이 부스트는 정확히 전미분만큼 라그랑지안을 바꾼다. ↩
“충돌은 국소 유체와 함께 움직이는 계에서 일어나야 한다”는 문장은 물리적으로 너무나 당연해서, 오히려 표준 LBM이 그걸 안 하고도 30년간 잘 굴러왔다는 사실이 더 신기하다. 답은 간단하다 — 다들 $Ma$ 를 0.1 아래로 눌러 놓고 썼다. ↩
그래서 “상류 차분은 물리적으로 자연스럽다”는 흔한 설명은 절반만 맞다. 정보 전파 방향을 존중한다는 점에서는 자연스럽고, 기준계를 고정한다는 점에서는 전혀 자연스럽지 않다. ↩
이 증상은 디버깅하기 좋은 편이다. 균일류를 넣고 돌려서 아무 일도 안 일어나면 통과, 뭐라도 움직이면 실패. 판정이 이렇게 명확한 시험은 CFD에 흔치 않다. ↩