Hai Cơ Chế Vận Chuyển Chất Trong Dòng Phản Ứng

Khi nước thải chảy qua một bình phản ứng hình trụ dài — kiểu bình phản ứng dòng lý tưởng (Plug Flow Reactor — PFR) — nồng độ chất ô nhiễm không đồng đều tại mọi điểm mà thay đổi liên tục theo vị trí dọc theo chiều dài bình. Sự phân bố nồng độ theo không gian này là kết quả của hai cơ chế vận chuyển chất hoạt động đồng thời:

Cơ chế thứ nhất — đối lưu (convection): Dòng chảy với vận tốc tuyến tính v mang chất ô nhiễm di chuyển theo chiều dài bình. Chất ô nhiễm được "đẩy" từ đầu vào về phía đầu ra bởi dòng chảy.

Cơ chế thứ hai — khuếch tán phân tử (molecular diffusion): Sự chênh lệch nồng độ dọc theo bình tạo ra gradient nồng độ, dẫn đến quá trình khuếch tán phân tử từ vùng nồng độ cao về vùng nồng độ thấp — ngược chiều dòng chảy. Cơ chế này được mô tả bởi định luật khuếch tán Fick.

GS. Lê Văn Cát trong "Xử Lý Nước Thải Giàu Hợp Chất Nitơ và Phốtpho" (NXB Khoa học Tự nhiên và Công nghệ, 2007) đã xây dựng phương trình mô tả cả hai cơ chế này qua PT 7-15 (định luật Fick thứ hai) và PT 7-16 (phương trình tổng hợp đối lưu-khuếch tán). Bài viết này phân tích chi tiết cơ sở vật lý và toán học của hai phương trình này — nền tảng để hiểu tại sao PFR hoạt động khác về bản chất so với CSTR (bể khuấy trộn đều).
 

Phần 1: Định Luật Fick Thứ Nhất — Khởi Điểm Lý Thuyết Khuếch Tán

1.1 Thông Lượng Khuếch Tán Và Gradient Nồng Độ

Định luật Fick thứ nhất (Adolf Fick, 1855) mô tả thông lượng khuếch tán phân tử J (lượng chất vận chuyển qua đơn vị diện tích trong đơn vị thời gian):

J = −D · ∂C/∂x

Trong đó:

• J: thông lượng khuếch tán (g/cm²·s)
• D: hệ số khuếch tán phân tử (cm²/s hay m²/s)
• ∂C/∂x: gradient nồng độ theo chiều x (g/cm³ per cm)
• Dấu âm: khuếch tán xảy ra ngược chiều gradient (từ cao xuống thấp)

Ý nghĩa vật lý: Khuếch tán phân tử là hệ quả của chuyển động nhiệt hỗn loạn (Brownian motion) của các phân tử. Các phân tử chất tan di chuyển ngẫu nhiên theo mọi hướng, nhưng do vùng nồng độ cao có nhiều phân tử hơn, thống kê xác suất cho thấy dòng net phân tử luôn từ nồng độ cao về nồng độ thấp.

1.2 Giá Trị D Trong Môi Trường Nước Thải

GS. Lê Văn Cát nhấn mạnh: hệ số khuếch tán phân tử D trong môi trường nước rất nhỏ:

D = 10⁻⁶ − 10⁻⁵ cm²/s = 10⁻¹⁰ − 10⁻⁹ m²/s

Để có cảm giác về tầm quan trọng của giá trị này, so sánh với các môi trường khác:

Môi trường	D điển hình (cm²/s)
Khí (không khí, 25°C)	10⁻¹ - 1
Nước lỏng (25°C)	10⁻⁶ - 10⁻⁵
Màng sinh học (biofilm)	10⁻⁷ - 10⁻⁶
Bông bùn (floc)	10⁻⁸ - 10⁻⁷
Đất bão hòa	10⁻⁹ - 10⁻⁸

D trong nước nhỏ hơn trong khí khoảng 100.000 lần — phản ánh sự cản trở lớn của các phân tử nước xung quanh đối với sự di chuyển của phân tử chất tan.

Phụ thuộc nhiệt độ: D tăng theo nhiệt độ theo phương trình Stokes-Einstein:

D = k_B·T / (3π·μ·d)

Trong đó k_B là hằng số Boltzmann, T là nhiệt độ tuyệt đối, μ là độ nhớt động lực, d là đường kính phân tử.

