Abstract:The analyses of interior penalty discontinuous Galerkin methods of any order k for solving elliptic and parabolic problems with Dirac line sources are presented. For the steady state case, we prove convergence of the method by deriving a priori error estimates in the L2 norm and in weighted energy norms. In addition, we prove almost optimal local error estimates in the energy norm for any approximation order. Further, almost optimal local error estimates in the L2 norm are obtained for the case of piecewise linear approximations whereas suboptimal error bounds in the L2 norm are shown for any polynomial degree. For the time-dependent case, convergence of semi-discrete and of backward Euler fully discrete scheme is established by proving error estimates in L2 in time and in space. Numerical results for the elliptic problem are added to support the theoretical results.
What problem does this paper attempt to address?
This paper attempts to address the challenges encountered in the numerical approximation using the interior penalty discontinuous Galerkin (dG) method for elliptic and parabolic problems with Dirac line sources. Specifically, the paper focuses on how to achieve effective numerical solutions in these special types of partial differential equations and ensure the convergence of solutions and error estimates.
### Main problems
1. **Treatment of non - smooth solutions**:
- The true solutions in elliptic and parabolic problems are not spatially smooth enough, especially in the case of including Dirac line sources. This causes the solutions to exhibit singular behavior in certain regions (such as near the line source), making it difficult for traditional finite element methods to obtain ideal error estimates.
2. **Global and local error estimates**:
- For elliptic problems, the author aims to prove the convergence in the global \(L^2\) and weighted energy norms and derive almost optimal local \(L^2\) error estimates (for linear polynomials) as well as sub - optimal error bounds for any polynomial order.
- For parabolic problems, the author hopes to prove the global convergence of the semi - discrete and fully - discrete schemes in the \(L^2(0,T;L^2(\Omega))\) norm for time - dependent problems.
3. **Validation of the effectiveness of numerical methods**:
- The paper validates the theoretical results through numerical experiments, especially for the numerical results of elliptic problems.
### Specific contributions
- **Elliptic problems**:
- Proved the convergence in the global \(L^2\) and weighted energy norms.
- In the region excluding the line source \(\Lambda\), derived almost optimal \(L^2\) error estimates (for linear polynomials) and sub - optimal error bounds for any polynomial order.
- Established an almost optimal error rate in the local energy norm for any polynomial order.
- **Parabolic problems**:
- Proved the global convergence of the semi - discrete and backward Euler fully - discrete schemes in the \(L^2(0,T;L^2(\Omega))\) norm.
### Mathematical formula representation
- The true solution \(u\) satisfies the following elliptic model problem:
\[
-\Delta u = f\delta_\Lambda, \quad \text{in } \Omega,
\]
\[
u = 0, \quad \text{on } \partial\Omega.
\]
where \(f\in L^2(\Lambda)\), and \(f\delta_\Lambda\) is defined as:
\[
\langle f\delta_\Lambda, v\rangle=\int_\Lambda fv \, ds, \quad \forall v\in L^\infty(\Omega).
\]
- The parabolic model problem is as follows:
\[
\frac{\partial u}{\partial t}-\Delta u = f\delta_\Lambda, \quad \text{in } \Omega\times(0, T],
\]
\[
u = 0, \quad \text{on } \partial\Omega\times(0, T],
\]
\[
u = u_0, \quad \text{in } \{0\}\times\Omega.
\]
Through these studies, the paper not only solves the numerical approximation problems of elliptic and parabolic problems with Dirac line sources but also provides a solid theoretical foundation for future research.