푸앵카레 부등식

1. 개요[편집]

푸앵카레 부등식 Poincaré Inequality
진술	‖u‖L² ≤ C ‖∇u‖L²
조건	경계에서 0 (프리드리히스) 또는 평균이 0 (비르팅거)
최적 상수	C = 1/√λ₁ (라플라시안 제1고유값)
차원	길이 — 영역 지름에 비례해 커진다
하는 일	강제성(coercivity) · 유일성 · 오차 추정
이산판	FV · DG의 안정성 증명 전제

푸앵카레 부등식은 함수의 크기를 그 기울기의 크기로 억제하는 부등식, 즉

$$\|u\|_{L^2(\Omega)} \;\le\; C_\Omega \,\|\nabla u\|_{L^2(\Omega)}$$

이 성립한다는 진술이다. 그냥 놔두면 이 부등식은 거짓이다 — $u \equiv 1$ 을 넣으면 좌변은 $|\Omega|^{1/2}$ 인데 우변은 0이다. 그래서 상수함수를 배제하는 조건이 반드시 붙는다. 조건을 어떻게 다느냐에 따라 두 판이 있고, 이 둘을 헷갈리는 것이 이 부등식에 얽힌 사고의 8할이다.

• 경계에서 0: $u \in H_0^1(\Omega)$. 흔히 프리드리히스 부등식(Friedrichs inequality)이라 부른다.
• 평균이 0: $u \in H^1(\Omega)$, $\bar u = \frac{1}{|\Omega|}\int_\Omega u = 0$. 푸앵카레-비르팅거 부등식(Poincaré–Wirtinger)이다.

겉보기엔 한 줄짜리 소박한 부등식인데, 이것 하나가 무너지면 유한요소법의 오차 이론도 타원형 문제의 유일성 증명도 통째로 무너진다. **“기울기 노름만으로 전체 $H^1$ 노름을 통제할 수 있다”**는 허가증이고, 그 허가증이 없으면 락스-밀그램 정리의 강제성 가정을 채울 방법이 없기 때문이다. 소볼레프 공간 문서가 이 부등식을 한 문단으로 스쳐 지나가므로, 여기서는 상수 $C_\Omega$ 의 정체와 이산 대응물까지 파고든다.

2. 두 부등식을 구분한다[편집]

용어가 문헌마다 섞여 있으니 정리해 두자. 좁은 의미의 “푸앵카레 부등식”은 평균 0 판을 가리키고, 경계 0 판은 프리드리히스라고 부르는 것이 정확하다. 다만 응용 문헌(특히 유한요소 쪽)에서는 둘 다 뭉뚱그려 “푸앵카레-프리드리히스 부등식”이라 쓰는 관행이 굳어져 있다. 어느 쪽이든 어떤 조건으로 상수함수를 죽였는가만 확인하면 된다.

경계 0 판은 증명이 세 줄이다. $u \in C_c^\infty(\Omega)$ 이고 $\Omega$ 가 $x_1$ 방향으로 폭 $d$ 인 띠 안에 들어 있다고 하자. 미적분의 기본정리와 코시-슈바르츠로

$$|u(x)|^2 = \Bigl|\int_{0}^{x_1} \partial_1 u \,dt\Bigr|^2 \le d \int_0^{d} |\partial_1 u|^2\,dt$$

을 얻고, 이것을 $\Omega$ 위에서 적분하면 $\|u\|_{L^2} \le d\,\|\nabla u\|_{L^2}$ 가 나온다. 밀도 논법으로 $H_0^1$ 전체로 확장하면 끝. 여기서 두 가지가 즉시 읽힌다 — **필요한 것은 유계성 전체가 아니라 “한 방향으로만 유계”**이고, 상수는 그 방향의 폭에만 의존한다. 무한 채널이나 무한 슬래브 위에서도 프리드리히스가 살아 있는 이유다.

