Chương 8: DẪN NHIỆT ( Khoa Công Nghệ Nhiệt Lạnh )

Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 66

Chương 8:

DẪN NHIỆT
8.1. Dẫn nhiệt qua vách phẳng (qv = 0)
8.1.1. Vách phẳng một lớp
Xét một vách phẳng đồng chất và đẳng
hướng, chiều dày δ và hệ số dẫn nhiệt λ; vách có
chiều rộng rất lớn so với chiều dày, nhiệt độ hai bề
mặt giữ không đổi là tw1 và tw2. Trong điều kiện này
nhiệt độ chỉ biến thiên theo phương vuông góc với bề
mặt. Nếu chọn trục toạ độ Ox như (hình 10.3) thì
nhiệt độ sẽ không thay đổi theo trục Oy và Oz, nghĩa
là:
∂t ∂t
=
=0
∂y ∂z
Khi các thông số vật lý λ, c, ρ là hằng số thì
phương trình vi phân dẫn nhiệt đối với vách phẳng
một lớp lúc này:
∂²t
(8.1)
=0
∂x ²
Điều kiện biên trong bài toán đang xét (điều kiện biên loại 1) có dạng:
Khi: x = 0 → t = tw1
x = δ → t = tw2
(8.2)

Tích phân (9.22) ta sẽ tìm được quy luật phân bố nhiệt độ trong vách .
Tích phân lần thứ nhất được:
∂t
(8.3)
= C1
∂x
Tích phân lần thứ hai được:
t = C1x + C2
(8.4)
Phương trình (8.4) cho ta thấy khi hệ số dẫn nhiệt không đổi, nhiệt độ trong
vách phân bố theo quy luật đường thẳng.
Hằng số tích phân C1 và C2 trong (10.25) được xác định từ điều kiện biên. Khi:
x = 0;
t = tw;
và: C2 = tw1
t _t
x=δ;
t = tw2;
và: C1 = – w1 w 2
δ
Thay các hằng số tích phân C 1 và C2 vào (8.4) sẽ tìm được phương trình đường
cong phân bố nhiệt độ có dạng:
t _t
t = tw1 – w1 w 2 x
(8.5)
δ
Để xác định mật độ dòng nhiệt theo phương Ox, ta dựa và định luật Fuorier:
Trường ĐHCN Tp. HCM

Khoa CN Nhiệt - Lạnh

Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 67

q = - λ

t _t
dt
dt
. Thay giá trị
= C1 = – w1 w 2 , vào biểu thức định luật
dx
dx
δ

Fuorier:
λ
(tw1 – tw2)
[W/m²]
(8.6)
δ
t _t
hoặc:
q = w1 w 2
[W/m²]
(8.7)
δ/λ
Tỷ số δ/λ được gọi là nhiệt trở dẫn nhiệt của vách một lớp.
Nhiệt lượng truyền qua bề mặt vách F sau khoảng thời gian τ được xác định

theo phương trình sau:
λ
Q = qFτ = (tw1 – tw2)Fτ
[J]
(8.8)
δ
t _t
q
Từ (8. 7) ta có: w1 w 2 =
δ
λ
q
Thay giá trị này vào (8. 6) ta được: t = tw1 – x
(8. 9)
λ
Biểu thức (8.6) và (8.7) tìm được với điều kiện giả thiết λ = const. Thực tế hệ số
dẫn nhiệt biến thiên theo nhiệt độ:
λ = λ0(1+bt)
(8.10)
λ0 – giá trị của hệ số dẫn nhiệt ở 0°C
b - hệ số xác định bằng thực nghiệm.
Theo định luật Fuorier:
dt
dt
q = - λ( t )
= - λ 0 (1 + bt )
(a)
dx
dx
Tách biến, sau đó tích phân trong khoảng từ x = 0 đến

x = δ; và từ t = tw1 đến t = tw2:
t w1 + t w 2
qδ = λ0 1+b
(tw1 – tw2)
(b)
2
q=

(

thừa số:λ0

)

(

t w1 + t w 2
1+b
2

)

t w1
1
λ ( t )dt = λ tb
= t w1 _ t w 2 t ∫
w2

