GIÁO TRÌNH THIẾT BỊ TRAO ĐỔI NHIỆT

pdf

31

4.7 (
9 lượt)

31447 KB

Nhấn vào bên dưới để tải tài liệu

Đang xem trước 10 trên tổng 31 trang, để tải xuống xem khá đầy đủ hãy nhấn vào bên trên

Chủ đề tương quan

Tài liệu tương tự

Nội dung

CHƯƠNG 1:

CÁC KHÁI NIỆM CHUNG

1.1. ĐỊNH NGHĨA VÀ PHÂN LOẠI CÁC THIẾT BỊ TRAO ĐỔI NHIỆT.
1.1.1. Các định nghĩa.
Thiết bị trao đổi nhiệt (TBTĐN) là thiết bị trong đó thực hiện sự trao đổi nhiệt giữa
chất cần gia công với chất mang nhiệt hoặc lạnh.
Chất mang nhiệt hoặc lạnh được gọi chung là môi chất có nhiệt độ cao hơn hoặc thấp
hơn chất gia công, dùng để nung nóng hoặc làm nguội chất gia công.
Chất gia công và môi chất thường ở pha lỏng hoặc hơi, gọi chung là chất lỏng. Các
chất này có nhiệt độ khác nhau.
Để phân biệt mỗi thông số ϕ là của chất lỏng nóng hay chất lỏng lạnh, đi vào hay ra
khỏi thiết bị, người ta quy ước:
– Dùng chỉ số 1 để chỉ chất lỏng nóng: ϕ1.
– Dùng chỉ số 2 để chỉ chất lỏng nóng: ϕ2.
– Dùng dấu “ ′ ” để chỉ thông số vào thiết bị: ϕ1′; ϕ2′.
– Dùng dấu “ ″ ” để chỉ thông số ra thiết bị: ϕ1″; ϕ2″.
Ví dụ:
Cl1 t1′
t”2

t1”
Cl 2 t’2

Hình 1.1. Sơ đồ khối của TBTĐN
1.1.2. Phân loại các TBTĐN.
1.1.2.1. Phân loại theo nguyên lý làm việc của TBTĐN.
1) TBTĐN tiếp xúc (hay hỗn hợp), là loại TBTĐN trong đó chất gia công và môi chất
tiếp xúc nhau, thực hiện cả quá trình trao đổi nhiệt và trao đổi chất, tạo ra một hỗn hợp. Ví
dụ bình gia nhiệt nước bằng cách sục 1 dòng hơi.
2) TBTĐN hồi nhiệt, là loại thiết bị TĐN có mặt trao đổi nhiệt được quay, khi tiếp
xúc chất lỏng 1 mặt nhận nhiệt, khi tiếp xúc chất lỏng 2 mặt toả nhiệt. Quá trình TĐN là
không ổn định và trong mặt trao đổi nhiệt có sự dao động nhiệt. Ví dụ: bộ sấy không khí
quay trong lò hơi nhà máy nhiệt điện.

3) TBTĐN vách ngăn, là loại TBTĐN có vách rắn ngăn cách chất lỏng nóng và chất
lỏng lạnh và 2 chất lỏng TĐN theo kiểu truyền nhiệt. Loại TBTĐN vách ngăn bảo đảm độ
kín tuyệt đối giữa hai chất, làm cho chất gia công được tinh khiết và vệ sinh, an toàn, do đó
được sử dụng rộng rãi trong mọi công nghệ.
4) TBTĐN kiểu ống nhiệt, là loại TBTĐN dùng ống nhiệt để truyền tải nhiệt từ chất
lỏng nóng đến chất lỏng lạnh. Môi chất trong các ống nhiệt nhân nhiệt từ chất lỏng 1, sôi
và hoá hơi thành hơi bão hoà khô, truyền đến vùng tiếp xúc chất lỏng 2, ngưng thành lỏng
rồi quay về vùng nóng để lặp lại chu trình. Trong ống nhiệt, môi chất sôi, ngưng và chuyển
động tuần hoàn, tải 1 lượng nhiệt lớn từ chất lỏng 1 đến chất lỏng 2.
a. Bình gia nhiệt hỗn hợp

b. Thùng gia nhiệt khí hồi nhiệt

c. Bình ngưng ống nước

d. Lò hơi ống nhiệt

