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ố 3: Khảo sát động học một hệ phản ứng dạng ống lý tưởng
lượt xem 3
download
Nội dung của tài liệu bao gồm: rời rạc hóa một hệ thống phản ứng đượ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; 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ình luận(0) Đăng nhập để gửi bình luận!
Nội dung Text: 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ố 3: Khảo sát động học một hệ phản ứng dạng ống lý tưởng
- 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Ố 3 KHẢO SÁT ĐỘNG HỌC MỘT HỆ PHẢN ỨNG DẠNG ỐNG LÝ TƯỞNG 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. Rời rạc hóa một hệ thống phản ứng đượ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ố. Tài liệu tham khảo của bài thực hành: [1] Martin Ruszkowski et al., Passivity based control of transport reaction systems, AIChE Journal, 2005. 1. Mô tả hệ thống phản ứng ống: Hệ khảo sát trong bài thực hành này là một hệ phản ứng dạng ống (tubular reactor) lý tưởng phân bố, được mô tả như trong hình dưới đây: Dòng A vào Dòng ra R AÆB 0 L x Hệ thống trên được giới hạn diễn ra ở các điều kiện đẳng nhiệt, T=const., và sự thay đổi của các đại lượng vật lý (nồng độ,…) chỉ theo phương truc x, bỏ qua sự thay đổi theo phương bán kính. Hệ xem xét được cung cấp ở lối vào bằng chất phản ứng A. Bên trong ống xảy ra phản ứng bậc 1 dạng AÆB với tốc độ phản ứng : σ = −kc (1) vói k là hằng số và c là nồng độ (cục bộ) của cấu tử hóa học A. Câu hỏi 1 (Mô hình hóa động học) : Nghiên cứu cân bằng vật chất của hệ đẳng nhiệt trên, chỉ ra rằng sự biến thiên của nông độ c được chi phối bởi phương trình vi phân đạo hàm riêng sau [1]: ∂c ∂ 2c ∂c (2) = D 2 −υ + σ ∂t ∂x ∂x Copyright © by Hoàng Ngọc Hà CuuDuongThanCong.com https://fb.com/tailieudientucntt
- 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 Trong phương trình (2) t ∈ [0, + ∞ ) , x ∈ [0, L] ; υ và D lần lượt là vận tốc dòng đối lưu và hệ số khuếch tán. □ 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 kiện biên: • Với điều kiện ban đầu : c( x,0) = c 0 ( x) (3) • Với điều kiện biên Hulburt (ĐK biên Direchlet + ĐK biên Neumann): ⎧c(0, t ) = c0 (t ) ⎪ (4) ⎨ ∂c( x, t ) ⎪ ∂x =0 ⎩ x= L Phương trình (2) với (3) và (4) miêu tả đầy đủ hệ thống. 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 : υ (cm/s) D (cm2/s) L (cm) k (1/s) R (cm) c0 (t ) (mol/cm3) c 0 (t ) (mol/cm3) (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 (2) 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ằng h = . Chúng ta ký hiệu xi là các điểm rời rạc : N x0 = 0, x1 = h, x2 = 2h,..., x N = L Các điểm này phân bố hình học như sau : x0 = 0 x1 = h x2 = 2h ... xN = L x x x123x x x x N +1 h Copyright © by Hoàng Ngọc Hà CuuDuongThanCong.com https://fb.com/tailieudientucntt
- 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 c0 = c(0), c1 = c(h),..., c N = c( Nh) . Tại các điểm trung gian x1 , x2 ,..., x N −1 chúng 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 : ∂ 2c ci +1 + ci −1 − 2ci ≅ , i = 1,..., N − 1 (5) ∂x 2 x = xi h2 ∂c ci +1 − ci −1 ≅ , i = 1,..., N − 1 (6) ∂x x = xi 2h Chúng ta đưa vào tập điểm rời rạc trên một điểm ảo x N +1 (và ký hiệu c N +1 = c(( N + 1)h) , xem sơ đồ phân bố các điểm bên trên) cho mục đích rời rạc hóa tại điểm biên x N , cụ thể : ∂ 2c c N +1 + c N −1 − 2c N ≅ (7) ∂x 2 x = xN h2 ∂c c N +1 − c N −1 ≅ (8) ∂x x = xN 2h Dùng quan hệ thứ nhất trong (4), (8) được xấp xỉ thành: ∂c c N +1 − c N −1 ≅ = 0 ⇔ c N +1 = c N −1 (9) ∂x x= xN 2h Sau cùng, để thuân lợi cho biểu diễn chúng ta ký hiệu : ⎡c(h, t ) ⎤ ⎡c1t ⎤ ⎢ c ( 2 h, t ) ⎥ ⎢ c ⎥ C (t ) = ⎢ ⎥ = ⎢ 2t ⎥ ⎢M ⎥ ⎢M ⎥ (10) ⎢ ⎥ ⎢ ⎥ ⎣c( Nh, t )⎦ ⎣c Nt ⎦ 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)Æ(6), chỉ ra rằng phương trình (4) tại điểm rời rạc i trở thành : dcit ⎛ D υ ⎞ ⎛ 2D ⎞ ⎛D υ ⎞ (11) = ⎜ 2 + ⎟c(i −1)t − ⎜ 2 + k ⎟cit + ⎜ 2 − ⎟c(i +1) t , i = 1,..., N dt ⎝ h 2h ⎠ ⎝h ⎠ ⎝h 2h ⎠ Copyright © by Hoàng Ngọc Hà CuuDuongThanCong.com https://fb.com/tailieudientucntt
- 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) Tính đến các điều kiên biên (4) và (9), chỉ ra rằng phương trình (11) có thể viết dưới dạng hệ phương trình vi phân thường sau : dC = AC + b dt (12) với C (t ) là vector được cho trong (10) và giá trị ban đầu C (t = 0) có được dùng (3). Xác định các ma trận A, b ? c) Tìm nghiệm của (12) với Matlab dùng lệnh ode. Biểu diễn biên dạng (profile) của nồng độ nhận được tại các điểm rời rạc. d) Viết biểu thức nghiệm c(x,t) dùng nội suy kiểu đa thức Lagrange từ các nút cit (12) đã có. Biểu diễn sự thay đổi của nồng độ với Matlab. Câu hỏi 3 (Câu hỏi mở rộng) 1. Biết rằng biểu thức số mol của A trong ống phản ứng được cho bởi : L N (t ) = πR ∫ c( x, t )dx 2 (13) 0 Dùng qui tắc hình thang, chứng minh rằng (13) có thể xấp xỉ như sau : N −1 (cit + c( i +1)t ) N (t ) = πR 2 h∑ i =0 2 (14) Lập trình tính giá trị số của (14) với giá trị các cit đã tìm thấy trong (12). 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à CuuDuongThanCong.com https://fb.com/tailieudientucntt
CÓ THỂ BẠN MUỐN DOWNLOAD
-
Bài giảng môn học Công nghệ sản xuất protein, axit amin và axit hữu cơ - Trương Thị Minh Hạnh (biên soạn)
66 p | 391 | 75
-
Luyện tập bài tập thực hành mô hình hóa môi trường
87 p | 319 | 56
-
Bài giảng môn lý thuyết ôtômát và ngôn ngữ hình thức - Chương 4
0 p | 247 | 54
-
CÁC BÀI THỰC HÀNH MÔN HỌC TÍN HIỆU VÀ HỆ THỐNG
43 p | 342 | 37
-
Bài giảng môn lý thuyết ôtômát và ngôn ngữ hình thức - Chương 7
0 p | 191 | 30
-
Sử dụng hệ thống bài tập thực hành thí nghiệm phần phi kim trong dạy học hóa học nhằm phát triển năng lực thực nghiệm cho học sinh lớp 11
6 p | 108 | 10
-
Một số ví dụ về Mô hình phân tích hành vi của doanh nghiệp
17 p | 138 | 9
-
Vận dụng mô hình 5E trong dạy học thực hành thí nghiệm môn Khoa học Tự nhiên ở trường trung học cơ sở
10 p | 125 | 9
-
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: Làm quen với phần mềm Matlab
9 p | 44 | 6
-
Sử dụng bài toán thực tiễn trong dạy học chủ đề tích phân lớp 12
6 p | 38 | 5
-
Đề xuất cấu trúc và đánh giá năng lực thực nghiệm cho học sinh thông qua môn khoa học tự nhiên cấp trung học cơ sở
10 p | 120 | 5
-
Một nghiên cứu về hoạt động thực hành nghiệp vụ của giáo viên Toán
6 p | 47 | 4
-
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ố 1: Khảo sát động học của một bình phản ứng lý tưởng khuấy liên tục
5 p | 66 | 3
-
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
5 p | 43 | 3
-
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ố 2: Khảo sát một hệ phản ứng dạng ống lý tưởng trạng thái dừng
5 p | 47 | 2
-
Thiết kế và tổ chức hoạt động thực hành thí nghiệm trong học tập trải nghiệm để dạy học phần Sinh học cơ thể trong chương trình trung học phổ thông nhằm phát triển năng lực sáng tạo cho học sinh
14 p | 23 | 2
-
Bài giảng Quản lý lưu vực (Watershed Management): Giới thiệu chương trình học - ThS. Nguyễn Duy Liêm
8 p | 5 | 1
Chịu trách nhiệm nội dung:
Nguyễn Công Hà - Giám đốc Công ty TNHH TÀI LIỆU TRỰC TUYẾN VI NA
LIÊN HỆ
Địa chỉ: P402, 54A Nơ Trang Long, Phường 14, Q.Bình Thạnh, TP.HCM
Hotline: 093 303 0098
Email: support@tailieu.vn