테일러-후드 요소

1. 개요[편집]

테일러-후드 요소 Taylor–Hood element
정의	속도를 압력보다 한 차수 높게 잡은 연속 압력 쌍
대표형	삼각형 P2–P1 · 사각형 Q2–Q1
제안	테일러 & 후드 (1973)
안정성	이산 inf-sup(LBB) 만족 — 증명은 1980년대
수렴차수	속도 H¹ O(h²) · 속도 L² O(h³) · 압력 L² O(h²)
2D 비용	절점당 미지수 9개 (P1–P1 은 3개)
약점	압력이 연속 → 요소별 질량보존 안 됨, 압력 견고성 없음

속도를 한 차수 올려라. 끝. 이 한 줄이 50년을 버텼다.

테일러-후드 요소는 속도와 압력을 모두 연속 다항식으로 보간하되 속도의 차수를 압력보다 정확히 하나 높게 잡은 유한요소 쌍($\mathrm{P}_k$–$\mathrm{P}_{k-1}$, $k\ge2$)이다. 가장 널리 쓰이는 것은 삼각형 위의 P2–P1(6절점 이차 속도 + 3절점 선형 압력)과 사각형 위의 Q2–Q1(9절점 이중이차 속도 + 4절점 이중선형 압력)이며, 스토크스 문제와 비압축성 나비에-스토크스 방정식의 유한요소 해석에서 사실상 기준선 노릇을 한다.

왜 속도-압력 쌍을 골라야 하는지, 잘못 고르면 무엇이 부서지는지, 안장점 계를 어떻게 푸는지는 혼합유한요소법 문서가 다룬다. 이 문서는 「그 쌍 하나」에 집중한다 — 정확히 왜 이 조합이 통과하고 등차수 조합이 실패하는지, 어떤 차수로 수렴하며 얼마를 지불하는지, 3차원에서 무엇이 달라지는지.

2. 요소의 생김새[편집]

쌍	속도	압력	비고
P2–P1 (2D)	삼각형 6절점, 성분당 6	꼭짓점 3절점	표준 중의 표준
Q2–Q1 (2D)	사각형 9절점	꼭짓점 4절점	이중이차. 8절점 세렌디피티가 아님
P2–P1 (3D)	사면체 10절점	꼭짓점 4절점	3차원 표준
Q2–Q1 (3D)	육면체 27절점	꼭짓점 8절점	비싸다
P3–P2 등	고차	한 차수 낮게	매끄러운 해에서 유리

구현상 요긴한 성질이 하나 있다. 압력 절점이 속도 절점의 부분집합(꼭짓점)이므로 격자를 두 벌 만들 필요가 없다. 같은 격자 위에서 속도는 꼭짓점 + 변 중점, 압력은 꼭짓점만 쓰면 끝이고, 이 편의성이 채택률의 절반쯤을 설명한다. 반면 압력이 연속이라는 사실은 뒤에서 볼 대가를 부른다.

이름은 테일러(C. Taylor)와 후드(P. Hood)가 1973년 나비에-스토크스 유한요소 논문에서 제시한 데서 왔다. 후드-테일러라고 쓴 문헌도 흔하다. 정작 제안 당시에는 안정성 증명이 없었다 — 되니까 쓴 것이고, 왜 되는지는 10년쯤 뒤에 정리되었다.

3. 왜 P1–P1 은 실패하는가[편집]

이 쌍을 이해하는 가장 빠른 길은 실패 사례를 정확히 해부하는 것이다. 속도-압력 결합항을 보자.

$$b(\mathbf{v},q)=-\int_\Omega q\,\nabla\!\cdot\!\mathbf{v}\;\mathrm{d}x$$

속도를 P1(조각별 선형)으로 잡으면 $\nabla\!\cdot\!\mathbf{v}_h$ 는 요소마다 상수다. 그러면 요소 $K$ 별로 묶어서

