Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

Size: px
Start display at page:

Download "Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer"

Transcription

1 Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here, Q/A = 35 W/m, T = 35 C, x = 0.08 m Thus, we get: k = 0.08 W/m K. The wattage of the heater s the same as the rate at whch energy s transferred from the nsde to the outsde, across the not-so-well nsulated wall. Ths rate of energy transfer s gven by: 3. We start off by drawng the system dagram, whch s as follows: We begn by consderng transfer of heat between the nterface of the nsulaton & ppe and the nsde wall of the ppe. For ths, the energy transferred s gven by: Here, Q/L = 0 W/m, k = W/m K, r = 0.04 m, r o = m, x = m Thus, we get: A lm = 0.59L m Substtutng these values, we get: T = 0.1 C Thus, Temperature at nsde surface of ppe = = 6.1 C We now consder transfer of heat between the nterface of the nsulaton & ppe and the outsde of the nsulaton. Usng the same expresson as earler for the energy transferred, we have: Q/L = 0 W/m, k = 0.0 W/m K, r = m, r o = m, x = 0.01 m Thus, we get: A lm = 0.97L m Substtutng these values, we get: T = 33.7 C Thus, Temperature at outsde of nsulaton = = 8.3 C We now consder transfer of heat between the nsde surface of the ppe and the bulk of the product. For ths, the energy transferred s gven by: Q = h A ( T) Here, Q/L = 0 W/m, r = 0.04 m, h = 00 W/m K Thus, we get: A = 0.51L m Substtutng these values, we get: T = 0.4 C

2 Thus, Bulk temperature of product = = 6.5 C We now consder transfer of heat between the outsde surface of the nsulaton and the bulk ar outsde. For ths, the energy transferred s gven by: Q = h A ( T) Here, Q/L = 0 W/m, r = m, h = 5 W/m K Thus, we get: A = 0.33L m Substtutng these values, we get: T = 1 C Thus, Bulk temperature of ambent ar = = 16.3 C 4. The system dagram s as shown below: We start off by consderng convecton between the wall of the heat exchanger and the bulk flud. The energy transferred s gven by: Q = h A ( T) wth A = r L Here, Q/L = 50 W/m, r = 0.03 m, T = = 5 C Thus, A = 0.189L m Hence, h = 53 W/m K We then consder transfer of heat across the metal ppe. The energy transferred s gven by: Q = k A lm ( T/ r) Here, Q/L = 50 W/m, A lm = 0.7L (based on r = 0.03 m, r 0 = 0.06 m), r = 0.03 m Thus, T = 0.6 C Hence, T = = 74.4 C 1 We then consder transfer of heat across the nsulaton. The energy transferred s gven by: Q = k A lm ( T/ r) Here, Q/L = 50 W/m, A lm = 0.5L (based on r = 0.06 m, r 0 = 1.1 m), r = 0.05 m Thus, T = 48.1 C Hence, T = = 6.3 C 5. Ths s a forced convecton problem wth the mean temperature of water beng: (0 + 60)/ = 40 C., L = 10 m Thus, the flow s lamnar. Now, (N Re)(N Pr)(d c/l) = 16.6 Usng equaton 4.64 on page 91 (3 ed: pg. 48), we get: N Nu = 5.3 = hd/k Thus, h = 67 W/m K

3 th 6. a. Lamnar flow: We need to use equaton 4.64 (4 ed: pg. 91; 3 ed: pg. 48) for ths stuaton. When the velocty of flow s doubled, Reynolds number s the only quantty that s affected n ether equaton. Let us examne equaton It can be seen from ths equaton that as velocty s doubled, 0.33 N Re s doubled, whch n turn results n a factor of (= 1.57) ncrease n N Nu. Thus, t can be estmated that doublng the velocty results n approxmately 5.7 % ncrease n the Nusselt number or the convectve heat transfer coeffcent under lamnar flow condtons. th b. Turbulent flow: We need to use equaton 4.67 (4 ed: pg. 91; 3 ed: pg. 49) for ths. 0.8 When the flow velocty s doubled, N Re s doubled and hence N Nu ncreases by (= 1.74). Thus, t can be estmated that doublng the velocty results n approxmately 74 % ncrease n the Nusselt number or the convectve heat transfer coeffcent under turbulent flow condtons. 7. The loss n energy as water flows through the ppe s gven by: Ths s the energy that s beng transferred across the ppe. Thus, Q = U A T = 1,594 W (1) lm lm and () From equatons (1), (), and (3), we get: (3) A = r L, A = r L r = 0.05 m, r = m, L = 100 m, A = 3.97 m x = m, h = 7 W/m K, k = 10 W/m K o o o lm (ppe) ppe o ppe Ths yelds: Thus, (4) where R o s the outsde radus ncludng the nsulaton. We now determne the nsde heat transfer coeffcent as follows: The bulk flud temperature s: ( )/ = 85 C

4 Now, We then use eqn 4.67 on pg. 91 (3 ed: pg. 49) to get: Thus, h = 358 W/m K. Usng ths n equaton (4), we get: Thus, the flow s turbulent. N = 53 = (h D)/k Nu Thus, any combnaton of R o and k ns that satsfes the above equaton wll be suffcent to mnmze the heat loss as per the requrements. For example, R = 0.06 m ( x = m), k = 0.1 W/m K o ns ns th 8. Ths s a case of flow over an mmersed object. We need to use eqn 4.69 (4 ed: pg. 9; 3 ed: pg. 49) The flm temperature n ths case s ( )/ = 50 C. The resultant velocty n ths case s 50 cm/s (= 0.5 m/s) N = 0.71 Thus, we get: N = 10 = hd/k = h (0.008)/0.07 Pr Nu h = 34 W/m K 9. The system dagram for the problem s as shown below: Consderng transfer of heat from the nsde surface of the ppe and the outsde surface of the nsulaton, we get: For ppe, x = 0.0 m, k = 15 W/m K For nsulaton, x = 0.04 m, k = 0.04 W/m K For ppe, r = 0.06 m, r = 0.04 m For nsulaton, r = 0.10 m, r = 0.06 m Ths yelds, (A ) = 0.31 L m and (A ) = 0.49 L m lm ppe o o lm ns