평균 0 판은 사정이 다르다. 위와 같은 직접 증명이 안 되고, 보통 모순 논법 + 렐리히-콘드라쇼프 콤팩트성으로 증명한다. 부등식이 거짓이라 가정하면 $\|u_n\|_{L^2}=1$, $\|\nabla u_n\|_{L^2}\to 0$, $\bar u_n = 0$ 인 수열이 존재하고, $H^1$ 유계이므로 $L^2$ 에서 강수렴하는 부분수열을 뽑을 수 있으며, 그 극한은 기울기가 0이라 (연결 영역에서) 상수인데 평균이 0이므로 0이다. 그런데 $L^2$ 노름이 1이어야 하니 모순. 이 증명 구조가 두 가지를 폭로한다.

• 연결성이 필수다. 영역이 두 조각으로 떨어져 있으면 각 조각에서 다른 상수를 갖는 함수가 전체 평균 0이면서 기울기 0일 수 있다. 반례가 즉시 만들어진다.
• 상수의 크기에 대한 정보를 하나도 주지 않는다. 모순 논법이라 $C$ 가 존재한다는 것만 알려 준다. 상수를 알아야 하는 실무에서는 이게 치명적이라, 아래처럼 별도의 정리들이 필요해진다.

3. 상수 C는 어디서 오는가[편집]

먼저 차원 분석이 절반을 말해 준다. $u_\lambda(x) = u(x/\lambda)$ 로 영역을 $\lambda$ 배 키우면 $\|u_\lambda\|_{L^2}$ 은 $\lambda^{n/2}$ 배, $\|\nabla u_\lambda\|_{L^2}$ 은 $\lambda^{n/2-1}$ 배가 된다. 따라서

$$C_{\lambda\Omega} = \lambda\, C_\Omega$$

즉 $C$ 는 길이 차원이고 영역 지름에 비례한다. 차원이 몇이든, $p$ 가 몇이든 이 스케일링은 변하지 않는다. 실무에서 “이 부등식의 상수가 얼마쯤인가”를 물으면 첫 답은 언제나 “영역 크기 정도”다.

나머지 절반, 즉 모양 의존성은 고유값이 쥐고 있다. 최적 상수는 정의상

$$\frac{1}{C_\Omega^2} \;=\; \inf_{u \in H_0^1(\Omega)\setminus\{0\}} \frac{\|\nabla u\|_{L^2}^2}{\|u\|_{L^2}^2} \;=\; \lambda_1(\Omega)$$

이고, 우변의 레일리 몫 최소화는 정확히 디리클레 라플라시안의 제1고유값 문제 $-\Delta \phi = \lambda \phi$, $\phi|_{\partial\Omega}=0$ 이다. 즉 최적 푸앵카레 상수 = 1/√(제1고유값). 평균 0 판에서는 같은 논리로 $C = 1/\sqrt{\mu_2}$ 이고, $\mu_2$ 는 노이만 라플라시안의 첫 번째 0이 아닌 고유값이다(노이만의 $\mu_1 = 0$ 은 상수함수이고, 우리가 평균 0 조건으로 죽인 것이 바로 그 모드다). 부등식의 최적 상수와 고유값 문제가 같은 물건이라는 사실은 이 부등식을 다루는 방식을 통째로 바꾼다 — 상수를 추정하는 문제가 스펙트럼 하한을 추정하는 문제가 되기 때문이다.

구체적 값 몇 개는 외워 둘 만하다.

영역	상수
폭 $d$ 인 띠 안의 임의 영역, $H_0^1$	$C \le d/\pi$ (길고 얇은 띠에서 등호에 접근)
변 $a$ 인 정사각형, $H_0^1$, 2차원	$\lambda_1 = 2\pi^2/a^2$ 이므로 $C = a/(\pi\sqrt2)$
지름 $d$ 인 볼록 영역, 평균 0	$C \le d/\pi$ (페인-와인버거)