là hệ số dẫn nhiệt trung bình trong khoảng nhiệt độ từ
tw1 đến tw2. Như vậy, mật độ dòng nhiệt được tính theo

biểu thức:
λ
q = tb (tw1 – tw2) [W/m²]
(8.11)
δ
Phương trình (a) sau khi tích phân từ x = 0 đến x, nhiệt độ tích phân từ t w1 đến t
và chỉnh lý lại, ta sẽ tìm được phương trình trường nhiệt độ có dạng sau:
1
1
2qx
t = – + ( + t w1 )² + _
b
b
λοb
Trường ĐHCN Tp. HCM

[°C]

(8.12)

Khoa CN Nhiệt - Lạnh


Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 68

Từ phương trình này ta thấy nhiệt độ trong vách phân bố theo đường cong. Khi
b > 0 bề lồi của đường cong hướng lên và khi b < 0 thì ngược lại.

8.1.2. Vách phẳng nhiều lớp

Giả thiết, có một vách gồm ba lớp hợp thành, chiều dày lớp thứ nhất là δ 1, lớp
thứ hai là δ2 và lớp thứ ba là δ3; hệ số dẫn nhiệt lần lượt bằng λ1, λ2và λ3. Biết nhiệt độ
hai bề mặt ngoài cùng không đổi bằng t w1 và tw4. Bởi vì giữa hai bề mặt có sự tiếp xúc
tốt nên tại đấy hai bề mặt có cùng một nhiệt độ, gọi là tw2 và tw3. (Hình 8.3)

Hình 8.3: Dẫn nhiệt qua vách phẳng nhiều lớp
Ở chế độ nhiệt ổn định dòng nhiệt qua các bề mặt đẳng nhiệt bất kỳ của vách
∂q
bằng nhau, nghĩa là:
=0
∂x
Mật độ dòng nhiệt dẫn qua các lớp:
λ1
q=
(tw1 – tw2)
δ1
q=

λ2
(tw2 – tw3)
δ2

q=

λ3
(tw3 – tw4)
δ3

(c)

Từ (c) ta xác định được độ chênh lệch nhiệt độ qua các lớp:
δ1
tw1 – tw2 = q
λ1
tw2 – tw3 = q

δ2
λ2

tw3 – tw4 = q

δ3
λ3

(d)

Cộng từng vế biểu thức (d) ta được:
Trường ĐHCN Tp. HCM

Khoa CN Nhiệt - Lạnh


Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 69

tw3 – tw4 = q(

δ1 δ 2 δ 3
+
+

)
λ1 λ 2 λ 3

t w1 _ t w 4
q = δ1 δ 2 δ 3
+
+
λ1 λ 2 λ 3

và:

(8.13)

Tương tự có thể suy ra cho vách phẳng nhiều lớp (n lớp)
t w1 _ t w ( n +1)
q=
i=n δ
[W/m²]
∑ i
i =1 λ i

(8.14)

i=n

δi
- là tổng nhiệt trở của n lớp và được gọi là nhiệt trở toàn phần
i =1 λ i
hoặc là nhiệt trở dẫn nhiệt của vách nhiều lớp.
Nếu đem giá trị của q trong (8.14) thay vào (d) ta tính được nhiệt độ của các bề

mặt tiếp xúc trung gian:
δ1
tw2 = tw1 – q
λ1
Trị số: ∑

tw3 = tw2 – q

δ2
δ1 δ 2
= tw1 – q ( + )
λ2
λ1 λ 2

(e)

Tổng quát:
i=k

tw(k – 1) = tw1 – q ∑

δi
λi

[°C]

(8.15)

Để đơn giản cho việc tính toán, ta có thể xem vách nhiều lớp như vách một lớp
i=n

có chiều dày ∆ = ∑ δ i , và nhiệt trở của vách tương đương ấy bằng:
i =1

i=n

∑ δi

i =1

λ td

i=n

δi
i =1 λ i

= ∑

i=n

Suy ra:

λ td =

∑ δi
i =1
i=n δ
i

[W/m độ]

(8.16)

∑

i =1 λ i

Mật độ dòng nhiệt: q =

λ td
(tw1 – tw(n+1)) [W/m²]
Δ

(8.17)

