완전정합층

전자기해석 전자기학 마지막 수정: 2026-07-13 04:14:02

1. 개요[편집]

완전정합층(Perfectly Matched Layer, PML)은 무한히 열린 공간을 유한한 격자로 잘라 쓸 때, 잘린 경계에서 파(wave)를 반사 없이 흡수하도록 설계한 인공 흡수 영역이다. FDTD나 유한요소 전자기 같은 시간영역·주파수영역 해석에서 안테나가 뿜어낸 전파나 산란파는 원래 무한 원방(far field)으로 퍼져 나가야 한다. 그런데 컴퓨터 메모리는 유한하니 격자를 어디선가 잘라야 하고, 아무 조치 없이 자르면 그 벽에서 파가 되튕겨 들어와 해석 도메인을 오염시킨다. PML은 이 벽에 “파를 빨아들이는 스펀지 층”을 덧대어, 마치 계산 영역이 무한히 이어지는 것처럼 속이는 기법이다.¹

핵심은 이름 그대로 완전 정합(perfectly matched)에 있다. 단순한 손실 매질을 갖다 붙이면 매질 경계에서 임피던스가 바뀌어 또 반사가 생긴다. PML은 흡수 매질을 원래 매질과 임피던스가 정확히 일치하도록, 그것도 입사각과 주파수에 무관하게 일치하도록 교묘하게 설계한다. 이론상 연속 방정식 수준에서는 반사 계수가 0이다. 1994년 프랑스의 장피에르 베렝제(Jean-Pierre Bérenger)가 발표한 이 아이디어는 전산전자기학에서 손꼽히는 히트작이 되었다.²

2. 왜 필요한가: 흡수 경계조건의 한계[편집]

PML 이전에도 열린 경계를 흉내 내는 방법은 있었다. 통칭 흡수 경계조건(Absorbing Boundary Condition, ABC)이라 부르는데, 무르(Mur) 조건이나 엔퀴스트-마이다(Engquist-Majda) 조건처럼 경계에서 파동방정식을 한쪽 방향으로만 전파하도록 근사하는 미분 연산자다.

문제는 이런 ABC가 경계에 수직으로 입사하는 파에만 잘 듣는다는 것이다. 비스듬히 들어오는 파나 넓은 대역의 파에 대해서는 반사가 눈에 띄게 남는다. 안테나처럼 사방으로 방사되는 문제에서는 이 각도 의존성이 치명적이다. 반사를 -40 dB 이하로 억누르려면 경계를 해석 대상에서 아주 멀찍이 떨어뜨려야 했고, 그만큼 격자가 낭비됐다.

PML은 발상을 뒤집었다. 경계에서 방정식을 근사하는 대신, 경계 바깥에 물리적으로 파를 감쇠시키는 층을 두되 그 층으로 들어가는 순간의 반사만 0으로 만든다. 층 안에서 파는 지수적으로 죽고, 설령 남은 파가 층 끝 벽에서 되튕겨도 돌아오는 길에 다시 감쇠되므로 이중으로 줄어든다.

3. 좌표 확장: PML의 수학적 심장[편집]

현대적으로 PML을 이해하는 가장 깔끔한 관점은 복소 좌표 확장(complex coordinate stretching)이다. 맥스웰 방정식을 주파수영역에서 쓴 뒤, 경계에 수직한 방향(예: $x$ )의 미분 연산자를 다음처럼 복소수로 늘려버린다.

\frac{\partial}{\partial x} \;\longrightarrow\; \frac{1}{s_x}\frac{\partial}{\partial x}, \qquad s_x = 1 + \frac{\sigma_x(x)}{j\omega \varepsilon_0}

여기서 $\sigma_x$ 는 PML 층 안에서만 0이 아닌 가상의 전도도(conductivity) 프로파일이고, $\omega$ 는 각주파수다. 이렇게 하면 PML 안으로 들어간 파는 $\exp(-\sigma_x x / (c\,\varepsilon_0))$ 꼴의 감쇠 인자를 얻어 지수적으로 줄어든다. 결정적으로, $s_x$ 를 $x$ 방향에만 걸고 접선 방향은 건드리지 않기 때문에 경계에서 임피던스가 보존되어 반사가 생기지 않는다.³

전도도는 층 안쪽 끝에서 0이었다가 바깥으로 갈수록 매끄럽게 커지도록 다항식이나 기하급수 프로파일로 준다. 흔히 쓰는 것은 차수 $m$ 의 다항식 등급(polynomial grading)이다.

\sigma_x(x) = \sigma_{\max}\left(\frac{x}{d}\right)^{m}

$d$ 는 PML 두께, $m$ 은 보통 24를 쓴다. 프로파일을 급격하게 주면 이산 격자에서 수치적 반사가 커지고, 너무 완만하게 주면 두꺼운 층이 필요하니 이 사이에서 타협한다. 실무에서는 층 두께를 파장의 절반한 파장, 셀 수로는 8~16층 정도 두는 것이 국룰이다.

