유한차분법

오류 제보
수치해석 이산화 기법 마지막 수정: 2026-10-19 19:00:00
○ 공개 검수 대기 AI 보조 작성 표시 기준 보기

1. 개요

유한차분법
Finite Difference Method
약칭FDM
분야수치해석 · 이산화 기법
핵심 아이디어테일러 급수 기반 미분 연산자 근사
주 사용처DNS, 학술 연구, 단순 형상 고정밀 해석

유한차분법(Finite Difference Method, FDM)은 미분 연산자를 격자점 값들의 차분 몫(difference quotient)으로 직접 근사하여 편미분방정식을 대수방정식으로 바꾸는 이산화 기법이다. 세 가지 주요 이산화 기법(유한체적법, 유한요소법, FDM) 중 가장 역사가 길고 개념적으로 가장 단순하다. 미분의 정의에서 극한 기호만 떼면 그대로 차분식이 되기 때문에, 수치해석 교과서가 예외 없이 FDM부터 시작하는 데는 이유가 있다.1 산업용 전산유체역학의 주류 자리는 유한체적법에 내주었지만, 고차 정확도가 생명인 DNS와 학술 연구에서는 여전히 현역 주전이다.

2. 테일러 급수와 차분식 유도

FDM의 모든 것은 테일러 급수 전개에서 나온다. 격자 간격이 Δx\Delta x일 때 이웃 격자점의 값은

ϕi+1=ϕi+Δx∂ϕ∂x∣i+Δx22∂2ϕ∂x2∣i+Δx36∂3ϕ∂x3∣i+⋯\phi_{i+1} = \phi_i + \Delta x \left.\frac{\partial \phi}{\partial x}\right|_i + \frac{\Delta x^2}{2}\left.\frac{\partial^2 \phi}{\partial x^2}\right|_i + \frac{\Delta x^3}{6}\left.\frac{\partial^3 \phi}{\partial x^3}\right|_i + \cdots

로 전개된다. 이 식을 1계 미분에 대해 정리하면 전방차분(forward difference)을 얻는다.

∂ϕ∂x∣i=ϕi+1−ϕiΔx+O(Δx)\left.\frac{\partial \phi}{\partial x}\right|_i = \frac{\phi_{i+1} - \phi_i}{\Delta x} + \mathcal{O}(\Delta x)

ϕi−1\phi_{i-1}을 쓰면 후방차분(backward difference)이 되고 역시 1차 정확도다. 두 전개식을 빼면 짝수 차 항이 상쇄되면서 2차 정확도의 중앙차분(central difference)이 나온다.

∂ϕ∂x∣i=ϕi+1−ϕi−12Δx+O(Δx2)\left.\frac{\partial \phi}{\partial x}\right|_i = \frac{\phi_{i+1} - \phi_{i-1}}{2\Delta x} + \mathcal{O}(\Delta x^2)

2계 미분의 표준 중앙차분도 같은 방식으로 유도된다.

∂2ϕ∂x2∣i=ϕi+1−2ϕi+ϕi−1Δx2+O(Δx2)\left.\frac{\partial^2 \phi}{\partial x^2}\right|_i = \frac{\phi_{i+1} - 2\phi_i + \phi_{i-1}}{\Delta x^2} + \mathcal{O}(\Delta x^2)

여기서 O(Δxp)\mathcal{O}(\Delta x^p)가 절단오차(truncation error)의 차수다. pp차 기법이라면 격자를 절반으로 줄일 때 오차가 2p2^p분의 1로 감소한다. 차수만 외우고 넘어가기 쉽지만, 테일러 전개를 한 항 더 붙들고 있으면 선행 오차항의 계수까지 그대로 나온다. 위 두 중앙차분의 경우

ϕi+1−ϕi−12Δx=∂ϕ∂x∣i+Δx26∂3ϕ∂x3∣i+O(Δx4)\frac{\phi_{i+1} - \phi_{i-1}}{2\Delta x} = \left.\frac{\partial \phi}{\partial x}\right|_i + \frac{\Delta x^2}{6}\left.\frac{\partial^3 \phi}{\partial x^3}\right|_i + \mathcal{O}(\Delta x^4) ϕi+1−2ϕi+ϕi−1Δx2=∂2ϕ∂x2∣i+Δx212∂4ϕ∂x4∣i+O(Δx4)\frac{\phi_{i+1} - 2\phi_i + \phi_{i-1}}{\Delta x^2} = \left.\frac{\partial^2 \phi}{\partial x^2}\right|_i + \frac{\Delta x^2}{12}\left.\frac{\partial^4 \phi}{\partial x^4}\right|_i + \mathcal{O}(\Delta x^4)