8.2. Dẫn nhiệt qua vách trụ (qv = 0)
8.2.1. Vách trụ một lớp:

Trường ĐHCN Tp. HCM

Khoa CN Nhiệt - Lạnh


Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 70

Xét vách trụ một lớp (ống tròn), đường kính
trong d1 = 2r1 và đường kính ngoài d2 = 2r2 (Hình 8.6).
Nhiệt độ bề mặt vách trong và vách ngoài không đổi

và bằng tw1 và tw2 (điều kiện biên loại 1). Trong khoảng
nhiệt độ cho, hệ số dẫn nhiệt λ có giá trị không đổi.
Phương trình vi phân dẫn nhiệt của vật rắn biểu diễn
trong hệ toạ độ trụ có dạng:
∂ ² t 1 ∂t 1 ∂ ² t ∂ ² t
∇² t =
+
+
+
=0
(8.18)
∂r ² r ∂r r ² ∂φ ² ∂z ²
Trục Oz trùng với trục ống. Trong trường hợp
chiếu dài ống rất lớn so với đường kính ống, nhiệt độ
chỉ thay đổi theo phương bán kính và trường nhiệt độ
là trường một chiều, cho nên:
∂t
∂² t
(a)
= 0 và
=0
∂z
∂z ²
Nhiệt độ bề mặt ngoài và bề mặt trong không thay đổi, nên mặt đẳng nhiệt là
những mặt trụ đồng trục với ống, khi ấy nhiệt độ không biến thiên theo φ; nghĩa là:
∂t
∂² t
= 0 và
=0
(b)

∂φ
∂φ ²
Lúc này (10.38) sẽ có dạng đơn giản hơn:
d ² t 1 dt
(8.19)
+
=0
dr ² r dr
Điều kiện biên:
r = r1→ t = tw1
r = r2→ t = tw2
(8.20)
Giải (8.19) kết hợp với điều kiện biên (8.20) sẽ tìm được phương trình nhiệt độ
trong vách trụ.
dt
Đặt u =
(c)
dr
d ² t du
1 dt u
Khi đó:
;
(d)
=
=
dr ² dr
r dr r
Thay (c) và (d) vào (8.19) được:
du 1
(8.21)

+ u=0
dr r
Tích phân (8.21) ta có:
lnu + lnr = lnC1 ; ur = C1
(e)
dt
dt
dr
thay: u =
→ r = C1 ; nên dt = C1
(f)
dr
dr
dt
Tích phân (f) ta được:
t = C1lnr + C2
(8.22)
Từ điều kiện biên (8.20) ta xac định được C1 và C2:
Trường ĐHCN Tp. HCM

Khoa CN Nhiệt - Lạnh


Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 71

r = r1→ t = tw1 = C1lnr1 + C2
r = r2→ t = tw2 = C1lnr2 + C2
Giải phương trình (g) tìm được C1 và C2:
t _t

ln r1
C1 = w1 w 2
r
;
C2 = tw1 – ( tw1 – tw2 ) r1
ln 1
ln
r2
r2

(g)

Thay C1 và C2 vào (8.22) ta tìm được phương trình đường cong phân bố nhiệt
độ trong vách trụ như sau:
ln(r / r1 )
t = tw1 – ( tw1 – tw2 )
[°C]
(8.23)
ln(r2 / r1 )
hoặc:

t = tw1 – ( tw1 – tw2 )

ln(d / d1 )
ln(d 2 / d1 )

[°C]

(8.24)

Để tính nhiệt lượng truyền qua mặt trụ F trong một đơn vị thời gian, áp dụng
định luật Fuorier:
2 πλL( t w1 _ t w 2 )
dt
Q=– λ F=
[W]
(8.25)
dr
ln(d 2 / d1 )
( t w1 _ t w 2 )
hoặc
Q= 1
[W]
(8.25’)
ln(d 2 / d1 )
2πλL
Trong đó:
L - chiều dài vách trụ.
d
1
ln 2 - nhiệt trở dẫn nhiệt của vách trụ một lớp có độ dài L.
2 πλL d1
Mật độ dòng nhiệt trên một đơn vị dài:
Q t _t
q L = = w1 w 2
d
1
L
ln 2
2 πλ d1

