Xem mẫu

  1. CHƯƠNG 3 TRAO ĐỔI NHIỆT ĐỐI LƯU 3.1. CÔNG THỨC NEWTON VÀ HỆ SỐ TỎA NHIỆT α 3.1.1. Định nghĩa và phân loại tỏa nhiệt Trao đổi nhiệt đối lưu hay toả nhiệt, là hiện tượng trao đổi nhiệt giữa một lớp chất lưu (chất lỏng hoặc chất khí) chuyển động với bề mặt của 1 vách tiếp xúc. Hiện tượng chất lưu chuyển động theo dòng tuần hoàn gọi là đối lưu. Theo nguyên nhân gây chuyển động đối lưu, người ta phân toả nhiệt ra 2 loại: 1) Toả nhiệt tự nhiên là hiện tượng toả nhiệt khi đối lưu được sinh ra một cách tự nhiên, do trọng lực tạo ra trong chất lưu không đồng nhất về nhiệt độ. 2) Toả nhiệt cưỡng bức là hiện tượng toả nhiệt khi đối lưu được tạo ra do một ngoại lực, ví dụ do bơm quạt. 3.1.2. Công thức Newton và hệ số toả nhiệt α Xét mô hình toả nhiệt gồm một bề mặt vách t nhiệt độ tW tiếp xúc một lớp chất lưu chuyển động có nhiệt độ giảm dần từ tW sát vách đến tf ở xa tw ω vách, như hình Hình 3.1. Lượng nhiệt toả ra chất 1 α 2 lưu qua mặt vách được tính theo công thức qui ước, F gọi là công thức Newton, có dạng: t f q=α(tW - tf), [W/m2] và Q = αF(tW - tf), [W] x trong đó tW là nhiệt độ mặt vách, tf là nhiệt độ chất Hình 3.1 Phân bố t(x) trong lưu ở xa vách, F là diện tích mặt tiếp xúc và chất lưu khi tỏa nhiệt q Q α= = , [W/m2K] gọi là hệ số toả nhiệt . ∆t F∆t Hệ số toả nhiệt α đặc trưng cho cường độ toả nhiệt, là ẩn số chính của mọi bài toán toả nhiệt. 28
  2. 3.2. PHƯƠNG TRÌNH TỔNGG QUÁT CỦA HỆ SỐ TOẢ NHIỆT 3.2.1. các thông số ảnh hưởng đến α Hệ số tỏa nhiệt α phụ thuộc vào các yếu tố gây ra đối lưu và dẫn nhiệt trong chất lưu, bao gồm các thông số chính sau đây: 1) Các thông số hình học của mặt toả nhiệt . Hình dạng, kích thước và vị trí của mặt toả nhiệt ảnh hưởng tới dòng chuyển động của chất lưu, do đó sẽ ảnh hưởng tới α. Đặc trưng hình học của bề mặt toả nhiệt có thể gồm nhiều kích thước khác nhau, nhưng khí tính α thường chọn một kích thước duy nhất, gọi là kích thước định tính lt. Kích thước định tính lt do người lập công thức tính α lựa chọn, theo qui tắc như sau ⎧chiãöu cao h cuía màût phàóng hoàûc truû thàóng âæïng lt = ⎨ ⎩âæåìng kênh tæång âæång cuía äúng nàòm ngang hoàûc bãö màût ngoaìi váût hæîu haûn Ống nằm ngang với diện tích và chu vi của tiết diện chứa chất lỏng là f và u, 4f sẽ có đường kính tương đương d td = . Vật hữu hạn với thể tích V, diện tích xung u quanh F sẽ có dtđ = GV , [m]. F ddt=4f/u ddt=6V/F u f F V Hình 3.2 dtd của vách phẳng và V hình trụ 2) Các thông số vật lý của chất lưu Các thông số vật lý trực tiếp ảnh hưởng đến α bao gồm: - Các thông số ảnh hưởng đến chuyển động là: khối lượng riêng ρ[kg/m3], hệ ∆V µ số nở nhiệt β = , [K −1 ] , độ nhớp động học ν = [m2/s] V∆T ρ 29
  3. - Các thông số ảnh hưởng tới dẫn nhiệt là: hệ số dẫn nhiệt λ[W/mK], hệ số λ khuếch tán nhiệt a = [m2/s]. ρC p Các thông số vật lý nới trên đều thay đổi theo nhiệt độ chất lưu . Để xác định giá trị các thông số vật lý khi tính α, người lập công thức qui định 1 trị số nào đó của nhiệt độ chất lưu, gọi là nhiệt độ định tính. Nhiệt độ định tính [t] có thể lấy một trong các giá trị sau: ⎧ ⎪t f laì nhiãût âäü cháút læu xa vaïch ⎪ [ t ] = ⎨t W laì nhiãût âäü cháút læu saït vaïch ⎪ 1 ⎪t m = ( t W + t f ) laì nhiãût âäü cháút læu åí giæîa låïp biãn nhiãût ⎩ 2 3) Các thông số đặc trưng cho cường độ đối lưu - Đối lưu tự nhiên là dòng đối lưu tự phát sinh trong chất lưu khi có độ chênh trọng lượng riêng giữa các lớp chất lưu. độ chênh trọng lượng riêng tỷ lệ thuận với gia tốc trọng lực g[m/s2], với hệ số nở nhiệt β[K-1] và với độ chênh nhiệt độ ∆t = tW = tf giưac chất lưu ở gần và ở xa vách. Do đó, cường độ đối lưu tự nhiên đặc trưng bằng tích số gβ∆t, [m/s] của chất lưu. -Đối lưu cưỡng bức gây ra bởi lực cưỡng bức của bơm, quạt được đặc trưng bằng tốc độ ω(m/s) của chất lưu. Khi đối lưu cững bức, nếu trong chất lưu có g≠0 và ∆t ≠ 0 thì luôn kèm theo đối lưu tự nhiên. 4) Các thông số đặc trư ng cho chế độ chuyển động của chất lưu . Khi chảy tầng, các phân tử chất lưu chuyển động song song vách nên hệ số α không cao. Khi tăng vận tốc ω đủ lớn dòng chảy rối sẽ xuất hiện. Lúc này các phân tử chất lưu xuất hiện các thành phần chuyển động rối loạn theo phương ngang, tăng cơ hội va đập lên vách, khiến cho hệ số α tăng cao. Chế độ chuyển động chất lưu được phân ra 3 dòng, đặc trưng bởi các thông số ω, l, ν , thông qua giá trị của vận tốc không thứ nguyên hay số Reynolds như sau: 30
  4. ⎧Re ≤ 2300 doìng chaíy táöng ωl ⎪ Re = = ⎨Re ∈ (2300 ÷ 10 4 ) - doìng quaï âäü ν ⎪ ⎩Re ≥ 10 - doìng chaíy räúi 4 3.2.2. phương trình tổng quát của hệ số tỏa nhiệt Phương trình tổng quát của α là phương trình chứa tất cả các thông số ảnh hưởng tới giá trị α, như đã phân tích trên đây, có dạng: α = f(l, ρ, ν, a, λ, g, β(tW- tf), ω), (1) Đây là phương trình tính α dạng tích phân tổng quát, phụ thuộc vào 10 biến số, mà dạng cụ thể của nó sẽ được tìm chủ yếu bằng thực nghiệm. 3.3. HỆ PHƯƠNG TRÌNH VI PHÂN TỎA NHIỆT 3.3.1.các phương trình cân bằng nhiệt–động lực học chất lưu t Quá trình trao đổi nhiệt đối lưu của phân tố chất lưu dV có các thông số ρ, Cp, a, µ, qv, p, t, ω, với mặt tw µρ tiếp xúc W có thể mô tả bằng 1 hệ các phương trình, ω a g p qv gồm phương trình cân bằng nhiệt, phương trình cân W qλ dV qω t bằng động lực học và phương trình cân bằng lượng f chất lưu như sau. 1) Phương trình cân bằng nhiệt Hình 3.3 Để lập hệ phương x Định luật bảo toàn năng lượng cho dV có nội dung trinh trình vi phân tỏa nhiệt là: độ tăng entanpi của dV = hiệu số dòng nhiệt(vào – ra )dV + lượng nhiệt tự phát sinh trong dV hay: ∂t ρdVC p = −divqdV + q v dV với ∂τ divq = div(ρωC p t − λgr adt ) = ρC p (ω.gr adt + tdivω) − λ∇ 2 t ở đây ω.gr adt là tích vô hướng của 2 vectow ω và gr adt λ Do đó, nếu ký hiệu a = và ρC 31
  5. ∂t ∂t ∂x ∂t ∂y ∂t ∂z ∂t dt + ω.gr adt = + + + = , [K/s] là đạo hàm toàn phần của nhiệt ∂τ ∂τ ∂τ ∂x ∂τ ∂y ∂τ ∂z dτ độ theo thời gian thì phương trình cân bằng nhiệt có dạng: dt q = a∇ 2 t − tdivω + v dτ ρC p ∂t q Nếu dV không chuyển động, ω = 0, thì có = a∇ 2 t + v là phương trình vi phân ∂τ ρC p dẫn nhiệt như nêu ở chương 2 2) Phương trình can bằng động lực họ cho dV có nội dung là: Lực quán tính của dV = trọng lực của dV + hiệu số áp lực lên dV + lực ma sát quanh dV, hay phương trình có dạng: dω ρ = ρg − gr adp + µ∇ 2 ω với µ∇ 2 ω = µ( i∇ 2 ω x + j ∇ 2 ω y + k ∇ 2 ω z ) [N/m3] là dτ tổng các lực nội ma sát ứng với 1m3 của dV. Phương trình trên còn được gọi là phương trình Naver – Stockes, là phương trình cơ bản của động lực học chất lưu 3) Phương trình liên tục Khi trong phân tố chất lưu dV không có điển rò hoặc điểm nguồn, thì hiệu số lưu luợng (vào - ra)dV bằng độ tăng khối lượng riêng của dV, hay: ∂ρ = −div(ρω) = − gr adρ.ω − ρdivω ∂τ ∂ρ dρ Suy ra + ω.gr adρ = =- ρdivω ∂τ dτ Với chất lỏng không chịu nén như nước hoặc dầu, ρ = const, thì phương trình liên tục có dạng divω = 0 3.3.2. Mô hình toán học của bài toán tỏa nhiệt Phát biểu toán học của bài toán toả nhiệt là tìm các hàm phân bố vận tốc ω , nhiệt độ t và hệ số toả nhiệt α thoả mãn hệ phương trình vi phân, gồm 3 phương trình cân bằng nhiệt - động lực học chất lưu và 4 loại điều kiện đơn trị như sau: 32
  6. ⎧ dt ⎪ dτ = a∇ t − tdivω + q v 2 ⎪ ⎪ ρd ω = pg − gr adp + µ∇ 2 ω ⎪ dτ ⎪ ∂ρ ⎪ = −div(ρω), våïi caïc ÂKÂT : ⎨ ∂τ ⎪− hçnh hoüc : hçnh daûng, kêch thæåïc l cuía W ⎪ ⎪- váût lyï : luáût xaïc âënh (ρ, C p , µ, a, λ, q v ) = f(t) ⎪- ban âáöu : t(∀M, τ = 0) = t(x, y, z) ⎪ ⎪- biãn W3 : - λgradt(M ∈ W3 ) = α[t(M ∈ W3 ) - t f ] ⎩ Đây là bài toán rất phức tạp, hiện nay chưa có lời giải tổng quát. Việc tính α chủ yếu dựa vào các số liệu và công thức thực nghiệm, như sẽ trình bày ở bài sau. 3.4. LÝ THUYẾT ĐỒNG DẠNG Lý thuyết đồng dạng là môn học nghiên cứu tính chất của các hiện tượng vật lý đồng dạng nhau, là một lý thuyết chỉ đạo cho công tác thực nghiệm. 3.4.1. Các khái niệm đồng dạng Các khái niệm đồng dạng được mở rộng từ đồng dạng hình học. 1) Đồng dạng hình học: Hai vật thể hình học được gọi là đồng dạng nếu tỷ sô 2 kích thước l z z' tương ứng bất kỳ là không đổi. l' y y' x' x Ví dụ: Hai hình hộp sẽ đồng dạng khi Hình 3.4 Đồng dạng x ' y' z ' l ' hình học = = = = C l = const. x y z l Khi 2 vật thể đồng dạng với tỷ số c l thì tỷ số các diện tích hoặc thể tích tương f' V' ứng là = Cl 2 = C3 l f V 2) Đồng dạng của 2 trường vật lý 33
  7. Hai trường của một đại lượng vật lý ϕ = ϕ(M( x, y, z), τ) được gọi là đồng dạng nhau nếu tỷ số 2 giá trị của ϕ tại 2 điểm tương ứng bất kỳ trong không gian, thời gian là ϕ' ϕ' ( M ' (C l x , C l y, C l z), C τ τ) không đổi: = = C ϕ = const , ∀(M, τ) ϕ ϕ(M ( x , y, z), τ) Để 2 trường vật lý đồng dạng, cần có 3 hằng số C l , C τ , C ϕ . 3) Đồng dạng 2 hiện tượng vật lý: Hai hiện tượng vật lý cùng được mô tả bởi phương trình F(ϕ1, ϕ1,…ϕi,…ϕn)=0 được gọi là đồng dạng nhau, nếu 2 trường của mỗi đại lượng vật lý cùng tên đồng dạng nhau: ϕi' ϕ' (M ' (C l x , C l y, C l z), C l τ) = =Cφi, ∀( M, τ) và ∀i = (1 ÷n). ϕi ϕ(M ( x , y, z), τ) Hai hiện tượng vật lý chỉ đồng dạng trong không gian đồng dạng về hình học, và cần (n+2) hằng số đồng dạng, là Cl , C τ và n là hằng số Cφi 3.4.2. Tiêu chuẩn đồng dạng và phương trình tiêu chuẩn. 1) Định nghĩa: Tiêu chuẩn đồng dạng là 1 biến số không thứ nguyên chứa một số đại lượng vật lý, được suy ra từ tính đồng nhất(hay tương đương) của phương trình mô tả 2 hiện tượng đồng dạng. 2) Ví dụ: khi 2 hiên tượng toả nhiệt đồng dạng, thì phương trình ∂t ∂t α∆t = −λ và α' ∆t ' = −λ' là đương nhau → thay các thông số ϕ i' = C ϕi ϕ i vào ∂n ∂n ' phương trình thứ 2 có C α αC t ∆t = −C λ λC t ∂t /(C l ∂n ) →Do phương trình này tương đương với phương trình đầu , nên có ⎛ α' ⎞⎛ l' ⎞ ⎛ λ' ⎞ C α C t = C λ C t / C l → C α C l / C λ = 1 = ⎜ ⎟⎜ ⎟ ⎜ ⎟ → Suy ra ⎝ α ⎠⎝ l ⎠ ⎝ λ ⎠ α ' l ' αl αl = = idem (như nhau) →tổ hợp không thứ nguyên Nu = có giá trị như nhau λ' λ λ cho 2 hiện tượng đồng dạng, được gọi là tiêu chuẩn Nusselt. Xét tương tự với các phương trình khác trong hệ phương trình vi phân toả nhiệt, có thể dẫn ra các tiêu chuẩn : 34
  8. ν ωl Tiêu chuẩn Prandtl Pr = , tiêu chuẩn Reynolds Re = và tiêu chuẩn a ν gl 3 ρ − ρ0 Galilei: Ga = . Nếu nhân Ga với = β ∆t sẽ được tiêu chuẩn Grashoff: ν 2 ρ0 gβ ∆tl 3 Gr = . ν2 3)Phân loại : các tiêu chuẩn đồng dạng được phân ra 2 loại : tiêu chuẩn xác định chỉ chứa các thông số đã cho trong điều kiện đơn trị (ví dụ Pr và Gr) và tiêu αl chuẩn chưa xác định có chứa 1 thông số chưa biết. Ví dụ: Nu = là chưa xác λ ωl định vì chứa α chưa biết. Re = là tiêu chuẩn xác định khi toả nhiệt cưỡng bức ν vì đã biết ω, và là tiêu chuẩn chưa xác định trong toả nhiệt tự nhiên, vì lức này chưa biết ω. 4) Phương trình tiêu chuẩn là phương trình liên hệ các tiêu chuẩn đồng dạng. Ví dụ, phương trình tiêu chuẩn suy ra từ hệ phương trinh vi phân toả nhiệt có dạng tổng quát Nu=f(Re, Gr, Pr). 3.4.3. Các định lý đồng dạng và ứng dụng của nó. 1) Định lý 1: hai hiện tượng đồng dạng thì trị số các tiêu chuẩn đồng dạng cùng tên bằng nhau. Định lý này cho phép xác định các thông số trong hiện tượng thực tế bằng cách đo các thông sô trong mô hình thí nghiệm đồng dạng với nó. 2) Định lý 2: Hệ phương trình vi phân và phương trình tiêu chuẩn suy ra từ hệ ấy là tương đương nhau. Định lý này chỉ ra cách giải 1 hệ phưông trình vi phân bằng thực nghiệm theo các bước sau: - Từ hệ phương trình vi phân suy ra các tiêu chuẩn đồng dạng. - Lập mô hình đồng dạng đo các đại lý có trong các tiêu chuẩn. - Tính và lập bảng tương ứng các tiêu chuẩn đồng dạng và xấp xỉ bảng này bằng một phương trình tiêu chuẩn. 35
  9. Phương trình tiêu chuẩn là một công thức thực nghiệm, tương đương với nghiệm của hệ phương trình vi phân đã cho. 3) Định lý 3: Điều kiện cần và đủ để 2 hiện tượng đồng dạng là các điều kiện đơn trị đồng dạng và các tiêu chuẩn xác định cùng tên bằng nhau. Định lý này chỉ ra cách lập mô hình thí nghiệm đồng dạng với hiện tượng cần khảo sát, bằng cách tạo các trường thông số cho trước đồng dạng và chọn các hệ số đồng dạng sao cho tiêu chuẩn xác định bằng nhau. Nó cũng cho hay phạm vi ứng dụng của một công thức thực nghiệm là lớp các hiện tượng đồng dạng với mô hình đã xét. 3.5. PHƯƠNG PHÁP PHÂN TÍCH THỨ NGUYÊN TÌM PHƯƠNG TRÌNH TIÊU CHUẨN. 3.5.1. Cơ sở và các bước áp dụng phương pháp PTTN - Phương pháp PTTN là phương pháp tìm các biến số độc lập không thứ nguyên ảnh hưởng tới 1 hiện tượng vật lý nào đó, thay vào quan hệ giữa các đại lượng có thứ nguyên đã biết. - Cơ sở của phương pháp PTTN là nguyên lý cho rằng, nội dung của phương trình mô tả một hiện tượng vật lý sẽ không đổi khi thay đổi đơn vị đo các đại lượng vật lý chứa trong phương trình. - Các bước áp dụng phương pháp PTTN gồm : 1) Phân tích thứ nguyên các đại lượng vật lý có ảnh hưởng tới hiện tượng để tìm hệ đơn vị cơ bản. 2) Thay đổi hệ đơn vị cơ bản theo các tỷ lệ thích hợp để khử các biến số phụ thuộc và tìm ra các biến độc lập không thứ nguyên. - Đặc điển của phương pháp PTTN là chỉ dựa vào nhận xét tổng quát về các đại lượng vật lý có ảnh hưởng tới hiện tượng, không cần lập hệ phương trình vi phân mô tả hiện tượng. Do đó phương pháp này sẽ tiện lợi khi nghiên cứu các hiện tượng mới. 3.5.2. Tìm phương trình tiêu chuẩn tỏa nhiệt bằng phương pháp PTTN 36
  10. Nhận xét tổng quát về các đại lượng ảnh hưởng tới α trong bài (3.2) đã nêu ra phương trình tổng quát của toả nhiệt là: α = f (l, ρ, ν, a , λ, gβ∆t , ω) (1) Phương pháp PTTN chuyển phương trình (1) về dạng tiêu chuẩn gồm 2 bước như sau: 1) Phân tích thức nguyên các đại lượng vật lý trong phương trình (1) gồm: [l] = [m], [ρ] = [kg / m 3 ], [ν] = [m 2 s −1 ] = [a ], [λ ] = [ Wm −1 K −1 ] = [kgms −3 K −1 ] , do W= Js-1 = Nms-1=kgm2s-3 , [gβ∆t]=[ms-2], [ω]=[ms-1] và [α] = [Wm-2K-1]=[kgs-3K-1]. Suy ra đơn vị cơ bản của phương trình (1) là hệ gồm 4 đơn vị sau [kg, m, s, K] 2) Dùng hệ đơn vị cơ bản mới, là [Gkg, Mm, Ss, DK] với G, M, S, D là các hằng số tỷ lệ khác không tùy ý, thì phương trình (1) có dạng (1’) như sau: GS-3D-1α = f(M l , GM-3ρ, M2S-1ν, M2S-1a, GMS-3D-1λ, MS-2gβ∆t, MS-1ω) Để khử các đại lượng phụ thuộc, cần chọn 4 hằng số G, M, S, D theo 4 điều kiện sau: Ml = 1 ⎫ ⎧M = l −1 ⎪ ⎪ GM −3 ρ = 1 ⎪ −3 −1 ⎪G = l ρ ⎬ → suy ra ⎨ M 2 S −1ν = 1 ⎪ ⎪S = νl −2 GMS −3 D −1λ = 1⎪ ⎭ ⎪D = λl 2 ρ −1ν −3 ⎩ Thay các giá trị của M, G, S, D vào (1’) sẽ có: α l /λ = f(1; 1’;1;a/ν;1; gβ∆t l 3/ν2, ω l /ν), (1”) Đây là phương trình tiêu chuẩn tổng quát, có dạng Nu = (Pr, Gr, Re). Theo nguyên lý nói trên, các phương trình (1’), (1”) và (1) là tương đương nhau. 3.5.3. Ý nghĩa của tiêu chuẩn đồng dạng Các tiêu chuẩn đồng dạng là các biến số độc lập không thứ nguyên(KTN), đặc trưng cho hiện tượng toả nhiệt. 37
  11. αl 1) Nu = là tiêu chuẩn Nusselt hay hệ số toả nhiệt KTN, đặc trưng cho cường λ độ toả nhiệt, là tiêu chuẩn chưa xác định. 2) Pr = ν/a là tiêu chuẩn Prandtl hay độ nhớt KTN đặc trưng cho tính chất vật lý của chất lưu, là tiêu chuẩn xác định, thường được coi như một thông số vật lý tổng hợp. 3) Gr= gβ (tw - tf ) l 3/ν2 là tiêu chuẩn Grashoff hay lực đẩy Archimede KTN, đặc trưng cho cường độ đối lưu tự nhiên, là tiêu chuẩn xác định. 4) Re = ω l /ν là tiêu chuẩn Reynolds hay vận tốc không thứ nguyên, đăc trưng cho cường độ đối lưu cưỡng bức và chế độ chuyển động của chất lưu, là tiêu chuẩn xác định khi đối lưu cưỡng bức , và là tiêu chuẩn chưa xác định phụ thuộc Gr khi đối lưu tự nhiên. 3.5.4. Các dạng đăc biệt của phương trình tiêu chuẩn Phương trình tiêu chuẩn tổng quát Nu = f(Re, Gr, Pr) sẽ có dạng đặc biệt trong các trường hợp sau đây. 1) khi đối lưu tự nhiên hoàn toàn, thì Re là ẩn số phụ thuộc Gr, nên phương trình tiêu chuẩn sau khi khử Re = f(Gr) sẽ có dạng Nu = f(Gr, Pr). Khi đối lưu tự nhiên trong chất lưu ở pha khí, có thể coi Pr = const, nên có Nu=f(Gr). 2) Khi đối lưu cưỡng bức mạnh, lúc Re>104 và chảy rối, thì có thể coi Gr = const, nên phương trình tiêu chuẩn dạng Nu = f(Re, Pr). Khi chất khí đối lưu cưõng bức mạnh, thì Nu = f(Re). 3.6. CÁCH XÁC LẬP 1 CÔNG THỨC THỰC NGHIỆM 3.6.1. các bước thực nghiệm Dể xác lập quan hệ cụ thể của phương trình tiêu chuẩn cho một hiện tượng toả nhiệt cần xét, gnười ta tiến hành các bước thực nghiệm như sau: 1)Lập mô hình thí nghiệm đồng dạng với hiện tượng cần xét 2) Đo mọi giá trị của tất cả các đại lượng vật lý trong các miền giá trị cần khảo sát 38
  12. 3) Lập bảng tính các giá trị tương ứng của các tiêu chuẩn đông dạng Re, Gr, Pr, Nu theo các số liệu đo được tại k điểm đo khác nhau, có dạng bảng (3.6.1.): Điểm ωl gβ∆tl 3 ν αl Re = Gr = pr = Nu = đo ν ν2 a λ 1 Re1 Gr1 Pr1 Nu1 2 Re2 Gr2 Pr2 Nu2 . . . . . . . . . . i Rei Gri Pri Nui . . . . . . . . . . k Rek Grk Prk Nuk Đây là quan hệ cụ thể của Nu=f(Re, Gr, Pr), được trình bày ở dạng bảng số, có thể dùng để xác định Nu và α. 4) Lập công thức thực nghiệm Nu = f(Re, Gr, Pr) tương ứng với bảng giá trị (Re, Gr, Pr, Nu) nói trên bằng phương pháp đồ thị . 3.6.2. phương pháp đồ thị tìm công thức thực nghiệm Từ bảng số liệu (Re, Gr, Pr, Nu), ta có thể tìm một công thức thực nghiệm tương ứng ở dạng Nu = CRemGrnPrp, bằng cách biểu lgNu diễn các điểm đo trên các đồ thị logarit và lần lượt xác định các hằng số n, m, p, C như sau Grk Grf Gri 1) Cố định Pr, Gr, tại các trị số Pri, Grj , biểu diễn k Prk θif Grk Grf Pri Gri điểm đo trên tọa độ (lgNu, lgRe) sẽ được k họ, mỗi họ Grk Grf Pr1 Gri lgRe k đường thẳng dạng lgNu = nlgRe + lg(CGrim.Prjp) với Re1 Rei Rek góc nghiêng γij , và tìm được số mũ n trung bình theo Hình 3.5 Để tìm n 1 k ⎛1 k ⎞ công thức n = ∑ ⎜ k ∑ tgγ ij ⎟ k i =1 ⎜ j=1 ⎟ ⎝ ⎠ 39
  13. Nu 2) Cố định Pr tại các Pri khác nhau, biểu diễn k điểm đo trên toạ độ mới (lg , Re n Nu lgGr) sẽ được 1 họ k đường thẳng dạng lg n = m lg Gr + lg(C Prip ) với góc Re 1 k nghiêng βi , và tìm được trị trung bình của số mũ m theo công thức m = ∑ tgβ i k i =1 lgNu n Re 3) Biểu diễn k điểm đo trên toạ độ ⎛ Nu ⎞ ⎜ lg n m , lg Pr ⎟ , xấp xỉ theo phương pháp bình Prk ⎝ Re Gr ⎠ βi Pri phương nhỏ nhất sẽ được 1 đường thẳng dạng Pr1 Nu lgGr lg = p lg Pr + lg C có góc nghiêng ϕ, và giao Gr1 Gri Grk Re n Gr m Hình 3.6 Để tìm với trục tung tại lgC, theo đó tìm được p=tgϕ và C = 1 k m= ∑ tgβ i k i =1 10lgC . Khi miền biến thiên của Re hoặc Gr khá rộng, để tăng độ chính xác của phép xấp xỉ, người ta chia miền đó ra các khoảng nhỏ, và tìm số n hoặc m ứng với mỗi khoảng. 3.7. CÁC CÔNG THỨC THỰC NGHIỆM TÍNH α 3.7.1. Bài toán tỏa nhiệt và các bước tính α Nu lgRenGrm Bài toán toả nhiệt thường được phát biểu như sau: Tìm hệ số tảo nhiệt α từ bề mặt có vị trí và hình Pri βi dạng cho trước, được đặt trưng bởi kích thước định tính lgC lgPr l , có nhiệt độ tw, đến môi trường chất lưu cho trước có Pr1 Pri Prk nhiệt độ tf và vận tốc ω, nếu có tác nhân cưỡng bức. Hình 3.6 Để tìm Lời giải của bài toán trên là α = λNu/ l , với Nu = p=tgϕ và f(Re, Gr, Pr) tính theo công thức thực nghiệm ứng với mô hình toả nhiệt đã cho, trong đó các giá trị (λ, ν, β, Pr ) được tra từ bảng thông số vật lý của chất lưu đã cho, tại nhiệt độ định tính qui định trong công thức. 40
  14. 1 Với chất lưu ở pha khí , có thể tính β theo công thức β = ,(1/K) [ t ] + 273 3.7.2. Công thức thực nghiệm khi tỏa nhiêt tự nhiên 3.7.2.1. Toả nhiệt tự nhiên trong không gian vô hạn Không gian vô hạn là không gian chứa chất lưu có chiều dày đủ lớn để có thể coi chất lưu chỉ trao đổi nhiệt với riêng mặt đang xét. Công thức thực nghiệm tính α cho bề mặt phẳng hoặc trụ và các mặt hữu hạn có dạng Num = C(GrPr)mn , trong đó nhiệt độ định tính là [t] = (tw + tf)/2 = tm, kích thước định tính là ⎧ ⎪chiãöu cao h cuía màût phàóng hoàûc truû thàóng âæïng ⎪ ⎪ f l = ⎨âæåìng kênh d = 4f / u cuía äúng truû nàòm ngang, coï diãûn têch/ chu vi tiãút diãûn laì ⎪ u ⎪âæåìng kênh d = 6V/F cuía váût hæîu haûn báút kyì, coï thãø têch/ diãûn têch xung quanh laì v ⎪ ⎩ f Các số C và n cho theo bảng sau: (GrPr)m C N 10-3 ÷ 5.102 1,18 1/8 5.102 ÷ 2.107 0,54 1/4 2.107 ÷ 1013 0,135 1/3 Với tấm phẳng nằm ngang thì lấy [ l ] = 4f/u, và khi toả nhiệt lên trên thì lấy α n ↑ = 1,3α l , khitoả nhiệt xuống dưới thì lấy α n ↓ = 0,7α l . 3.7.2.2. Toả nhiệt tự nhiên trong không gian hữu hạn Không gian hữu hạn được hiểu là 1 khe hẹp phẳng hoặc trụ chứa chất lưu có chiều dày δ nhỏ, giữa 2 mặt có nhiệt độ khác nhau tw1 > tw2, khiến cho chất lưu nhận nhiệt từ mặt nóng nhả nhiệt cho mặt lạnh. Lượng nhiệt truyền từ mặt nóng đến mặt lạnh qua chất lưu được tính theo công thức dẫn nhiệt qua vách chất lưu có 2 biên loại 1 và hệ số dẫn nhiệt tương đương λtđ cho theo công thức thực nghiệm: 41
  15. λtđ = λmC(GrPr)mn với [t] = tm = (tw1+ tw2)/2 , ⎧chiãöu daìy δ cuía khe phàóng [l] = ⎨ ⎩chiãöu daìy δ = (d 2 - d 1 )/2 cuía khe truû và C, n cho theo bảng sau: tw1 tw2 (Gr.Pr)m C n δ q tw2
  16. bởi bước ngang S1, bước dọc S2, đường kính ống d, số hàng ống theo phương dòng chảy n. Hệ số tỏa nhiệt α trung bình giữa chất lỏng và mặt ống có thể tính theo công thức sau: 1 0 ,15 ⎛ Pr ⎞ ⎛ d ⎞ - Khi chùm ống song song α= n − 0,5 .0,26. Re 0 , 65 . Prf0, 33 .⎜ f λ 4 fd ⎜ Pr ⎟ ⎟ .⎜ ⎜S ⎟ ⎟ . n ⎝ w ⎠ ⎝ 2 ⎠ d - Khi chùm sole với S1/S2
  17. ε1 là số hiệu chỉnh theo chiều dài, ε1 = f ⎡ 1 , Ref ⎤ cho theo bảng ở phần phụ ⎢d ⎣ ⎥ ⎦ lục. Nếu ống cong với bán kính cong R, chẳng hạn tại đoạn cút hoặc ống xoắn ruột gàm, thì hệ số tỏa nhiệt trong ống cong là: d1 α R = α t .ε R = α t .(1 + 1,77. ) R trong đó αt là hệ số tỏa nhiệt khi ống thẳng tính theo các công thức trên 3.8. TỎA NHIỆT KHI CHẤT LƯU CHUYỂN PHA . 3.8.1. Toả nhiệt khi sôi Sôi là hiện tượng hoá hơi xảy ra trong thể tích chất lỏng. Ở áp suất p = const chất lỏng chi sôi khi đạt tời nhiệt độ sôi ts f(p) , và trong suốt quá trình sôi luôn có t = ts = const. Trên bề mặt nóng có nhiệt độ tw > ts , tuỳ theo cách tạo bọt hơi, hiện tượng sôi được chia ra 2 dạng: Sôi bọt là hiện tượng bọt hơi được tạo ra rời nhau trên mặt nóng, xảy ra khi ∆t = tw – ts không lớn lắm, lúc ∆t < 120C. Sôi màng là hiện tượng bọt hơi tạo ra một màng ngăn cách mặt nóng với chất lỏng, xảy ra khi ∆t lớn, khiến tw tăng cao rất nguy hiểm cho thiết bị và làm giảm α. Hệ số toả nhiệt α [W/m2K] khi nước sôi bọt tại áp suất p ∈ (1 ÷40)bar được tính theo phụ tải nhiệt qua mặt nóng q[W/m2] bởi các công thức thực nghiệm sau: α = 3,14q 0,7 p 0,15 hoặc α = 46( t w − t s ) 2,33 p 0,5 3.8.2. Tỏa nhiệt khi ngưng 3..8.2.1. Mô tả quá trình ngưng Ngưng là hiện tượng hoá lỏng của hơi bảo hoà, xảy ra khi tiếp xúc với vách có nhiệt độ tw
  18. Trên mặt mặt tiếp xúc có nhiệt độ tw < ts, tuỳ theo cách phân bố nước ngưng, hiện tượng ngưng được chia ra 2 dạng. Ngưng giọt là hiện tượng nước ngưng tạo thành các giọt rời nhau rồi rơi xuống, xảy ra khi mặt vách không dính ướt chất lỏng ngưng – Ngưng màng là hiện tượng nước ngưng tạo ra 1màng có độ dày tăng dần theo phương chuyển động xuống dưới. Màng này ngăn cách hơi và vách, làm hệ số α giảm. Ngưng màng xảy ra khi vách dính ướt chất lỏng ngưng, là hiện tượng thường gặp trong thực tế. 3.8.2.2. Công thức Nusselt tính α khi ngưng màng. 1) Phát biểu bài toán Tính hệ số toả nhiệt α khi hơi bảo hoà nhiệt độ ts ngưng màng trên vách đứng cao h, nhiệt độ tw < ts. z 2) Các giả thiết: Để có thể giải phương trình Navier – Stockes tìm o trường vận tốc ω, ta giả thiết rằng chất lỏng ngưng y rλδµ chuyển động đều xuống dưới theo phương x, bỏ tw 1m δx qua hiệu số áp lực. Khi đó phương trình Navier - x dx dy 1m Stockes có dạng h ωx(y) ts d 2ωx ρg − ρg 2 x//g O = ρg + µ → ω x ( y) = ∫∫ dy 2 = y + C1 y + C 2 Hình 3.10 Để tính α dy 2 µ 2µ với C1,C2 tìm theo điều kiện biên ωx(y=0) = C2 = 0, và lực nội ma sát với hơi tại y dωx = δx, là s = µ y =δ x =0 dy ρ.g ρ.g Suy ra µ(C1 - .δ x )=0 → C1 = .δ x µ µ ρ.g → ωx(y) = .(δ x .y − y 2 ) , có dạng đường parabol 2µ 3. Tìm phân bố độ dày lớp nước ngưng δx 45
  19. Tại độ cao x, tốc độ nước trung bình theo y є (0÷δx) là δx 1 ρ.g y3 σ ρ.g 2 ωx = . ∫ ωx ( y).dy = .(δ x .y − ) 0 = 2 x .δ x σx 0 2.µ.δ x 3 3.µ Lưu lượng nước ngưng qua mặt ( δ x x1) tại x là ρ2 .g 3 Gx = ωx .ρ(δ x x1) = δx 3.µ → Lưu lượng nước ngưng thêm trên mặt (dx.1) là ρ2 .g 2 dGx = .δ x .dδ x µ Nhiệt do dGx ngưng tỏa ra bằng nhiệt dẫn qua vách nước dày δ x nên có (t s − t w ) r.dGx = qx(dx.1)= .dx δx λ δ x .r.dG x ρ2 .g.r → dx = = .δ3 .dδ x → tích phân phương trình này ta được λ.(t s − t w ) µ.λ.(t s − t w ) x ρ2 .g.r δ 4 x= . x + C1 với C1 tìm theo điều kiện x( δ x =0) = C1 = 0. Do đó phân µ.λ.∆t 4 1 ⎛ 4.µ.λ.∆t ⎞ 4 bố độ dày lớp nước theo x có dạng δ x = ⎜ 2 .x ⎟ ⎝ ρ .g.r ⎠ 4. Xác định hệ số tỏa nhiệt khi ngưng trên vách đứng Hệ số tỏa nhiệt αx trên độ cao x xác định theo phương trình cân bằng nhiệt λ qx = αx.(ts – tw) = .(t s − t w ) → δx 1 λ ⎛ ρ2 .g.r.λ 3 ⎞ 4 αx = =⎜ ⎟ δ x ⎝ 4.µ.∆t.x ⎠ Hệ số tỏa nhiệt trung bình theo độ cao của vách đưng là 1 1 h ⎛ ρ .g.r.λ ⎞ 1 ⎛ 4 3/ 4 ⎞ 2 3 4 α= h ∫0 α x .dx = ⎜ 4.µ.∆t.x ⎟ . h .⎜ 3 .h ⎟ hay ⎝ ⎠ ⎝ ⎠ 46
  20. 1 1 2 2 ⎛ ρ .g.r.λ ⎞ 2 3 2 2 ⎛ ρ.g.r.λ ⎞ 4 2 3 4 α= .⎜ ⎟ = .⎜ ⎟ ,[W/m K] 3 ⎝ µ.∆t.h ⎠ 3 ⎝ γ.∆t.h ⎠ với ρ,γ,µ= ρ.γ, λ, r là khối lượng riêng, độ nhớt động học, độ nhớt động lực, hệ số dẫn nhiệt, nhiệt ngưng tụ ( hay hóa hơi) của nước bão hòa ở nhiệt độ tm = ( ts+tw)/2 5. tỏa nhiệt khi ngưng trên vách nghiêng hoặc ống ngang Khi vách có chiều cao H đặt dϕ nghiêng so với phương trọng lực 1 ϕ ϕ αϕ góc φ, thì thành phần trọng lực gây ts tw tw αn ρg ts d chuyển động dọc vách là ρ.g.cosφ. Sau tính toán như trên sẽ tìm được hệ Hình 3.11 Để tính αφ trên vách nghiêng số tỏa nhiệt trung bình trên vách và αm trên ống nằm ngang nghiêng là 1 2 2 ⎛ ρ2 .g.r.λ 3 .cosϕ ⎞ 4 α φ= .⎜ ⎟ 3 ⎝ γ.∆t.H ⎠ - Khi ngưng màng trên ống ngang đường kính ngoài d, hệ số tỏa nhiệt trung 1 z 1 π/ 2 ⎛ ρ2 .g.r.λ 3 ⎞ 4 bình αn= .∫− π / 2 α(ϕ)dϕ = 0,72. ⎜ ⎟ π ⎝ γ.∆t.d ⎠ 1m o y Các công thức trên có đủ độ chính xác cho các tw tính toán trong kỹ thuật. Khi có sự khai khác so với x δx ts giả thiết, người ta nhân thêm các hệ số hiệu chỉnh, dx rλδµ được xác định bằng thực nghiệm. h ωx x//g Hình 3.12 Để lập công thức tính Nusselt 47
nguon tai.lieu . vn