이고, 두 식 모두 홀수 차 항이 대칭성으로 상쇄되기 때문에 오차가 Δx2\Delta x^2부터 시작한다. 이 선행항을 원래 방정식에 되돌려 넣은 것이 수정 방정식(modified equation)이며, “이 도식이 실제로 푸는 방정식”을 보여준다 — 1계 미분의 선행 오차가 홀수 차 미분이면 수치확산이 아니라 분산(dispersion)으로 나타나는 이유가 여기서 보인다. 격자점을 더 많이 동원하면(넓은 스텐실) 4차, 6차 등 임의의 고차 도식을 체계적으로 만들 수 있으며, 이 확장 용이성이 FDM의 대표 강점이다.

ϕ=sin⁡x\phi = \sin x 를 x=π/2x = \pi/2 에서 2계 중앙차분하면 차수가 말로 끝나지 않는다. 참값은 ϕ′′=−1\phi'' = -1 이다. Δx=0.2\Delta x = 0.2 이면

sin⁡(π/2+0.2)−2sin⁡(π/2)+sin⁡(π/2−0.2)0.22=−0.99667.\frac{\sin(\pi/2+0.2) - 2\sin(\pi/2) + \sin(\pi/2-0.2)}{0.2^2} = -0.99667.

오차는 0.0033290.003329 다. 선행항 (Δx2/12)ϕ′′′′(\Delta x^2/12)\phi'''' 은 0.0033330.003333 이라, 남은 차이는 그다음 항이다. Δx=0.1\Delta x = 0.1 이면 오차는 0.0008330.000833 으로 거의 정확히 4분의 1이 된다. 격자를 반으로 줄였는데 오차가 4분의 1이 아니면, 그 구간은 2차가 아니다.

이론 차수를 믿지 말고 실제로 재보는 것이 검증 및 확인의 기본기다. 격자를 rr배 조밀하게 한 세 해의 오차비로 관측 수렴 차수를 뽑는 절차는 격자 수렴 지수가 다루고, 엄밀해가 없는 방정식에서 오차를 정의하려면 제조해를 쓴다.

3. 일관성 · 안정성 · 수렴성

차분식을 만들었다고 끝이 아니다. 수치해가 엄밀해에 다가간다는 보장이 필요하며, 이는 세 가지 개념으로 정리된다.

  • 일관성(consistency): Δx,Δt→0\Delta x, \Delta t \to 0일 때 차분방정식이 원래 미분방정식으로 수렴하는가. 절단오차가 0으로 가면 성립.
  • 안정성(stability): 계산 도중 발생한 오차(반올림 오차 포함)가 시간 전진 과정에서 증폭되지 않는가.
  • 수렴성(convergence): 격자를 조밀하게 할수록 수치해가 엄밀해에 다가가는가. 정작 우리가 원하는 것.

이 셋을 묶는 것이 락스 동치 정리(Lax equivalence theorem)다: 적절히 제기된(well-posed) 선형 초기값 문제에서, 일관성을 갖춘 기법이 수렴할 필요충분조건은 안정성이다. 수렴성을 직접 증명하는 것은 어렵지만 일관성과 안정성은 각각 따로 확인할 수 있으므로, 이 정리는 “어려운 문제 하나”를 “할 만한 문제 둘”로 바꿔주는 수치해석의 대들보다.2

안정성 판별의 표준 도구는 폰 노이만 안정성 해석으로, 오차를 푸리에 모드 eikxe^{\mathrm{i}kx}로 분해한 뒤 각 모드의 증폭 인자 ∣G∣|G|가 1 이하인지 검사한다. 명시적 시간 전진 기법에 이 해석을 적용하면 그 유명한 CFL 조건이 안정성 한계로 튀어나온다. 락스 동치 정리가 선형 문제에 대한 정리라는 점은 계속 기억해야 한다 — 쿠랑·프리드리히스·레비의 1928년 원논문도, 락스·리히트마이어의 1956년 정리도 전제는 선형이다.

4. 구조 격자의 필요성

FDM의 아킬레스건. 차분식이 (i,j,k)(i, j, k) 인덱스로 이웃을 지목하는 구조라서, 격자점이 규칙적으로 배열된 정렬(구조) 격자가 사실상 필수다. 복잡한 형상 앞에서는 다음과 같은 우회로를 동원한다.

  • 좌표 변환: 물리 공간의 곡면 격자를 계산 공간의 직교 격자로 사상(mapping)한다. 변환 야코비안 항들이 방정식에 주렁주렁 붙는 대가를 치른다.
  • 중첩 격자(overset/Chimera): 형상별로 만든 격자를 겹쳐 놓고 경계에서 보간한다. 보간 과정에서 보존성이 깨질 수 있다.
  • 침입 경계법(immersed boundary): 격자는 직교로 두고 물체를 힘 항으로 표현한다. 경계층 해상도가 고민거리.

전부 훌륭한 기법이지만, “비정렬 격자를 그냥 쓰면 되는” 유한체적법과 비교하면 태생적 제약을 땜질하는 인상을 지우기 어렵다. 산업 CFD 시장에서 FDM이 밀려난 근본 원인이다.

5. 컴팩트 차분 — 좁은 스텐실로 고차 내기

스텐실을 넓혀 차수를 올리는 고전적 방식은 경계 근처에서 쓸 격자점이 모자라고 병렬 통신량이 늘어난다는 대가가 있다. 컴팩트(Padé) 차분은 발상을 뒤집어, 미분값 자체를 미지수로 두고 좌변에도 이웃 미분값을 섞는다.

α ϕi−1′+ϕi′+α ϕi+1′=a ϕi+1−ϕi−12Δx+b ϕi+2−ϕi−24Δx\alpha\,\phi'_{i-1} + \phi'_i + \alpha\,\phi'_{i+1} = a\,\frac{\phi_{i+1} - \phi_{i-1}}{2\Delta x} + b\,\frac{\phi_{i+2} - \phi_{i-2}}{4\Delta x}

테일러 전개를 차수별로 맞추면 계수가 a=23(α+2)a = \tfrac{2}{3}(\alpha + 2), b=13(4α−1)b = \tfrac{1}{3}(4\alpha - 1)로 묶인다. α=14, b=0\alpha = \tfrac{1}{4},\,b=0이면 3점 스텐실만으로 4차(a=3/2a = 3/2), α=13\alpha = \tfrac{1}{3}이면 5점 스텐실로 6차(a=14/9, b=1/9a = 14/9,\ b = 1/9)다. 명시적 도식이 같은 차수를 내려면 각각 5점·7점이 필요하니, 스텐실 두 칸을 아낀 셈이다. 대가는 미분값을 얻을 때마다 삼중대각 연립방정식을 풀어야 한다는 것 — 토마스 알고리즘이면 O(N)O(N)이라 비싸지는 않지만, 방향마다 전 격자선을 훑어야 해서 영역 분할 병렬화가 까다로워진다.

값진 것은 차수가 아니라 파수 분해능이다. 수정 파수 k′Δxk'\Delta x가 실제 kΔxk\Delta x를 얼마나 오래 따라가는지로 재면, 6차 컴팩트 도식은 같은 격자로 명시적 6차보다 훨씬 높은 파수까지 위상 오차 없이 버틴다. “spectral-like resolution”이라는 이름이 여기서 나왔다.3

6. DNS와 학술 연구에서의 위상

산업을 내준 대신 FDM은 연구의 최전선을 지키고 있다. 난류의 모든 스케일을 모델 없이 직접 해상하는 DNS는 격자당 정확도가 곧 생명인데, 여기서 FDM의 고차 확장성이 빛난다. DNS의 계산량은 레이놀즈수에 대해 살인적으로 증가하므로4, 같은 격자로 더 정확한 답을 주는 기법이 곧 더 높은 레이놀즈수를 의미한다. 채널 유동, 평판 경계층처럼 형상이 단순한 정준 문제(canonical problem)가 DNS의 주 무대라는 점도 구조 격자 제약을 무의미하게 만든다.

다만 정준 DNS의 주력이 곧 FDM이라고 말하면 과장이다. 주기 방향에 푸리에, 벽수직 방향에 체비셰프를 쓰는 스펙트럴 방법 이 채널·평판 DNS의 사실상 표준이고(김·모인·모저의 1987년 Reτ=180Re_\tau = 180 채널 DNS부터 그렇다), FDM·컴팩트 차분은 스펙트럴이 쓸 수 없는 곳에서 주전이 된다 — 압축성 유동, 비주기·복잡 형상, 가변 물성, 충격파가 섞인 문제. 실제로 널리 쓰이는 코드 상당수는 방향별로 둘을 섞은 혼합 구성이다.

7. 유한체적법과의 비교

항목유한차분법 (FDM)유한체적법 (FVM)
근사 대상미분 연산자 자체검사체적 적분 보존식
보존성자동 보장되지 않음국소 보존 보장
격자구조 격자 필수비정렬·다면체 가능
고차 확장체계적이고 용이2차 이상은 번거로움
주 무대DNS, 학술 연구산업용 CFD

요약하면 FDM은 “단순한 형상에서 극한의 정확도”, FVM은 “복잡한 형상에서 견고한 보존성”으로 역할이 분담되어 있다. 어느 쪽이 우월하냐는 질문은 스포츠카와 트럭 중 뭐가 좋은 차냐는 질문과 같다.

8. 이 다음 30분

  • 폰 노이만 안정성 해석 — 차수가 맞아도 시간 전진에서 오차가 커지는 조건.
  • CFL 조건 — 64칸 펄스가 C=1.2C = 1.2 에서 30스텝 만에 2500 을 넘는 계산.
  • 검증 및 확인 — 위의 사인 테스트가 코드 검증의 가장 작은 형태인 이유.
  • 제조해 — 엄밀해가 없는 방정식에서 오차를 정의하는 방법.
  • 격자 수렴 지수 — 격자 세 벌로 관측 차수를 뽑는 절차.

9. 관련 문서

10. Footnotes

  1. 오일러가 1768년 『적분학 원론』에서 이미 차분 근사로 미분방정식을 풀었으니, 컴퓨터보다 약 180년 선배인 셈이다. 물론 오일러는 격자 수렴성 연구를 손으로 했다. ↩

  2. Lax, P. D. & Richtmyer, R. D. (1956), Comm. Pure Appl. Math. 9, 267-293. 단, 어디까지나 “선형” 문제에 대한 정리다. 나비에-스토크스 같은 비선형 방정식에는 엄밀히 적용되지 않지만, 다들 마음의 평화를 위해 성립하는 셈 치고 쓴다. ↩

  3. 파수 분해능이라는 말이 거창해 보이지만 뜻은 단순하다. 차분 연산자를 푸리에 모드에 먹여 보면 ik\mathrm{i}k 대신 ik′\mathrm{i}k'가 나오는데, 이 k′k'가 진짜 kk를 얼마나 따라오는지가 전부다. 2차 중앙차분은 kΔxk\Delta x가 1만 넘어가도 슬슬 게을러지기 시작한다 — 격자 8칸짜리 파동을 “대충 이 정도 기울기?” 하고 뭉개는 셈이다. ↩

  4. 지수가 문헌마다 다른데, 어떤 난류를 재느냐가 다르기 때문이다. 균질 등방성 난류를 콜모고로프 스케일까지 해상하면 격자점 수가 대략 Re9/4Re^{9/4}(약 2.25제곱)로 늘고, 평판 경계층처럼 벽이 있는 유동은 Choi & Moin (2012)의 갱신치로 N∼ReLx37/14N \sim Re_{L_x}^{37/14}(약 2.64제곱)다. 후자 기준이면 레이놀즈수를 10배 올릴 때 격자점만 102.64≈44010^{2.64} \approx 440배 — 시간 스텝까지 줄어드니 총 비용은 그보다 더 나빠진다. 같은 논문이 벽해상 LES 를 ReLx13/7Re_{L_x}^{13/7}, 벽모형 LES 를 ReLxRe_{L_x}로 놓는다. DNS 연구자들이 슈퍼컴퓨터 지원사업 공고에 민감한 이유. ↩

검증 자료

이 문서의 주요 출처

  1. Lax & Richtmyer (1956), Survey of the stability of linear finite difference equations, Comm. Pure Appl. Math. 9, 267-293
  2. Courant, Friedrichs & Lewy (1928), Über die partiellen Differenzengleichungen der mathematischen Physik, Math. Ann. 100, 32-74
  3. Lele (1992), Compact finite difference schemes with spectral-like resolution, J. Comput. Phys. 103, 16-42
  4. Choi & Moin (2012), Grid-point requirements for large eddy simulation: Chapman's estimates revisited, Phys. Fluids 24, 011702
  5. Kim, Moin & Moser (1987), Turbulence statistics in fully developed channel flow at low Reynolds number, J. Fluid Mech. 177, 133-166
  6. Euler (1768), Institutionum calculi integralis, vol. I (Sectio II, Caput VII — 오일러 방법의 원전)
이 문서가 문제를 푸는 데 도움이 됐나요? 한 번의 클릭이 다음 검수 순서를 정합니다.

의견을 반영할게요. 감사합니다.

다음으로 읽기 나비에-스토크스 방정식 약 5분 더 탐색 →