[W/m]

(8.26)

8.2.2. Vách trụ nhiều lớp:
Tương tự như vách phẳng nhiều lớp, vách trụ
nhiều lớp cũng được tạo nên bởi nhiều lớp vật liệu
khác nhau, dòng nhiệt truyền qua vách cũng được
xác định bởi tỷ số giữa độ chênh lệch nhiệt độ toàn
phần và tổng nhiệt trở.
Giả sử một vách trụ gồm 3 lớp (Hình 10.7),
bán kính tương ứng là r1, r2, r3, và r4, hệ số dẫn nhiệt
các lớp λ1, λ2, λ3 là hằng số, nhiệt độ mặt trong cùng
là tw1 và bề mặt ngoài cùng là t w4 không thay đổi (tw1
> tw4 ), nhiệt độ các lớp tiếp xúc chưa biết gọi là t w2 ,
tw3. Ta cần tìm dòng nhiệt truyền qua vách q L và nhiệt
độ tại các lớp tiếp xúc.
Trường ĐHCN Tp. HCM

Khoa CN Nhiệt - Lạnh


Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 72

Trong điều kiện ổn định nhiệt, mật độ dòng nhiệt trên một đơn vị dài q L dẫn qua
∂q L
các lớp đều bằng nhau, nghĩa là:
=0

∂r
Nhiệt lượng dẫn qua các lớp được tính:
2 π( t w1 _ t w 2 )
qL =
1 d2
ln
λ1 d1
qL =

qL =

2π( t w 2 _ t w 3 )
1 d3
ln
λ2 d2

(a)

2π( t w 3 _ t w 4 )
1 d4
ln
λ3 d3

Từ (a) ta tìm được độ chênh lệch nhiệt độ qua các lớp:
qL 1 d2
ln
tw1 – tw2 =
2 π λ1 d 1
tw2 – tw3 =

q L 1 d3
ln
2π λ 2 d 2

tw3 – tw4 =

qL 1 d4
ln
2π λ 3 d 3

(b)

Cộng từng vế của phương trình (b), ta được:
qL 1 d2 1 d3 1 d4
( ln
+
ln
+ ln )
tw1 – tw4 =
2 π λ1 d1 λ 2 d 2 λ 3 d 3
Mật độ dòng nhiệt trên một đơn vị độ dàiqua vách trụ nhiều lớp:
π( t w1 _ t w 4 )
qL =
d
d
d
1
1
1
[W/m]

ln 2 +
ln 3 +
ln 4
2 λ 1 d1 2 λ 2 d 2 2 λ 3 d 3
Tổng quát cho vách có n lớp:
π( t w1 _ t w ( n +1) )
qL = i=n
d
1
∑
ln i +1
di
i =1 2 λ i

[W/m]

(c)

(8.27)

(8.28)

Trong đó:
d
1
ln i +1 có đơn vị [m độ/W] gọi là nhiệt trở đường của mỗi lớp.
2λ i
di
i=n

d
1
ln i +1 gọi là nhiệt trở đường toàn phần của vách nhiều lớp.
di
i =1 2 λ i
∑

Hệ số dẫn nhiệt tương đương đối với vách trụ nhiều lớp bằng nhiệt trở toàn
phần của vách nhiều lớp:
Trường ĐHCN Tp. HCM

Khoa CN Nhiệt - Lạnh


Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 73
i=n 1
d
1
d
ln n +1 = ∑
ln i +1
2λ td
d1
di
i =1 2 λ i

d n +1
d1
λtđ = i = n

1 d
∑ ln i +1
di
i =1 λ i
ln

do đó:

[W/m độ]

(8.29)

Từ (10.49) ta thấy rằng hệ số dẫn nhiệt tương đương λ tđ không những phụ thuộc
vào tính chất vật lý của các lớp (λ 1, λ2, ..., λn ) mà còn phụ thuộc vào chiều dày của
vách.
Nhiệt độ các lớp tiếp xúc có thể tính như sau:
qL 1
d
ln 2
tw2 = tw1 –
π 2 λ1 d 1
tw3 = tw1 –

d
qL 1
d
1
(
ln 2 +
ln 3 )

