Tổng quan nghiên cứu

Phản ứng hợp nước xúc tác axit pha lỏng từ 2,3-epoxy-1-propanol (glycidol) tạo glycerol là một quá trình nhiệt động tỏa lượng nhiệt đặc biệt lớn với biến thiên entanpi đạt mức âm 87,7 nhân 10 mũ 3 J/mol (-87,7 kJ/mol). Trong thực tế sản xuất công nghiệp, glycidol là hợp chất kém bền nhiệt với điểm sôi ở mức 440,15 K và có tính nguy hại cao, trong khi sản phẩm glycerol lại có giá trị kinh tế vượt trội và thân thiện với môi trường. Tuy nhiên, việc vận hành thiết bị phản ứng khuấy trộn liên tục (CSTR) cho hệ phản ứng này thường xuyên đối mặt với hiện tượng đa trạng thái dừng (multiple steady states), khiến hệ thống rất dễ rơi vào trạng thái vận hành bất ổn định hoặc xảy ra sự cố tăng nhiệt mất kiểm soát.

Vấn đề cốt lõi đặt ra là các mô hình động lực học truyền thống xây dựng từ phương trình vi phân thường (ODEs) dựa trên cân bằng vật chất và năng lượng chỉ mô tả được sự biến thiên nồng độ và nhiệt độ, mà chưa thể biểu diễn tường minh bản chất phân tán năng lượng cũng như mức độ bất thuận nghịch sinh ra trong hệ thống. Luận văn tập trung vào mục tiêu thiết lập mô hình Hamilton dựa trên cổng (Port-Controlled Hamiltonian) có xét đến thành phần phân tán năng lượng, từ đó định tính và định lượng hóa nguyên nhân gây tiêu tán năng lượng của quá trình.

Nghiên cứu được triển khai thực hiện trong thời gian 6 tháng (từ ngày 04/07/2016 đến ngày 04/12/2016) tại Phòng thí nghiệm Trọng điểm Điều khiển số và Kỹ thuật Hệ thống (DCSELAB), Trường Đại học Bách Khoa - ĐHQG TP. Hồ Chí Minh. Ý nghĩa then chốt của công trình là cung cấp cơ sở toán học chuẩn xác giúp giải mã 3 trạng thái cân bằng dừng của thiết bị CSTR, hỗ trợ việc thiết kế các thuật toán điều khiển phi tuyến nhằm tối ưu hóa hiệu suất năng lượng, kiểm soát nhiệt độ thấp hơn ngưỡng sôi tới hạn 440,15 K và giảm thiểu mức tổn hao năng lượng bất thuận nghịch từ 15% đến 25% trong quá trình vận hành thực tế.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên sự giao thoa của 2 khung lý thuyết nền tảng: lý thuyết truyền vận kinh điển và lý thuyết hệ thống Hamilton dựa trên cổng mở rộng cho nhiệt động lực học không cân bằng.

Lý thuyết truyền vận cung cấp 3 phương trình cốt lõi để mô tả trường thủy động lực học, trường nhiệt độ và trường nồng độ trong không gian phản ứng: phương trình Navier-Stokes mô tả bảo toàn xung lượng của lưu chất không chịu nén, phương trình vi phân dẫn nhiệt đối lưu Fourier-Kirchhoff mô tả trường nhiệt độ, và phương trình Fick II mô tả quy luật phân bố nồng độ cấu tử khuếch tán.

Lý thuyết hệ thống Hamilton dựa trên cổng (PCHS) phát triển bởi Maschke và Van Der Schaft được ứng dụng để cấu trúc hóa động lực học của hệ thống qua hàm lưu trữ Hamilton H(x), ma trận kết nối phản đối xứng J(x) thỏa mãn J(x) = -J(x)^T và ma trận phân tán đối xứng bán xác định dương R(x) thỏa mãn R(x) = R(x)^T >= 0. Vì nội năng tổng cộng luôn được bảo toàn theo nguyên lý 1 nhiệt động học, hàm lưu trữ Hamilton trong hệ hóa học được lựa chọn là đại lượng âm của entropy (-S) hoặc bình phương của ái lực hóa học (A^2). Các khái niệm chính bao gồm: biến trạng thái quảng tính (nội năng, số mol cấu tử), biến cường tính (nhiệt độ, hóa thế), và tốc độ sinh entropy bất thuận nghịch chia làm 3 thành phần riêng biệt gồm khuấy trộn, trao đổi nhiệt qua vỏ áo và phản ứng hóa học.

Phương pháp nghiên cứu

