Abstract:Two main aims of this paper are to develop a numerical method to solve an inverse source problem for parabolic equations and apply it to solve a nonlinear coefficient inverse problem. The inverse source problem in this paper is the problem to reconstruct a source term from external observations. Our method to solve this inverse source problem consists of two stages. We first establish an equation of the derivative of the solution to the parabolic equation with respect to the time variable. Then, in the second stage, we solve this equation by the quasi-reversibility method. The inverse source problem considered in this paper is the linearization of a nonlinear coefficient inverse problem. Hence, iteratively solving the inverse source problem provides the numerical solution to that coefficient inverse problem. Numerical results for the inverse source problem under consideration and the corresponding nonlinear coefficient inverse problem are presented.
What problem does this paper attempt to address?
The problems that this paper attempts to solve are mainly divided into two aspects:
1. **Inverse Source Problem**: Researchers have developed a numerical method to solve the inverse source problem of parabolic equations. Specifically, they hope to reconstruct the source term from external observation data. This problem is very important in practical applications, such as in the identification of pollution sources in rivers or lakes, electroencephalogram analysis in electrophysiology, and in fields such as ultrasonic imaging and photoacoustic tomography.
2. **Inverse Problem of Non - linear Coefficients**: By solving the above - mentioned inverse source problem, researchers further apply it to solve an inverse problem of non - linear coefficients. Specifically, they use an iterative method to gradually approximate the solution of the inverse problem of non - linear coefficients by repeatedly solving the inverse source problem.
### Problem Description
#### Inverse Source Problem
Given a positive number \(T\), assume that the function \(f(x,t)\) is known on the region \(\Omega\times[0, T]\), and \(f(x,0)\) is not zero at any point in \(\Omega\). The goal is to determine the function \(p(x)\) from the following measurement data:
\[G(x,t)=\partial_{\nu}u(x,t)\quad\text{for all}\;x\in\partial\Omega\;\text{and}\;t\in[0, T]\]
where \(u(x,t)\) satisfies the following initial - value problem:
\[
\begin{cases}
u_{t}(x,t) = Au(x,t)+f(x,t)p(x)&x\in\Omega,\;t > 0\\
u(x,t)=0&x\in\partial\Omega,\;t > 0\\
u(x,0)=0&x\in\Omega
\end{cases}
\]
#### Inverse Problem of Non - linear Coefficients
By solving the above - mentioned inverse source problem, researchers apply it to a more complex inverse problem of non - linear coefficients. Specifically, they consider the following problem:
\[F(x,t)=\partial_{\nu}u(x,t)\quad\text{for all}\;x\in\partial\Omega\;\text{and}\;t\in[0, T]\]
where \(u(x,t)\) satisfies the following parabolic equation:
\[
\begin{cases}
u_{t}(x,t)=\Delta u(x,t)+c(x)u(x,t)&x\in\Omega,\;t > 0\\
u(x,t)=g_{1}(x,t)&x\in\partial\Omega,\;t > 0\\
u(x,0)=g(x)&x\in\Omega
\end{cases}
\]
The goal is to determine the coefficient \(c(x)\).
### Solution Method
1. **Numerical Method for Inverse Source Problem**:
- Researchers first derive a new equation that contains only one unknown.
- Solve this equation by the quasi - reversibility method. The quasi - reversibility method is an effective method for solving ill - posed problems.
2. **Iterative Method for Inverse Problem of Non - linear Coefficients**:
- The initial guess \(c_{0}(x)\) is set as a constant function.
- Gradually approximate the solution of the inverse problem of non - linear coefficients by repeatedly solving the inverse source problem.
- The specific steps are as follows:
1. Set \(c_{0}\) as the background constant and calculate the solution \(u_{0}\) when \(c_{0}\) replaces \(c\).
2. Assume that \(c_{n}(x)\) and \(u_{n}(x,t)\) are known, and calculate the new Neumann data \(G_{n}(x,t)=F(x,t)-\partial_{n}u_{n}\)