SIMWIKI DOC // 수치해석

유한차분법

편집 역사 토론
수치해석 이산화 기법 CONVERGED last update 2026-06-28 22:05:31
유한차분법
Finite Difference Method
약칭FDM
분야수치해석 · 이산화 기법
핵심 아이디어테일러 급수 기반 미분 연산자 근사
주 사용처DNS, 학술 연구, 단순 형상 고정밀 해석

1. 개요[편집]

유한차분법(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로 감소한다. 격자점을 더 많이 동원하면(넓은 스텐실) 4차, 6차 등 임의의 고차 도식을 체계적으로 만들 수 있으며, 이 확장 용이성이 FDM의 대표 강점이다.

3. 일관성 · 안정성 · 수렴성[편집]

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

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

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

안정성 판별의 표준 도구는 von Neumann 안정성 해석으로, 오차를 푸리에 모드 eikxe^{\mathrm{i}kx}로 분해한 뒤 각 모드의 증폭 인자 ∣G∣|G|가 1 이하인지 검사한다. 명시적 시간 전진 기법에 이 해석을 적용하면 그 유명한 CFL 조건이 안정성 한계로 튀어나온다.

4. 구조 격자의 필요성[편집]

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

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

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

5. DNS와 학술 연구에서의 위상[편집]

산업을 내준 대신 FDM은 연구의 최전선을 지키고 있다. 난류의 모든 스케일을 모델 없이 직접 해상하는 DNS는 격자당 정확도가 곧 생명인데, 여기서 FDM의 고차 확장성이 빛난다. 특히 고차 컴팩트 도식(compact scheme)은 좁은 스텐실로 스펙트럴 방법에 근접하는 파수 분해능을 달성한다. DNS의 계산량은 레이놀즈수에 대해 살인적으로 증가하므로3, 같은 격자로 더 정확한 답을 주는 기법이 곧 더 높은 레이놀즈수를 의미한다. 채널 유동, 평판 경계층처럼 형상이 단순한 정준 문제(canonical problem)가 DNS의 주 무대라는 점도 구조 격자 제약을 무의미하게 만든다. 실제로 널리 인용되는 난류 DNS 데이터베이스 상당수가 FDM 혹은 스펙트럴-FDM 혼합 코드의 산물이다.

6. 유한체적법과의 비교[편집]

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

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

7. 관련 문서[편집]

8. Footnotes[편집]

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

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

  3. 격자점 수 기준으로 대략 레이놀즈수의 37/14제곱(약 2.64제곱)으로 증가한다는 것이 고전적 추정이다. 레이놀즈수를 10배 올리려면 컴퓨터를 400배쯤 키워야 한다는 뜻으로, DNS 연구자들이 슈퍼컴퓨터 지원사업 공고에 민감한 이유다. ↩