Vol.8, No.2, 2008, pp.53-66, doi:10.3970/cmc.2008.008.053
OPEN ACCESS
ARTICLE
A Time-Marching Algorithm for Solving Non-Linear Obstacle Problems with the Aid of an NCP-Function
  • Chein-Shan Liu1
Department of Mechanical and Mechatronic Engineering, Computation and Simulation Center, Taiwan Ocean University, Keelung, Taiwan. E-mail: csliu@mail.ntou.edu.tw
Abstract
Proposed is a time-marching algorithm to solve a nonlinear system of complementarity equations: Pi(xj) ≥ 0, Qi(xj) ≥ 0 , Pi(xj)Qi(xj) = 0, i, j = 1,...,n, resulting from a discretization of nonlinear obstacle problem. We transform the above nonlinear complementarity problem (NCP) into a nonlinear algebraic equations (NAEs) system: Fi(xj) = 0 with the aid of the Fischer-Burmeister NCP-function. Such NAEs are semi-smooth, highly nonlinear and usually implicit, being hard to handle by the Newton-like method. Instead of, a first-order system of ODEs is derived through a fictitious time equation. The time-stepping equations are obtained by applying a numerical integration on the resultant ODEs, which are derivative-free and do not need the inverse of any matrix. The computational cost is thus greatly reduced. The numerical examples of Bratu, von Karman and other elliptic equations are used to demonstrate that the new fictitious time integration method (FTIM) is highly efficient to calculate the obstacle problems.
Keywords
Nonlinear Obstacle Problem, Nonlinear complementarity problem, Nonlinear algebraic equations, Iterative method, Elliptic equations, Fictitious time integration method (FTIM)
Cite This Article
Liu, C. (2008). A Time-Marching Algorithm for Solving Non-Linear Obstacle Problems with the Aid of an NCP-Function. CMC-Computers, Materials & Continua, 8(2), 53–66.