Scribd
Upload
25 views

2d Transient Flow Analysis

This document describes a Matlab program that uses a finite difference method to calculate 2D transient heat conduction through a rectangular solid. The program takes user inputs like material properties, boundary temperatures and convection coefficients, and initial temperature to calculate the temperature profile over time. It works by generating a matrix of temperature values at nodes in the solid. The program is limited to rectangular geometry and can have errors for large time steps or high convection coefficients due to memory limitations. Improvements could be made to the methods used to calculate time steps and grid spacing.

Uploaded by

Anil Kumar Rout
Copyright
© © All Rights Reserved
Available Formats
Download as PPT, PDF, TXT or read online on Scribd
25 views

2d Transient Flow Analysis

This document describes a Matlab program that uses a finite difference method to calculate 2D transient heat conduction through a rectangular solid. The program takes user inputs like material properties, boundary temperatures and convection coefficients, and initial temperature to calculate the temperature profile over time. It works by generating a matrix of temperature values at nodes in the solid. The program is limited to rectangular geometry and can have errors for large time steps or high convection coefficients due to memory limitations. Improvements could be made to the methods used to calculate time steps and grid spacing.

Uploaded by

Anil Kumar Rout
Copyright
© © All Rights Reserved
Available Formats
Download as PPT, PDF, TXT or read online on Scribd
You are on page 1/ 13

2D Transient

Conduction
Calculator
Using Matlab
Greg Teichert
Kyle Halgren
Assumptions
 Use Finite Difference Equations shown in
table 5.2
 2D transient conduction with heat transfer
in all directions (i.e. no internal corners as
shown in the second condition in table
5.2)
 Uniform temperature gradient in object
 Only rectangular geometry will be
analyzed
Program Inputs
 The calculator asks for
 Length of sides (a, b) (m)
 Outside Temperatures (T_inf 1-T_inf 4) (K)
 Temperature of object (T_0) (K)
 Thermal Convection Coefficient (h1-h4)
(W/m^2*K)
 Thermal Conduction Coefficient (k) (W/m*K)
 Density (ρ) (kg/m^3)
 Specific Heat (Cp) (J/kg*K)
 Desired Time Interval (t) (s)
Transient Conduction
 Example problem
 suppose we have an object with rectangular cross-
section with these boundary conditions:

Origin T_inf 2, h2

T_inf 1, h1 a T_inf 3, h3

b
T_inf 4, h4
Conditions
%Userdefined h values
h(1) = 10;
h(2) = .1;
h(3) = 10;
h(4) = .1;

%Boundary conditions
%Userdefined T infinity values in kelvin
T_inf(1) = 293;
T_inf(2) = 293;
T_inf(3) = 353;
T_inf(4) = 353;

%Initial condition (assume uniform initial temperature)


%Userdefined initial temperature value
T_0 = 573;

%Material properties
%Userdefined material values
k = .08;
rho = 7480;
c_p = .460;

%Userdefined physical variables


a = 1; %height of cross section
b = 1.3; %width of cross section
t = 3600; %time at which results are given
Time Step (Δt)
 We assumed a value of Δx = Δy = gcd(a,
b)
 Using each of the conditions (except the
second) in the table 5.2, we calculate the
Δt and choose the smallest value
 Using that Δt we calculate Fo
 Our outputs for delta_x, delta_t, Fo
respectively
 0.0500, 3.7078, 0.0345
Method
 Using the Finite Difference Method,
matlab generates a matrix of
temperature values that are represented
in the graph shown on the next slide
 This method allows for the calculation of
every node in any 2D direction
Results
Transient conduction (the origin of the plot is the top left corner of the cross section) Transient conduction (the origin of the plot is the top left corner of the cross section)

550

550
500

500
450

Temperature (K)
450 400
Temperature (K)

400 350

300
350

250
300 1.5
0 1

250 0.5 0.5 0.2 0


0.6 0.4
0 0.2 0.4 0.6 0 0.8
0.8 1 1.2 1 1
1.4 b (m)
a (m)
a (m)
b (m)
Solution to different Problem
%Userdefined h values
h(1) = 0;
h(2) = 1000;
h(3) = 1000;
h(4) = 100;
Transient conduction (the origin of the plot is the top left corner of the cross section)

%Boundary conditions
%Userdefined T infinity values in kelvin
T_inf(1) = 273; 600

T_inf(2) = 150; 500


T_inf(3) = 590;

Temperature (K)
T_inf(4) = 273; 400

300
%Initial condition (assume uniform
initial temperature) 200
%Userdefined initial temperature value
T_0 = 250; 100
1.5
1
%Material properties 1
0.8
0.6
%Userdefined material values 0.5 0.4
k = .8; 0 0
0.2
b (m) a (m)
rho = 1000;
c_p = .460;

%Userdefined physical variables


a = 1; %height of cross section
b = 1.3; %width of cross section
t = 20; %time at which results are given
Conclusion and
Recommendations
 Works only in rectangular geometry
 High values of h and t>1 causes errors to
occur due to lack of memory
 Use a better method to find Δx and Δt
Appendix-References
 Incropera, Frank P. DeWitt, DaviD P.
Fundamentals of Heat and Mass Transfer
Fifth Edition, R. R. Donnelley & Sons
Company. 2002 John Wiley & Sons, Inc
Appendix-hand work
Appendix-hand work

You might also like