[논문 리뷰] Inverse problems for the stationary transport equation in the diffusion scaling

Ru-Yu Lai, Qin Li외 1명|arXiv (Cornell University)|2018. 08. 06.

Numerical methods in inverse problems참고 문헌 22인용 수 5

한 줄 요약

이 논문은 확산 스케일링 하에서 복사 전달 방정식(RTE)의 광학적 매개변수 복원 역문제에 대한 안정성 추정을 수립한다. 이는 문제의 해가 리프시츠 안정성을 보이지만, 쿤델 수치(Kn)가 감소할수록 리프시츠 상수가 지수적으로 증가함을 보여주며, 이는 소규모-Kn 영역에서 확산 방정식 역문제의 심각한 불안정성 원인을 설명한다. 이는 약간의 불안정성을 보이는 RTE와 심각한 불안정성을 보이는 확산 방정식 역문제 사이의 격차를 메운다.

ABSTRACT

We consider the inverse problem of reconstructing the optical parameters of the radiative transfer equation (RTE) from boundary measurements in the diffusion limit. In the diffusive regime (the Knudsen number $\mathsf{Kn}\ll 1$), the forward problem for the stationary RTE is well approximated by an elliptic equation. However, the connection between the inverse problem for the RTE and the inverse problem for the elliptic equation has not been fully developed. This problem is particularly interesting because the former one is mildly ill-posed , with a Lipschitz type stability estimate, while the latter is well known to be severely ill-posed with a logarithmic type stability estimate. In this paper, we derive stability estimates for the inverse problem for RTE and examine its dependence on $\mathsf{Kn}$. We show that the stability is Lipschitz in all regimes, but the coefficient deteriorates as $e^{\frac{1}{\mathsf{Kn}}}$, making the inverse problem of RTE severely ill-posed when $\mathsf{Kn}$ is small. In this way we connect the two inverse problems. Numerical results agree with the analysis of worsening stability as the Knudsen number gets smaller.

연구 동기 및 목표

• 소규모-Kn 영역에서 복사 전달 방정식(RTE)의 역문제와 확산 방정식 간의 관계를 명확히 하기.
• RTE 역문제의 안정성과 쿤델 수치(Kn)에 대한 의존성을 분석하기.
• RTE 역문제가 약간의 불안정성(리프시츠 안정성)으로 보이지만, 동일한 조건에서 확산 방정식 역문제가 심각한 불안정성(로그 안정성)으로 보이는 이유를 설명하기.
• Kn→0에 따라 리프시츠 상수가 지수적으로 악화됨을 드러내는 정량적 안정성 추정을 유도하기.
• 소규모-Kn에서 안정성 악화가 증가함을 보여주는 수치 실험을 통해 이론적 결과를 검증하기.

제안 방법

• 확산 스케일링 하에서 RTE의 총 흡수 계수 σ_Kn 복원 역문제에 대한 안정성 추정을 유도하기.
• 알베도 연산자 커널의 특이 분해를 사용하여 경계 측정치의 구조 분석하기.
• RTE 역문제에 대해 리프시츠 유형의 안정성 추정을 수립하며, 리프시츠 상수가 1/Kn에 대해 지수적으로 의존함을 보장하기.
• Kn→0에서 알베도 연산자의 점근적 행동을 분석하여, 주요 특이 부분(A₁)이 지수적으로 감쇠됨을 보여주기.
• 알베도 연산자 및 그 특이 부분 A₁의 연산자 노름을 수치적으로 계산하여 이론적 감쇠 및 폭발 행동 검증하기.
• Kn, σ_s 및 편미분 크기 z를 변화시키는 수치 실험을 수행하여 이론적 1/Kn 및 z에 대한 의존성 검증하기.

실험 결과

연구 질문

• RQ1 확산 스케일링 영역에서 RTE 역문제의 안정성은 쿤델 수치 Kn에 어떻게 의존하는가?
• RQ2 두 모델이 점근적으로 동일한데도 불구하고 RTE 역문제가 약간의 불안정성(리프시츠)으로 보이지만, 확산 방정식 역문제가 심각한 불안정성(로그)으로 보이는 이유는 무엇인가?
• RQ3 안정성 상수의 Kn에 대한 정확한 의존성은 무엇이며, 이는 불안정성 전이를 어떻게 설명하는가?
• RQ4 Kn→0에서 알베도 연산자의 특이 부분 A₁의 지수 감쇠는 수치적으로 검증될 수 있는가?
• RQ5 이론적으로 예측한 것처럼, 두 알베도 연산자 간의 차이 ∥A - Ã∥_* 는 소규모 Kn에서 e^{-c/ Kn} 비율로 감소하는가?

주요 결과

• RTE 역문제는 리프시츠 안정성을 보이지만, 리프시츠 상수가 e^{1/ Kn} 비율로 증가하여 소규모-Kn(확산) 영역에서 심각한 불안정성이 발생한다.
• 수치 결과는 알베도 연산자의 특이 부분 A₁의 연산자 노름이 1/ Kn 증가에 따라 지수적으로 감쇠됨을 확인하며, 이는 e^{-0.1/ Kn}과 일치한다.
• 고정된 편미분 크기 z=0.025일 때, 두 알베도 연산자 간의 차이 ∥A - Ã∥_* 는 e^{-0.05/ Kn} 비율로 감소함을 확인하여 이론적 지수 의존성을 검증한다.
• Kn→0일 때, 전체 알베도 연산자 A와 그 특이 부분 A₁ 간의 격차가 커지며, A₁이 A에 기여도가 줄어들게 되어 매개변수 분리가 어려워짐을 나타낸다.
• 수치 실험 결과, Kn이 1로 고정된 상태에서 ∥A - Ã∥_* 는 편미분 크기 z에 대해 선형 의존성을 보이며, 리프시츠 연속성 결과를 지지한다.
• 이론적 안정성 상수의 상한은 수치적으로 검증되었으며, 유한한 Kn 및 z 값 범위 내에서의 계산 제약에도 불구하고 성립한다.