부분 경계 측정을 이용한 유전율 및 전도도 재구성 변분 최적화 접근법

초록

본 논문은 제한된 경계 전기장 관측을 활용해 시간 의존 맥스웰 방정식의 유전율과 전도도 두 계수를 동시에 복원하는 변분 최적화 방법을 제시한다. 라그랑지안의 약한 형태와 최적성 조건을 유도하고, 인접 문제의 안정성 추정 및 정칙화된 티크노프 함수의 프레셰 미분 가능성을 증명한다. 2·3 차원 수치 실험에서는 공액 기울기(CGA)와 적응형 공액 기울기(ACGA) 알고리즘을 적용해 악성 흑색종 검출 등 실제 의료 영상에 대한 유효성을 확인한다.

상세 분석

이 연구는 전도성 매질에서 시간 의존 맥스웰 방정식을 비자기성 가정 하에 정규화된 형태(식 2.8)로 변형하고, 유전율 ε(x)와 전도도 σ(x)를 공간 함수로 두어 부분 경계 데이터만을 이용한 역문제(IP)를 설정한다. 기존의 스칼라 파동 방정식 기반 방법과 달리 전기장 E 자체를 직접 다루어 복합 물성(복소 유전율) 복원을 가능하게 한다. 핵심 아이디어는 정칙화 티크노프 함수 F(ε,σ)와 제약식(전방 방정식)을 라그랑지안 L(E,λ,ε,σ)에 결합해 최적화 문제를 전역적으로 풀 수 있게 하는 것이다.

논문은 먼저 L의 약한 형태(Lemma 3.1)를 도출하고, 이를 기반으로 변분 미분(Fréchet derivative)을 구해 최적성 조건(식 3.8‑3.11)을 제시한다. 특히 ε와 σ에 대한 파라미터 미분은 정칙화 항 γ_ε‖ε−ε₀‖², γ_σ‖σ−σ₀‖²와 결합돼 안정적인 역문제 해를 보장한다. 인접 문제는 라그랑지안에서 λ에 대한 변분을 통해 얻어지며, Lemma 4.2에서 제시된 에너지 추정과 Theorem 4.3의 안정성 추정은 데이터 노이즈 δ에 대한 연속성을 확보한다.

프레셰 미분 가능성은 함수 공간 U=H¹_E×H¹_λ×C_ε×C_σ에서 L이 2차 연속 미분 가능함을 증명함으로써, 수치 최적화 알고리즘(공액 기울기, 적응형 공액 기울기)의 수렴 이론적 기반을 제공한다. 또한 정칙화 파라미터 선택에 대한 가이드라인과, 비정칙화와 정칙화 티크노프 함수 모두 최소점 존재성을 보이는 Theorem 6.1을 제시해 전역 최소화 가능성을 확보한다.

수치 구현에서는 도메인 분할(Ω_FDM, Ω_FEM)과 유한 요소(FEM) 기반 공간 이산화를 적용하고, 시간 전진-후진 스키마(전방-인접) 루프를 공액 기울기 방식으로 최적화한다. 2‑D 실험에서는 단순 원형·다각형 물체의 ε와 σ를 정확히 복원했으며, 3‑D 실험에서는 실제 악성 흑색종 모델(6 GHz, 복합 유전율) 데이터를 사용해 물체 형태와 물성값을 높은 정밀도로 재구성했다. 결과는 관측 영역이 제한적이고 노이즈가 존재함에도 불구하고 재구성 오차가 5 % 이하로 수렴함을 보여, 의료 마이크로파 영상 및 폭발물 탐지와 같은 실용 분야에 적용 가능함을 시사한다.