Luận án tiến sĩ: Mô hóa kênh xúc tác và lò phản ứng monolith

Hệ phương trình cơ bản bao gồm cân bằng khối lượng cho pha khí và pha rắn. Phương trình bảo toàn năng lượng tính đến dẫn nhiệt trục và bức xạ nhiệt. Mô hình sử dụng giả thiết chuẩn ổn định (quasi-steady) cho pha khí. Độ dày lớp biên động lượng và nhiệt được xác định qua số Reynolds. Cân bằng vật liệu và năng lượng áp dụng cho từng thể tích hữu hạn dọc theo kênh.

Cơ chế phản ứng dị thể trên bề mặt platinum bao gồm hấp phụ và giải hấp các chất. Tốc độ phản ứng phụ thuộc vào phân số vị trí bề mặt của CO(s) và O(s). Mô hình sử dụng cơ chế Deutschmann cho phản ứng CO/O2 khô và ẩm. Điều kiện ban đầu bề mặt thường giả định bao phủ hoàn toàn bởi oxy. Tích số ZCO(s)·ZO(s) tỷ lệ thuận với tốc độ chuyển hóa nhiên liệu.

Rời rạc hóa không gian sử dụng phương pháp thể tích hữu hạn. Tích phân theo thời gian áp dụng cho pha rắn với dung sai lỗi tuyệt đối và tương đối. Thuật toán giải lặp qua các bước thời gian cho đến đạt trạng thái ổn định. Sơ đồ giải bao gồm cập nhật tính chất vận chuyển và kiểm tra hội tụ. Chương trình máy tính tối ưu hóa hiệu suất tính toán cho hệ thống lớn.

Dẫn nhiệt trục dọc theo thành kênh platinum ảnh hưởng đến phân bố nhiệt độ. Số Peclet nhiệt rắn đánh giá tầm quan trọng của quá trình này. Với giá trị Peclet lớn, dẫn nhiệt trục có thể bỏ qua. Tuy nhiên, trong vùng phản ứng mạnh, gradient nhiệt độ cao làm tăng vai trò dẫn nhiệt. Mô hình đầy đủ giữ lại số hạng này để chính xác hơn.

Trao đổi nhiệt giữa dòng khí và thành kênh xảy ra qua lớp biên nhiệt. Số Nusselt đặc trưng cho hiệu quả trao đổi nhiệt này. Giá trị Nu=4.36 thường dùng cho kênh tròn phát triển đầy đủ. Nghiên cứu kiểm tra ảnh hưởng của việc thay đổi Nu đến kết quả mô phỏng. Hệ số truyền nhiệt tính từ Nu, độ dẫn nhiệt khí và đường kính thủy lực.

Bức xạ nhiệt trở nên quan trọng ở nhiệt độ cao trên 800K. Mô hình tính đến trao đổi bức xạ giữa các phần tử bề mặt kênh. Độ phát xạ của platinum ảnh hưởng đến cường độ bức xạ. Tại vùng đầu vào lạnh, bức xạ giúp làm ấm pha rắn. Ở vùng phản ứng nóng, bức xạ hỗ trợ phân tán nhiệt dọc trục.

Số Sherwood định lượng hiệu quả truyền khối đối lưu so với khuếch tán. Giá trị Sh tính từ tương tự nhiệt-khối với điều chỉnh theo số Schmidt. Hệ số truyền khối liên hệ trực tiếp với Sh, hệ số khuếch tán và kích thước đặc trưng. Nghiên cứu thay đổi Sh để kiểm tra độ nhạy của mô hình. Kết quả cho thấy Sh ảnh hưởng mạnh đến chuyển hóa nhiên liệu.

Lớp biên nồng độ hình thành gần thành kênh do phản ứng bề mặt. Độ dày lớp biên phụ thuộc vào vận tốc dòng chảy và hệ số khuếch tán. Phân số khối lượng tại thành (YkW) khác với giá trị dòng chính. Gradient này tạo ra lực kéo khuếch tán vận chuyển chất phản ứng. Mô hình tính toán YkW cho từng loài tại mỗi vị trí trục.

