WesterntechVN – Khi ∂C/∂t = 0 — Bài Toán PFR Có Nghiệm Đóng Đẹp Nhất

Phương trình PT 7-18 (dòng đẩy không khuếch tán) là phương trình vi phân riêng (PDE) với hai biến độc lập x và t. Trong điều kiện vận hành thực tế, kênh dài và bể lọc sinh học XLNT hầu như luôn ở trạng thái ổn định (steady state) — nồng độ tại mỗi vị trí x cố định theo thời gian, chỉ thay đổi theo vị trí dọc chiều dài.

Khi áp đặt điều kiện ∂C/∂t = 0 vào PT 7-18 và giả thiết phản ứng bậc nhất (f(C) = C), bài toán PDE phức tạp giản lược thành một phương trình vi phân thường (ODE) đơn giản — PT 7-19 — có nghiệm tường minh tuyệt đẹp: C = C₀·e^(-k·θ) (PT 7-21).

GS. Lê Văn Cát chỉ ra sự tương đồng sâu sắc: nghiệm này giống hệt với nghiệm của phản ứng gián đoạn (batch) PT 7-3 — chỉ khác ở chỗ biến thời gian t trong batch được thay bởi thời gian lưu thủy lực θ = V/Q trong PFR. Đây là một trong những kết quả đẹp nhất của lý thuyết kỹ thuật phản ứng: bình phản ứng liên tục dòng đẩy ở trạng thái ổn định tương đương với bình phản ứng gián đoạn trong thời gian θ.

Phần 1: Dẫn Xuất PT 7-19 — Từ PDE Đến ODE

1.1 Điều Kiện Trạng Thái Ổn Định ∂C/∂t = 0

Trong thực tế vận hành PFR, trạng thái ổn định được xác lập khi:

• Lưu lượng Q không đổi
• Nồng độ đầu vào C₀ không đổi
• Điều kiện môi trường (nhiệt độ, pH, DO) ổn định
• Quần xã vi sinh vật đã thích nghi

Khi đó, nồng độ tại mỗi điểm x cố định: ∂C/∂t = 0, mặc dù C vẫn thay đổi theo x.

1.2 Từ PT 7-18 Đến PT 7-19

Thay ∂C/∂t = 0 và f(C) = C (bậc nhất) vào PT 7-18:

−0 = v·∂C/∂x + k·C

v·dC/dx + k·C = 0 (PT 7-19)

(Viết ∂C/∂x thành dC/dx vì C chỉ còn là hàm một biến x)

PT 7-19 là ODE tuyến tính bậc nhất đơn giản, dạng chuẩn:

dC/dx = −(k/v)·C

Đây là dạng hàm mũ cơ bản: tốc độ thay đổi nồng độ theo không gian tỷ lệ với nồng độ hiện tại tại vị trí đó.

1.3 Ý Nghĩa Vật Lý PT 7-19

PT 7-19 nói lên điều gì về vật lý của dòng đẩy ở trạng thái ổn định?

Mỗi điểm x trong bể có vai trò độc lập: Tại mỗi vị trí x, tốc độ tiêu thụ cơ chất = k·C(x). Cơ chất bị tiêu thụ khi đi qua từng đoạn nhỏ dx của bình, với lượng tiêu thụ tỷ lệ với nồng độ còn lại.

Không có tương tác ngược chiều: Nồng độ tại x chỉ phụ thuộc vào nồng độ tại các vị trí x’ < x (phía thượng lưu), không phụ thuộc vào x’ > x (phía hạ lưu). Đây là đặc tính căn bản phân biệt PFR với CSTR (nơi mọi điểm trong bể tương tác với nhau).

Phần 2: Tích Phân PT 7-19 — PT 7-20 Và Nghiệm PT 7-21

2.1 Tách Biến Và Tích Phân — PT 7-20

Từ PT 7-19: dC/C = −(k/v)·dx

Tách biến và viết dưới dạng tích phân:

∫dC/C = −(k/v)·∫dx (PT 7-20)

Lấy tích phân định thức từ đầu vào (x = 0, C = C₀) đến vị trí bất kỳ (x = x, C = C(x)):

∫[từ C₀ đến C] dC/C = −(k/v)·∫[từ 0 đến x] dx

ln(C) − ln(C₀) = −(k/v)·x

ln(C/C₀) = −k·x/v

C(x) = C₀·e^(−k·x/v)

2.2 Nghiệm Tại Đầu Ra — PT 7-21

Tại đầu ra x = L:

C(L) = C₀·e^(−k·L/v)