Nguồn dữ liệu của nghiên cứu xuất phát từ các thông số hóa lý thực nghiệm của phản ứng hợp nước glycidol, bao gồm năng lượng hoạt hóa Ea, hằng số tốc độ phản ứng bậc một k0, nhiệt phản ứng tỏa nhiệt -ΔH = 87.700 J/mol, và nhiệt dung riêng Cp của hỗn hợp phản ứng.

Phương pháp phân tích chủ đạo là mô hình hóa giải tích lý hóa kết hợp phân tích hệ thống. Hệ phương trình vi phân phi tuyến bậc 4 được thiết lập từ phương trình cân bằng vật chất cho 4 cấu tử (A: 2,3-epoxy-1-propanol, B: H2O, C: glycerol, D: xúc tác H2SO4) cùng phương trình cân bằng năng lượng không đẳng nhiệt. Lý do lựa chọn mô hình giải tích PCHS là nhằm liên kết trực tiếp các quy luật bảo toàn với hàm suy giảm Lyapunov, giúp phân tích ổn định hệ thống một cách chặt chẽ. Để khảo sát đặc tính tĩnh, phương pháp giản đồ Van Heerden được áp dụng nhằm xác định các giao điểm giữa đường cong sinh nhiệt phi tuyến dạng Arrhenius và đường thẳng tiêu thụ nhiệt qua áo làm mát.

Toàn bộ thuật toán mô hình hóa được cài đặt và kiểm chứng số thông qua phần mềm MATLAB & SIMULINK trong suốt chu kỳ nghiên cứu 6 tháng, sử dụng các khối giải tích tích phân thời gian thực và trích xuất dữ liệu để đánh giá tốc độ sinh entropy tại từng chế độ vận hành.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, nghiên cứu đã thiết lập thành công mô hình toán học dạng Hamilton phân tán năng lượng cho thiết bị phản ứng CSTR với 4 cấu tử, chứng minh một cách tường minh rằng tốc độ biến thiên hàm Hamilton dH/dt luôn bị triệt tiêu một phần năng lượng bởi đại lượng phân tán d = (dH/dx)^T * R(x) * (dH/dx) >= 0.

Thứ hai, khảo sát đặc tính tĩnh thông qua giản đồ Van Heerden đã phát hiện chính xác 3 trạng thái cân bằng dừng riêng biệt (ký hiệu là P1, P2, P3) của hệ thống phản ứng CSTR:

  • Trạng thái dừng P1: Đặc trưng cho vùng nhiệt độ thấp, tốc độ phản ứng chậm với độ chuyển hóa của cấu tử A chỉ đạt dưới mức 15%, nhiệt lượng sinh ra không đáng kể.
  • Trạng thái dừng P2: Tương ứng với điểm nhiệt độ trung gian (điểm bốc cháy), là trạng thái cân bằng dừng không ổn định của hệ thống động học.
  • Trạng thái dừng P3: Đặc trưng cho vùng nhiệt độ cao, tốc độ phản ứng diễn ra cực mạnh với độ chuyển hóa cấu tử A đạt trên 95%, giải phóng nhiệt lượng cực đại.

Thứ ba, kết quả mô phỏng trên MATLAB & SIMULINK tại 4 bộ điều kiện ban đầu (C1, C2, C3, C4) cho thấy tốc độ sản sinh entropy bất thuận nghịch do phản ứng hóa học chiếm tỷ trọng áp đảo, lên tới hơn 70% tổng mức năng lượng phân tán của toàn hệ thống khi phản ứng vận hành ở vùng nhiệt độ cao, trong khi tổn thất do truyền nhiệt qua vỏ áo chiếm khoảng 20% đến 25% và do khuấy trộn cơ học chiếm dưới 5%.

Thảo luận kết quả

Nguyên nhân căn bản dẫn đến sự phân tán năng lượng và tính bất thuận nghịch trong thiết bị CSTR bắt nguồn từ hiệu ứng tỏa nhiệt mãnh liệt (-ΔH = 87,7 kJ/mol) kết hợp cùng độ chênh lệch nhiệt độ lớn giữa khối chất lỏng phản ứng với lưu chất làm mát trong vỏ áo có hệ số trao đổi nhiệt xác định.

Khi so sánh với các hệ thống vật lý kinh điển, sự phân tán trong mạch điện RLC đơn thuần do điện trở R tiêu tán dưới dạng nhiệt lượng tỏa ra môi trường với công suất Ri^2, và trong hệ cơ học lò xo - giảm chấn là do lực ma sát c*(dx/dt)^2. Ngược lại, trong hệ quá trình hóa học, sự phân tán năng lượng mang bản chất nhiệt động học sâu sắc, gắn liền với định luật hai nhiệt động học và sự tăng trưởng liên tục của entropy nội tại.

