Application of Fem and FDM in Solving 2D Irregular Geometry Heat Transfer Problem
Application of Fem and FDM in Solving 2D Irregular Geometry Heat Transfer Problem
Application of Fem and FDM in Solving 2D Irregular Geometry Heat Transfer Problem
Faculty of Science
JUNE 2014
iii
Dr Yeak Su Hoe
Thank you
iv
ACKNOWLEDGEMENT
In the name of Allah The Most Beneficent and The Most Merciful. The ultimate praises
is only to Him for giving me this chance to complete my dissertation.
I would love to dedicate my greatest indebtedness to the most important person who is
my supervisor for being the backbone of this dissertation. His advice and guidance are the things
that keep me going until the very end. Also my gratitude to all my lecturers whose deeds and
kindness I will never forget.
Million thanks to my family for the endless support and care, for always lending me a
helpful hand and for staying steadfastly by my side no matter what.
Last but not least a bouquet of appreciation to all my thoughtful friends for the
encouragement and assistance throughout the ups and downs of this challenging yet inspiring
journey.
.
v
ABSTRACT
The study will focus on the application of finite element method (FEM) in solving two
dimensional irregular geometry heat transfer problem. FEM is known as one of the numerical
technique for finding approximate solutions to boundary value problems for differential
equation. Therefore this method is a perfect application to solve a parabolic partial differential
equation such as heat equation that describes the distribution of heat in a given region over time.
In order to minimize an error function and produce a stable solution, FEM uses variation method
or also known as the calculus of variations. All the solutions to this method related to the given
problem will be compared with the solution from another popular numerical method which is
finite difference method (FDM). FDM is also very useful in approximating the solutions to
differential equation. However the application is limited to regular geometry and simple irregular
geometry problems. The comparison of the solutions from these two methods will confirm that
FEM is a better choice in solving two dimensional irregular geometry problems involving heat
transfer.
vi
ABSTRAK
Kajian ini akan bertumpu kepada penggunaan kaedah unsur terhingga untuk
menyelesaikan masalah yang melibatkan perpindahan haba dalam geometri dua dimensi yang
tidak sekata. Kaedah unsur terhingga adalah salah satu kaedah berangka bagi mencari
penyelesaian anggaran untuk masalah nilai sempadan bagi persamaan pembezaan. Maka dengan
itu kaedah ini adalah pilihan yang tepat untuk menyelesaikan masalah berkaitan perpindahan
haba kerana persamaannya turut melibatkan nilai sempadan untuk persamaan pembezaan. Bagi
meminimumkan kesalahan fungsi and mendapatkan penyelesaian yang lebih stabil, kaedah unsur
terhingga ini menggunakan kaedah variasi atau turut dikenalai sebagai kalkulus variasi. Semua
penyelesaian terhadap masalah yang diberi akan dibandingkan dengan penyelesaian daripada
satu lagi kaedah berangka yang popular iaitu kaedah perbezaan terhingga. Kaedah perbezaan
terhingga juga sangat berguna untuk mendapatkan anggaran jawapan kepada persamaan
pembezaan. Tetapi aplikasi kaedah ini limit hanya untuk dimensi yang sekata dan dimensi tidak
sekata yang mudah sahaja. Perbandingan antara dua kaedah ini akan membuktikan bahawa
kaedah unsur terhingga merupakan pilihan ynag lebih baik untuk menyelesaikan masalah
perpindahan haba dalam geometri dua dimensi yang tidak sekata.
vii
TABLE OF CONTENTS
DECLARATION ii
DEDICATION iii
ACKNOWLEDGEMENT iv
ABSTRACT v
ABSTRAK vi
LIST OF TABLES xi
LIST OF SYMBOLS xv
1 INTRODUCTION
1.0 Introduction 1
2 LITERATURE REVIEW
2.0 Introduction 5
3 METHOD OF SOLUTION
3.0 Introduction 17
23
3.3 2D Simple Irregular Geometry Heat Transfer Problem
4.0 Introduction 35
5.0 Introduction 46
5.1 Conclusion 46
5.2 Recommendation 47
REFERENCES
APPENDIX
xi
LIST OF TABLES
problem.
FDM.
LIST OF FIGURES
4.15 Heat flux for simple irregular geometry problem using FEM 43
when f = 10.
4.17 Heat flux for simple irregular geometry problem using FEM 43
when f = -10.
4.19 Heat flux for simple irregular geometry problem using FEM 44
when f = 100.
4.21 Heat flux for simple irregular geometry problem using FEM 45
when f = -100.
xv
LIST OF SYMBOLS
delta
Te displacement vector
Fe load vector
Q heat generation
Ei truncation error
pi
n normal flux
T temperature
ST specified temperature
Sc convection
xvi
LIST OF ABBREVIATIONS
2D Two Dimensional
INTRODUCTION
1.0 Introduction
Finite Element Method (FEM) and Finite Difference Method (FDM) are the examples
of numerical methods that can be used to solve this kind of problems. FEM is a numerical
method for solving problems of engineering and mathematical physics. It is useful for
problems with complicated geometries, loadings and material properties where analytical
solutions cannot be obtained. It is a numerical technique in finding approximate solutions to
boundary value for differential equations problems of engineering and mathematical physics.
Heat transfer is the study of thermal energy transport within a medium or among
neighbouring media by molecular interaction, fluid motion, and electro-magnetic waves,
resulting from a spatial variation in temperature (Long and Seyma, 2009). This variation in
temperature is governed by the principle of energy conservation, which when applied to a
control volume or a control mass, states that the sum of the flow of energy and heat across the
system, the work done on the system, and the energy stored and converted within the system,
is zero. Heat transfer finds application in many important areas, namely design of thermal and
2
nuclear power plants including heat engines, steam generators, condensers and other heat
exchange equipment, catalytic convertors, heat shields for space vehicles, heaters, electronic
equipment, internal combustion engines, refrigeration and air conditioning units, design of
cooling systems for electric motors generators and transformers, heating and cooling of fluids
in chemical operations, construction of dams and structures, minimization of building heat
losses using improved insulation techniques, thermal control of space vehicles, heat treatment
of metals and dispersion of atmospheric pollutants.
Basically, the two-dimensional heat equation can be solved theoretically and also
numerically by using numerical method such as the finite difference method. Heat equation is
one of the most important partial differential equation which describes the distribution of heat
(or variation in temperature) in a given region over time.
Grid or mesh based numerical methods such as finite difference methods (FDM) and
finite element methods (FEM) have been widely applied to various areas such as fluid
dynamics and solid mechanics. These methods are currently the dominant methods in solving
problems of engineering and science.
Problems related to heat transfer can be solve by using finite difference method
(FDM). However this approach will cause a few flaws in the findings mainly because the
method itself is low in accuracy. Although FDM is easier to compute as well as to code, some
problems that involved complex geometry is either difficult or completely cannot be solved
3
by FDM. Generally FDM is a simple method to use for common problems defined on regular
geometries.
The two sources of error in FDM are round-off error which is the loss of precision due
to computer rounding of decimal quantities and the local truncation error typically expressed
using Big-O notation that refers to the error from a single application of the method. Besides,
the majority of heat transfer problem is related to engineering and the problem domain is in
irregular geometry. It is not ideal to solve this kind of problems using FDM since it is only
applicable for problems that involve regular geometry and simple irregular geometry.
Even though FEM is a complex method but it can be used to solve complex geometry
problems. Industrial or engineering heat transfer problems normally involve irregular
geometry domain. The need for numerical methods arises from the fact that analytical
solutions does not exist for most practical problems.
On the other hand FEM uses variation methods to minimize an error function and
produce a stable solution. Therefore this method is mostly ideal for problems with
complicated geometries when analytical solutions cannot be obtained. For that reason, this
research will focus on how to solve heat transfer problems in twodimensional irregular
geometry by using the application of FEM.
2. To associate the application of FEM and FDM in solving problems related to heat
transfer.
This research will focus on the application of FEM with the approach of Galerkins
method to solve heat transfer problems in simple 2D irregular geometries while comparing
the solutions with FDM.
The results of this research can be used for further research in related areas. One of
the goals of this research is that, this research will gives a clearer view of FEM, where it will
finally lead to further understanding on applying the FEM in solving related problems
especially in the science and engineering fields.
In this report of dissertation, there will be five chapters which are introduction,
literature review, method of solution, results and discussion and conclusion and
recommendation.
As for Chapter 1 there will be seven subtopics which are introduction, background of
study, problem statement, objective of the study, scope of the study, significant of the study
and organization of report.
Chapter 2 is going to introduce a literature review on FEM and FDM also its
application in solving heat transfer problem.
The formulation and the methods that will be used are discussed in Chapter 3. This
chapter provides the modeling and derivation of FEM, FDM and heat transfer equation.
Then all the results that are obtained from Matlab will be discussed in Chapter 4. The
results will be presented in figures and table.
Finally Chapter 5 is for the conclusion of this dissertation and recommendation for
further research.
REFERENCES
Reimer, A. S. and Cheviakov, A. F. (2012). A Matlab-Based Finite Difference Solver for the
Poisson Problem with Mixed Dirichlet-Neumann Boundary Conditions. Department of
Mathematics and Statistics, University of Saskatchewan, Saskatoon S7N 5E6 Canada.
Roos C. (2008). Principles of Heat Transfer. Washington State University Extension Energy
Program.
Pepper, D. W. and Heinrich, J. C. (1992). The Finite Element Method Basic Concepts and
Applications. Hemisphere Publishing Corporation.
Peiro, J. and Sherwin, S. (2010). Finite Difference, Finite Element and Finite Volume
Methods for Partial Difference Equations. Department of Aeronautics, Imperial College.
London. UK.
Larson, G. M. and Bengzon, F. (2010). The Finite Element Method: Theory, Implementation,
and Practice. Springer.
49
Nor Hafizah Ahmad Kailani (2014), The Application of Finite Element Method in 2D Heat
Distribution Problem for Irregular Geometry. Thesis Master. Universiti Teknologi Malaysia.
John C. Chai, Girija Parthasarathy, HaeOk S. Lee, Suhas V. Patankar (1995). Finite Volume
Radiative Heat Transfer Procedure for Irregular Geometries. Journal of Thermophysics and
Heat Transfer, Vol. 9, No. 3 (1995), pp. 410-415.
Susane C. Brenner, L. Ridgway Scott (2000). The Mathematical Theory of Finite Element
Methods: Third Edition. New York: Springer Science + Business Media.
Zabaras, N. (2012). Finite Element Analysis for Mechanical and Aerospace Design. New
York: Cornell University.
Elden, L. (1997). Solving Inverse Heat Conduction Problem by a Method of Lines. Journal
of Heat Transfer 119(3), 406-412 J. Heat Transfer 119(3), 406-412.
Ash, J., Baldwin, J., Hirt, J. Lance, A. (2005). Derivation of the Heat Conduction Equation.
www.math.vt.edu.
Nayroles, B., Touzot, G., Villon, P. (1992). Generalizing the finite element method: Diffuse
approximation and diffuse elements. Computational Mechanics Springer.
Logan, D.L. (2007). A First Course in the Finite Element Method Fourth Edition. Canada
Thomson.