적응형 웨이브렛 격자법을 이용한 시간 의존 맥스웰 방정식 시뮬레이션

초록

본 논문은 전자기 파동을 기술하는 맥스웰 방정식을 풀기 위해, 매 시간 단계마다 필드의 웨이브렛 분해를 이용해 격자를 동적으로 재구성하는 적응형 웨이브렛 콜로케이션(time‑domain) 방법을 제안한다. 고밀도 전자기 현상이 발생하는 영역에는 격자점을 집중시키고, 필드가 거의 없는 영역에서는 격자밀도를 크게 낮춘다. 이렇게 얻어진 적응 격자 위에서 고차 정확도의 공간 차분식과 명시적 시간 적분법을 적용함으로써, 압축률이 높고 연산 비용이 크게 감소하는 시뮬레이션이 가능함을 보인다. 특히 가이드파형 광학 소자와 같은 복잡한 구조의 전파 해석에 유리하다.

상세 분석

이 연구는 웨이브렛 기반 적응 격자 생성 메커니즘을 맥스웰 방정식의 시간 전진 스키마와 결합함으로써, 전통적인 고정 격자 방식이 갖는 공간 해상도와 계산량 사이의 트레이드오프를 근본적으로 해소한다는 점에서 의미가 크다. 웨이브렛 변환은 신호의 지역적 특성을 다중 해상도 형태로 표현할 수 있어, 계수의 절댓값이 일정 임계값 이하인 경우 해당 영역을 자동으로 제거한다. 논문에서는 Daubechies‑4와 같은 정규 직교 웨이브렛을 선택해, 전기장·자기장 각각에 대해 2‑차원 다중 스케일 분해를 수행한다. 이때, 세부 레벨(l)과 위치(k) 인덱스를 이용해 동적으로 활성화된 격자 포인트 집합 Gⁿ을 정의하고, Gⁿ⁺¹을 생성하기 위해 현재 시각 n에서의 필드 계수와 전파 속도에 기반한 예측 보정 절차를 도입한다.

공간 차분 연산자는 적응 격자에 맞게 비균일 스텝을 허용하도록 설계되었으며, 중앙 차분 형태의 curl 연산자를 유지하면서도 각 격자점에서 필요한 이웃 포인트를 웨이브렛 재구성을 통해 자동으로 확보한다. 시간 적분은 CFL 조건을 만족하도록 제한된 Δt를 사용한 2차 정확도 명시적 스킴을 적용했으며, 이는 기존 FDTD와 동일한 안정성 한계를 갖지만, 격자 포인트 수가 크게 감소함에 따라 실제 연산 시간은 획기적으로 단축된다.

압축 효율을 정량화하기 위해, 동일한 물리적 도메인에 대해 고정 격자 FDTD와 비교했을 때, 평균 격자 점 수가 10% 이하로 감소하면서도 L₂ 오차는 10⁻⁴ 수준으로 유지되는 것을 실험적으로 확인했다. 또한, 파동이 급격히 변하는 경계층이나 비선형 매질 내부에서도 웨이브렛 기반 적응이 충분히 세밀한 해상도를 제공함을 보여준다. 이러한 결과는 적응형 웨이브렛 콜로케이션이 고주파 광학 디바이스, 광섬유 모드 해석, 그리고 플라즈마 시뮬레이션 등 복잡한 전자기 현상을 효율적으로 다룰 수 있는 강력한 도구임을 시사한다.