5 Thus, Q/L = 17.1 W/m a. Consderng heat transfer from the bulk of the flud to the nsde surface of the ppe, we get: Q = h (A ) (90-75) wth A = r L where r = 0.04 m Thus, 17.1 = h ( )(0.04)(15) Ths yelds h = 5 W/m K b. Consderng transfer of heat between the nsde and outsde of the ppe, we get: Q = ka( T/ x) = 15 (0.31 L) (75 - T)/(0.0) Ths yelds, T = 74.9 C c. Consderng transfer of heat between the outsde of the nsulaton and the surroundng ar, we get: Q = h (A ) (40 - T ) wth A = r L where r = 0.10 m o o a o o o Thus, 17.1 = 10 ( )(0.10)(40 - T ) a Ths yelds T = 37.3 C 10. Ths s a free convecton problem wth D = 0.3 m and T flm = ( )/ = 50 C a N Pr = th N Ra = (N Gr)(N Pr) = x 10 Thus, from Table 4. (4 ed: pg ; 3 ed: pg ), a = 0.59, m = 0.5 m Thus, N = a (N x N ) = 64 = hd/k = h(0.3)/0.07 Thus, h = 6 W/m K Nu Gr Pr 11. Ths s a forced convecton problem and hence we begn by computng the Reynolds number Snce Reynolds number s greater than 10,000, the flow s turbulent N = 8.1 (at 15 C) Pr Thus, -6 wth = x 10 Pa s w Thus, N Nu = 90 = hd/k = h (0.05)/(0.587) Thus, h = 3405 W/m K 1. The system dagram s as shown below: Ths s a forced convecton problem. So we begn by determnng the Reynolds number. Here, -6 = 1 kg/s, D = m, = x 10 Pa s (at 40 C) We determned the vscosty of water at the bulk temperature [T bulk = ( )/ = 40 C]

6 Thus, N Re = 76,179 Snce the flow s turbulent, we use the followng equaton to determne N Nu: N Nu = 0.03 (N Re) (N Pr) ( b/ w) -6-6 Here, N Re = 76179, N Pr = 4.3 (at 40 C), b = x 10 Pa s, w = 9.38 x 10 Pa s Thus, we get: N Nu = 335 = hd/kf wth D = m, k f = W/m K Solvng, we get: h = 8349 W/m K 13. The system dagram s as follows: a. Ths s a forced convecton problem. Thus, we begn by determnng the Reynolds number. Here, -6 = 0.1 kg/s, D = 0.0 m, = 9.38 x 10 Pa s (at 95 C) Thus, N Re = 1,784 Snce the flow s turbulent, we use the followng equaton to determne N Nu: N Nu = 0.03 (N Re) (N Pr) ( b/ w) -6-6 Here, N Re = 1,784, N Pr = 1.84 (at 95 C), b = 9.38 x 10 Pa s, w = x 10 Pa s Thus, we get: N Nu = 8.5 = hd/kf wth D = 0.0 m, k f = 0.68 W/m K Solvng, we get: h = 805 W/m K b. Q = ha (95-90) = hoa o(87-13) wth A = rl, A o = rol, r = 0.01 m, r o = m, L = 0.9 m Solvng, we get: h = 137 W/m K o 14. In ths problem, all we need to do s determne the logarthmc mean temperature dfference. We note that T lm = ( T 1 - T )/[ln( T 1/ T)] Snce the steam s saturated steam, the temperature of the steam s the same at the nlet and outlet. Based on the pressure of kpa, we determne the steam temperature to be 10 C (from steam tables) Thus, T 1 = 10-0 = 100 C, T = = 50 C Thus, T = 7.1 C lm

7 15. Crtcal radus (r c) = k/h = 0.1/5 = 0.0 m Snce the radus of the tube s greater than the crtcal radus, addng any amount of nsulaton wll decrease the heat loss from the ppe U = Q/(A T lm) = [15.34 x 10 ] / [{( ) (0.05)(10)}{50.8}] = 157 W/m K 17. Consderng heat transfer from the hot water to the ambent ar, Q = 500 W, T = 85 C, h o = 0 W/m K A = rl = 1.6 m, A o = rol =.0 m ( x/ka) ppe = T/Q = 5/500 = 0.0 K/W Solvng, we get h = 49 W/m K 18. Snce we have a fnte cylnder, we consder to be a combnaton of an nfnte cylnder and an nfnte slab. For the nfnte cylnder, N B = hd/k = (40)(0.05)/50 Thus, N B = 0.04 For the nfnte slab, N B = hd/k = (40)(0.1)/50 Thus, N B = 0.08 Snce N < 0.1, the lumped parameter analyss can be used. B Here, T = 35 C, T = 0 C, T = 5 C 3 h = 40 W/m K, = 8314 kg/m, c p = 4190 J/kg K r = 0.05 m,l = 0. m 3 A = r L = 0.06 m V = r L = m Substtutng these values, we get: t = 16,10 s (= 4.5 hrs) a 19. Snce one of the dmensons s very large n comparson wth the other two dmensons, the product wll be assumed to nfnte n one drecton. The resultng product can be obtaned by the ntersecton of a product 4 cm thck n one dmenson (and nfnte n the other two drectons) and another product 10 cm thck n one dmenson (and nfnte n the other two dmensons). Thus, the temperature rato (TR) of the fnte object can be obtaned as the product of the TR s for the two nfnte objects. -7 = k/( c ) = 1.3 x 10 m /s,t = 3600 s p Frst slab: D = 0.0 m Thus, N B = hd/k = 4 (thus, k/hd = 0.5) and N Fo = 1. TR = 0.19 (From Hesler chart -- pg. 90)

