Track new results across pure and applied mathematics.
Search papers, labs, and topics across Lattice.
The authors analyze the numerical approximation of steady-state heat conduction governed by a hemivariational inequality with non-smooth boundary conditions parametrized by a penalization factor $\alpha > 0$. They establish strong double-limit convergence of the finite element solutions to the constrained limit state as the mesh size $h \to 0$ and penalization parameter $\alpha \to \infty$ independently and simultaneously. Additionally, the analysis yields an optimal-order error estimate for the discrete limit problem while proving why uniform error bounds across arbitrary penalization values are fundamentally unachievable.
Mesh refinement and boundary penalization can be driven to their limits completely independently without sacrificing strong convergence in non-smooth thermal contact problems.
In this paper we study the numerical approximation and convergence analysis of a steady-state heat conduction problem with mixed boundary conditions. The physical model is governed by a hemivariational inequality depending on a heat transfer parameter $α> 0$. We consider the finite element approximation of this penalized problem, as well as its corresponding limit problem with a prescribed constant temperature on a part of the boundary. The main theoretical result establishes the strong convergence of the discrete solutions. Specifically, we prove the double-limit convergence of the finite element approximations to the limit solution as the mesh size $h$ tends to zero and the penalization parameter $α$ tends to infinity, independently and simultaneously. We further complement the convergence analysis with an error estimate of optimal order for the discrete limit problem, and we explain why an estimate uniform in the penalization parameter cannot be expected. The theoretical results are illustrated by numerical simulations on four examples verified by the method of manufactured solutions.