Professional Documents
Culture Documents
Ss T Conductmnbvjkb
Ss T Conductmnbvjkb
340
NUMERICAL METHODS
341
CHAPTER 5
(a)
(b) (d)
(e) (f)
(c)
(g)
FIGURE 5–56
The SS-T-CONDUCT Input window
cengel_ch05.qxd 1/7/10 10:21 AM Page 342
342
NUMERICAL METHODS
(a)
(a)
FIGURE 5–57
The SS-T-CONDUCT Tabular and Graphical Output windows
cengel_ch05.qxd 1/5/10 3:05 PM Page 343
343
CHAPTER 5
Discussion The limitation on the time step of the explicit method can be
explored by simply changing the time step in item a of Fig. 5–56. The stabil-
ity criterion for this problem requires a maximum time step of 15.5 s. By
using time steps much larger than 15.5 s, the solution for the nodal temper-
atures will oscillate and at times become negative (violating the second law of
thermodynamics). One can also see that no matter what time step is used for
the implicit method, the solution remains stable.
h, T⬁
1 2 3
344
NUMERICAL METHODS
(the primary coefficient) to be greater than or equal to zero for all nodes. The
smallest primary coefficient in the nine equations here is the coefficient of
T1i in the T1i 1 expression, and thus the stability criterion for this problem can
be expressed as
hl 1 l2
1 4t 4t 0 → t → t
k 4(1 hl/k) 4(1 hl/k)
(0.1 m)2
t 179 s
4(12 10 6 m2/s)[1 (45 W/m2. K)(0.1 m)/(28 W/m.K)]
Therefore, any time step less than 179 s can be used to solve this problem. For
convenience, let us choose the time step to be t 60 s. Then the mesh
Fourier number becomes
6
t (12 10 m /s)(60 s)
2
t 2
0.072 (for t 60s)
l (0.1 m)2
Using the specified initial condition as the solution at time t 0 (for i 0),
sweeping through the 3 finite difference equations will give the solution at a
time step of 60 s.
The solution at the center node (node 5) can also be determined using the
SS-T-CONDUCT software. On the SS-T-CONDUCT Input window, select the
2-Dimensional Transient Problem tab (item (a) Fig. 5–59). Enter the problem
(c)
(e)
FIGURE 5–59
The SS-T-CONDUCT Input window for 2-Dimensional Transient Problem
cengel_ch05.qxd 1/7/10 4:46 PM Page 345
345
CHAPTER 5
(b)
(a)
FIGURE 5–60
The SS-T-CONDUCT Tabular and Graphical Output windows
cengel_ch05.qxd 1/7/10 10:27 AM Page 346
346
NUMERICAL METHODS
FIGURE 5–61
The SS-T-CONDUCT material properties list
parameters in the appropriate text boxes listed in the Problem Parameters box
(item (b) in Fig. 5–59).
To compute the nodal temperatures using the explicit method, the explicit
scheme radio button is selected (item (c) in Fig. 5–59). The boundary conditions
for this problem are specified in the Boundary Conditions box (item (d) in
Fig. 5–59). Once all the necessary inputs are entered, the results are computed
by clicking on the Calculate Temperature button (item (e) in Fig. 5–59).
The computed results are presented in the SS-T-CONDUCT Tabular Output
window (Fig. 5–60a.) The nodal temperatures are tabulated as indicated by
item (a) in Fig. 5–60a. The temperature at the center of the bar after 20 min
(or 1200 s) is 379.31°C. The solution at different times can be viewed by se-
lecting from the drop-down list (item (b) in Fig. 5–60a). In doing so, the tem-
perature at the center of the bar, after steady conditions are established, is
determined to be 1023°C, after about 6 hours.
To view the graphical output, select the Graphical Output tab and the tem-
perature contour for this problem is presented (Fig. 5–60b). User may view the
temperature contour at different time by controlling the slider bar.
Discussion The SS-T-CONDUCT software also has a built-in feature for user to
choose the properties (thermal conductivity and thermal diffusivity) of various
materials or explore the influence of properties of different materials on the
nodal temperatures. To select a material, click on the Properties menu on the
menu bar on the SS-T-CONDUCT Input window (Fig. 5–59) to see the available
list of materials (Fig. 5–61).