Sử dụng v = Q/A (A là diện tích mặt cắt ngang) và V = L·A (thể tích bình):

L/v = L/(Q/A) = L·A/Q = V/Q = θ

C = C₀·e^(−k·L/v) = C₀·e^(−k·θ) (PT 7-21)

Trong đó θ = V/Q = L·A/Q là thời gian lưu thủy lực (HRT) của bình phản ứng.

2.3 Sự Tương Đồng Với Phản Ứng Gián Đoạn PT 7-3

So sánh hai phương trình:

PT 7-3 (Batch): C(t) = C₀·e^(−k·t) PT 7-21 (PFR steady state): C = C₀·e^(−k·θ)

Hai phương trình đồng dạng hoàn toàn — chỉ khác nhau ở biến:

• Batch: biến là thời gian phản ứng t (giờ)
• PFR: biến là thời gian lưu thủy lực θ = V/Q (giờ)

Ý nghĩa vật lý sâu sắc: Trong PFR ở trạng thái ổn định, một “phần tử chất lỏng” di chuyển từ đầu vào đến đầu ra mất đúng θ = V/Q giây. Trong θ giây đó, nó trải qua phản ứng bậc nhất y hệt như trong bể gián đoạn với thời gian phản ứng t = θ. Trạng thái ổn định không gian trong PFR tương đương với quá trình thời gian trong batch.

GS. Lê Văn Cát kết luận chính xác: “So với lời giải trong trường hợp tĩnh (batch), sự phân bố nồng độ từ đầu đến cuối bình phản ứng tương tự như mối quan hệ của nó theo thời gian.”

Phần 3: Hiệu Suất PFR Ở Trạng Thái Ổn Định

3.1 Biểu Thức Hiệu Suất E Của PFR

Từ PT 7-21 và định nghĩa hiệu suất E = 1 − C/C₀:

E_PFR = 1 − e^(−k·θ)

So sánh với hiệu suất CSTR trạng thái ổn định (PT 7-14):

E_CSTR = 1 − 1/(1 + k·θ)

(Với cùng k và θ = τ)

3.2 So Sánh Định Lượng PFR vs CSTR — Ai Vượt Trội?

Bảng so sánh E_PFR và E_CSTR theo Da = k·θ:

PFR luôn vượt trội CSTR với cùng k và θ. Chênh lệch lớn nhất ở Da ≈ 2-4 (vùng thiết kế phổ biến).

3.3 So Sánh Thể Tích Cần Thiết Cho Cùng E

Ngược lại, để đạt cùng một hiệu suất E mục tiêu:

• PFR: θ_PFR = −ln(1−E)/k
• CSTR: θ_CSTR = E/[k·(1−E)]

Tỷ số: θ_CSTR/θ_PFR = E/[(1−E)·(−ln(1−E))]

Ví dụ với E = 90%: θ_CSTR/θ_PFR = 0,9/[0,1 × ln(10)] = 0,9/(0,1 × 2,303) = 3,91 lần

Ý nghĩa kinh tế: CSTR cần thể tích V = Q·θ_CSTR = 3,91 × Q·θ_PFR lớn hơn 3,91 lần so với PFR để đạt E = 90%. Chọn PFR tiết kiệm 74% thể tích bể so với CSTR đơn!

Phần 4: Bảng So Sánh Ba Hệ Thống — Batch, PFR, CSTR

4.1 Bảng Tổng Hợp Toàn Diện

Nhận xét quan trọng: Batch và PFR cho cùng hiệu suất với cùng k và thời gian (t hoặc θ) — đây là hệ quả trực tiếp của sự tương đồng toán học giữa PT 7-3 và PT 7-21. Cả hai đều ưu việt hơn CSTR về hiệu suất với cùng thể tích.

4.2 Tam Giác Hiệu Suất — Gradient Nồng Độ Và Tốc Độ Phản Ứng

Lý do căn bản tại sao PFR và Batch tốt hơn CSTR:

Trong PFR/Batch: Tốc độ phản ứng thay đổi từ v_max = k·C₀ (đầu vào/ban đầu) → v_min = k·Ce (đầu ra/cuối). Toàn bộ quá trình diễn ra ở nồng độ cao nhất có thể tại mỗi thời điểm/vị trí.

Trong CSTR: Toàn bộ bể ở nồng độ Ce (thấp nhất) → tốc độ phản ứng = k·Ce (thấp nhất). Phần lớn thể tích bể “lãng phí” do làm việc ở nồng độ thấp.

