삼각형 요소를 이용한 거의 정확한 BEM 솔루션

초록

본 논문은 평면 삼각형 및 직사각형 요소에 균일하게 분포된 특이원(전하, 전위 등)의 영향을 정확히 계산하는 C 함수 라이브러리(ISLES)를 개발하고, 이를 기반으로 한 neBEM 솔버의 성능을 검증한다. 삼각형 요소에 대한 폐곡선 적분식을 유도하고, 복잡한 식의 계산 시간과 전통적인 수치 적분 방식과의 효율성을 비교한다. 또한 단위 전압을 가한 정사각형 플레이트의 정전용량을 고정밀 결과와 비교하여 정확성을 입증하고, 코너 근처 전하밀도 분포가 이산화 수준에 따라 매끄럽게 수렴함을 보여준다.

상세 분석

이 연구는 경계요소법(BEM)의 근본적인 한계인 ‘특이점 분포의 근사’ 문제를 근본적으로 해결하고자 한다. 기존 BEM에서는 요소 내부에 점특이점(예: 점전하)이나 평균값을 가정해 근사적으로 적분을 수행한다. 그러나 이러한 접근은 요소 경계, 특히 코너와 같은 급격한 기하학적 변화 부근에서 큰 오차를 초래한다. 저자들은 균일 표면 전하(또는 전위) 분포를 정확히 적분한 폐곡선 해를 삼각형 요소에 대해 유도하고, 이를 C 함수 형태로 구현한 ISLES 라이브러리를 제공한다. 핵심은 삼각형의 세 변과 내부 각을 파라미터화하여, 로그와 아크탄젠트 함수들의 조합으로 전위와 플럭스(전기장)의 정확한 폐곡선 표현식을 얻는 것이다.

식 자체는 길고 복잡하지만, 저자들은 이를 고정밀 부동소수점 연산과 메모리 효율적인 루프 구조로 최적화함으로써, 동일한 정확도를 요구하는 수치 적분(예: 가우스‑레전드)보다 현저히 짧은 실행 시간을 달성한다. 특히, 요소당 연산량이 O(1)인 반면, 수치 적분은 적분점 수에 비례해 O(N)으로 증가한다는 점에서 큰 장점을 가진다.

또한, neBEM 솔버에 ISLES를 통합함으로써 전통적인 BEM에서 발생하던 ‘전하밀도 발산’ 문제를 완화한다. 정사각형 플레이트의 전하밀도는 코너 근처에서 급격히 증가하는 경향이 있는데, 균일 분포에 대한 정확한 영향 함수를 사용하면 이 현상이 물리적으로 타당한 형태로 매끄럽게 수렴한다. 이는 요소 크기를 점점 작게 하여도 전하밀도 프로파일이 수렴함을 실험적으로 확인함으로써 입증된다.

마지막으로, 저자들은 단위 전압을 가한 정사각형 플레이트의 정전용량을 계산하고, 기존 문헌에 보고된 고정밀 결과(예: 전위 적분법, 고차원 적분법)와 비교한다. neBEM이 제공하는 값은 0.3668 pF 정도로, 5~6자리까지 일치한다. 이는 복잡한 3차원 전자기 문제에서도 ISLES 기반 BEM이 높은 정확도와 효율성을 동시에 제공한다는 강력한 증거이다.

요약하면, 이 논문은 삼각형 요소에 대한 정확한 폐곡선 적분식을 수학적으로 도출하고, 이를 실용적인 C 라이브러리로 구현함으로써 BEM의 근본적인 근사오차를 제거한다. 계산 효율성, 정전용량 정확도, 코너 근처 전하밀도 수렴성 등 다각적인 검증을 통해 neBEM의 우수성을 입증하였다. 이러한 접근은 정전기뿐 아니라 전자기, 열전도, 유체역학 등 다양한 분야의 경계값 문제에 적용 가능하며, 고정밀 시뮬레이션이 요구되는 최신 공학·과학 연구에 큰 파급 효과를 기대한다.