Giải bi toán dẫn nhiệt không ổn định tổng quát
bằng phơng pháp phân ly biến số

ThS. Nguyễn đức huy
Bộ môn Kỹ thuật nhiệt - ĐH GTVT
Tóm tắt: Bi báo trình by phơng phân ly biến số để giải bi toán dẫn nhiệt không ổn
định dạng tổng quát.
Summary: The article presents the separation-of-variables method of solving the problem
of unsteady conduction in general form.
Nh đã biết, phơng trình đặc trng của hiện tợng dẫn nhiệt là phơng trình vi phân Fourier:


t
= a
2
t (1)
ở đây: a =


c
- hệ số khuếch tán nhiệt, (m
2
/s);

2
t - toán tử Laplace của nhiệt độ, trong hệ toạ độ Descartes bằng:

2
2
x
t

a
1
(3)
Bằng phơng pháp phân ly biến số, coi t có thể đợc biểu diễn dới dạng tích của 2 hàm
số độc lập:
X = X(x) và T = T() tức là:
t(x, ) =
X(x).T() (4)
ta sẽ tìm đợc nghiệm tổng quát của (1) là:
t = (C
1
cos kx + C
2
sin kx) exp(- k
2
a) (5)
Nếu nhiệt độ phụ thuộc vào 2 hay 3 biến không gian, cho đến nay vẫn phải sử dụng các
phơng pháp gần đúng (sai phân hữu hạn, phần tử hữu hạn) để tìm phơng trình trờng nhiệt độ.
Thực ra, phơng pháp phân ly biến số vẫn có thể đợc sử dụng để giải phơng trình vi
phân Fourier trong trờng hợp tổng quát, khi mà nhiệt độ phụ thuộc cả 3 biến không gian:
t = f (x, y, z, ) (6)

Nhiều tác giả nh Frank Kreith & Mark S. Bohn [1], André B. De Vriendt [2] đã chứng
minh rằng t có thể đợc coi là tích của 3 hàm: X



= X(x, ); Y



),z,y,x(


=
0
),x(


0
),y(


0
),z(


=
X,

.
Y,

.
Z,

(8)
ở đây:

X,



trong hệ toạ độ Descartes. ở thời điểm đầu, nhiệt độ tại mọi điểm trong vật đều bằng t
0
, sau đó
vật đợc làm nguội trong môi trờng chất lỏng có nhiệt độ không đổi t
L
(t
L
< t
0
). Hệ số toả nhiệt
giữa môi trờng và bề mặt vật là (W/m
2
.độ). Cần xác định phân bố nhiệt độ trong vật tại thời
điểm tuỳ ý, nói cách khác, cần tìm nghiệm của phơng trình (1) ở dạng (6).
Với điều kiện nh trên của bài toán, trờng nhiệt độ sẽ đối xứng qua tâm hình hộp, vì thế có
thể chọn tâm hình hộp làm gốc toạ độ.
2L
x
2L
y
2L
z
Phơng trình (1) có dạng:


t
= a
2
t hoặc



= (x, y, z, )
x
Lx m=y
Ly
y
),z,y,x(
m=








= (x, y, z, )

y
Ly m=

z
Lz
z
),z,y,x(
m=

y
),z,y,x(
=








= 0;
0z
z
),z,y,x(
=








= 0
Theo kết quả nhận đợc khi giải phơng trình vi phân Fourier:
2
2
x
t

i
2

2
b
a
)
=


=
1i
i
A

cos
b
i

x. exp(-
i
2

2
b
a
)
với
iii
i

cos.A.A.A cos
y
L
y
y,j

cos
z
L
z
z,k

.
.exp[-(
2
x
2
x,i
L

+
2
y
2
y,j
L

+
2
z

L
z
z,1

.exp[-(
2
x
2
x,1
L

+
2
y
2
y,1
L

+
2
z
2
z,1
L

)a]
Các hệ số A
1,x
,


=


L
y
Bi
z
=


L
z
Hiện nay, các bảng tính sẵn hoặc các đồ thị cho trớc trong [2], [3] có thể giúp ta đơn giản
hóa rất nhiều các bớc giải bài toán dẫn nhiệt không ổn định tổng quát bằng cách cho phép tìm
đợc ngay giá trị của

X,

=
0
),x(


,
Y,

=
0
),y(


C, hệ số dẫn nhiệt = 2,2
W/m.độ, nhiệt dung c = 1930 J/kg.độ, khối lợng riêng
= 913 kg/m
3
. Khối nớc đá tiếp xúc với môi trờng
không khí có nhiệt độ 22
0
C, hệ số toả nhiệt giữa
không khí và bề mặt khối nớc đá = 30 W/m
2
độ.
Cần tính nhiệt độ tại bề mặt khối nớc đá sau 1 phút.
2L
x
= 0, 2 m
2L
y
= 0, 06 m
2L
z
= 0, 1 m
Ta có quan hệ:

=
X,

.
Y,

.

-6
m
2
/ s
Fo
x
=
2
x
L
.a
=
2
6
1,0
60.10.25,1

= 0,0075
- Tra bảng: đợc

X,

= 0,885
Xác định

y,

:
- Tính Bi
y

z,

:
- Tính Bi
z
=


z
L.
=
2,2
05,0.30
= 0,68
- Tính Fo
z
=
2
z
L
.a
=
2
6
05,0
60.10.25,1

= 0,03
- Tra bảng: đợc


= 0,685 t = - 3,3
0
C
Bằng phơng pháp trên, có thể tính đợc nhiệt độ tại mọi điểm của khối nớc đá tại thời
điểm bất kỳ bằng cách thay các giá trị toạ độ và thời gian tơng ứng.
Cũng có thể giải bài toán ngợc (cho trớc nhiệt độ, yêu cầu tính thời gian) bằng cách tính
hoàn toàn tơng tự.
Cần lu ý rằng với dữ kiện nh trờng hợp trên, ta không thể dùng phơng pháp nhiệt độ
đồng nhất (phơng pháp quy tụ) để tính nhiệt độ của khối nớc đá bởi lý do nh sau:
Do chiều dài quy ớc L của khối nớc đá bằng: L =
S
V
= 0,015 m nên tiêu chuẩn Biot của
nó có giá trị: Bi =

L.
=
2,2
015,0.30
= 0,2 . Giá trị này không đáp ứng điều kiện của phơng pháp
nhiệt độ đồng nhất (Bi 0,1) nên không đợc phép áp dụng phơng pháp đó.
Phơng pháp phân ly biến số kết hợp với việc sử dụng các bảng tính và đồ thị cho sẵn sẽ
giúp cho việc giải bài toán dẫn nhiệt không ổn định dạng tổng quát trở nên đơn giản và chính
xác hơn nhiều so với các phơng pháp gần đúng đã biết.

Tài liệu tham khảo
[1]. Frank Kreith & Mark S. Bohn, "Principles of heat transfer", West Publishing Company, 1993.
[2]. André B. De Vriendt, "La transmission de la chaleur", Gaetan Morin éditeur, 1990.
[3]. D. Sivoukhine, "Téplopérédatra", Vshaia Skola, 1988.
[4]. Nguyễn Đức Huy, "Bài giảng môn học Kỹ thuật nhiệt", 2002