Nguyên lý tổng quát: Để cực đại hiệu quả phản ứng, nên duy trì nồng độ cao nhất có thể trong suốt quá trình → PFR (gradient không gian) hay Batch (gradient thời gian) đều thực hiện điều này. CSTR không thể vì tính chất khuấy trộn đều.

Phần 5: Thiết Kế PFR Từ PT 7-21 — Ứng Dụng Thực Tế Tại Việt Nam

5.1 Tính L Hoặc V Từ Yêu Cầu Ce

Từ PT 7-21: C = C₀·e^(-k·θ)

Đảo ngược: θ = −ln(C/C₀)/k = ln(C₀/C)/k

θ_min = ln(C₀/Ce) / k

V_min = Q · θ_min = Q · ln(C₀/Ce) / k

Hay nếu biết v và diện tích A:

L_min = v · θ_min = (Q/A) · ln(C₀/Ce) / k

5.2 Ví Dụ Thiết Kế Kênh Oxy Hóa Thẳng

Bài toán: Thiết kế kênh oxy hóa thẳng (straight oxidation channel) xử lý nước thải sinh hoạt khu đô thị mới tại Huế.

Thông số đầu vào:

• Q = 3.000 m³/ngày = 125 m³/giờ
• C₀_BOD = 200 mg/L
• Yêu cầu: Ce_BOD ≤ 20 mg/L (QCVN 14:2008 Cột A)
• Nhiệt độ trung bình T = 24°C

Tính E và θ:

• E = 1 − 20/200 = 0,90 (90%)
• k₂₄ = 0,18 × (1,04)^(24-20) = 0,18 × 1,170 = 0,211 giờ⁻¹
• θ_min = ln(200/20)/0,211 = ln(10)/0,211 = 2,303/0,211 = 10,9 giờ
• Với SF = 1,2: θ_thiết kế = 1,2 × 10,9 = 13,1 giờ ≈ 13 giờ

Tính thể tích và kích thước:

• V = Q·θ = 125 × 13 = 1.625 m³
• Chọn chiều rộng b = 4 m, chiều sâu h = 3 m → A = 12 m²
• Chiều dài: L = V/A = 1.625/12 = 135 m
• Vận tốc dòng: v = Q/A = 125/(12 × 3600) = 0,00289 m/s = 0,29 cm/s

Kiểm tra điều kiện dòng đẩy (Pe_L): Với v nhỏ (0,29 cm/s) và L lớn (135 m), cần đảm bảo không có phân tán đáng kể. Tốc độ thực tế cần khuấy trộn ngang (transverse mixing) để tránh phân tầng, nhưng không khuấy dọc chiều kênh.

Kiểm tra mùa đông (T = 18°C): k₁₈ = 0,18 × (1,04)^(18-20) = 0,18 × 0,925 = 0,166 giờ⁻¹ θ_đông = ln(10)/0,166 = 2,303/0,166 = 13,9 giờ < θ_thiết kế = 13 giờ ✗

→ Không đạt mùa đông! Cần thiết kế lại với θ = 14 giờ hoặc chọn T_thiết kế = 18°C: θ_thiết kế_an toàn = 1,15 × 13,9 = 16 giờ → V = 125 × 16 = 2.000 m³, L = 2.000/12 ≈ 167 m

5.3 So Sánh Chi Phí: PFR Kênh Dài vs CSTR Aerotank

Với cùng Q = 3.000 m³/ngày, C₀ = 200 mg/L, Ce = 20 mg/L, k = 0,211 giờ⁻¹ tại T = 24°C:

PFR (kênh dài): V_PFR = Q·ln(C₀/Ce)/k = 125 × 10,9 = 1.363 m³ (chưa SF)

CSTR đơn: V_CSTR = Q·(C₀/Ce − 1)/k = 125 × (10 − 1)/0,211 = 125 × 42,7 = 5.337 m³ (chưa SF)

Tỷ lệ: V_CSTR/V_PFR = 5.337/1.363 = 3,91 lần

CSTR cần thể tích gấp 3,91 lần PFR — chi phí xây dựng bê tông cốt thép lớn hơn ~3,91 lần. Với giá xây dựng bể bê tông tại Huế khoảng 3 triệu đồng/m³, chênh lệch chi phí:

Δ Chi phí = (5.337 − 1.363) × 3.000.000 = 3.974 × 3.000.000 ≈ 12 tỷ đồng

Đây là lý do kỹ thuật và kinh tế thuyết phục cho việc ưu tiên thiết kế kênh dài (PFR) thay vì bể vuông khuấy trộn (CSTR) trong các dự án XLNT đô thị quy mô vừa tại Việt Nam.

