1. An iron rod of length L = 0.5 m, diameter D = 0.02 m, thermal conductivity k = 50W/(m 0C) protruding from wall is exposed to an ambient at T∞ = 30 0C. The heat transfer coefficient between the ambient and the rod surface is h= 100 W/(m2 0C) and the base of the rod is kept at a constant temperature T0 = 330 0C. Assuming a one-dimensional steady-state heat flow, the mathematical formulation of this problem is given by the following equations:
The finite difference approximation of this problem using 10 equal sub-regions (n = 10), is given by:
Where,
Now, find out the temperature differences at each locations (i.e. find out θ1, θ2, θ3, θ4, and θ5 ) solving the finite difference equations using Gauss Elimination method. Also calculate the temperatures at each locations (i.e. calculate T1, T2, T3, T4, T5...........) where Ti = θi + T∞ Plot the temperature (T) versus location (x).
2. Repeat the above analysis for 20 equal sub-regions. Plot T vs. x. Compare (1) and (2).
3. For n=20 case, extract the subdomain 4:10 rows and print the matrix for subdomain.