$$b(\mathbf{v}_h,q_h)=-\sum_K \left(\nabla\!\cdot\!\mathbf{v}_h\right)\big|_K \int_K q_h\,\mathrm{d}x =-\sum_K \left(\nabla\!\cdot\!\mathbf{v}_h\right)\big|_K\,|K|\,\bar q_K$$

결합항이 압력의 요소별 평균 $\bar q_K$ 밖에 보지 못한다. 압력을 P1 으로 잡아 요소 안에서 기울어질 수 있게 해 줬는데, 그 기울기 정보는 방정식에 아예 안 들어간다는 뜻이다. 요소 평균이 모두 0이면서 0이 아닌 압력장이 존재하면 그것은 완전한 가짜 모드이고, 정확히 0이 아니더라도 그런 장에 가까운 방향이 있으면 이산 inf-sup 상수 $\beta_h$ 가 격자를 조일수록 0으로 떨어진다. 결과가 인접 절점마다 부호가 뒤집히는 압력 체커보드다.¹

자유도를 세어 보면 사정이 더 분명하다. 2차원 삼각형 격자에서 절점 수 $V$ 에 대해 요소 수는 대략 $2V$ 다. P1 속도의 벡터 자유도는 $V$ 개인데, 위 식이 요구하는 “요소별 평균을 제어하라”는 조건은 $2V$ 개다. 제어해야 할 것이 제어 수단보다 두 배 많다.

같은 렌즈로 P2–P1 을 보면 처방이 자명해진다. 속도를 P2 로 올리면 $\nabla\!\cdot\!\mathbf{v}_h$ 가 요소마다 선형이 되어, 압력의 요소 내 기울기까지 결합항이 볼 수 있게 된다. 그리고 P2 의 변 중점 자유도가 요소당 평균 1.5개(변 수 $\approx3V$, 요소 수 $\approx2V$)이므로 요소당 하나씩 필요한 제약을 감당할 여유가 생긴다.

4. 안정성의 정확한 근거 — 포르탱 연산자[편집]

위의 논증은 감각이고, 증명은 포르탱 기준(Fortin’s criterion)으로 한다. 연속 문제의 inf-sup 조건이 성립할 때, 다음을 만족하는 사상 $\Pi_h:\mathbf{V}\to\mathbf{V}_h$ 가 있으면 이산 inf-sup 조건이 따라온다.

$$b(\mathbf{v}-\Pi_h\mathbf{v},\,q_h)=0\quad\forall q_h\in Q_h, \qquad \lVert\Pi_h\mathbf{v}\rVert_1\le C\lVert\mathbf{v}\rVert_1$$

핵심은 첫 조건을 테일러-후드에서 아주 구체적으로 만들 수 있다는 것이다. $q_h$ 가 연속이고 $\mathbf{v}\in\mathbf{H}^1_0$ 이므로 부분적분하면 경계항이 사라져

$$b(\mathbf{v},q_h)=\int_\Omega \nabla q_h\cdot\mathbf{v}\;\mathrm{d}x =\sum_K \nabla q_h\big|_K\cdot\int_K \mathbf{v}\,\mathrm{d}x$$

$q_h\in\mathrm{P}_1$ 이면 $\nabla q_h$ 는 요소별 상수다. 따라서 모든 요소에서 $\int_K(\mathbf{v}-\Pi_h\mathbf{v})\,\mathrm{d}x=\mathbf{0}$ 이기만 하면 첫 조건이 성립한다. 그리고 P2 의 변 중점 자유도는 요소 내부 적분값을 조절하기에 충분하다 — 클레망/스콧-쟝 보간으로 $\mathbf{H}^1$ 유계성을 확보한 뒤, 변 중점 값만 살짝 고쳐 요소별 평균을 맞추면 된다. 요소당 벡터 제약 1개, 쓸 수 있는 변 자유도 1.5개. 앞 절의 자유도 세기가 여기서 증명의 뼈대로 그대로 재등장한다.

