Tổng quan về luận án

Sự cố ô nhiễm dầu trên biển và các thủy vực nội địa là hiểm họa môi trường nghiêm trọng bắt nguồn từ hoạt động khai thác hàng hải, dàn khoan ngoài khơi và các cơ sở công nghiệp ven bờ. Theo số liệu thống kê trong luận án, "tổng trọng tải của đội tàu dầu hiện nay khoảng 270 triệu DWT (Dead weight), chiếm tỷ trọng khoảng 30% tổng trọng tải của đội tàu biển thế giới" và "số lượng dầu mỏ thấm qua đất và lan truyền ra nước biển ước tính 3 triệu tấn mỗi năm". Đồng thời, quá trình xả nước dằn tàu ước tính phát tán hơn 300.000 tấn dầu khoáng mỗi năm ra đại dương. Mặc dù Tổ chức Hàng hải Quốc tế (IMO) áp dụng quy chuẩn nghiêm ngặt MARPOL 73/78 (nồng độ dầu xả thải không vượt quá 15 ppm), việc xử lý các dòng hỗn hợp dầu - nước lưu lượng lớn trong các vụ tràn dầu khẩn cấp vẫn đối mặt với giới hạn kỹ thuật nghiêm trọng.

                    +-------------------------------------------------------------+
                    |     KHÔNG GIAN DÒNG CHẢY HAI PHA TRONG THIẾT BỊ LY TÂM       |
                    +-------------------------------------------------------------+
                                                   |
                     [Vào: Hỗn hợp Dầu - Nước (alpha = 300 - 600 ppm)]
                                                   |
                                                   v
                         +-----------------------------------+
                         | Cụm Bầu & 9 Cánh Dẫn (Beta = 13°) |
                         | Tăng tốc tiếp tuyến và tạo xoáy   |
                         +-----------------------------------+
                                                   |
                                                   v
            +-----------------------------------------------------------------+
            |   Ống quay ly tâm (Do = 100 mm, Lo = 2000 mm, n = 5500-6000 rpm)|
            +-----------------------------------------------------------------+
                       /                                             \
                      / (Lực ly tâm cưỡng bức)                        \ (Lực hướng tâm pha nhẹ)
                     v                                                 v
      +-----------------------------+                   +-----------------------------+
      | Pha Nặng: Nước (sát vách)   |                   | Pha Nhẹ: Dầu (trục tâm)     |
      | rho_n = 1067.8 kg/m3        |                   | rho_d = 881 kg/m3           |
      +-----------------------------+                   +-----------------------------+
                     |                                                 |
                     v                                                 v
      [Thoát ra ngoài qua vách ống]                     [Thu gom qua ống tâm d = 50 mm]

Khoảng trống nghiên cứu (Research Gap) cốt lõi xuất phát từ sự bất cập của các công nghệ hiện hành: phương pháp lắng đọng trọng trường đòi hỏi dung tích bồn chứa khổng lồ do vận tốc lắng quá nhỏ ($v_g \sim 10^{-4}\text{ m/s}$); màng lọc polymer SQS-1 hay vật liệu kỵ nước siêu việt dễ bị tắc nghẽn bề mặt ở lưu lượng cao; máy phân ly đĩa nón (100–110 đĩa, khe hở 0,6–0,8 mm, tốc độ 7.000–8.000 vòng/phút) chỉ thích hợp tinh lọc nhiên liệu có tỷ lệ dầu cao chứ không thể xử lý nước nhiễm dầu lưu lượng lớn; trong khi các ống xoáy lốc tĩnh (hydrocyclone) bị suy giảm xoáy nhanh dọc trục và bất ổn định vùng lõi chảy ngược khi lưu lượng biến thiên ($20\div60\text{ m}^3/\text{h}$).

Để giải quyết triệt để bài toán này, luận án tiến hành giải quyết hai câu hỏi nghiên cứu trung tâm:

  • Câu hỏi nghiên cứu 1 ($Q_1$): Quy luật thủy động lực học nào chi phối sự hình thành và phân tách của giọt dầu kích thước micro ($10\div100\ \mu\text{m}$) trong trường ly tâm cưỡng bức tạo bởi ống quay đồng bộ với cụm cánh - bầu dẫn dòng?
  • Câu hỏi nghiên cứu 2 ($Q_2$): Tập hợp tối ưu các thông số hình học (đường kính bầu $D_b$, số cánh $Z$, góc đặt cánh $\beta$, biên dạng Hermite) và chế độ vận hành (tốc độ quay $n$, nồng độ đầu vào $\alpha$) để đạt hiệu suất tách dầu cao nhất là gì?

Hệ thống giả thuyết khoa học được kiểm chứng bao gồm:

  • Giả thuyết 1 ($H_1$): Việc tích hợp ống quay ly tâm tốc độ cao ($5.500\div6.000\text{ vòng/phút}$) với cụm cánh xoắn biên dạng thủy động học tối ưu sẽ duy trì trường vận tốc xoáy tiếp tuyến ổn định dọc toàn bộ chiều dài ống ($L_o = 2.000\text{ mm}$), ngăn chặn hiện tượng suy thoái xoáy thường gặp ở thiết bị tĩnh.
  • Giả thuyết 2 ($H_2$): Sự chênh lệch khối lượng riêng giữa nước ngọt/nước biển ($\rho_n = 1.067,8\text{ kg/m}^3$) và dầu diesel ($\rho_d = 881\text{ kg/m}^3$) dưới gia tốc ly tâm cực đại sẽ rút ngắn thời gian di chuyển hướng tâm ($\tau$) của giọt dầu $50\ \mu\text{m}$ về vùng lõi trục ($r_d \to 0$) xuống dưới 1 giây.

Khung lý thuyết của nghiên cứu dựa trên sự mở rộng định luật bồi lắng Stokes trong hệ quy chiếu quay, tích hợp với hệ phương trình Navier-Stokes đa pha và lý thuyết thủy lực cánh dẫn dòng chảy xoáy. Nghiên cứu thực hiện trên quy mô mô hình thực nghiệm với ống phân ly đường kính trong $D_o = 100\text{ mm}$, chiều dài $L_o = 2.000\text{ mm}$, khảo sát nồng độ dầu đầu vào từ 300 ppm đến 600 ppm ở dải lưu lượng công nghiệp, mang lại giải pháp công nghệ có khả năng nội địa hóa cao cho ngành cơ khí động lực và bảo vệ môi trường biển Việt Nam.

Literature Review và Positioning

Lịch sử phát triển của công nghệ phân ly chất lỏng hai pha bằng trường xoáy đã trải qua nhiều giai đoạn với các đóng góp học thuật quan trọng:

+----------------------------------------------------------------------------------------------------+
|                                 TIẾN TRÌNH NGHIÊN CỨU DÒNG XOÁY PHÂN TÁCH                         |
+----------------------------------------------------------------------------------------------------+
  [1965-1969] Kreith, Sonju, Backshall, Landis: Nghiên cứu băng xoắn, quy luật phân rã xoáy (Re ~ 10^4)
       |
       v
  [1974-1991] Baker & Sayre (1974), Kitoh (1991): Cánh cố định, đo trường vận tốc xoáy bằng LDV & Hot-wire
       |
       v
  [1989-1998] Nejad (1989), Parchen & Steenbergen (1998): Động học phân tán hạt, đo dòng hai pha
       |
       v
  [2021] Luận án Nguyễn Việt Đức: Tích hợp Ống quay + Cụm bầu cánh quay đồng bộ (n = 6000 rpm)
         --> Xóa bỏ lớp biên ma sát vách, duy trì xoáy cưỡng bức liên tục dọc trục Lo = 2000 mm.
+----------------------------------------------------------------------------------------------------+

Giai đoạn đầu tập trung vào các phần tử tạo xoáy thụ động trong ống cố định. Kreith và Sonju (1965) cùng Backshall và Landis (1969) đã khảo sát sự phân bố vận tốc dọc trục và vận tốc tiếp tuyến qua dải băng xoắn ở dải chuẩn số Reynolds $Re_D = 10.000 \div 60.000$, chỉ ra sự phân rã nhanh chóng của động năng xoáy do ma sát nhớt tại vách ống. Tiếp đó, Kitoh (1991) cùng Baker và Sayre (1974) sử dụng kỹ thuật đo gió dây nóng (Hot-wire anemometry) và đầu dò áp suất Pitot để đánh giá trường xoáy qua cánh dẫn cố định ở $Re_D = 12.500 \div 40.000$. Young và Rao (1978), Nejad và cộng sự (1989), Parchen và Steenbergen (1998) đã nâng tầm phương pháp nghiên cứu thông qua kỹ thuật đo vận tốc Doppler Laser hai thành phần (LDV) ở chuẩn số Reynolds cao lên tới $Re_D = 125.000$, xác nhận cấu trúc xoáy Rankin hỗn loạn và hiện tượng mất ổn định tâm xoáy (vortex core precessing).

Một cuộc tranh luận học thuật sâu sắc tồn tại giữa hai trường phái:

  • Trường phái 1 (Thiết bị xoáy tĩnh / Hydrocyclone): Đại diện bởi các nghiên cứu của Colman, Thew và cộng sự. Trường phái này ủng hộ việc sử dụng ống xoáy lốc tĩnh hoặc cánh cố định vì ưu thế không có chi tiết chuyển động cơ học, độ tin cậy cơ khí cao. Tuy nhiên, điểm yếu cốt tử là lực ly tâm bị tiêu tán mạnh theo chiều dài ống; gradient vận tốc cắt lớn tại lớp biên tĩnh gây vỡ các giọt dầu thành các giọt siêu nhỏ ($d_d < 10\ \mu\text{m}$), gây nhũ hóa thứ cấp và làm sụp đổ dòng lõi ngược khi nồng độ dầu dao động.
  • Trường phái 2 (Thiết bị tạo xoáy cưỡng bức bằng cơ năng quay): Đại diện bởi Rocklage-Marliani, Weske, Sturov (1974), Imao, Itoh, Harada (1989) và Kocherscheidt, Schmidts (1996) với bó ống quay ở chuẩn số $Re_D = 218.000 \div 280.000$. Nhóm này chứng minh rằng việc truyền chuyển động quay trực tiếp cho thành ống giúp thiết lập trường xoáy đồng đều, triệt tiêu ứng suất cắt tại vách và duy trì lực ly tâm trên toàn bộ chiều dài làm việc.

So sánh với hai nghiên cứu quốc tế tiêu biểu làm đối chứng trực tiếp trong luận án:

  1. Nghiên cứu của nhóm tác giả sử dụng cánh tạo xoáy cố định trong hỗn hợp dầu diesel - nước muối ($Re_D = 1,81 \times 10^5, S_w = 6,72$): Nghiên cứu này ghi nhận hạt dầu $50\ \mu\text{m}$ mất khoảng cách di chuyển $z = 0,6\text{ m}$ mới tiến gần tâm $0,01\text{ m}$ trong điều kiện không phân tán hỗn loạn, nhưng khi xuất hiện dao động dòng, giọt dầu bị trôi dạt ra vách và hiệu suất tách sụt giảm nghiêm trọng do không có năng lượng bổ sung duy trì xoáy.
  2. Nghiên cứu của Đại học Michigan về màng lọc kỵ dầu cải tiến: Sử dụng cơ chế thấm ướt chọn lọc với cấu trúc silsesquioxan florodexyl liên kết chéo. Dù đạt hiệu quả phân tách cao ở trạng thái tĩnh, màng lọc này bị giới hạn nghiêm trọng về lưu lượng đỉnh (tối đa $250\text{ m}^3/\text{giờ}$ trên $1\text{ m}^2$ bề mặt nhưng nhanh chóng suy giảm khi dầu có độ nhớt cao bám dính).

Luận án của NCS Nguyễn Việt Đức đã định vị chính xác khoảng trống: kết hợp ưu điểm của cơ cấu cánh dẫn định hình dòng xoáy với công nghệ ống quay cưỡng bức đồng bộ tốc độ cao ($5.500\div6.000\text{ rpm}$). Thiết kế này xóa bỏ hoàn toàn lớp biên ma sát tĩnh, triệt tiêu sự hình thành dòng xoáy thứ cấp có hại, đảm bảo dòng chảy ổn định ở lưu lượng lớn phục vụ ứng phó sự cố tràn dầu.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng nền tảng lý thuyết phân ly đa pha trong cơ học chất lưu thông qua việc mô hình hóa chuyển động của giọt chất lỏng nhẹ phân tán trong trường lực ly tâm cực đại. Xuất phát từ cân bằng động học giữa lực ly tâm cưỡng bức và lực cản thủy động Stokes, luận án làm sáng tỏ công thức xác định vận tốc tách hướng tâm $v_{lt}$ của pha phân tán:

$$v_{lt} = \frac{d_d^2 \omega^2 r_d (\rho_n - \rho_d)}{18 \mu_n}$$

Trích dẫn công thức trực tiếp từ tài liệu luận án (Công thức 1.4):

"Vận tốc tách của dầu ra khỏi nước nhờ lực ly tâm: $v_{lt} = \frac{d_d^2 \omega^2 r_d (\rho_n - \rho_d)}{18 \mu_n}$, trong đó $v_{lt}$ là vận tốc tách của dầu ra khỏi nước (m/s); $\omega$ là vận tốc góc (rad/s); $r_d$ là bán kính xuyên tâm quay của giọt dầu (m); $\rho_n, \rho_d$ lần lượt là khối lượng riêng của nước và dầu $(\text{kg/m}^3)$; $\mu_n$ là độ nhớt động lực của nước $(\text{Pa}\cdot\text{s})$."

Đóng góp lý thuyết của nghiên cứu thể hiện qua việc chứng minh rằng: trong hệ quy chiếu quay đồng bộ (ống quay cùng vận tốc góc $\omega$ với dòng chất lỏng), vận tốc trượt hướng tâm của giọt dầu $v_r = \frac{dr_d}{dt} = -v_{lt}$ là hàm số bậc nhất theo bán kính $r_d$ và bậc hai theo tần số góc $\omega$. Nghiên cứu đã xác lập quan hệ giải tích về thời gian lưu chuyển $\tau$ để giọt dầu di chuyển từ thành ống ($R = D_o/2$) vào bán kính lõi thu gom ($R_{core} = d/2$):

$$\tau = \int_{R}^{R_{core}} \frac{dr_d}{-v_{lt}} = \frac{18 \mu_n}{d_d^2 \omega^2 (\rho_n - \rho_d)} \ln\left(\frac{D_o}{d}\right)$$

Phương trình này đóng vai trò dẫn đường để tính toán chiều dài tối thiểu của ống quay $L_o = V_0 \cdot \tau$ (với $V_0$ là vận tốc chảy dọc trục), chứng minh sự phụ thuộc có tính quy luật của quá trình tách pha vào tổ hợp thông số hình học và chế độ động học.

                      KHUNG PHÂN TÍCH TỔNG HỢP CỦA LUẬN ÁN
+-----------------------------------------------------------------------------------+
| 1. Lý thuyết Thủy lực Cánh dẫn & Máy Thủy lực Thể tích                            |
|    - Thiết lập số vòng quay đặc trưng ns = n.sqrt(Q) / H^(3/4)                    |
|    - Tỷ số đường kính bầu d_b = Db/Do = f(kH, ns)                                 |
|    - Độ mau lưới profin (L/T)tb tại tiết diện trung bình                          |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| 2. Hình học Giải tích Tối ưu hóa (Hermite Polynomials)                            |
|    - Nội suy bậc 3 (P1 - P4) cho đường trục nhân profil cánh                      |
|    - Nội suy bậc 5 (Q1 - Q6) cho biên dạng đầu/đuôi bầu quay thuôn mượt           |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| 3. Cơ học Chất lưu Tính toán Đa pha (Multiphase CFD Modeling)                      |
|    - Hệ phương trình liên tục và Navier-Stokes phân pha                           |
|    - Bộ lọc không gian Fourier / Gaussian / Tophat trong mô phỏng dòng xoáy      |
|    - Đánh giá trường nồng độ alpha(r, z) và trường áp suất p(r, z)                |
+-----------------------------------------------------------------------------------+

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp ba trụ cột khoa học:

  1. Lý thuyết Thủy lực cánh dẫn (Turbomachinery Cascade Theory): Ứng dụng quan hệ thực nghiệm giữa số vòng quay đặc trưng $n_s = \frac{3,65 n \sqrt{Q}}{H^{3/4}}$, hệ số cột áp $k_H$, hệ số lưu lượng $k_Q$ để xác định tỷ số đường kính bầu $\bar{d}b = \frac{D_b}{D_o}$, bước cánh $t = \frac{2\pi R{tb}}{Z}$ và độ dày profin $\delta_{max}$.
  2. Kỹ thuật trơn hóa biên dạng bằng toán tử Hermite: Sử dụng các đa thức Hermite bậc ba ($P_1 \div P_4$) để xây dựng đường nhân cánh và đa thức Hermite bậc năm ($Q_1 \div Q_6$) để thiết lập đường bao bầu thuôn mượt ở cả đầu vào ($\Delta z_n$) và đầu ra ($\Delta z_t$). Kỹ thuật này triệt tiêu hoàn toàn các điểm kỳ dị hình học, ngăn ngừa hiện tượng phân tách dòng (flow separation) và giảm thiểu tổn thất động áp lực $\Delta h_h$.
  3. Mô hình hóa số động lực học dòng chảy hai pha (Multiphase CFD): Thiết lập bài toán mô phỏng chuyển động xoáy đối xứng trục trong phần mềm ANSYS Fluent, kết hợp các hàm lọc không gian thực và không gian Fourier (Tophat, Gaussian filters) để giải quyết tương tác giữa pha dầu phân tán và pha nước liên tục.

Điều kiện biên xác định cho khung phân tích bao gồm: đường kính giọt dầu $d_d = 10 \div 100\ \mu\text{m}$, khối lượng riêng pha lỏng $\rho_d = 869 \div 881\text{ kg/m}^3, \rho_n = 1.000 \div 1.067,8\text{ kg/m}^3$, độ nhớt $\mu_d = (16,6 \div 20,5) \times 10^{-3}\text{ Pa}\cdot\text{s}$, nồng độ dầu đầu vào $\alpha \le 1.000\text{ ppm}$, dải vận tốc dòng dọc trục $V_0 = 1,5 \div 2,5\text{ m/s}$.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường nhận thức luận thực chứng (Positivism), kết hợp chặt chẽ giữa tính toán giải tích lý thuyết, mô phỏng số động lực học chất lưu (CFD) và kiểm chứng thực nghiệm vật lý đa cấp độ. Thiết kế nghiên cứu tiến hành tuần tự qua 4 bước: (1) Tính toán xác định kích thước thủy lực hình học cơ bản; (2) Mô hình hóa 3D và mô phỏng số trường thủy động lực học trên ANSYS Fluent; (3) Chế tạo mẫu thử nghiệm vật lý theo các thông số tối ưu từ mô phỏng; (4) Thử nghiệm đối chứng trên hệ thống đo đạc quang học tiêu chuẩn.

+----------------------------------------------------------------------------------------------------+
|                                    QUY TRÌNH NGHIÊN CỨU TOÀN DIỆN                                  |
+----------------------------------------------------------------------------------------------------+
  [BƯỚC 1: TÍNH TOÁN THỦY LỰC GIẢI TÍCH]
  Xác định: Do = 100 mm, Lo = 2000 mm, ns, db = 0.8, Z = 9 cánh, Beta = 13°, c, t, (L/T)tb
                                   |
                                   v
  [BƯỚC 2: MÔ HÌNH HÓA VÀ MÔ PHỎNG SỐ TRÊN ANSYS FLUENT]
  Khảo sát đa thông số: n (5000-6000 rpm), Z (7-10), Beta (11°-19°), Db (60-80 mm), alpha (300-600 ppm)
  --> Đánh giá trường pha, áp suất, vận tốc xoáy, bán kính hội tụ giọt dầu
                                   |
                                   v
  [BƯỚC 3: CHẾ TẠO MẪU VẬT LÝ VÀ LẮP RÁP HỆ THỐNG]
  Ống inox đúc Do = 100 mm, cụm bầu-cánh Z = 9, bộ làm kín cơ khí trục quay, gối đỡ ổ bi
                                   |
                                   v
  [BƯỚC 4: THỰC NGHIỆM ĐỐI CHỨNG VÀ PHÂN TÍCH QUANG HỌC]
  Đo nồng độ dầu bằng máy OMD-21, kiểm soát lưu lượng, đối chiếu kết quả mô phỏng số và thực nghiệm
+----------------------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ các tiêu chuẩn nghiêm ngặt về kỹ thuật thủy khí:

  • Thiết kế hình học và chia lưới CFD: Tạo lập mô hình 3D trên SolidWorks, xây dựng miền tính toán không gian dòng chảy trong ống quay. Lưới tính toán được chia mịn đặc biệt tại các vùng sát bề mặt cánh dẫn, bề mặt bầu và vách ống quay để thu giữ chính xác gradient vận tốc trong lớp biên.
  • Thiết lập điều kiện biên mô phỏng: Biên đầu vào (Velocity Inlet) với phân bố đều vận tốc dọc trục $V_0 = 2,0\text{ m/s}$ và nồng độ thể tích pha dầu; Biên đầu ra (Pressure Outlet) phân tách độc lập thành cửa xả pha dầu tại lõi trung tâm và cửa xả nước sạch tại vành khăn sát vách; Biên vách ống và cụm bầu cánh được gán điều kiện chuyển động quay đồng bộ (Moving Wall / Rotating Frame) với tốc độ góc $\omega = 575,96 \div 628,32\text{ rad/s}$ (tương ứng $5.500 \div 6.000\text{ rpm}$).
  • Hệ thống đo lường thực nghiệm: Sử dụng thiết bị đo nồng độ dầu chuyên dụng OMD-21 (Oil Mist Detector/Oil Content Meter sử dụng nguyên lý tán xạ quang học theo tiêu chuẩn IMO MEPC.107(49)), kết hợp hệ thống đo lưu lượng dòng điện từ và hệ thống lấy mẫu đẳng tốc tại các đầu xả.

Data và phân tích

Các thông số kích thước vật lý và đặc tính lưu chất được chuẩn hóa chính xác trong các bảng dữ liệu của luận án:

  • Kích thước hình học thiết bị chế tạo: Đường kính trong ống quay $D_o = 100\text{ mm}$; Chiều dài làm việc ống quay $L_o = 2.000\text{ mm}$; Đường kính bầu cánh $D_b = 80\text{ mm}$; Số lượng cánh dẫn $Z = 9$; Góc đặt cánh $\beta = 13^\circ$; Đường kính ngoài ống thu dầu bẩn trung tâm $d = 50\text{ mm}$ (bề dày ống $2\text{ mm}$); Ống chế tạo bằng thép không gỉ inox đúc, gối đỡ sử dụng vòng bi tốc độ cao cùng cụm làm kín cơ khí chuyên dụng (Mechanical Face Seal) chống rò rỉ ở áp suất cao.
  • Tính chất vật lý của lưu chất thử nghiệm:
    • Pha liên tục (Nước ngọt pha muối): Khối lượng riêng $\rho_n = 1.067,8\text{ kg/m}^3$; Độ nhớt động lực $\mu_n = 1,183 \times 10^{-3}\text{ Pa}\cdot\text{s}$.
    • Pha phân tán (Dầu Diesel DO thương mại): Khối lượng riêng $\rho_d = 881\text{ kg/m}^3$; Độ nhớt động lực $\mu_d = 19,4 \times 10^{-3}\text{ Pa}\cdot\text{s}$.
    • Nồng độ dầu thử nghiệm tại đầu vào: $\alpha_1 = 300\text{ ppm}$ và $\alpha_2 = 600\text{ ppm}$.

Phát hiện đột phá và implications

Những phát hiện then chốt

Kết quả mô phỏng số chuyên sâu và thực nghiệm đối chứng đã mang lại 5 phát hiện then chốt:

+------------------------------------------------------------------------------------------------------+
|                     TỔNG HỢP CÁC KẾT QUẢ KHẢO SÁT THÔNG SỐ CÔNG NGHỆ TỐI ƯU                         |
+----------------------+--------------------------+----------------------------------------------------+
| Thông số khảo sát    | Dải giá trị thử nghiệm   | Giá trị tối ưu & Cơ chế thủy động học              |
+----------------------+--------------------------+----------------------------------------------------+
| 1. Tốc độ quay (n)   | 5.000 - 6.000 rpm        | 5.500 - 6.000 rpm: Gia tốc ly tâm cực đại, kéo giọt|
|                      |                          | 50 µm về lõi chỉ sau z = 0,6 - 1,7 m.              |
| 2. Số cánh dẫn (Z)   | 7, 8, 9, 10 cánh         | Z = 9 cánh: Đạt phân bố xoáy tối ưu, tránh nghẽn   |
|                      |                          | thủy lực và giảm thiểu tổn thất áp suất.           |
| 3. Góc đặt cánh (Beta)| 11°, 13°, 15°, 17°, 19°  | Beta = 13°: Đảm bảo góc dạt dòng trơn mượt, triệt  |
|                      |                          | tiêu xoáy cuộn phân tách cục bộ tại mép cánh.      |
| 4. Đường kính bầu (Db)| 60 mm, 70 mm, 80 mm      | Db = 80 mm (db = 0,8): Tăng tốc độ dòng chảy qua   |
|                      |                          | khe và định hình lõi xoáy tập trung ở tâm ống.     |
| 5. Nồng độ dầu vào   | 300 ppm & 600 ppm        | Hiệu suất tách đạt đỉnh, nồng độ dầu cửa ra nước   |
|                      |                          | sạch giảm sâu, đạt chuẩn quy chuẩn hàng hải quốc tế|
+----------------------+--------------------------+----------------------------------------------------+
  1. Quy luật ảnh hưởng của vận tốc vòng quay ống ($n$): Khảo sát biến thiên tốc độ từ $5.000$ đến $6.000\text{ rpm}$ cho thấy lực ly tâm tỷ lệ thuận với bình phương tốc độ góc ($\omega^2$). Ở mức $5.500 \div 6.000\text{ rpm}$, các giọt dầu đường kính $50\ \mu\text{m}$ phân tán ngẫu nhiên sát vách ngoài ống ($r = 50\text{ mm}$) bị gia tốc hướng tâm đẩy mạnh vào trục tâm, hoàn tất pha hội tụ tại bán kính thu gom $r \le 25\text{ mm}$ ở khoảng cách dọc trục $z \le 1.700\text{ mm}$.
  2. Quy luật tối ưu số lượng cánh dẫn ($Z$): Thử nghiệm các phương án $Z = 7, 8, 9, 10$ cánh trên bầu chỉ ra rằng: với $Z = 7$ và $8$, động năng xoáy tiếp tuyến truyền vào dòng chất lỏng chưa đủ mạnh, khiến vùng phân ly bị kéo dài vượt quá chiều dài ống ($>2.000\text{ mm}$); ngược lại, với $Z = 10$, tổn thất ma sát cục bộ tăng vọt và xuất hiện hiện tượng nghẹt dòng tại rãnh cánh. Phương án $Z = 9$ cánh tạo ra cấu trúc xoáy đồng nhất và tối ưu hóa chuyển đổi động năng.
  3. Ảnh hưởng phi tuyến của góc đặt cánh ($\beta$): Phân tích các góc $\beta = 11^\circ, 13^\circ, 15^\circ, 17^\circ, 19^\circ$ chứng minh mối quan hệ phi tuyến giữa góc nghiêng cánh và độ ổn định dòng chảy. Tại $\beta = 11^\circ$, góc tạo xoáy quá nhỏ làm giảm gia tốc tiếp tuyến. Tại các góc lớn ($\beta \ge 15^\circ$), xuất hiện hiện tượng tách dòng (flow separation) tại mặt hút của cánh, sinh ra các xoáy thứ cấp làm vỡ giọt dầu thành các hạt mịn ($d_d < 10\ \mu\text{m}$) khó phân ly. Góc $\beta = 13^\circ$ là điểm tối ưu tuyệt đối, giúp dòng chảy đi vào ống quay với độ xoáy cực đại và tổn thất cột áp nhỏ nhất.
  4. Tác động của tỷ số đường kính bầu ($\bar{d}_b$): Khảo sát các kích thước bầu $D_b = 60\text{ mm}, 70\text{ mm}, 80\text{ mm}$ (tương ứng $\bar{d}_b = 0,6; 0,7; 0,8$). Bầu đường kính $D_b = 80\text{ mm}$ làm tăng vận tốc qua tiết diện thu hẹp cục bộ, đẩy nhanh quá trình định hình lõi xoáy tập trung ngay tại đầu vào ống quay.
  5. Tính nhất quán giữa mô phỏng và thực nghiệm: Thử nghiệm đối chứng trên mô hình vật lý với mẫu nước lẫn dầu diesel ở nồng độ $300\text{ ppm}$ và $600\text{ ppm}$ cho thấy sự tương thích xuất sắc giữa dữ liệu đo đạc bằng máy OMD-21 và kết quả tính toán mô phỏng trên ANSYS Fluent (sai số tương đối dưới 8%), khẳng định độ tin cậy của mô hình toán học đã xây dựng.

Implications đa chiều

  • Về mặt học thuật: Cung cấp phương pháp luận hoàn chỉnh tích hợp giữa lý thuyết cánh dẫn máy thủy lực và mô phỏng số động lực học dòng đa pha CFD cho các thiết bị phân ly ly tâm tốc độ siêu cao.
  • Về mặt kỹ thuật công nghệ: Mở ra giải pháp chế tạo thiết bị tách dầu ly tâm dạng ống quay có kích thước nhỏ gọn ($D_o = 100\text{ mm}, L_o = 2.000\text{ mm}$), khối lượng nhẹ, công suất xử lý lớn, vượt trội hoàn toàn so với các bồn lắng trọng lực cồng kềnh.
  • Về mặt ứng dụng thực tiễn: Cung cấp bộ thông số công nghệ chuẩn để các nhà máy đóng tàu và cơ sở cơ khí hàng hải Việt Nam tự chủ thiết kế, chế tạo hệ thống xử lý nước la canh tàu thủy, thiết bị thu gom dầu tràn trên biển và xử lý nước thải nhiễm dầu công nghiệp.
  • Về mặt chính sách và môi trường: Đóng góp công cụ kỹ thuật đắc lực giúp các đội tàu vận tải và cơ sở cảng biển tuân thủ triệt để Công ước quốc tế MARPOL 73/78, giảm thiểu rủi ro xử phạt quốc tế và bảo vệ nguồn lợi thủy hải sản, hệ sinh thái biển ven bờ.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận các giới hạn kỹ thuật và phạm vi hiệu lực của đề tài:

  1. Giới hạn về đặc tính lưu chất: Thử nghiệm mới chỉ tập trung vào hệ dầu Diesel (DO) và nước ngọt/nước muối chuẩn; chưa mở rộng khảo sát đối với các loại dầu thô nặng (Heavy Fuel Oil - HFO), dầu nhờn bôi trơn có độ nhớt cực cao ($\mu > 100\text{ mPa}\cdot\text{s}$) hoặc các hỗn hợp nhũ tương hóa chứa hàm lượng chất hoạt động bề mặt cao.
  2. Giới hạn cơ khí vận hành tốc độ cao: Vận hành ống quay ở tốc độ $6.000\text{ vòng/phút}$ đòi hỏi độ chính xác gia công cơ khí, cân bằng động khắt khe và hệ thống làm mát/bôi trơn liên tục cho cụm làm kín cơ khí (Mechanical Seal), hạn chế khả năng hoạt động liên tục dài ngày trong điều kiện dằn xóc trên tàu biển nhỏ.
  3. Mô hình hóa giọt dầu cố định: Mô phỏng CFD hiện tại giả định đường kính giọt dầu đồng nhất ($d_d = 50\ \mu\text{m}, 100\ \mu\text{m}$), chưa tính toán đến phổ phân bố kích thước hạt ngẫu nhiên (Droplet Size Distribution - DSD) cũng như hiện tượng kết tụ (coalescence) và vỡ hạt (breakup) trong trường cắt hỗn loạn.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) được đề xuất theo 4 hướng:

  • Hướng 1: Nghiên cứu tích hợp mô hình cân bằng quần thể giọt (Population Balance Model - PBM) vào CFD để dự báo chính xác động học kết tụ và phân rã của nhũ tương dầu.
  • Hướng 2: Thử nghiệm với các loại dầu có độ nhớt cao, dầu thô và khảo sát ảnh hưởng của nhiệt độ dòng chảy ($20^\circ\text{C} \div 60^\circ\text{C}$) đến hiệu suất phân ly.
  • Hướng 3: Nghiên cứu tối ưu hóa vật liệu chế tạo (hợp kim titan, composite gia cường sợi carbon) và ứng dụng đệm từ trường (magnetic bearings) nhằm tăng tuổi thọ làm việc của gối trục ở dải vòng quay trên $10.000\text{ rpm}$.
  • Hướng 4: Phát triển hệ thống điều khiển tự động thích ứng (Adaptive Control) tự điều chỉnh tốc độ quay $n$ và khe hở ống thu dầu $d$ theo tín hiệu thời gian thực từ cảm biến nồng độ quang học OMD.

Tác động và ảnh hưởng

Đề tài luận án tiến sĩ của NCS Nguyễn Việt Đức tạo ra những tác động khoa học và kinh tế - xã hội rõ nét:

  • Tác động học thuật: Đặt nền móng cho hướng nghiên cứu thiết bị phân ly chất lỏng đa pha bằng phương pháp ly tâm ống quay tại Việt Nam; tạo tài liệu tham khảo chuyên sâu cho các chương trình đào tạo sau đại học ngành Kỹ thuật Cơ khí Động lực, Khai thác Bảo trì Tàu thủy.
  • Tác động kinh tế và công nghiệp hàng hải: Cho phép các doanh nghiệp đóng tàu (Tổng công ty Công nghiệp Tàu thủy SBIC, các nhà máy đóng tàu quân đội và dân sự) tự chủ sản xuất thiết bị phân ly dầu nước nội địa với giá thành ước tính chỉ bằng 40–50% so với thiết bị nhập khẩu từ châu Âu hay Nhật Bản.
  • Tác động môi trường và xã hội: Nâng cao năng lực ứng phó sự cố tràn dầu khẩn cấp tại các vùng biển nhạy cảm, vịnh kín, khu bảo tồn biển của Việt Nam; kiểm soát ô nhiễm nước thải la canh từ đội tàu biển nội địa, góp phần phát triển bền vững kinh tế biển theo Chiến lược Biển Việt Nam đến năm 2030, tầm nhìn 2045.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giảng viên ngành Kỹ thuật Thủy khí / Cơ khí Động lực: Tiếp cận thuật toán tính toán thủy lực cánh dẫn hoàn chỉnh, mô hình mô phỏng CFD dòng xoáy đa pha và bộ số liệu thực nghiệm đối chứng chi tiết.
  • Kỹ sư Thiết kế Đóng tàu và Công trình Biển: Sở hữu hướng dẫn thiết kế hình học chuẩn xác (biên dạng Hermite, tỷ số bầu, góc đặt cánh) để phát triển các hệ thống xử lý nước la canh, nước dằn tàu đạt chuẩn IMO.
  • Doanh nghiệp Khai thác Cảng biển và Giàn khoan Dầu khí: Ứng dụng giải pháp thiết bị phân ly ly tâm công suất lớn để xử lý nước thải nhiễm dầu ven bờ và nước tách khai thác giàn khoan ngoài khơi với chi phí đầu tư và vận hành tối ưu.
  • Cơ quan Quản lý Nhà nước (Cục Hàng hải Việt Nam, Chi cục Bảo vệ Môi trường Biển): Có thêm căn cứ khoa học và công cụ kỹ thuật để kiểm soát, giám sát và thẩm định công nghệ xử lý chất thải nhiễm dầu trên các phương tiện thủy nội địa và quốc tế.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và mở rộng lý thuyết nào?
Trả lời: Đóng góp lý thuyết độc đáo nhất là việc mở rộng định luật bồi lắng Stokes cổ điển sang trường gia tốc ly tâm cưỡng bức trong hệ quy chiếu ống quay đồng bộ tốc độ cao ($5.500 \div 6.000\text{ rpm}$), tích hợp với lý thuyết cánh dẫn máy thủy lực thông qua việc giải tích hóa biên dạng cánh và bầu bằng toán tử Hermite bậc ba và bậc năm. Công trình đã thiết lập được phương trình giải tích xác định chính xác thời gian và quỹ đạo di chuyển hướng tâm của giọt dầu kích thước $10 \div 100\ \mu\text{m}$ dưới tác dụng đồng thời của lực ly tâm, lực cản nhớt và trường áp suất phân bố dọc trục.

2. Đổi mới về mặt phương pháp nghiên cứu so với các công trình quốc tế trước đây là gì?
Trả lời: So với các nghiên cứu của Kitoh (1991) hay Baker & Sayre (1974) (chỉ khảo sát phần tử tạo xoáy cố định trong ống tĩnh dẫn đến suy giảm xoáy nhanh) và các nghiên cứu ống xoáy lốc của Colman & Thew (bị mất ổn định vùng lõi ngược), luận án đã kết hợp đột phá giữa: (1) Cụm bầu cánh quay dẫn dòng sơ cấp góc đặt $\beta = 13^\circ$ và (2) Ống quay ly tâm thứ cấp duy trì động năng xoáy trên toàn bộ chiều dài $L_o = 2.000\text{ mm}$. Toàn bộ quy trình được kiểm chứng chéo giữa mô phỏng số 3D đa pha ANSYS Fluent và thực nghiệm quang học OMD-21.

3. Phát hiện bất ngờ nhất (counter-intuitive result) thu được từ dữ liệu là gì?
Trả lời: Phát hiện bất ngờ nhất là việc gia tăng góc đặt cánh ($\beta > 15^\circ$) hoặc tăng số lượng cánh dẫn ($Z > 9$) không làm tăng hiệu quả phân ly như suy luận hình học thông thường, mà ngược lại làm sụt giảm nghiêm trọng hiệu suất tách dầu. Nguyên nhân do góc đặt lớn gây hiện tượng tách dòng (separation) và xoáy cục bộ tại mặt sau của cánh, làm vỡ các giọt dầu $50\ \mu\text{m}$ thành các giọt micro mịn ($d_d < 10\ \mu\text{m}$), gây nhũ hóa thứ cấp và cản trở giọt dầu di chuyển về lõi trục.

4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) hoàn chỉnh không?
Trả lời: Có. Luận án cung cấp bộ quy trình chế tạo và thử nghiệm chi tiết từ kích thước hình học chính xác ($D_o = 100\text{ mm}, L_o = 2.000\text{ mm}, D_b = 80\text{ mm}, Z = 9, \beta = 13^\circ$), quy chuẩn ống inox, thông số vòng bi, bản vẽ kết cấu bộ làm kín cơ khí trục quay (Hình 4.9), thuật toán chương trình tính toán thủy lực (Chương 2), điều kiện biên chia lưới CFD (Phụ lục 2) đến quy trình 6 bước tháo lắp và vận hành thử nghiệm hệ thống trong phòng thí nghiệm (Hình 4.10).

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
Trả lời: Lộ trình 10 năm tới bao gồm: (1) Giai đoạn 1-3 năm: Tích hợp mô hình cân bằng quần thể hạt PBM vào CFD và thử nghiệm với các dòng dầu nặng HFO, dầu thô; (2) Giai đoạn 4-6 năm: Nghiên cứu ứng dụng gối đỡ từ trường (active magnetic bearings) và vật liệu siêu nhẹ để nâng dải tốc độ quay lên $10.000 \div 15.000\text{ rpm}$, giảm kích thước thiết bị xuống $50%$; (3) Giai đoạn 7-10 năm: Thương mại hóa các trạm phân ly tự động hóa tích hợp IoT/AI điều khiển thích ứng nồng độ dầu xả thải phục vụ đội tàu biển thông minh và giàn khoan khai thác nước sâu.

Kết luận

Luận án tiến sĩ kỹ thuật của NCS Nguyễn Việt Đức đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 6 đóng góp mang tính bước ngoặt:

  1. Xây dựng thành công thuật toán và chương trình tính toán thủy lực xác định chính xác các thông số hình học cơ bản của thiết bị tách dầu ly tâm dạng ống quay dựa trên số vòng quay đặc trưng $n_s$ và tỷ số đường kính bầu $\bar{d}_b$.
  2. Thiết lập mô hình giải tích biên dạng cánh và bầu quay tối ưu bằng cách ứng dụng đa thức Hermite bậc 3 ($P_1 \div P_4$) và bậc 5 ($Q_1 \div Q_6$), triệt tiêu điểm kỳ dị hình học và tối thiểu hóa tổn thất động áp lực.
  3. Mô hình hóa số hoàn chỉnh quá trình phân ly hai pha lỏng - lỏng trong ống quay ly tâm trên phần mềm ANSYS Fluent, làm sáng tỏ quy luật phân bố trường vận tốc xoáy, áp suất và nồng độ thể tích pha dầu.
  4. Xác định chính xác bộ thông số công nghệ tối ưu của thiết bị: Đường kính trong ống quay $D_o = 100\text{ mm}$, chiều dài $L_o = 2.000\text{ mm}$, đường kính bầu $D_b = 80\text{ mm}$ ($\bar{d}_b = 0,8$), số cánh dẫn $Z = 9$, góc đặt cánh $\beta = 13^\circ$, số vòng quay làm việc $n = 5.500 \div 6.000\text{ vòng/phút}$.
  5. Chế tạo thành công mẫu thiết bị thử nghiệm vật lý hoàn chỉnh tích hợp hệ thống gối đỡ và bộ làm kín cơ khí tốc độ cao chịu áp lực.
  6. Thực nghiệm đối chứng khẳng định độ chính xác cao của mô hình lý thuyết, chứng minh khả năng phân tách triệt để các giọt dầu $50\ \mu\text{m}$ ở nồng độ đầu vào $300 \div 600\text{ ppm}$, mở ra hướng phát triển các sản phẩm cơ khí hàng hải công nghệ cao do Việt Nam hoàn toàn làm chủ thiết kế và chế tạo.