Ở 20°C → 30°C: D tăng khoảng 25-30% — đáng kể nhưng vẫn trong phạm vi nhỏ so với tốc độ đối lưu điển hình.

Phần 2: Định Luật Fick Thứ Hai — PT 7-15

2.1 Từ Fick Thứ Nhất Đến Fick Thứ Hai

Định luật Fick thứ nhất cho biết tại một thời điểm và vị trí xác định, thông lượng khuếch tán là bao nhiêu. Nhưng để mô tả sự thay đổi nồng độ theo thời gian tại một điểm cố định trong không gian — cần định luật Fick thứ hai.

Áp dụng bảo toàn khối lượng cho một phần tử thể tích dV = A·dx tại vị trí x:

Tốc độ tích lũy = Thông lượng vào − Thông lượng ra

A·dx·(∂C/∂t) = J(x)·A − J(x+dx)·A = −[J(x+dx) − J(x)]·A

A·dx·(∂C/∂t) = −(∂J/∂x)·A·dx

∂C/∂t = −∂J/∂x

Thay J = −D·∂C/∂x:

∂C/∂t = D · ∂²C/∂x² (PT 7-15)

2.2 Phân Tích Toán Học PT 7-15

PT 7-15 là phương trình khuếch tán (diffusion equation) hay phương trình nhiệt (heat equation) — một trong những phương trình vi phân riêng (PDE) căn bản nhất trong vật lý toán.

Chiều biến thiên:

• ∂²C/∂x² > 0 (hàm lõm lên): ∂C/∂t > 0 → nồng độ tăng theo thời gian (vùng trũng của phân bố nồng độ)
• ∂²C/∂x² < 0 (hàm lõm xuống): ∂C/∂t < 0 → nồng độ giảm theo thời gian (vùng đỉnh của phân bố nồng độ)

Ý nghĩa: PT 7-15 mô tả quá trình san bằng (smoothing) phân bố nồng độ theo thời gian. Bất kỳ phân bố không đều nào cũng dần dần được làm phẳng bởi khuếch tán.

