The Scalar Homotopy Method for Solving Non-Linear Obstacle Problem
Resource
Computers, Materials, & Continua, 15(1), 67-86
Journal
Computers, Materials, & Continua
Journal Volume
15
Journal Issue
1
Pages
67-86
Date Issued
2010
Date
2010
Author(s)
Fan, Chia-Ming
Liu, Chein-Shan
Yeih, Weichung
Chan, Hsin-Fang
Abstract
In this study, the nonlinear obstacle problems, which are also known as the nonlinear free boundary problems, are analyzed by the scalar homotopy method (SHM) and the finite difference method. The one- and two-dimensional nonlinear obstacle problems, formulated as the nonlinear complementarity problems (NCPs), are discretized by the finite difference method and form a system of nonlinear al- gebraic equations (NAEs) with the aid of Fischer-Burmeister NCP-function. Addi- tionally, the system of NAEs is solved by the SHM, which is globally convergent and can get rid of calculating the inverse of Jacobian matrix. In SHM, by introduc- ing a scalar homotopy function and a fictitious time, the NAEs are transformed to the ordinary differential equations (ODEs), which can be integrated numerically to obtain the solutions of NAEs. Owing to the characteristic of global convergence in SHM, the restart algorithm is adopted to fasten the convergence of numerical integration for ODEs. Several numerical examples are provided to validate the ef- ficiency and consistency of the proposed scheme. Besides, some factors, which might influence on the accuracy of the numerical results, are examined by a series of numerical experiments.
Type
journal article
File(s)![Thumbnail Image]()
Loading...
Name
251.pdf
Size
23.21 KB
Format
Adobe PDF
Checksum
(MD5):dbe00225008383e93f97efa682b1ce88