Transcript 열전달 수업자료
( Heat Transfer ) 2013년도 2학기 김 내 현 교수 Homepage : http://lion.incheon.ac.kr/~knh0001 열전달의 중요성 - 에어컨 QH 응축기 압축기 팽창장치 W 증발기 QL QL 을 구현하기 위한 증발기 설계 QH 4 Q L QL mh W -> 열전달 지식 필요 열전달의 중요성 - 전자 장치 냉각 c h i p • 칩의 온도(T)를 제어하기 위한 냉각장치 설계 -> 열전달 지식 필요 4 - 기타 에너지 관련기기 (발전소,자동차, 연료전지등)의 설계에 필요 열전달 메카니즘 전도(conduction) : 고체 내부 (자유전자, 격자 떨림 등) 열전달 대류(convection) : 유체유동에 의한 열전달 복사(radiation) : 빛 ( 전자파 또는 photon)에 의한 열전달 4 전 도(conduction) • Fourier 의 열전도 법칙 ( 면적 )( 온도차 ) 전도열전달량 두 께 T1 T2 T Q kA kA x x dT q Q / A k dx 4 열전도도 (Thermal conductivity) 4 Figure 1-6 Figure 1-7 열 확 산 율(Thermal diffusivity) Heat conducted k Heat stored C p 4 (m2 / s) 대 류(convection) • Newton 의 냉각 법칙 Q hAs Ts T 4 Table 1-5 복 사(radiation) • Stephan-Boltzmann 의 법칙 Q AsT 4 s • 실제 표면 Q AsT 4 s 4 예 제 1-1 Q kA T1 T2 L 0 (15 4) C 0 2 (0.8W / m C )(48m ) 0.25m 1690W 4 예 제 1-2 1 1 We (44W ) 22W 2 2 1 1 A D 2 (0.05m) 2 0.00196m2 4 4 Q QL (22W )(0.03m) k AT (0.00196m 2 )(15 oC ) 22.4W / m C o 4 예 제 1-4 Q VI (60V )(1.5 A) 90W As DL (0.003m)(2m) 0.01885m2 Q 90W 2o h 34.9 W / m C 2 o As (Ts T ) 4 (0.01885m )(152 15) C 예 제 1-4 Qw As (Ts 4 Tw4 ) (0.95)(5.67 108W / m2 K 4 )(1.4m2 ) [(30 273)4 (10 273)4 ]K 4 152W 4 예 제 1-6 Qconv hAs (Ts T ) (6W / m2 oC)(1.6m2 )(29 20) oC 86.4W Qrad As (Ts 4 T4 ) (0.95)(5.67 108W / m2 K 4 )(1.6m2 ) [(29 273)4 (20 273)4 ]K 4 81.7W 4 Qtotal Qconv Qrad (86.4 81.7)W 168.1W 예 제 1-7 (a)공기층 Qcond T1 T2 kA 219W L Qrad A(T14 T24 ) 368W Qtotal 587W (b)진공 Qrad 368W (c)우레탄 Qcond kA 4 T1 T2 260W L Qcond 0 Qrad 0 (d )초단열재 Qcond kA T1 T2 0.2W L Qrad 0 예 제 1-9 As qsolar hAs (Ts T ) qsolar 0.6 (700W / m2 ) Ts T 25 h 50W / m2 oC 33.4 oC 4 2-1 평면 벽에서 열전도 Qcond dT kA const dx T2 L Q cond 0 Qcond 4 dx kAdT T1 T1 T2 kA L 열저항(thermal resistance) 개념 Qcond R L : 전도저항 kA Qconv 4 R T1 T2 R Ts T R 1 : 대류저항 hAs 열저항 회로 Rtotal Rconv ,1 Rwall Rconv ,2 Q T Rtotal 4 1 L 1 h1 A kA h2 A Q UAT U : 열관류율( W / m2 oC ) 다층 평면 벽 4 예제 2-3 Ri 1 1 0.08333 oC / W 2 o 2 h1 A (10W / m C )(1.2m ) R1 R3 R2 L1 0.004m o 0.00427 C /W 2o 2 k1 A (0.78W / m C )(1.2m ) L2 0.01m o 0.3205 C /W 2o 2 k2 A (0.026W / m C )(1.2m ) 1 1 o Ro 0.02083 C /W h2 A (40W / m2 oC )(1.2m2 ) Rtotal Ri R1 R2 R3 Ro 0.4332 C / W 4 T1 T 2 [20 (10)] oC Q 69.2W o Rtotal 0.4332 C / W 예제 2-21 Rtotal Ri R1 R2 Rmid R6 Ro 1 1 1 1 Rmid R3 R4 R5 4 Q kA 원형관 r2 r1 dT dr T 2 Q dr kdT A T1 A 2 rL Q 2 Lk R ln( r2 / r1 ) 2 Lk Rtotal ln(r2 r1 ) 1 1 (2 r1 L)h1 2 Lk (2 r2 L)h2 구의경우 4 T1 T2 ln( r2 / r1 ) A=4 r 2 R r2 r1 4 r1r2 k Rtotal r2 r1 1 1 (4 r12 )h1 4 r1r2 k (4 r2 2 )h2 예제 2-7 A1 D12 (3m)2 28.3m2 A2 D22 (3.04m)2 29.0m2 Qrad (T24 T24 ) (T24 T24 )(T2 T2 )(T2 T2 ) 표면온도 5℃가정 hrad (1)(5.67 108W / m2 K 4 )[(295K )2 (278K )2 ][(295 278) K ] 5.34W / m2 K Ri 1 1 0.000442o C / W 2o 2 h1 A1 (80W / m C )(28.3m ) R1 r2 r1 0.000047o C / W 4 kr1r2 1 0.00345o C / W h2 A2 1 0.00646o C / W hrad A2 R0 Rrad 4 Rtotal 0.00274o C / W T 2 T1 (22 0)o C Q 8029W Ttotal 0.00274o C / W 예제 2-7 표면온도 재 계산 T 2 T1 Q Requiv 1 Requiv 1 1 444.7W oC R0 Rrad T2 T2 QRequiv 4C 24시간동안 열전달량 (8.029KJ / s)(24 3600) 673,700KJ 4 mice Q 673, 700KJ 2079kg hif 333.7 KJ 임계 단열 반지름 Q T1 T T1 T ln( r2 / r1 ) 1 Rims Rconv 2 Lk h(2 r2 L) dQ 0 dr2 구형관 rcr 4 rcr k h 2k h 예제 2-9 Rconv Rcond Q 1 1 0.76m2 C 2 2 hA2 (12W / m C )(0.110m ) ln(r2 / r1 ) ln(3.5 /1.5) 0.18m2 C 2 2 2 kL 2 (0.15W / m C )(5m ) T1 T Rtot T1 105C 단열이 없을 때 Rconv 1 1 2 1.76 m C 2 2 hA1 (12W / m C )(0.0471m ) T1 T QRTOT 30C (80W )(1.76C / W ) 4 170.8C 휜으로부터 열전달 Qcond , x Qcond , xx Qconv Qconv h( Px)(T T ) Qcond , x x Qcond , x x hP(T T ) 0 dQcond hP (T T ) 0 dx Qcond kAc dT dx d dT kA c hP(T T ) 0 dx dx 휜으로부터 열전달 단면과 열전도도가 일정한 경우 d 2 a 2 2 dx T T 일반해 경계조건 a2 hP kAc ( x) C1eax C2eax (o) b Tb T 휜으로부터 열전달 1. 무한히 긴 휜 ( L) 0 beax Q kAc dT dx 또는 T T e ax Tb T hPkAc (Tb T ) x 0 2. 휜 끝에서 열손실 무시 d dx 0 xL T T cosh a( L x) Tb T cosh aL Q kAc dT dx hPkAc (Tb T ) tanh aL x 0 3. 대류가 있는경우 : Lc L Ac P 휜 효율 f 실제 전열량 최대 전열량 hPkAc (Tb T ) tanh mL tanh mL hAf (Tb T ) mL 휜 유효도 f Q fin Qunfin f hAf (Tb T ) hAb (Tb T ) Af Ab f 휜 달린 표면 Qtotal Qfin Qunfin f hAf (Tb T ) hAuf (Tb T ) 예제 2-11 Q R T R T (90 30)C 1.0C / W Q 60W 예제 2-12 휜효율 계산 r2 r1 1 t 2 2.07 1 L t 2 h 0.207 kt f 0.95 Af 2 (r22 r12 ) 2 r2t 0.00462m2 Qfin f hAf (Tb T ) 25.0W Auf D1S 0.000283m2 Qunfin hAuf (Tb T ) 1.6W Qtotal n(Qfin Qunfin ) 200(25.0 1.6)W 5320W 전도형상계수 Q Sk (T1 T2 ) S 2 L 62.9m ln(4 z / D) Q Sk (T1 T2 ) (62.9m)(0.9W / mC )(80 10)C 3963W 열전도 방정식 ◎ ( 왼쪽 면을 통한 전도 열 ) + ( 요소 안의 발생 열 ) – ( 오른쪽 면을 통한 전 도열) = ( 내부 열 축척 ) ◎ 왼쪽 면 전도 열 : T q x kA x ◎ 오른쪽 면 전도 열 : qx dx T kA x x dx T T A k k x x x 4 dx 열전도 방정식 ◎ 요소안의 발생 열 : ◎ 내부 열축척 : qAdx cA T dx T T T T kA qAdx A k k dx dx cA x x x x T T k q c x x 4 3차원 방정식 T k x x T k y y T k z z T q c ◎ K가 상수일 경우 2T 2T 2T q 1 T 2 2 2 x y z k k : 열 확산 계수 c ◎ 원통 좌표계 2T 1 T 4 1 2T 2T q 1 T 2 2 2 2 r r r r z k 평판의 2차원 열전도 2T 2T 2 0 2 x y ◎ 경계조건 y x x y = 0, = 0, = w, = H, T = T1 T = T1 T = T1 x T Tm sin T1 w ◎해 T T1 2 (1) n 1 1 n x sinh( n y / W ) sin 4 T2 T1 n w sinh(n H / W ) n 1 T x 수치해석 T x T y T y T x 2 2 2T y 2 m,n m,n T x 1 m ,n 2 T x 1 m ,n 2 x Tm4 , m 1 Tm , n 1 2Tm , n y 2 1 m ,n 2 1 m ,n 2 m,n 1 2 1 m,n 2 Tm 1, n Tm , n x Tm , n Tm 1, n x Tm , n 1 , Tm, n y Tm, n Tm, n 1 y Tm 1, n Tm 1, n 2Tm , n ( x ) 2 Tm1,n Tm1,n 2Tm,n (x) 2 Tm,n1 Tm,n1 2Tm,n (y) 2 0 ◎ △x = △y 인 경우 Tm1,n Tm1,n Tm,n1 Tm,n1 4Tm,n 0 4 100 + 500 + T2 + T3 - 4T1 = 0 T1 + 500 + 100 + T4 - 4T2 = 0 100 + T1 + T4 + 100 - 4T3 = 0 T1 = T2, T3 = T4 이므로 100 + 500 + T3 - 3T1 = 0 100 + T1 + 100- 3T3 = 0 T1 = T2 = 250 , T3 = T4 = 150 4 대류 경계조건 k y Tm,n Tm1,n x hy(Tm,n T ) k x Tm,n Tm,n 1 x Tm,n Tm,n 1 k 2 y 2 y 1 hx hx Tm,n 2 T (2Tm1,n Tm,n 1 Tm,n 1 ) 0 2 k k y Tm ,n Tm 1,n x Tm ,n Tm ,n 1 k k 2 x 2 y x y h Tm, n T h Tm,n T 2 2 4 hx hx 2Tm,n 1 2 T Tm1,n Tm,n 1 0 k k 예제 3-5 ◎1, 2, 4, 5 절 : Tm1,n Tm1,n Tm,n1 Tm,n1 4Tm,n 0 500 + 100 + T2 + T4 – 4T1 = 0 500 + T1 + T5 + T3 - 4T2 = 0 T1 + 100 + T5 + T7 - 4T4 = 0 T2 + T4 + T6 + T8 - 4T5 = 0 1 hx hx 2 T (2Tm1,n Tm,n 1 Tm,n 1 ) 0 ◎ 3, 6, 7, 8 절 : Tm,n k 2 k 2T2 2T5 2T5 2T5 + + + + T6 T3 T7 T7 + + + + 567 - 4.67T3 = 0 T9 + 67 – 4.67T6 = 0 T9 + 67 - 4.67T8 = 0 T9 + 67 - 4.67T8 = 0 4 ◎ 9절 : T6+T8+67-2.67T9=0 4 1 0 1 0 0 0 0 0 1 4 0 1 2 4.67 0 0 1 0 0 1 0 0 0 0 0 0 4 1 0 0 1 0 4 1 0 2 0 0 0 1 4 2 0 2 0 0 T1 600 0 T2 500 1 0 0 0 T3 567 0 0 0 0 T4 100 1 1 1 0 T5 0 4.67 0 0 1 T6 67 0 1 1 0 T7 167 0 4.67 4.67 1 T8 67 1 1 1 2.67 T9 67 0 0 0 0 0 0 ◎해 4 절점 온도 1 280.67 2 330.30 3 309.38 4 192.38 5 231.15 6 217.19 7 157.70 8 184.71 9 175.62 4 집중계 해석 hAs (T T )dt mc p dT hAs d (T T ) dt T T VC p hAs T T ln t Ti T VC p T T e Ti T hAs t VC p 집중계 해석 기초 hLc Bi k V Lc As Bi 0.1인 경우 T- T오차 5%미만 예제 3-1 1 D3 V 6 1 1 4 Lc = D (0.001 m ) 1.67 10 m 2 As D 6 6 hLc (210W / m2 C )(1.67 104 m) Bi = 0.001 0.1 k 35W / mC T T 0.01 Ti T hAs h 210W / m2 C b VC p C p Lc (8500kg / m3 )(320 J / kg C )(1.67 104 m) 0.462 S 1 0.01 e0.462t t 10 (s) 비정상 열전도 T T Ti T x X L hL Bi ( Biot Number ) k t 2 L ( Fourier Number ) 예제 3-3 hro (1200W / m2 C )(0.025m) Bi 47.8 k 0.627W / mC To T 70 95 0.278 Ti T 5 95 표에서 =0. 209 ro 2 (0.209)(0.025m)2 t 865s 2 0.15110 6m / s 반무한 고체 예제 3-6 h t k 1 T T 0 (10) 1 0.6 Ti T 15 (10) 표에서 = x 2 t 0.36 t (90days)(24h / day)(3600s / h) 7.78 106 s x 2 t 2 0.36 (0.15 106 m 2 / s )(7.78 106 s) 0.77m 다차원 비정상 열전도 T T T T T T T T T T T T wall i cyl i i Q Qmax Q Q Q 1 Qmax 1 Qmax 2 Qmax 1 예제 3-7 평판 t (3.39 105 m2 / s )(900s ) 2 8.48 L (0.06) 2 1 k 110W / mC 30.6 Bi hL (60W / m2 C )(0.06m) To T 0.8 T T plate i 원통에서도 마찬가지로 To T 0.5 T T cyl i 따라서 To T 0.8 0.5 0.4 Ti T To T 0.4(Ti T ) 25 0.4(120 25) 63C 예제 3-8 m ro 2 L (8530kg / m3 ) (0.05m)2 (0.06m) 4.02kg Qmax mCp (Ti T ) (4.02kg )(0.380kJ / kgC)(120 25)C 145.1kJ 평판 Bi 1 0.0327 30.6 Bi 2 (0.0327)2 (8.48) 0.0091 표에서 Q / Qmax 0.23 원통 Q / Qmax 0.47 Q Q Q Q 1 Q max Qmax Qmax 1 Qmax 2 1 0.23 0.47(1 0.23) 0.592 Q 0.592Qmax 0.592 145.1kJ 85.9kJ 제 5장 대류의 원리 층류 : Re x u x u x 5 105 경계층 : 자유흐름속도의 99%까지 거리 난류 : Re D um d 2300 열 경계층 열 경계층 : 자유흐름속도의 99%까지 거리 열경계층 속도경계층 ( )와 열경계층 (t )의 비 C p Pr k t 액체금속(Pr <<1) t 오일(Pr 1) t 물, 공기 Pr 1 열경계층 국소 열전달계수 hxx 층류 : Nu x = 0.332 Re x 0.5 Pr1/ 3 k 난류: Nux =0.0296Rex0.8 Pr1/ 3 1 평균 열전달계수 h L L 0 hx dx hL 0.664 Re L 0.5 Pr1/ 3 k 난류: Nu=0.037 ReL0.8 Pr1/ 3 층류 : Nu= 임계 Reynolds 수 Rec 5 105 층류 영역을 고려한 경우 1 h L xcr 0 hx ,la min ar dx L xcr hx ,turbulent dx Nu (0.037 ReL0.8 871) Pr1/3 균일 열유속 층류 : Nu x =0.453Rex 0.5Pr1/3 난류 : Nu x =0.0308Rex 0.8Pr1/3 예 5-4 27℃ 1기압의 공기가 2m/s의 속도로 평판위를 흐르고 있다. 평판 전체가 60℃ 로 가열되었다고 하자. 평판의 선단으로부터 (a) 20cm와 (b) 40cm 되는 곳까지 의 열전달을 계산하여라. Tf 27 60 43.5C 316.5 K 110.3 F 2 물성치 17.36 106 m 2 / s 1.87 10 4 ft 2 / s k 0.02749W / m C 0.0159 Btu / h ft F Pr 0.7 C p 1.006kJ / kg C 0.24 Btu / lbm F 예 5-4 (a) x=20cm Re x u x (2)(0.2) 23, 041 17.36 106 1 1 1 1 hx x Nu x 0.332 Re x 2 Pr 3 (0.332)(23, 041) 2 (0.7) 3 44.74 k k (44.74)(0.02749) hx Nu x 0.2 x 6.15W / m 2 C 1.083Btu / h ft 2 F h (2)(6.15) 12.3W / m 2 C 2.17 Btu / h ft 2 F q hA(Tw T ) q (12.3)(0.2)(60 27) 81.18W 277 Btu / h 예 5-4 (b) x=40cm Re x u x (2)(0.4) 46, 082 17.36 106 1 1 1 1 hx x 3 2 2 Nu x 0.332 Re x Pr (0.332)(46, 082) (0.7) 3 63.28 k k (63.28)(0.02749) hx Nu x 0.4 x 4.349W / m 2 C h (2)(4.349) 8.698W / m 2 C 1.53Btu / h ft 2 F q hA(Tw T ) q (8.698)(0.4)(60 27) 114.8W 392 Btu / h 5.2 운동방정식 - 연속 방정식 uy ( ) xy y x y x x ( u ) uy yx x ( u) ( ) xy xy 0 x y ( u ) ( ) 0 x y 가 일정하면 u x y 0 -운동량 방정식 Newton의 제 2법칙 d (mV ) F dt i 방향의 i 방향의 i 방향의 운동량 변화 체적력 표면력 1) 운동량 변화 ( u ) dxdy x u dy ( uu ) dxdy y uudy uudy u dx ( uu ) ( u ) M out M in dxdy dxdy x y u u u dxdy y x 2) 체적력 및 표면력 yx dx Px dy ( yx ) y dydx Fx Px dy yx dx yx u y P 2u 표면력+체적력= Fx dxdy 2 x y P dxdy x u u P 2u x 방향 : ( u ) Fx 2 x y x y P 2 y 방향: ( u ) Fy x y y x 2 원통좌표계: ( u u u P u ) Fx r x r x r r r - 에너지 방정식 대류에 의한 에너지 입력 체적력에 의한 표면력에 의한 전도에 의한 에너지 입력 에너지입력 에너지입력 T T 2T cP u k 2 y y x 5.4 강제유동 F f (V , D, , ) M M ML L f , L, 3 , 2 L LT T T (길이- L, 시간- T, 질량- M , 온도- ) 질량제거 f1 V , D, F L2 T L f , L , 1 T L2 T 시간제거 F V f 2 D, V 1 L f L , 2 L 길이제거 VD F f3 f 3 Re VD 5.6 강제대류 h f ( D, V , , , c p , k ) L M M L2 ML M T 3 f L, T , L3 , LT , T 2 , T 3 질량제거 k f1 D, V , , cP , h L3 L L2 L2 L4 T 3 f1 L, T , T , T 2 , T 2 온도제거 h k f 2 D, V , , cP cP L L2 L2 L T f 2 L, T , T , T 시간제거 h V V f 3 D, , cPV 1 1 0 f 3 L, , L L 길이제거 h VD VD f4 , cPV St f 4 (Re, Pe) St hD : prandtl 수, : Nusselt 수 k Nu , Pe Re Pr Re Pr 또는 Nu=f 4 (Re, Pr) 5.7 자연대류 d 부력 ~ g ~ g ~ g T dT 1 d : 열팽창계수 dT h=f ( D, g T , , , cP , k ) L M M L2 ML M T 3 f L, T 2 , T 3 , LT , T 2 , T 3 질량제거 h =f1( D, g T , , , cP ) k k k L T 3 T 2 L2 1 L f1 L, T 2 , L4 , L2 , T 2 온도제거 h =f 2 ( D, g T , , Pr) k L L2 1 f L , , , 0 2 T2 T L 시간제거 h g T =f 3( D, , Pr) k 2 1 1 f L , , 0 3 3 L L 길이제거 hD g TD 3 =f 4 ( , Pr) k 2 Nu f ( Ra , Pr) g TD 3 Ra : Rayleigh 수 v g TD 3 v g TD 3 / 2 v 강제대류 Qconv hAs (Ts T ) qconv h(Ts T ) k fluid h T y y 0 k fluid T / y y 0 Ts T 예제 4-1 물성치는 막(film)온도에서 구한다 876kg / m3 , Pr=2870 k=0.144W/mC, =242 10-6 m2 /s Re L VL 2m / s 5m 0.242 105 m2 / s 4.13 104 Nu 0.664 ReL0.5 Pr1/ 3 0.664 (4.13104 )0.5 (2870)1/ 3 1918 h k 0.144W / mC Nu (1918) 55.2W / m 2 C L 5m Q hAs (T Ts ) 55.2W / m2 C)(5 1m2 )(60 20) 11.040W 원통과 구를 지나는 유동 5/8 hD 0.62 Re Pr Re Nucyl 0.3 1 1/ 4 2/ 3 k 1 0.4 / Pr 282,000 1/ 2 1/ 3 4/ 5 1/ 4 Nusph hD 2 0.4 Re1/ 2 0.06 Re 2 / 3 Pr 0.4 k s 예제 4-4 막온도 60℃의 물성치 k 0.02808W / mC , Pr=0.72 =1.896 10-5 m2 /s Re Nu 0.3 h VD 8m / s 0.1m 5 1.896 10 m / s 1/ 2 1/ 3 0.62 Re Pr 2/ 3 1/ 4 1 0.4 / Pr 2 4.219 104 Re 5/8 1 282,000 4/ 5 124 k 0.02808W / mC Nu 124 34.8W / m2 C D 0.1m As DL 0.1m1m 0.314m2 Q hAs (Ts T ) 34.8W / m 2 C 0.314m 2 110 10 1093W 파이프유동 Vm D Re 층류: Re 2300 난류 : Re>2300 C T( V2 r dr ) T= p m 2 V( R )C p m 2 Vm R 2 R TVrdr 0 입구영역 층류 : Lh =0. 05ReD Lt 0.05Re Pr D 난류: Lh Lt 10D 일정 열유속 Q qs As mC p (Te Ti ) Te Ti qs As m Cp dTs dTm T const X dX dx mC p dTm qs (2 R)dx dTm qs (2 R ) dX m Cp 일정 온도 m C p dTm h(Ts Tm )dAs d (Ts Tm ) hP dx Ts Tm mC p ln Ts Te hAs Ts Ti m Cp Te Ts (Ts Ti ) exp hAs / m C p Q mC p (Te Ti ) hAs T ln T ln Ti Te ln Ts Te / Ts Ti Te Ti ln Te / Ti 예제 4-6 물성치: Cp =4187J / KgC Q mC p (Te Ti ) (0.3kg / s)(4.187 KJ / Kg C )(115C 15C ) 125.6kW Tln Te Ti ln Te Ti 5 105 32.85C 5 ln 105 Q 125.6kW 2 As 4.78 m hTln (0.8kW / m2 C )(32.85C ) As 4.78m2 L 61m D (0.025m) 관내 층류 유동 mCpTx mCpTxdx Qr Qr dr 0 m V (2 rdr ) C pV Tx dx Tx Qr dr Qr 1 dx 2 rdr dr V T 1 Q x 2 C p rdx r Q T T ( k 2 rdx ) 2 kdx( r ) r r r r V 4 T T (r ) x r r r k C p 일정 표면 열유속 Ts Tm T 2q const x x x VmC p R r2 V 2Vm (1 2 ) R 4q r2 1 d dT (1 2 ) (r ) kR R r dr dr q r2 2 T (r ) C1 R C2 2 kR 4R T =0 r r R T =Ts B.C. : r 0 4 Q q 2 Rdx qR 3 r 2 r4 T Ts ( 2 ) 4 k 4 R 4R 2 Tm Vm R 2 Tm R 0 TVrdr 11 qR Ts 24 k q h(Ts Tm )이므로 24 k 48 k k h 4.36 11 R 11 D D Nu 4.36 4 관내 층류 유동의 속도 분포 ` (2 rdrP) x (2 rdrP) xdx (2 rdr )r (2 rdr )r dr 0 r Px dx Px ( r ) r dr ( r ) x r 0 dx dr dP d (r ) r r 0 dx dr d r dx (r du dP ) dr dx 2 w dP dx R 1 dP u (r ) ( ) C1 ln r C2 4 dx vavg R 2 dP r2 u (r ) ( )(1 2 ) C1 ln r C2 4 4 dx R 2 R 2 R R 2 dP r2 R 2 dP 2 u (r )rdr 2 ( )(1 2 )rdr ( ) 0 0 R R 4 dx R 8 dx 2 r u (r ) 2Vavg (1 2 ) R 비원형관내 층류 유동 4 입구영역에서 발달중인 층류유동 • 원형관 0.065( D / L) Re Pr Nu 3.66 1 0.04[( D / L) Re Pr]2/3 • 평행 평판 0.03( Dh / L) Re Pr Nu 7.54 1 0.016[( Dh / L) Re Pr]2/3 4 Table 1-5 예4-7 물성 치 =888kg/ m3 901106 m2 / s k 0.145W / m℃ c p 1880J / kgK Pr 10, 400 Re Vm Dh (2m / s)(0.3m) 666 901106 m2 / s Lt 0.05Re Pr D (0.05)(666)(10, 400)(0.3m) 104, 000m hD 0.065( D / L) Re Pr 3.66 k 1 0.04[( D / L) Re Pr)]2 / 3 37.3 Nu h k 0.145W / m Nu (37.3) 18.0W / m℃ D 0.3m As DL 188.5m2 m AcVm (888kg / m3 )[ (0.3m) 2 ](2m / s) 4 125.5kg / s Te Ts (Ts Ti )exp(hAs / mcp ) 19.71℃ Tlm Ti Te 19.85℃ Ts Te ln Ts Ti Q hAs Tlm (18.0W / m℃) (188.5m2 )(19.85℃) 6.74 104W f 64 0.0961 Re L Vm 2 P f 1.14 105 N / m2 D 2 Wpump mP (125.5kg / s)(1.14 105 N / m2 ) 16.1kw 888kg / m3 관내 난류 유동 f (0.79ln Re1.64)2 Nu ( f / 8)(Re 1000) Pr 1.0 12.7 f / 8(Pr 2/ 3 1) 또는 Nu 0.023Re0.8 Pr n (가열 : n=0.4) (냉각 : n=0.3) 액체금속 T 일정 : Nu 4.8 0.0156 Re0.85 Prs 0.93 q 일정 : Nu 6.3 0.0167 Re0.85 Prs 0.93 예4-8 물성치 =992. 1kg / m3 cp 4179J / kg℃ k =0.631W / m℃ Pr 4.32 0.658 106 m2 / s Q mcp (Te Ti ) 34.6kw Re Vm D Nu hD 0.023Re0.8 Pr 0.4 69.5 k 10, 760 h 1462W / m2℃ qs 73, 460W / m2 Ts Tm 65℃+ 115℃ h 1462W / m2℃ 자연대류 검사체적의 운동량 법칙 m ax Fx ax du x dx u dy u u u v dt u dt y dt x y m (dx dy) P Fx dy dx 1 dx (dy 1) g (dx dy 1) x y 2u P 2 g (dx dy 1) x y u y u u 2u P u v g 2 x y y x 외부의 정지한 유체 P g x P P 이므로 P P g x x u u 2u u v 2 ( ) g y y x u u 2u u v v g (T T ) x y y 2 (T T ) 체적팽창계수 ideal =- 1 T 1 P T P 지배방정식 무차원화 x* T T x y u v , y* , u* , v* , T * Lc Lc V V Ts T * g (Ts T ) Lc 3 T u * 1 2u * * u u v * * 2 2 x y v Re L y*2 Re L * g (Ts T ) Lc 3 GrL v2 GrL 109 : 난류 수직판의경우 자연대류와 강제대류의 상대적 중요성 GrL / Re L 2 1 : 강제대류 지배적 GrL / Re L 2 1 : 자연대류 지배적 GrL / Re L 2 1 : 두 모드 모두 중요 자연대류 상관식 N u=CRa Ln g ( Ts - T ) Lc 3 Ra L =GrLPr = Pr 2 v 예제 5-1 막온도에서의 물성치 k 0.0270W / mC Pr 0.724 v 1.75 105 m2 / s 특성길이 1 1 Tf 318 Lc =D=0. 08m g (Ts T ) D3 Pr 6 RaD 1.869 10 v2 0.387 RaD1/ 6 Nu 0.6 9 /16 8 / 27 1 (0.559 / Pr) h k Nu 5.869W / mC D As DL 1.508m2 Q hAs (Ts T ) 443W 복사 열전달 Q As (Ts 4 Tsurr 4 ) 553W 2 17.40 휜달린 표면의 자연대류 냉각(Ts=일정) Ts constant : hS 576 2.873 Nu= k ( Ras S / L) 2 ( Ras S / L)0.5 g (Ts T )S 3 Ras Pr, 2 v Ts constant : S Sopt : g (Ts T ) L3 L3 Ra L = Pr Ras 3 2 v S S3 L Sopt =2. 714 Ra s Nu= 0.5 hSopt 1.307 k 0.25 2.714 L RaL 0.25 수직 PCB의 자연대류 냉각(q=일정) hL S 48 2.51 Nu * * 0.4 k Ra S / L ( Ra S / L ) L s 0.5 4 g q S s Ras* Pr 2 kv q s constant : S L Sopt =2. 12 * Ra s 4 0.2 예제 5-3 물성치 k 0.02772W / mC 5 v 1.846 10 m / s 2 Pr 0.7215 1 1 Tf 328 g (Ts T ) L3 Pr RaL 1.846 107 2 v Sopt 2.714 n L 7.45mm RaL 0.25 W 15 fins st h Nu k 0.2012W / m2 C Sopt Q hAs (Ts T ) 1.30W 밀폐공간 Nu 0.159 RaL1/ 4 104 <Ra L <4 105 Nu 0.068RaL1/ 3 4 105 <Ra L <107 1708 RaL1/ 3 Nu 1 1.44 1 1 Ra 18 L 1708 Nu 1 1.44 1 RaL cos Ra L <108 1708(sin1.8 )1.6 ( RaL cos )1/ 3 1 1 RaL cos 18 Pr Nu 0.18 RaL 0.2 Pr 0.28 Pr H Nu 0.22 RaL 0.2 Pr L Q 2 keff ln( Do / Di ) 1 H / L 2 any Prandtl number 0.29 RaL Pr/(0.2 Pr) 103 2 H / L 10 any Prandtl number 1/ 4 RaL 1010 (Ti To ) 1/ 4 Pr 0.386 K 0.861 Pr keff ( Fcyl RaL )1/ 4 ln( Do / Di ) 4 Fcyl Lc 3 ( Di 3/ 5 Do 3/ 5 )5 Di Do Q keff Lc (Ti To ) 1/ 4 Pr 0.74 K 0.861 Pr keff Fsph ( Fsph RaL )1/ 4 Lc ( Di Do )4 ( Di 7 / 5 Do 7 / 5)5 예제 5-4 물성치 k 0.02416W / mC 5 v 1.399 10 m / s 2 Pr 0.734 1 1 Tave 280 g (T1 T2 ) L3 4 RaL 1.051 10 v2 H Nu 0.42RaL1/ 4 Pr 0.012 L 0.3 1.401 As 1.6m2 Q hAs (T1 T2 ) kNuAs T1 T2 27.1W L 비등열전달 풀비등 상관식 g l v l h fg 1/ 2 qnucleate C p Ts Tsat n Csf h fg Prl 3 최대 열유속 qmax Ccr h fg g v l v 2 1/ 4 최소 열유속 1/ 4 qmin g l v 0.09 v h fg 2 l v 막비등 gkv3 v l v h fg 0.4C pv Ts Tsat C film v D Ts Tsat 1/ 4 q film qrad Ts 4 Tsat 4 qtotal q film 3 qrad 4 Ts Tsat 예제 6-1 100C 물 0.0589 N / m l 957.9kg / m3 v 0.6kg / m3 h fg 2257 KJ / kg l 0.282 10 3kg m / s Prl 1.75 Cpl 4217 J / kg C Csf=0.013 n=1.0 Rohsenow 관계식 g l q l h fg 1/ 2 Cpl Ts Tsat n C h Pr sf fg l 7.2 104 W / m 2 Q q A 5093W mevapo Q 2.26 103 kg / s h fg 3 유동비등 막응축 du bdx v g y bdx dy l g y bdx l du g l v y dy l y 0 에서 y=y 까지적분 g l v y2 u( y) y l 2 m( x) l u ( y )dA A l u ( y )bdy y 0 gbl l v 3 m( x) 3l gbl l v 2 d dm dx l dx 막응축 Tsat Ts dQ h fg dm kl bdx 3 d l kl Tsat Ts dx g l l v h fg dm kl b Tsat Ts dx h fg 4 k T Ts x ( x) l l sat g h l l v fg 1/ 4 q( x) hx Tsat Ts kl Tsat Ts hx kl ( x) g l l v h fg kl 3 hx 4 l Tsat Ts x 1/ 4 h have 1 L 3 h dx 4 h x 0 g l l v h fg kl 3 0.943 l Tsat Ts L 1/ 4 L xL g l l v h fg *kl 3 0.943 T T L l sat s 1/ 4 hvert hfg* hfg 0.68Cpl Tsat Ts (W / m2 C ), 0 Re 30 액막의 Reynolds 수 Re Dh l vl l 4 Ac l vl 4 l vl 4m pl l p l Qcond hAs Tsat Ts mhfg* 4 As h Tsat Ts 4Qcond Re * pl h fg pl h fg * 4 g l l v 3 4 g l 2 Re 2 3l 3l 2 1/ 3 g hvert 1.47 kl Re 1/ 3 2 vl kl hx L 4g 3 3 h / 4 vert 0 Re 30 v l 3 파형 층류유동 ( Re>30) 30 Re 1800 1/ 3 hvert , wavy g Re kl 1.08 Re1.22 5.2 vl 2 Re vert , wavy v 1/ 3 3.70 Lkl Tsat Ts g 4.81 2 h * l fg vl l 0.820 , v l 난류유동 ( Re>1800) 1/ 3 hvert ,turbulent g Re kl 2 8750 58 Pr 0.5 Re0.75 253 vl Re 1800 v Revert ,turb l 0.0690 Lk Pr 0.5 T T g 1/ 3 l sat s 0.5 2 151Pr 253 h * l fg vl 3/ 4 예제 6-4 물성치 l =965. 3kg/ m3 Cpl =4206J / kgC l =0. 315 10- 3 kg/ m s 2 vl =l / l =0. 326 10 - 6 m/ s kl =0. 675W / mC h fg * h fg 0.68C pl Tsat Ts 2314 103 J / kg 1/ 3 3.70 Lkl Tsat Ts g Re 4.81 2 l h fg * vl 1287 0.820 1/ 3 h g Re kl 1.08 Re1.22 5.2 vl 2 5848W / m 2 C 예제 6-6 물성치 l =994kg/ m3 l =0. 720 10 - 3 kg/ m s Cpl =4178J / kg C kl =0. 623W / mC hfg* hfg 0.68Cpl Tsat Ts 2435 103 J / kg g l l v h fg *kl 3 h 0.725 l T T D sat s 1/ 4 = 9292W / m2 C 1/ 4 g l 2 h fg *kl 3 0.725 T T D l sat s 복사 Stephan Boltsman 법칙 Eb ,0 (T ) Eb ( , T )d 0 Eb ( , T ) C1 5 exp c2 / T 1 4 C1 2 hco 2 3.742 108 W m/ m2 C2 hco / k 1.439 104 m K ( T )max T power =2897. 8 m K Eb (T ) Eb (, T )d T 4 0 W/m 2 W/m W/m 2 2 예제 7-1 Eb T 4 23.2kW / m 2 5분동안 방사한 총 복사량 Qr ad =Eb As T=876KJ 파장 3 m에서의 흑체 방사도 C1 Eb = 3848W / m 2 m C 5 exp 2 1 T Eb,0 (T ) Eb ( , T )d W/m 2 0 E ( , T ) d b f (T ) 0 T 4 W/ m 2 예제 7-2 1T 1000 m K 2T 1900 m K f 1 2 f 1 f 2 0.0527 방사도가 가장 높은 파장 ( T ) max 2897.8 m K max 1.16 m 방사율 = 흡수율 E Ebl ack body 반사율 투과율 Gabs G G = r ef G Gt r = G = 1 형상계수 Fij : 표면 i를 떠나서표면 j에부딪히는 복사의분율 Ai Fij A j F ji N F j 1 ij 1 (밀폐공간) 예제 8-1 A1 F12 =A2 F21 A 4 r F21 = 1 F12 A2 4 r 2 1 2 2 r1 r2 F21 F22 1 r1 F22 1 F21 1 r2 2 2 예제 8-3 F12 F13 F14 F15 F ij 1 F12 F13 F14 F15 0.25 복사열전달 : 흑체 Q12 A1 Eb1 F12 A2 Eb 2 F21 = A1F12 T14 -T2 4 예제 8-6 Q12 A1F12 T14 -T2 4 F11 F12 F13 1 F13 1 F11 F12 1 0 0.2 0.8 Q13 394 103W Q12 1329 103W Q1 Q11 Q12 Q13 0 1319 394 925( kW ) 복사열전달 : 회체 Radiation emitted Radiation reflected Ji by surface i by surface i i Ebi i Gi i Ebi (1 i )Gi Ji Ebi Ti 4 복사열전달 : 회체 Q12,one shield Eb1 Eb 2 1 3,1 1 2 1 1 1 A11 A1 F12 A3 3,1 A2 2 복사열전달 : 회체 Qi j 표면j에 부딪치는 표면i에 부딪치는 완전한표면 i를 완전한표면 j를 떠나는 복사 떠나는 복사 Ai J i Fi j Ai J i Fj i Qi j Ai Fi j J i J j Qi j Ri j Q12 Ji J j Ri j 1 Ai Fi j (T14 T2 4 ) 1 1 1 2 1 A11 A1 F12 A2 2 예제 8-7 Q12 (T14 T2 4 ) q12 3625W / m 2 1 1 A 1 1 2 복사차폐 Q12,one shield Q12,one shield Q12, N shield Q12, N shield Eb1 Eb 2 1 1 1 1 3,1 1 3,2 1 1 2 A11 A1F12 A3 3,1 A3 3,2 A3 F32 A2 2 (T14 T2 4 ) 1 1 1 1 1 1 1 2 3,1 3,2 (T14 T2 4 ) 1 1 1 1 1 1 1 1 1 1 2 3,1 3,21 N ,1 N ,2 (T14 T2 4 ) 1 Q12,no shield N 1 1 1 N 1 1 1 2 예제 8-11 q12 (T14 T2 4 ) 1 1 1 1 1 1 1 2 3,1 3,2 806W / m2 예제 8-12 q h T f Tth th Tth 4 Tw 4 T f Tth th Tth 4 Tw 4 715K h 열교환기 열관류율 (총괄열전달계수) R Rtotal Ri Rwall Ro ln Do / Di 1 1 hi Ai 2 kL ho Ao 대수평균온도차법 Q mh C ph dTh Q mc C pc dTc dTh Q mh C ph dTc Q mc C ph dTh dTc d Th Tc 1 1 Q m m h C ph c C pc Q U Th Tc dAs 로그평균온도차법 d Th Tc 1 1 UdAs Th Tc mh C ph mc C pc ln Th,out Tc ,out Th,in Tc ,in mh ch 1 1 UAs m m h C ph c C pc q Th1 Th 2 mc cc Q UAs Tlm Tlm T1 T2 ln T1 / T2 q Tc1 Tc 2 수정계수의 사용 Q UAF Tlm,CF 예제 9-3 T1 Th ,in Tc ,out 8C T2 Th ,out Tc ,in 16C T1 T2 Tlm 11.5C T1 ln T2 Q UATlm 1087 kW 예제 9-3 Q mC p Tin Tout 62.93kW water Ai n Di L 0.408m2 Ui Q Ai F Tlm,CF Tlm,CF P T1 T2 47.6C T ln 1 T2 t2 t1 0.36, T1 t1 R T1 T2 0.8 t2 t1 F 0.97 Ui Q 3341W / m2 C Ai F Tlm,CF 유용도-NTU법 Q Q Qmax Cmin (Th ,in Tc ,in ) UA NTU Cmin c Cmin Cmax func( NTU , c) 예제 9-8 Q m C p Tout Tin 301kW water Ch m h C ph 8.62kW / C Cc m c C pc 5.02kW / C c C min Cmax 5.02 0.583 8.62 Q Cmin (Th ,in Tc ,in ) 0.428 1 1 ln 0.651 c 1 c 1 NTU Cmin A 5.11m 2 U A A DL L 108m D NTU