π 2λ1 d1 2λ 2 d 2

Tổng quát:
qL n 1
d
(∑
ln i +1 )
tw(i+1) = tw1 –
π i =12λ i
di

[°C] (8.30)

8.3 Dẫn nhiệt qua thanh có tiết diện không đổi
8.3.1. Phương trình vi phân dẫn nhiệt và nghiệm của phương trình:
Chúng ta nghiên cứu sự phân bố nhiệt độ dọc theo một thanh dài có tiết diện
không đổi. Gọi tiết diện ngang của thanh là f, chu vi của tiết diện là U, thanh được đặt
trong môi trường có nhiệt độ không đổi là t f, hệ số toả nhiệt trên bề mặt thanh không
đổi bằng α1. Giả thiết vật liệu làm thanh có hệ số dẫn nhiệt lớn, đồng thời tiết diện
ngang của thanh bé, nên nhiệt độ chỉ thay đổi theo phương dọc trục; gọi θ là nhiệt độ
thừa của thanh. Ta có:
θ = t – tf
[°C]
(8.31)
với: tf - nhiệt độ môi trường xung quanh thanh
t - nhiệt độ ở một tiết diện bất kỳ của thanh.
Nếu biết nhiệt độ ở gốc thanh là t1 thì nhiệt độ thừa ở gốc thanh là:
θ1 = t1 – tf
[°C]
(8.32)

Cách gốc thanh một khoảng x ta lấy một phần tử thanh có chiều dài dx thì
phương trình cân bằng nhiệt của phần tử thanh được biểu diễn dưới dạng:
Qx – Qx+dx = dQ
(a)
ở đây: Qx - nhiệt lượng dưa vào bề mặt bên trái của phần tử thanh trong một đơn vị
thời gian.
Qx+dx - nhiệt lượng đưa ra khỏi bề mặt đối diện trong một đơn vị thời gian
dQ - nhiệt lượng toả ra môi trường qua bề mặt xung quanh của phần tử thanh.

Trường ĐHCN Tp. HCM

Khoa CN Nhiệt - Lạnh


Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 74

Theo định luật Fourier: Q x = _ λ
và: Q x + dx = _ λ

dθ
f
dx

d
dθ
(θ + dx )f
dx
dx

d ²θ
(b)
dx ²
mặt khác theo định luật Newton – Ricman có:
dQ = α1θUdx
(c)
Cân bằng vế phải của phương trình (b) và (c), ta sẽ tìm được phương trình vi
phân dẫn nhiệt biểu thị sự phân bố nhiệt độ dọc theo thanh:
d ²θ α1U
(8.33)
=
θ = m ²θ
dx ²
λf
do đó:

Q x _ Q x + dx = λf

α1U
[1/m]
(d)
λf
Nếu α1 và λ không thay đổi trong khoảng nhiệt dộ khảo sát thì m = const.
Nghiệm của phương trình (8.33) có dạng:
θ = C1emx + C2e- mx
(8.34)
Hằng số tích phân C1 và C2 được xác định theo điều kiện biên.
ở đây:

m=±

10.3.4.2. Thanh dài vô hạn
Nhiệt độ ở gốc thanh không đổi, nghĩa là khi: x = 0; θ = θ 1; nếu chiều dài của
thanh: H = ∞, thì: θ = 0 (điều này có nghĩa là khi θ = tx – tf là rất nhỏ), thay điều kiện
biên này vào phương trình (8.34) được:
x = 0 → θ = C1 + C2
x = ∞ → C1e∞ = 0
Từ đẳng thức hai ta rút được C1 = 0, như vậy C2 = θ1, thay C1 và C2 vào phương
trình (8.34) được phương trình đường cong phân bố nhiệt độ dọc theo thanh có dạng:
θ = θ1e- mx
[°C]
(8.35)
Trường ĐHCN Tp. HCM

Khoa CN Nhiệt - Lạnh


Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 75

