Nội dung của bài thực hành 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; tìm nghiệm số của hệ ODE này. Xấp xỉ nghiệm của bài toán dùng phương pháp nội suy Lagrange; kết luận về ảnh hưởng của số điểm nút chọn trên các nghiệm số. | 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 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 độ T 1 có thể thay đổi theo thời gian. 0 1 T T ( x, t ) T (t ) R 0 x x+dx 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 ) ∂ 2T ( x , t ) ρc =λ ∂t ∂x 2 (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 ) ∂ 2T ( x, t ) =D ∂t ∂x 2 (2) Copyright © by Hoàng Ngọc Hà .
