ALE 기법

수치해석 시뮬레이션 전산유체역학 마지막 수정: 2026-07-17 04:23:15

1. 개요[편집]

ALE 기법(Arbitrary Lagrangian-Eulerian method, 임의 라그랑주-오일러 기법)은 계산 격자를 물질에 완전히 붙이지도(라그랑주), 공간에 완전히 고정하지도(오일러) 않고 제3의 임의 속도로 따로 움직이게 하는 이산화 기법이다. 격자에게 독립적인 인생을 허락한 것이 핵심이며, 그 대가로 방정식에는 대류항이 하나 더 붙는다.

이름의 “임의(Arbitrary)“가 전부를 말한다. 격자 속도 $\hat{\mathbf{u}}$ 를 물질 속도 $\mathbf{u}$ 와 같게 두면 라그랑주, 0으로 두면 오일러가 되므로 ALE는 두 고전 기술법을 특수 경우로 품는 상위 개념이다. 즉 새로운 물리가 아니라 **새로운 시점(kinematic description)**이다.¹

2. 왜 필요한가[편집]

연속체역학에는 원래 두 가지 시점이 있었고, 둘 다 결정적인 약점이 있었다.

라그랑주 기술 — 격자가 물질을 따라다닌다. 자유표면과 물질 경계가 자동으로 추적되고, 이력 의존적인 소성 구성방정식을 쓰기에 자연스럽다. 구조해석이 기본으로 채택하는 시점. 문제는 변형이 커지면 격자가 뒤집히거나 찌그러져 자코비안이 음수가 되고 해석이 그냥 죽는다는 것.
오일러 기술 — 격자가 공간에 못 박혀 있다. 아무리 격렬한 유동이 지나가도 격자는 멀쩡하다. 전산유체역학의 기본값. 대신 물질 경계가 격자 사이를 헤엄쳐 다니므로 계면이 수치 확산으로 뭉개진다.

ALE는 이 둘의 장점만 취하려는 시도다. 물질 경계는 라그랑주처럼 따라가되, 내부 격자는 스스로 재배치해서 품질을 지킨다. 유체-구조 연성에서 ALE가 사실상 표준인 이유가 이것이다. 구조물 표면에는 격자가 딱 붙어 있어야 벽면 전단응력을 제대로 뽑는데, 그 구조물이 펄럭이면 주변 유체 격자도 같이 움직여줘야 한다.

3. 지배 방정식[편집]

임의 속도로 움직이는 관찰자에 대한 물질 도함수를 다시 쓰는 것에서 출발한다. 격자 기준 시간미분과 물질 도함수의 관계는 다음과 같다.

\left. \frac{\partial f}{\partial t} \right|_{\mathbf{X}} = \left. \frac{\partial f}{\partial t} \right|_{\boldsymbol{\chi}} + \mathbf{c} \cdot \nabla f, \qquad \mathbf{c} \equiv \mathbf{u} - \hat{\mathbf{u}}

여기서 $\mathbf{c}$ 가 **대류 속도(convective velocity)**로, 물질이 격자를 상대로 얼마나 빠르게 스쳐 지나가는지를 뜻한다. 라그랑주에서는 $\mathbf{c} = 0$ , 오일러에서는 $\mathbf{c} = \mathbf{u}$ 다. 이걸 나비에-스토크스 방정식에 대입하면 ALE 형태의 운동량 방정식이 나온다.

\rho \left( \left. \frac{\partial \mathbf{u}}{\partial t} \right|_{\boldsymbol{\chi}} + (\mathbf{u} - \hat{\mathbf{u}}) \cdot \nabla \mathbf{u} \right) = -\nabla p + \nabla \cdot \boldsymbol{\tau} + \rho \mathbf{g}

유한체적법으로 움직이는 검사체적에 적분할 때는 반드시 기하 보존 법칙(Geometric Conservation Law, GCL)을 만족시켜야 한다.

\frac{d}{dt} \int_{V(t)} dV - \int_{S(t)} \hat{\mathbf{u}} \cdot \mathbf{n} \, dS = 0

말로 풀면 “균일한 유동장은 격자가 아무리 움직여도 균일하게 유지되어야 한다”는 지극히 당연한 요구다. 그런데 이걸 어기면 격자가 흔들리는 것만으로 없던 질량이 생기거나 사라진다. GCL 위반은 ALE 코드에서 가장 악명 높은 버그 원인이며, 증상이 “천천히 이상해지는” 형태라 잡기가 극도로 성가시다.²