Hình 1.2. Các loại TBTDN phân theo nguyên lý làm việc.
1.1.2.2. Phân loại TBTĐN theo sơ đồ chuyển động chất lỏng, với loại TBTĐN có vách
ngăn.
a. Sơ đồ song song cùng chiều.
b. Sơ đồ song song ngược chiều.
c. Sơ đồ song song đổi chiều.
d. Sơ đồ giao nhau 1 lần.
e. Sơ đồ giao nhau nhiều lần.
Hình 1.3. Các sơ đồ chuyển động chất lỏng trong TBTDN.
1.1.2.3. Phân loại TBTĐN theo thời gian.
– Thường phân ra 2 loại: Thiêt bị liên tục (ví dụ bình ngưng, calorife) và thiết bị làm việc
theo chu kỳ (nồi nấu, thiết bị sấy theo mẻ).
1.1.2.4. Phân loại TBTĐN theo công dụng.
– Thiết bị gia nhiệt dùng để gia nhiệt cho sản phẩm (Ví dụ nồi nấu, lò hơi).
– Thiết bị làm mát để làm nguội sản phẩm đến nhiệt độ môi trường (Ví dụ tháp giải nhiệt
nước, bình làm mát dầu)
– Thiết bị lạnh để hạ nhiệt độ sản phẩm đến nhiệt độ nhỏ hơn môi trường (Ví dụ tủ cấp
đông, tủ lạnh).
1.2. CÁC YÊU CẦU KỸ THUẬT CHUNG CHO MỌI TBTĐN.
1.2.1. Các yêu cầu kỹ thuật chung cho mọi TBTĐN.

Khi thiết kế chế tạo hoặc lựa chọn trang bị, các TBTĐN cần đạt các yêu cầu kỹ thuật
chính sau đây.
−1

⎛ 1 δ
1 ⎞
⎟⎟ cần phải lớn, để tăng cường công suất TĐN
1) Hệ số truyền nhiệt k = ⎜⎜ + +
⎝ α1 λ α 2 ⎠

Q = kF ∆t. Muốn tăng k, cần tăng λ, α1, α2, nhất là tăng min (α1, α2) và giảm chiều dày δ
của vách, không làm vách nhiều lớp.
2) Giảm trở kháng thuỷ lực trên dòng chảy các môi chất ∆p1, ∆p2, để giảm công suất
bơm quạt p = ∆pV/η. Muốn vậy cần giảm độ nhớt của chất lỏng, giảm tốc độ ω, giảm các
tổn thất cục bộ đến mức có thể.
3) Tăng diện tích mặt trao đổi nhiệt, là mặt có 2 phía tiếp xúc trực tiếp chất lỏng nóng
và chất lỏng lạnh để tăng công suất Q = kF ∆t .
4) Bảo đảm an toàn tại áp suất và nhiệt độ làm việc cao nhất và có tuổi thọ cao. Muốn
vậy phải chọn kim loại đủ bền ở p, t làm việc, tính toán độ dày δ theo các quy tắc sức bền.
5) Bảo đảm độ kín giữa 2 chất lỏng với nhau và với môi trường bên ngoài, để gữ độ
tính nhiệt của sản phẩm và vệ sinh an toàn cho môi trường.
6) Cấu tạo đơn giản, gọn nhẹ, dễ vận chuyển, dễ lắp ráp, vận hành, dễ kiểm tra, điều
khiển và dễ vệ sinh, bảo dưỡng.
1.2.2. Các nguyên tắc lựa chọn môi chất.
Môi chất là chất trung gian dùng để gia nhiệt hay làm lạnh sản phẩm trong TBTĐN.
Môi chất được phân loại theo mục đích sử dụng (Môi chất tải nhiệt như hơi nước, môi chất
tải lạnh như dung dịch NaCl, môi chất lạnh như NH3), theo pha khi làm việc (1 pha, 2 pha,
3 pha), theo nhiệt độ làm việc ∆tlv = (tmin ÷ tmax) (nhiệt độ rất cao, cao, trung bình, thấp, rất
thấp).
Việc lựa chọn môi chất cần đạt các yêu cầu sau:
1) Chọn môi chất có ρ, c, λ, r lớn để có d, k lớn nhằm tăng cường trao đổi nhiệt.
2) Chất có nhiệt độ nóng chảy tnc, nhiệt độ sôi ts và có pha thích hợp với ∆tlàm việc và
áp suất làm việc.
3) Chất có độ nhớt ν nhỏ để giảm ∆p.
4) Chất không gây cháy nổ, ít độc hại, ít ăn mòn, không chứa tạp chất (cặn, bụi).

