About this Journal Submit a Manuscript Table of Contents
Advances in Numerical Analysis
Volume 2013 (2013), Article ID 967342, 10 pages
http://dx.doi.org/10.1155/2013/967342
Research Article

Numerical Method for Solving Matrix Coefficient Elliptic Equation on Irregular Domains with Sharp-Edged Boundaries

1Department of Mathematics, College of Science, China University of Petroleum, Beijing 102249, China
2Department of Mathematics and Statistics, Louisiana Tech University, Ruston, LA 71272, USA

Received 2 June 2013; Accepted 21 November 2013

Academic Editor: Yanping Lin

Copyright © 2013 Liqun Wang and Liwei Shi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

We present a new second-order accurate numerical method for solving matrix coefficient elliptic equation on irregular domains with sharp-edged boundaries. Nontraditional finite element method with non-body-fitting grids is implemented on a fictitious domain in which the irregular domains are embedded. First we set the function and coefficient in the fictitious part, and the nonsmooth boundary is then treated as an interface. The emphasis is on the construction of jump conditions on the interface; a special position for the ghost point is chosen so that the method is more accurate. The test function basis is chosen to be the standard finite element basis independent of the interface, and the solution basis is chosen to be piecewise linear satisfying the jump conditions across the interface. This is an efficient method for dealing with elliptic equations in irregular domains with non-smooth boundaries, and it is able to treat the general case of matrix coefficient. The complexity and computational expense in mesh generation is highly decreased, especially for moving boundaries, while robustness, efficiency, and accuracy are promised. Extensive numerical experiments indicate that this method is second-order accurate in the norm in both two and three dimensions and numerically very stable.