The quasi-reversibility method to numerically solve an inverse source problem for hyperbolic equations (Q2031868)

The inverse source problem of computing the initial condition of the following hyperbolic equations from the measurements of Cauchy data are considered \[ \begin{gathered} \mathbf{a}(\mathbf{x})u_{tt} (\mathbf{x},t)+a(\mathbf{x})u_t (\mathbf{x},t)=\Delta u(\mathbf{x},t)+ \mathbf{b}\cdot\nabla u(\mathbf{x},t)+ \mathbf{c}u(\mathbf{x},t), \quad (\mathbf{x},t) \in \Omega \times (0,T),\\ u_t (\mathbf{x},0)=0, \quad \mathbf{x}\in \Omega, \quad u(\mathbf{x},0)=p(\mathbf{x}), \quad \mathbf{x}\in \Omega, \end{gathered} \] with either homogeneous Dirichlet \(u(\mathbf{x},t)=0\), \((\mathbf{x},t)\in\partial\Omega\times [0,T]\) or Neumann \(u_ν (\mathbf{x},t)=0\), \((\mathbf{x},t)\in\partial\Omega\times [0,T]\) boundary conditions. Here \(\Omega\) is a smooth and bounded domain in \(\mathbb{R}^d\), \(d\geq 1\), and \(T\) is a positive number; \(1\leq \mathbf{a}\in C^1(\Omega)\), \(a, \mathbf{c}\in L^\infty (\Omega)\), \(b\in L^\infty (\Omega, \mathbb{R}^d)\). The authors consider problems of solving \(u(\mathbf{x},t)\in H^2 (\Omega\times (0, T))\) generated by the source \(p(\mathbf{x})\in C_0^2 (\Omega)\) and subjected to the above-mentioned Dirichlet or Neumann boundary condition. The paper is organized as follows. Section 1 is an introduction. In Section 2, the inverse problems under consideration are stated and an approximation model whose solution directly yields their solutions is derived. In Section 3, some auxiliary results are introduced and the Carleman estimate is proved. In Section 4, the quasi-reversibility method to solve the system of elliptic equations is implemented and the convergence of the solution as the noise level tends to 0 is proved. Numerical experiments with figures are given in Section 5. Section 6 is to provide some numerical results for the 1D problem and to show the significance of the used orthonormal basis. Finally, some conclusions are fixed in Section 7.