편미분방정식

수치해석 물리 마지막 수정: 2026-07-14 04:19:03

1. 개요[편집]

편미분방정식 Partial Differential Equation
약칭	PDE
분야	해석학 · 수리물리 · 수치해석
미지수	둘 이상 변수의 함수
2계 분류	타원형 · 포물형 · 쌍곡형
대표 수치기법	유한요소법, 유한차분법, 유한체적법

자연은 미분방정식으로 말하고, 우리는 그걸 못 알아들어서 컴퓨터에 통역을 맡긴다.

편미분방정식(Partial Differential Equation, PDE)은 두 개 이상의 독립변수를 갖는 미지 함수와 그 편도함수들 사이의 관계를 기술하는 방정식이다. 상미분방정식(ODE)이 시간 하나에 대한 변화를 다룬다면, PDE는 공간과 시간에 걸친 장(field)의 변화를 다룬다. 열이 어떻게 퍼지는지, 파동이 어떻게 전파되는지, 정상 상태에서 퍼텐셜이 어떻게 분포하는지 — 물리학의 지배 방정식은 거의 예외 없이 PDE다.

나비에-스토크스 방정식, 맥스웰 방정식, 슈뢰딩거 방정식, 탄성체의 평형 방정식 — 시뮬레이션이라는 분야가 통째로 하는 일이 결국 “어떤 PDE를 어떤 경계 조건 아래에서 푸느냐”로 요약된다. 그래서 이 문서는 심위키 여러 문서가 결국 돌아오는 지배 방정식 허브다. 문제는 극소수의 예쁜 경우를 빼면 손으로는 안 풀린다는 것. 그래서 인류는 공간을 잘게 쪼개 대수방정식으로 바꾸는 이산화(discretization)에 인생을 갈아 넣게 되었다.

2. 계수와 차수[편집]

PDE를 구분하는 첫 번째 잣대는 차수(order)와 선형성(linearity)이다. 차수는 방정식에 등장하는 최고 편도함수의 계수를 뜻한다. 물리에서 압도적으로 자주 나오는 것은 2계 PDE인데, 뉴턴의 운동 법칙이 2계이고 확산·파동·퍼텐셜이 모두 2계 라플라시안( $\nabla^2$ )을 품고 있기 때문이다.

선형 PDE는 미지 함수와 그 도함수들이 1차로만 등장하고 서로 곱해지지 않는 얌전한 부류다. 중첩 원리(superposition)가 성립해서 해를 조각조각 더할 수 있고, 스펙트럴 방법이나 그린 함수 같은 우아한 도구가 먹힌다. 반면 비선형 PDE는 미지 함수끼리 곱해지는 순간부터 지옥이 열린다. 나비에-스토크스 방정식의 대류항 $(\mathbf{u}\cdot\nabla)\mathbf{u}$ 가 바로 그 비선형의 원흉이고, 난류라는 인류 미해결 난제의 뿌리다.¹

3. 2계 PDE의 분류: 타원·포물·쌍곡[편집]

2변수 2계 선형 PDE의 일반형은 다음과 같이 쓸 수 있다.

A\,u_{xx} + 2B\,u_{xy} + C\,u_{yy} + (\text{저계항}) = 0

놀랍게도 이 방정식의 성격은 최고계 계수 $A, B, C$ 가 만드는 판별식(discriminant) $B^2 - AC$ 의 부호 하나로 갈린다. 원뿔곡선의 판별식과 완전히 똑같은 구조라서 이름도 거기서 따왔다.

$B^2 - AC < 0$ → 타원형(elliptic)
$B^2 - AC = 0$ → 포물형(parabolic)
$B^2 - AC > 0$ → 쌍곡형(hyperbolic)

이 분류가 단순한 수학 놀음이 아닌 이유는, 각 유형이 근본적으로 다른 물리와 다른 수치 전략을 요구하기 때문이다. 정보가 어떻게 전파되는지(특성곡선), 어떤 경계 조건이 적절한지, 어떤 이산화가 안정한지가 전부 이 부호에서 결정된다. 유형을 잘못 읽고 엉뚱한 솔버를 붙이면 수렴성은 신에게 맡기게 된다.

4. 세 유형의 물리적 의미[편집]

