One Dimensional Heat Transfer through a Uniform Plane Wall by Using Finite Volume Method

In this paper, finite volume method has been used to investigate one dimensional conductive heat transfer throw a uniform plane wall. Then the step by step procedures of this numerical solution are described and implemented in a real world problem where tri-diagonal matrix algorithm and Gaussian elimination matrix method are applied to solve the system of our discretized algebraic system of equations. Finally, to check the accuracy of our method, a comparison between the numerical solution obtained by finite volume techniques and exact solution is presented which show a minimum error compared to other existing methods.


INTRODUCTION:
Nowadays, to investigate the mass, momentum, and heat transfer phenomena, computational fluid dynamics (CFD) study numerical solutions which is associated with the phenomena of thermodynamics, turbo machinery, chemical manufacturing, power generation, weather simulation, biological engineering, meteorology, aerospace, reaction chemistry, predicting fluid flow and heat transfer . Though the fundamental source of almost all CFD problems is Navier-Stokes equations which describe many single-phase (liquid or gas, but not both) fluid flows.
Recently, due to the revolution of computer technology, abundant computational grid techniques have been developed which is very efficient to solve numerous engineering problems (Cheniguel and Reghioua, 2013;Chuathong and Toutip, 2011). Among these numerical grid techniques, finite difference method (FDM) is used as a common numerical technique to solve numerous engineering problems (Lau and Kuruganty, 2010). Again, finite element method (FEM) is another kind of commonly used numerical technique which has been applied to solve many heat transfer problems (Patil and Prasad, 2013). A generalized transfer equation for a dependent variable which can be mass, concentration, heat and momentum is given by Patankar (1980), where discretization technique is applied for CFD analysis (Islam and Hossain, 2019).
Moreover, the finite volume method (FVM) is becoming most popular numerical technique which is commonly used in commercial CFD software such as COMSOL Multiphysics (Versteeg and Malalasekera, 2007). Though to solve any problem these three methods have its own merits besides its demerits, FVM technique is one of the utmost flexible and multipurpose techniques to solve CFD (Uddin et al., 2020). In this paper, we have described an engineering problem by the point of view of FVM. Due to its wide-spread popularity in CFD, this method has been applied previously in different mechanical engineering problem as like thermo elastic and linear elastic problems (Berezovski and Maugin, 2001;Demirdžić and Muzaferija, 1994;and Jasak and Weller, 2000).
The rest of our paper is prepared as follows. A short review of FVM with the help of Tri-Diagonal Matrix Algorithm TDMA is given in Section 2. In Section 3, the numerical solutions obtained by this technique is given where a comparison between exact and our numerical solutions is also described. Finally, in Section 5 we have concluded this paper.

FINITE VOLUME METHOD:
Consider the steady state diffusion of a property in a one-dimensional domain defined in Fig 1. The process is governed by Where is the thermal connectivity, is the temperature and is the source term (i.e., the rate of heat generation per unit volume).

Grid Generation
In the finite volume method first of all, we have to divide the domain into discrete control volumes. Let us consider a number of nodal points in the domain space from to . The boundaries (or faces) of control volumes are positioned mid-way between adjacent nodes. Thus each node is surrounded by a control volume or cell. Today it is common practice to set up control volumes near to the edge of the domain in such a way where the physical boundaries coincide with the control volume boundaries. At this point it is appropriate to establish a system of notation that can be used in future developments. A general nodal point is identified by and its neighbor's nodes to the north and south are represented by and respectively. The south side control volume face is mentioned by and the north side face of the control volume is mentioned to by . The distances between the nodes and S, and between nodes N and are represented by and respectively. In the same way, distances between point and face s is denoted by and distances between face n and point is denoted by .