Như vậy: θ = Ф(θ1, x, α1, λ, f, U)
α1 U
ta thấy m tỷ lệ thuận với hệ số toả nhiệt, với chu vi
λf
của tiết diện và tỷ lệ nghịch với λf là những nhân tố xác định sự dẫn nhiệt dọc theo
thanh.
Nhiệt lượng của thanh truyền cho môi trường sẽ bằng nhiệt lượng dẫn qua gốc
thanh. Dòng nhiệt dẫn qua gốc thanh được tính:
dθ
[W]

(8.36)
Q = _ λf ( ) x = 0
dx
Đạo hàm phương trình (10.59) ta có:
_
dθ
(8.37)
( ) x = 0 = _ me mx θ1 x = 0 = _ mθ1
dx
Thay giá trị của phương trình (8.37) vào (8.36), ta tìm được biểu thức tính nhiệt
lượng của thanh toả ra môi trường xung quanh:
Từ biểu thức: m =

Q = λfmθ1 = θ1 α1Uλf

[W]

(8.38)

10.3.4.3. Thanh dài hữu hạn
1. Đối với thanh dài hữu hạn phương trình vi phân (8.33) và nghiệm của nó
(10.58) vẫn đúng, chỉ có điều kiện biên sẽ khác đi.
x = 0 → θ = θ1
dθ
)
x = H → _ λ(
(8.39)
= θ H α1
dx x = H
( dθ ) x = H = _ α1 θ H

hoặc:
(8.39’)
dx
λ
trong đó:
θH - nhiệt độ thừa cuối thanh
α1 - hệ số toả nhiệt trên bề mặt cuối thanh
Ở tại bề mặt cuối thanh (x = H) ta thấy có sự cân bằng nhiệt lượng do dẫn nhiệt
với nhiệt lượng toả ra cho môi trường xung quanh.
2. Nếu toả nhiệt trên bề mặt cuối thanh có thể bỏ qua (diện tích tiết diện đỉnh
thanh thông thường rất bé so với diện tích xung quanh), thì điều kiện biên (8.39) có thể
viết dưới dạng:
Khi:
x = 0 → θ = θ1
dθ
)
x=H → (
(8.40a)
=0
dx x = H
Phương trình (8.40a) kết hợp với (8.34), ta có:
Khi:
x = 0 → θ = C1 + C2
dθ
)
x=H → (
(8.40b)
=0
dx x = H
do đó:

C1emH + C2e- mH = 0
Giải (8.40b) sẽ xác định được các hằng số tích phân C1 và C2:
Trường ĐHCN Tp. HCM

Khoa CN Nhiệt - Lạnh


Chương 8: Khái Niệm Cơ Bản Về Dẫn Nhiệt
Trang 76

C1 = θ1

_

e

mH
_

e mH

C 2 = θ1

;

_

e mH + e mH
e
+e

Thay C1 và C2 vào (8.34) ta tìm được phương trình đường cong phân bố nhiệt
độ dọc theo thanh:
mH

θ = θ1

mH

e m(H
e

_

x)

mH

+e
+e

_

m(H _ x )

_

mH

_

x
x
ex _ e
Ta biết rằng: e + e
,
và
= ch ( x )
2
2
Khi đó phương trình (8.41) sẽ trở thành:

θ = θ1

ch[m(H _ x )]
ch (mH)

Nhiệt độ thừa cuối thanh là:
θ1
θx=H =
ch (mH)

(8.41)
_

x

= sh ( x )

[°C]

(8.42)

[°C]

(8.42’)

Nhiệt lượng của thanh toả ra môi trường xung quanh bằng nhiệt lượng dẫn qua
dθ
gốc thanh: Q = _ λf ( ) x = 0
[W]
dx
Từ phương trình (8.42) tìm được:
dθ
sh (mH) _
( ) x = 0 = _ θ1m
= mθ1th (mH)
(8.43)
dx
ch (mH)
và:
thay: m =

Q = λ.f.θ1.m.th(mH)

[W]

(8.44)

[W]

(8.45)

α1U
vào (8.44) được:
λf
Q = θ1 α1.U.λ.f th (mH)

Trong trường hợp chiều dài thanh rất lớn thì ch(mH) → ∞ còn th(mH) ≈ 1 và
θx = H = 0; phương trình (8.45) sẽ trở thành (8.38)

Trường ĐHCN Tp. HCM

Khoa CN Nhiệt - Lạnh