세 번째 줄이 실무에서 가장 유용하다. 페인-와인버거 부등식은 볼록 영역에 대해 $\mu_2 \ge \pi^2/d^2$ 를 주며, 볼록성만 요구하고 모양에는 무관하다. 유한요소 요소 하나(삼각형·사면체)는 볼록이므로, 국소 푸앵카레 상수를 요소 지름 $h_K$ 의 $1/\pi$ 배로 바로 잡을 수 있다. 아래에서 볼 사후 오차 추정자의 $h_K$ 가중치가 여기서 나온다.

거꾸로 상수가 나빠지는 경우도 알아 두어야 한다. 부피를 고정하면 파버-크란 부등식에 의해 공이 $\lambda_1$ 을 최소화하므로, 같은 부피에서 푸앵카레 상수가 가장 큰(가장 불리한) 영역은 공이다. 얇고 길쭉한 영역은 디리클레 조건 아래에서는 오히려 유리하다 — 한 방향으로 얇으면 $d$ 가 작기 때문. 반면 평균 0 판에서는 정반대라, 아령 모양(dumbbell)의 목이 가늘어질수록 $\mu_2 \to 0$ 이고 상수가 폭발한다. 질량이 두 방으로 나뉘어 각각 다른 값을 갖는 함수가 기울기 비용을 거의 안 내고 만들어지기 때문이며, 이건 등주부등식·체거 상수와 직결되는 현상이다.

4. 강제성과 유일성 — 타원형 문제의 시동키[편집]

푸앵카레 부등식이 실제로 하는 일은 딱 하나다. 에너지 반노름을 진짜 노름으로 승격시킨다.

$$\|u\|_{H^1}^2 = \|u\|_{L^2}^2 + \|\nabla u\|_{L^2}^2 \le (1 + C_\Omega^2)\,\|\nabla u\|_{L^2}^2$$

이므로, 푸아송 문제 $-\nabla\cdot(\kappa\nabla u) = f$ 의 약형식에 나오는 쌍선형형식 $a(u,v) = \int \kappa\nabla u\cdot\nabla v$ 에 대해 ($\kappa \ge \kappa_0 > 0$ 일 때)

$$a(u,u) \;\ge\; \kappa_0 \|\nabla u\|_{L^2}^2 \;\ge\; \frac{\kappa_0}{1+C_\Omega^2}\,\|u\|_{H^1}^2$$

가 성립한다. 이것이 락스-밀그램 정리가 요구하는 강제성(coercivity)이고, 강제성 상수는 $\alpha = \kappa_0/(1+C_\Omega^2)$ 다. 여기서부터 도미노가 넘어간다.

유일성 증명이 두 줄로 끝난다. 라플라스 방정식 $-\Delta u = f$, $u|_{\partial\Omega} = g$ 의 해가 둘이라 하고 차 $w$ 를 잡으면 $w \in H_0^1$ 이고 $\int|\nabla w|^2 = 0$ 이다. 푸앵카레로 $\|w\|_{L^2} \le C \cdot 0 = 0$, 따라서 $w = 0$. 푸앵카레 없이 여기서 얻는 것은 ”$w$ 가 상수”뿐이고, 경계조건으로 그 상수를 죽이려면 어차피 자취 정리가 필요하다. 라플라스 방정식·포아송 방정식 교과서가 “에너지 방법에 의한 유일성”이라 부르는 논증의 정체가 이것이다.

순수 노이만 문제의 정체도 여기서 설명된다. $-\Delta u = f$ 에 경계 전체가 노이만이면 해가 상수만큼 결정되지 않고, 이산화하면 강성행렬이 특이해져 상수벡터가 영공간에 앉는다. 반복 솔버가 발산하는 그 흔한 사고의 원인이 “이 문제에는 프리드리히스가 없다”는 것이다. 처방은 정석적이다 — 적합성 조건 $\int_\Omega f = 0$ 을 확인하고, 평균 0 부분공간으로 상(quotient)을 취하고, 그 위에서 푸앵카레-비르팅거로 강제성을 회복한다. 압력이 상수만큼 미정인 비압축성 유동에서 압력 기준점을 하나 못 박거나 평균을 0으로 고정하는 것도 정확히 같은 조작이다.