Tỷ lệ đương lượng tại thành (φW) khác với giá trị dòng chính do truyền khối chọn lọc. Các loài có hệ số khuếch tán khác nhau di chuyển với tốc độ không đồng đều. Điều này tạo ra sự khác biệt về thành phần hỗn hợp tại bề mặt. Giá trị φW ảnh hưởng đến động học phản ứng bề mặt. Mô hình theo dõi φW dọc theo chiều dài kênh để đánh giá hiệu ứng này.

Thực nghiệm Khitrin và Solovyeva sử dụng ống platinum ở ba vận tốc khác nhau. Mô hình tính toán với Nu=4.36 cho dòng phát triển đầy đủ. Kết quả cho thấy sự phù hợp tốt về chuyển hóa CO dọc theo kênh. Tại vận tốc 34 m/s, mô hình dự đoán chính xác vị trí hoàn thành phản ứng. Phân tích cho thấy sự quan trọng của truyền khối đến bề mặt xúc tác.

Cấu trúc monolith gồm nhiều kênh song song với diện tích bề mặt lớn. Mô hình áp dụng cho một kênh đại diện với điều kiện biên đối xứng. Trao đổi nhiệt giữa các kênh qua thành ngăn được xem xét. Kết quả nhiệt độ khớp với đo đạc hồng ngoại dọc theo monolith. Mô hình nắm bắt được điểm bắt lửa và vùng phản ứng mạnh.

Thí nghiệm ống platinum ngang quan sát lan truyền mặt trận phản ứng. Mô hình tạm thời theo dõi sự phát triển nhiệt độ theo thời gian. Tốc độ lan truyền sóng phụ thuộc vào cân bằng nhiệt khí-rắn. Kết quả cho thấy sự di chuyển vùng nóng dọc theo kênh. Dẫn nhiệt trục trong pha rắn ảnh hưởng đến tốc độ lan truyền.

Thang thời gian đối lưu khí tính từ chiều dài kênh và vận tốc trung bình. Giá trị điển hình trong khoảng mili-giây cho kênh cỡ cm. Thang thời gian khuếch tán khí dựa trên hệ số khuếch tán và kích thước đặc trưng. Tỷ số giữa hai thang thời gian này cho số Peclet khối. Kết quả cho thấy đối lưu chiếm ưu thế trong vận chuyển khí.

Thang thời gian khuếch tán nhiệt rắn phụ thuộc vào độ dày thành và độ khuếch tán nhiệt. Giá trị thường lớn hơn thang thời gian khí hàng trăm lần. Điều này tạo ra sự tách biệt rõ ràng giữa động lực học hai pha. Tỷ số nhiệt dung thể tích khí-rắn thường nhỏ hơn 0.01. Pha rắn hoạt động như bể chứa nhiệt lớn cho hệ thống.

Số Peclet nhiệt khí đánh giá tầm quan trọng dẫn nhiệt trục khí. Giá trị lớn cho phép bỏ qua dẫn nhiệt trục trong pha khí. Số Peclet nhiệt rắn nhỏ hơn, đặc biệt trong vùng gradient cao. Điều này yêu cầu giữ lại số hạng dẫn nhiệt trục rắn. Số Peclet khối tương tự đánh giá khuếch tán trục các loài.

Mở rộng sang mô hình 2D hoặc 3D để nắm bắt biến đổi hướng kính. Tích hợp hiệu ứng bức xạ phức tạp hơn giữa các kênh lân cận. Xem xét sự không đồng nhất của lớp phủ xúc tác. Bao gồm ảnh hưởng của lão hóa xúc tác và nhiễm độc. Mô hình hóa chi tiết hơn vùng đầu vào với dòng chảy phát triển.

Đo đạc chi tiết hơn về phân bố nhiệt độ bề mặt. Phân tích thành phần khí tại nhiều vị trí dọc kênh. Kiểm soát chính xác hơn điều kiện biên nhiệt độ và lưu lượng. Đặc tính hóa đầy đủ tính chất vật liệu xúc tác. Thí nghiệm với nhiều loại nhiên liệu và điều kiện vận hành khác nhau.

Cải thiện thuật toán giải để giảm thời gian tính toán. Áp dụng kỹ thuật song song hóa cho hệ thống đa kênh. Tối ưu hóa lưới tính toán dựa trên gradient địa phương. Sử dụng phương pháp bậc cao hơn để tăng độ chính xác. Phát triển công cụ hậu xử lý tự động cho phân tích kết quả.