2D Conduction Calculator Using Matlab

2D Conduction Calculator Using Matlab

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 crosssection with these boundary conditions: Origin T_inf 2, h2 T_inf 1, h1 a b T_inf 4, h4 T_inf 3, h3

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 Temperature (K)

500 Temperature (K) 450 400 400 350 300 350 250 1.5 300 0 250 0

450 1 0.5 0.5 0.2 0.4 0.6 0.8 1 b (m) 1.2 1.4

1 a (m) b (m) 0 1 0.6 0.8 a (m) 0.4 0.2 0 Solution to different Problem

%Userdefined h values h(1) = 0; h(2) = 1000; h(3) = 1000; h(4) = 100; %Initial condition (assume uniform initial temperature) %Userdefined initial temperature value T_0 = 250; %Material properties %Userdefined material values k = .8; 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 600

500 Temperature (K) %Boundary conditions %Userdefined T infinity values in kelvin T_inf(1) = 273; T_inf(2) = 150; T_inf(3) = 590; T_inf(4) = 273; Transient conduction (the origin of the plot is the top left corner of the cross section) 400 300 200 100 1.5 1 0.5 b (m)

0 0 0.2 0.4 a (m) 0.6 0.8 1 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

Recently Viewed Presentations

  • Lesson 13 Figurative Lang. & Foils

    Lesson 13 Figurative Lang. & Foils

    Dr. Seus wrote in iambic tetrameter. If he had split his lines differently, he would have used enjambment:I do not eat green eggs. And ham, I will not eat them, Sam I am Internal Rhyme * A rhyme created by...
  • Data Analytics: The Data Mining Process

    Data Analytics: The Data Mining Process

    The Data Mining Process. Core Problems & Approaches . Problems: ... Establish Atlantic Canada as a leader. Education . Develop and deliver KDD seminar material at the academic and industry level . ... (ANNs, IDTs, Regression models)
  • The History of the Measurement of the Muon

    The History of the Measurement of the Muon

    Based on carefully measured muon rates at different altitudes, with and without absorbers, Rossi and Hall obtained a lifetime of 2.4 ± 0.3 µs. Direct measurement of the muon lifetime. Scenes from the 1962 movie "Time Dilation". Direct measurement of...
  • Cold War & The Postwar World

    Cold War & The Postwar World

    Cold War & The New Red Scare SS.A.3.4.9; SS.A3.4.10 Communism In Asia 1949: Mao Zedong wins Chinese Revolution China turns to Communism Chinese Nationalists flee to exile in Taiwan Japan now seen as best chance to defend against Communism in...
  • SERA PROPOSAL ON NATURAL DISASTERS CHRISTA D. COURT,

    SERA PROPOSAL ON NATURAL DISASTERS CHRISTA D. COURT,

    Community informed, science driven results to measure progress towards resiliency at multiple scales. Tool and information for disaster preparation efforts. Generalizable methodology that can be extended to other geographies, sectors, disasters (e.g. man-made disasters), etc.
  • Global review of gender statistics

    Global review of gender statistics

    Global review of gender statistics Key findings in UNECE Region Objectives of the survey Review gender statistics in national statistical systems Understand how gender statistics are mainstreamed in NSO, in all areas of statistics Identify good practices and challenges in...
  • Early Postoperative Ambulation Pricilla Puente University of South

    Early Postoperative Ambulation Pricilla Puente University of South

    Ineffective peripheral tissue perfusion r/t circulatory stasis, prolonged immobility aeb fatigue . Acute pain r/t inflammation in surgical area aeb patient rating pain a 7/10. Urinary retention r/t anesthesia, pain, unfamiliar surroundings aeb urine output of 20cc in 3 hours.
  • Chapter 4 -- Modular Combinational Logic

    Chapter 4 -- Modular Combinational Logic

    Chapter 4 -- Modular Combinational Logic Decoders Decoder Realization More complex decoders Example 4.1 -- Realize f(Q,X,P) = m(0,1,4,6,7) = M(2,3,5) Example 4.1 (concluded) SN 74138 Decoder Module Four-to-one multiplexer design Use a 74151A multiplexer to Realize f(x1,x2,x3) = m(0,2,3,5)...