Abstract:In this work, we derive a reliable and efficient residual-typed error estimator for the finite element approximation of a 2d cathodic protection problem governed by a steady-state diffusion equation with a nonlinear boundary condition. We propose a standard adaptive finite element method involving the Dörfler marking and a minimal refinement without the interior node property. Furthermore, we establish the contraction property of this adaptive algorithm in terms of the sum of the energy error and the scaled estimator. This essentially allows for a quasi-optimal convergence rate in terms of the number of elements over the underlying triangulation. Numerical experiments are provided to confirm this quasi-optimality.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in the cathodic protection problem, how to effectively use the Adaptive Finite Element Method (AFEM) to solve two - dimensional problems described by the steady - state diffusion equation and nonlinear boundary conditions, and ensure the quasi - optimal convergence rate of the algorithm.
Specifically, the paper focuses on the following partial differential equation system:
\[
-\nabla \cdot (\sigma \nabla u) = 0 \quad \text{in } \Omega,
\]
\[
\sigma \frac{\partial u}{\partial n} =
\begin{cases}
0 & \text{on } \Gamma_0, \\
g & \text{on } \Gamma_A, \\
-f(u) & \text{on } \Gamma_C,
\end{cases}
\]
where $\Omega$ is a two - dimensional bounded polygonal region, and the boundary $\Gamma$ is divided into three mutually disjoint parts: $\Gamma_0$, $\Gamma_A$ and $\Gamma_C$. The conductivity $\sigma$ is a piecewise $W^{1,\infty}$ function, and satisfies $\sigma_1 \leq \sigma \leq \sigma_2$. The boundary conditions include insulation conditions, current density conditions, and nonlinear conditions in the corrosion process.
Due to the sudden change of boundary conditions and the discontinuity of conductivity, singularities in local solutions are caused. Traditional uniform mesh refinement will lead to low computational efficiency. Therefore, the paper proposes an Adaptive Finite Element Method, which improves the accuracy of numerical results by locally refining the mesh and minimizes the number of degrees of freedom as much as possible.
The main contributions include:
1. A reliable a posteriori error estimator is proposed.
2. The contraction property is established, that is, the sum of the energy error and the scaled estimator will decrease after each iteration.
3. The quasi - optimal complexity of the algorithm is proved, that is, as the number of elements increases, the error converges at the optimal rate.
These results are verified by theoretical analysis and numerical experiments, indicating that the proposed Adaptive Finite Element Method is efficient and accurate in dealing with cathodic protection problems.