1. Dòng đẩy lý tưởng khác cách trộn toàn bồn
Trong mô hình PFR, thành phần thay đổi theo chiều đi của dòng và không có trộn ngược dọc trục. Ở mỗi mặt cắt, mô hình một chiều thường dùng trạng thái đại diện cho toàn tiết diện. Đây là giả thiết dòng chảy, không phải kết luận rằng bất kỳ ống dài nào cũng là PFR lý tưởng.
Tài liệu Kitchin Research Group trình bày phương trình vi phân theo thể tích phản ứng và nghiệm bậc nhất của PFR.[1] Bài dùng phản ứng giả định A → B, không thuận nghịch, đẳng nhiệt, k không đổi, pha lỏng không nén được và không có hạn chế truyền khối. Những điều kiện này phải được giữ khi so sánh với CSTR.
2. Tích phân tốc độ dọc thể tích
Gọi F_A là dòng mol A, Q là lưu lượng thể tích không đổi và v là tọa độ thể tích tính từ cửa vào. Nồng độ C_A = F_A/Q. Cân bằng trên lát mỏng cho:
dF_A/dv = −kC_A = −(k/Q)F_A
F_A(v) = F_A0 exp(−kv/Q)
X = 1 − exp(−kV/Q)
Với Q = 2 m³/h, V = 5 m³, k = 0,8 h⁻¹ và C_A0 = 1,2 kmol/m³, Da = kV/Q = 2. Kết quả X = 0,864665 và C_A ra = 0,162402 kmol/m³. Dòng A vào 2,4 kmol/h, còn lại 0,324805 kmol/h và đã phản ứng 2,075195 kmol/h.
Ở cửa vào, tốc độ tiêu thụ là 0,96 kmol/(m³·h); ở cửa ra chỉ còn 0,129922 kmol/(m³·h). Nhân tốc độ cửa ra với toàn thể tích sẽ không cho lượng phản ứng đúng, vì tốc độ đã giảm liên tục dọc ống. Cân bằng phải tích phân theo v.
3. So sánh với CSTR trên cùng đầu vào
CSTR có cùng Q, V, k cho X = Da/(1 + Da) = 0,666667. PFR trong ví dụ chuyển hóa cao hơn 19,799805 điểm phần trăm. Sự khác biệt xuất phát từ phân bố nồng độ: phần đầu PFR vẫn tiếp xúc với nồng độ A cao, còn CSTR trộn toàn bộ về nồng độ ra.
Phương trình thiết kế theo chuyển hóa trong tài liệu Fogler cho phép viết thể tích dưới dạng tích phân F_A0/(-r_A) đối với PFR.[2] Với hệ bậc nhất và Q không đổi, hai biểu thức đảo là:
V_PFR = −(Q/k) ln(1 − X)
V_CSTR = (Q/k) X/(1 − X)
Tại X = 0,90, PFR cần 5,756463 m³, còn CSTR cần 22,5 m³; tỷ số thể tích là 3,908650. Không được chuyển tỷ số này thành tỷ số giá thành: thiết bị còn khác hình học, truyền nhiệt, trở lực, vật liệu và phương thức bảo trì.
Kết luận về chuyển hóa cũng không nên mở rộng sang mọi phản ứng. Phản ứng nhiều bước, phản ứng thuận nghịch, hiệu ứng nhiệt và mục tiêu chọn lọc có thể làm bài toán chọn thiết bị khác hẳn ví dụ một phản ứng này.
4. Dãy bồn nhỏ tiến dần tới dòng đẩy trong mô hình
Chia tổng thể tích 5 m³ thành N CSTR bằng nhau nối tiếp, không có dòng phụ, cùng nhiệt độ và k. Mỗi bồn có Da/N và tỷ số nồng độ ra/vào bằng 1/(1 + Da/N). Qua N bồn:
C_A,N/C_A0 = (1 + Da/N)^(−N)
X_N = 1 − (1 + Da/N)^(−N)
Với Da = 2, một bồn cho 0,666667; hai bồn cho 0,75; năm bồn cho 0,814066; mười bồn cho 0,838494. Khi N tăng không giới hạn, tỷ số nồng độ tiến về exp(−Da), trùng nghiệm PFR.
Đây là giới hạn toán học của mô hình bồn lý tưởng nối tiếp, không phải lời khẳng định phải lắp thật vô số bồn. Nếu có lưu lượng cấp thêm, nhiệt độ khác nhau hoặc k thay đổi theo bậc, không thể dùng nguyên công thức lũy thừa trên.
Với hai bồn có tổng thể tích cố định, tích mẫu số là (1 + kV_1/Q)(1 + kV_2/Q). Khi V_1 + V_2 không đổi, tích V_1V_2 lớn nhất ở hai thể tích bằng nhau. Trong riêng mô hình này, chia đều làm nồng độ ra nhỏ nhất. Lập luận không áp tự động cho động học phi tuyến khác.
5. Thứ tự một CSTR và một PFR có thể không làm đổi kết quả
Xét CSTR 3 m³ và PFR 2 m³ nối tiếp, vẫn Q = 2 m³/h và k = 0,8 h⁻¹. Nếu CSTR đi trước, tỷ số còn lại của A qua cả hai thiết bị là:
C_A,out/C_A0 = [1/(1 + 0,8 × 3/2)] exp(−0,8 × 2/2)
Nếu đổi thứ tự, hai thừa số chỉ đổi vị trí, nên kết quả không đổi. Chuyển hóa bằng khoảng 0,795760. Nguyên nhân là hai phép biến đổi đều tuyến tính theo nồng độ vào và các hệ số giữ cố định.
Nếu thay bằng phản ứng bậc hai, nhiều phản ứng, hoặc để nhiệt độ thiết bị sau phụ thuộc nhiệt tỏa ở thiết bị trước, tính giao hoán này không còn được bảo đảm. Bài toán thứ tự cần xét phương trình thật, không suy từ tên loại thiết bị.
6. Khí giãn nở làm mất giả thiết Q không đổi
Xét một kịch bản độc lập: khí lý tưởng A tinh khiết phản ứng A → 2B, đẳng nhiệt, áp suất không đổi và không có khí trơ. Khi chuyển hóa X, tổng mol tăng theo 1 + X, nên Q = Q_0(1 + X), còn C_A = C_A0(1 − X)/(1 + X).
Với động học tiêu thụ A bậc nhất, thay vào cân bằng PFR rồi tích phân:
V = (Q_0/k) integral[(1 + X)/(1 − X)dX]
V = (Q_0/k)[−2 ln(1 − X) − X]
Chọn Q_0 = 2 m³/h, k = 0,8 h⁻¹ và X = 0,90, thể tích cần là 9,262925 m³. Công thức Q không đổi cho 5,756463 m³ sẽ đánh giá thấp trong bộ giả thiết này. Dòng đã giãn nở làm nồng độ và tốc độ khác với trường hợp lỏng ban đầu.
Nếu áp suất giảm dọc lớp xúc tác, nhiệt độ đổi hoặc khí không lý tưởng, cần bổ sung các quan hệ tương ứng. Tốc độ tính trên kg xúc tác cũng không được nhân trực tiếp với m³ thiết bị; phải đổi cơ sở hoặc dùng phương trình theo khối lượng xúc tác.
7. Cùng thời gian trung bình vẫn có thể khác chuyển hóa
Giả định dòng được chia đều thành hai nhánh PFR riêng, một nhánh có thời gian tiếp xúc 1 h, nhánh kia 4 h. Trung bình theo lượng dòng là 2,5 h, bằng thời gian của PFR cơ sở. Với k = 0,8 h⁻¹, phần A còn lại sau khi trộn hai nhánh là 0,5 exp(−0,8) + 0,5 exp(−3,2).
Chuyển hóa tổng chỉ 0,754954, thấp hơn 0,864665 của mọi phần tử cùng tiếp xúc 2,5 h. Nhánh nhanh để lại nhiều A; phần lưu lâu hơn ở nhánh chậm không bù tương ứng do quan hệ hàm mũ. Đây là ví dụ tự xây dựng cho thấy chỉ kiểm V/Q chưa đủ để mô tả phân phối dòng.
8. Từ nghiệm đẳng nhiệt tới mô hình thiết bị
ReactorLab của Richard K. Herz mô tả một PFR có lớp hạt xúc tác và áo trao đổi nhiệt, kèm các giả thiết về khối lượng riêng, nhiệt dung và trao đổi nhiệt.[3] Việc ghi các giả thiết ấy cho thấy mô hình cần thêm cấu trúc khi nhiệt độ và tích lũy nhiệt tham gia.
Trước khi dùng thể tích tính được, cần giữ hồ sơ động học, cơ sở tốc độ, lưu lượng tại điều kiện thực và mô hình nhiệt. Một phép so sánh có ích là chỉ ra điều kiện nào tạo ra ưu thế về chuyển hóa, đồng thời giữ các tiêu chí truyền nhiệt, chọn lọc và khả năng vận hành cho bước đánh giá thiết bị tiếp theo.
Nguồn tham khảo
[1] Single Unit: The PFR — Differentiable Flowsheets — Kitchin Research Group, Carnegie Mellon University. Truy cập ngày 09/11/2026.
[2] Elements of Chemical Reaction Engineering — Chapter2: Design Equations — H. Scott Fogler, University of Michigan. Truy cập ngày 09/11/2026.
[3] Plug Flow Reactor with heat exchange jacket — Web Lab7 information — Richard K. Herz, ReactorLab. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Trong mô hình PFR, thành phần thay đổi theo chiều đi của dòng và không có trộn ngược dọc trục. Ở mỗi mặt cắt, mô hình một chiều thường dùng trạng thái đại diện cho toàn tiết diện. Đây là giả thiết dòng chảy, không phải kết luận rằng bất kỳ ống dài nào cũng là PFR lý tưởng.
Tài liệu Kitchin Research Group trình bày phương trình vi phân theo thể tích phản ứng và nghiệm bậc nhất của PFR.[1] Bài dùng phản ứng giả định A → B, không thuận nghịch, đẳng nhiệt, k không đổi, pha lỏng không nén được và không có hạn chế truyền khối. Những điều kiện này phải được giữ khi so sánh với CSTR.
2. Tích phân tốc độ dọc thể tích
Gọi F_A là dòng mol A, Q là lưu lượng thể tích không đổi và v là tọa độ thể tích tính từ cửa vào. Nồng độ C_A = F_A/Q. Cân bằng trên lát mỏng cho:
dF_A/dv = −kC_A = −(k/Q)F_A
F_A(v) = F_A0 exp(−kv/Q)
X = 1 − exp(−kV/Q)
Với Q = 2 m³/h, V = 5 m³, k = 0,8 h⁻¹ và C_A0 = 1,2 kmol/m³, Da = kV/Q = 2. Kết quả X = 0,864665 và C_A ra = 0,162402 kmol/m³. Dòng A vào 2,4 kmol/h, còn lại 0,324805 kmol/h và đã phản ứng 2,075195 kmol/h.
Ở cửa vào, tốc độ tiêu thụ là 0,96 kmol/(m³·h); ở cửa ra chỉ còn 0,129922 kmol/(m³·h). Nhân tốc độ cửa ra với toàn thể tích sẽ không cho lượng phản ứng đúng, vì tốc độ đã giảm liên tục dọc ống. Cân bằng phải tích phân theo v.
3. So sánh với CSTR trên cùng đầu vào
CSTR có cùng Q, V, k cho X = Da/(1 + Da) = 0,666667. PFR trong ví dụ chuyển hóa cao hơn 19,799805 điểm phần trăm. Sự khác biệt xuất phát từ phân bố nồng độ: phần đầu PFR vẫn tiếp xúc với nồng độ A cao, còn CSTR trộn toàn bộ về nồng độ ra.
Phương trình thiết kế theo chuyển hóa trong tài liệu Fogler cho phép viết thể tích dưới dạng tích phân F_A0/(-r_A) đối với PFR.[2] Với hệ bậc nhất và Q không đổi, hai biểu thức đảo là:
V_PFR = −(Q/k) ln(1 − X)
V_CSTR = (Q/k) X/(1 − X)
Tại X = 0,90, PFR cần 5,756463 m³, còn CSTR cần 22,5 m³; tỷ số thể tích là 3,908650. Không được chuyển tỷ số này thành tỷ số giá thành: thiết bị còn khác hình học, truyền nhiệt, trở lực, vật liệu và phương thức bảo trì.
Kết luận về chuyển hóa cũng không nên mở rộng sang mọi phản ứng. Phản ứng nhiều bước, phản ứng thuận nghịch, hiệu ứng nhiệt và mục tiêu chọn lọc có thể làm bài toán chọn thiết bị khác hẳn ví dụ một phản ứng này.
4. Dãy bồn nhỏ tiến dần tới dòng đẩy trong mô hình
Chia tổng thể tích 5 m³ thành N CSTR bằng nhau nối tiếp, không có dòng phụ, cùng nhiệt độ và k. Mỗi bồn có Da/N và tỷ số nồng độ ra/vào bằng 1/(1 + Da/N). Qua N bồn:
C_A,N/C_A0 = (1 + Da/N)^(−N)
X_N = 1 − (1 + Da/N)^(−N)
Với Da = 2, một bồn cho 0,666667; hai bồn cho 0,75; năm bồn cho 0,814066; mười bồn cho 0,838494. Khi N tăng không giới hạn, tỷ số nồng độ tiến về exp(−Da), trùng nghiệm PFR.
Đây là giới hạn toán học của mô hình bồn lý tưởng nối tiếp, không phải lời khẳng định phải lắp thật vô số bồn. Nếu có lưu lượng cấp thêm, nhiệt độ khác nhau hoặc k thay đổi theo bậc, không thể dùng nguyên công thức lũy thừa trên.
Với hai bồn có tổng thể tích cố định, tích mẫu số là (1 + kV_1/Q)(1 + kV_2/Q). Khi V_1 + V_2 không đổi, tích V_1V_2 lớn nhất ở hai thể tích bằng nhau. Trong riêng mô hình này, chia đều làm nồng độ ra nhỏ nhất. Lập luận không áp tự động cho động học phi tuyến khác.
5. Thứ tự một CSTR và một PFR có thể không làm đổi kết quả
Xét CSTR 3 m³ và PFR 2 m³ nối tiếp, vẫn Q = 2 m³/h và k = 0,8 h⁻¹. Nếu CSTR đi trước, tỷ số còn lại của A qua cả hai thiết bị là:
C_A,out/C_A0 = [1/(1 + 0,8 × 3/2)] exp(−0,8 × 2/2)
Nếu đổi thứ tự, hai thừa số chỉ đổi vị trí, nên kết quả không đổi. Chuyển hóa bằng khoảng 0,795760. Nguyên nhân là hai phép biến đổi đều tuyến tính theo nồng độ vào và các hệ số giữ cố định.
Nếu thay bằng phản ứng bậc hai, nhiều phản ứng, hoặc để nhiệt độ thiết bị sau phụ thuộc nhiệt tỏa ở thiết bị trước, tính giao hoán này không còn được bảo đảm. Bài toán thứ tự cần xét phương trình thật, không suy từ tên loại thiết bị.
6. Khí giãn nở làm mất giả thiết Q không đổi
Xét một kịch bản độc lập: khí lý tưởng A tinh khiết phản ứng A → 2B, đẳng nhiệt, áp suất không đổi và không có khí trơ. Khi chuyển hóa X, tổng mol tăng theo 1 + X, nên Q = Q_0(1 + X), còn C_A = C_A0(1 − X)/(1 + X).
Với động học tiêu thụ A bậc nhất, thay vào cân bằng PFR rồi tích phân:
V = (Q_0/k) integral[(1 + X)/(1 − X)dX]
V = (Q_0/k)[−2 ln(1 − X) − X]
Chọn Q_0 = 2 m³/h, k = 0,8 h⁻¹ và X = 0,90, thể tích cần là 9,262925 m³. Công thức Q không đổi cho 5,756463 m³ sẽ đánh giá thấp trong bộ giả thiết này. Dòng đã giãn nở làm nồng độ và tốc độ khác với trường hợp lỏng ban đầu.
Nếu áp suất giảm dọc lớp xúc tác, nhiệt độ đổi hoặc khí không lý tưởng, cần bổ sung các quan hệ tương ứng. Tốc độ tính trên kg xúc tác cũng không được nhân trực tiếp với m³ thiết bị; phải đổi cơ sở hoặc dùng phương trình theo khối lượng xúc tác.
7. Cùng thời gian trung bình vẫn có thể khác chuyển hóa
Giả định dòng được chia đều thành hai nhánh PFR riêng, một nhánh có thời gian tiếp xúc 1 h, nhánh kia 4 h. Trung bình theo lượng dòng là 2,5 h, bằng thời gian của PFR cơ sở. Với k = 0,8 h⁻¹, phần A còn lại sau khi trộn hai nhánh là 0,5 exp(−0,8) + 0,5 exp(−3,2).
Chuyển hóa tổng chỉ 0,754954, thấp hơn 0,864665 của mọi phần tử cùng tiếp xúc 2,5 h. Nhánh nhanh để lại nhiều A; phần lưu lâu hơn ở nhánh chậm không bù tương ứng do quan hệ hàm mũ. Đây là ví dụ tự xây dựng cho thấy chỉ kiểm V/Q chưa đủ để mô tả phân phối dòng.
8. Từ nghiệm đẳng nhiệt tới mô hình thiết bị
ReactorLab của Richard K. Herz mô tả một PFR có lớp hạt xúc tác và áo trao đổi nhiệt, kèm các giả thiết về khối lượng riêng, nhiệt dung và trao đổi nhiệt.[3] Việc ghi các giả thiết ấy cho thấy mô hình cần thêm cấu trúc khi nhiệt độ và tích lũy nhiệt tham gia.
Trước khi dùng thể tích tính được, cần giữ hồ sơ động học, cơ sở tốc độ, lưu lượng tại điều kiện thực và mô hình nhiệt. Một phép so sánh có ích là chỉ ra điều kiện nào tạo ra ưu thế về chuyển hóa, đồng thời giữ các tiêu chí truyền nhiệt, chọn lọc và khả năng vận hành cho bước đánh giá thiết bị tiếp theo.
Nguồn tham khảo
[1] Single Unit: The PFR — Differentiable Flowsheets — Kitchin Research Group, Carnegie Mellon University. Truy cập ngày 09/11/2026.
[2] Elements of Chemical Reaction Engineering — Chapter2: Design Equations — H. Scott Fogler, University of Michigan. Truy cập ngày 09/11/2026.
[3] Plug Flow Reactor with heat exchange jacket — Web Lab7 information — Richard K. Herz, ReactorLab. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