Nhiệt độ làm việc, áp suất làm việc và khả năng trao đổi nhiệt của 1 số môi chất
thông dụng được giới thiệu ở bảng 1 và bảng 2.
Bảng 1 – Khoảng nhiệt độ và áp suất làm việc của các môi chất.
tlv [0C]

plv tuyệt đôi [bar]

Khí H2

≥ 273

≤ 10

Khí O2, N2, không khí

≥ 210

≤ 200

Khí metal CH4

-160 ÷ -100

≤ 40

Khí etal, etylen, freon

-150 ÷ -70

≤ 40

Freon 12, 22, NH3, CO2

-70 ÷ 0

≤ 15

Nước muối (dung dịch NaCl)

-50 ÷ 0

≤3

Freon 11, 12, 113, 114

-10 ÷ 0

≤3

Dầu

0 ÷ 215

≤2

Nước H2O

0 ÷ 374

1 ÷ 225

Hơi nước

0 ÷ 650

1 ÷ 300

Hỗn hợp difenyl

260 ÷ 350

1÷6

Thuỷ ngân

350 ÷ 500

1÷9

Khói nóng

450 ÷ 1000

≤1

Chất rắn (samot)

≤ 1500

≤1

Plasma t0 thấp

≤ 3500

≤1

Môi chất

Bảng 2 – Khả năng trao đổi nhiệt của các môi chất.
Quá trình TĐN
Đốt nóng hoặc làm nguội

Sôi bọt

Môi chất

α

[W/m2K]

α min

α max

Khí

1

60

Hơi quá nhiệt

20

120

Dầu

60

1.700

Nước

200

10.000

Chất lỏng hữu cơ

600

10.000

Nước
Ngưng màng

Hơi chất hữu cơ
Hơi nước

6.000

50.000

600

2.500

5.000

20.000

1.2.3. Chọn sơ đồ chuyển động của 2 chất lỏng.
Các kết quả thực nghiệm cho biết, hệ số toả nhiệt α khi dòng chất lỏng cắt ngang ống
lớn hơn, khi dòng chảy dọc ống, αn > αd, còn trở kháng thuỷ lực thì ∆pn > ∆pd.
Qua phân tích, Berman cho biết:
1) Với chất lỏng, khi Nu/Pr < 61 thì nên cho chảy dọc ống (ưu tiên ngược chiều, đảo chiều).; khi Nu/Pr > 61 nên cho chảy cắt ngang ống (ưu tiên giao nhiều lần).
2) Với chất khí, khi Re ∈ [4.103 ÷ 4.104] nên cho chảy cắt ngang ống.
1.2.4. Các nguyên tắc chọn chất lỏng chảy trong ống.
Khi cần chọn 1 chất lỏng cho đi trong ống thì ưu tiên cho:
1) Chất lỏng có lưu lượng thể tích V (m3/s) nhỏ hơn, để giảm vận tốc ω = V/ρ, do đó
giảm ∆p và công suất bơm.
2) Chất lỏng có độ nhớt cao hơn để để tăng ∆p lúc bơm.
3) Chất lỏng có (p, t)lv lớn để vỏ thiết bị không chịu (p, t) cao, thiết bị sẽ nhẹ và rẻ
hơn.
4) Chất lỏng độc hại, bẩn, gây ăn mòn, để dễ làm kín, dễ vệ sinh và ít tốn vật liệu bị
ăn mòn hóa chất.
1.2.5. Chọn tốc độ dòng môi chất.
Khi tốc độ ω tăng thì α, k tăng, làm TĐN tốt hơn, nhưng cũng làm tăng ∆p và công
suất tiêu hao cho bơm quạt. Do đó, cần chọn một vận tốc hợp lý để giảm chi phí vận hành,
tăng hiệu quả kinh tế. Bài toán tối ưu cho biết, nên chọn ω hợp lý cho môi chất theo bảng 3
sau đây:
Bảng 3. Khoảng giá trị hợp lý của vận tốc môi chất.
Môi chất

ωtư (m/s)

