오픈소스 FDTD 솔버 비교 Elecode gprMax MEEP의 번개 전자기장 시뮬레이션 적용성

초록

본 논문은 오픈소스 FDTD 솔버인 gprMax, Elecode, MEEP을 이용해 번개 전자기장 전파를 재현하고, 완전 전도성 및 손실성 지면 위에서의 기준 해와 비교하였다. 모든 솔버가 대체로 만족스러운 정확도를 보였지만, 공간 격자 크기와 경계 조건 설정에 따라 인공적인 수치 분산 및 반사 현상이 발생한다는 점을 강조한다. 각 솔버의 특징, 성능, 제한점 및 사용상의 장단점을 정리하고, 재현 가능한 스크립트를 공개한다.

상세 요약

이 연구는 번개 전자기장 전파라는 고주파, 급격한 전압 상승을 포함하는 복합 현상을 FDTD 방식으로 모델링할 때, 오픈소스 솔버가 상업용 도구와 견줄 만한 신뢰성을 제공하는지를 체계적으로 검증한다. 첫 번째 핵심은 공간 격자 해상도이다. 번개 전류 파형은 수십 MHz에서 수백 MHz까지 넓은 스펙트럼을 포함하므로, 최소 10 점/파장 이상의 샘플링이 필요하다. 논문에서는 Δx = 0.5 m, Δt = Δx/(2c) 정도의 설정이 수치 분산을 최소화하면서도 계산 비용을 합리적으로 유지한다는 결론을 제시한다. 두 번째는 경계 조건이다. PML(Perfectly Matched Layer) 두께와 손실 프로파일을 부적절하게 설정하면, 실제 물리적 반사와 구분이 어려운 인공 반사가 발생한다. 특히 지면을 완전 전도성으로 모델링할 때는 지면 경계 자체가 반사면이 되므로, PML을 지면 위에 배치하거나 지면을 전도성 경계 조건으로 직접 지정하는 방법을 비교한다. 세 솔버 모두 PML 구현 방식에 차이가 있다. gprMax는 1‑D 스택형 PML을 기본 제공하지만, 복합 지형에서는 사용자 정의가 제한적이다. Elecode는 전자기‑전기 회로 결합 모델링에 강점이 있어, 지면 임피던스를 직접 입력할 수 있다. MEEP은 고차원 PML과 복합 물질 모델링이 가장 유연하지만, 스크립트 작성이 다소 복잡하다. 세 번째는 물질 모델링 정확도이다. 손실성 토양을 복소 유전율 ε = ε′ − jσ/ω 로 표현할 때, 각 솔버가 주파수 의존성을 어떻게 처리하는지가 결과에 큰 영향을 미친다. gprMax는 고정된 복소 유전율을 사용하므로 저주파 영역에서 과소평가가 발생하고, Elecode는 주파수‑시간 변환을 통해 동적 손실을 구현한다. MEEP은 디지털 필터를 이용해 손실을 주파수‑의존적으로 적용할 수 있다. 마지막으로 성능 측면에서는 병렬화와 GPU 가속 지원이 차별화 포인트다. gprMax와 MEEP은 MPI와 CUDA를 통한 확장이 용이하지만, Elecode는 현재 CPU 전용이며 대규모 3‑D 시뮬레이션에서 메모리 한계에 직면한다. 종합적으로, 논문은 각 솔버가 제공하는 기능과 제한을 명확히 구분하고, 사용자가 시뮬레이션 목표(정확도 vs. 계산 시간)와 물리적 조건(지면 특성, 파형 스펙트럼)에 따라 적절한 도구를 선택하도록 가이드한다.