5.4 Áp Dụng PT 7-21 Cho Nitrat Hóa Trong Kênh Dài

PT 7-21 cũng áp dụng cho quá trình nitrat hóa trong kênh hiếu khí dài:

NH₄-N đầu vào: C₀_N = 45 mg/L Yêu cầu đầu ra: Ce_N = 5 mg/L (QCVN 14 Cột A) k_N,25°C = 0,15 giờ⁻¹ (Nitrosomonas, ở nồng độ thấp)

θ_min_N = ln(45/5)/k_N = ln(9)/0,15 = 2,197/0,15 = 14,6 giờ

So sánh với θ_BOD = 10,9 giờ: nitrat hóa yêu cầu HRT dài hơn → θ_thiết kế phải đủ cho cả hai quá trình.

Điều này cho thấy: trong hệ thống kênh dài xử lý tổng hợp BOD và N, chiều dài kênh phải thiết kế theo yêu cầu khó hơn (nitrat hóa với θ = 14,6 giờ), không phải theo BOD.

Phần 6: Mở Rộng PT 7-21 — Xấp Xỉ Cho Phản Ứng Phi Tuyến

6.1 PFR Với Phản Ứng Michaelis-Menten

Khi động học tuân theo Michaelis-Menten (không phải đơn giản bậc nhất):

v = Vmax·C/(Ks + C)

Phương trình PFR trạng thái ổn định:

v·dC/dx = −Vmax·C/(Ks + C)

Tách biến: (Ks + C)·dC/C = −(Vmax/v)·dx

Tích phân: Ks·ln(C₀/C) + (C₀ − C) = (Vmax/v)·L = Vmax·θ

Đây là phương trình ẩn (implicit) — không có nghiệm tường minh dạng C = f(θ). Phải giải số.

Xấp xỉ bậc nhất hiệu dụng (khi C << Ks):

v ≈ (Vmax/Ks)·C = k_eff·C → C ≈ C₀·e^(-k_eff·θ) (PT 7-21 áp dụng gần đúng)

Với C₀ << Ks: k_eff = Vmax/Ks — đây là điều kiện vận hành bình thường của bể xử lý nước thải sinh hoạt (nồng độ BOD << Ks_BOD ≈ 100-200 mg/L).

6.2 Điều Chỉnh PT 7-21 Theo Nhiệt Độ

Kết hợp PT 7-21 với hiệu chỉnh nhiệt độ:

k_T = k₂₀ × (1,04)^(T-20) (bùn hoạt tính hiếu khí)

C_T = C₀·e^(-k_T·θ) = C₀·e^(-k₂₀·(1,04)^(T-20)·θ)

Để duy trì cùng Ce = C_T khi T thay đổi theo mùa, cần điều chỉnh θ:

θ_T = ln(C₀/Ce)/k_T = θ₂₀/(1,04)^(T-20)

Tại Hà Nội mùa đông (T = 15°C): θ₁₅ = θ₂₀/(1,04)^(-5) = θ₂₀/0,822 = 1,22·θ₂₀

→ Mùa đông cần HRT dài hơn 22% so với mùa hè/thu ở 20°C. Thiết kế phải theo θ₁₅ để đảm bảo đạt tiêu chuẩn quanh năm.

Kết Luận

Chuỗi phương trình PT 7-19 → PT 7-20 → PT 7-21 thể hiện một trong những kết quả thanh lịch nhất của lý thuyết kỹ thuật phản ứng: bình phản ứng dòng đẩy lý tưởng ở trạng thái ổn định — với giả thiết D = 0, ∂C/∂t = 0, và phản ứng bậc nhất — cho nghiệm C = C₀·e^(-k·θ), hoàn toàn tương đồng với bình phản ứng gián đoạn trong thời gian θ = V/Q.

Sự tương đồng này không chỉ có giá trị lý thuyết. Trong thực tiễn thiết kế hệ thống XLNT tại Việt Nam, nó cho phép kỹ sư dùng cùng một công thức cho cả hai loại hệ thống: đối với hệ batch (SBR), thiết kế thời gian phản ứng t; đối với hệ liên tục dòng đẩy (kênh dài, bể lọc sinh học), thiết kế thời gian lưu thủy lực θ = V/Q. Và quan trọng hơn cả: cả PFR và batch đều ưu việt hơn CSTR khoảng 3-4 lần về thể tích cần thiết để đạt E = 90% — lợi thế kinh tế hàng chục tỷ đồng trong các dự án XLNT quy mô vừa và lớn.