4. 종류: 분할형·단축형·CFS-PML[편집]

베렝제의 원조 PML은 각 전자기장 성분을 두 개로 쪼개 방정식을 다시 쓴 분할장(split-field) 형태였다. 개념은 명쾌하지만 맥스웰 방정식의 원래 형태를 깨뜨려 구현이 지저분하다는 단점이 있었다.

이후 좌표 확장 관점이 정립되면서 단축형 PML(Uniaxial PML, UPML)이 등장했다. 이건 PML을 이방성(anisotropic) 유전율·투자율 텐서를 가진 진짜 물질처럼 취급하는 방식이라, 맥스웰 방정식을 원형 그대로 유지한 채 물성만 바꿔주면 되어 유한요소 전자기 코드에 얹기 좋다.

가장 널리 쓰이는 현대판은 CFS-PML(Complex Frequency Shifted PML) 혹은 컨볼루션 PML(CPML)이다. 좌표 확장 인자에 항을 하나 더 붙인다.

s_x = \kappa_x + \frac{\sigma_x}{\alpha_x + j\omega \varepsilon_0}

$\alpha_x$ 라는 주파수 이동 항 덕분에 저주파 성분과 격자에 거의 나란히 스치는 에바네센트파(evanescent wave)까지 잘 흡수한다. 원조 PML이 취약했던 이 두 경우를 CFS-PML이 메꿔주기 때문에, 오늘날 상용·오픈소스 FDTD 코드의 기본 경계는 사실상 CPML이다.⁴

5. 이산 격자에서의 현실[편집]

연속 방정식에서 PML의 반사 계수는 이론상 0이지만, 정작 컴퓨터가 푸는 것은 유한차분법이나 유한요소법으로 이산화한 방정식이다. 이산화하는 순간 완벽했던 정합이 깨져 수치적 반사가 생긴다. 이 반사를 억누르는 것이 PML 튜닝의 전부라고 해도 과언이 아니다.

격자 급변 금지: PML 안에서 셀 크기가 갑자기 바뀌면 그 자체가 반사원이 된다. PML 영역은 균일 격자로 두는 것이 안전하다.
최적 $\sigma_{\max}$ : 너무 작으면 흡수가 부족하고, 너무 크면 첫 층에서 수치 반사가 커진다. $\sigma_{\max}$ 에는 셀 크기와 다항식 차수로 결정되는 경험적 최적값 공식이 알려져 있다.
모서리 처리: 3차원에서 도메인 모서리와 꼭짓점은 두 방향, 세 방향의 PML이 겹치는데, 각 방향 감쇠를 독립적으로 걸어야 한다. 여기서 실수하면 코너에서만 파가 되튕겨 나온다.

PML은 안테나 해석, 광도파로, 레이더 산란 단면적(RCS) 계산, 심지어 음향·탄성파 해석까지 무한 영역을 다루는 거의 모든 파동 시뮬레이션에 이식되었다. “무한을 유한 상자에 담는 마법”이라는 별명이 붙을 만하다.⁵

6. 관련 문서[편집]

비유하자면 무반향실(anechoic chamber)의 벽에 붙은 뾰족한 흡음 폼과 정확히 같은 역할이다. 다만 실제 폼은 임피던스가 완벽히 정합되지 않아 저주파에서 반사가 남는 반면, PML은 연속 방정식 수준에서 반사가 0이라는 점이 다르다. ↩
Bérenger, J.-P. (1994). “A perfectly matched layer for the absorption of electromagnetic waves.” Journal of Computational Physics, 114(2), 185-200. 인용 횟수가 수만 회에 달하는, 전산전자기학 역사상 최다 피인용 논문 중 하나다. ↩
이 좌표 확장 해석은 1994~96년 첸(W. C. Chew)과 웨드랜드(W. H. Weedon) 등이 정립했다. 베렝제의 다소 인위적으로 보였던 분할장 트릭이 사실은 “복소 공간으로 좌표를 늘리는 것”이라는 깔끔한 기하학적 의미를 가진다는 게 밝혀진 순간, 다들 무릎을 쳤다. ↩
Roden, J. A., & Gedney, S. D. (2000). “Convolution PML (CPML): An efficient FDTD implementation of the CFS-PML.” Microwave and Optical Technology Letters, 27(5), 334-339. 컨볼루션 형태로 시간영역에 직접 구현해 메모리·연산 효율까지 챙긴 판본이라 실무 표준이 되었다. ↩
다만 PML도 만능은 아니다. 특정 이방성·분산성 매질이나 특정 주기 구조에서는 PML이 도리어 불안정(지수적으로 발산)해지는 사례가 보고되어 있다. “무한을 담는 마법”에도 결계가 필요한 조건이 있다는 뜻. ↩