All instructions in file
Project 1 Project 1 Numerical Methods I: MIE334 H1S Department of Mechanical and Industrial Engineering University of Toronto Winter 2020 1 Project 1 MIE334: Numerical Methods I Due Date: Tuesday, March 3rd, 2020, 11:59 pm Problem Statement: In heat transfer, fins are used to increase the effective surface area of a heat exchanger and, thus, enhance the rate of convective heat transfer. As an example, fins are extensively used in electronic circuit boards to dissipate the generated heat by transistors to the environment. Figure 1. A finned heatsink used to cool an integrated circuit. Here, we are dealing with a single fin with triangular cross-section attached to a base object as shown in Figure 2. A constant heat flux of ?? ′′ = 3.2 ?/??2 enters the fin through its base. Also, the fin is being cooled with a stream of air having a temperature of ?∞ = 300 ?. The convective heat transfer coefficient between the fin and cooling air is assumed to be ℎ? = 18 ?/? 2? at the fin base and decrease linearly towards the tip of the fin where ℎ? = 6 ?/?2?. The fin is made of aluminum for which the properties―density, specific heat capacity, conductivity, and melting point―are listed in Table 1 along with the geometrical dimensions of the fin. The fin and the base are designed in such a way that ? is much longer than ? (i.e., ? ≫ ?). Figure 2. Schematics of the triangular fin. Table 1. Geometrical dimensions of the fin along with properties of aluminum. ? (??) ? (??) ? (??/?3) ?? (?/?? ?) ? (?/??) ??????? ????? (?) 1 7 2720 871 202 930 2 We are trying to obtain the temperature profile inside the fin as a function of ?. For this purpose, a control volume approach can be used which is explained in the following section. Solution Approach As mentioned above, the width ? of the fin is much bigger than its length ?. Therefore, temperature variations along the width of the fin can be neglected and ? can be set to unity (i.e., ? = 1). Also, since the fin is relatively thin (due to its thickness being smaller than its length) and the conductivity of aluminum is high, temperature gradients in the ? direction can also be neglected without losing much accuracy. Based on these assumptions, the problem can be solved one-dimensionally only in the ? direction. To further simplify the problem, radiation from the fin to the surrounding is neglected and the fin is considered to be in a steady state; meaning that the temperature profile does not change with time. Figure 3. Schematics of the solution domain divided into ? control volumes. The length of all the control volumes in the ?-direction is ∆? except the first and last ones for which the length is ∆? 2 . As seen in Figure 3, the fin is divided into ? control volumes. The first and the last control volumes include the boundary conditions. Due to the triangular shape of the fin, the cross- section of the fin changes along the x-direction which should be carefully accounted for in the solution procedure. ??,? ′′ ? ?− 1 2 ′′ ? ?+ 1 2 ′′ ??,? ′′ Figure 4. Geometrical parameters of a control volume (left) and heat fluxes entering/exiting the control volume (right). Note that ?? is NOT equal to ∆?. 3 For the control volume at node ?, the energy equation can be written by balancing the energies entering and exiting the boundaries of the control volume as follows: ∑(�̇�?? − �̇�???) = 0 → ? ?− 1 2 ′′ × ? ?− 1 2 − ? ?+ 1 2 ′′ × ? ?+ 1 2 − 2 × ??,? ′′ × ?? = 0 (1) where ? ?− 1 2 ′′ and ? ?+ 1 2 ′′ show the rate of heat conduction per unit area entering/exiting the left (? − 1 2 ) and right (? + 1 2 ) boundary surfaces of the control volume ?, respectively, and ??,? ′′ shows the rate of heat convection exiting from the upper and lower boundaries to the environment per unit area. All the fluxes entering the control volume are taken to be positive (adding to the energy content) while the exiting fluxes are considered to be negative (reducing the energy content). From heat transfer you know that the rate of heat conduction per unit area (flux) can be obtained using the Fourier’s law: ?′′ = −? ?? ?? (2) Whereas, the convective heat transfer per unit area can be obtained using the Newton’s cooling law: ?? ′′ = ℎ(? − ?∞) (3) Using these relation, Equation 1 can be rewritten as: −?? ?− 1 2 ( ?? ?? ) ?− 1 2 + ?? ?+ 1 2 ( ?? ?? ) ?+ 1 2 − 2 × ℎ(?) × ??(?? − ?∞) = 0 (4) In the above relation the two derivatives ?? ?? can be approximated using the Taylor’s series expansion―at the left and right boundaries of the control volume―as you have learnt in this course. For ( ?? ?? ) ?− 1 2 a 2nd order accurate central approximation can be developed as follows: ( ?? ?? ) ?− 1 2 = ?? − ??−1 ∆? (5) where ∆? is the distance between the centers of two adjacent control volumes. A similar relation can also be developed for ( ?? ?? ) ?+ 1 2 : ( ?? ?? ) ?+ 1 2 = ??+1 − ?? ∆? (6) Consequently, Equation 4 will be transformed into a system of equations which should be solved to obtain the values of temperature at every control volume. 4 Project Requirements: Online Report Submissions: You should upload your report in the designated place on Quercus. You have to submit a single *.pdf file containing all parts. 1. First, using the Taylor’s series expansion show why the numerical approximation developed for ( ?? ?? ) ?− 1 2 is central and 2nd order accurate. Hint: Notice that ( ?? ?? ) ?− 1 2 is the temperature derivative at point (? − 1 2 ) NOT at (?). Submission: Provide the complete details of your solution for this part along with some discussion to support your answer. 2. Substitute the derived relations for ?? ?? (Equations 5 and 6) back into Equation 4, factor out the terms containing ?, and rearrange the obtained relation to reach the following form: ????−1 + ???? + ????+1 = ?? (7) which clearly results in a tridiagonal system of equations. Submission: Show all the relations you obtain for ??, ??, ??, and ??. In each case include the complete details of the derivation procedure. 3. Equation 7 can only be used for the intermediate control volumes (i.e. for ? = 2, 3, … , ? − 1). For the first and last control volumes the situation is slightly different. Using the figure below, develop an expression for the first control volume in the same manner as in Equation 1 by starting with energy balancing. In this procedure, you should finally find relations for ?1, ?1, ?1, and ?1. ??,1 ′′ ?? ′′ ?1→2 ′′ ??,1 ′′ Now, follow a similar procedure to find the relations for ??, ??, ??, and ??. 5 Submission: Show the all relations you obtained for ?1, ?1, ?1, and ?1 as well as ??, ??, ??, and ??. In each case include the complete details of the derivation procedure. 4. In the final step you are required to develop a program in MATLAB. Using the relations you obtained for ??, ??, and ?? (and the boundary conditions), form the tridiagonal matrix of coefficients. Also form the matrix of constants using the formula obtained for ??. Subsequently, employ the Tridiagonal Matrix Algorithm (also known as the Thomas Algorithm) to solve the system of equations and obtain the values of ?? for all the control volumes. Submission: Plot the temperature variations against ? for different number of control volumes (? = 5, 10, 20, 40, and 80) all in the same figure. The figures should be complete in terms of axis labels, units, legends, title, line style, etc. Comment on your results. How does, in your opinion, increasing the number of control volumes affect your results? Also, report the values of temperature at the tip and base of the fin for different values of ?. Which one of these values you think represents the true value? Explain your answer. 5. Assume that the stream of cooling air is suddenly turned off. This causes the convective heat transfer coefficient to fall down to 2 ?/?2? on the whole surface of the fin. What would happen in this situation? Submission: Plot the temperature variation resulted in this case together with the one you obtained in Part 4 in the same figure. Use only ? = 40 for both cases. Write down your observations based on the results. Online Code Submissions: You should upload your m-files in the designated place on Quercus. You have to submit a single *.zip file containing all of your m-files. The name of the *.zip file should be: lastname_studentID_project1.zip For example: einstein_9999999999_project1.zip To test your program, the TAs will first declare ?, an array containing the values of q''?, ?∞, ??, ??, ? and another array containing ? and ?. Then they will call your program in the Command Window. For example: >> n = 20; >> properties = [32000; 300; 18; 6; 202]; >> dimensions = [0.07; 0.01]; >> einstein_9999999999_project1(properties, dimensions, n) You should define the name of your main function to make possible such an action in MATLAB. Once called, you program should give out an array containing the values of ?, an array containing the values of ?, and a figure where the temperature is plotted against ?. It is strongly recommended that you check your program in the same 6 manner to avoid any problems. It is your responsibility to make sure that your code works