Abstract

In this paper, we consider two methods, the Second order Central Difference Method (SCDM) and the Finite Element Method (FEM) with P1 triangular elements, for solving two dimensional general linear Elliptic Partial Differential Equations (PDE) with mixed derivatives along with Dirichlet and Neumann boundary conditions. These two methods have almost the same accuracy from theoretical aspect with regular boundaries, but generally Finite Element Method produces better approximations when the boundaries are irregular. In order to investigate which method produces better results from numerical aspect, we apply these methods into specific examples with regular boundaries with constant step-size for both of them. The results which obtained confirm, in most of the cases, the theoretical results.

Highlights

  • Finite Difference schemes and Finite Element Methods are widely used for solving partial differential equations [1]

  • Finite Element methods are more complicated than Finite Difference schemes because they use various numerical methods such as interpolation, numerical integration and numerical methods for solving large linear systems

  • We contacted a numerical study with Matlab R2015a. We apply these methods into specific elliptical problems, in order to test which of these methods produce better approximations when the Dirichlet and Neumann boundary conditions are imposed

Read more

Summary

Introduction

Finite Difference schemes and Finite Element Methods are widely used for solving partial differential equations [1]. We will describe the Second order Central Difference Scheme and the Finite Element Method for solving general second order elliptic partial differential equations with regular boundary conditions on a rectangular domain. For both of these methods, we consider the Dirichlet and Neumann Boundary conditions, along the four sides of the rectangular area. We contacted a numerical study with Matlab R2015a We apply these methods into specific elliptical problems, in order to test which of these methods produce better approximations when the Dirichlet and Neumann boundary conditions are imposed.

Second Order Central Difference Scheme
B2 D2 O O
Finite Element Method
Bilinear Interpolation in P1
Gauss Quadrature
Finite Element Algorithm
Error Analysis
Numerical Study
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call