혼합 경계조건에서도 살아남는다. 경계의 일부 $\Gamma_D \subset \partial\Omega$ 가 양의 측도를 갖고 거기서만 $u=0$ 이면 프리드리히스가 그대로 성립한다(상수는 $\Gamma_D$ 가 작아질수록 커진다). 실무의 구조·열해석 모델이 대개 이 꼴이고, “고정 경계를 한 점에만 걸었더니 강성행렬이 특이하더라”는 사고가 점은 측도 0이라는 사실의 직접적 귀결이다.

5. 유한요소 오차해석에서의 역할[편집]

강제성이 확보되면 세아(Céa) 보조정리가 나온다.

$$\|u - u_h\|_{H^1} \;\le\; \frac{M}{\alpha}\, \inf_{v_h \in V_h}\|u - v_h\|_{H^1}, \qquad \alpha = \frac{\kappa_0}{1+C_\Omega^2}$$

즉 푸앵카레 상수가 곧 오차 상수의 분모다. 영역이 커지거나(대형 구조물) 이방성이 심해지면 $C_\Omega$ 가 커지고 $M/\alpha$ 가 함께 커진다. 이론적 상수라 실무 수치와 직결되지는 않지만, “같은 격자인데 도메인만 키웠더니 오차 상수가 나빠진다”는 관측의 근거는 여기에 있다.

더 직접적인 등장은 두 곳이다.

국소 근사와 브램블-힐버트. $h^k$ 수렴 차수를 만드는 브램블-힐버트(데니-리옹) 보조정리는 사실상 다항식으로 나눈 몫공간 위의 푸앵카레 부등식이다. “함수에서 그 평균을 뺀 것은 기울기로 통제된다”를 “함수에서 그 $k$ 차 다항식 근사를 뺀 것은 $k+1$ 계 도함수로 통제된다”로 일반화한 것이고, 요소 지름 $h_K$ 의 거듭제곱이 붙는 방식이 앞서 본 스케일링 $C_{\lambda\Omega} = \lambda C_\Omega$ 그대로다. 참조요소에서 부등식을 한 번 증명하고 아핀 사상으로 옮기면 $h_K$ 승수가 자동으로 따라온다.

사후 오차 추정자의 가중치. 잔차형 추정자

$$\eta_K^2 = h_K^2 \,\|f + \nabla\!\cdot\!(\kappa\nabla u_h)\|_{L^2(K)}^2 \;+\; \tfrac{1}{2}\sum_{F \subset \partial K} h_F \,\bigl\|\,[\kappa\nabla u_h\cdot n\,]\,\bigr\|_{L^2(F)}^2$$

에서 $h_K$ 와 $h_F$ 라는 가중치가 왜 하필 그 거듭제곱인지를 설명하는 것이 국소 푸앵카레 부등식이다. 잔차와 시험함수를 짝지을 때 클레망(Clément) 보간의 국소 근사오차를 쓰는데, 그 추정이 요소 패치 위의 푸앵카레 부등식이고 상수가 $h_K/\pi$ 급이다(요소가 볼록이라 페인-와인버거를 그대로 쓸 수 있다). 적응 격자 세분화의 세분 기준이 요소마다 $\eta_K$ 를 비교하는 방식으로 굴러가는 이상, 이 가중치가 틀리면 세분 판정 자체가 틀린다.

덤으로 조건수의 $h^{-2}$ 법칙도 푸앵카레의 몫이다. 강성행렬 $A$ 의 최소 고유값 하한이 푸앵카레(정확히는 강제성)에서, 최대 고유값 상한이 역부등식(inverse inequality)에서 나오고, 둘을 나누면 $\kappa(A) = O(h^{-2})$ 다. 다중격자법이나 영역 분할법의 수렴률 해석에서 부분영역별 푸앵카레 상수가 명시적으로 등장하는 것도 같은 이유이며, 이때는 부분영역 지름 $H$ 가 상수에 들어와 $H/h$ 꼴의 인자가 만들어진다.