4. 격자 이동 전략[편집]

격자 속도 $\hat{\mathbf{u}}$ 를 뭘로 줄 것인가 — 여기서 ALE의 실력이 갈린다. 경계에서의 이동량은 물리가 정해주지만, 내부 절점은 우리가 정해야 한다.

의사 탄성체(pseudo-solid) — 격자를 가상의 탄성체로 보고 라플라스 방정식이나 선형 탄성 방정식을 풀어 내부 절점을 배분한다. 작은 요소일수록 강성을 크게 주는 트릭(jacobian-based stiffening)을 얹으면 경계층 격자를 지켜낼 수 있다.
RBF 보간 — 경계 변위를 방사기저함수로 내부에 뿌린다. 연결 정보가 필요 없어 견고하지만 절점 수가 많아지면 밀집행렬 때문에 비싸진다.
라플라스 평활화 + 재격자 — 격자가 감당 못 할 만큼 찌그러지면 포기하고 메시 생성을 다시 한 뒤 해를 옮겨 심는다. 이 옮겨 심는 과정이 리맵(remap)이다.

연산자 분리형 ALE(operator-split ALE)는 이 발상을 아예 알고리즘으로 굳혔다. 한 스텝을 (1) 라그랑주 단계로 물리를 풀고, (2) 격자를 예쁘게 재배치하고, (3) 해를 새 격자로 보존적으로 리맵하는 셋으로 쪼갠다. 명시적 동해석 코드인 LS-DYNA와 폭발·관통 해석용 하이드로코드가 이 구조를 쓴다.

5. 적용 분야와 한계[편집]

ALE가 밥값을 하는 곳은 대체로 “물질 경계가 크게 움직이는데 그 경계가 중요한” 문제다.

분야	대표 사례	격자가 따라가는 대상
유체-구조 연성	심장 판막, 날개 플러터	고체 표면
자유표면 유동	슬로싱, 조파 저항	액면
금속 성형	압출, 절삭	공구-소재 계면
충격·폭발	관통, 조류 충돌	물질 계면
터보기계	회전-정지부 계면	블레이드

한계도 뚜렷하다. 위상(topology)이 바뀌는 순간 ALE는 무너진다. 물방울이 두 개로 갈라지거나 파도가 부서지며 공기를 가두는 순간, 격자를 계면에 붙여둔다는 전제 자체가 붕괴한다. 이런 문제는 VOF 방법이나 레벨셋 방법 같은 오일러 계면 포착법, 혹은 침지경계법이나 격자가 아예 없는 SPH로 가는 게 맞다. 그리고 리맵을 반복할수록 수치 확산이 쌓여 결국 오일러 해석과 비슷한 정확도로 수렴해간다는 점도 잊으면 안 된다.³

6. 현업에서의 현실[편집]

해석이 죽었을 때 물리를 의심하기 전에 격자 자코비안부터 찍어보게 된다. 열에 아홉은 요소 하나가 뒤집혔다.
격자 이동 방정식을 푸는 시간이 정작 유동을 푸는 시간보다 긴 경우가 흔하다. 배보다 배꼽.
“그냥 오일러로 하면 안 되나요?” — 벽면 마찰이나 응력과 변형률 이력이 필요 없다면 대체로 된다. 그러면 ALE 안 쓰는 게 이득이다.
상용 코드에서 ALE 옵션을 켜는 건 클릭 한 번이지만, 격자 이동 파라미터를 잡는 건 며칠이다. 수렴성은 여전히 신에게 맡긴다.

7. 관련 문서[편집]

그래서 “ALE 방정식”이라는 별도의 물리 법칙은 존재하지 않는다. 같은 보존 법칙을 움직이는 관찰자 좌표에서 다시 쓴 것뿐이다. 물리는 그대로인데 좌표계를 바꿨을 뿐이라는 점에서, 상대성이론을 배울 때의 그 킹받는 기분을 살짝 재현해준다. ↩
GCL을 만족시키는 가장 확실한 방법은 격자 속도를 해석적으로 주는 게 아니라, 실제로 쓸어낸 체적(swept volume)을 기하학적으로 계산해서 역산하는 것이다. 코드가 지저분해지지만 균일 유동 테스트를 통과시키는 유일한 길인 경우가 많다. ↩
리맵은 결국 보간이고, 보간은 결국 뭉개는 일이다. 차분 도식에서 상류차분이 수치점성을 만드는 것과 근본적으로 같은 이야기다. 공짜 점심은 여기에도 없다. ↩