Các dữ liệu nghiên cứu có thể được trực quan hóa tối ưu thông qua giản đồ Van Heerden, nơi giao điểm giữa đường cong hàm mũ sinh nhiệt Arrhenius và đường thẳng tải nhiệt tuyến tính theo hệ số góc Cp làm nổi bật rõ ranh giới an toàn nhiệt độ. Đồng thời, biểu đồ mô phỏng theo trục thời gian thực từ SIMULINK cho thấy rõ xu hướng tiến tới cân bằng ổn định của tốc độ sinh entropy tại các điểm P1 và P3, minh chứng rằng nếu không có sự can thiệp của bộ điều khiển, hệ thống sẽ tự động trượt khỏi điểm cân bằng tối ưu P2 và có nguy cơ vượt qua ngưỡng nhiệt độ sôi an toàn 440,15 K.

Đề xuất và khuyến nghị

Thứ nhất, kỹ sư thiết kế và vận hành quy trình công nghệ cần áp dụng ngay phương pháp điều khiển thụ động gán kết nối và giảm chấn (IDA-PBC) dựa trên biểu diễn Hamilton, nhằm ổn định hóa hệ thống CSTR tại điểm làm việc tối ưu P2 trong dải nhiệt độ an toàn từ 360 K đến 390 K, giúp loại bỏ hoàn toàn nguy cơ quá nhiệt trong thời gian triển khai từ 3 đến 6 tháng.

Thứ hai, bộ phận kỹ thuật nhà máy sản xuất hóa chất cần điều chỉnh tối ưu lưu lượng và nhiệt độ dòng làm mát vỏ áo, hướng tới mục tiêu giảm thiểu từ 18% đến 22% tốc độ sinh entropy do truyền nhiệt, nâng cao hiệu quả thu hồi năng lượng trong 6 tháng đầu năm vận hành.

Thứ ba, các nhóm nghiên cứu tại các trường đại học và viện chuyên ngành tự động hóa nên mở rộng thuật toán mô hình Hamilton phân tán năng lượng này cho các hệ phản ứng dị thể đa pha hoặc thiết bị phản ứng dạng ống (PFR), hoàn thiện bộ thư viện khối chức năng mô phỏng chuyên dụng trong thời hạn 12 tháng tới.

Thứ tư, ban quản lý an toàn công nghệ hóa chất cần xây dựng quy chuẩn vận hành nghiêm ngặt và lắp đặt hệ thống cảm biến nhiệt độ đa điểm có thời gian đáp ứng dưới 0,5 giây tại các thiết bị phản ứng tỏa nhiệt cao, nhằm ngăn chặn kịp thời hiện tượng quá nhiệt đột ngột (runaway reaction) trong quý tới.

Đối tượng nên tham khảo luận văn

Nhóm 1: Kỹ sư công nghệ hóa chất và thiết kế quá trình sản xuất. Luận văn cung cấp toàn bộ mô hình toán học giải tích cân bằng vật chất và năng lượng của thiết bị CSTR, giúp kỹ sư ứng dụng trực tiếp vào việc tính toán kích thước bồn phản ứng, thiết kế bề mặt truyền nhiệt vỏ áo và xác định chính xác các điểm cân bằng dừng cho các phản ứng tỏa nhiệt cao trên 80 kJ/mol.

Nhóm 2: Chuyên gia điều khiển tự động và kỹ sư kỹ thuật hệ thống. Luận văn mở ra phương pháp tiếp cận hiện đại thông qua mô hình Hamilton dựa trên cổng (PCHS), hỗ trợ thiết lập hàm Lyapunov phi tuyến chuẩn xác để thiết kế các bộ điều khiển phi tuyến nâng cao cho các hệ thống hóa học có động học phức tạp.

Nhóm 3: Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Hóa học và Tự động hóa. Tài liệu là nguồn tham khảo học thuật giá trị cao về cách kết hợp giữa nhiệt động lực học không cân bằng và lý thuyết điều khiển hệ thống, đi kèm hướng dẫn chi tiết về kỹ thuật mô phỏng trên MATLAB & SIMULINK.

Nhóm 4: Cán bộ quản lý an toàn công nghệ và kiểm toán năng lượng tại các nhà máy hóa chất, dược phẩm. Luận văn cung cấp phương pháp định lượng sự phân tán năng lượng và tốc độ sinh entropy, giúp xác định các điểm tổn thất năng lượng chính và thiết lập ranh giới an toàn cho hệ thống phản ứng.

Câu hỏi thường gặp

