›
›

Bài thực hành môn Mô hình hóa, mô phỏng và tối ưu hóa các quá trình hóa học - Bài thực hành số 4: Khảo sát một quá trình truyền nhiệt

Download
Bài thực hành môn học Mô hình hóa, mô phỏng và tối ưu hóa các quá trình hóa học  
BÀI THỰC HÀNH SỐ 4  
KHẢO SÁT MỘT QUÁ TRÌNH TRUYỀN NHIỆT  
Mục đích của bài thực hành này là mô phỏng, giải quyết các bài toán trong CNHH. Cụ thể :  
1. Mô hình hóa, rời rạc hóa quá trình dẫn nhiệt trong một thanh kim loại được mô tả  
bằng phương trình vi phân đạo hàm riêng PDE dùng phương pháp sai phân hữu hạn để  
nhận được một hệ ODE.  
2. Tìm nghiệm số của hệ ODE này.  
3. Xấp xỉ nghiệm của bài toán dùng phương pháp nội suy Lagrange.  
4. Kết luận về ảnh hưởng của số điểm nút chọn trên các nghiệm số.  
1. Mô tả quá trình truyền nhiệt:  
Mục đích của bài thực hành này là để khảo sát động học của nhiệt độ T(x,t) phân bố trong  
một thanh kim loại. Chúng ta giả sử rằng nhiệt độ tại x = 0 của thanh kim loại được duy trì ở  
nhiệt độ môi trường bên ngoài T 0 và đầu kia được duy trì ở nhiệt độ T1 có thể thay đổi theo  
thời gian.  
T1(t)  
T 0  
R
T(x,t)  
x
x+dx  
0
L
x
Giả sử truyền nhiệt chỉ xảy ra theo chiều của trục x. Qua bài thực hành này, chúng ta biết cách  
mô phỏng tiến triển của biên dạng nhiệt độ T(x,t) trong thanh kim loại (giả sử thanh đồng  
nhất và bỏ qua sự giản nở nhiệt).  
Câu hỏi 1 (Mô hình hóa toán học): Dùng định luật Fourier và với cân bằng năng lượng được  
viết cho phẩn tử thể tích vô cùng bé giũa x và x+dx, chỉ ra rằng động học của nhiệt độ trong  
thanh kim loại được chi phối bởi phương trình đạo hàm riêng sau:  
∂T(x,t)  
∂t  
∂2T(x,t)  
∂x2  
ρc  
= λ  
(1)  
Với ρ (g/cm3) là khối lượng thể tích ; c (J/g/K) vàλ (W/cm/K) tương ứng là nhiệt dung riêng  
khối lượng và hệ số dẫn nhiệt của thanh kim loại. □  
Phương trình (1) có thể viết lại tương đương như sau :  
∂T(x,t)  
∂t  
∂2T(x,t)  
∂x2  
= D  
(2)  
Copyright © by Hoàng Ngọc Hà  
Bài thực hành môn học Mô hình hóa, mô phỏng và tối ưu hóa các quá trình hóa học  
λ
với D =  
. Nghiệm nếu có của phương trình (2) cần bổ sung các điều kiện ban đầu và điều  
ρc  
kiện biên:  
• Với điều kiện ban đầu  
T(x,0) = T init (x)  
• Với điều kiện biên  
(3)  
(4)  
0
⎧
T(0,t) = T (t)  
⎪
⎨
1
⎪
T(L,t) = T (t)  
⎩
Trong phương trình (1) t ∈  
đủ hệ thống.  
[0, + ∞  
)
, x∈  
[
0, L  
]
. Phương trình (2) với (3) và (4) miêu tả đầy  
Bảng dưới đây cho các dữ liệu các tham số của hệ phản ứng nghiên cứu :  
D (cm2/s)  
L (cm)  
R (cm)  
T0 (x) (K)  
T 0 (t) (K)  
T1(t) (K)  
(Sinh viên tùy chọn giá trị các tham số sao cho phù hợp. Có điểm ưu tiên cho việc chọn lựa  
tốt, sáng tạo)  
2. Nghiệm số dùng phương pháp sai phân hữu hạn bước trung tâm:  
Trước tiên chúng ta xác định biểu diễn đại số của phương trình vi phân (1) mà chúng ta sẽ  
nhận được bởi rời rạc hóa nó dùng phương pháp sai phân hữu hạn tại từng điểm nút.  
Gọi h là bước rời rạc không gian sao cho khoảng  
[
0, L được chia thành N khoảng con có  
]
L
cùng chiều dài h, có nghĩa rằngh = . Chúng ta ký hiệu xi là các điểm rời rạc :  
N
x0 = 0, x1 = h, x2 = 2h,..., xN = L  
Các điểm này phân bố hình học như sau :  
xN = L  
x2 = 2h  
x = h  
...  
x0 = 0  
x
x
1x  
123  
x
x
h
Copyright © by Hoàng Ngọc Hà  
Bài thực hành môn học Mô hình hóa, mô phỏng và tối ưu hóa các quá trình hóa học  
Tiếp theo, để đơn giản cách trình bày chúng ta sẽ kí  
hiệu  
T0 (t) = T(0,t),T (t) = T(h,t),...,TN (t) = T(Nh,t) . Tại các điểm trung gian x1, x2 ,..., xN−1 chúng  
1
ta sử dụng các biểu thức sau để xấp xỉ đạo hàm bậc hai và bậc 1 như sau :  
∂2T(x,t)  
Ti+1(t) +Ti−1(t) − 2Ti (t)  
≅
, i =1,..., N −1  
(5)  
∂x2  
h2  
x=xi  
và  
Ti+1(t) −Ti−1(t)  
2h  
∂T  
∂x  
≅
, i =1,..., N −1  
x=xi  
Cuối cùng chúng ta diễn giải điều kiện biên (4) tại các điểm đầu-cuối như sau :  
T0 (t) = T 0 (t) và TN (t) = T1(t)  
(6a)  
(6b)  
và điều kiện ban đầu (3) tại các điểm rời rạc như sau  
T(ih,0) = T init (ih), i = 0...N  
Sau cùng, để thuân lợi cho biểu diễn chúng ta ký hiệu :  
T
⎡
⎤
⎥
⎥
⎥
⎥
T(h,t)  
⎡
⎤
⎥
⎥
⎥
⎥
⎦
1t  
⎢
⎢
T2t  
T(2h,t)  
⎢
⎢
X (t) =  
=
(7)  
⎢
⎢
M
M
⎢
⎢
T
T((N −1)h,t)  
⎢
⎥
⎦
(N−1)t  
⎣
⎣
3. Yêu cầu :  
Chọn số điểm nút (N>2) N= ?  
Câu hỏi 2 :  
a) Dùng xấp xỉ (5), chỉ ra rằng phương trình (2) tại điểm rời rạc i trở thành :  
dTit  
dt  
Ti+1(t) +Ti−1(t) − 2Ti (t)  
(8)  
= D  
, i =1,..., N −1  
h2  
b) Tính đến các điều kiên biên (6a), chỉ ra rằng phương trình (8) có thể viết dưới dạng hệ  
phương trình vi phân thường sau :  
dX (t)  
= AX (t) + b  
dt  
(9)  
với X (t) là vector được cho trong (7) và giá trị ban đầu X (t = 0) có được dùng (6b). Xác  
định các ma trận A, b ?  
Copyright © by Hoàng Ngọc Hà  
Bài thực hành môn học Mô hình hóa, mô phỏng và tối ưu hóa các quá trình hóa học  
c) Tìm nghiệm của (9) với Matlab dùng lệnh ode. Biểu diễn sự thay đổi của nhiệt độ  
nhận được tại các điểm rời rạc dùng lệnh plot.  
d) Viết biểu thức nghiệm T(x,t) dùng nội suy kiểu đa thức Lagrange từ các nút giá trị T  
it  
đã có. Biểu diễn sự thay đổi của biên dạng (profile) nhiệt độ với Matlab.  
Câu hỏi 3 (Câu hỏi mở rộng)  
1. Giả sử rằng biểu thức của năng lượng nội năng U(t) của thanh kim loại được  
cho bởi :  
L
U(t) = πR2 u(x,t)dx  
(10)  
∫
0
với u(x,t) = ρcT là mật độ thể tích năng lượng. Dùng qui tắc hình thang, chứng  
minh rằng (13) có thể xấp xỉ như sau :  
N−1 (Tit + T(i+1)t  
)
U(t)  
ρc  
= πR2h  
∑
2
i=0  
(11)  
Lập trình tính giá trị số của (11) với giá trị các T đã tìm thấy.  
it  
2. Quan sát ảnh hưởng của số điểm nút (tăng/giảm N) trên các nghiệm số bằng  
mô phỏng ? Kết luận.  
Copyright © by Hoàng Ngọc Hà  
pdf 4 trang baolam 29/04/2022 4780
Download
Bạn đang xem tài liệu "Bài thực hành môn Mô hình hóa, mô phỏng và tối ưu hóa các quá trình hóa học - Bài thực hành số 4: Khảo sát một quá trình truyền nhiệt", để tải tài liệu gốc về máy hãy click vào nút Download ở trên

File đính kèm:

  • pdfbai_thuc_hanh_mon_mo_hinh_hoa_mo_phong_va_toi_uu_hoa_cac_qua.pdf