Indian Journal of Science and Technology
Year: 2016, Volume: 9, Issue: 36, Pages: 1-6
S. M. Afzal Hoq1,3*, Erwin Sulaeman1 and Abdurahim Okhunov2
1 Department of Mechanical Engineering, International Islamic University, Malaysia; [email protected]
2 Department of Science and Engineering, International Islamic University, Malaysia; [email protected]
3 Department of Civil Engineering, Southern University Bangladesh, Bangladesh; [email protected]
*Author for correspondence
Department of Mechanical Engineering
The present work explores an error analysis of Galerkin finite element method (GFEM) for computing steady heat conduction in order to show its convergence and accuracy. The steady state heat distribution in a planar region is modeled by two-dimensional Laplace partial differential equations. A simple three-node triangular finite element model is used and its derivation to form elemental stiffness matrix for unstructured and structured grid meshes is presented. The error analysis is performed by comparison with analytical solution where the difference with the analytical result is represented in the form of three vector norms. The error analysis for the present GFEM for structured grid mesh is tested on heat conduction problem of a rectangular domain with asymmetric and mixed natural-essential boundary conditions. The accuracy and convergence of the numerical solution is demonstrated by increasing the number of elements or decreasing the size of each element covering the domain. It is found that the numerical result converge to the exact solution with the convergence rates of almost O(h²) in the Euclidean L2 norm, O(h²) in the discrete perpetuity L∞norm and O(h1 ) in the 1 H norm.
Keywords: Error Analysis, Finite Element Method, Galerkin’s Weight Residual Approach, Heat Conduction, Laplace Equation, Partial Differential Equation
Subscribe now for latest articles and news.