Discretization
To discretize the governing equation, we will integrate the governing equation (1) over the control volume at its nodal point P.
Then we get, Where the cross-sectional area of the control volume face is, ∆ is the volume and is the source term over the control volume. This is a very attractive feature of this finite volume method where the discretized equation contains a clear physical interpretation. Equation 2 states that the diffusive flux of leaving the south face minus the diffusive flux of entering the north face is equal to the generation of over the control volume. In order to derive useful forms of the discretized equations, the interface diffusion coefficient and the gradient at west and east are required. And the diffusive flux terms are evaluated as - For practical situations, the source term may be a function of the dependent variable. By substituting the values of equations (3) and (4) into equation (2), we get By rearranging equation (5), we can write By representing the coefficient of , and as , and respectively then the equation (6) can be written as Where (here = = ) But there will be a little difference between the boundary conditions at first node (at point ) and for last node (at point ). In this case, linear approximation will be use between the adjacent nodal point and the boundary. For this reasons, we will introduce the linear approximation for temperatures between and as like: This can be rewrite as, Where Similarly, we will introduce the linear approximation for temperatures between and as like: This can be rewrite as, Where 0 + − − /2 ( /2 ) + When we will get the all discretized equations from equation (7), (9) and (11), to solve a problem we have to set up at each of the nodal points. The resulting system of linear algebraic equations will be solved to obtain the distribution of the property at each nodal point. There are different types of matrix solution technique (Gaussian Elimination Method, Tri-Diagonal Matrix Algorithm, Gauss-Jordan Elimination, etc.) exist which can be used to solve the system of linear algebraic equations. From this system, we will use Gaussian Elimination (GE) Method (Sasaki and Murao, 1982) and Tri-Diagonal Matrix Algorithm (TDMA) (El-Mikkawy, 2004) to find the values of temperature distributions .

RESULT AND DISCUSSION:
Consider a plane wall whose thickness is 1 (Fig  2) and this has a source effect in this heat conduction system. The two ends of this wall maintained at constant temperature as like 300° and 800° respectively. The thermal conductivity of this wall is equals to 0.5 / . , heat generation = 1000 / 3 , and the cross-sectional area is 10 2 . The main goal of this work is to calculate the temperature distributions of this wall by using our control volume method where we will consider the dimensions of the -and -directions are so large. For this reason the temperature gradients has significant change along the X -direction only.

Numerical Solution by FVM
First of all, we have to divide the whole length of the wall into some control volume (along onlydirections). In this work, we will calculate the temperature distribution in five internal points. That is, the difference between two control volumes is 0.2 (since total length of the wall is 1 , so = 1/5 =0.2 ). For our given problem (one-dimensional heat conduction) governing equation is, Where, is the thermal conductivity, is the temperature and is the source term. To get the discretization equation, we have to integrate equation (12) over the control volume. After simplifying (which has already described in Section 2) we get as like equation (6), By representing the coefficient of , and as , and respectively then the equation (13) can be written as Where (here = = ) In this problem, the given data are / .
By using these values, we get = = 25, = 50, and = 2000. Again by using the value of , , and into equation (14), we get the discretized equation for the nodal points 2, 3, and 4 as,  (14) can be written as, Which is the discretized equation for the nodal point 5.
Therefore, the set of discretized equations are:

Exact Solution
Our governing equation is ( ) + = 0 where boundary conditions, = 300° at = 0 , and = 800° , = 1 . By simplifying this we get our general solution, Where and are integral constant. By applying the boundary condition we get = 1500 and = 300, so the complete solution is, For different values of , we will get the temperature distribution at different nodal points ( Table 1) But if we divide the domain into more and more control volume, the error will be minimized than the before.

CONCLUSION:
In this work, to investigate one dimensional conductive heat transfer in a uniform plane wall, we have used the finite volume method (FVM). Obtaining discretized equation by FVM, we have applied Tri-Diagonal Matrix Algorithm (TDMA) and Gaussian elimination (GE) method. Finally, with the help of MATLAB R2014a, we represent a tabular and a graphical comparison between our numerical solution and exact solutions. The FVM has given outstanding results which have few errors with respect to exact solution but this error can be minimized by taking more and more grid points. That is, this method is very effective, accurate, reliable and easier to appliance in MATLAB or any other programming languages compared to the other exorbitant methods.

ACKNOWLEDGEMENT:
We are so much obligated to Md. Sirajul Islam for his patience and thoughtful guidance as a complete supervisor to continue this research work.

CONFLICTS OF INTEREST:
The authors declare that they have no competing interests with respect to the research.

Nomenclature
Control Volume boundary Area (