1. Mô hình khuấy trộn lý tưởng bắt đầu từ nồng độ trong bồn
CSTR là mô hình thiết bị phản ứng khuấy trộn liên tục, trong đó thành phần dòng ra bằng thành phần pha đã trộn trong thiết bị. Vì vậy tốc độ phản ứng phải được tính tại nồng độ trong bồn, không phải nồng độ dòng cấp. Tài liệu của Kitchin Research Group trình bày cân bằng mol và giả thiết trộn đều theo cách này.[1]
Bài xét phản ứng giả định A → B, không thuận nghịch, bậc nhất theo A, pha lỏng có khối lượng riêng gần không đổi. Nhiệt độ giữ cố định nên hằng số tốc độ k không đổi; thể tích pha phản ứng V và lưu lượng thể tích Q vào ra cũng không đổi. Các số liệu được tự đặt, chưa đại diện cho phản ứng hydrocarbon hay xúc tác cụ thể.
2. Từ cân bằng mol tới chuyển hóa
Gọi C_A0 là nồng độ cấp, C_A là nồng độ trong bồn và dòng ra. Tốc độ tiêu thụ dương là −r_A = kC_A. Ở trạng thái ổn định, không còn tích lũy:
0 = Q C_A0 − Q C_A − VkC_A
C_A = C_A0/(1 + kV/Q)
X = 1 − C_A/C_A0 = Da/(1 + Da)
Da = k tau; tau = V/Q
V tính bằng m³, Q bằng m³/h, k bằng h⁻¹ thì Da không có đơn vị. Nếu k được cho theo phút, phải đổi thời gian trước khi nhân. Tên “hằng số” không cho phép bỏ đơn vị của k.
Chọn Q = 2 m³/h, C_A0 = 1,2 kmol/m³, V = 5 m³ và k = 0,8 h⁻¹. Ta có tau = 2,5 h, Da = 2, C_A = 0,4 kmol/m³ và X = 0,666667. Dòng A cấp 2,4 kmol/h; dòng A ra 0,8 kmol/h; lượng tiêu thụ 1,6 kmol/h.
Kiểm độc lập bằng tốc độ phản ứng: VkC_A = 5 × 0,8 × 0,4 = 1,6 kmol/h. Nếu dùng nhầm C_A0 trong tốc độ, sẽ tính 4,8 kmol/h, lớn hơn cả A được cấp; đây là dấu hiệu vi phạm cân bằng rất rõ.
3. Thể tích tăng nhanh khi chuyển hóa tiến gần một
Giải ngược phương trình ổn định:
V = (Q/k) X/(1 − X)
Với Q và k như trên, X = 0,90 cần V = 22,5 m³; X = 0,95 cần 47,5 m³. Tăng 5 điểm phần trăm ở mức chuyển hóa cao làm thể tích tăng 111,111111%, không phải 5%.
Khi X tiến đến 1, mẫu số tiến về 0 và thể tích trong mô hình tăng không giới hạn. Điều này không có nghĩa mọi công nghệ đều cần bồn vô hạn để đạt một chỉ tiêu sản phẩm; nó là giới hạn của mô hình CSTR đơn, phản ứng bậc nhất đã chọn.
Đạo hàm dV/dX = Q/[k(1 − X)²] tại X = 0,90 bằng 250 m³ trên một đơn vị X. Biến thiên nhỏ 0,01 cho ước lượng tuyến tính thêm 2,5 m³, trong khi tính đúng từ 0,90 tới 0,91 cho thêm 2,777778 m³. Độ nhạy tăng nhanh khiến xấp xỉ tuyến tính kém dần khi khoảng thay đổi rộng.
4. Nồng độ cấp và tỷ lệ chuyển hóa là hai câu hỏi khác nhau
Trong mô hình này, X chỉ phụ thuộc kV/Q. Nếu tăng C_A0 gấp đôi mà giữ Q, V, k và nhiệt độ không đổi, X vẫn bằng 0,666667, nhưng C_A ra tăng thành 0,8 kmol/m³ và A tiêu thụ thành 3,2 kmol/h.
Vì thế “chuyển hóa không đổi” không đồng nghĩa tải phản ứng hoặc tải nhiệt không đổi. Nếu phản ứng tỏa nhiệt, việc giữ k cố định còn cần một cân bằng nhiệt phù hợp. Mô hình CSTR của MathWorks ghép phương trình nồng độ với nhiệt độ để mô tả tác động của phản ứng và trao đổi nhiệt.[3] Bài hiện tại chưa giải bài toán nhiệt đó.
Khi dữ liệu thực cho thấy chuyển hóa phụ thuộc mạnh nồng độ cấp, cần xem các giả thiết bậc nhất, nhiệt độ không đổi và thể tích thực tham gia còn phù hợp hay không. Không nên chỉnh k tùy ý để che mọi sai khác của mô hình.
5. Suy k từ phép đo cần trạng thái ổn định
Nếu Q, V, C_A0 và C_A đã được đo ở cùng trạng thái ổn định, có thể suy:
k = (Q/V)(C_A0/C_A − 1)
Giả sử C_A đo được 0,42 kmol/m³ thay vì 0,40, còn các dữ liệu khác giữ nguyên. k suy ra là 0,742857 h⁻¹. Đây là giá trị biểu kiến theo mô hình, chưa tự chứng minh động học đã thay đổi.
Tại C_A = 0,40, đạo hàm theo nồng độ ra là −QC_A0/(VC_A²) = −3 theo đơn vị tương ứng. Một độ lệch nhỏ 0,01 kmol/m³ có thể làm k suy ra lệch khoảng 0,03 h⁻¹ nếu các biến khác giữ cố định. Cần lưu phương pháp lấy mẫu, độ ổn định và thể tích pha thực khi diễn giải con số.
6. Trong khởi động, phần thiếu ở cửa ra còn có thể đang tích lũy
Tài liệu tham khảo của Fogler xử lý khởi động CSTR bằng nồng độ và cân bằng có tích lũy.[2] Với các giả thiết thể tích và lưu lượng không đổi của bài:
dC_A/dt = (Q/V)(C_A0 − C_A) − kC_A
C_A(t) = C_A,ss + [C_A(0) − C_A,ss] exp[−(Q/V + k)t]
Giả định bồn ban đầu đã chứa đủ 5 m³ dung môi nhưng chưa có A, nên C_A(0) = 0. Đây không phải bồn rỗng đang được làm đầy. Với bộ số cơ sở, C_A,ss = 0,4 và Q/V + k = 1,2 h⁻¹. Sau một giờ, C_A = 0,279522 kmol/m³.
Tại thời điểm đó, A vào trừ A ra bằng 1,840955 kmol/h. Nhưng A thực bị phản ứng chỉ 1,118089 kmol/h; phần 0,722866 kmol/h còn lại đang tích lũy. Nếu lấy toàn bộ chênh vào–ra làm lượng đã phản ứng, kết quả sẽ sai ngay cả khi mọi phép đo đều chính xác.
Thời gian để sai khác nồng độ so với trạng thái cuối còn 5% của sai khác ban đầu là −ln(0,05)/1,2 = 2,496444 h. Với mức 1%, thời gian là 3,837642 h. Các mốc này được suy từ phương trình cụ thể, không phải quy tắc luôn chờ một số cố định lần V/Q cho mọi phản ứng.
7. Chia bồn song song cần chia lưu lượng phù hợp
Nếu thay bồn 5 m³ bằng hai bồn 2,5 m³ song song, mỗi bồn nhận 1 m³/h cùng nồng độ cấp, từng bồn vẫn có tau = 2,5 h và X = 0,666667. Tổng thể tích bằng nhau chỉ cho kết quả này khi phân phối dòng tương ứng.
Nếu vẫn hai bồn 2,5 m³ nhưng dòng chia thành 1,5 và 0,5 m³/h, chuyển hóa từng nhánh thành 0,571429 và 0,80. Chuyển hóa của dòng gộp phải lấy trọng số theo lượng A cấp từng nhánh; vì nồng độ cấp bằng nhau, kết quả là (1,5 × 0,571429 + 0,5 × 0,80)/2 = 0,628571.
Trung bình không trọng số sẽ cho 0,685714 và còn cao hơn trường hợp chia đều, dẫn tới kết luận sai. Đối với hai nhánh khác nồng độ cấp, trọng số phải là dòng mol A, không chỉ lưu lượng thể tích.
8. Dùng kết quả trong đúng phạm vi
Bảng tính cần nêu rõ phản ứng, cơ sở tốc độ, nhiệt độ, thể tích pha, dòng cấp và điều kiện trộn. Tự kiểm ba cách: khép mol A, kiểm đơn vị của Da và kiểm giới hạn V tiến về 0 hoặc tăng lớn.
Một nghiệm CSTR bậc nhất là nền tảng để học cách liên hệ động học với thể tích. Thiết bị thực còn có thể chịu giới hạn truyền nhiệt, truyền khối, phân bố dòng và suy giảm xúc tác; các yếu tố ấy phải được đưa vào mô hình tương ứng trước khi dùng kết quả để chọn hay thay đổi thiết bị.
Nguồn tham khảo
[1] Single Unit: The CSTR — Differentiable Flowsheets — Kitchin Research Group, Carnegie Mellon University. Truy cập ngày 09/11/2026.
[2] Professional Reference Shelf — Startup of a CSTR — H. Scott Fogler, University of Michigan. Truy cập ngày 09/11/2026.
[3] CSTR Model — Nonlinear Model — MathWorks official documentation. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
CSTR là mô hình thiết bị phản ứng khuấy trộn liên tục, trong đó thành phần dòng ra bằng thành phần pha đã trộn trong thiết bị. Vì vậy tốc độ phản ứng phải được tính tại nồng độ trong bồn, không phải nồng độ dòng cấp. Tài liệu của Kitchin Research Group trình bày cân bằng mol và giả thiết trộn đều theo cách này.[1]
Bài xét phản ứng giả định A → B, không thuận nghịch, bậc nhất theo A, pha lỏng có khối lượng riêng gần không đổi. Nhiệt độ giữ cố định nên hằng số tốc độ k không đổi; thể tích pha phản ứng V và lưu lượng thể tích Q vào ra cũng không đổi. Các số liệu được tự đặt, chưa đại diện cho phản ứng hydrocarbon hay xúc tác cụ thể.
2. Từ cân bằng mol tới chuyển hóa
Gọi C_A0 là nồng độ cấp, C_A là nồng độ trong bồn và dòng ra. Tốc độ tiêu thụ dương là −r_A = kC_A. Ở trạng thái ổn định, không còn tích lũy:
0 = Q C_A0 − Q C_A − VkC_A
C_A = C_A0/(1 + kV/Q)
X = 1 − C_A/C_A0 = Da/(1 + Da)
Da = k tau; tau = V/Q
V tính bằng m³, Q bằng m³/h, k bằng h⁻¹ thì Da không có đơn vị. Nếu k được cho theo phút, phải đổi thời gian trước khi nhân. Tên “hằng số” không cho phép bỏ đơn vị của k.
Chọn Q = 2 m³/h, C_A0 = 1,2 kmol/m³, V = 5 m³ và k = 0,8 h⁻¹. Ta có tau = 2,5 h, Da = 2, C_A = 0,4 kmol/m³ và X = 0,666667. Dòng A cấp 2,4 kmol/h; dòng A ra 0,8 kmol/h; lượng tiêu thụ 1,6 kmol/h.
Kiểm độc lập bằng tốc độ phản ứng: VkC_A = 5 × 0,8 × 0,4 = 1,6 kmol/h. Nếu dùng nhầm C_A0 trong tốc độ, sẽ tính 4,8 kmol/h, lớn hơn cả A được cấp; đây là dấu hiệu vi phạm cân bằng rất rõ.
3. Thể tích tăng nhanh khi chuyển hóa tiến gần một
Giải ngược phương trình ổn định:
V = (Q/k) X/(1 − X)
Với Q và k như trên, X = 0,90 cần V = 22,5 m³; X = 0,95 cần 47,5 m³. Tăng 5 điểm phần trăm ở mức chuyển hóa cao làm thể tích tăng 111,111111%, không phải 5%.
Khi X tiến đến 1, mẫu số tiến về 0 và thể tích trong mô hình tăng không giới hạn. Điều này không có nghĩa mọi công nghệ đều cần bồn vô hạn để đạt một chỉ tiêu sản phẩm; nó là giới hạn của mô hình CSTR đơn, phản ứng bậc nhất đã chọn.
Đạo hàm dV/dX = Q/[k(1 − X)²] tại X = 0,90 bằng 250 m³ trên một đơn vị X. Biến thiên nhỏ 0,01 cho ước lượng tuyến tính thêm 2,5 m³, trong khi tính đúng từ 0,90 tới 0,91 cho thêm 2,777778 m³. Độ nhạy tăng nhanh khiến xấp xỉ tuyến tính kém dần khi khoảng thay đổi rộng.
4. Nồng độ cấp và tỷ lệ chuyển hóa là hai câu hỏi khác nhau
Trong mô hình này, X chỉ phụ thuộc kV/Q. Nếu tăng C_A0 gấp đôi mà giữ Q, V, k và nhiệt độ không đổi, X vẫn bằng 0,666667, nhưng C_A ra tăng thành 0,8 kmol/m³ và A tiêu thụ thành 3,2 kmol/h.
Vì thế “chuyển hóa không đổi” không đồng nghĩa tải phản ứng hoặc tải nhiệt không đổi. Nếu phản ứng tỏa nhiệt, việc giữ k cố định còn cần một cân bằng nhiệt phù hợp. Mô hình CSTR của MathWorks ghép phương trình nồng độ với nhiệt độ để mô tả tác động của phản ứng và trao đổi nhiệt.[3] Bài hiện tại chưa giải bài toán nhiệt đó.
Khi dữ liệu thực cho thấy chuyển hóa phụ thuộc mạnh nồng độ cấp, cần xem các giả thiết bậc nhất, nhiệt độ không đổi và thể tích thực tham gia còn phù hợp hay không. Không nên chỉnh k tùy ý để che mọi sai khác của mô hình.
5. Suy k từ phép đo cần trạng thái ổn định
Nếu Q, V, C_A0 và C_A đã được đo ở cùng trạng thái ổn định, có thể suy:
k = (Q/V)(C_A0/C_A − 1)
Giả sử C_A đo được 0,42 kmol/m³ thay vì 0,40, còn các dữ liệu khác giữ nguyên. k suy ra là 0,742857 h⁻¹. Đây là giá trị biểu kiến theo mô hình, chưa tự chứng minh động học đã thay đổi.
Tại C_A = 0,40, đạo hàm theo nồng độ ra là −QC_A0/(VC_A²) = −3 theo đơn vị tương ứng. Một độ lệch nhỏ 0,01 kmol/m³ có thể làm k suy ra lệch khoảng 0,03 h⁻¹ nếu các biến khác giữ cố định. Cần lưu phương pháp lấy mẫu, độ ổn định và thể tích pha thực khi diễn giải con số.
6. Trong khởi động, phần thiếu ở cửa ra còn có thể đang tích lũy
Tài liệu tham khảo của Fogler xử lý khởi động CSTR bằng nồng độ và cân bằng có tích lũy.[2] Với các giả thiết thể tích và lưu lượng không đổi của bài:
dC_A/dt = (Q/V)(C_A0 − C_A) − kC_A
C_A(t) = C_A,ss + [C_A(0) − C_A,ss] exp[−(Q/V + k)t]
Giả định bồn ban đầu đã chứa đủ 5 m³ dung môi nhưng chưa có A, nên C_A(0) = 0. Đây không phải bồn rỗng đang được làm đầy. Với bộ số cơ sở, C_A,ss = 0,4 và Q/V + k = 1,2 h⁻¹. Sau một giờ, C_A = 0,279522 kmol/m³.
Tại thời điểm đó, A vào trừ A ra bằng 1,840955 kmol/h. Nhưng A thực bị phản ứng chỉ 1,118089 kmol/h; phần 0,722866 kmol/h còn lại đang tích lũy. Nếu lấy toàn bộ chênh vào–ra làm lượng đã phản ứng, kết quả sẽ sai ngay cả khi mọi phép đo đều chính xác.
Thời gian để sai khác nồng độ so với trạng thái cuối còn 5% của sai khác ban đầu là −ln(0,05)/1,2 = 2,496444 h. Với mức 1%, thời gian là 3,837642 h. Các mốc này được suy từ phương trình cụ thể, không phải quy tắc luôn chờ một số cố định lần V/Q cho mọi phản ứng.
7. Chia bồn song song cần chia lưu lượng phù hợp
Nếu thay bồn 5 m³ bằng hai bồn 2,5 m³ song song, mỗi bồn nhận 1 m³/h cùng nồng độ cấp, từng bồn vẫn có tau = 2,5 h và X = 0,666667. Tổng thể tích bằng nhau chỉ cho kết quả này khi phân phối dòng tương ứng.
Nếu vẫn hai bồn 2,5 m³ nhưng dòng chia thành 1,5 và 0,5 m³/h, chuyển hóa từng nhánh thành 0,571429 và 0,80. Chuyển hóa của dòng gộp phải lấy trọng số theo lượng A cấp từng nhánh; vì nồng độ cấp bằng nhau, kết quả là (1,5 × 0,571429 + 0,5 × 0,80)/2 = 0,628571.
Trung bình không trọng số sẽ cho 0,685714 và còn cao hơn trường hợp chia đều, dẫn tới kết luận sai. Đối với hai nhánh khác nồng độ cấp, trọng số phải là dòng mol A, không chỉ lưu lượng thể tích.
8. Dùng kết quả trong đúng phạm vi
Bảng tính cần nêu rõ phản ứng, cơ sở tốc độ, nhiệt độ, thể tích pha, dòng cấp và điều kiện trộn. Tự kiểm ba cách: khép mol A, kiểm đơn vị của Da và kiểm giới hạn V tiến về 0 hoặc tăng lớn.
Một nghiệm CSTR bậc nhất là nền tảng để học cách liên hệ động học với thể tích. Thiết bị thực còn có thể chịu giới hạn truyền nhiệt, truyền khối, phân bố dòng và suy giảm xúc tác; các yếu tố ấy phải được đưa vào mô hình tương ứng trước khi dùng kết quả để chọn hay thay đổi thiết bị.
Nguồn tham khảo
[1] Single Unit: The CSTR — Differentiable Flowsheets — Kitchin Research Group, Carnegie Mellon University. Truy cập ngày 09/11/2026.
[2] Professional Reference Shelf — Startup of a CSTR — H. Scott Fogler, University of Michigan. Truy cập ngày 09/11/2026.
[3] CSTR Model — Nonlinear Model — MathWorks official documentation. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