6. 이산 푸앵카레 부등식[편집]

여기서부터가 수치해석 고유의 이야기다. 적합(conforming) 유한요소는 $V_h \subset H_0^1$ 이므로 연속판 부등식을 그냥 물려받는다. 문제는 이산 공간이 $H^1$ 안에 없는 방법들 — 유한체적법과 불연속 갤러킨법 — 이다. 조각별 상수함수는 요소 내부에서 기울기가 0인데 $L^2$ 노름은 0이 아니므로, 요소별 기울기만으로 만든 반노름에 대한 푸앵카레 부등식은 그냥 거짓이다.

처방은 하나다. 점프를 반노름에 집어넣는다.

• 유한체적. 셀 평균 $u_K$ 로 이루어진 조각별 상수 공간에 이산 $H^1$ 반노름 $$|u_h|_{1,h}^2 = \sum_{\sigma = K|L} \frac{|\sigma|}{d_{KL}} \,\bigl(u_K - u_L\bigr)^2$$ 을 주면 $\|u_h\|_{L^2(\Omega)} \le C\,|u_h|_{1,h}$ 가 성립하고, 결정적으로 $C$ 는 격자에 무관하고 $\operatorname{diam}(\Omega)$ 에만 의존한다. 에마르-갈루에-에르뱅의 유한체적 수렴 이론이 이 부등식 위에 서 있으며, 상수가 격자와 무관하다는 것이 곧 $h \to 0$ 에서의 균등 안정성이다.
• DG. 조각별 다항식 공간 위의 깨진(broken) 푸앵카레 부등식은 $$\|v\|_{L^2(\Omega)}^2 \;\le\; C\Bigl( \sum_K \|\nabla v\|_{L^2(K)}^2 \;+\; \sum_{F} h_F^{-1}\,\bigl\|[v]\bigr\|_{L^2(F)}^2 \Bigr)$$ 꼴이다.¹ 오른쪽 두 번째 항이 바로 IPDG의 벌점항이고, 벌점 상수 $\sigma$ 가 충분히 커야 강제성이 나온다는 DG 문서의 그 조건은 “깨진 푸앵카레 부등식을 쓸 수 있을 만큼 점프에 값을 매겨라”는 요구다. 벌점을 빼면 조각별 상수가 커널에 앉아 시스템이 특이해진다 — 순수 노이만 문제에서 상수함수가 앉는 것과 구조가 똑같다.

같은 이야기의 이산 극한판이 그래프 라플라시안이다. 정점 값 벡터의 평균을 0으로 두면 $\sum_v x_v^2 \le C \sum_{(u,v)\in E}(x_u-x_v)^2$ 꼴이 되고, 최적 상수는 $1/\lambda_2$ 즉 스펙트럼 갭의 역수다. 그래프가 병목으로 두 덩어리로 나뉘려 할수록 $\lambda_2$ 가 작아지고 상수가 커지는데, 이는 연속판의 아령 영역과 정확히 같은 현상이며 양쪽 다 체거 부등식이 정량화한다. 스펙트럴 클러스터링, 다중격자법의 대수적 조립(AMG), 마르코프 연쇄의 혼합시간이 전부 이 한 부등식의 변주다.

7. 성립하지 않는 경우[편집]

부등식이 깨지는 조건을 아는 것이 상수를 아는 것만큼 중요하다.