역사적으로는 볼런드-니콜라이데스(1983)의 매크로요소 논법이 표준 도구가 되었다. 격자를 몇 개 요소씩 묶은 패치 위에서 국소 조건을 확인하면 전역 안정성이 따라온다는 것으로, 이 논법 덕에 “어떤 격자에서 성립하는가”가 국소적으로 판정 가능해졌다. 브레치-포크(1991)가 고차 $\mathrm{P}_k$–$\mathrm{P}_{k-1}$ 와 3차원까지 정리했다.

격자 조건이 붙는다는 사실은 알아 둘 만하다. 병리적으로 성긴 격자(요소가 몇 개뿐이거나 특정 퇴화 패치가 있는 경우)에서는 조건이 깨질 수 있다. 실무 격자에서는 사실상 언제나 만족하므로 걱정할 일이 아니지만, “테일러-후드는 무조건 안정”이라는 문장은 엄밀히는 틀렸다.

5. 수렴차수와 비용[편집]

매끄러운 해에 대해 $\mathrm{P}_k$–$\mathrm{P}_{k-1}$ 은 다음을 준다.

$$\lVert\nabla(\mathbf{u}-\mathbf{u}_h)\rVert_0=O(h^{k}),\qquad \lVert p-p_h\rVert_0=O(h^{k}),\qquad \lVert\mathbf{u}-\mathbf{u}_h\rVert_0=O(h^{k+1})$$

$k=2$ 면 속도 에너지 노름과 압력이 모두 2차, 속도 $L^2$ 가 3차다. 속도와 압력의 정확도가 같은 차수로 균형 잡힌다는 것이 이 쌍의 설계 의도이고, 그래서 “속도만 좋고 압력은 못 쓰는” 불안정 쌍과 대비된다.

주의할 조건이 둘 있다. 첫째, 속도 $L^2$ 의 $O(h^3)$ 는 곡선 경계를 직선 변으로 근사하면 그 자리에서 죽는다. 경계가 곡면이면 등매개변수 요소로 이차 사상을 써야 하며, 이건 선택이 아니라 차수를 지키기 위한 필수 조건이다. 둘째, 각힌 구석이나 재진입 모서리가 있으면 해의 정칙성이 부족해 위 차수가 안 나온다 — 요소 탓이 아니라 문제 탓이다.

비용은 정직하게 세어 보는 것이 좋다. 2차원 삼각형 격자($E\approx3V$, $T\approx2V$)와 3차원 사면체 격자($E\approx7V$, $T\approx6V$)에서 절점 수 $V$ 당 총 미지수는

차원	P2–P1	P1–P1(안정화)	비율
2D	$2(V+3V)+V=9V$	$2V+V=3V$	3배
3D	$3(V+7V)+V=25V$	$3V+V=4V$	6배

3차원에서 여섯 배라는 숫자가 이 문서에서 가장 실무적인 한 줄이다. 게다가 P2 는 절점 간 연결이 조밀해 행렬의 0 아닌 원소 수가 자유도 수보다 더 빠르게 늘고, 조립·구적 비용도 함께 오른다. 대규모 3차원 산업 해석이 테일러-후드 대신 등차수 + SUPG/PSPG 안정화로 기우는 이유는 이론적 우월성이 아니라 이 표다.

6. 대안들과의 자리 배치[편집]

