Scribd
Upload
300 views9 pages

Fin Cooling Analysis with MATLAB

The document describes a complex engineering problem of determining the temperature distribution of fins used to cool a hot surface. The surface is initially at 90°C and the fins are made of copper, are 6 cm long and have a cross-sectional area of 5x10-6 m2. The problem is modeled using finite elements in MATLAB to calculate the temperature along the fin length. The temperature distribution is determined by varying the heat transfer coefficient from 10-50 W/m2°C and fin length from 6-12 cm in steps of 2 cm each time.

Uploaded by

Ali Raza
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
300 views9 pages

Fin Cooling Analysis with MATLAB

The document describes a complex engineering problem of determining the temperature distribution of fins used to cool a hot surface. The surface is initially at 90°C and the fins are made of copper, are 6 cm long and have a cross-sectional area of 5x10-6 m2. The problem is modeled using finite elements in MATLAB to calculate the temperature along the fin length. The temperature distribution is determined by varying the heat transfer coefficient from 10-50 W/m2°C and fin length from 6-12 cm in steps of 2 cm each time.

Uploaded by

Ali Raza
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

COMPLEX ENGINEERING PROBLEM

PROBLEM STATEMENT:

A hot surface is cooled by attaching fins (called pin fins) to it, as shown in Figure 1. The surface of the
plate (left end of the pin) is 90 °C. The fins are 6 cm long and have a cross-sectional area of 5x 10-6 m2
with a perimeter of 0.006 m.
The fins are made of copper (k = 400 W/m-°C). The temperature of the surrounding air is T infinite =
20 °C with heat transfer coefficient on the surface (including the end surface) of h = 10 W/ m2-°C. A
model of the typical fin is also shown in Figure 2.

Figure 1

Figure 2
a. Use three elements in your finite element model to determine the temperatures along the fin
length.
b. Use PDE MATLAB tool, determine the temperature distribution around the fins.

Show the contours of temperature distribution and defines the maximum temperature. c.
Compare the results with changing values of
• Heat transfer coefficient from 10 W/ m2-°C to 50 W/ m2-°C with a step of 10 W/ m2-°C (5
iterations)
• Length of fins from 6 cm to 12 cm with a step of 2 cm (5 iterations)
STEP 3: FORCE VECTOR CALCULATION:
As there is no internal heat generation so Q=0 and also no flux so q*=0. WE only have the
following considerations for force vector calculations:
CODE FILE GENERATED BY PDE TOOLBOX:
pdetool('appl_cb',9);
set(ax,'DataAspectRatio',[3.25 1 64]);
set(ax,'PlotBoxAspectRatio',[1 0.49129592909413505 0.49129592909413505]);
set(ax,'XLim',[-0.022739955357143044 0.078822544642856956]); set(ax,'YLim',[-
0.011288139763779043 0.019961860236220957]); set(ax,'XTickMode','auto');
set(ax,'YTickMode','auto');
% Geometry description: pderect([0 0.059999999999999998
0.0025000000000000001 0],'FIN');
set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String','FIN'
)
% Boundary
conditions:
pdetool('changemode',0)
pdesetbd(4,...
'dir',...
1,...
'1',... '90')
pdesetbd(3,...
'neu',...
1,...
'10',... '0')
pdesetbd(2,...
'neu',...
1,...
'10',... '0')
pdesetbd(1,...
'neu',...
1,...
'10',...
'0')
coefficients:
pdeseteq(1,...
'400',...
'10',...
'(0)+(10).*(20)',...
'(1.0).*(1.0)',...
'0:10',...
'0.0',...
'0.0',... '[0 100]')
setappdata(pde_fig,'currparam',...
['1.0';... '1.0';...
'400';...
'0 ';...
'10 ';...
'20 '])
% Solve parameters:
Graph 3: Temperature change with length of fin:

PART C

HEAT TRANSFER COEFFICIENT VARIATIONS:

For h=30 W/m2 0C

You might also like