[논문 리뷰] Lipschitz stability in an inverse problem for the wave equation

Lucie Baudouin|arXiv (Cornell University)|2011. 06. 08.

Numerical methods in inverse problems참고 문헌 12인용 수 9

한 줄 요약

이 논문은 파동 방정식에 대한 역문제에서 리프시츠 안정성을 확립한다. $ \partial_{tt}u - \Delta u + p(x)u = 0 $에서 잠재력 $ p(x) $는 $ \partial\Omega $의 부분집합인 $ \Gamma_0 \subset \partial\Omega $에서 정규 도함수 $ \partial_\nu u $의 측정값으로 유일하고 안정적으로 재구성될 수 있음을 증명한다. 핵심 결과는 전역 카르레만 추정을 통해 도출되며, 이는 $ \Gamma_0 $와 관측 시간 $ T $에 대한 기하 조건 하에서 안정성이 확보됨을 보여주며, 데이터 정규성과 초기 데이터의 하한에 대한 명시적 의존성을 포함한다.

ABSTRACT

We are interested in the inverse problem of the determination of the potential $p(x), x\inΩ\subset\mathbb{R}^n$ from the measurement of the normal derivative $\partial_νu$ on a suitable part $Γ_0$ of the boundary of $Ω$, where $u$ is the solution of the wave equation $\partial_{tt}u(x,t)-Δu(x,t)+p(x)u(x,t)=0$ set in $Ω imes(0,T)$ and given Dirichlet boundary data. More precisely, we will prove local uniqueness and stability for this inverse problem and the main tool will be a global Carleman estimate, result also interesting by itself.

연구 동기 및 목표

• 파동 방정식에서 잠재력 $ p(x) $를 경계에서 정규 도함수 $ \partial_\nu u $의 측정값으로부터 복원하는 역문제에 대해 局부 리프시츠 안정성을 확립하는 것.
• 관측 영역 $ \Gamma_0 $와 관측 시간 $ T $에 대한 기하 조건 하에서 유일성 및 안정성 결과를 증명하는 것. 이는 재구성에 필요한 충분한 정보를 확보하기 위함이다.
• 이전 연구에서 사용된 컴팩트성-유일성 추론을 피하고, 전역 카르레만 추정을 직접 활용하여 안정성의 증명을 제공하는 것.
• 해의 정규성에 대한 가정을 약화시키고, 전체 경계가 아닌 부분집합 $ \Gamma_0 $에서의 측정값으로도 충분함을 보이는 것.

제안 방법

• 해의 정규성에 따라 $ |x - x_0|^2 $와 시간 $ t $에 의존하는 가중 함수 $ \varphi $를 사용하여 파동 연산자 $ \partial_{tt} - \Delta + q(x) $에 대한 전역 카르레만 추정을 유도하는 것. 여기서 $ x_0 \notin \Omega $이다.
• 역문제의 선형화된 형태에 카르레만 추정을 적용한다. 이 경우 해 $ u $는 동차 초깃값 및 경계 조건을 만족하며, 소스 항 $ f(x)R(x,t) $에 의해 구동된다.
• 공간과 시간에 걸쳐 해를 국소화하기 위해 컷오프 함수 $ \chi $를 사용하여 에너지 및 플럭스 항을 카르레만 부등식 내에서 제어할 수 있도록 한다.
• 카르레만 추정을 활용하여, 경계 플럭스 $ \partial_\nu u $가 $ \Gamma_0 \times (0,T) $에서 관측된 결과로 소스 항 $ f $의 가중 $ L^2 $-노름을 추정하는 관측 유형의 부등식을 수립한다.
• 에너지 추정과 가중 함수의 성질을 이용하여 저차항을 흠집내고, 최종적으로 $ \|f\|_{L^2} $와 경계 측정값 $ \|\partial_\nu u\|_{L^2(\Gamma_0 \times (0,T))} $ 사이의 관계를 유도한다.
• 마지막으로 리프시츠 안정성 추정을 도출한다: $ \|f\|_{L^2(\Omega)} \leq C \|\partial_\nu u\|_{L^2(\Gamma_0 \times (0,T))} $, 여기서 $ C $는 $ m, T, K, r $에 의존한다.

실험 결과

연구 질문

• RQ1 관측 시간 $ T $가 충분할 경우, 파동 방정식에서 잠재력 $ p(x) $는 $ \partial\Omega $의 부분집합인 $ \Gamma_0 \subset \partial\Omega $에서 정규 도함수 $ \partial_\nu u $의 측정값으로 유일하게 결정될 수 있는가?
• RQ2 역문제에 대해 리프시츠 안정성 추정을 확립할 수 있는가? 이 경우 $ L^2 $-노름으로서의 차 $ q - p $의 크기를 경계 플럭스의 차 $ \partial_\nu y[q] - \partial_\nu y[p] $의 $ L^2 $-노름으로 유계화할 수 있는가?
• RQ3 전역 카르레만 추정을 통해 더 약한 정규성 조건과 전체 경계가 아닌 부분집합 $ \Gamma_0 $에서의 측정만으로도 안정성 결과를 확보할 수 있는가?
• RQ4 $ \Gamma_0 $와 $ T $에 대해 어떤 기하적 및 시간적 조건에서 파동 방정식의 역문제가 국부적으로 리프시츠 안정적인가?
• RQ5 안정성 상수는 데이터 정규성과 초기 데이터 하한에 대해 명시적으로 추정할 수 있는가?

주요 결과

• 논문은 국부적 리프시츠 안정성을 증명한다: $ q, p \in L^\infty_{\leq m}(\Omega) $일 때, 차 $ \|q - p\|_{L^2(\Omega)} $는 $ \|\partial_\nu y[q] - \partial_\nu y[p]\|_{H^1(0,T;L^2(\Gamma_0))} $의 상수배로 유계화되며, 상수 $ C = C(m, T, K, r) $는 $ L^\infty $-상한 $ m $, 시간 $ T $, 해의 정규성 $ K $, 初기 데이터 하한 $ r $ 에 의존한다.
• 안정성 결과는 다음 기하 조건 하에서 성립한다: $ x_0 \notin \Omega $이 존재하여 $ \Gamma_0 \supset \{x \in \partial\Omega : (x - x_0) \cdot \nu(x) \geq 0\} $ 이고, $ T > \sup_{x \in \Omega} |x - x_0| $ 이다.
• 증명은 파동 연산자에 대한 전역 카르레만 추정에 기반하며, 이는 전체 경계가 아닌 부분집합 $ \Gamma_0 $에서의 측정값을 사용할 수 있음을 허용한다.
• 이전 연구에서 사용된 컴팩트성-유일성 추론을 피하고, 카르레만 추정에서 직접 안정성 증명을 제공한다.
• 카르레만 추정은 초깃값 및 경계 조건이 0인 선형화된 문제에 적용되며, 여기서 소스 항 $ f(x)R(x,t) $는 경계 플럭스 $ \partial_\nu u $로부터 재구성된다.
• 최종 안정성 추정은 $ \|f\|_{L^2(\Omega)} \leq C \|\partial_\nu u\|_{L^2(\Gamma_0 \times (0,T))} $ 형태로 유도되며, $ C $는 $ m, T, K, r $에 의존한다. 이 결과는 선형화를 통해 비선형 역문제로 확장된다.