Finite element method solution of Poisson's equation is studied for three dimensional semiconductor devices. Newton-Raphson iteration method is used to solve nonlinear equation, while linear matrix equation for each iteration step is solved by Jacobi Conjugate Gradient Method(JCGM) for numerical efficiency. Also automatic mesh refinements are incorporated with local error estimation. Delaunay tetrahedral algorithm is used for initial mesh generation and following mesh refinements.