Numerical solution of obstacle problem by using the mesh-free radial point interpolation method

Document Type : Research Article

Authors

1 Department of Mathematics, Isfahan (Khorasgan) branch, Islamic Azad University, Isfahan, Iran

2 Department of Mathematics, College of Computer and Mathematics, University of Thi-Qar, Iraq

Abstract

The obstacle problem is a well-known type of elliptic variational inequality that originally arose from contact issues in solid mechanics. The solutions to the obstacle problem often have irregularities along a free boundary, which can be effective in designing an appropriate numerical method for these problems.  In this paper, we present a  mesh-free method based on the radial point interpolation method for numerically solving an obstacle problem. In the proposed method, the  radial point interpolating shape functions are utilized in the global weak form of the obstacle problem, based on the element-free Galerkin method. This approach is combined with an active set strategy to address the obstacle problem. One of the key benefits of the proposed method is its independence from any mesh of the computing domain, along with its straightforward implementation and  high numerical stability. To ensure the efficiency of the presented method, we have investigated the convergence of  the proposed method.  The  obtained numerical results confirm the theoretical achievements and demonstrate the method's effectiveness and accuracy.

Keywords

Main Subjects

Journal of Mathematical Modeling
Volume 13, Issue 4
December 2025
Pages 929-942

Files

Share

How to cite

Statistics