Scribd
Upload
4 views3 pages

Chapter3_Prob40

The document outlines a method for solving a nonlinear system of equations related to heat transfer in a coating process using fixed-point iteration. It provides the iteration functions and a MATLAB program to carry out 100 iterations, starting with initial values, resulting in final temperatures and radiosities. The final computed values are T_c = 481 K, J_c = 6222 W/m2, T_h = 671 K, and J_h = 10504 W/m2.

Uploaded by

josa.palaganas.up
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
4 views3 pages

Chapter3_Prob40

The document outlines a method for solving a nonlinear system of equations related to heat transfer in a coating process using fixed-point iteration. It provides the iteration functions and a MATLAB program to carry out 100 iterations, starting with initial values, resulting in final temperatures and radiosities. The final computed values are T_c = 481 K, J_c = 6222 W/m2, T_h = 671 K, and J_h = 10504 W/m2.

Uploaded by

josa.palaganas.up
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
You are on page 1/ 3

1

3.40 A coating on the panel surface is cured by radiant energy


from a heater. The temperature of the coating is determined by
Ambient at
radiative and convective heat transfer processes. If the radiation is 298 K
treated as diffuse and gray, the following nonlinear system of
simultaneous equations determine the unknowns J h , T h , J c , T c : coating
–8 4 panel
5.67 × 10 T c + 17.41T c – J c = 5188.18
J c – 0.71J h + 7.46T c = 2352.71
–8 4
5.67 × 10 T h + 1.865T h – J h = 2250
J h – 0.71J c + 7.46T h = 11093
where J h and J c are the radiosities of the heater and coating surfaces, respectively, and T h and T c are the
respective temperatures.
(a) Show that the following iteration function can be used for solving the nonlinear system of equations
with the fixed-point iteration method:
J c – 17.41T c + 5188.18 1 ⁄ 4 2250 + J h – 1.865T h 1 ⁄ 4
T c = -------------------------------------------------------
–8
T h = ------------------------------------------------
–8
5.67 × 10 5.67 × 10
J c = 2352.71 + 0.71J h – 7.46T c J h = 11093 + 0.71J c – 7.46T h
(b) Solve the nonlinear system of equations with the fixed-point iteration method using the iteration func-
tions from part (a). Use the following initial values: T h = T c = 298 K, J c = 3000 W/m2, and
J h = 5000 W/m2. Carry out 100 iterations, and plot the respective values to observe their convergence.
The final answers should be: T c = 481 K, J c = 6222 W/m2, T h = 671 K, J h = 10504 W/m2.
Solution
(a) The iteration functions are derived by solving the first equation for T c , the second equation for J c ,
the third equation for T h , and the fourth equation for J h .

(b) The following MATLAB program, written in a script file, solves the nonlinear system of equations
with the fixed-point iteration method using the iteration functions from part (a). The following initial val-
ues are used: T h = T c = 298 K, J c = 3000 W/m2, and J h = 5000 W/m2. The program carries out 100
iterations, display the value of the solutions after the last iteration, and plots the respective values to
observe their convergence.

clear all
Tc(1)=298; Th(1)=298; Jc(1)=3000; Jh(1)=5000; IT(1)=1;

Excerpts from this work may be reproduced by instructors for distribution on a not-for-profit basis
for testing or instructional purposes only to students enrolled in courses for which the textbook
has been adopted. Any other reproduction or translation of this work beyond that permitted by
Sections 107 or 108 of the 1976 United States Copyright Act without the permission of the
copyright owner is unlawful.
2

for i=2:100
IT(i)=i;
Tc(i)=((Jc(i-1)-17.41*Tc(i-1)+5188.18)/5.67E-8)^(1/4);
Th(i)=((2250+Jh(i-1)-1.865*Th(i-1))/5.67E-8)^(1/4);
Jc(i)=2352.71+0.71*Jh(i-1)-7.46*Tc(i-1);
Jh(i)=11093+0.71*Jc(i-1)-7.46*Th(i-1);
end
i
TcLast=Tc(i)
ThLast=Th(i)
JcLast=Jc(i)
JhLast=Jh(i)
subplot(2,2,1)
plot(IT,Tc)
xlabel('Iteration'); ylabel('Tc')
subplot(2,2,2)
plot(IT,Th)
xlabel('Iteration'); ylabel('Th')
subplot(2,2,3)
plot(IT,Jc)
xlabel('Iteration'); ylabel('Jc')
subplot(2,2,4)
plot(IT,Jh)
xlabel('Iteration'); ylabel('Jh')

When the program is executed, the following is displayed in the Command Window.

i =
100
TcLast =
481.0273
ThLast =
671.1240
JcLast =
6.2222e+003

Excerpts from this work may be reproduced by instructors for distribution on a not-for-profit basis
for testing or instructional purposes only to students enrolled in courses for which the textbook
has been adopted. Any other reproduction or translation of this work beyond that permitted by
Sections 107 or 108 of the 1976 United States Copyright Act without the permission of the
copyright owner is unlawful.
3

JhLast =
1.0504e+004

The following figure is displayed in the Figure Window.

600 700

550
600
500

450 500

Th
Tc

400 400
350
300
300

250 200
0 20 40 60 80 100 0 20 40 60 80 100
Iteration Iteration

7000 12000

11000
6000
10000

9000
Jh
Jc

5000
8000

7000
4000
6000

3000 5000
0 20 40 60 80 100 0 20 40 60 80 100
Iteration Iteration

Excerpts from this work may be reproduced by instructors for distribution on a not-for-profit basis
for testing or instructional purposes only to students enrolled in courses for which the textbook
has been adopted. Any other reproduction or translation of this work beyond that permitted by
Sections 107 or 108 of the 1976 United States Copyright Act without the permission of the
copyright owner is unlawful.

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy