Một vận tốc lắng chưa phải hiệu suất tách
Trong tính sơ bộ bình tách dầu–nước, công thức Stokes thường được đưa vào bảng tính ngay sau đường kính giọt. Tuy nhiên, kết quả chỉ là vận tốc tương đối của một đối tượng lý tưởng với chất lỏng xung quanh. Để suy liệu giọt có tới mặt phân cách trước cửa ra hay không, còn cần hướng dòng, quãng đường, phân bố kích thước và mô hình tiếp xúc bề mặt.
Bài này dựng các ví dụ giả định để sàng lọc tính nhất quán. Không lấy một đường kính duy nhất làm cam kết chất lượng nước hoặc dầu sau tách. Ta xét pha phân tán nặng hơn pha liên tục; đối với giọt nhẹ nổi lên, đổi chiều chuyển động và dùng độ lớn chênh khối lượng riêng phù hợp.
Bắt đầu từ lực, xác định đúng đường kính
Với một cầu đường kính d chuyển động chậm trong chất lỏng Newton, lực cản Stokes có độ lớn F_D=3πμdv. Trọng lượng trừ lực nổi bằng (ρ_d−ρ_c)gπd³/6. Khi vận tốc đã ổn định, hai lực cân bằng, suy ra:
v_t = (ρ_d−ρ_c)g d²/(18μ).
ρ_d và ρ_c lần lượt là khối lượng riêng pha phân tán và liên tục, kg/m³; μ là độ nhớt động lực của pha liên tục, Pa·s; d bằng m; g bằng m/s²; v_t bằng m/s. Không đưa độ nhớt động học, đơn vị m²/s, vào vị trí μ mà thiếu nhân khối lượng riêng.
Tài liệu ASPECT trình bày cân bằng này theo bán kính r, với v_t=2Δρgr²/(9μ), đồng thời nêu khác biệt giữa cầu cứng, miền vô hạn và mô hình số có biên hữu hạn. [1] Hai cách viết bằng nhau khi d=2r. Nhập bán kính vào công thức đường kính làm vận tốc nhỏ đi bốn lần; đó là lỗi định nghĩa biến, không phải hệ số bảo thủ có chủ ý.
Giọt lỏng có thể khác cầu cứng
COMSOL phân biệt mô hình Stokes với Hadamard–Rybczynski cho giọt hoặc bọt rất sạch: chuyển động bên trong và trạng thái bề mặt làm thay đổi lực cản. Với chất hoạt động bề mặt, ứng xử có thể gần Stokes hơn. [2] Vì vậy ta chỉ dùng xấp xỉ cầu cứng hoặc bề mặt gần bất động cho bài toán dưới đây; không khẳng định mọi giọt nước trong dầu đều đáp ứng giả thiết ấy.
Giả thiết bổ sung là giọt gần hình cầu, hệ đủ loãng để bỏ tương tác giữa các giọt, không vỡ hay nhập giọt trên đường xét, không sát thành và chưa có dao động rối làm chi phối quỹ đạo. Nếu những điều này không được dữ liệu hỗ trợ, phép tính vẫn có ích để xác định đại lượng cần đo nhưng không đủ xác lập năng lực bình.
Ba kích thước trong một chất lỏng giả định
Cho ρ_c=850 kg/m³, ρ_d=1.000 kg/m³, μ=0,020 Pa·s và g=9,81 m/s². Chênh khối lượng riêng chỉ bằng 150 kg/m³; dùng riêng 1.000 kg/m³ ở tử sẽ bỏ lực nổi và làm sai rất lớn.
Với d=100 µm, tức 0,000100 m, v_t=0,000040875 m/s. Để lắng xuống 0,30 m trong chất lỏng đứng yên cần 7.339,449541 s, khoảng 122,32 phút. Với d=300 µm, v_t=0,000367875 m/s và thời gian còn 815,494393 s, khoảng 13,59 phút. Với d=500 µm, v_t=0,001021875 m/s, thời gian 293,577982 s.
Tăng đường kính ba lần đã tăng vận tốc chín lần trong mô hình. Vì vậy giá trị trung bình kích thước không mô tả được tất cả phần giọt nhỏ. Nếu phần lớn thể tích nằm ở giọt lớn nhưng một số lượng lớn giọt nhỏ còn đi theo dầu ra, việc báo cáo một đường kính trung bình mà không nêu cách lấy trọng số có thể làm diễn giải sai kết quả.
Kiểm tra Reynolds bằng chính vận tốc vừa tính
Đại lượng kiểm tra là Re_d=ρ_c v_t d/μ, dùng vận tốc tương đối chứ không dùng tùy tiện tốc độ trung bình của dòng trong ống cấp. Ba trường hợp lần lượt cho Re_d=0,000173719; 0,004690406 và 0,021714844. Chúng đều nhỏ so với 1, phù hợp với việc dùng Stokes như xấp xỉ dòng chảy nhớt rất chậm.
Đây là phép kiểm sau tính. Không thể gán sẵn Re nhỏ chỉ vì giọt được gọi là “micromet”. Nếu giảm độ nhớt, tăng đường kính hoặc tăng chênh khối lượng riêng, vận tốc và Reynolds cùng đổi. Với các tính chất giữ cố định, thay v_t vào Re_d cho thấy Re_d tỷ lệ d³; tăng đường kính hai lần làm Reynolds tăng tám lần.
Thêm dòng ngang và dòng đi lên
Xét giọt 300 µm, cần đi xuống 0,30 m. Nếu vận tốc ngang đều u_x=0,006 m/s, không có chuyển động thẳng đứng của pha liên tục, giọt đi ngang x=u_x×0,30/v_t=4,892966 m trước khi chạm mặt phân cách. Đây là chiều dài quỹ đạo lý tưởng của một giọt xuất phát tại vị trí đã chọn, chưa cộng vùng phân phối dòng hay vật cản.
Trong kịch bản khác, pha liên tục chuyển động đều đi lên với u_z=0,000200 m/s. Tốc độ xuống tương đối so với thành bình là v_t−u_z=0,000167875 m/s; thời gian đi 0,30 m tăng thành 1.787,043931 s. Dùng v_t trực tiếp sẽ bỏ gần một nửa thời gian cần thiết trong ví dụ này.
Đặt v_t=u_z cho đường kính tới hạn trong mô hình: d_c=√[18μu_z/(Δρg)]=221,200505 µm. Giọt nhỏ hơn giá trị này bị dòng đi lên cuốn lên nếu các giả thiết khác còn giữ. Ngay tại d_c, thời gian lắng theo phương đứng không hữu hạn. Do đó vùng gần cân bằng hai vận tốc rất nhạy, không thể dùng một hệ số hiệu chỉnh thời gian cố định cho mọi cỡ giọt.
Một trường hợp Stokes tự bác bỏ chính nó
Giữ hai khối lượng riêng và d=500 µm, nhưng đổi μ thành 0,001 Pa·s. Công thức Stokes cho v_S=0,0204375 m/s, kéo theo Re_S=8,685938. Điều kiện Reynolds rất nhỏ đã mất, nên không tiếp tục dùng vận tốc này như kết quả hợp lệ.
Để minh họa việc đổi luật cản cho cầu, dùng Schiller–Naumann trong tài liệu COMSOL: C_D=(24/Re)[1+0,15Re^0,687]. [2] Cân bằng lực dẫn tới phương trình v[1+0,15(ρ_cvd/μ)^0,687]=v_S. Giải lặp được v=0,013615522 m/s, Re=5,786597. Vận tốc Stokes cao hơn kết quả hiệu chỉnh khoảng 50,10%. Kết quả mới vẫn chỉ thuộc giả thiết cầu đã chọn; đổi luật cản không tự giải quyết biến dạng hoặc động học bề mặt giọt.
Độ nhạy và dữ liệu quyết định
Trong miền Stokes, vi phân tương đối là dv_t/v_t=dΔρ/Δρ+2dd/d−dμ/μ. Tăng đường kính 20% làm vận tốc tăng chính xác 44%, còn tăng độ nhớt gấp đôi làm vận tốc giảm một nửa. Nếu hai khối lượng riêng gần nhau, sai số nhỏ của từng phép đo có thể tạo sai số tương đối lớn trong chênh Δρ.
Nghiên cứu của Whorton và cộng sự đo hạt lắng trong chất lỏng giảm nhớt theo tốc độ cắt và có tính đàn nhớt; các công thức Stokes hiệu chỉnh thành thông thường dự đoán kém trong tập thí nghiệm đó. [3] Kết quả này hỗ trợ việc kiểm tra tính chất lưu biến trước khi áp dụng, không cung cấp hệ số sẵn dùng cho bình dầu–nước đang xét.
Hồ sơ cần bổ sung gồm phân bố cỡ giọt tại đúng vị trí vào, tính chất ở nhiệt độ vận hành, trạng thái bề mặt, trường vận tốc và khả năng biến đổi giọt trong bình. Khi đã có chúng, mô hình quỹ đạo mới có cơ sở để nối với phân bố thời gian lưu và kết quả thử tách. Phép tính một giọt giúp kiểm logic vật lý và phát hiện độ nhạy; hiệu suất toàn thiết bị còn đòi hỏi tích hợp trên toàn bộ dòng phân tán.
Nguồn tham khảo
[1] The “Stokes’ law” benchmark — ASPECT 3.1.0 — Juliane Dannberg; ASPECT and Computational Infrastructure for Geodynamics. Truy cập ngày 09/11/2026.
[2] Particle Motion in a Fluid — COMSOL 6.4 Particle Tracing Module — COMSOL AB. Truy cập ngày 09/11/2026.
[3] Particle settling in a shear-thinning, viscoelastic fluid in the presence of wall effects — Jodie Whorton, Thomas Jones and Kelly Russell; Scientific Reports 15, article 4482 (2025), Lancaster University research directory. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Trong tính sơ bộ bình tách dầu–nước, công thức Stokes thường được đưa vào bảng tính ngay sau đường kính giọt. Tuy nhiên, kết quả chỉ là vận tốc tương đối của một đối tượng lý tưởng với chất lỏng xung quanh. Để suy liệu giọt có tới mặt phân cách trước cửa ra hay không, còn cần hướng dòng, quãng đường, phân bố kích thước và mô hình tiếp xúc bề mặt.
Bài này dựng các ví dụ giả định để sàng lọc tính nhất quán. Không lấy một đường kính duy nhất làm cam kết chất lượng nước hoặc dầu sau tách. Ta xét pha phân tán nặng hơn pha liên tục; đối với giọt nhẹ nổi lên, đổi chiều chuyển động và dùng độ lớn chênh khối lượng riêng phù hợp.
Bắt đầu từ lực, xác định đúng đường kính
Với một cầu đường kính d chuyển động chậm trong chất lỏng Newton, lực cản Stokes có độ lớn F_D=3πμdv. Trọng lượng trừ lực nổi bằng (ρ_d−ρ_c)gπd³/6. Khi vận tốc đã ổn định, hai lực cân bằng, suy ra:
v_t = (ρ_d−ρ_c)g d²/(18μ).
ρ_d và ρ_c lần lượt là khối lượng riêng pha phân tán và liên tục, kg/m³; μ là độ nhớt động lực của pha liên tục, Pa·s; d bằng m; g bằng m/s²; v_t bằng m/s. Không đưa độ nhớt động học, đơn vị m²/s, vào vị trí μ mà thiếu nhân khối lượng riêng.
Tài liệu ASPECT trình bày cân bằng này theo bán kính r, với v_t=2Δρgr²/(9μ), đồng thời nêu khác biệt giữa cầu cứng, miền vô hạn và mô hình số có biên hữu hạn. [1] Hai cách viết bằng nhau khi d=2r. Nhập bán kính vào công thức đường kính làm vận tốc nhỏ đi bốn lần; đó là lỗi định nghĩa biến, không phải hệ số bảo thủ có chủ ý.
Giọt lỏng có thể khác cầu cứng
COMSOL phân biệt mô hình Stokes với Hadamard–Rybczynski cho giọt hoặc bọt rất sạch: chuyển động bên trong và trạng thái bề mặt làm thay đổi lực cản. Với chất hoạt động bề mặt, ứng xử có thể gần Stokes hơn. [2] Vì vậy ta chỉ dùng xấp xỉ cầu cứng hoặc bề mặt gần bất động cho bài toán dưới đây; không khẳng định mọi giọt nước trong dầu đều đáp ứng giả thiết ấy.
Giả thiết bổ sung là giọt gần hình cầu, hệ đủ loãng để bỏ tương tác giữa các giọt, không vỡ hay nhập giọt trên đường xét, không sát thành và chưa có dao động rối làm chi phối quỹ đạo. Nếu những điều này không được dữ liệu hỗ trợ, phép tính vẫn có ích để xác định đại lượng cần đo nhưng không đủ xác lập năng lực bình.
Ba kích thước trong một chất lỏng giả định
Cho ρ_c=850 kg/m³, ρ_d=1.000 kg/m³, μ=0,020 Pa·s và g=9,81 m/s². Chênh khối lượng riêng chỉ bằng 150 kg/m³; dùng riêng 1.000 kg/m³ ở tử sẽ bỏ lực nổi và làm sai rất lớn.
Với d=100 µm, tức 0,000100 m, v_t=0,000040875 m/s. Để lắng xuống 0,30 m trong chất lỏng đứng yên cần 7.339,449541 s, khoảng 122,32 phút. Với d=300 µm, v_t=0,000367875 m/s và thời gian còn 815,494393 s, khoảng 13,59 phút. Với d=500 µm, v_t=0,001021875 m/s, thời gian 293,577982 s.
Tăng đường kính ba lần đã tăng vận tốc chín lần trong mô hình. Vì vậy giá trị trung bình kích thước không mô tả được tất cả phần giọt nhỏ. Nếu phần lớn thể tích nằm ở giọt lớn nhưng một số lượng lớn giọt nhỏ còn đi theo dầu ra, việc báo cáo một đường kính trung bình mà không nêu cách lấy trọng số có thể làm diễn giải sai kết quả.
Kiểm tra Reynolds bằng chính vận tốc vừa tính
Đại lượng kiểm tra là Re_d=ρ_c v_t d/μ, dùng vận tốc tương đối chứ không dùng tùy tiện tốc độ trung bình của dòng trong ống cấp. Ba trường hợp lần lượt cho Re_d=0,000173719; 0,004690406 và 0,021714844. Chúng đều nhỏ so với 1, phù hợp với việc dùng Stokes như xấp xỉ dòng chảy nhớt rất chậm.
Đây là phép kiểm sau tính. Không thể gán sẵn Re nhỏ chỉ vì giọt được gọi là “micromet”. Nếu giảm độ nhớt, tăng đường kính hoặc tăng chênh khối lượng riêng, vận tốc và Reynolds cùng đổi. Với các tính chất giữ cố định, thay v_t vào Re_d cho thấy Re_d tỷ lệ d³; tăng đường kính hai lần làm Reynolds tăng tám lần.
Thêm dòng ngang và dòng đi lên
Xét giọt 300 µm, cần đi xuống 0,30 m. Nếu vận tốc ngang đều u_x=0,006 m/s, không có chuyển động thẳng đứng của pha liên tục, giọt đi ngang x=u_x×0,30/v_t=4,892966 m trước khi chạm mặt phân cách. Đây là chiều dài quỹ đạo lý tưởng của một giọt xuất phát tại vị trí đã chọn, chưa cộng vùng phân phối dòng hay vật cản.
Trong kịch bản khác, pha liên tục chuyển động đều đi lên với u_z=0,000200 m/s. Tốc độ xuống tương đối so với thành bình là v_t−u_z=0,000167875 m/s; thời gian đi 0,30 m tăng thành 1.787,043931 s. Dùng v_t trực tiếp sẽ bỏ gần một nửa thời gian cần thiết trong ví dụ này.
Đặt v_t=u_z cho đường kính tới hạn trong mô hình: d_c=√[18μu_z/(Δρg)]=221,200505 µm. Giọt nhỏ hơn giá trị này bị dòng đi lên cuốn lên nếu các giả thiết khác còn giữ. Ngay tại d_c, thời gian lắng theo phương đứng không hữu hạn. Do đó vùng gần cân bằng hai vận tốc rất nhạy, không thể dùng một hệ số hiệu chỉnh thời gian cố định cho mọi cỡ giọt.
Một trường hợp Stokes tự bác bỏ chính nó
Giữ hai khối lượng riêng và d=500 µm, nhưng đổi μ thành 0,001 Pa·s. Công thức Stokes cho v_S=0,0204375 m/s, kéo theo Re_S=8,685938. Điều kiện Reynolds rất nhỏ đã mất, nên không tiếp tục dùng vận tốc này như kết quả hợp lệ.
Để minh họa việc đổi luật cản cho cầu, dùng Schiller–Naumann trong tài liệu COMSOL: C_D=(24/Re)[1+0,15Re^0,687]. [2] Cân bằng lực dẫn tới phương trình v[1+0,15(ρ_cvd/μ)^0,687]=v_S. Giải lặp được v=0,013615522 m/s, Re=5,786597. Vận tốc Stokes cao hơn kết quả hiệu chỉnh khoảng 50,10%. Kết quả mới vẫn chỉ thuộc giả thiết cầu đã chọn; đổi luật cản không tự giải quyết biến dạng hoặc động học bề mặt giọt.
Độ nhạy và dữ liệu quyết định
Trong miền Stokes, vi phân tương đối là dv_t/v_t=dΔρ/Δρ+2dd/d−dμ/μ. Tăng đường kính 20% làm vận tốc tăng chính xác 44%, còn tăng độ nhớt gấp đôi làm vận tốc giảm một nửa. Nếu hai khối lượng riêng gần nhau, sai số nhỏ của từng phép đo có thể tạo sai số tương đối lớn trong chênh Δρ.
Nghiên cứu của Whorton và cộng sự đo hạt lắng trong chất lỏng giảm nhớt theo tốc độ cắt và có tính đàn nhớt; các công thức Stokes hiệu chỉnh thành thông thường dự đoán kém trong tập thí nghiệm đó. [3] Kết quả này hỗ trợ việc kiểm tra tính chất lưu biến trước khi áp dụng, không cung cấp hệ số sẵn dùng cho bình dầu–nước đang xét.
Hồ sơ cần bổ sung gồm phân bố cỡ giọt tại đúng vị trí vào, tính chất ở nhiệt độ vận hành, trạng thái bề mặt, trường vận tốc và khả năng biến đổi giọt trong bình. Khi đã có chúng, mô hình quỹ đạo mới có cơ sở để nối với phân bố thời gian lưu và kết quả thử tách. Phép tính một giọt giúp kiểm logic vật lý và phát hiện độ nhạy; hiệu suất toàn thiết bị còn đòi hỏi tích hợp trên toàn bộ dòng phân tán.
Nguồn tham khảo
[1] The “Stokes’ law” benchmark — ASPECT 3.1.0 — Juliane Dannberg; ASPECT and Computational Infrastructure for Geodynamics. Truy cập ngày 09/11/2026.
[2] Particle Motion in a Fluid — COMSOL 6.4 Particle Tracing Module — COMSOL AB. Truy cập ngày 09/11/2026.
[3] Particle settling in a shear-thinning, viscoelastic fluid in the presence of wall effects — Jodie Whorton, Thomas Jones and Kelly Russell; Scientific Reports 15, article 4482 (2025), Lancaster University research directory. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