8 Second slab: D = 0.05 m Thus, N B = hd/k = 10 (thus, k/hd = 0.1) and N Fo = 0.19 TR = 0.85 (From Hesler chart -- pg. 90) Thus, TR for fnte object = 0.19 x 0.85 = 0.16 = (T a - T) / (T a - T ) Thus, T = 19 C 0. a. The Bot number n ths case s: N B = hd/k = (75)(0.04)/0.4 = 7.5 k/hd = Thus, we use the Hesler chart. = k/( c p) = 1.35 x 10 m /s TR = 0.17 N = = t/d t = 4444 s (= 74 mn) Fo b. We need to use eqn on pg. 34 (3 ed: eqn 4.17 on pg. 84). 3 3 A = 4 r = 0.0 m, V = (4/3) r = m Substtutng these nto eqn (or eqn 4.17 n 3 ed), we get: Thus, t = 937 s (= 15.6 mn) c. The reason for the bg dfference between the answers n parts a and b are because the lumped parameter analyss s not vald. 1. a. A fnte cylnder s obtaned by the ntersecton of an nfnte cylnder and an nfnte slab. For nfnte cylnder: D = 0.03 m N B = hd/k = 1.5 Thus, k/hd = 0.67 For nfnte slab: D = 0.04 m N = hd/k =.0 Thus, k/hd = Thus, we use the Hesler chart. = k/( c ) = 1.17 x 10 m /s N = t/d = 0.47 for the nfnte cylnder and N = t/d = 0.6 for the nfnte slab Fo p Fo Thus, TR = 0.45 for the nfnte cylnder and TR = 0.75 for nfnte slab Thus, TR for fnte cylnder = (0.45) (0.75) = 0.34 = (T - T) / (T - T ) B a a T = 45.4 C b. Here, the TR for fnte cylnder s gven: (T a - T) / (T a - T ) = We know k/hd for the nfnte cylnder; but we do not know ts TR or N Fo Let us guess a value of t. Ths wll enable us to determne the TR for the nfnte cylnder. We then use ths t for the nfnte slab and thus determne the TR for the nfnte slab. We then multply the TR for the two nfnte objects and compare wth If t s not equal to 0.083, we attempt another guess value for t and contnue the ITERATIVE procedure tll the product of the two TR s s The Bot number n ths case s: N B = hd/k = (1000)(0.015)/0.4 = 37.5 k/hd = 0.03 We thus use the Hesler chart. A fnte cube of sde 0.03 m s obtaned by the ntersecton of three nfnte slabs, each havng a thckness of 0.03 m 1/3 TR for fnte cube = 0. Thus, TR for nfnte slab = (0.) = = k/( c p) = 1.17 x 10 m /s For one such nfnte slab: TR = 0.6 Thus, we get: N = 0. Thus, t = 385 s 3. a. The propertes requred for solvng the problem are those of the product b. h s the convectve heat transfer coeffcent for transfer of heat between the surface of the product and water c. Snce ths s an unsteady state problem, we begn by determnng the Bot number N = hd/k = (650)(0.05)/(0.5) = 3.5 B Fo

9 Snce the Bot number les between 0.1 and 40, we use the Hesler chart The temperature rato for the fnte cube s (85-50)/(85-5) = 0.44 We now note that a fnte cube can be obtaned as the ntersecton of 3 nfnte slabs, each of thckness equal to that of the length of the cube. 3 1/3 Thus, (TR nfnte slab) = TRfnte cube Thus, TR nfnte slab = (0.44) = 0.76 We thus use the Hesler chart for an nfnte slab wth TR = 0.76 and k/hd = 0.03 to obtan a Fourer number of Thus, ( t/d ) = 0. = k/( c p) = 0.5/(950)(3700) = 1.4 x 10 m /s D = 0.05 Thus, we get t = 880 s = 14.7 mn 4. Ths s an unsteady state heat transfer problem and hence we begn by computng Bot number. Snce ths s a fnte cylnder, we consder t to be a combnaton of an nfnte cylnder and an nfnte slab. For nfnte cylnder, N B = hd/k = (5)(0.04)/(10) = 0.0 For nfnte slab, N B = hd/k = (5)(0.05)/(10) = 0.05 th Snce Bot number s less than 0.1, we use the lumped parameter analyss (4 ed: eqn on pg. 34; 3 ed: eqn on pg. 84). Thus, A = rl + r = m 3 V = r L = m t = 1800 s Solvng, we get: (T - 5)/(90-5) = 0.85 Thus, T = 77.3 C 5. A fnte brck-shaped object can be obtaned by the ntersecton of 3 nfnte slabs. Neglect the convectve resstance mples that k/hd = 0 (h 0; n fact h s very hgh) d = 0.03 m, d = m, d = 0.04 m 1 3 Thus, (T - 5) / (-5-5) = T = 3.6 C 6. N = (0 x 0.0) / (5) = Snce N < 0.1, use the lumped parameter analyss. B B Thus, t = 7 s = 1 mn 7. a. A fnte cylnder s obtaned by the ntersecton of an nfnte cylnder and an nfnte slab. For nfnte cylnder: D = 0.04 m N B = hd/k =.0 Thus, k/hd = 0.5 For nfnte slab: D = 0.05 m N B = hd/k =.5 Thus, k/hd = Thus, we use the Hesler chart. = k/( c p) = 1.3 x 10 m /s N Fo = t/d = 0.5 for the nfnte cylnder and N Fo = t/d = 0.33 for the nfnte slab Thus, TR = 0.4 for the nfnte cylnder and TR = 0.9 for nfnte slab Thus, TR for fnte cylnder = (0.4) (0.9) = 0.36 = (T - T ) / (T - T ) T = 45. C b. Here, the TR for fnte cylnder s gven: (T - T ) / (T - T ) = (40-0)/(90-0) = 0.9 We know k/hd for the nfnte cylnder; but we do not know ts TR or N. Let us guess a Fo

