Chương 9: TỎA NHIỆT ĐỐI LƯU ( Khoa Công Nghệ Nhiệt Lạnh )

Chương 9: Tỏa nhiệt đối lưu Trang 77 Chương 9: TỎA NHIỆT ĐỐI LƯU 9.1 TỎA NHIỆT ĐỐI LƯU TỰ NHIÊN 9.1.1 Đặc trưng chuyển động tự nhiên Chuyển động tự nhiên chất lỏng dọc theo bề mặt xảy tồn trao đổi nhiệt Thông thường nhiệt độ thay đổi lớp mỏng gần sát bề mặt trao đổi nhiệt Còn trường tốc độ lớp mỏng có dạng (hình 11.1), bề mặt tốc độ 0, phần lớp mỏng nhiệt độ không thay đổi nên tốc độ Tốc độ lớp mỏng (giữa hai khoảng cách ấy) có thay đổi tương đối lớn, dạng cụ thể đường cong phân bố tốc độ phụ thuộc vào trường nhiệt độ tính chất vật lý chất lỏng Khi trao đổi nhiệt mãnh liệt, biến thiên nhiệt độ lớn chuyển động chất lỏng mạnh Chất lỏng chuyển động tự nhiên có hai chế độ chuyển động bản: chảy tầng chảy rối (Hình 11.2), đoạn phía chất lỏng chảy tầng đoạn phía chảy rối Điểm mà chất lỏng chuyển từ chế độ chảy tầng sang chế độ chảy rối gọi điểm tới hạn Khi chất lỏng chảy tầng, nhiệt trở hoàn toàn chiều lớp mỏng dày lớp mỏng chảy tầng định (lớp biên nhiệt) Từ phía bề mặt trao đổi nhiệt ta thấy chiều cao tăng dần, chiều dày lớp mỏng chảy tầng tăng lên, hệ số toả nhiệt cục giảm dần theo chiều cao Khi chất lỏng trạng thái độ chuyển từ chảy tầng sang chảy rối, tác dụng xoáy mãnh liệt lớp chảy rối, hệ số toả nhiệt cục tăng lên Lúc độ cao đạt đến khu vực chảy rối mãnh liệt, hệ số toả nhiệt cục không thay đổi Khi đạt đến độ cao này, hệ số toả nhiệt trung bình tiến tới không đổi Đối với dây kim loại nhỏ, diện tích bề mặt toả nhiệt bé nên nhiệt lượng truyền nhỏ, điều kiện độ chênh lệch nhiệt độ tương đối lớn tồn chế độ chảy tầng Khi độ chênh lệch nhiệt độ nhỏ, xung quanh dây hình thành lớp mỏng không khí gần không chuyển động, chế độ gọi toả nhiệt trạng thái màng 9.1.2 Toả nhiệt đối lưu không gian vô hạn Quá trình toả nhiệt đối lưu tự nhiên có giá trị thực dụng rât lớn, chẳng hạn sinh hoạt hàng ngày: nấu nước, sưởi ấm , kỹ thuật: làm nguội bề mặt lò hơi, lò công nghiệp, làm mát cánh quạt tua bin khí, Hiện có xu hướng sâu vào phạm vi chảy rối cao độ phạm vi tiêu chuẩn Pr bé Chọn hệ toạ độ chiều x trùng với chiều lực nâng Phương trình vi phân toả nhiệt đối lưu tự nhiên có dạng: ∂ω x ∂ω x ∂ω x ∂ ² ω x ∂ ²ω x ∂ ²ω x ρ(ω x + ωy + ωz ) = µ( + + ) + ρgβΔt (9.1) ∂x ∂y ∂z ∂x ² ∂y ² ∂z ² ρ(ω x ∂ω y ∂x + ωy ∂ω y ∂y + ωz ∂ω y ∂z ) = µ( ∂ ²ω y ∂x ² + ∂ ²ω y ∂y ² + ∂ ²ω y ∂z ² ) (9.2) Phương trình chuyển động theo hướng z: Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh Chương 9: Tỏa nhiệt đối lưu Trang 78 ∂t ∂t ∂t ∂²t ∂²t ∂²t ωx + ωy + ωz ) = a( + + ) ∂x ∂y ∂z ∂x ² ∂y ² ∂z ² α=_ Từ (9.1) được: λ ∂t Δt ∂y C ρ C² ω C1 = (9.4) Cµ C ω C ²1 = CρCgCβC∆t (9.5) Từ (9.3) được: C ω C1 =1 Ca (9.6) Từ (9.4) được: C α C1 αl = → = Nu Cλ λ (9.7) C ρ C² ω C1 Thay vào (9.6) ta được: đem: (9.3) Cµ C ω C ²1 = CρCgCβC∆t → = Cµ C ω C ²1 → Cω = Cν ν = → = Pr Ca a Cω = Cv C1 (9.8) (9.9) C g Cβ C Δt C ²1 Cµ Cρ C² ν = CgCβC∆t → gl βΔt = Gr (9.10) C ²1 ν² Theo định lý đồng dạng thứ hai, phương trình tiêu chuẩn trường hợp toả nhiệt đối lưu tự nhiên không gian vô hạn có dạng: Nu = f (Gr, Pr) (9.11) Mikheyev tiến hành thí nghiệm chất lỏng khác bề mặt trao đổi nhiệt khác nhau, kết số liệu sau tổng hợp biểu diễn đồ thị (Hình 11.3) dạng phương trình tiêu chuẩn chỉnh lý sau: thay vào: (9.8) ta được: Nu = C(Gr Pr) nm = C.Ra nm (9.12) đó: αl - tiêu chuẩn Nusselt λ - Ra = (Gr.Pr) gọi tiêu chuẩn Rayleigh - Các thông số vật lý công thức (9.12) chọn theo nhiệt độ tính toán - Nu = tm - Nhiệt độ xác định (tm) chọn nhiệt độ trung bình vách (t w)và môi trường (tf) tm = 0,5(tf + tw) (9.13) - Đối với ống nằm ngang, kích thước xác định đường kính ống: L = D - Đối với ống đứng vách đặt đứng kích thước xác định chiều cao:L = H Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh Chương 9: Tỏa nhiệt đối lưu Trang 79 - Ký hiệu (m) phương trình (9.12) tiêu chuẩn chọn theo nhiệt độ tm Hình 11.3: Toả nhiệt vật khác chuyển động tự nhiên Các đại lượng C, n phương trình (9.12) hệ số thực nghiệm Tùy thuộc vào giá trị (Gr.Pr)m mà C n có giá trị theo bảng (9.1): Bảng 11.1: Trị số C n theo (Gr.Pr)m công thức (11.12) Trạng thái chuyển động (Gr.Pr)m C n < 10 –3 0.5 Quá độ màng sang tầng 10 -3 – 5.10² 1.18 1/8 Chảy tầng 5.10² – 2.107 0.54 1/4 Chảy màng Chảy rối 2.107 – 1013 0.135 1/3 Công thức (9.12) dùng để tính toán toả nhiệt vật có hình dạng khác hình trụ cụt đặt nằm đứng, phẳng hẹp đặt đứng, phẳng đặt nằm ngang, kích thước tính toán dùng kích thước tương đương sau: 1 = + (9.14) L td L ng L dg Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh Chương 9: Tỏa nhiệt đối lưu Trang 80 Hình 9.4: Sự chuyển động tự nhiên chất lỏng bề mặt vật a) Tấm phẳng hẹp bề mặt nóng hướng lên b) Tấm phẳng rộngbề mặt nóng hướng lên c) Tấm phẳng hẹp bề mặt nóng hướng xuống d) Ống ngang Do đó: - Tấm hẹp, dài: chọn cạnh bé làm kích thước xác định - Đối với phẳng đặt nằm ngang, kích thước tính toán chọn theo chiều hẹp, trị số C n chọn theo bảng 9.2: Bảng 9.2 C n Ra Bề mặt nóng Bề mặt nóng hướng hướng lên xuống < 0,001 0,65 0,35 0,001 ÷ 500 1,53 0,83 1/8 500 ÷ x 10 0,70 0,38 1/4 13 x 10 ÷ 10 0,176 0,095 1/3 Ngoài chất lỏng chất khí, toả nhiệt đối lưu tự nhiên quanh ống ngang, dùng công thức sau: Nuf = 0,5Grf1/4Prf1/4(Prf/Prw)1/4 (9.15) Đối với không khí (11.15) có dạng đơn giản hơn: Nuf = 0,46Grf1/4 (9.16) 9.1.3 Toả nhiệt đối lưu tự nhiên không gian hữu hạn Khi đối lưu tự nhiên phát sinh không gian nhỏ, kín khe hẹp chuyển động chất lỏng (khí) phức tạp (hình 9.5) phụ thuộc vào vị rí tương hổ độ chênh lệch nhiệt độ hai bề mặt, hình dạng, kích thước không gian Để đơn giản việc chỉnh lý số liệu thí nghiệm sử dụng thuận tiện, ta dùng kháo niệm “hệ số dẫn nhiệt tương đương” λtđ ; xem trình trao đổi nhiệt phức tạp tượng nhiệt mà dẫn nhiệt qua lớp chất lỏng với hệ số dẫn nhiệt tương đương tính theo công thức sau: Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh Chương 9: Tỏa nhiệt đối lưu Trang 81 Hình9.5: Tính chất chảy vòng tự nhiên chất lỏng không gian hẹp vật kín λ td (9.17) ( t w1 _ t w ) δ Trong điều kiện, lượng nhiệt truyền qua đơn dẫn nhiệt, q= thì: qλ = λ (t _ t ) δ w1 w (9.18) λ td ( t w1 _ t w ) λ q Do đó: (9.19) = δ = td = ε td λ qλ λ _ (t t ) δ w1 w λ Suy ra: ε td = td = Nu ' ; ε td = f (Gr Pr) f (9.20) λ εtđ - hệ số đối lưu, thứ nguyên Trong phương trình (9.20) chọn: - Chiều dày khe hẹp δ làm kích xác định - Nhiệt độ trung bình chất lỏng: tf = 0,5(tw1 + tw2) làm nhiệt độ xác định - Khi (Gr.Pr)f < 103 εtđ = - Khi 103 < (Gr.Pr)f < 106 ε td = 0,105(Gr Pr) 0f,3 (9.21) - Khi 106 < (Gr.Pr)f < 1010 ε td = 0,40(Gr Pr) 0f,2 (9.22) - Theo Mikheyev (Gr.Pr)f > 103 thay công thức (9.21) (9.22) bằng: Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh Chương 9: Tỏa nhiệt đối lưu Trang 82 0, 25 ε td = 0,18(Gr Pr) f (9.23) 9.2 TỎA NHIỆT ĐỐI LƯU KHI CHẤT LỎNG CHUYỂN ĐỘNG CƯỠNG BỨC 9.2.1 Đặc trưng chuyển động ống Chảy tầng chảy rối Tính chất chuyển động chất lỏng ống rãnh phân làm hai chế độ bản: chảy tầng chảy rối Khi tốc độ nhỏ, chất lỏng chảy tầng tốc độ lớn chất lỏng chảy rối, tốc độ của điểm phân giới chuyển chuyển từ chảy tầng sang chảy rối gọi tốc độ tới hạn Tiêu chuẩn Reynolds định chế độ chuyển động chất lỏng Căn vào tính chất này, dùng trị số Re để phân biệt điểm tới hạn, trị số Re tương ứng có chế độ chảy tầng chuyển sang chảy rối gọi số Reynolds tới hạn (Re tới hạn = 2200) - Khi Re < 2200 chất lỏng chảy tầng - Khi Re > 2200 chất lỏng chảy rối Sự phân bố tốc độ ống tiêu chuẩn Reynolds định Ở trạng thái chảy tầng (Hình 9.6a) phân bố tốc độ tiết diện ngang có dạng đường parabol, trạng thái chảy rối phân bố tốc độ tiết diện ngang có dạng đường parabol cụt (hình 9.6b) Hình9.6: a) Sự phân bố tốc độ chất lỏng chảy tầng, Re < Rettb b) Sự phân bố tốc độ chất lỏng chảy rối, Re > Rettb Trong tính toán thực tế thường cần biết tốc độ trung bình mặt cắt đủ Tốc độ trung bình tỷ số lưu lượng thể tích diện tích mặt cắt mà chất lỏng chuyển động qua V (9.24) ω = ∫ ωdf = f f ω - tốc độ trung bình chất lỏng mặt cắt, [m/s] đó: V – lưu lượng thể tích chất lỏng qua mặt cắt, [m3/s] f - diện tích mặt cắt, [m2] Tính chất chuyển động chất lỏng mặt cắt chỗ vào Ở chỗ vào có tiếp xúc đột ngột chất lỏng với bề mặt nên quy luật phân bố tốc độ có thay đổi lớn Do phải thời gian chuyển động ống phân bố tốc độ đạt đến ổn định Xét chất lỏng chuyển động từ không gian lớn vào ống tròn Tại mặt cắt chỗ vào, tốc độ không thay đổi, sau vào ống tác dụng ma sát chất lỏng Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh Chương 9: Tỏa nhiệt đối lưu Trang 83 bề mặt, tốc độ chất lỏng gần sát vách giảm dần xuống, lưu lượng chất không thay đổi, tốc độ phần mặt cắt tăng dần lên (Hình 9.7) Hình: 9.7: Sự phân bố tốc độ chất lỏng chuyển động từ không gian rộng vào ống Qua giai đoạn định, chuyển động chất lỏng đạt đến trạng thái ổn định, lớp biên hợp chung thành (Hình 9.8) Trong trường hợp tốc độ chất lỏng lớn, không cần đạt đến trạng thái ổn định, lớp biên trở thành chảy rối Trong lớp biên chảy rối sát vách có lớp chảy tầng gọi lớp đệm tầng, độ dày δl khu vực ổn định hàm số Reynolds (đường chấm chấm Hình 9.9), cuối lớp biên chảy rối hợp lại thành một, trạng thái chuyển động ổn định chảy rối Trên quãng đường chuyển động tiến tới ổn định chiều dày lớp biên chỗ vào không sau tăng dần lên, đạt đến trạng thái ổn định chiều dày lớp biên bán kính ống Chiều dày lớp đệm tầng tính theo công thức: d δ l = 64,2 / (9.25) Re trị số tuyệt đối, chiều dày lớp mỏng giảm dần Re tăng: Re = 104 Re = 105 Re = 106 δl δl δl 1 = = = d 466 d 3600 d 23400 Bởi đoạn ống chuyển động tiến tới ổn định, chiều dày lớp biên có thay đổi lớn, chế độ toả nhiệt thay đổi theo Tính chất trường nhiệt độ mặt cắt chỗ vào Tương tự lớp biên thuỷ động, chỗ vào lớp biên nhiệt bắt đầu hình thành dày dần lên Ở gần chỗ vào chiều dày lớp biên nhiệt tương đối mỏng gradien nhiệt độ lớp biên lớn, xa chỗ vào chiều dày lớp biên tăng lên gradien nhiệt độ lớp biên giảm xuống ổn định Tương ứng với thay đổi này, hệ số toả nhiệt lớn gần chỗ vào sau giảm dần cuối tiến tới trị số giới hạn xac định, gọi đạt đến “trạng thái ổn định nhiệt” Ảnh hưởng toả nhiệt đối lưu tự nhiên đến phân bố tốc độ Khi xảy trình trao đổi nhiệt, nhiệt độ mặt cắt không đồng đều, biến thiên nhiệt độ ảnh hưởng đến độ nhớt chất lỏng làm cho phân bố tốc độ mặt cắt khác với chuyển động đẳng nhiệt Tuỳ theo phương hướng dòng nhiệt mà thay đổi khác (hình 9.8) Thực nghiệm dùng hệ số (Prf/Prw)0,25 để hiệu chỉnh đến phương hướng dòng nhiệt Ngoài ra, nhiệt độ mặt cắt không đồng tạo nên đối lưu tự nhiên phần, ảnh hưởng đến phân bố tốc độ Những điều kiện sau làm ảnh hưởng đối lưu tự nhiên phân bố tốc độ: Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh Chương 9: Tỏa nhiệt đối lưu Trang 84 a) b) Hình 9.8: Sự phân bố tốc Hình 9.9: Sự phân bố nhiệt độ tốc độ độ chất lỏng mặt cắt chất lỏng bị làm nguội - Chảy đẳng nhiệt a) Không xét đến ảnh hưởng chuyển - Chảy chất lỏng bị làm nguội động tự nhiên - Chảy chất lỏng gia nhiệt b) Có xét đến ảnh hưởng chuyển động tự nhiên - Cách đặt ống (đứng hay nằm) - Chuyển động tự nhiên hay chuyển động cưỡng chiều hay ngược chiều Khi xảy chuyển động tự nhiên ống đứng mà sát vách có chiều trùng với hướng chuyển động cưỡng bức, phân bố tốc độ nhiệt độ có dạng Hình 11.9a chuyển dạng hình 9.9b 9.2.2 Tỏa nhiệt chảy rối (Re > 104) Khi chảy rối, truyền lượnggiữa phần tử chất lỏng hỗn hợp thân chất lỏng Trong chế độ chảy rối mãnh liệt Re > 10 4, trình xáo trộn mạnh, làm cho nhiệt độ lõi chảy rối không thay đổi, có lớp đệm tầng xuất sư biến thiên nhiệt độ cách rõ ràng Trên mặt cắt ngang, phân bố nhiệt độ không gây nên chuyển động tự nhiên ảnh hưởng chuyển động cưỡng chất lỏng bỏ qua Phương trình tiêu chuẩn tìm trường hợp có dạng: Nu = f(Re, Pr) (9.26) Xét đến việc chọn nhiệt độ xác định loại trừ ảnh hưởng phương hướng dòng nhiệt, phương trình hiệu chỉnh thành: Pr Nu f = f (Re f Prf )( f ) 0, 25 (9.27) Prw M Mikheyev tìm phương trình tiêu chuẩn sau: 0,80 0, 43 Pr Nu f = 0,021(Re f Prf )( f ) 0, 25 ε l ε R Prw (9.28) đó: Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh Chương 9: Tỏa nhiệt đối lưu Trang 85 - Nhiệt độ xác định: nhiệt độ trung bình chất lỏng tf - Prf – tiêu chuẩn Pr chọn theo nhiệt độ trung bình chất lỏng tf - Prw – tiêu chuẩn Pr chọn theo nhiệt độ vách ống tw - Kích thước xác định tính theo đường kính tương đương: d td = = 4F U - F- diện tích tiết diện ngang dòng chất lỏng [m²] - U – chu vi ướt (chu vi mà chất lỏng tiếp xúc với bề mặt trao đổi nhiệt) [m] - εl - trị số hiệu chỉnh ảnh hưởng đoạn đầu ống, phụ thuộc vào trị số Re l/d (εl = f(Re, l/d) - l, d chiều dài đường kính ống [m] + Nếu l/d > 50 εl = + Nếu l/d < 50 εl tra bảng 9.3 Bảng 9.3 εl = f(Ref, l/d) l/d Ref 10 15 20 30 40 50 1÷10 1,65 1,50 1,34 1,23 1,17 1,13 1,07 1,03 2÷10 1,51 1,40 1,27 1,18 1,13 1,10 1,05 1,02 5÷10 1,34 1,27 1,18 1,13 1,10 1,08 1,04 1,02 1÷10 1,28 1,22 1,15 1,10 1,08 1,06 1,03 1,02 1÷106 1,14 1,11 1,08 1,05 1,04 1,03 1,02 1,01 * Ảnh hưởng ống cong: Khi chất lỏng chuyển động ống cong, tác động lực ly tâm gây nên chảy vòng phụ mặt cắt Bán kính cong R lớn tác dụng lực ly tâm nhỏ, R = ∞ ống thẳng Thực nghiệm với nước không khí lập nên quan hệ sau: d (9.29) ε R = + 1,77 R R- bán kính cong ống xoắn, [m] d - đường kính ống, [m] Đối vói không khí Pr ≈ const, dùng công thức đơn giản sau: 0,80 Nu f = 0,018 Re f (11.30) 9.2.3 Tỏa nhiệt chảy tầng (Re < 2200) Khi trị số Re < 2200 chất lỏng chảy tầng Khi chảy tầng, theo hướng bán kính chất lỏng tác dụng hỗn hợp nhiệt trở lớn khiến cho hệ số toả nhiệt có trị số nhỏ chảy rối nhiều Đa số thiết bị nhiệt có chất lỏng chuyển động cưỡng luôn phạm vi chảy rối, trường hợp gặp phải chế độ chảy tầng Khi phân tích phương trình vi phân trao đổi nhiệt, ta phương trình vi phân chuẩn có dạng sau: Nu = f(Re, Pr, Gr) (9.31) Xét đến việc chọn nhiệt độ xác định ảnh hưởng phương hướng dòng nhiệt, phương trình chuyển thành: Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh Chương 9: Tỏa nhiệt đối lưu Trang 86 Nu f = f (Re f , Grf , Prf )( Prf 0, 25 ) Prw (9.32) Từ kết thực nghiệm với nhiều chất lỏng khác ống có kích thước hình dạng khác nhau, người ta hiệu chỉnh phương trình tiêu chuẩn sau: Nu f = 0.15 Re 0f.33 Prf0.43 Grf0.1 ( Prf 0,25 ) ε l Prw (9.33) đó: - Nhiệt độ xác định: nhiệt độ trung bình chất lỏng tf - Prf – tiêu chuẩn Pr chọn theo nhiệt độ trung bình chất lỏng tf - Prw – tiêu chuẩn Pr chọn theo nhiệt độ vách ống tw - Kích thước xác định tính theo đường kính d ; đối vơí rãnh có tiết diện ngang hình tròn dùng đường kính tương đương làm kích thước 4F xác định: d td = = U - F- diện tích tiết diện ngang dòng chất lỏng [m²] - U – chu vi ướt (chu vi mà chất lỏng tiếp xúc với bề mặt trao đổi nhiệt) [m] - εl - trị số hiệu chỉnh ảnh hưởng đoạn đầu ống, phụ thuộc vào trị số Re l/d (εl = f(Re, l/d) - l, d chiều dài đường kính ống [m] + Nếu l/d > 50 εl = + Nếu l/d < 50 εl tra bảng 9.4 Bảng 9.4: εl = f(l/d) l/d 10 15 20 30 40 ≥ 50 εl 1.9 1.7 1.44 1.28 1.18 1.13 1.05 1.02 Trong trường hợp toả nhiệt loại dầu có độ nhớt cao chuyển động ống, tác dụng đối lưu tự nhiên bị ức chế gần chảy tầng hoàn toàn (không bị nhiẽu loạn), lúc tiêu chuẩn Gr (9.33) bỏ qua Đối với không khí phương trình tiêu chuẩn có dạng đơn giản sau: Nu f = 0,13 Re 0f ,33 Gr f0,1ε l (9.34) 9.2.4 Tỏa nhiệt trạng thái độ (2200 < Re < 104 ) Khi trị số 2200 < Re < 10 chất lỏng chảy độ Khi phương trình tiêu chuẩn có dạng: Pr Nu = K ο Pr 0, 43 ( f ).ε l (9.35) Prw trị số K0 = f(Ref) tra bảng 11.5 Bảng 9.5 K0 = f(Ref ) Ref 103 2,2 2,3 2,5 3,0 3,5 K0 1,9 3,2 4,0 6,8 9,5 Trường ĐHCN Tp HCM 4,0 11 5,0 16 6,0 19 7,0 24 8,0 27 9,0 30 10 33 Khoa CN Nhiệt - Lạnh Chương 9: Tỏa nhiệt đối lưu Trang 87 εl - chọn theo trường hợp chảy tầng Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh ... 0,105(Gr Pr) 0f,3 (9 .2 1) - Khi 106 < (Gr.Pr)f < 1010 ε td = 0,40(Gr Pr) 0f,2 (9 .2 2) - Theo Mikheyev (Gr.Pr)f > 103 thay công thức (9 .2 1) (9 .2 2) bằng: Trường ĐHCN Tp HCM Khoa CN Nhiệt - Lạnh Chương. .. thì: qλ = λ (t _ t ) δ w1 w (9 .1 8) λ td ( t w1 _ t w ) λ q Do đó: (9 .1 9) = δ = td = ε td λ qλ λ _ (t t ) δ w1 w λ Suy ra: ε td = td = Nu ' ; ε td = f (Gr Pr) f (9 .2 0) λ εtđ - hệ số đối lưu, thứ.. .Chương 9: Tỏa nhiệt đối lưu Trang 78 ∂t ∂t ∂t ∂²t ∂²t ∂²t ωx + ωy + ωz ) = a( + + ) ∂x ∂y ∂z ∂x ² ∂y ² ∂z ² α=_ Từ (9 . 1) được: λ ∂t Δt ∂y C ρ C² ω C1 = (9 . 4) Cµ C ω C ²1 = CρCgCβC∆t (9 . 5) Từ