Chất lỏng có ν nhỏ (H2O, glycol)

0.5 ÷3

Chất lỏng nhớt cao (dầu, dd NaCl)

0.2 ÷1

Khí + bụi ở pk ( khói, khí bụi)

6 ÷ 10

Khí sạch ở pk ( không khí )

12 ÷ 16

Khí nén ở p > pk ( khí nén)

15 ÷ 30

Hơi bảo hoà

30 ÷ 50

Hơi quá nhiệt

30 ÷ 75

1.3. PHƯƠNG TRÌNH CƠ BẢN CỦA TBTĐN.
Khi tính toán các thiết bị trao đổi nhiệt, người ta luôn dựa vào 2 phương trình sau
đây gọi là phương trình cơ bản của TBTĐN.
1.3.1. Phương trình cân bằng nhiệt (CBN).
Phương trình cân bằng nhiệt là phương trình mô tả định luật bảo toàn và biến hoá
năng lượng cho TBTĐN.
1.3.1.1. Phương trình cân bằng nhiệt tổng quát.
Phương trình cân bằng nhiệt tổng quát liên hệ các hệ số entanpi ra vào thiết bị với
nhiệt truyền qua vỏ thiết bị ra môi trường và biến thiên nội năng của thiết bị:
(Hiệu entanpi ra – vào của chất lỏng 1) + (Hiệu entanpi ra- vào của chất lỏng 2) +
(Nhiệt truyền qua vỏ thiết bị ra môi trường) + (Biến thiên nội năng của thiết bị) = 0.
Ở dạng tích phân, phương trình cân bằng nhiệt tổng quát có dạng:
∑Q = ( ∆I1 + ∆I2 + Qk) τ + ∆U = 0,
Trong đó:
∆I1 = G1( i1// − i1/ ) = G1Cp1( t1// − t1/ ) < 0, do chất lỏng 1 toả nhiệt. ∆I2 = G2( i2// − i2/ ) = G2Cp2( t 2// − t 2/ ) > 0, do chất lỏng 2 thu nhiệt.
Qk = ∑kiFi( t – tf), [W] là nhiệt truyền từ chất lỏng có nhiệt độ t qua các diện tích Fi
của vỏ thiết bị ra môi trường nhiệt độ tf. Với thiết bị gia nhiệt thường có: t > tf nên Qk > 0
tức môi trường nhận nhiệt. Với thiết bị làm lạnh, thường t < tf nên Qk < 0 tức môi trường toả nhiệt vào thiết bị. τ [s] là thời gian từ khi khởi động thiết bị ở nhiệt độ t0 đến nhiệt độ tτ nào đó. ∆U = ∑ViρiCi( tτ - t0), [J] là biến thiên nội năng của các chi tiết tạo ra thiết bị. Trong thiết bị gia nhiệt, thường tτ > t0 nên ∆U > 0; trong thiết bị làm lạnh, thường tτ < t0 nên ∆U < 0. Nếu tính từ khi thiết bị đã làm việc ổn định, thì ∆U = 0. Nếu đặt W = GCp, [W/K] là đương lượng nước của chất lỏng thì liên hệ giữa W, lưu lượng G(kg/s); khối lượng riêng ρ[kg/m3], nhiệt dung riêng Cp[J/kgK], vận tốc ω[m/s] của chất lỏng với tiết diện dòng chất lỏng f sẽ có dạng: W = GCp = ρωfCp Trong đó: V = ωf [m3/s] được gọi là lưu lượng thể tích. Phương trình CBN tích phân tổng quát, liên hệ các thông số như trên sẽ có dạng: [ρ1ω1f1( i1// − i1/ ) + ρ2ω2f2( i2// − i2/ ) + ∑kiFi( t - tf)] τ + ∑ρiViCi( tτ - t0). Phương trình này cho phép tìm được một đại lượng chưa biết nào đó, ví dụ thời gian τ để khởi động thiết bị khi có thể xác định tất cả các đại lượng còn lại. Khi xét cân bằng nhiệt qua 1 vi phân dF của diện tích trao đổi nhiệt của thiết bị thì phương trình cân bằng nhiệt tổng quát có dạng vi phân sau: ρ1ω1f1di1 + ρ2ω2f2di2 + ∑kiFi( t - tf)dFi + ∑ρiViCi dt =0. dτ Đây là phương trình vi phân cân bằng công suất nhiệt trao đổi qua diện tích dF của TBTĐN. Nó cho phép tìm được luật biến thiên theo thời gian τ của nhiệt độ các chất lỏng, thông qua di = Cpdt. 1.3.1.2. Các phương trình cân bằng nhiệt đặc biệt. 1) Khi thiết bị cách nhiệt tốt với môi trường: Coi Qk = 0, (∆I1 + ∆I2) τ + ∆U = 0. 2) Khi TBTĐN làm việc ổn định, coi ∆U = 0, ∆I1 + ∆I2 + Qk = 0. 3) Khi thiết bị được cách nhiệt, làm việc ổn định thì: ∆I1 + ∆I2 = 0 hay G1( i1/ − i1// ) = G2Cp2( i2// − i2/ ) G1Cp1( t1/ − t1// ) = G2Cp2( t 2// − t 2/ ) hay W1( t1/ − t1// ) = W2( t 2// − t 2/ ) Dạng vi phân của phương trình cân bằng nhiệt khi đó là W1dt1 = W2dt2. 4) Khi ∆U = 0, Qk = 0 nếu các chất lỏng có sự chuyển pha trong TBTĐN, từ chất lỏng Cp đến sôi ở ts nhận nhiệt r, rồi quá nhiệt đến hơi có nhiệt dung riêng Cph, thì phương trình cân bằng nhiệt có dạng: G1 [C p1h (t1/ − t s1 ) + r1 + C p1 (t s1 − t1// )] = G2 [C p 2 (t s 2 − t 2/ ) + r2 + C p 2 (t 2// − t s 2 )] Hình 1.4. Phân bố nhiệt độ các chất lỏng khi chuyển pha trong TBTĐN cùng chiều. Ví dụ: + Phương trình cân bằng nhiệt trong lò hơi: G1Cp1 (t1/ − t1// ) = G2 [C p 2 (t s 2 − t 2/ ) + r2 + C p 2 h (t 2// − t s 2 )] với: 1- khối nóng, 2- H2O. + Phương trình cân bằng nhiệt cho bình ngưng: 1- hơi ngưng, 2- nước làm mát. G1 [C p1h (t1/ − t s1 ) + r1 + C p1 (t s1 − t1// )] = G2C p 2 (t 2/ − t 2// ) . 1.3.2. Phương trình truyền nhiệt. Phương trình truyền nhiệt là những phương trình mô tả lượng nhiệt trao đổi giữa 2 chất lỏng qua mặt TĐN bằng phương thức truyền nhiệt. 1) Dạng vi phân. Lượng nhiệt δQ truyền từ chất lỏng nóng nhiệt độ t1 qua diện tích dFx của mặt TĐN đến chất lỏng lạnh nhiệt độ t2 là: δQ = k(t1 – t2)dFx = k∆txdFx, W Trong đó: ⎛ 1 δ 1 +∑ i k = ⎜⎜ + λi ⎝ α1 α 2 −1 ⎞ ⎟⎟, [W/m2K] là hệ số truyền nhiệt qua vách thường được coi là ⎠ không đổi trong mặt F. ∆tx = t1(x) – t2(x) = f(Fx) là độ chênh nhiệt độ của 2 chất lỏng hai bên mặt dFx, phụ thuộc vị trí Fx. 2) Dạng tích phân. Lượng nhiệt Q truyền từ chất lỏng 1 qua diện tích TĐN F đến chất lỏng 2 là: F Q = ∫ k∆t x (Fx )dFx = k ∫ ∆t x (Fx )dFx = kF ∆t, [W]. 0 Với: ∆t = 1 ∆t x (Fx )dFx, gọi là độ chênh trung bình trên mặt F của nhiệt độ 2 chất lỏng. F∫ 1.4. XÁC ĐỊNH ĐỘ CHÊNH NHIỆT ĐỘ TRUNG BÌNH ∆t . Giá trị ∆t phụ thuộc vào t 1/, t 1//, t 2/, t 2// và loại sơ đồ chuyển động của 2 chất lỏng. 1.4.1. Sơ đồ song song ngược chiều. Phương trình cân bằng nhiệt và truyền nhiệt qua dFx của TBTĐN song song ngược chiều, theo hình 1.5 có dạng: ⎧δQ = − W1dt1 = − W2 dt 2 ⎨ ⎩δQ = k∆t x dFx ⎛ 1 1 ⎞ ⎟⎟δQ − W W 2 ⎠ ⎝ 1 Theo đó có: dt1-dt2 = - ⎜⎜ ⎛ 1 1 ⎞ ⎟⎟, [k/W]. − ⎝ W1 W2 ⎠ hay d ∆t x = - mk∆txdFx, với m = ⎜⎜ Nếu m và k không đổi thì: Do đó: ∆t x F x d∆t x ∆t = − mk ∫∆t ∆t x ∫0 dFx hay ln ∆tX0 = −mkFx . 0 ∆tx(Fx) = ∆t0exp(-mkFx). F Theo định nghĩa ∆t : ∆t ∆t 0 1 ∆t x ( Fx )dFx = 0 ∫ exp(− mkFx )dFx = e −mkF − 1 . ∫ F0 F − mkF ∆t = ( Thay ∆tF = ∆t0exp(-mkF) vào trên sẽ được: ∆t = ∆t 0 ∆t ln F ∆t 0 ) ⎛ ∆t F ⎞ ∆t − ∆t 0 ⎜⎜ − 1⎟⎟ = F ∆t ⎝ ∆t0 ⎠ ln F ∆t0 ⎧⎪∆t s = t1/ − t 2// với ⎨ ⎪⎩∆t F = t1// − t 2/ t t'1 ∆t o t1 t"2 C1 ∆t x dt 2 λ t2 dt 1 t"1 C2 ∆t F t'2 d Fx O Fx F Fx Hình 1.5. Sơ đồ trao đổi nhiệt 2 chất lỏng song song ngược chiều. 1.4.2. Sơ đồ song song cùng chiều. ⎧δQ = − W1dt1 = W2 dt 2 Phương trình cân bằng nhiệt và truyền nhiệt dFx là: ⎨ ⎩δQ = −k∆t x dFx ⎛ 1 1 ⎞ ⎟⎟k∆tdFx = -mk∆txdFx và biến đổi như trên sẽ Sau khi đưa về dạng: d∆tx= - ⎜⎜ + ⎝ W1 W2 ⎠ thu được: ⎧⎪∆t0 = t1/ − t 2/ ∆t F − ∆t 0 ∆t = với ⎨ ∆t F ⎪⎩∆t F = t1// − t 2// ln ∆t 0 t t'1 t1 C1 dt 1 ∆t x C2 ∆t o t2 t'2 O t"1 ∆t F dt 2 t"2 d Fx Fx Fx F Hình 1.6. Sơ đồ trao đổi nhiệt 2 chất lỏng song song cùng chiều. Các công thức trên dùng khi: ∆t0 ≠ ∆tF ≠ 0. Các công thức đặc biệt khác có thể tính ∆t theo: ⎧∆t khi ∆t0 = ∆t F ≠ 0 ⎪1 ⎪ ∆t = ⎨ (∆t 0 + ∆t F ) ⎪2 ⎪⎩0 khi ∆t0 = ∆t F = 0 1.4.3. Các sơ đồ khác. Để tính ∆t cho các sơ đồ khác (song song đổi chiều, giao nhau 1 hay n lần), ta tính ∆t theo sơ đồ song song ngược chiều rồi nhân với hệ số ε∆t, được xác định bằng thực nghiệm và cho ở dạng đồ thị. ε∆t = f( p = t 2// − t 2/ t1/ − t1// , R = , loại sơ đồ): ∆t = ∆t ↑↓. ε∆t (P, R,loại sơ đồ). t1/ − t 2/ t 2// − t 2/ 1.5. CÁC CHỈ TIÊU CHẤT LƯỢNG CỦA TBTĐN . Để đánh giá chất lượng của TBTĐN, người ta dựa vào các chỉ tiêu sau đây: 1.5.1. Chỉ tiêu về năng lượng. Để đặc trưng cho một công suất nhiệt thu được ứng với 1kW điện tiêu hao khi vận hành bơm quạt của thiết bị, người ta dùng chỉ tiêu năng lượng E0, được định nghĩa: E0 = Công suất nhiệt sản phẩm thu được từ môi chất Tổng công suất để bơm quạt sản phẩm và môi chất

Source: https://vvc.vn
Category : Đồ Cũ