10 value of t. Ths wll enable us to determne the TR for the nfnte cylnder. We then use ths t for the nfnte slab and determne the TR for the nfnte slab. We then multply the TR for the two nfnte objects. If ths s not equal to 0.9, we attempt another guess value for t and contnue the ITERATIVE procedure tll the product of the two TR s s The system dagram s as follows: Snce ths s a forced convecton problem, we begn by determnng the Reynolds number. Here, = 0.5 kg/s, D = 0.06 m, = Pa s (at 85 C) Thus, N Re = 10,610 Snce the flow s turbulent, we use the followng equaton to determne N Nu: N Nu = 0.03 (N Re) (N Pr) ( b/ w) Here, N Re = 10,610, N Pr = cp /k = 8 Snce we do not know the wall temperature, we assume that the vscosty correcton factor ( b/ w) s neglgble. Thus, we get: N Nu = 76 = hd/kf wth D = 0.06 m, k f = 0.5 W/m K Solvng, we get: h = 633 W/m K The heat loss n the un-nsulated case s gven by: = 0.5 (4000) (10) = 0,000 Thus heat loss n the nsulated case = 10,000 W (Ths means that the ext temperature of the water wll be 85 C) For the nsulated case, Q = = UA lm Tlm Thus, 1/UA lm = T lm/q = (1/hA ) + ( r/ka lm) + (1/hoA o) Here, T lm = 77.5 C (based on T 1 = = 80 C and T = = 75 C) A = rl = 3.77 m (based on r = 0.03 m) A o = rol = 4.4 m (based on r o = m) A lm = 4.08 m (based on A and A o) Also, h = 633 W/m K, h o = 100 W/m K Thus, 77.5/10000 = [1/(633)(3.77)] + [0.005/(k)(4.08)] + [1/(100)(4.4)] Smplfyng, we get: 198 = (4.08 k)/0.005 Solvng, we get: k = 0.4 W/m K 9. The Bot number n ths stuaton s gven by: N B = hd c/k Here, h = 10 W/m K, d = 0.05 m, k = 3.5 W/m K c

11 Thus, N B = 0.07 Snce N < 0.1, the lumped parameter analyss can be used. B Here, T a = 150 C, T = 5 C, T = 80 C, h = 10 W/m K, A = 6 (0.05) = m = 950 kg/m, c = 3700 J/kg K, V = (0.05) = 15 x 10 m Thus, TR = 0.483, ha/ cpv = s -1 Solvng, we get: t =,134 s p 30. a. The system dagram s as follows: Snce ths s a forced convecton problem, we begn by determnng the Reynolds number. Here, = 0. kg/s, D = 0.04 m, = Pa s Thus, N Re = 173 Snce the flow s lamnar, we use the followng equaton to determne N Nu: N Nu = 1.86 [N Re N Pr (d c/l)] ( b/ w) Here, N Re = 173, N Pr = cp /k = 41 Snce we do not know the nsde wall temperature, we assume that the vscosty correcton factor ( b/ w) s neglgble. Thus, we get: N Nu = 10.8 = hd/kf wth D = 0.04 m, k f = 0.5 W/m K Solvng, we get: h = 135 W/m K Consderng heat transfer between the outsde wall and the bulk of the flud, we get: Here, T = = 5 C h = 135 W/m K, k = 8 W/m K A = rl, A o = rol, r = 0.0 m, r o = 0.05 m, L = 10 m, x = = m Thus, A = 1.6 m, A o = 1.57 m, A lm = (A o - A )/[ln(a o/a )] = 1.41 m Solvng, we get: Q = 791 W

12 Now, Q = hoa o (T surface - T a) Here, Q = 791 W, A o = 1.57 m, T surface = 80 C, T a = 10 C Solvng, we get: h = 7. W/m K o Note: You can also solve ths problem consderng free convecton at the surface b. Here, r c = k/h o = 0.1/7. = m In ths case, r = 0.05 m; Snce r > r, addng any amount of nsulaton wll dec. heat loss o o c 31. a. The overall heat transfer coeffcent can be determned usng the followng equaton: Q = UA T (1) lm In ths case, lm = 0. (400)(70-5) = 54,600 W Snce there are no heat losses, ths amount of energy s the energy lost by hot water. Thus, 54,600 = = 0.4 (400) (90 - T ) Thus, T h(o) = 57.5 C Now, T lm = ( T 1 - T )/[ln( T 1/ T)] Here, T 1 = = 0 C, T = = 5.5 C Thus, T lm = 33.7 C Now, A lm = [ L (r o - r )]/ln (r o/r )] Here, L = 10 m, r o = 0.04 m, r = 0.05 m Substtutng, we get: A = 1.46 m lm Substtutng these values n equaton 1, we get: U = 1,110 W/m K b. When a rapd ntal rate of coolng s desred, a co-current confguraton s preferable snce the drvng force (temperature dfference) s hgher for ths confguraton. 3. When the lumped parameter analyss s vald, we can use the followng equaton to determne the temperature as a functon of tme: h(o) If the ntal and ambent condtons are the same n both cases, the temperature after a gven amount of tme would depend on the value of A/V. 3 For the cube, A = 6D, V = D ; thus, A/V = 6/D 3 For the sphere, A = 4 R, V = (4/3) R ; thus, A/V = 3/R = 6/D Thus, t would take the same amount of tme for the cube and the sphere to heat up (or cool down) to the same temperature.

