Flash cần hai loại thông tin: cân bằng vật chất và quan hệ cân bằng pha
Khi một dòng nguyên liệu phân thành pha lỏng và pha hơi, không thể xác định lượng mỗi pha chỉ bằng việc biết tổng lưu lượng. Thành phần nguyên liệu và xu hướng phân bố từng cấu tử giữa hai pha đều quan trọng. Phương trình Rachford–Rice giúp tìm tỷ phần hơi khi các tỷ số cân bằng K đã biết; nó không tự tạo ra K đúng cho mọi nhiệt độ và áp suất.
Bài này chọn hỗn hợp hai cấu tử giả định A và B, với K được cho trước và giữ cố định. Đây là bài tập đại số và kiểm tra cân bằng. Các cấu tử không đại diện cho methane, propane hoặc một phân đoạn dầu cụ thể; không suy ra nhiệt độ, áp suất hay tính chất thực từ các số K đã đặt. Cách tách rõ phần cân bằng vật chất này giúp kiểm tra đầu ra của phần mềm mô phỏng mà không nhầm nó với toàn bộ mô hình nhiệt động.
Cơ sở mol và các ký hiệu
Gọi F là lưu lượng mol nguyên liệu; L là lưu lượng mol lỏng; V là lưu lượng mol hơi. Trong ví dụ chúng đều tính bằng kmol/h. z_i, x_i và y_i lần lượt là phần mol cấu tử i trong nguyên liệu, lỏng và hơi. Tất cả là đại lượng không thứ nguyên, có tổng bằng một trong từng dòng. Tỷ số K_i = y_i/x_i và tỷ phần hơi β = V/F cũng không có đơn vị.
F = L + V; Fz_i = Lx_i + Vy_i
Thay L = F(1−β), V = Fβ và y_i = K_ix_i vào cân bằng từng cấu tử, rút ra:
x_i = z_i/[1 + β(K_i − 1)]; y_i = K_ix_i
Từ điều kiện tổng thành phần hai pha, phương trình cần giải là:
f(β) = Σ z_i(K_i − 1)/[1 + β(K_i − 1)] = 0
Tài liệu Chemicals trình bày công thức trên theo cơ sở mol và mô tả các bộ giải Rachford–Rice. [1] Tài liệu MOOSE cũng sử dụng dạng phương trình flash, nhưng định nghĩa các biến theo khối lượng trong ngữ cảnh mô hình của họ. [2] Hai quy ước không được trộn trong một bảng tính. Nếu z là phần khối lượng mà K lại là tỷ số phần mol, kết quả không có cơ sở dù bộ giải vẫn trả về một con số.
Bộ số giả định có thể kiểm tra bằng tay
Chọn F = 100 kmol/h, z_A = 0,40 và z_B = 0,60. Đặt K_A = 2,50, K_B = 0,50. Khi ấy phương trình trở thành:
0,60/(1 + 1,50β) − 0,30/(1 − 0,50β) = 0
Nhân chéo trong miền các mẫu số dương: 0,60 − 0,30β = 0,30 + 0,45β. Do đó β = 0,40. Lượng hơi V = 40 kmol/h; lượng lỏng L = 60 kmol/h. Đây là tỷ phần mol: không được kết luận 40% thể tích bình là hơi hoặc 40% khối lượng dòng là hơi nếu chưa có khối lượng mol và mật độ pha.
Tính thành phần pha lỏng: x_A = 0,40/(1+0,40×1,50) = 0,25; x_B = 0,60/(1−0,40×0,50) = 0,75. Thành phần pha hơi là y_A = 2,50×0,25 = 0,625 và y_B = 0,50×0,75 = 0,375. Cả tổng x và tổng y đều bằng một.
Đóng cân bằng từng cấu tử
Nguyên liệu chứa 40 kmol/h A. Dòng lỏng mang 60×0,25 = 15 kmol/h A; dòng hơi mang 40×0,625 = 25 kmol/h A. Tổng đầu ra đúng 40 kmol/h. Với B, dòng lỏng mang 45 kmol/h và dòng hơi mang 15 kmol/h, tổng đúng 60 kmol/h.
Ba phép kiểm tra khác nhau đã được thỏa mãn: tổng lưu lượng mol, tổng thành phần mỗi pha và cân bằng từng cấu tử. Chỉ kiểm tra F = L + V chưa đủ. Một bộ số sai thành phần vẫn có thể cộng tổng lưu lượng đúng; một chương trình chuẩn hóa x và y riêng rẽ cũng có thể làm tổng thành phần đẹp nhưng phá cân bằng cấu tử nếu β chưa đúng.
Tài liệu của nhóm Kitchin tại Carnegie Mellon minh họa bước tính lại x và y sau khi giải phân số hơi, cùng cách giữ nghiệm trong các giới hạn thích hợp. [3] Trong một bảng tính kiểm tra, nên lưu phần dư cân bằng thay vì chỉ hiện trạng thái “đã hội tụ”. Mức dung sai cần phù hợp độ chính xác dữ liệu và mục đích tính, không tự chọn một ngưỡng để biến kết quả chưa rõ thành đạt.
Kiểm tra khả năng tồn tại nghiệm trong khoảng 0 đến 1
Với bộ số này, f(0) = 0,30 và f(1) = −0,36. Hàm đổi dấu giữa hai đầu. Đạo hàm của mô hình K cố định là f′(β) = −Σz_i(K_i−1)²/[1+β(K_i−1)]². Trong miền mẫu số hợp lệ, đạo hàm không dương và ở ví dụ là âm; nghiệm vì vậy có thể được khoanh bằng phương pháp chia đôi.
Một vòng lặp số có thể bắt đầu tại β = 0,5, tính dấu phần dư rồi thu hẹp khoảng. Không cần dùng một thuật toán phức tạp để kiểm chứng ví dụ hai cấu tử này. Sự đơn giản có lợi: kết quả số có thể đối chiếu trực tiếp với nghiệm đại số β = 0,4. Nếu hai cách khác nhau, cần xem lại đầu vào, dấu của K_i−1 hoặc quy ước cơ sở.
Nếu cả hai giá trị đầu mút không khoanh một nghiệm hai pha, không nên ép β vào giữa 0 và 1 rồi tiếp tục tính như thể hai pha chắc chắn tồn tại. Trong mô hình nhiệt động đầy đủ, số pha và tính ổn định cần được kiểm tra. Một nghiệm của phương trình vật chất với K giả định không tự chứng minh trạng thái ấy là cân bằng vật lý của một hỗn hợp thực.
Độ nhạy với K và thành phần đầu vào
Giữ z_A = 0,40, z_B = 0,60 và K_B = 0,50, tăng K_A từ 2,50 lên 2,80. Với hai cấu tử, có thể rút β = −[z_A(K_A−1)+z_B(K_B−1)]/[(K_A−1)(K_B−1)]. Thay số cho β xấp xỉ 0,4667; lượng hơi tăng lên khoảng 46,667 kmol/h. Sự thay đổi không chỉ ảnh hưởng lượng hơi mà còn thay đổi thành phần của cả hai pha.
Ở một kịch bản khác, giữ K như ban đầu nhưng tăng z_A lên 0,45, đồng thời giảm z_B xuống 0,55. Nghiệm β trở thành khoảng 0,5333. Không được chỉ tăng z_A mà giữ z_B = 0,60 vì tổng phần mol sẽ thành 1,05. Độ nhạy thành phần luôn phải xét cùng điều kiện tổng bằng một.
Các kịch bản này không có nghĩa nhiệt độ tăng một mức nào đó chắc chắn làm K_A đổi từ 2,50 lên 2,80. Muốn nối thay đổi K với nhiệt độ hoặc áp suất cần một mô hình cân bằng pha và dữ liệu tính chất cụ thể. Ở đây chúng chỉ lượng hóa tác động của sự bất định hoặc thay đổi đầu vào toán học.
Giới hạn khi chuyển sang thiết bị thực
Một bình tách thực còn có động học tách giọt, thời gian lưu, entrainment, tổn thất áp suất và trao đổi nhiệt. Cân bằng pha lý tưởng không bảo đảm thiết bị đạt đúng thành phần cân bằng tại mọi chế độ. Nếu bài toán là flash đoạn nhiệt, cần thêm cân bằng enthalpy để tìm nhiệt độ; giữ K cố định rồi suy ra nhiệt độ ra là thiếu phương trình.
Hồ sơ có thể kiểm tra nên ghi nguồn mô hình K, cơ sở mol hay khối lượng, điều kiện nhiệt độ–áp suất, phương pháp xác định số pha, phần dư cân bằng và độ nhạy. Ví dụ này cung cấp một trường hợp chuẩn nhỏ để thử cách triển khai công thức, đồng thời nhắc rằng “bộ giải hội tụ” và “mô hình đại diện đúng hệ thực” là hai bước đánh giá cần thực hiện riêng.
Nguồn tham khảo
[1] Rachford-Rice Equation Solvers (chemicals.rachford_rice) — Chemicals documentation — Caleb Bell and contributors. Truy cập ngày 09/11/2026.
[2] Compositional flash calculations — Idaho National Laboratory — MOOSE PorousFlow. Truy cập ngày 09/11/2026.
[3] Solvers and Utilities — Differentiable Flowsheets — Carnegie Mellon University — Kitchin Research Group / Differentiable Flowsheets contributors. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Khi một dòng nguyên liệu phân thành pha lỏng và pha hơi, không thể xác định lượng mỗi pha chỉ bằng việc biết tổng lưu lượng. Thành phần nguyên liệu và xu hướng phân bố từng cấu tử giữa hai pha đều quan trọng. Phương trình Rachford–Rice giúp tìm tỷ phần hơi khi các tỷ số cân bằng K đã biết; nó không tự tạo ra K đúng cho mọi nhiệt độ và áp suất.
Bài này chọn hỗn hợp hai cấu tử giả định A và B, với K được cho trước và giữ cố định. Đây là bài tập đại số và kiểm tra cân bằng. Các cấu tử không đại diện cho methane, propane hoặc một phân đoạn dầu cụ thể; không suy ra nhiệt độ, áp suất hay tính chất thực từ các số K đã đặt. Cách tách rõ phần cân bằng vật chất này giúp kiểm tra đầu ra của phần mềm mô phỏng mà không nhầm nó với toàn bộ mô hình nhiệt động.
Cơ sở mol và các ký hiệu
Gọi F là lưu lượng mol nguyên liệu; L là lưu lượng mol lỏng; V là lưu lượng mol hơi. Trong ví dụ chúng đều tính bằng kmol/h. z_i, x_i và y_i lần lượt là phần mol cấu tử i trong nguyên liệu, lỏng và hơi. Tất cả là đại lượng không thứ nguyên, có tổng bằng một trong từng dòng. Tỷ số K_i = y_i/x_i và tỷ phần hơi β = V/F cũng không có đơn vị.
F = L + V; Fz_i = Lx_i + Vy_i
Thay L = F(1−β), V = Fβ và y_i = K_ix_i vào cân bằng từng cấu tử, rút ra:
x_i = z_i/[1 + β(K_i − 1)]; y_i = K_ix_i
Từ điều kiện tổng thành phần hai pha, phương trình cần giải là:
f(β) = Σ z_i(K_i − 1)/[1 + β(K_i − 1)] = 0
Tài liệu Chemicals trình bày công thức trên theo cơ sở mol và mô tả các bộ giải Rachford–Rice. [1] Tài liệu MOOSE cũng sử dụng dạng phương trình flash, nhưng định nghĩa các biến theo khối lượng trong ngữ cảnh mô hình của họ. [2] Hai quy ước không được trộn trong một bảng tính. Nếu z là phần khối lượng mà K lại là tỷ số phần mol, kết quả không có cơ sở dù bộ giải vẫn trả về một con số.
Bộ số giả định có thể kiểm tra bằng tay
Chọn F = 100 kmol/h, z_A = 0,40 và z_B = 0,60. Đặt K_A = 2,50, K_B = 0,50. Khi ấy phương trình trở thành:
0,60/(1 + 1,50β) − 0,30/(1 − 0,50β) = 0
Nhân chéo trong miền các mẫu số dương: 0,60 − 0,30β = 0,30 + 0,45β. Do đó β = 0,40. Lượng hơi V = 40 kmol/h; lượng lỏng L = 60 kmol/h. Đây là tỷ phần mol: không được kết luận 40% thể tích bình là hơi hoặc 40% khối lượng dòng là hơi nếu chưa có khối lượng mol và mật độ pha.
Tính thành phần pha lỏng: x_A = 0,40/(1+0,40×1,50) = 0,25; x_B = 0,60/(1−0,40×0,50) = 0,75. Thành phần pha hơi là y_A = 2,50×0,25 = 0,625 và y_B = 0,50×0,75 = 0,375. Cả tổng x và tổng y đều bằng một.
Đóng cân bằng từng cấu tử
Nguyên liệu chứa 40 kmol/h A. Dòng lỏng mang 60×0,25 = 15 kmol/h A; dòng hơi mang 40×0,625 = 25 kmol/h A. Tổng đầu ra đúng 40 kmol/h. Với B, dòng lỏng mang 45 kmol/h và dòng hơi mang 15 kmol/h, tổng đúng 60 kmol/h.
Ba phép kiểm tra khác nhau đã được thỏa mãn: tổng lưu lượng mol, tổng thành phần mỗi pha và cân bằng từng cấu tử. Chỉ kiểm tra F = L + V chưa đủ. Một bộ số sai thành phần vẫn có thể cộng tổng lưu lượng đúng; một chương trình chuẩn hóa x và y riêng rẽ cũng có thể làm tổng thành phần đẹp nhưng phá cân bằng cấu tử nếu β chưa đúng.
Tài liệu của nhóm Kitchin tại Carnegie Mellon minh họa bước tính lại x và y sau khi giải phân số hơi, cùng cách giữ nghiệm trong các giới hạn thích hợp. [3] Trong một bảng tính kiểm tra, nên lưu phần dư cân bằng thay vì chỉ hiện trạng thái “đã hội tụ”. Mức dung sai cần phù hợp độ chính xác dữ liệu và mục đích tính, không tự chọn một ngưỡng để biến kết quả chưa rõ thành đạt.
Kiểm tra khả năng tồn tại nghiệm trong khoảng 0 đến 1
Với bộ số này, f(0) = 0,30 và f(1) = −0,36. Hàm đổi dấu giữa hai đầu. Đạo hàm của mô hình K cố định là f′(β) = −Σz_i(K_i−1)²/[1+β(K_i−1)]². Trong miền mẫu số hợp lệ, đạo hàm không dương và ở ví dụ là âm; nghiệm vì vậy có thể được khoanh bằng phương pháp chia đôi.
Một vòng lặp số có thể bắt đầu tại β = 0,5, tính dấu phần dư rồi thu hẹp khoảng. Không cần dùng một thuật toán phức tạp để kiểm chứng ví dụ hai cấu tử này. Sự đơn giản có lợi: kết quả số có thể đối chiếu trực tiếp với nghiệm đại số β = 0,4. Nếu hai cách khác nhau, cần xem lại đầu vào, dấu của K_i−1 hoặc quy ước cơ sở.
Nếu cả hai giá trị đầu mút không khoanh một nghiệm hai pha, không nên ép β vào giữa 0 và 1 rồi tiếp tục tính như thể hai pha chắc chắn tồn tại. Trong mô hình nhiệt động đầy đủ, số pha và tính ổn định cần được kiểm tra. Một nghiệm của phương trình vật chất với K giả định không tự chứng minh trạng thái ấy là cân bằng vật lý của một hỗn hợp thực.
Độ nhạy với K và thành phần đầu vào
Giữ z_A = 0,40, z_B = 0,60 và K_B = 0,50, tăng K_A từ 2,50 lên 2,80. Với hai cấu tử, có thể rút β = −[z_A(K_A−1)+z_B(K_B−1)]/[(K_A−1)(K_B−1)]. Thay số cho β xấp xỉ 0,4667; lượng hơi tăng lên khoảng 46,667 kmol/h. Sự thay đổi không chỉ ảnh hưởng lượng hơi mà còn thay đổi thành phần của cả hai pha.
Ở một kịch bản khác, giữ K như ban đầu nhưng tăng z_A lên 0,45, đồng thời giảm z_B xuống 0,55. Nghiệm β trở thành khoảng 0,5333. Không được chỉ tăng z_A mà giữ z_B = 0,60 vì tổng phần mol sẽ thành 1,05. Độ nhạy thành phần luôn phải xét cùng điều kiện tổng bằng một.
Các kịch bản này không có nghĩa nhiệt độ tăng một mức nào đó chắc chắn làm K_A đổi từ 2,50 lên 2,80. Muốn nối thay đổi K với nhiệt độ hoặc áp suất cần một mô hình cân bằng pha và dữ liệu tính chất cụ thể. Ở đây chúng chỉ lượng hóa tác động của sự bất định hoặc thay đổi đầu vào toán học.
Giới hạn khi chuyển sang thiết bị thực
Một bình tách thực còn có động học tách giọt, thời gian lưu, entrainment, tổn thất áp suất và trao đổi nhiệt. Cân bằng pha lý tưởng không bảo đảm thiết bị đạt đúng thành phần cân bằng tại mọi chế độ. Nếu bài toán là flash đoạn nhiệt, cần thêm cân bằng enthalpy để tìm nhiệt độ; giữ K cố định rồi suy ra nhiệt độ ra là thiếu phương trình.
Hồ sơ có thể kiểm tra nên ghi nguồn mô hình K, cơ sở mol hay khối lượng, điều kiện nhiệt độ–áp suất, phương pháp xác định số pha, phần dư cân bằng và độ nhạy. Ví dụ này cung cấp một trường hợp chuẩn nhỏ để thử cách triển khai công thức, đồng thời nhắc rằng “bộ giải hội tụ” và “mô hình đại diện đúng hệ thực” là hai bước đánh giá cần thực hiện riêng.
Nguồn tham khảo
[1] Rachford-Rice Equation Solvers (chemicals.rachford_rice) — Chemicals documentation — Caleb Bell and contributors. Truy cập ngày 09/11/2026.
[2] Compositional flash calculations — Idaho National Laboratory — MOOSE PorousFlow. Truy cập ngày 09/11/2026.
[3] Solvers and Utilities — Differentiable Flowsheets — Carnegie Mellon University — Kitchin Research Group / Differentiable Flowsheets contributors. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
