[논문 리뷰] Inverse Problem of Finding the Time-dependent Coefficient of Heat Equation from Integral Overdetermination Condition Data

Mansur I. Ismailov, Fatma Kanca|arXiv (Cornell University)|2010. 04. 30.

Numerical methods in inverse problems참고 문헌 7인용 수 14

한 줄 요약

이 논문은 비국소 경계 조건과 적분 과다조건을 갖는 1차원 열방정식에서 시간에 따라 변하는 열확산계수 계수 $ a(t) $와 온도 분포 $ u(x,t) $를 식별하는 역문제를 다룬다. 푸리에 방법과 고정점 원리를 사용하여 특정 부드러움 및 호환성 조건 하에서 고전적 해의 존재성과 유일성을 증명하고, 크랭크-니콜슨 기반의 수치적 스킴과 예측-보정 반복법을 제안하여 역문제를 해결한다. 수치 예제를 통해 높은 정확도로 검증되었다.

ABSTRACT

In this paper we consider the problem of simultaneously determining the time-dependent thermal diffusivity and the temperature distribution in one-dimensional heat equation in the case of nonlocal boundary and integral overdetermination conditions. We establish conditions for the existence and uniqueness of a classical solution of the problem under considerations. We present some results on the numerical solution with an example.

연구 동기 및 목표

• 비국소 경계 조건과 적분 과다조건을 갖는 1차원 열방정식에서 $ a(t) $와 $ u(x,t) $를 식별하는 역문제에 대해 고전적 해의 존재성과 유일성을 보장하는 조건을 설정하는 것.
• 해의 연속적 의존성(문제 데이터에 대한 민감도)을 분석하는 것.
• 크랭크-니콜슨 스킴과 반복적 예측-보정 방법을 기반으로 한 수치 알고리즘을 개발하는 것.
• 정확한 해가 알려진 수치 예제를 통해 방법을 검증하고 수렴성과 정확도를 입증하는 것.

제안 방법

• 공간 영역에서 고유함수 전개를 활용하여 문제를 고정점 문제로 변환한다.
• 비국소 경계 조건과 관련된 스펙트럼 문제를 분석하고, 고유함수와 관련 함수의 이중수체계를 구성한다.
• 해는 고유함수 기반의 급수 전개로 표현되며, 시간에 따라 변하는 계수 $ a(t) $는 적분 과다조건 $ \int_0^1 u(x,t)dx = E(t) $ 로부터 유도된다.
• 예측-보정 반복법을 설계한다: 계수 $ a^{j+1} $ 는 $ E(t) $ 의 시간 도함수와 $ u $ 의 공간 기울기로부터 예측되고, 이후 크랭크-니콜슨 스킴을 통해 보정된다.
• 크랭크-니콜슨 스킴은 공간과 시간에 대해 두 번째 차수 정확도를 갖도록 적용되며, 가우스 소거법을 사용하여 반복적으로 해를 구한다.
• 수치적 구현은 $ h = 0.005 $, $ \tau = h/4 $ 인 유한차분 이산화를 사용하고, 수렴성은 연속된 반복 간의 차이를 모니터링하여 확인한다.

실험 결과

연구 질문

• RQ1 비국소 경계 조건과 적분 과다조건을 갖는 1차원 열방정식에서 시간에 따라 변하는 열확산계수 $ a(t) $와 온도 $ u(x,t) $를 식별하는 역문제가 고전적 해를 갖는 조건은 무엇인가?
• RQ2 해의 존재성과 유일성을 보장하기 위해 초기 조건, 소스 항, 과다조건 함수에 대해 필요한 부드러움 조건과 호환성 조건는 무엇인가?
• RQ3 $ \int_0^1 u(x,t)dx = E(t) $ 조건으로부터 시간에 따라 변하는 계수 $ a(t) $ 를 어떻게 재구성할 수 있으며, 이 재구성의 안정성은 어떻게 평가할 수 있는가?
• RQ4 비국소 조건이 있는 역문제를 해결할 때 두 번째 차수 정확도와 안정성을 보장하는 수치적 스킴은 무엇인가?
• RQ5 제안된 반복 수치 방법이 정확한 해와 비교하여 $ a(t) $ 와 $ u(x,t) $ 를 얼마나 정확하게 복원하는가?

주요 결과

• 가정 (A1)–(A3) 하에, 즉 데이터에 대한 부드러움, 호환성, 양성 조건이 만족될 경우, 역문제는 $ Q_T $ 에서 고전적 해를 갖는다.
• 해는 어떤 $ T_0 < T $ 에 대해 $ Q_{T_0} $ 에서 유일하게 결정되며, $ T_0 $ 는 문제의 데이터에 따라 달라진다.
• 크랭크-니콜슨과 반복적 예측-보정 스킴을 기반으로 한 수치 방법은 $ a(t) $ 와 $ u(x,t) $ 를 높은 정확도로 성공적으로 복원한다.
• 테스트 예제에서 $ T = 1/4 $, $ h = 0.005 $, $ \tau = h/4 $ 일 때, 정확한 해 $ a(t) = \frac{1}{(2\pi)^2} + \exp(4t) $ 와 $ u(x,t) = (1-x)\sin(2\pi x)\exp(-t) $ 와의 일치가 높게 유지된다.
• 반복 스킴은 수렴하며, 각 시간 단계에서의 $ a(t) $ 와 $ u(x,t) $ 의 계산값은 연속된 반복 간의 차이가 사전에 정해진 허용 오차 이내일 경우 수락된다.
• 표 1과 2의 수치 결과로 확인되듯이, 방법은 안정성과 두 번째 차수 정확도를 보여준다.