Numerical Method for Solving Matrix Coefficient Elliptic Equation on Irregular Domains with Sharp-Edged BoundariesReportar como inadecuado

Numerical Method for Solving Matrix Coefficient Elliptic Equation on Irregular Domains with Sharp-Edged Boundaries - Descarga este documento en PDF. Documentación en PDF para descargar gratis. Disponible también para leer online.

Advances in Numerical AnalysisVolume 2013 2013, Article ID 967342, 10 pages

Research Article

Department of Mathematics, College of Science, China University of Petroleum, Beijing 102249, China

Department 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.


We present a new second-order accurate numerical method for solving matrix coefficient elliptic equation on irregular domains withsharp-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 expensein mesh generation is highly decreased, especially for moving boundaries, while robustness, efficiency, and accuracy are promised. Extensive numericalexperiments indicate that this method is second-order accurate in the norm in both two and three dimensions and numerically very stable.

Autor: Liqun Wang and Liwei Shi



Documentos relacionados