• 모든 방향으로 무한한 영역. $\Omega = \mathbb{R}^n$ 에서는 $H_0^1$ 판도 거짓이다($n \ge 3$ 에서는 대신 소볼레프 부등식 $\|u\|_{L^{2^*}} \le C\|\nabla u\|_{L^2}$ 이 성립하는데, 좌변 지수가 다르다는 점이 핵심이다). 외부 유동·무한 영역 문제에서 에너지 방법이 곧바로 안 통하고 가중 소볼레프 공간이나 무한요소가 등장하는 이유다.
• 비연결 영역(평균 0 판). 위에서 본 대로 반례가 즉시 나온다.
• 첨점(cusp)·병리적 경계. 외부 첨점을 가진 영역이나 “방과 복도(rooms and corridors)” 구조에서는 렐리히-콘드라쇼프 콤팩트성 자체가 깨져 평균 0 판이 무너진다. 립시츠 경계 가정이 정리 진술에 늘 붙어 있는 것이 장식이 아니다.
• 측도 0인 곳에만 건 경계조건. 2·3차원에서 한 점(또는 3차원에서 한 선)에만 디리클레를 걸면 $H^1$ 함수가 그 조건을 “느끼지” 못한다. 수학적으로는 부등식이 없고, 수치적으로는 강성행렬이 거의 특이해져 반복 솔버가 수렴하지 않는다.

한 줄로 요약하면 — 푸앵카레 부등식은 “이 문제에 유일한 해가 있는가”를 묻는 가장 값싼 검사지다. 부등식이 성립하는 조건을 못 채우는 정식화는 이산화 전에 이미 병들어 있고, 격자를 아무리 곱게 깔아도 낫지 않는다.² 편미분방정식 문제를 새로 세울 때 “상수함수(또는 강체운동 모드)가 어디서 죽는가”를 먼저 확인하는 습관이 사고를 절반으로 줄인다.³

8. 관련 문서[편집]

• 소볼레프 공간 · 르베그 적분 · 변분법
• 락스-밀그램 정리 · 갤러킨 방법 · 유한요소법
• 유한체적법 · 불연속 갤러킨법 · 사후 오차 추정
• 고유값 문제 · 레일리 몫 · 등주부등식
• 라플라스 방정식 · 포아송 방정식 · 편미분방정식
• 조건수 · 다중격자법 · 영역 분할법 · 적응 격자 세분화

9. Footnotes[편집]

1. 면 점프항의 가중치가 왜 하필 $h_F^{-1}$ 인지도 스케일링으로 확인된다. 요소를 $\lambda$ 배 키우면 $\|\nabla v\|_{L^2(K)}^2$ 은 $\lambda^{n-2}$ 배, 면 위의 $\|[v]\|_{L^2(F)}^2$ 은 $\lambda^{n-1}$ 배가 되므로 둘을 같은 차원으로 맞추려면 뒤쪽을 $h_F$ 로 한 번 나눠야 한다. 벌점 상수가 $\sigma/h_F$ 꼴로 들어가는 것은 취향이 아니라 차원 맞추기이고, 여기에 $p^2$ 이 곱해지는 것은 다항식 차수에 대한 역부등식에서 온다. ↩

2. 구조해석에서 이 병은 강체운동 모드라는 이름으로 매일 나타난다. 구속을 충분히 걸지 않은 3차원 구조물의 강성행렬은 병진 3 + 회전 3 = 6차원 영공간을 갖고, 솔버는 “pivot is zero” 또는 “singular matrix”를 뱉는다. 관성 릴리프(inertia relief)나 약한 스프링을 붙여 넘기는 관행이 있는데, 이건 문제를 푼 게 아니라 상수함수를 인공적으로 죽여 푸앵카레 부등식을 억지로 성립시킨 것이다. 결과를 볼 때 그 사실을 기억해야 한다. ↩

3. 이름의 주인 앙리 푸앵카레는 이 부등식을 1890년대 디리클레 문제 연구 중에 썼고, 쿠르트 프리드리히스는 1927년 판(경계 0 조건)을 정식화했다. 두 사람 이름이 붙은 부등식이 사실상 같은 물건인데 문헌마다 다르게 불리는 통에, 논문에서 “Poincaré–Friedrichs inequality”라고 붙여 쓰는 절충안이 굳어졌다. 이 바닥에서 정치적 해법이 표기법으로 나타나는 드문 사례. ↩