Nghiệm cơ bản (Green's function): Nếu tại t = 0 đặt lượng chất M vào một điểm x = 0 trong môi trường vô hạn, nghiệm của PT 7-15 là phân bố Gaussian:

C(x,t) = M/(2A·√(π·D·t)) · e^(-x²/(4·D·t))

Phân bố này lan rộng theo thời gian với độ lệch chuẩn σ = √(2·D·t). Đây là nền tảng để hiểu hiện tượng phân tán (dispersion) trong dòng chảy.

2.3 Chiều Dài Khuếch Tán Đặc Trưng L_diff

Từ nghiệm Gaussian, chiều dài khuếch tán đặc trưng trong thời gian t:

L_diff = √(2·D·t)

Với D = 10⁻⁵ cm²/s (nước):

Thời gian t	L_diff (khuếch tán)
1 giây	0,0045 cm = 0,045 mm
1 phút	0,035 cm = 0,35 mm
1 giờ	0,27 cm = 2,7 mm
1 ngày	1,3 cm = 13 mm
1 tuần	3,5 cm = 35 mm

Trong 8 giờ HRT điển hình của bể aerotank, khuếch tán phân tử chỉ vận chuyển chất trong phạm vi khoảng 6-9 mm — hoàn toàn không đáng kể so với chiều dài bể hàng chục mét.

Phần 3: Phương Trình Đối Lưu-Khuếch Tán PT 7-16

3.1 Tổng Hợp Hai Cơ Chế Vận Chuyển

Khi xét đồng thời cả đối lưu (convection) và khuếch tán phân tử (molecular diffusion) mà không có phản ứng hóa học, phương trình bảo toàn khối lượng cho một phần tử thể tích dV = A·dx:

Tích lũy = Vào − Ra (do đối lưu) + Vào − Ra (do khuếch tán)

• Đối lưu mang vào: v·C(x)·A dt
• Đối lưu mang ra: v·C(x+dx)·A dt → Net đối lưu = −v·(∂C/∂x)·A·dx·dt
• Khuếch tán vào: J(x)·A·dt = −D·(∂C/∂x)|_x · A·dt
• Khuếch tán ra: J(x+dx)·A·dt → Net khuếch tán = +D·(∂²C/∂x²)·A·dx·dt

Kết hợp:

A·dx·dC = [−v·(∂C/∂x) + D·(∂²C/∂x²)]·A·dx·dt

Chia cho A·dx·dt:

−∂C/∂t = v·∂C/∂x − D·∂²C/∂x² (PT 7-16, dấu âm bên trái biểu diễn suy giảm nồng độ nguyên liệu)

3.2 Phân Tích Vai Trò Từng Số Hạng

Số hạng đối lưu v·∂C/∂x: Mô tả sự vận chuyển nồng độ bởi dòng chảy. Nếu ∂C/∂x < 0 (nồng độ giảm theo chiều dòng chảy — trường hợp điển hình), số hạng này âm → đóng góp vào việc giảm nồng độ tại điểm xét theo thời gian.

Số hạng khuếch tán D·∂²C/∂x²: Mô tả tác động san bằng của khuếch tán. Luôn có xu hướng làm phẳng gradient nồng độ — giảm bớt sự chênh lệch giữa đầu vào (nồng độ cao) và đầu ra (nồng độ thấp).

3.3 Số Péclet — Đánh Giá Tương Đối Hai Cơ Chế

Tỷ lệ giữa vận chuyển đối lưu và khuếch tán được đặc trưng bởi số Péclet (Pe):

Pe = v·L / D

Trong đó:

• v: vận tốc dòng chảy tuyến tính (m/s)
• L: chiều dài đặc trưng (chiều dài bình, m)
• D: hệ số khuếch tán phân tử (m²/s)

Khi Pe >> 1: Đối lưu thống trị — khuếch tán không đáng kể → Dòng đẩy lý tưởng (Plug Flow)

Khi Pe << 1: Khuếch tán thống trị — bình phản ứng hành xử như bể khuếch tán đồng nhất

Khi Pe ≈ 1: Cả hai cơ chế quan trọng → cần giải PT 7-16 đầy đủ

Ví dụ tính toán — kênh oxy hóa dài 50m tại TP.HCM:

• v = 0,3 m/s (vận tốc điển hình kênh oxy hóa Carrousel)
• L = 50 m
• D = 10⁻⁹ m²/s (trong nước ở 25°C)
• Pe = 0,3 × 50 / 10⁻⁹ = 1,5 × 10¹⁰

Pe cực kỳ lớn → khuếch tán phân tử hoàn toàn không đáng kể. Điều này xác nhận quan sát của GS. Lê Văn Cát: "D rất nhỏ, sự đóng góp của khuếch tán không lớn, có thể bỏ qua."

Phần 4: Phân Tán Trục (Axial Dispersion) — Vượt Ra Ngoài Khuếch Tán Phân Tử

4.1 Sự Khác Biệt Giữa Khuếch Tán Phân Tử Và Phân Tán

Trong khi khuếch tán phân tử (molecular diffusion) có D rất nhỏ và có thể bỏ qua, phân tán trục (axial dispersion) trong dòng chảy thực là hiện tượng khác và không thể bỏ qua:

Khuếch tán phân tử (PT 7-15): Do chuyển động nhiệt ngẫu nhiên của từng phân tử. D = 10⁻⁹ - 10⁻⁸ m²/s trong nước.

Phân tán trục: Do sự kết hợp của khuếch tán phân tử với sự phân bố vận tốc dòng chảy (velocity profile). Trong ống có dòng tầng (laminar), vận tốc tại tâm ống lớn hơn tại thành ống → các phần tử chất tan ở các vị trí bán kính khác nhau di chuyển với tốc độ khác nhau → "kéo dài" phân bố nồng độ theo chiều dọc.

Hệ số phân tán trục D_L >> D_phân tử:

• Dòng tầng trong ống (Taylor dispersion): D_L = v²·d²/(192·D_m) ∝ v²/D_m — tỷ lệ nghịch với D_m!
• Dòng rối: D_L ≈ (10-100) × v × d (thực nghiệm)

Ví dụ: Ống nhựa d = 5 cm, v = 0,1 m/s, D_m = 10⁻⁹ m²/s: D_L = (0,1)²·(0,05)²/(192·10⁻⁹) = 0,01·0,0025/(1,92×10⁻⁷) ≈ 1,3 × 10⁻⁴ m²/s

D_L lớn hơn D_m khoảng 100.000 lần! Phân tán trục không thể bỏ qua trong thiết kế kênh dài.

4.2 Số Péclet Trục (Axial Péclet Number)

Thực tế thiết kế dùng số Péclet trục Pe_L dựa trên D_L thay vì D_m:

Pe_L = v·L / D_L

Đây là thông số quyết định xem bình phản ứng có hành xử gần với dòng đẩy lý tưởng hay không:

• Pe_L > 100: dòng đẩy lý tưởng (sai số < 5%)
• 10 < Pe_L < 100: phân tán đáng kể, cần hiệu chỉnh
• Pe_L < 10: phân tán mạnh, hành xử gần CSTR

4.3 Ứng Dụng: Đánh Giá Thủy Lực Kênh Xử Lý Tại Việt Nam

Kênh oxy hóa Carrousel (oxidation ditch) tại trạm XLNT Đà Lạt: hình elip, chiều dài hiệu dụng ~120 m, v ≈ 0,3 m/s, D_L ≈ 0,05 m²/s: Pe_L = 0,3 × 120 / 0,05 = 720 → Dòng đẩy tốt, có thể dùng PT 7-21

Bể aerotank dạng vuông 20×20 m với sục khí phân tán: v ≈ 0,05 m/s (dòng tuần hoàn), D_L cao do khuấy trộn mạnh ≈ 0,5 m²/s: Pe_L ≈ 0,05 × 20 / 0,5 = 2 → Gần CSTR, phải dùng PT 7-14 thay vì PT 7-21

Phần 5: Điều Kiện Biên Của PT 7-16 — Từ Toán Học Đến Thực Tiễn Thiết Kế

5.1 Điều Kiện Biên Danckwerts

PT 7-16 là PDE cần điều kiện biên (BC) tại cả hai đầu của bình phản ứng. Điều kiện biên Danckwerts (1953) được chấp nhận rộng rãi:

Đầu vào (x = 0): v·C₀ = v·C(0⁺) − D·∂C/∂x|_(x=0⁺)

Điều kiện này nói rằng: thông lượng tổng (đối lưu + khuếch tán) tại đầu vào bằng v·C₀ (thông lượng của dòng nạp).

Đầu ra (x = L): ∂C/∂x|_(x=L) = 0

Điều kiện này nói rằng: tại đầu ra, không có khuếch tán ngược chiều ra khỏi bình — chất chỉ rời khỏi bình bằng đối lưu.

5.2 Nghiệm Phân Tích Trong Trường Hợp Đơn Giản

Với Pe_L >> 1 (bỏ qua khuếch tán), PT 7-16 đơn giản hóa thành:

∂C/∂t + v·∂C/∂x = 0

Đây là phương trình sóng thuần túy với nghiệm:

C(x,t) = C₀(t − x/v) — Sóng nồng độ truyền với vận tốc v, không biến dạng

Một lượng cơ chất nạp vào đầu vào tại t = 0 sẽ xuất hiện nguyên vẹn ở đầu ra tại t = L/v = θ (thời gian lưu thủy lực). Đây là đặc trưng lý tưởng của PFR — mỗi "plug" (phần tử chất lỏng) di chuyển như một khối thống nhất từ đầu vào đến đầu ra.

5.3 Phân Tích RTD (Residence Time Distribution) Của PFR Lý Tưởng

Phân phối thời gian lưu (RTD) của PFR lý tưởng (Pe_L → ∞) là hàm delta Dirac:

E(t) = δ(t − θ)

Tất cả phần tử chất lỏng lưu lại đúng θ giây, không sớm không muộn. Đây là phân bố hẹp nhất có thể — không có "ngắn mạch" (short-circuiting) và không có "vùng chết" (dead zones).

Thực nghiệm xác định RTD bằng thí nghiệm dấu vết (tracer test):

• Bơm xung (impulse) dung dịch NaCl hoặc thuốc nhuộm rhodamine vào đầu vào tại t = 0
• Đo nồng độ tại đầu ra theo thời gian
• PFR lý tưởng → xuất hiện peak sắc nét tại t = θ
• PFR thực → peak tù hơn do phân tán (độ tù ∝ 1/Pe_L)
• CSTR → đường cong mũ giảm dần E(t) = (1/θ)·e^(-t/θ)

Phần 6: Ứng Dụng Thực Tiễn PT 7-15, 7-16 Trong Đánh Giá Kênh Xử Lý Tại Việt Nam

6.1 Kênh Oxy Hóa — Đánh Giá Tính Xấp Xỉ Dòng Đẩy

Kênh oxy hóa (oxidation ditch) là hệ thống phổ biến tại nhiều đô thị nhỏ và vừa ở Việt Nam (Đà Lạt, Bảo Lộc, Hội An, Tam Kỳ). Việc đánh giá xem kênh có hành xử gần dòng đẩy lý tưởng hay gần CSTR cần thiết cho thiết kế chính xác.

Thông số điển hình kênh oxy hóa:

• Chiều dài toàn kênh: L_tổng = 200 m (vòng elip, chu vi 80 m × 2)
• Vận tốc dòng: v = 0,3 m/s (đảm bảo bùn không lắng)
• Thời gian lưu thủy lực: θ = V/Q = 8 giờ = 28.800 giây

Dòng chảy trong kênh oxy hóa là tuần hoàn nhiều vòng: số vòng tuần hoàn trong thời gian θ: N_vòng = v·θ/L_tổng = 0,3 × 28.800 / 200 = 43 vòng

Với 43 vòng tuần hoàn, mỗi "plug" nước thải đi qua toàn bộ kênh 43 lần trước khi ra. Hiệu quả khuếch tán trục tích lũy qua nhiều vòng làm cho Pe_L hiệu dụng giảm đáng kể. Thực tế, kênh oxy hóa hành xử giống CSTR hơn là PFR — đây là lý do cần tính toán bằng PT 7-14 (CSTR) thay vì PT 7-21 (PFR) cho loại công trình này.

6.2 Mương Oxy Hóa Thẳng — Gần PFR Hơn

Ngược lại, mương oxy hóa thẳng (straight oxidation channel) hoặc hệ thống AO/AOB hình chữ nhật dài hành xử gần PFR hơn:

• Chiều dài L = 50 m, chiều rộng 5 m, sâu 4 m
• Thể tích V = 50 × 5 × 4 = 1.000 m³, Q = 2.000 m³/ngày = 83 m³/h
• θ = V/Q = 12 giờ
• v = Q/(5×4) = 83/(20×3600) = 0,00115 m/s = 0,115 cm/s

Với v thấp này, D_L (Taylor dispersion) nhỏ, Pe_L cao → hành xử gần PFR. Đây là thiết kế phù hợp khi muốn tận dụng lợi thế thể tích nhỏ của PFR so với CSTR.

6.3 Bể UASB — Dòng Đẩy Từ Dưới Lên

Bể UASB (Upflow Anaerobic Sludge Blanket) phổ biến trong XLNT công nghiệp thực phẩm tại Bình Dương và Tiền Giang là ví dụ thực tế của PFR theo phương đứng:

• Nước thải vào từ đáy, chảy lên đều qua lớp bùn hạt (granular sludge)
• Phân bố vận tốc: đồng đều theo mặt cắt ngang (do bộ phân phối đều)
• Gradient nồng độ theo chiều đứng: COD cao ở đáy, thấp ở đỉnh
• Áp dụng PT 7-16 với chiều x là chiều đứng từ dưới lên

Số Péclet trục của UASB: v = Q/(A) ≈ 1-5 m/h lên (vận tốc bề mặt), D_L ≈ 10⁻⁵ - 10⁻⁴ m²/s trong lớp bùn hạt, L = chiều cao bể 3-6 m:

Pe_L = v·L/D_L = (1/3600) × 4 / 10⁻⁴ ≈ 11

Pe_L ≈ 11 — phân tán đáng kể nhưng vẫn gần PFR hơn CSTR. Trong thiết kế UASB, thường dùng mô hình PFR có hiệu chỉnh phân tán (dispersion model) thay vì PFR lý tưởng hoàn toàn.

Kết Luận

Định luật Fick thứ hai (PT 7-15) và phương trình đối lưu-khuếch tán (PT 7-16) cung cấp nền tảng lý thuyết đầy đủ để hiểu sự phân bố nồng độ trong bình phản ứng dòng chảy. Hệ số khuếch tán phân tử D rất nhỏ trong nước (10⁻⁹ - 10⁻⁸ m²/s) và số Péclet Pe = v·L/D cực kỳ lớn trong mọi thiết bị phản ứng thực tế — điều này cho phép bỏ qua khuếch tán phân tử và đơn giản hóa phương trình trong hầu hết ứng dụng XLNT.

Tuy nhiên, phân tán trục (D_L >> D_m) không thể bỏ qua trong đánh giá chính xác hiệu suất PFR thực. Số Péclet trục Pe_L = v·L/D_L quyết định mức độ xấp xỉ dòng đẩy lý tưởng đạt được trong thực tế — thông số quan trọng mà kỹ sư thiết kế và vận hành hệ thống XLNT tại Việt Nam cần nắm vững để phân biệt khi nào áp dụng PT 7-21 (PFR) và khi nào cần dùng PT 7-14 (CSTR).