• MINI 요소(P1+버블–P1). P1–P1 의 부족분을 요소 내부 버블 하나로 정확히 메운다. 버블은 정적 축약으로 요소 단계에서 지워지므로 전역 자유도가 거의 안 늘어, 2차원에서 절점당 3~4개 수준으로 끝난다. 대신 수렴차수는 1차로 내려간다 — 싸지만 정확하지 않다.
• P2–P0, Q2–P1 불연속. 압력을 불연속으로 두면 요소별 질량보존이 정확히 성립한다. Q2–P1 불연속은 상용 CFD 코드가 즐겨 쓰는 조합이고, 대신 압력 자유도가 요소당 3개라 테일러-후드보다 비싸다.
• 안정화 등차수 쌍(PSPG·GLS·국소 사영). 자료구조가 단순하고 3차원 비용이 압도적으로 싸지만 매개변수 $\tau$ 라는 노브가 생긴다. 테일러-후드의 진짜 장점은 정확도보다도 “튜닝할 것이 없다”는 데 있다.
• 스콧-보글리우스 $\mathrm{P}_k$–$\mathrm{P}_{k-1}$ 불연속. 속도가 정확히 발산 0이 되는 쌍이다. 일반 2차원 격자에서는 $k\ge4$ 가 필요하고, 무게중심 세분(Alfeld) 격자에서는 더 낮은 차수로도 성립한다. 아래에서 볼 압력 견고성 문제의 정공법이지만 격자 생성에 제약이 붙는다.

7. 압력 견고성 — 테일러-후드의 약한 고리[편집]

압력이 연속이라는 대가가 여기서 청구된다. 테일러-후드의 이산 속도는 점별로도 요소별로도 발산이 0이 아니다. 약형식에서만 질량보존이 성립하고, 그 결과 스토크스 문제의 오차 상한이 다음 모양이 된다.

$$\lVert\nabla(\mathbf{u}-\mathbf{u}_h)\rVert_0\ \le\ C\left(\inf_{\mathbf{v}_h}\lVert\nabla(\mathbf{u}-\mathbf{v}_h)\rVert_0 +\frac{1}{\nu}\inf_{q_h}\lVert p-q_h\rVert_0\right)$$

속도 오차 안에 압력 근사 오차가 $1/\nu$ 배로 들어온다. 압력 기울기가 크고 점성이 작은 문제 — 부력이 지배하는 자연대류, 큰 원천항이 있는 흐름 — 에서는 이 항이 속도해를 통째로 오염시킨다. 물리적으로 회전이 없는 힘은 압력에만 흡수되어야 하는데, 이산 공간이 그 분리를 정확히 못 해내는 것이다. 이 성질이 없는 이산화를 압력 견고(pressure-robust)하다고 부르며, 스콧-보글리우스 같은 정확히 발산 없는 쌍이나 발산 없는 재구성 연산자가 그 처방이다.

값싼 완화책이 grad-div 안정화다. 운동량 식에 $\gamma\,(\nabla\!\cdot\!\mathbf{u}_h,\nabla\!\cdot\!\mathbf{v}_h)$ 항을 더해 발산 오차에 벌점을 매기는 것으로, 위 상한의 $1/\nu$ 계수를 실질적으로 눌러 준다. 대가는 행렬이 조밀해지고 $\gamma$ 라는 노브가 하나 생긴다는 것 — 안정화를 피하려고 테일러-후드를 골랐는데 결국 안정화를 하나 붙이게 되는 흔한 아이러니다.

8. 구현할 때 실제로 물리는 것들[편집]