13 33. The system dagram s as follows: In oer to determne T, we consder heat transfer between the nterface of the two layers of nsulaton and the outsde ar and wrte down the expresson for heat transfer between these two ponts as follows: Solvng, we get T = 8.7 C 34. Ths s a free convecton problem. Hence we begn by determnng N Gr. Assumng that the entre object s hot and that t s placed wth the 0 cm face on the counter-top, d = 0. m and T = (40 + 4)/ = C c flm N Pr = 0.71 N Ra = (N Gr)(N Pr) = 3.1 x 10 7 th Thus, from Table 4. (page 55 of 3 ed; pg. 99 of 4 ed), a = 0.59, m = 0.5 m Thus, N Nu = a (N Gr x N Pr) = 44.1 = hd c/k = h(0.)/0.05 Thus, h = 5.6 W/m K 35. Ths s a forced convecton problem. So, we begn by computng the Reynolds number. = (4)(0.35)/[( )(0.05)(0.0015)] = 11,883.6 Snce N Re > 10,000, the flow s turbulent and we use the followng equaton to determne h: N Nu = 0.03 (N Re) (N Pr) ( b/ w) Here, N Re = 11,883.6, N Pr = cp /k = (4100)(0.0015)/(0.55) = 11., b = w = Pa s Thus, we get: N Nu = 9.9 = hd c/kf wth d c = 0.05 m, k f = 0.55 W/m K Solvng, we get: h =,044 W/m K

14 36. A fnte cylnder s obtaned by the ntersecton of an nfnte cylnder and an nfnte slab. For nfnte cylnder: d c = 0.04 m, k = 0.54 W/m K N B = hd c/k = 1.85 Thus, k/hd c = 0.54 For nfnte slab: d c = 0.05 m, k = 0.54 W/m K N B = hd c/k =.31 Thus, k/hd c = Thus, we use the Hesler chart. = k/( c p) = 1.54 x 10 m /s t = 90 mn = 5400 s N Fo = t/d c = 0.5 for the nfnte cylnder and N Fo = t/d c = 0.33 for the nfnte slab Thus, TR = 0.4 for the nfnte cylnder and TR = 0.9 for nfnte slab Thus, TR for fnte cylnder = (0.4) (0.9) = 0.36 = (T - T ) / (T - T ) T = 45. C 37. a. Snce the tubular heat exchanger s well-nsulated, all the heat lost by hot water s ganed by the product. Thus, Solvng, we get: The energy transferred from hot water to the product s gven by: Q = 0. (400) (90-70) = W Also, Q = U A lm T lm (1) T lm = ( T 1 - T )/[ln( T 1/ T )] () Here, T 1 = = 60 C, T = = 30 C Substtutng these n equaton, we get: T lm = 43.3 C Also, A lm = (A o - A )/[ln(a o/a )] (3) Here, A = rl, A o = rol wth r = m, r o = m, L = 50 m Thus, A = 4.4 m, A o = 6.6 m Substtutng these n equaton 3, we get: A lm = 5.4 m Substtutng ths n equaton (1), we get: U = 7 W/m K b. Ths s a forced convecton problem. So, we begn be determnng the Reynolds number., L = 0 m Thus, the flow s lamnar. So, we use equaton 4.64 wth the vscosty correcton factor assumed to be 1.0. Thus, equaton 4.64 reduced to: N Nu = 1.86 {(N Re) (N Pr) (d c/l)} 0.33 Here, N Pr = (c p)( )/k = (4000) (0.005)/0.4 = 50 Also, d c = 0.06 m, L = 0 m Thus, N Nu = 6.9 Also, N Nu = hd c/kf Substtutng N = 6.9, d = 0.06 m, k = 0.4 W/m K yelds: h = 106 W/m K Nu c f 38. Ths s an unsteady state heat transfer problem. So, we have to determne the Bot number. However, we do not know the convectve heat transfer coeffcent. We determne that by performng a free convecton analyss.

15 d c = ½ D wth D = 0.01 m Thus, we get: d c = m We begn by determnng the N as follows: Gr We determned prop. of water at the flm temperature (= avg. of 15 C and 45 C = 30 C) Also, N Pr = 5.4 N Ra = (N Gr)(N Pr) =.91 x 10 6 Substtutng ths n the Nusselt number equaton n Table 4. (pg. 99), we get: N Nu = 4 = hd c/kf wth d c = m and k f = W/m K Thus, h = 940 W/m K We now begn the unsteady state heat transfer analyss as follows: The Bot number n ths case s: N B = hd c/k s = (940)(0.005)/0.4 = Snce N B > 0.1, we use the Hesler chart for a sphere. 1/N B = 1/11.75 = 0.09 TR = (T a - T)/ (T a - T ) = (45-30)/(45-15) = 0.5 From the Hesler chart, we get: N Fo = 0.15 N Fo = t/d c = Here, = k/( c p) = 1.5 x 10 m /s and d c = m Solvng, we get: t = 30 s

1. Nusselt number and Biot number are computed in a similar manner (=hd/k). What are the differences between them? When and why are each of them used?

1. Nusselt number and Biot number are computed in a similar manner (=hd/k). What are the differences between them? When and why are each of them used? 1. Nusselt number and Biot number are computed in a similar manner (=hd/k). What are the differences between them? When and why are each of them used?. During unsteady state heat transfer, can the temperature

More information

Principles of Food and Bioprocess Engineering (FS 231) Exam 2 Part A -- Closed Book (50 points)