Tại sao nội năng không thể được chọn làm hàm lưu trữ Hamilton cho hệ phản ứng hóa học? Theo nguyên lý 1 của nhiệt động lực học, tổng năng lượng trong vũ trụ luôn được bảo toàn, do đó nội năng không thể tự suy giảm trong suốt quá trình biến đổi nội tại. Để phản ánh tính bất thuận nghịch tự nhiên và sự phân tán năng lượng theo nguyên lý 2 nhiệt động học, hàm lưu trữ Hamilton bắt buộc phải được chọn từ đại lượng âm của entropy (-S) hoặc bình phương ái lực hóa học (A^2).

Hiện tượng đa trạng thái dừng trong thiết bị CSTR đe dọa an toàn vận hành như thế nào? Giản đồ Van Heerden chỉ ra sự tồn tại của 3 điểm cân bằng P1, P2, P3. Trong đó, điểm P2 là trạng thái không ổn định; chỉ cần một dao động nhiệt nhỏ, tốc độ sinh nhiệt sẽ vượt qua tốc độ làm mát, đẩy hệ thống chuyển dịch nhanh chóng về trạng thái bốc cháy P3 hoặc vượt quá ngưỡng sôi 440,15 K của glycidol, gây nguy cơ quá áp và nổ thiết bị.

Ma trận phân tán R(x) trong biểu diễn Hamilton thể hiện bản chất vật lý gì? Ma trận R(x) là ma trận đối xứng bán xác định dương, lượng hóa toàn bộ các yếu tố gây mất mát năng lượng hữu ích trong hệ thống. Trong thiết bị CSTR, ma trận này đặc trưng trực tiếp cho trở lực ma sát dòng chảy khuấy trộn, sự truyền nhiệt bất thuận nghịch qua vỏ áo và tốc độ sản sinh entropy do phản ứng hóa học tạo glycerol.

Phần mềm MATLAB & SIMULINK đóng vai trò gì trong việc giải quyết bài toán của luận văn? MATLAB & SIMULINK cung cấp môi trường giải tích số mạnh mẽ để giải hệ phương trình vi phân phi tuyến bậc 4 mô tả 4 cấu tử phản ứng. Thông qua các khối chức năng như Integrator, Math Operations và Scope, phần mềm cho phép mô phỏng chính xác động học phản ứng và vẽ đồ thị tốc độ sinh entropy tại 4 điều kiện đầu (C1 đến C4) theo thời gian thực.

Mô hình Hamilton phân tán năng lượng có thể mở rộng cho các phản ứng hóa học khác không? Mô hình hoàn toàn có thể áp dụng cho mọi hệ phản ứng hóa học khác trong công nghiệp. Do được xây dựng từ các phương trình truyền vận tổng quát và phương trình năng lượng Gibbs, cấu trúc không gian trạng thái Hamilton và ma trận phân tán R(x) có thể tùy biến dễ dàng cho các phản ứng dị thể, phản ứng nối tiếp hoặc phản ứng xúc tác đa pha.

Kết luận

  • Luận văn đã thiết lập thành công mô hình toán học giải tích mô tả chính xác động học phi tuyến của thiết bị phản ứng CSTR cho phản ứng thủy phân 2,3-epoxy-1-propanol tạo glycerol với 4 cấu tử tham gia.
  • Ứng dụng thành công biểu diễn Hamilton dựa trên cổng (PCHS) với hàm lưu trữ âm entropy, chỉ ra tường minh nguồn gốc phân tán năng lượng thông qua ma trận cấu trúc R(x) bán xác định dương.
  • Xác định và phân tích thấu đáo 3 trạng thái cân bằng dừng P1, P2, P3 trên giản đồ Van Heerden, làm sáng tỏ cơ chế bất ổn định nhiệt động học tại điểm vận hành tối ưu P2.
  • Định lượng chính xác tốc độ sản sinh entropy bất thuận nghịch bằng công cụ mô phỏng số MATLAB & SIMULINK tại 4 bộ điều kiện đầu, khẳng định phản ứng hóa học tỏa nhiệt (-ΔH = 87,7 kJ/mol) là nguyên nhân gây phân tán năng lượng chủ đạo.
  • Mở ra nền tảng lý thuyết vững chắc cho việc thiết kế các thuật toán điều khiển phi tuyến thụ động (PBC) nhằm đảm bảo an toàn tuyệt đối và nâng cao hiệu quả năng lượng cho các thiết bị phản ứng công nghiệp.

Trong giai đoạn 6 đến 12 tháng tới, các kỹ sư và nhà nghiên cứu nên ứng dụng ngay mô hình này để hoàn thiện thuật toán điều khiển ổn định hóa vòng kín, đồng thời mở rộng thử nghiệm trên các hệ thống phản ứng đa cấu tử quy mô công nghiệp thực tế.