• 압력의 상수 모드. 속도 경계조건이 전부 디리클레면 압력은 상수만큼 부정이라 계가 특이하다. 절점 하나를 고정하는 방식이 가장 쉽지만 그 절점 주변에 국소 오차를 남기고 조건수를 나쁘게 만든다. 평균 0 제약(라그랑주 승수 한 개를 추가하거나 반복 솔버에서 매 스텝 사영)이 정석이다.
• 구적 차수. 스토크스 항만 보면 P2 기울기가 P1 이라 낮은 차수 구적으로도 정확하지만, 나비에-스토크스의 대류항과 곡선 등매개 사상이 들어오면 피적분함수 차수가 훌쩍 오른다. 여기서 수치적분 차수를 아끼면 안정성이 조용히 갉아먹힌다.
• inf-sup 조건은 대류를 다뤄 주지 않는다. 테일러-후드는 압력에 대해 안정할 뿐, 높은 페클레 수·레이놀즈수에서는 여전히 대류 안정화(SUPG/GLS)나 충분히 조밀한 격자가 필요하다. “안정한 요소를 썼는데 진동한다”의 대부분이 이 혼동이다.
• 솔버 궁합. 안정한 쌍을 쓰면 슈어 보수가 압력 질량행렬과 스펙트럼 동등해져, 압력 질량행렬을 전처리기로 쓴 우자와 알고리즘·MINRES의 반복 수가 격자 크기에 무관해진다. 이산화 안정성이 곧 솔버 성능으로 번역되는 대목이며, $\beta_h\to0$ 인 쌍에서는 이 좋은 성질이 함께 사라진다.
• 자기 검사. 새 코드를 짰으면 수치 inf-sup 검사(일반화 고유값 문제의 최소 0 아닌 고윳값이 격자를 조여도 상수로 남는지)를 한 번 돌려 보는 것이 좋다. 절차는 혼합유한요소법 문서에 있다. 그리고 제조해(MMS)로 위 표의 차수 $2,2,3$ 이 실제로 나오는지 확인하면 구적·경계 사상 실수가 대부분 잡힌다.²

대부분의 오픈소스 유한요소 프레임워크(FEniCS, Firedrake, deal.II, MFEM 등)가 테일러-후드를 스토크스 예제의 기본값으로 싣고 있는 것은 우연이 아니다. 한 줄로 설명되고, 튜닝 노브가 없고, 차수가 예측대로 나오는 쌍이라 교재와 회귀시험에 이보다 나은 후보가 없다. 산업 현장의 3차원 대규모 해석에서는 안 쓰이더라도, 어떤 코드가 옳게 짜였는지 판정하는 기준선 자리는 여전히 이쪽이다.³

9. 관련 문서[편집]

• 혼합유한요소법 · 유한요소법 · 갤러킨 방법 · 형상함수
• 등매개변수 요소 · 요소 잠김 · 불연속 갤러킨법 · 수치적분
• 스토크스 유동 · 나비에-스토크스 방정식 · SIMPLE 알고리즘 · 페클레 수
• 슈어 보수 · 우자와 알고리즘 · 안장점 · GMRES · 전처리기
• SUPG · 정적 축약 · FEniCS · 다중격자법

10. Footnotes[편집]

1. 여기서 흔한 오해 하나. 체커보드는 “격자가 성겨서” 생기는 것이 아니라 이산화가 압력을 유일하게 정하지 못해서 생긴다. 그래서 격자를 조여도 안 없어지고, 오히려 진동의 파장만 짧아진다. 격자를 두 배로 조였는데 그림이 똑같이 지저분하면 수렴 문제가 아니라 요소 선택 문제라는 신호다. ↩

2. P2 삼각형을 처음 구현하는 사람이 반드시 한 번은 겪는 버그가 변 자유도의 방향 문제다. 인접한 두 요소가 공유하는 변에 대해 국소 번호가 서로 반대로 매겨지면 그 변의 속도가 두 요소에서 다른 값으로 조립되고, 결과는 “수렴은 하는데 차수가 1차로 나오는” 미묘한 실패다. 발산은 안 하고 그림도 그럴듯해서 제조해 수렴차수 시험을 돌리기 전까지는 아무도 모른다. 그래서 차수 시험은 선택이 아니다. ↩

3. 재미있는 역사적 사실 하나. 테일러와 후드는 1973년 논문에서 이 조합이 “잘 작동한다”고 보고했을 뿐 왜 그런지는 몰랐다. LBB 조건 자체가 정리된 것이 1971~1974년이고, 이 쌍이 그 조건을 만족한다는 증명은 1980년대에야 나왔다. 즉 실무가 이론보다 10년 앞서 정답을 골랐다. 수치해석에서 이런 순서는 생각보다 흔하고, 문제는 같은 방식으로 고른 다른 조합들 중 상당수가 틀렸다는 것이다. 되는 걸 고르는 것과 왜 되는지 아는 것의 차이가 여기 있다. ↩