Principles of Food and Bioprocess Engineering (FS 231) Exam 2 Part A -- Closed Book (50 points) Principles of Food and Bioprocess Engineering (FS 231) Exam 2 Part A -- Closed Book (50 points) 1. Are the following statements true or false? (20 points) a. Thermal conductivity of a substance is a measure

More information

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850) hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In

More information

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00 ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental

More information

MODULE 2: Worked-out Problems

MODULE 2: Worked-out Problems MODUE : Worked-out Problems Problem : he steady-state temperature dstrbuton n a one dmensonal wall of thermal conductvty 5W/m and thckness 5 mm s observed to be ( C) abx, where a C, B- c/ m, and x n meters

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

( ) 1/ 2. ( P SO2 )( P O2 ) 1/ 2.

( ) 1/ 2. ( P SO2 )( P O2 ) 1/ 2. Chemstry 360 Dr. Jean M. Standard Problem Set 9 Solutons. The followng chemcal reacton converts sulfur doxde to sulfur troxde. SO ( g) + O ( g) SO 3 ( l). (a.) Wrte the expresson for K eq for ths reacton.

More information

EXAMPLES of THEORETICAL PROBLEMS in the COURSE MMV031 HEAT TRANSFER, version 2017

EXAMPLES of THEORETICAL PROBLEMS in the COURSE MMV031 HEAT TRANSFER, version 2017 EXAMPLES of THEORETICAL PROBLEMS n the COURSE MMV03 HEAT TRANSFER, verson 207 a) What s eant by sotropc ateral? b) What s eant by hoogeneous ateral? 2 Defne the theral dffusvty and gve the unts for the

More information

Principles of Food and Bioprocess Engineering (FS 231) Problems on Heat Transfer

Principles of Food and Bioprocess Engineering (FS 231) Problems on Heat Transfer Principles of Food and Bioprocess Engineering (FS 1) Problems on Heat Transfer 1. What is the thermal conductivity of a material 8 cm thick if the temperature at one end of the product is 0 C and the temperature

More information

Numerical Transient Heat Conduction Experiment

Numerical Transient Heat Conduction Experiment Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use

More information

Outline. Unit Eight Calculations with Entropy. The Second Law. Second Law Notes. Uses of Entropy. Entropy is a Property.

Outline. Unit Eight Calculations with Entropy. The Second Law. Second Law Notes. Uses of Entropy. Entropy is a Property. Unt Eght Calculatons wth Entropy Mechancal Engneerng 370 Thermodynamcs Larry Caretto October 6, 010 Outlne Quz Seven Solutons Second law revew Goals for unt eght Usng entropy to calculate the maxmum work

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Chemical Engineering Department University of Washington

Chemical Engineering Department University of Washington Chemcal Engneerng Department Unversty of Washngton ChemE 60 - Exam I July 4, 003 - Mass Flow Rate of Steam Through a Turbne (5 onts) Steam enters a turbne at 70 o C and.8 Ma and leaves at 00 ka wth a qualty

More information

Digital Signal Processing

Digital Signal Processing Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over

More information

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018 MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

Physics 114 Exam 2 Fall 2014 Solutions. Name:

Physics 114 Exam 2 Fall 2014 Solutions. Name: Physcs 114 Exam Fall 014 Name: For gradng purposes (do not wrte here): Queston 1. 1... 3. 3. Problem Answer each of the followng questons. Ponts for each queston are ndcated n red. Unless otherwse ndcated,

More information

The Finite Element Method

The Finite Element Method The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of

More information

CIRCLE YOUR DIVISION: Div. 1 (9:30 am) Div. 2 (11:30 am) Div. 3 (2:30 pm) Prof. Ruan Prof. Naik Mr. Singh

CIRCLE YOUR DIVISION: Div. 1 (9:30 am) Div. 2 (11:30 am) Div. 3 (2:30 pm) Prof. Ruan Prof. Naik Mr. Singh Frst CIRCLE YOUR DIVISION: Dv. 1 (9:30 am) Dv. (11:30 am) Dv. 3 (:30 m) Prf. Ruan Prf. Na Mr. Sngh Schl f Mechancal Engneerng Purdue Unversty ME315 Heat and Mass ransfer Eam #3 Wednesday Nvember 17 010

More information

between standard Gibbs free energies of formation for products and reactants, ΔG! R = ν i ΔG f,i, we

between standard Gibbs free energies of formation for products and reactants, ΔG! R = ν i ΔG f,i, we hermodynamcs, Statstcal hermodynamcs, and Knetcs 4 th Edton,. Engel & P. ed Ch. 6 Part Answers to Selected Problems Q6.. Q6.4. If ξ =0. mole at equlbrum, the reacton s not ery far along. hus, there would

More information

PART I: MULTIPLE CHOICE (32 questions, each multiple choice question has a 2-point value, 64 points total).

PART I: MULTIPLE CHOICE (32 questions, each multiple choice question has a 2-point value, 64 points total). CHEMISTRY 123-07 Mdterm #2 answer key November 04, 2010 Statstcs: Average: 68 p (68%); Hghest: 91 p (91%); Lowest: 37 p (37%) Number of students performng at or above average: 58 (53%) Number of students

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Lecture 5.8 Flux Vector Splitting

Lecture 5.8 Flux Vector Splitting Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form

More information

The Analysis of Convection Experiment

The Analysis of Convection Experiment Internatonal Conference on Appled Scence and Engneerng Innovaton (ASEI 5) The Analyss of Convecton Experment Zlong Zhang School of North Chna Electrc Power Unversty, Baodng 7, Chna 469567@qq.com Keywords:

More information

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton

More information

Convection Heat Transfer. Textbook: Convection Heat Transfer. Reference: Convective Heat and Mass Transfer. Convection Heat Transfer

