[논문 리뷰] A phase-field approach for detecting cavities via a Kohn-Vogelius type functional

Andrea Aspri|arXiv (Cornell University)|2022. 02. 24.

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

한 줄 요약

이 논문은 선형 탄성 매질 내의 구멍을 복원하기 위한 새로운 단계장 방법을 제안한다. 이 방법은 둘레 정규화로 정규화된 Kohn-Vogelius 유형의 기능을 사용한다. 둘레의 Modica-Mortola 완화와 소규모 탄성 텐서를 사용하는 가짜 재료 접근법을 결합함으로써, 형태 최적화 문제를 H¹(Ω)의 볼록 부분집합 위에서 매끄럽고 미분 가능한 최적화 문제로 변환한다. 이는 반복적 원-이중 알고리즘을 통해 안정적이고 효율적인 수치적 복원을 가능하게 하며, 노이즈가 있는 데이터 조건에서도 정확한 구멍 형상을 수렴한다.

ABSTRACT

We deal with the geometrical inverse problem of the shape reconstruction of cavities in a bounded linear isotropic medium by means of boundary data. The problem is addressed from the point of view of optimal control: the goal is to minimize in the class of Lipschitz domains a Kohn-Vogelius type functional with a perimeter regularization term which penalizes the perimeter of the cavity to be reconstructed. To solve numerically the optimization problem, we use a phase-field approach, approximating the perimeter functional with a Modica-Mortola relaxation and modeling the cavity as an inclusion with a very small elastic tensor. We provide a detailed analysis showing the robustness of the algorithm through some numerical experiments.

연구 동기 및 목표

• 선형 탄성 매질 내의 구멍을 경계에서의 응력 및 변위 측정치로부터 식별하는 불안정한 역문제를 해결하는 것.
• 약한 안정성 추정치가 존재하는 상황에서도 안정성과 수렴성을 보장하는 구멍 복원을 위한 강력한 수치 알고리즘을 개발하는 것.
• 기존에는 이 분야에 적용된 바가 없는 Kohn-Vogelius 유형 기능에 단계장 방법을 확장하는 것.
• 단계장 완화와 둔변 정규화를 통해 형태 복원 문제를 수치적으로 안정적이고 미분 가능한 형태로 제시하는 것.

제안 방법

• 두 개의 경계값 문제의 해 간의 L² 에너지 격차로 정의된 Kohn-Vogelius 기능의 제약 최적화 문제로 역구멍 탐지 문제를 수립한다.
• 불안정 문제를 안정화하고 복잡한 구멍 경계를 방지하기 위해 둔변 정규화 항을 도입한다.
• 구멍을 나타내기 위해 스칼라 단계 변수 v ∈ [0,1]을 도입함으로써 단계장 접근법을 적용한다. 여기서 C1 = δC0 (δ > 0이 매우 작음)이며, 구멍을 낮은 강성의 포함체로 모델링한다.
• 에너지원 ε|∇v|² + (1/ε)v(1−v)을 사용하여 둘레를 Modica-Mortola 기능으로 근사함으로써, 구멍 경계 부근에 너비 ε의 분산 경계를 생성한다.
• Kohn-Vogelius 에너지와 단계장 둔변 근사의 조합으로, Frechét 미분 가능성이 보장되는 이완 기능 Jδ,ε(v)를 유도한다. 이는 H¹(Ω) 위에서 정의된다.
• 이완 기능의 1차 최적성 조건에 기반한 반복적 원-이중 활성집합 방법을 구현하여 수치적 복원을 수행한다.

실험 결과

연구 질문

• RQ1 Kohn-Vogelius 유형 기능에 대해 단계장 방법을 선형 탄성에서의 구멍 복원에 성공적으로 적용할 수 있는가?
• RQ2 Modica-Mortola 둔변 완화와 가짜 재료 모델링의 조합이 형태 복원에 대해 안정적이고 수렴하는 수치적 스킴을 제공하는가?
• RQ3 기존의 표준 오차 기능에 비해, 노이즈가 있는 경계 측정치 조건에서 제안된 방법은 어떻게 성능을 발휘하는가?
• RQ4 비볼록 형상에 대해서도 단계장 완화가 진정한 구멍의 기하적 특성을 어느 정도 유지하는가?

주요 결과

• Kohn-Vogelius 기능의 단계장 완화와 둔변 정규화는 H¹(Ω) 위에서 Frechét 미분 가능성을 보장하는 기능을 생성하여 효율적인 수치 최적화를 가능하게 한다.
• 수치 실험 결과, δ와 ε가 충분히 작을 경우, 이완 기능 Jδ,ε의 최소값은 원래 Jreg의 최소값을 높은 정확도로 근사한다.
• 원형, 타원형, 직사각형 등의 볼록 구멍의 경우, 경계 측정치에 최대 6.5%의 노이즈가 존재하더라도 정확한 복원을 달성한다.
• 두 개의 원형 포함체와 같은 다중 구멍의 경우, 5% 노이즈 조건에서 n = 35,318회 반복 후 수렴이 관찰된다.
• 비볼록 영역의 복원 결과는 Modica-Mortola 둔변 근사로 인해 볼록화되는 경향을 보이며, 이는 비볼록 기하학에 대한 한계를 시사한다.
• 오차 기능(Jmisfitδ,ε)과의 비교 결과, 다중 포함체의 경우 Kohn-Vogelius 기반 방법이 다소 더 우수한 성능을 보였지만, 경계 아티팩트가 더 많이 발생하였다.