Microsoft Word 00 a loinoidau TV (moi thang1 2016) docx 42 Mai Thị Đoan Thanh, Nguyễn Đình Lâm, Đoàn Quang Vinh XÂY DỰNG MÔ HÌNH THỬ NGHIỆM THIẾT BỊ PHẢN ỨNG LIÊN TỤC THỦY PHÂN ACETIC ANHYDRIDE PHỤC V[.]
42 Mai Thị Đoan Thanh, Nguyễn Đình Lâm, Đồn Quang Vinh XÂY DỰNG MƠ HÌNH THỬ NGHIỆM THIẾT BỊ PHẢN ỨNG LIÊN TỤC THỦY PHÂN ACETIC ANHYDRIDE PHỤC VỤ NGHIÊN CỨU ĐIỀU KHIỂN AN EXPERIMENTAL MODEL OF CONTINUOUS STIRED TANK REACTOR WITH ACETIC ANHYDRIDE HYDROLYSIS FOR CONTROL RESEARCH Mai Thị Đoan Thanh1, Nguyễn Đình Lâm2, Đồn Quang Vinh3 Trường Cao đẳng Nghề Đà Nẵng; maithidoanthanh@gmail.com Trường Đại học Bách khoa, Đại học Đà Nẵng; ndlam@dut.udn.vn Đại học Đà Nẵng; dqvinh@ac.udn.vn Tóm tắt - Một phương pháp tạo acetic acid thủy phân hóa Acetic Anhydride bình phản ứng khuấy trộn liên tục (CSTR– Continuous Stirred Tank Reactor) Trong báo này, tác giả xây dựng mơ hình thiết bị phản ứng CSTR, với mục đích nghiên cứu động học điều khiển thiết bị phản ứng Xuất phát từ động học trình phản ứng, thơng số cơng nghệ phản ứng thủy phân hóa Acetic Anhydride, nội dung báo thiết kế mơ hình, tính tốn thơng số vận hành Mơ hình thiết kế kiểm tra qua mơ hai giai đoạn: giai đoạn khởi động giai đoạn vận hành, đánh giá đặc tính thiết bị Mơ hình điều khiển thiết bị phản ứng hệ đa biến, có đặc tính phi tuyến tác động xen kênh Với thiết kế điều khiển tuyến tinh phản hồi đầu dùng điều khiển PID, cho thấy chất lượng sản phẩm bị ảnh hưởng lớn nhiễu đầu vào, cần phải thiết kế bù nhiễu dùng điều khiển phi tuyến Abstract - One of the methods that generate acetic acid is Anhydride Acetic hydrolysis reactor in Continuous Stirred Tank Reactor (CSTR) For this reason, the article aims at modeling CSTR bioreactor with a view to studying the kinetics and control of the reactor The paper deals with designing an experimental model based on dynamics and technological parameters of Acetic Anhydride reaction The model is verified by numerical simulations via starting and operating phases, rating characteristics of the device The controlling CSTR model is multivariable, nonlinear and coupled The output feedback controller with PID action is greatly affected by noisy input; therefore there is a need to design noise offset or use nonlinear control Từ khóa - thiết bị phản ứng liên tục; điều khiển phi tuyến; điều khiển trình; cân bằng; phản ứng hydro hóa; Key words - Continuous Stirred Tank Reactor; nonlinear control; process control; balance; hydrolysis reaction; Đặt vấn đề Acetic Acid (hay gọi axit etanoic) chất hữu dạng lỏng, hợp chất hữu điển hình dãy đồng đẳng axit mono cacboxylic, có công thức CH3COOH, sản phẩm hữu cơ quan trọng sử dụng rộng rãi cơng nghệ tổng hợp hữu hóa dầu Ngày xưa người biết sử dụng Acetic Acid làm dấm ăn Hiện với phát triển khoa học kỹ thuật, Acetic Acid ứng dụng lĩnh vực khác công nghiệp như: sản xuất muối axetat, axit cloaxetic, hợp chất nhựa, dược phẩm, thuốc nhuộm thuốc trừ sâu, hóa chất ngành ảnh, thực phẩm phụ gia cho thực phẩm (ví dụ: dấm ăn, nước dầm hoa quả, nước sốt), ngồi dùng để điều chế nhiều loại hợp chất hữu khác Do có tầm quan trọng nên nhu cầu sử dụng Acetic Acid lớn, đòi hỏi ngành sản xuất Acetic Acid phải có bước phát triển với nhu cầu Người ta tạo Acetic Acid phản ứng thủy phân Acetic Anhydride thiết bị phản ứng khuấy trộn liên tục (CSTR) Trên Hình sơ đồ nguyên lý chung thiết bị phản ứng khuấy trộn liên tục Phần thiết bị phản ứng CSTR bao gồm: Bình phản ứng tích V chứa dung dịch phản ứng, khuấy trộn cánh khuấy quay động điện, dung dịch phản ứng cấp vào bình có lưu lượng F1, nồng độ CA1 nhiệt độ T1 sản phẩm hoàn thành lấy có lưu lượng F2, nồng độ CA2 nhiệt độ T, điều khiển lưu lượng dùng van bơm; Vỏ thiết bị được cấu tạo kiểu Jacket chứa môi chất gia nhiệt, nhiệt truyền qua thành bình phản ứng, điều chỉnh cơng suất nhiệt cấp cho phản ứng thơng qua lưu lượng mơi chất (có thể dùng van điều chỉnh bơm) Thiết bị phản ứng xảy bình có pha lỏng nên gọi phản ứng đồng thể Hình Thiết bị phản ứng khuấy trộn liên tục CSTR Nhằm tạo sản phẩm đầu có nồng độ đảm bảo chất lượng suất theo thiết kế, ta cần phải điều khiển q trình hóa lý theo u cầu cơng nghệ Vì vậy, ta cần tìm hiểu q trình hóa lý thiết bị phản ứng quan hệ phi tuyến, tính xen kênh, động học q trình mơ hình điều khiển, từ hiệu chỉnh lại thiết bị công nghệ Nội dung báo gồm ba phần: Nghiên cứu hóa lý phản ứng thủy phân Acetic Anhydride; Thiết kê mơ hình thiết bị phản ứng; Nghiên cứu điều khiển tuyến tính phản hồi đầu dùng PID đánh giá chất lượng hệ Phản ứng thủy phân Acetic Anhydride Xét phản ứng thủy phân Acetic Anhydride với phương trình phản ứng sau [1]: (1) (CH 3CO)2 O + H 2O ⇔ 2CH 3COOH Cơ chế phản ứng thủy phân mơ tả Hình 2, cặp điện tử tự nguyên tử Oxygen phân tử nước cơng vào ngun tử Carbon nhóm chức ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 11(108).2016, Quyển Andehyde tạo hai phân tử Acetic Acid Bảng Thông số ban đầu thiết bị phản ứng Hình Cấu trúc phản ứng hydro hóa metyl este axit peraxetic Stt Thơng số Ý nghĩa V=0.508.10-3 m3 Thể tích dung dịch phản ứng S=0.00635 m2 CA1=0,5 mol/l =51 kg/m3=100% Đây phản ứng đồng pha phản ứng hydro hóa este hay cịn gọi thủy phân este Cặp electron nguyên tử oxi nước bẻ gãy (cracking) liên kết C-O gốc axit CH3COO- gốc ancol CH3CO+ Phân tử H2O tách thành OH- H+, OH- liên kết với gốc COCH3 H+ liên kết với gốc CH3COO tạo thành hai phân tử CH3COOH Để đảm bảo phản ứng hoàn thành cần cấp nhiệt lượng đủ khuấy trộn đồng dung dịch Về thiết bị, dùng thiết bị theo mẻ hay liên tục [2], [3] Động học thiết bị phản ứng liên tục có khuấy trộn (CSTR – continuous stirred tank reactor) Theo [4], [5] ta có: - Cân thành phần dC A V = F1C A1 − F2 C A − k1 V C A (2) dt - Cân lượng cho phản ứng dT ρVC p = (ρC p FT 1 − ρC p F2T ) −ΔHk1VC A + Q (3) dt (4) Q = K T A(T − TJ ) - Cân lượng cho Jacket dT j ρ jV j C Pj = ρ j C pj F j (T j1 − T j ) − Q dt - Cân khối lượng dh Vmax = F1 − F2 dt V h= % Vmax (5) k1 =0,00791 s-1 Hệ số tốc độ phản ứng nhiệt độ 313K T1=278K Nhiệt độ dung dịch đầu vào T = 313K Nhiệt độ phản ứng E=53408 J/mol R=8,314 J/mol.K Hằng số chất khí lý tưởng 10 ∆H=142000 J/kg Nhiệt phản ứng 11 Cp=4200 J/kg.K Nhiệt dung riêng chất phản ứng 12 ρ=1000 kg/m3 Khối lượng riêng môi chất phản ứng 13 H=0,08 m Mức dung dịch thiết bị phản ứng 14 A= 0,029 m2 15 Cpj =4200 J/kgK 16 Vj=2,5411.10-3 m3 17 ρj =1000kg/m3 Năng lượng phản ứng Diện tích truyền nhiệt Nhiệt dung riêng nước gia nhiệt Thể tích nước gia nhiệt Khối lượng riêng nước gia nhiệt Từ thông số ban đầu thiết bị phản ứng CSTR thủy phân Acetic Anhydride ta thiết kế cấu trúc mơ hình thiết bị phản ứng trình bày Hình 40 (6) 70 (7) tốc độ phản ứng: k1 = k0 e KT hệ số truyền nhiệt (W/m2K), A diện tích truyền nhiệt (m2), T(K) nhiệt độ dung dịch bình phản ứng, ρ (kg/ m3) khối lượng riêng chất phản ứng, Cp (J/kg.K) nhiệt dung riêng chất phản ứngVj (m3) thể tích Jaket, ρj (kg/ m3) khối lượng riêng môi chất gia nhiệt, Cpj (J/kg.K) nhiệt dung riêng môi chất gia nhiệt, Tj1, Tj2 (K) nhiệt độ vào Jacket, Fj (m3/s) lưu lượng môi chất cấp cho Jacket, Q công suất truyền nhiệt từ môi chất đến phản ứng (xem 2.2), h mức dung dich bình phản ứng (tính theo đơn vị %) Tính tốn thơng số lị phản ứng 4.1 Thơng số ban đầu Diện tích đáy thiết bị Nồng độ dung dịch nguyên liệu (đầu vào) CA2=0,047 mol/l = Nồng độ dung dịch sản 4,794 kg/m3=9,3962% phẩm (đầu ra) 190 Trong đó: V (m3) thể tích thiết bị, CA1 CA2 (mol/l, kg/m3 tính theo %) nồng độ đầu vào chất phản ứng, T1(K) nhiệt độ đầu vào chất phản ứng, T(K) nhiệt độ phản ứng trì khơng đổi 313(K), F1 F2 (m3/s) lưu lượng vào chất phản ứng, k1(s) E − RT 43 80 140 40 90 170 Hình Cấu trúc mơ hình thiết bị phản ứng Trên Hình mơ tả mơ hình thiết bị phản ứng liên tục CSTR gồm bình phản ứng có chiều cao 150 mm, bề ngang 90 mm, chứa 0,508.10-3 m3 dung dịch, bao bọc lớp jacket có chiều cao 190 mm, bề ngang 170 mm, chứa 2,5411.10-3 m3 nước gia nhiệt 4.2 Tính tốn thơng số vận hành a Tính lưu lượng hoạt chất cấp vào thiết bị F1 lưu lượng sản phẩm thiết bị F2 44 Mai Thị Đoan Thanh, Nguyễn Đình Lâm, Đồn Quang Vinh Từ phương trình cân thành phần chế độ ổn định ta có: dC A V = F1C A1 − F2 C A − k1 V C A = (8) dt Nếu đảm bảo cân F1 = F2 ta tính được: F1 = F2 = k1VC A C A1 − C A (9) b Tính cơng suất nhiệt cấp cho thiết bị Từ phương trình cân lượng cho phản ứng ta có: dT ρVCp = (ρCp FT 1 − ρCp F2T ) −ΔHk1VCA + Q = (10) dt Công suất nhiệt yêu cầu cho phản ứng tính: Q = (ρC p F1T1 − ρC p F2T ) −ΔHk1VC A (11) c Tính tốn thơng số jacket Phương trình cân lượng Jacket phương trình truyền nhiệt là: dT ρ jV j C pj J = ρ j C pj Fj (T j1 − T j ) − Q = (12 dt Tính lưu lượng nước nóng cấp cho Jacket Fj, chọn Tj1=343K Tj2=323K ta tính lưu lượng Fj Fj = Q ρ j C pj (T j1 − T j ) a Giai đoạn khởi động Phương trình cân thành phần giải sau: dC A dt dC A CA2 = − k1 C A (16) = − k1 dt (17) C A2 = C A0 e − k1t ; k1 = k e − E RT (18) Phương trình cân nhiệt phản ứng ρVC p dT dt = −Vk0 e − E RT ΔH CA2 + kT A(Tj − T ) (19) Cân nhiệt cho jacket ρ jV j C pj dTj dt = ρ j C pj Fj (Tj1 − Tj ) − k T A(Tj − T ) (20) Trên Hình mơ tả đáp ứng nhiệt độ phản ứng T giai đoạn khởi động, nhiệt độ phản ứng T tăng từ nhệt độ ban đầu đến giá trị cân (13) Từ phương trình (4), ta tính KT (W/m2K) sau: Q (14) KT = A(T − T j ) Hình Đáp ứng nhiệt độ phản ứng T giai đoạn khởi động d Tính tốn hệ số tốc độ phản ứng k0, ta có: k1 k0 = e − (15) E RT Trong k1 tốc độ phản ứng (1/s); Thơng số tính tốn trình bày Bảng Trên Hình mơ tả đáp ứng nồng độ đầu CA2 giai đoạn khởi động, nồng độ CA2 giảm từ nồng độ ban đầu giá trị cân Bảng Thông số thiết bị phản ứng CSTR Stt Thông số Ý nghĩa F1= F2= 4,169.10-7 m3/s Lưu lượng dung dịch đầu vào KT=0,002712 (W/m2K) Hệ số truyền nhiệt k0 = 6477839 s-1 Hệ số tốc độ phản ứng Q= 64W Công suất truyền nhiệt Tj1=343 K Nhiệt độ đầu vào nước gia nhiệt Tj2=323 K Nhiệt độ đầu nước gia nhiệt Fj=6,55.10-7 m3/s Lưu lượng đầu vào nước gia nhiệt 4.3 Kiểm tra thiết kế thiết bị phản ứng Để kiểm tra thông số thiết bị thiết kế mục 4.1, ta tiến hành mô kiểm tra vận hành thiết bị (theo [6], [7]) Quá trình mơ có hai giai đoạn: - Giai đoạn khởi động: Nạp hóa chất vào đầy bình, khóa tất van cho gia nhiệt cho nhiệt độ bình tăng nhiệt từ 50C lên 400C kết thúc giai đoạn khởi động - Giai đoạn phản ứng: Khi kết thúc giai đoạn khởi động ta cho thiết bị phản ứng vận hành trạng thái cân Hình Đáp ứng nồng độ sản phẩm đầu CA2 giai đoạn khởi động Trên Hình Hình đánh giá đáp ứng giai đoạn khởi động: Nồng độ sản phẩm đầu CA2 nhiệt độ phản ứng T sau 1500s đạt giá trị cân thiết kế (CA2 đạt 4,794 kg/m3, T đạt 313K) b Giai đoạn phản ứng Mơ hình mơ cân khối lượng (6) theo Hình 6: Hình Mơ hình mơ cân khối lượng ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 11(108).2016, Quyển Trong h mức dung dịch thiết bị phản ứng đại lượng cần điều khiển (tính theo đơn vị tương đối), lưu lượng F2 biến điều khiển (giả thiết điều khiển qua bơm), F1 đại lượng nhiễu Mơ hình mơ cân thành phần (2) theo Hình 7: 45 dựng từ mơ hình mơ từ Hình đến Hình Trên Hình 10 mơ hình mơ thiết bị Xây dựng mơ hình điều khiển 5.1 Mơ hình cấu trúc điều khiển thiết bị phản ứng Mơ hình điều khiển thiết bị CSTR trình bày dạng sơ đồ khối Hình 11 Mơ hình điều khiển lị phản ứng có hai đại lượng cần điều khiển: mức h (cân khối lượng) nồng độ đầu CA2 Tuy nhiên nồng độ CA2 khơng đo trực tuyến được, cần điều khiển qua đại lượng nhiệt độ phản ứng T Hình Mơ hình mơ cân thành phần Mơ hình mơ cân lượng (3) (5) theo Hình Hình 9: Hình 11 Mơ hình cấu trúc điều khiển thiết bị phản ứng Hình Mơ hình mơ q trình cân lượng thiết bị phản ứng Hình Mơ hình mơ q trình cân lượng Jacket thiết bị phản ứng Mạch vòng điều khiển mức (h) ta chọn đại lượng tác động F2 Mạch vòng điều khiển nhiệt độ (nồng độ) chọn đại lượng tác động lưu lượng nước làm nóng Fj Các đại lượng lại xác định nhiễu: F1, T1, CA1, Tj1 Biến tự khơng kiểm sốt Tj2 Từ mơ hình điều khiển tính giá trị điểm cân sau: Khi nồng độ đầu vào CA10 = 51 kg/m3, nhiệt độ đầu vào T10 = 278K; để có nồng độ đầu C A 20 = 4, 794kg / m ứng với nhiệt độ đầu T = 313K ; lưu lượng cấp cho thiết bị xác định F10 = 4,169.10−7 m3 / s ; mức dung dịch thiết bịH = 0.08 m ; công suất nước gia nhiệt Q = 61,2W ; nhiệt độ nước gia nhiệt đầu vào T j10 = 343K ; tính lưu lượng nước gia nhiệt F j = 6,55.10−7 m3 / s 5.2 Chọn thông số làm việc thiết bị sau Để đảm bảo điều khiển cân khối lượng, chọn lưu lượng cực đại đầu thiết bị F1max = 8, 338.10−7 m3 / s (gấp lần giá trị làm việc ổn định) Để đảm bảo cân lượng chọn lưu lượng nước gia nhiệt cực đạiFjmax =13,1.10−7 m3 / s, (gấp lần giá trị làm việc ổn định) Nhận xét: Mơ hình điều khiển lị phản ứng có tính phi tuyến, cụ thể: Phi tuyến vào-ra: hệ số tốc độ phản ứng hàm mũ nhiệt độ: k1 = k e Hình 10 Mơ hình mơ hệ thống thiết bị phản ứng Mơ hình mơ hệ thống thiết bị phản ứng xây − E RT (21) Phi tuyến cấu trúc: Có tích hai biến cần điều khiển nhiệt độ nồng độ (3) Tác động xen kênh: Đại lượng cần điều khiển h lại có quan hệ với đại lượng điều khiển cân thành phần cân lượng 46 Mai Thị Đoan Thanh, Nguyễn Đình Lâm, Đồn Quang Vinh Từ phân tích trên, ta có sơ đồ cấu trúc điều khiển thiết bị phản ứng với hai mạch vòng điều khiển (mức h nhiệt độ phản ứng T) Hình 12: Nhận xét: Khi hệ làm việc ổn định, thời điểm 5000s, ta cho C A1 giảm 10%, T1 tăng 10%, nồng độ đầu C A có giảm lượng nhỏ, nhiên tác động điều khiển PID, đáp ứng nồng độ đầu C A ổn định sai lệch 12,6% b Nhiễu lưu lượng đầu vào F1 ±10% ta có đáp ứng CA2 Hình 15 Hình 12 Cấu trúc điều khiển thiết bị phản ứng với hai mạch vòng điều khiển 5.3 Hệ điều khiển phản hồi PID cho thiết bị phản ứng với mơ hình phi tuyến Trên Hình 13 mơ hình mơ điều khiển thiết bị CSTR, với điều khiển tuyến tính phản hồi PID Hình 15 Đáp ứng nồng độ đầu C A nhiễu F1 ± 10% (CA2 sai lệch 10% giảm sai lệch 8% tăng) Trên Hình 16 đáp ứng mức dung dịch h thay đổi lưu lượng đầu vào F1 ±10%: Hình 16 Đáp ứng mức dung dịch h F1 ± 10% tác động Nhận xét: Mạch vòng điều khiển mức h tác động nhanh khoảng 200s, ổn định Nồng độ CA2 sai lệch 10% Hình 13 Mơ hình mơ điều khiển phản hồi PID thiết bị phản ứng CSTR cấu trúc phi tuyến Chất lượng hệ đánh giá nồng độ sản phẩm phải đảm bảo ≤±5% Khi tất thông số đầu vào không đổi thiết kế nồng độ đầu đảm bảo lượng đặt Tuy nhiên có thay đổi nhiễu đầu vào như: nồng độ, nhiệt độ lưu lượng đầu vào; nhiệt độ đầu vào dung dịch gia nhiệt Sau ta xét ảnh hưởng nhiễu: a Nhiễu C A1 −10%; T1 + 10% tác động, ta có đáp ứng CA2 Hình 14 Hình 14 Đáp ứng nồng độ đầu C A nhiễuC A1 − 10%; T1 + 10%(CA2 đạt 4,192 kg/m3) c Nhiễu đồng thời ta có đáp ứng Hình 17 Nhận xét: Khi hệ làm việc ổn định, thời điểm 5000s, cho nhiễu C A1 giảm 10%, nhiễu T1 tăng 10% Tj1 tăng 10%, F1 giảm 10% nồng độ đầu C A có sai lệch 15% Đến 7500s ta thay đổi T1 giảm 10% F1 tăng 10%, nồng độ đầu C A có sai lệch 30% Hình 17 Đáp ứng nồng độ sản phẩm đầu CA2 nhiễu tác động lúc ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 11(108).2016, Quyển 47 Nhận xét chung: Hệ điều khiển thiết bị phản ứng với điều khiển PID cho mơ hình phi tuyến đảm bảo ổn định nhiễu đầu vào tác động, nhiên tồn sai lệch tĩnh giá trị cho phép Tác động xen kênh mạch vòng mức tới mạch vòng nhiệt độ ngược lại khơng đáng kể, mạch vịng mức tác động nhanh so với mạch vòng nhiệt độ chất lượng sản phẩm đầu đáp ứng yêu cầu Tuy nhiên có nhiễu đầu vào, đặc biệt trường hợp nhiễu tác động đồng thời, chất lượng sản phẩm có sai lệch lớn ngồi giá trị cho phép Điều cho thấy để đảm bảo cho hệ bám lượng đặt bền vững với nhiễu, cần phải nghiên cứu phương pháp thiết kế phản hồi phi tuyến phù hợp Kết luận Bài báo nghiên cứu động học q trình phản ứng thơng số cơng nghệ phản ứng thủy phân hóa Acetic Anhydride, từ thiết kế mơ hình, tính tốn thông số vận hành thiết bị phản ứng CSTR Mơ hình thiết kế kiểm tra qua mơ hai giai đoạn: giai đoạn khởi động giai đoạn vận hành, đánh giá đặc tính thiết bị phản ứng CSTR Mơ hình thử nghiệm lị phản ứng thiết kế theo tính tốn có thơng số đáp ứng yêu cầu đặt ra, kiểm chứng đảm bảo cân phản ứng sản phẩm đầu Với việc phân tích mơ hình cấu trúc điều khiển thiết bị phản ứng cho thấy mơ hình điều khiển thiết bị phản ứng hệ đa biến tác động xen kênh có tính phi tuyến Khi dùng điều khiển tuyến tính phản hồi đầu dùng PID, với thông số đầu vào không đổi, hệ ổn định cho kết [1] N Bin, Các trình, thiết bị cơng nghệ hóa chất thực phẩm vol Tập 5: Nhà xuất Khoa học Kỹ thuật, 2007 TÀI LIỆU THAM KHẢO [2] N A Amenaghawon, E I Osagie, S O Osemwengie, S E Ogbeide, and C O Okieimen, "Modelling and simulation of the batch hydrolysis of acetic anhydride to produce acetic acid”, Nigerian Journal of Technology (NIJOTECH), vol Vol 32 No 3, p 7, 2013 [3] S Haji and C Erkey, "Kinetics of hydrolysis of acetic anhydride by in-situ FTIR spectroscopy - An experiment for the Undergraduate Laboratory”, Chemical Engineering Education, p 6, 2005 [4] B.Q.Khánh, P.Q.Đăng, N H Phương, and V.T.Nguyên, Điều khiển trình, Nhà xuất Khoa học Kỹ thuật, 2014 [5] B Roffel and B Betlem, Process Dynamics and Control Modeling for Control and Prediction: John Wiley & Son, Ltd, 562 [6] J Vojtesek and P Dostal, "Use of MATLAB Environment for Simulation and Control of CSTR”, International journal of mathematics and computers in simulation, vol Volume 5, p 8, 2011 [7] W L Luyben, Process modeling simulation and control for chemical engineers McGraw-Hill 1996 (BBT nhận bài: 23/9/2016, phản biện xong: 12/11/2016) ... Cặp electron nguyên tử oxi nước bẻ gãy (cracking) liên kết C-O gốc axit CH3COO- gốc ancol CH3CO+ Phân tử H2O tách thành OH- H+, OH- liên kết với gốc COCH3 H+ liên kết với gốc CH3COO tạo thành... [4], [5] ta có: - Cân thành phần dC A V = F1C A1 − F2 C A − k1 V C A (2) dt - Cân lượng cho phản ứng dT ρVC p = (ρC p FT 1 − ρC p F2T ) −ΔHk1VC A + Q (3) dt (4) Q = K T A(T − TJ ) - Cân lượng cho...ISSN 185 9-1 531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 11(108).2016, Quyển Andehyde tạo hai phân