Convection Heat Transfer. Textbook: Convection Heat Transfer. Reference: Convective Heat and Mass Transfer. Convection Heat Transfer Convecton Heat Transfer Tetbook: Convecton Heat Transfer Adran Bean, John Wley & Sons Reference: Convectve Heat and Mass Transfer Kays, Crawford, and Wegand, McGraw-Hll Convecton Heat Transfer Vedat S.

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Flow Induced Vibration

Flow Induced Vibration Flow Induced Vbraton Project Progress Report Date: 16 th November, 2005 Submtted by Subhrajt Bhattacharya Roll no.: 02ME101 Done under the gudance of Prof. Anrvan Dasgupta Department of Mechancal Engneerng,

More information

HEAT TRANSFER THROUGH ANNULAR COMPOSITE FINS

HEAT TRANSFER THROUGH ANNULAR COMPOSITE FINS Journal of Mechancal Engneerng and Technology (JMET) Volume 4, Issue 1, Jan-June 2016, pp. 01-10, Artcle ID: JMET_04_01_001 Avalable onlne at http://www.aeme.com/jmet/ssues.asp?jtype=jmet&vtype=4&itype=1

More information

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems:

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems: TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases

More information

Estimation of the composition of the liquid and vapor streams exiting a flash unit with a supercritical component

Estimation of the composition of the liquid and vapor streams exiting a flash unit with a supercritical component Department of Energ oltecnco d Mlano Va Lambruschn - 05 MILANO Eercses of Fundamentals of Chemcal rocesses rof. Ganpero Gropp Eercse 8 Estmaton of the composton of the lqud and vapor streams etng a unt

More information

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient 58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

MTH 263 Practice Test #1 Spring 1999

MTH 263 Practice Test #1 Spring 1999 Pat Ross MTH 6 Practce Test # Sprng 999 Name. Fnd the area of the regon bounded by the graph r =acos (θ). Observe: Ths s a crcle of radus a, for r =acos (θ) r =a ³ x r r =ax x + y =ax x ax + y =0 x ax

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors Multple Lnear and Polynomal Regresson wth Statstcal Analyss Gven a set of data of measured (or observed) values of a dependent varable: y versus n ndependent varables x 1, x, x n, multple lnear regresson

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

y i x P vap 10 A T SOLUTION TO HOMEWORK #7 #Problem

y i x P vap 10 A T SOLUTION TO HOMEWORK #7 #Problem SOLUTION TO HOMEWORK #7 #roblem 1 10.1-1 a. In order to solve ths problem, we need to know what happens at the bubble pont; at ths pont, the frst bubble s formed, so we can assume that all of the number

More information

Math1110 (Spring 2009) Prelim 3 - Solutions

Math1110 (Spring 2009) Prelim 3 - Solutions Math 1110 (Sprng 2009) Solutons to Prelm 3 (04/21/2009) 1 Queston 1. (16 ponts) Short answer. Math1110 (Sprng 2009) Prelm 3 - Solutons x a 1 (a) (4 ponts) Please evaluate lm, where a and b are postve numbers.

More information

Sampling Theory MODULE V LECTURE - 17 RATIO AND PRODUCT METHODS OF ESTIMATION

Sampling Theory MODULE V LECTURE - 17 RATIO AND PRODUCT METHODS OF ESTIMATION Samplng Theory MODULE V LECTURE - 7 RATIO AND PRODUCT METHODS OF ESTIMATION DR. SHALABH DEPARTMENT OF MATHEMATICS AND STATISTICS INDIAN INSTITUTE OF TECHNOLOG KANPUR Propertes of separate rato estmator:

More information

TURBULENT FLOW A BEGINNER S APPROACH. Tony Saad March

TURBULENT FLOW A BEGINNER S APPROACH. Tony Saad March TURBULENT FLOW A BEGINNER S APPROACH Tony Saad March 2004 http://tsaad.uts.edu - tsaad@uts.edu CONTENTS Introducton Random processes The energy cascade mechansm The Kolmogorov hypotheses The closure problem

More information

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Problem Points Score Total 100

Problem Points Score Total 100 Physcs 450 Solutons of Sample Exam I Problem Ponts Score 1 8 15 3 17 4 0 5 0 Total 100 All wor must be shown n order to receve full credt. Wor must be legble and comprehensble wth answers clearly ndcated.

More information

GeoSteamNet: 2. STEAM FLOW SIMULATION IN A PIPELINE

GeoSteamNet: 2. STEAM FLOW SIMULATION IN A PIPELINE PROCEEDINGS, Thrty-Ffth Workshop on Geothermal Reservor Engneerng Stanford Unversty, Stanford, Calforna, February 1-3, 010 SGP-TR-188 GeoSteamNet:. STEAM FLOW SIMULATION IN A PIPELINE Mahendra P. Verma

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve

More information

where I = (n x n) diagonal identity matrix with diagonal elements = 1 and off-diagonal elements = 0; and σ 2 e = variance of (Y X).

where I = (n x n) diagonal identity matrix with diagonal elements = 1 and off-diagonal elements = 0; and σ 2 e = variance of (Y X). 11.4.1 Estmaton of Multple Regresson Coeffcents In multple lnear regresson, we essentally solve n equatons for the p unnown parameters. hus n must e equal to or greater than p and n practce n should e

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017) Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube A Numercal Study of Heat ransfer and Flud Flow past Sngle ube ZEINAB SAYED ABDEL-REHIM Mechancal Engneerng Natonal Research Center El-Bohos Street, Dokk, Gza EGYP abdelrehmz@yahoo.com Abstract: - A numercal

More information

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω

More information

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of

More information

How Differential Equations Arise. Newton s Second Law of Motion

