FDTD

전자기해석 수치해석 마지막 수정: 2026-07-03 14:27:05

시간영역 유한차분법 Finite-Difference Time-Domain
약칭	FDTD
제안	케인 이(Kane S. Yee), 1966
지배 방정식	맥스웰 방정식 (시간영역)
기반 기법	유한차분법
안정 조건	쿠랑 한계 (CFL 조건)

1. 개요[편집]

펄스 한 방 쏘고, 나머지는 컴퓨터가 알아서 한다.

FDTD(Finite-Difference Time-Domain, 시간영역 유한차분법)는 맥스웰 방정식의 회전(curl) 방정식을 시간과 공간 모두에서 차분으로 근사하여, 전자기장이 시간에 따라 어떻게 퍼져나가는지를 직접 계산하는 유한차분법 기반의 수치 기법이다. 1966년 케인 이(Kane S. Yee)가 제안한 격자 배치 아이디어에서 출발했으며¹, 그래서 지금도 그 격자를 Yee 격자라 부른다.

주파수 영역 기법이 특정 주파수 하나하나를 따로 푸는 것과 달리, FDTD는 시간 파형을 그대로 굴린다. 그래서 넓은 대역의 응답을 시뮬레이션 한 번으로 뽑아낼 수 있고, 전자기파가 구조물을 통과하고 반사되는 과정을 애니메이션처럼 눈으로 볼 수 있다는 것이 큰 매력이다.

2. Yee 격자 — 엇갈림의 미학[편집]

FDTD의 핵심 발명품은 전기장 $\mathbf{E}$ 와 자기장 $\mathbf{H}$ 를 같은 지점에 두지 않고 공간적으로 반 칸 엇갈리게(staggered) 배치하는 Yee 격자다. 전기장 성분은 격자 셀의 모서리(edge) 중앙에, 자기장 성분은 면(face) 중앙에 놓이는 식으로, 각 전기장을 자기장 성분들이 고리처럼 둘러싸고 그 반대도 성립한다.

이 배치가 왜 천재적인가 하면, 맥스웰 방정식의 회전 연산자가 요구하는 “어떤 장의 회전이 옆에 있는 다른 장을 만든다”는 구조를 격자가 그대로 흉내 내기 때문이다. 덕분에 중심 차분(central difference)이 자연스럽게 성립해 2차 정확도를 얻고, 발산 조건( $\nabla \cdot \mathbf{B} = 0$ 등)이 별도 강제 없이도 저절로 지켜진다.

3. Leapfrog 시간 전진[편집]

시간 방향으로도 $\mathbf{E}$ 와 $\mathbf{H}$ 는 반 스텝 엇갈린 시각에 갱신된다. 이를 도약 개구리(leapfrog) 방식이라 부른다.

\mathbf{E}^{\,n+1} = \mathbf{E}^{\,n} + \frac{\Delta t}{\varepsilon} \, \nabla \times \mathbf{H}^{\,n+1/2}

\mathbf{H}^{\,n+3/2} = \mathbf{H}^{\,n+1/2} - \frac{\Delta t}{\mu} \, \nabla \times \mathbf{E}^{\,n+1}

전기장은 정수 시각에서, 자기장은 반정수 시각에서 갱신되며 서로가 서로를 밟고 뛰어넘는다. 각 스텝이 직전에 계산된 값만으로 다음 값을 구하는 명시적(explicit) 갱신이라, 큰 행렬을 풀 필요 없이 곱셈과 덧셈만 반복하면 된다. 구현이 단순하고 병렬화가 극단적으로 쉬운 이유가 여기에 있다.

4. 안정 조건 (쿠랑 한계)[편집]

공짜 점심은 없다. 명시적 시간 전진의 대가로 FDTD는 시간 간격 $\Delta t$ 를 함부로 키울 수 없다. 파동이 한 스텝 동안 격자 한 칸 이상을 건너뛰면 정보가 새어나가 해가 폭발하기 때문이다. 3차원에서의 안정 조건, 즉 쿠랑 한계(Courant limit)는 다음과 같다.

\Delta t \le \frac{1}{c \sqrt{\dfrac{1}{\Delta x^2} + \dfrac{1}{\Delta y^2} + \dfrac{1}{\Delta z^2}}}

이것은 본질적으로 CFL 조건의 전자기 버전이다. 격자를 절반으로 촘촘하게 만들면 시간 간격도 절반으로 줄여야 하므로, 해상도를 2배 올리면 계산량은 (3차원 격자 8배) × (시간 스텝 2배) = 대략 16배로 뛴다. FDTD 엔지니어가 격자 셀 수를 늘리기 전에 심호흡을 하는 이유다.

5. 흡수 경계조건 (PML)[편집]

무한히 열린 공간에서 벌어지는 방사 문제를 유한한 격자에 담으려면, 격자 끝에 도달한 파동을 반사 없이 삼켜버리는 장치가 필요하다. 이 역할을 하는 것이 완전 정합층(Perfectly Matched Layer, PML)이다.² 1994년 베렝제(Berenger)가 제안한 이 인공 흡수층은 격자 가장자리를 둘러싸 마치 공간이 무한히 이어지는 것처럼 파동을 감쇠시킨다. PML이 등장하기 전의 흡수 경계조건들은 반사가 커서 FDTD의 신뢰도를 갉아먹었는데, PML이 이 문제를 사실상 해결하면서 FDTD의 실용성이 급상승했다.

6. 강점과 약점[편집]

FDTD의 장단점을 정리하면 다음과 같다.

구분	내용
강점	시뮬레이션 한 번으로 광대역 응답 확보 (펄스 → 푸리에 변환)
강점	시간 응답과 과도 현상을 직접 관찰 (시각화에 최적)
강점	알고리즘이 단순하고 병렬화·GPU 가속이 쉬움
약점	정렬 격자라 곡면이 계단 근사(staircasing)됨
약점	격자 분산(numerical dispersion)으로 파동 속도가 미세하게 왜곡
약점	큰 빈 공간까지 전부 격자로 채워야 함 (자유공간 산란은 모멘트법이 유리)

계단 근사는 곡면 안테나나 매끈한 렌즈처럼 형상이 중요한 문제에서 오차의 주범이 되고, 격자 분산은 파동이 격자를 오래 지날수록 위상이 어긋나는 형태로 누적된다. 두 약점 모두 격자를 촘촘하게 하면 완화되지만, 그러면 위에서 본 쿠랑 한계가 계산 비용으로 복수를 해온다.

7. 관련 문서[편집]

Yee, K. S. (1966). “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media”. IEEE Trans. Antennas Propag. 당시엔 크게 주목받지 못하다가 1970~80년대 컴퓨터 성능이 따라오면서 뒤늦게 재발견되어 전자기해석의 주력 기법으로 등극했다. 논문 한 편으로 격자에 이름을 남긴 케이스. ↩
이름부터가 “완벽하게(Perfectly)“라니 자신감이 대단하지만, 실제로는 이산 격자에서 완벽하진 않고 아주 약간의 반사가 남는다. 그래도 이전 기법 대비 반사가 수십 dB 낮아서 실무에서는 충분히 “완벽하다”고 쳐준다. ↩