고해상도 모멘트법을 이용한 정육면체와 중공 원통의 정전용량 계산

초록

본 논문은 전하 분포를 고해상도로 이산화한 모멘트법(MOM)을 적용해 정육면체와 중공 원통의 정전용량을 계산한다. 정육면체는 면당 최대 600 × 600개의 서브 영역으로 나누고, 대칭성 및 병렬 연산을 활용해 계산 효율을 높였다. 면당 90 × 90 분할에서 정전용량이 73.519 pF로 최대값을 보이며, 그 이후에는 이산화 오차로 감소한다. 중공 원통에 대해서는 원통 표면을 원형 링과 정사각형 서브 영역으로 세분화하고, 레크너(Lekner)의 행렬 전개와 실험값을 비교해 좋은 일치를 확인하였다.

상세 분석

본 연구는 전통적인 모멘트법을 고해상도 이산화와 대칭성 활용, 그리고 병렬 컴퓨팅으로 확장한 점이 가장 큰 특징이다. 정육면체의 경우, 기존 문헌에서는 6 × 6, 20 × 20, 48 × 48 등 비교적 낮은 해상도로만 다루었으나, 저자들은 면당 600 × 600까지 확장하면서도 대칭성을 이용해 실제 계산량을 10 × 10 행렬(P)만으로 축소하였다. 이는 각 서브 영역이 동일한 전위(1 V) 조건에서 전하 밀도 σ를 구하고, 전체 전하 Q를 합산해 정전용량 C = Q/φ를 얻는 전통적인 MOM 절차를 유지하면서도 메모리와 연산 시간을 크게 절감한다는 의미다.

특히, 전위 근사식 A2(σ/(4πϵ₀d))가 실제 전위보다 과대평가되는 점을 지적하고, 이로 인해 고해상도일수록 전하 밀도가 과소평가되어 정전용량이 감소하는 현상을 설명한다. 이는 전위 근사의 1 V 고정값이 실제 전위 분포와 차이를 만들기 때문에 발생하는 시스템적 오차이며, “Thomson 정리”를 통해 전기장 에너지 최소화와 정전용량 증가 사이의 관계를 이론적으로 뒷받침한다.

정육면체 결과는 90 × 90 분할에서 73.519 pF라는 피크값을 보이며, 이는 기존 0.660 a(cm) 수준(≈73.5 pF)과 일치한다. 그러나 600 × 600까지 확대하면 오히려 73.07 pF 수준으로 감소한다. 이는 이산화 오차가 지배적으로 작용함을 의미한다. 저자들은 이 현상을 “이산화 오류가 무한히 작은 전하 분포를 방해한다”는 점에서 MOM의 근본적 한계로 제시한다.

중공 원통 부분에서는 원통을 L개의 원형 링으로 나누고, 각 링을 K개의 정사각형 서브 영역으로 세분화한다. 전위와 전하 밀도를 동일하게 1 V 고정하고, 링 간 거리 d와 자체 전위 항을 이용해 전위 방정식을 구성한다. 여기서 K는 링의 둘레(2πR)와 동일하게 설정되어, 원통의 축대칭성을 완전히 활용한다. 계산된 C/R 비율을 레크너의 행렬 전개(6차 truncation)와 비교했을 때, 대부분의 L/K 조합에서 오차가 0.01 % 이하로 매우 정확함을 확인한다. 다만 L/K 비율이 극단적으로 작을 때(예: 5000/1571, 5000/524)에는 계산 자원의 제한으로 인한 오차가 약간 커진다.

실험적 검증으로는 18세기 Cavendish의 전기계 측정값을 인용했으며, 이론값과 실험값 사이에 존재하는 차이는 측정 장비의 한계와 실험 환경(전하 누설, 절연 불완전성 등) 때문이라고 설명한다. 전반적으로, 고해상도 MOM이 정전용량 계산에 있어 높은 정확도를 제공하지만, 이산화 오차와 근사식 선택이 결과에 큰 영향을 미친다는 점을 강조한다.

이 논문은 전산 전자기학에서 대규모 이산화와 대칭성 활용, 그리고 병렬 처리 기법을 결합한 실용적인 접근법을 제시함으로써, 복잡한 3차원 도체의 정전용량을 고정밀도로 예측하고, 기존 문헌과 실험값을 통합 검증하는 데 기여한다.