How Differential Equations Arise. Newton s Second Law of Motion page 1 CHAPTER 1 Frst-Order Dfferental Equatons Among all of the mathematcal dscplnes the theory of dfferental equatons s the most mportant. It furnshes the explanaton of all those elementary manfestatons

More information

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY The Bot-Savart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot- Savart formula

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO

FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO ISTP-,, PRAGUE TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO Mohammad Rahnama*, Seyed-Mad Hasheman*, Mousa Farhad**

More information

Mathematical modeling for finding the thermal conductivity of solid materials

Mathematical modeling for finding the thermal conductivity of solid materials Mathematcal modelng for fndng the thermal conductvty of sold materals Farhan Babu 1, Akhlesh Lodwal 1 PG Scholar, Assstant Professor Mechancal Engneerng Department Dev AhlyaVshwavdyalaya, Indore, Inda

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Turbulence and its Modelling

Turbulence and its Modelling School of Mechancal Aerospace and Cvl Engneerng 3rd Year Flud Mechancs Introducton In earler lectures we have consdered how flow nstabltes develop, and noted that above some crtcal Reynolds number flows

More information

Aerodynamics. Finite Wings Lifting line theory Glauert s method

Aerodynamics. Finite Wings Lifting line theory Glauert s method α ( y) l Γ( y) r ( y) V c( y) β b 4 V Glauert s method b ( y) + r dy dγ y y dy Soluton procedure that transforms the lftng lne ntegro-dfferental equaton nto a system of algebrac equatons - Restrcted to

More information

Statistics Chapter 4

Statistics Chapter 4 Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment

More information

13. One way of expressing the power dissipated by a resistor is P = ( V)

13. One way of expressing the power dissipated by a resistor is P = ( V) Current and esstance 9. One way of expressng the power dsspated by a resstor s ( ). Thus, f the potental dfference across the resstor s doubled, the power wll be ncreased by a factor of 4, to a value of

More information

% & 5.3 PRACTICAL APPLICATIONS. Given system, (49) , determine the Boolean Function, , in such a way that we always have expression: " Y1 = Y2

% & 5.3 PRACTICAL APPLICATIONS. Given system, (49) , determine the Boolean Function, , in such a way that we always have expression:  Y1 = Y2 5.3 PRACTICAL APPLICATIONS st EXAMPLE: Gven system, (49) & K K Y XvX 3 ( 2 & X ), determne the Boolean Functon, Y2 X2 & X 3 v X " X3 (X2,X)", n such a way that we always have expresson: " Y Y2 " (50).

More information

G4023 Mid-Term Exam #1 Solutions

G4023 Mid-Term Exam #1 Solutions Exam1Solutons.nb 1 G03 Md-Term Exam #1 Solutons 1-Oct-0, 1:10 p.m to :5 p.m n 1 Pupn Ths exam s open-book, open-notes. You may also use prnt-outs of the homework solutons and a calculator. 1 (30 ponts,

More information

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram Adabatc Sorpton of Ammona-Water System and Depctng n p-t-x Dagram J. POSPISIL, Z. SKALA Faculty of Mechancal Engneerng Brno Unversty of Technology Techncka 2, Brno 61669 CZECH REPUBLIC Abstract: - Absorpton

More information

Implicit Integration Henyey Method

Implicit Integration Henyey Method Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure

More information

Exercises of Fundamentals of Chemical Processes

Exercises of Fundamentals of Chemical Processes Department of Energ Poltecnco d Mlano a Lambruschn 4 2056 MILANO Exercses of undamentals of Chemcal Processes Prof. Ganpero Gropp Exercse 7 ) Estmaton of the composton of the streams at the ext of an sothermal

More information

CHEMICAL ENGINEERING

CHEMICAL ENGINEERING Postal Correspondence GATE & PSUs -MT To Buy Postal Correspondence Packages call at 0-9990657855 1 TABLE OF CONTENT S. No. Ttle Page no. 1. Introducton 3 2. Dffuson 10 3. Dryng and Humdfcaton 24 4. Absorpton

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity 1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum

More information

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists *

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists * How Strong Are Weak Patents? Joseph Farrell and Carl Shapro Supplementary Materal Lcensng Probablstc Patents to Cournot Olgopolsts * September 007 We study here the specal case n whch downstream competton

More information

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Georgia Tech PHYS 6124 Mathematical Methods of Physics I Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends

More information

Chapter Newton s Method

Chapter Newton s Method Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve

More information

8.592J: Solutions for Assignment 7 Spring 2005

8.592J: Solutions for Assignment 7 Spring 2005 8.59J: Solutons for Assgnment 7 Sprng 5 Problem 1 (a) A flament of length l can be created by addton of a monomer to one of length l 1 (at rate a) or removal of a monomer from a flament of length l + 1

More information

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

DECOUPLING THEORY HW2

DECOUPLING THEORY HW2 8.8 DECOUPLIG THEORY HW2 DOGHAO WAG DATE:OCT. 3 207 Problem We shall start by reformulatng the problem. Denote by δ S n the delta functon that s evenly dstrbuted at the n ) dmensonal unt sphere. As a temporal

More information

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11) Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

ME 440 Aerospace Engineering Fundamentals

ME 440 Aerospace Engineering Fundamentals Fall 006 ME 440 Aerosace Engneerng Fundamentals roulson hrust Jet Engne F m( & Rocket Engne F m & F ρ A - n ) ρ A he basc rncle nsde the engne s to convert the ressure and thermal energy of the workng

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY. Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR

More information

8.6 The Complex Number System

8.6 The Complex Number System 8.6 The Complex Number System Earler n the chapter, we mentoned that we cannot have a negatve under a square root, snce the square of any postve or negatve number s always postve. In ths secton we want

More information