[논문 리뷰] Discrete Geometric Structures in Homogenization and Inverse Homogenization with application to EIT

Mathieu Desbrun, Roger Donaldson외 1명|ArXiv.org|2009. 04. 17.

Numerical methods in inverse problems참고 문헌 55인용 수 3

한 줄 요약

이 논문은 미세구조의 전도도 균질화를 볼륨 평균 또는 볼륨 평균을 통한 선형 보간으로 재구성하는 기하학적 프레임워크를 제안한다. 이는 가중치가 부여된 델로아 삼각형 메esh를 통해 안정적이고 최적의 균질화를 가능하게 한다. 2차원 전기 임피던스 단층촬영(EIT)에 적용하여 경계 데이터로부터 안정적이고 비등방성의 전도도 재구성을 가능하게 하며, 볼록 함수 공간에서 해를 매개변수화함으로써 이전까지 달성되지 못했던 미세구조의 해상도 향상을 이룬다.

ABSTRACT

We introduce a new geometric approach for the homogenization and inverse homogenization of the divergence form elliptic operator with rough conductivity coefficients $σ(x)$ in dimension two. We show that conductivity coefficients are in one-to-one correspondence with divergence-free matrices and convex functions $s(x)$ over the domain $Ω$. Although homogenization is a non-linear and non-injective operator when applied directly to conductivity coefficients, homogenization becomes a linear interpolation operator over triangulations of $Ω$ when re-expressed using convex functions, and is a volume averaging operator when re-expressed with divergence-free matrices. Using optimal weighted Delaunay triangulations for linearly interpolating convex functions, we obtain an optimally robust homogenization algorithm for arbitrary rough coefficients. Next, we consider inverse homogenization and show how to decompose it into a linear ill-posed problem and a well-posed non-linear problem. We apply this new geometric approach to Electrical Impedance Tomography (EIT). It is known that the EIT problem admits at most one isotropic solution. If an isotropic solution exists, we show how to compute it from any conductivity having the same boundary Dirichlet-to-Neumann map. It is known that the EIT problem admits a unique (stable with respect to $G$-convergence) solution in the space of divergence-free matrices. As such we suggest that the space of convex functions is the natural space in which to parameterize solutions of the EIT problem.

연구 동기 및 목표

• 미세구조 매개변수에 의존하지 않는, 불규칙한 전도도 계수를 가진 타원형 연산자의 기하학적 균질화 프레임워크를 개발하는 것.
• EIT에서의 불안정한 역균질화 문제를 선형 불안정 문제와 잘 정의된 비선형 문제로 분해하여 해결하는 것.
• 전도도를 볼록 함수 공간에서 매개변수화함으로써 안정적이고 강건한 EIT 재구성 방법을 제공하고, 비등방성 미세구조의 해상도를 향상시키는 것.
• 이산 볼록 함수 보간을 위한 가중치가 부여되고 Q-적응형 메쉬로 확장된 델로아 삼각형 메쉬의 최적성 특성을 확장하는 것.

제안 방법

• 전도도 σ(x)를 발산이 없는 행렬 Q와 볼록 함수 s(x) 사이의 일대일 대응을 통해 매개변수화함으로써, 균질화를 선형 보간 또는 체적 평균으로 재표현하는 것.
• 이산 딜라클 에너지를 최소화하는 데 있어 전역 최적성을 보장하는 최적의 가중치가 부여된 델로아 삼각형 메쉬를 사용하여 볼록 함수 s(x)를 보간하는 것.
• 삼각형 스텐실에서 이차 보간을 사용하여 이산 딜라클-노이만(DtN) 매ap Λ_sh를 s^h_i에서 구성하며, 경계 안정성을 확보하기 위해 가상 정점(ghost vertices)을 포함하는 것.
• 데이터 오차와 tr(Q^h)의 총변동성 정규화를 최소화하는 최적화를 통해 역 EIT 문제를 해결하며, q^h_ij ≥ 0 조건을 만족시키는 것.
• IpOpt를 사용한 효율적 최적화를 위해 데이터 오차 항의 야코비안을 원시-수반 방법으로 계산하는 것.
• 계산의 실현 가능성을 확보하기 위해 tr(Q^h)를 정규화 항으로 사용하며, 실험적으로 티호노프 파라미터 α를 조정하는 것.

실험 결과

연구 질문

• RQ1 불규칙한 전도도 계수의 균질화가 비선형성을 피하는 기하학적 공간에서 선형 연산자로 재구성될 수 있는가?
• RQ2 이산 볼록 함수 보간을 위한 최적의 메쉬를 어떻게 구성할 수 있으며, 이는 균질화의 안정성과 정확성 보장을 위한가?
• RQ3 역 EIT 문제는 잘 정의된 비선형 문제와 선형 불안정 문제로 분해될 수 있는가? 이는 안정적인 재구성을 가능하게 하는가?
• RQ4 전도도를 볼록 함수로 매개변수화하는 것이, EIT에서 미세구조의 비등방성 특성을 더 효과적으로 해소하는가?
• RQ5 이 프레임워크는 안정적인 해상도 이하의 피치를 가진 레이어형 미세구조의 상향 스케일링 표현을 재구성할 수 있는가?

주요 결과

• 전도도 균질화가 볼록 함수 s(x) 공간에서 표현될 경우 선형 보간 연산자로, 발산이 없는 행렬 Q 공간에서는 체적 평균 연산자로 표현된다.
• 볼록 함수에 기반한 가중치가 부여된 델로아 삼각형 메쉬는 이산 딜라클 에너지를 최소화하는 데 있어 전역 최적성을 달성하며, 기존 델로아 삼각형 메쉬의 성질을 확장한다.
• 이 방법은 EIT에서 비등방성 전도도 패tern을 성공적으로 재구성하며, 비등방성의 방향성과 강도를 Q의 고유값 차이를 통해 시각화할 수 있다.
• 도식 2.2에 나타난 레이어형 미세구조의 재구성 결과는 저자들에 따르면 이전까지 달성되지 못했던 수준의 세부 정보를 포착하였으며, 해상도 이하의 피치를 가진 상태에서도 상향 스케일링된 비등방성 특성을 유지한다.
• 수치 실험을 통해 tr(Q^h)를 사용한 정규화가 det(Q^h)보다 계산 효율성이 뛰어나면서도 재구성 품질을 손상시키지 않음을 확인하였다.
• 일부 표준 방법들과 달리, 이 방법은 전도도 재구성에서 다이내믹 레인지의 과소평가를 방지하며, 발산이 없는 행렬 공간에서 G-수렴성에 대해 안정성을 유지한다.