타원형 — 정상 상태의 균형. 대표 주자는 라플라스 방정식 $\nabla^2 u = 0$ 과 그 소스 버전인 포아송 방정식 $\nabla^2 u = f$ 다. 시간이 사라진 방정식이라 “이미 안정된 상태”를 기술한다. 정전기 퍼텐셜, 정상 열분포, 비압축성 포텐셜 유동이 전부 여기 속한다. 특징은 도메인 안의 한 점 값이 경계 전체와 즉각적으로 얽혀 있다는 것. 그래서 전역적으로 큰 연립방정식을 한 번에 풀어야 하고, 크리로프 부분공간법 같은 반복 해법이 필수다.

포물형 — 확산과 망각. 대표 주자는 열방정식 $u_t = \alpha \nabla^2 u$ 와 대류-확산 방정식이다. 시간이 한 방향으로만 흐르고, 초기의 뾰족한 정보가 시간이 지나며 부드럽게 뭉개진다. 커피가 식고 잉크가 퍼지는 그 비가역적 과정. 정보가 무한 속도로 전파되는(수학적으로) 특성 때문에 안정성 조건이 까다롭고, 시간 적분에 음해법(implicit)을 자주 쓴다.²

쌍곡형 — 파동과 기억. 대표 주자는 파동방정식 $u_{tt} = c^2 \nabla^2 u$ 과 압축성 오일러 방정식이다. 정보가 유한한 속도 $c$ 로 특성곡선(characteristic)을 따라 전파되며, 초기의 파형이 흐트러지지 않고 멀리 실려 간다. 소리·빛·충격파가 여기 속한다. 유한 전파 속도 때문에 시간 스텝이 격자 간격에 묶이는 CFL 조건을 반드시 지켜야 하고, 충격파 근처의 진동을 잡기 위해 유한체적법의 정교한 차분 도식이 동원된다.

정리하면, 같은 2계 방정식이라도 판별식 부호 하나로 정보의 전파 방식(즉각/확산/파동)과 시간의 성격(없음/비가역/가역)이 완전히 갈린다. 이 세 원형(archetype)을 몸으로 익히는 것이 시뮬레이션 입문의 절반이다.

5. 수치 해법과 경계 조건[편집]

해석해가 없으니 결국 이산화다. 대표적인 세 갈래는 유한차분법(격자점에서 도함수를 차분으로 근사), 유한체적법(검사체적의 플럭스 보존), 유한요소법(약형식과 기저함수)이다. 어떤 방법을 쓰든 PDE만으로는 해가 유일하게 정해지지 않는다. **경계 조건**과, 시간 문제라면 초기 조건이 반드시 함께 주어져야 한다.

경계 조건은 크게 값을 지정하는 디리클레(Dirichlet), 기울기(플럭스)를 지정하는 노이만(Neumann), 둘을 섞은 로빈(Robin)으로 나뉜다. 여기서 아주 중요한 규칙 하나. PDE 유형마다 잘 맞는(well-posed) 경계 조건이 다르다는 것이다. 타원형은 도메인 전체를 감싸는 경계값이 필요하고, 쌍곡형·포물형은 시간의 출발선(초기값)과 함께 특정 방향에서의 조건이 어울린다. 이 궁합을 무시하면 해가 존재하지 않거나 무수히 많거나 데이터에 미친 듯이 민감해지는(ill-posed) 참사가 벌어진다.³ 잘 세워진 문제(well-posed problem)의 세 조건 — 존재·유일·안정 — 을 아다마르가 못 박은 이유가 여기 있다.

6. 관련 문서[편집]

3차원 나비에-스토크스의 해가 항상 매끄럽게 존재하는지는 클레이 연구소 밀레니엄 문제(상금 100만 달러)로 아직 미해결이다. 비선형 PDE 앞에서는 인류 전체가 아직 학부생이다. ↩
포물형에서 정보가 “무한 속도”로 퍼진다는 건 수학적 이상화다. 실제 열은 그렇게 안 퍼진다. 그래서 상대론과 충돌하는 이 성질을 고친 쌍곡형 열전도 모델(카타네오 방정식)도 있다. 물리학자들은 언제나 자기 방정식이 마음에 안 든다. ↩
대표적 함정이 쌍곡형인 파동방정식에 타원형처럼 사방을 둘러싼 경계값을 주는 것. 수치적으로 폭발하거나 말도 안 되는 해가 튀어나온다. “일단 돌려” 정신은 아름답지만, 유형과 경계 조건의 궁합만은 돌리기 전에 확인하자. ↩