Tổng quan về luận án

Trong cơ học chất lỏng tính toán (Computational Fluid Dynamics - CFD) và kỹ thuật vi lỏng (Microfluidics), hiện tượng tương tác giữa các cấu trúc giọt lỏng với bề mặt kênh dẫn đóng vai trò cốt lõi trong việc điều khiển dòng chảy sinh hóa ở quy mô micro và nano. Luận án tiến sĩ chuyên ngành Cơ học chất lỏng và chất khí (Mã số: 9 44 01 08) với đề tài "Nghiên cứu tương tác của hạt chất lỏng với bề mặt rắn dưới tác động của mao dẫn nhiệt", do nghiên cứu sinh Nguyễn Đình Hòe thực hiện tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam, dưới sự hướng dẫn khoa học của PGS. Vũ Văn Trường và TS. Nguyễn Hồng Phan (bảo vệ năm 2023), là công trình tiên phong giải quyết bài toán động lực học của hạt chất lỏng đa lớp (compound droplets / double emulsions) dưới tác động đồng thời của lực căng bề mặt phụ thuộc nhiệt độ và cấu trúc hình học của biên rắn.

Khoảng trống nghiên cứu (Research Gap) then chốt mà luận án nhắm tới là sự thiếu vắng các mô hình động lực học hoàn chỉnh cho cấu trúc giọt đa lớp chịu hiệu ứng mao dẫn nhiệt (thermocapillarity / Marangoni effect) khi tương tác với bề mặt rắn. Trong khi các nghiên cứu quốc tế kinh điển như công trình của Sui (2014) hay Fath & Bothe (2015) chỉ tập trung vào hạt chất lỏng đơn pha (single droplet) trên bề mặt đẳng nhiệt hoặc gradient nhiệt phẳng, thì hành vi thủy động lực học của hệ đa pha lồng ghép (hạt lõi bên trong - vỏ chất lỏng bên ngoài) khi va chạm hoặc chuyển động dọc bề mặt rắn dính ướt/không dính ướt vẫn chưa được giải mã tường minh. Cụ thể, các câu hỏi nghiên cứu (Research Questions) và giả thuyết cốt lõi (Hypotheses) được xác lập:

  1. RQ1: Sự hiện diện của hạt chất lỏng bên trong và tỷ lệ bán kính lõi - vỏ ($R_{io} = R_i / R_o$) điều biến như thế nào đến ứng suất biến dạng và độ lõm bề mặt phía sau của giọt khi đi qua kênh vi lỏng thắt hẹp hình sin?
    • H1: Lực cản thủy động và gradient áp suất cục bộ qua đoạn thắt tỷ lệ thuận với số mao dẫn ($Ca$), kích hoạt hiện tượng phân tách hoặc phá vỡ cấu trúc giọt ngoài khi $R_{io}$ vượt ngưỡng giới hạn.
  2. RQ2: Gradient nhiệt độ theo phương thẳng đứng tác động như thế nào đến khả năng nảy và trạng thái cân bằng động của hạt chất lỏng đa lớp trên bề mặt siêu kỵ nước (không dính ướt)?
    • H2: Lực mao dẫn nhiệt hướng tâm sinh ra do gradient nhiệt cục bộ làm giảm áp suất bề mặt dưới của hạt, tạo lực nâng nhiệt hỗ trợ hạt nảy lên và duy trì trạng thái lơ lửng ổn định tại khoảng cách cân bằng $H_s$.
  3. RQ3: Trên bề mặt rắn dính ướt có gradient nhiệt nằm ngang, sự tương tác giữa trường xoáy Marangoni nội tại và góc tiếp xúc động ($\theta_D$) quyết định xu hướng chuyển động của hạt đa lớp theo chiều hướng nóng hay hướng lạnh?
    • H3: Chiều chuyển động không chỉ bị chi phối bởi góc tiếp xúc tĩnh ($\theta_e$) mà còn có thể bị đảo ngược hoàn toàn khi thay đổi tỷ số bán kính hạt lõi và tỷ số độ nhớt giữa các pha.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của hệ phương trình Navier-Stokes đa pha không nén được, phương trình bảo toàn năng lượng với giả thiết Boussinesq cho sức căng bề mặt $\sigma(T) = \sigma_0 - \sigma_T(T - T_0)$, mô hình động lực học góc tiếp xúc động Muradoglu - Tasoglu (2014) cải biên từ tương quan Hoffman - Kistler, và phương pháp số theo dấu biên (Front-Tracking Method - FT). Phạm vi nghiên cứu bao quát các dải thông số không thứ nguyên rộng: số Reynolds ($Re \in [0.1, 100]$), số Weber ($We \in [0.1, 100]$), số Marangoni ($Ma \in [0.1, 50]$), số Bond ($Bo$), số Ohnesorge ($Oh = 0.2$), góc tiếp xúc tĩnh $\theta_e$ biến thiên từ $60^\circ$ đến $150^\circ$, cùng các tỷ số độ nhớt và bán kính đa dạng.

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy động lực học của dòng chảy vi lỏng chứa giọt được chia thành hai nhánh chính: tương tác dạng cản trở (hydrodynamic constriction interaction) và tương tác dạng trực tiếp (direct surface contact interaction).

Trong nhánh tương tác dạng cản trở, Bretherton (1961) đã đặt nền móng lý thuyết với công thức xác định chiều dày màng chất lỏng bôi trơn giữa bọt khí và thành ống mao dẫn thẳng: $$W = 1.29(3\mu U / \sigma)^{2/3}$$ áp dụng khi số mao dẫn rất nhỏ ($\mu U/\sigma < 0.003$). Tiếp đó, Schwartz et al. (1986) mở rộng công thức thành: $$W = \frac{0.643}{\kappa}(3Ca)^{2/3}$$ với độ cong $\kappa = 1/R$ trong dải $Ca \in [10^{-7}, 10^{-4}]$. Về mặt thực nghiệm và mô phỏng số trên giọt đơn lớp, Olbricht & Kung (1992) quan sát thấy bề mặt phía sau của giọt lỏng bị lõm khi $Ca = 0.277$ và bị phân rã hoàn toàn khi $Ca = 4.08$, tỷ số độ nhớt $\lambda = 0.0013$. Tsai et al. (2011) mô phỏng giọt đơn qua kênh thắt hình sin, chỉ ra xu hướng tách hạt ở $\lambda < 0.01$ và $Ca = 0.1$. Đối với hệ giọt đa lớp, Borthakur et al. (2018) đã phân tích biến dạng giọt đôi trong ống thẳng nhưng chỉ dừng lại ở điều kiện đẳng nhiệt; Vu et al. (2020) khảo sát giọt đa lớp trong ống côn nhưng chưa tích hợp hiệu ứng mao dẫn nhiệt.

Trong nhánh tương tác dạng trực tiếp, tranh luận học thuật kéo dài xoay quanh nghịch lý kỳ dị ứng suất cắt tại đường tiếp xúc ba pha (contact line singularity) khi áp dụng điều kiện biên không trượt (no-slip condition). Ford & Nadim (1994) áp dụng thuyết bôi trơn với điều kiện biên trượt Navier để tính toán vận tốc giọt dưới gradient nhiệt: $$U = \frac{\sigma_T G}{\mu} + \frac{\sigma_T G b}{\mu} J - \frac{\sigma_B \cos(\theta_B) - \sigma_A \cos(\theta_A)}{2\mu} J$$ tuy nhiên mô hình này chỉ chuẩn xác với góc tiếp xúc tĩnh nhỏ ($\theta_e < 45^\circ$). Về thực nghiệm, Dai et al. (2017a, 2017b) đã theo dõi chuyển động của giọt dầu hỏa trên bề mặt nhẵn và bề mặt có vi rãnh dưới gradient nhiệt độ $3^\circ\text{C/mm}$, chứng minh giọt giãn ra rồi co lại khi trôi về phía nguồn lạnh.

Về mặt mô phỏng số, tồn tại hai quan điểm đối nghịch rõ rệt:

  • Quan điểm thứ nhất từ Nguyen et al. (2013) và Le et al. (2015) cho rằng trên bề mặt ưa nước ($\theta_e < 90^\circ$), sự mất cân bằng giữa hai cấu trúc xoáy Marangoni nội tại luôn đẩy giọt di chuyển về phía nhiệt độ thấp hơn.
  • Quan điểm thứ hai từ Sui (2014) khi mở rộng góc tiếp xúc đến $150^\circ$ phát hiện rằng khi $\theta_e > 105^\circ$, giọt đảo chiều chuyển động sang phía có nhiệt độ cao hơn; tuy nhiên nghiên cứu của Sui bị chỉ trích do thiếu dữ liệu kiểm chứng độc lập ở góc tiếp xúc lớn. Fath & Bothe (2015) cải tiến mô hình với góc tiếp xúc động nhưng vẫn chỉ giới hạn ở cấu trúc hạt đơn lớp.

Luận án của NCS Nguyễn Đình Hòe đã định vị chính xác vào giao điểm chưa từng được khai phá: mở rộng toàn diện động lực học tương tác từ hạt đơn lớp sang hạt đa lớp (compound droplets), đồng thời giải quyết triệt để sự tương tác phi tuyến giữa lực mao dẫn nhiệt Marangoni, lực mao dẫn hình học, biến dạng mặt phân cách kép và động lực học góc tiếp xúc trên các loại bề mặt có độ dính ướt đối nghịch nhau.

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

graph TD
    A["Lực mao dẫn nhiệt (Marangoni Stress): sigma(T) = sigma_0 - sigma_T(T - T_0)"] --> D["Hệ phương trình Navier-Stokes Đa pha (Front-Tracking Method)"]
    B["Động lực học góc tiếp xúc (Muradoglu - Tasoglu Model)"] --> D
    C["Cấu trúc hạt đa pha (Pha lõi, Pha vỏ, Pha môi trường)"] --> D
    D --> E["Tương tác dạng cản trở (Kênh thắt hẹp hình sin đẳng nhiệt)"]
    D --> F["Tương tác trực tiếp trên bề mặt kỵ nước (Gradient nhiệt đứng)"]
    D --> G["Tương tác trực tiếp trên bề mặt dính ướt (Gradient nhiệt ngang)"]
    E --> H["Hiện tượng lõm bề mặt sau & Biểu đồ chuyển tiếp trạng thái giọt"]
    F --> I["Hiện tượng nảy & Cân bằng lơ lửng ổn định (Hs)"]
    G --> J["Hiện tượng phân cực xoáy & Đảo chiều chuyển động của giọt ngoài"]

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

Luận án mang lại bước chuyển dịch mô thức (paradigm shift) từ mô hình dòng chảy vi giọt đơn giản sang động lực học cấu trúc màng kép phức hợp, với 4 đóng góp lý thuyết căn bản:

  1. Mở rộng lý thuyết bôi trơn và mao dẫn của Bretherton - Schwartz cho hệ giọt lồng đa lớp: Thiết lập mối quan hệ hàm số giữa độ lõm bề mặt phía sau ($\Delta D_{dimple}$), tỷ số bán kính hạt lõi/vỏ $R_{io} = R_i / R_o$, và tỷ số sức căng bề mặt $\sigma_{io} = \sigma_i / \sigma_o$. Chứng minh rằng hạt lõi bên trong đóng vai trò như một "vật cản thủy động lực học cục bộ", làm biến đổi trường áp suất bên trong hạt ngoài, dẫn đến các chế độ biến dạng phi đối xứng khi vượt qua đoạn kênh thắt.
  2. Khám phá cơ chế nâng nhiệt mao dẫn (Thermocapillary Levitation Mechanism): Luận án trích dẫn trực tiếp cơ chế vật lý: "Sự chênh lệch sức căng bề mặt dọc theo biên phân cách của hạt chất lỏng tạo ra lực mao dẫn nhiệt. Khi lơ lửng bên trong chất lỏng khác, hạt chất lỏng luôn di chuyển từ vùng có nhiệt độ thấp hơn đến vùng có nhiệt cao hơn dưới tác động của lực mao dẫn nhiệt." Tuy nhiên, khi hạt tiếp xúc với bề mặt không dính ướt dưới gradient nhiệt độ đứng, lực mao dẫn nhiệt sinh ra mô-men ứng suất kéo căng giọt lên phía trên, triệt tiêu lực bám dính và trọng trường, tạo ra chiều cao ổn định $H_s$ mà lý thuyết hạt đơn lớp truyền thống không thể dự báo.
  3. Phát triển định luật ghép cặp xoáy Marangoni đa pha: Chứng minh rằng sự hiện diện của hạt lõi phân chia xoáy Marangoni thành hai hệ thống xoáy thứ cấp lồng nhau. Luận án ghi nhận: "Hạt chất lỏng bên trong có thể ảnh hưởng rất lớn đến động lực học của hạt chất lỏng bên ngoài, tương tác giữa hạt chất lỏng đa lớp với bề mặt rắn dưới tác động của mao dẫn nhiệt có thể phức tạp, khó đoán, và khác biệt so với hạt chất lỏng đơn lớp."
  4. Xác lập ranh giới ổn định cấu trúc nhũ tương đa lớp: Xây dựng hệ thống bất đẳng thức giới hạn phá vỡ giọt (emulsion breakup boundary conditions) dựa trên tổ hợp các số không thứ nguyên $Ca_{crit}, Ma_{crit}, We_{crit}$.

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

Khung phân tích của luận án tích hợp liền mạch ba trụ cột lý thuyết:

  • Thủy động lực học Navier-Stokes đa pha: Mô hình hóa tương tác giữa ba chất lỏng bất khả hòa tan (chất lỏng trong 1, chất lỏng trung gian 2, chất lỏng môi trường 3) thông qua một phương trình động lượng duy nhất với trường thuộc tính vật lý trơn hóa: $$\rho(\mathbf{x}) = I_i \rho_i + I_o \rho_m + (1 - I_i - I_o)\rho_o$$ $$\mu(\mathbf{x}) = I_i \mu_i + I_o \mu_m + (1 - I_i - I_o)\mu_o$$ trong đó $I_i, I_o$ là các hàm chỉ thị Heaviside xấp xỉ liên tục.
  • Lý thuyết nhiệt mao dẫn Marangoni: Ứng suất tiếp tuyến mặt phân cách được mô hình hóa chính xác qua lực khối Dirac: $$\mathbf{f}s = \int{\Gamma} \left( \sigma \kappa \mathbf{n} + \nabla_s \sigma \right) \delta(\mathbf{x} - \mathbf{x}_f) ds$$ với $\nabla_s \sigma = -\sigma_T \nabla_s T$.
  • Lý thuyết góc tiếp xúc động (Dynamic Contact Line Theory): Khử hoàn toàn điểm kỳ dị ứng suất tại thành rắn bằng cách áp dụng mô hình Muradoglu & Tasoglu (2014), tính toán góc tiếp xúc động $\theta_D$ thông qua hàm thực nghiệm Hoffman: $$f_{Hoff}(x) = \arccos\left{ 1 - 2\tanh\left[ 5.16 \left( \frac{x}{1 + 1.31x^{0.99}} \right)^{0.706} \right] \right}$$ với vận tốc đường tiếp xúc $V_{cl}$ được trích xuất tại khoảng cách giả định $h_{th} = 2\Delta x$ từ biên rắn.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái Thực chứng luận quy nạp - diễn dịch tính toán (Computational Positivism), kết hợp chặt chẽ giữa xây dựng thuật toán số nguyên thủy (in-house solver), kiểm chứng chéo với các nghiệm giải tích kinh điển và đối chuẩn (benchmarking) với các bộ dữ liệu thực nghiệm quốc tế độ chính xác cao.

Mô hình tính toán được triển khai trên hệ tọa độ Đề-các hai chiều ($2D$) và đối xứng trục ($2D$-axisymmetric), rời rạc hóa trên lưới so le (Staggered Grid / MAC grid). Áp suất $p$, nhiệt độ $T$, các thuộc tính $\rho, \mu$ được lưu trữ tại tâm ô lưới $(i, j)$, trong khi các thành phần vận tốc $u_{i+1/2, j}$ và $v_{i, j+1/2}$ được đặt tại các mặt biên ô lưới nhằm triệt tiêu dao động áp suất răng cưa (checkerboard pressure oscillations).

   (i, j+1) -------- v(i, j+1/2) -------- (i+1, j+1)
      |                                        |
      |                                        |
  u(i-1/2, j)           p, T, rho, mu         u(i+1/2, j)
      |                    (i, j)              |
      |                                        |
    (i, j)  --------  v(i, j-1/2) --------  (i+1, j)

Quy trình nghiên cứu rigorous

  1. Thuật toán theo dấu biên (Front-Tracking Scheme): Mặt phân cách $\Gamma$ được biểu diễn dưới dạng lưới di động Lagrange bao gồm các phần tử đoạn thẳng nối các điểm mốc (marker points). Mật độ điểm mốc được kiểm soát nghiêm ngặt từ 2 đến 4 điểm trên mỗi ô lưới Euler cố định. Quy trình tái cấu trúc biên (restructuring) tự động thêm điểm khi khoảng cách phần tử $\Delta s > 0.8\Delta x$ và xóa điểm khi $\Delta s < 0.2\Delta x$, đảm bảo độ trơn và bảo toàn thể tích tuyệt đối của các pha lỏng.
  2. Quy trình giải tách biến áp suất - vận tốc (Projection Method): Phương trình động lượng được giải qua hai bước:
    • Bước 1: Tính toán trường vận tốc trung gian $\mathbf{u}^$ loại bỏ gradient áp suất bằng lược đồ sai phân tiến bậc 1 theo thời gian kết hợp sai phân trung tâm bậc 2 cho số hạng khuếch tán và đối lưu: $$\frac{\mathbf{u}^ - \mathbf{u}^n}{\Delta t} = -\mathbf{A}^n + \frac{1}{\rho^n}\left(\mathbf{D}^n + \mathbf{f}_s^n\right) + \mathbf{g}$$
    • Bước 2: Giải phương trình Poisson áp suất: $$\nabla \cdot \left( \frac{1}{\rho^n} \nabla p^{n+1} \right) = \frac{\nabla \cdot \mathbf{u}^*}{\Delta t}$$ bằng phương pháp lặp quá giãn liên tiếp (Successive Over-Relaxation - SOR) với tiêu chí hội tụ thắt chặt: sai số tương đối $\epsilon < 10^{-6}$.
    • Bước 3: Cập nhật vận tốc không phân kỳ: $\mathbf{u}^{n+1} = \mathbf{u}^* - \frac{\Delta t}{\rho^n}\nabla p^{n+1}$.
  3. Giải phương trình cân bằng năng lượng: Trường nhiệt độ $T^{n+1}$ được cập nhật tường minh: $$T^{n+1} = T^n + \Delta t \left[ -\mathbf{u}^n \cdot \nabla T^n + \frac{1}{\rho c_p}\nabla \cdot (k \nabla T^n) \right]$$ với điều kiện biên nhiệt loại 1 (Dirichlet) sử dụng các nút nhiệt độ ảo $T_{2, j} = 2T_w - T_{1, j}$ tại vách rắn.
  4. Điều kiện ổn định bước thời gian (CFL & Viscous Limits): $$\Delta t \le \min \left( \frac{h}{2|\mathbf{u}|{\max}}, \frac{\rho{\min} h^2}{4\mu_{\max}}, \sqrt{\frac{(\rho_1+\rho_2)h^3}{4\pi\sigma}} \right)$$

Data và phân tích kiểm chứng (Validation Triangulation)

Mã nguồn được lập trình bằng ngôn ngữ FORTRAN, trải qua 4 bước kiểm chứng độ chính xác nghiêm ngặt:

Nghiên cứu đối chuẩn Dạng bài toán Thông số kỹ thuật kiểm chứng Kết quả & Sai số
Borthakur et al. (2018) Giọt đa lớp di chuyển trong kênh thẳng $Ca = 0.05 - 0.5$, $\lambda = 1.0$, $R_{io} = 0.5$ Trùng khớp biên dạng hạt trong và ngoài tại mọi thời điểm
Qu & Hu (2012) Va chạm giọt đơn trên bề mặt không dính ướt $We = 10 - 80$, $Re = 100 - 500$, $\theta_e = 160^\circ$ Khớp hoàn hảo hệ số giãn rộng cực đại và tỷ lệ hình dạng $D_h(t)$
Young et al. (1959); Kalichetty et al. (2019) Chuyển động mao dẫn nhiệt của giọt trong khe hẹp $Ma = 0.1, 1.0, 10.0$; đối lưu nhiệt dừng Vận tốc hạt tiệm cận chính xác nghiệm giải tích YGB kinh điển
Šikalo et al. (2005) Giọt va chạm bề mặt sáp dính ướt $Re = 36$, $We = 93$, $Eo = 0.285$, $\theta_e = 93^\circ$ Biến thiên bán kính dính ướt $R_w(t)$ tương đồng thực nghiệm; sai số góc tiếp xúc tĩnh $\le 1.3%$

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

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

  1. Hiệu ứng lõm bề mặt kép và sự phá vỡ cấu trúc giọt ngoài trong kênh thắt: Khi hạt đa lớp đi qua kênh thắt hình sin trong điều kiện đẳng nhiệt, bề mặt phía sau của giọt ngoài xuất hiện độ lõm sâu phụ thuộc chặt chẽ vào số mao dẫn $Ca$. Khi $Ca > 0.3$ kết hợp với tỷ số bán kính $R_{io} > 0.6$, hạt lõi bên trong gây ra sự tập trung ứng suất cắt cực đại tại vùng đuôi giọt, dẫn đến hiện tượng hạt trong đâm xuyên phá vỡ màng chất lỏng giọt ngoài (core penetration and outer droplet rupture).
  2. Cơ chế nảy và cân bằng lơ lửng ổn định trên bề mặt siêu kỵ nước: Dưới tác động của gradient nhiệt độ hướng thẳng đứng (nhiệt độ cao ở biên trên $T_{hot}$, bề mặt rắn lạnh $T_{cold}$), hạt chất lỏng đa lớp va chạm với bề mặt không dính ướt không chỉ nảy lên mà còn đạt được trạng thái lơ lửng cân bằng bền tại độ cao ổn định $H_s$. Lực nâng Marangoni sinh ra do biến thiên sức căng bề mặt $\Delta \sigma = \sigma_T (T_{hot} - T_{cold})$ cân bằng hoàn toàn với trọng lực hữu hiệu và lực cản nhớt.
  3. Hiện tượng nghịch lý về hướng chuyển động trên bề mặt dính ướt: Đối với hạt đơn lớp trên bề mặt dính ướt có gradient nhiệt ngang, luận án tái khẳng định: "Khi góc tiếp xúc tĩnh nhỏ hơn 90°, hướng chuyển động của hạt chất lỏng phụ thuộc vào độ lớn của hai xoáy xuất hiện bên trong hạt chất lỏng. Sự chênh lệch độ lớn của hai xoáy này gây ra chuyển động sang phía có nhiệt độ thấp hơn của hạt chất lỏng." Tuy nhiên, với hạt đa lớp, khi thay đổi tỷ số độ nhớt $\mu_{mo}$ và tỷ số bán kính tương đương $R_{io}$, hệ xoáy Marangoni bị phân cực bất đối xứng, cho phép điều khiển giọt di chuyển ngược về phía nhiệt độ cao hơn ngay cả khi $\theta_e < 90^\circ$.
  4. Biểu đồ trạng thái động lực học toàn diện (Regime Phase Diagrams): Xây dựng thành công 3 biểu đồ pha trạng thái:
    • Chế độ biến dạng - lõm - tách hạt trong kênh thắt theo cặp thông số $(Ca, R_{io})$.
    • Chế độ nảy - dao động - lơ lửng trên bề mặt không dính ướt theo cặp $(Ma, We)$.
    • Chế độ trôi trái - trôi phải - đứng yên trên bề mặt dính ướt theo cặp $(\theta_e, \mu_{mo})$.

Implications đa chiều

graph LR
    A["Phát hiện của Luận án"] --> B["Học thuật / Lý thuyết"]
    A --> C["Công nghệ Vi lỏng & Y sinh"]
    A --> D["Sản xuất Công nghiệp"]
    B --> B1["Hoàn thiện mô hình Marangoni đa pha"]
    B --> B2["Chuẩn hóa thuật toán Front-Tracking"]
    C --> C1["Hệ phân tích Lab-on-a-Chip không bơm"]
    C --> C2["Đóng gói & Phóng thích thuốc trúng đích (DDS)"]
    D --> D1["Thiết bị tản nhiệt vi kênh cho vi chip điện tử"]
    D --> D2["Chế tạo vi cầu vật liệu polymer đa lớp"]
  • Về mặt học thuật: Cung cấp cơ sở dữ liệu số chuẩn xác cao cho cơ học chất lỏng đa pha phi đẳng nhiệt, đóng góp trực tiếp vào lý thuyết bôi trơn vi mô và thủy động lực học bề mặt phân cách.
  • Về mặt công nghệ vi lỏng (Microfluidics & Lab-on-a-Chip): Mở ra giải pháp điều khiển hạt lỏng sinh hóa hoàn toàn thụ động (passive manipulation) bằng nhiệt độ, loại bỏ nhu cầu sử dụng vi van, vi bơm cơ học phức tạp dễ gây biến tính mẫu sinh học.
  • Về mặt dược học và y sinh: Tối ưu hóa quy trình tổng hợp hạt nhũ tương kép (W/O/W hoặc O/W/O) mang thuốc trúng đích, kiểm soát chính xác độ dày màng bao bọc để ngăn ngừa vỡ thuốc sớm khi đi qua các mao mạch hẹp trong cơ thể.
  • Về mặt công nghệ tản nhiệt vi điện tử: Ứng dụng hiện tượng nảy và tự dịch chuyển của giọt dưới gradient nhiệt để thiết kế các bề mặt làm mát tự thích ứng (self-cooling thermal management systems) cho chip vi xử lý mật độ công suất cao.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu nội tại mang tính phương pháp luận:

  1. Giới hạn về không gian tính toán ($2D$/Trục đối xứng): Mô hình số hiện tại được triển khai trên hệ tọa độ hai chiều phẳng và đối xứng trục. Mặc dù nắm bắt chuẩn xác các cơ chế vật lý nền tảng, mô hình chưa thể hiện đầy đủ các hiệu ứng mất ổn định không gian ba chiều ($3D$ azimuthal instabilities) như hiện tượng cuộn xoáy 3D dọc chu vi giọt hoặc sự đứt gãy không đối xứng của sợi chất lỏng.
  2. Giới hạn về màng tiền bôi trơn nano (Precursor Film): Phương pháp theo dấu biên lưới vĩ mô chưa tích hợp phương trình màng mỏng phân tử (Disjoining pressure) để mô phỏng lớp màng chất lỏng nano siêu mỏng còn lưu lại phía sau khi giọt chuyển động trên bề mặt thực tế.
  3. Giả thiết tính chất nhiệt - vật lý tuyến tính: Luận án giả thiết độ nhớt $\mu$, khối lượng riêng $\rho$, nhiệt dung $c_p$ và hệ số dẫn nhiệt $k$ không đổi theo nhiệt độ, chỉ có sức căng bề mặt $\sigma(T)$ biến thiên tuyến tính. Trong các môi trường chịu gradient nhiệt cực lớn, tính chất phi tuyến của độ nhớt nhiệt độ có thể gây ra sai số nhỏ.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Mở rộng thuật toán Front-Tracking sang không gian 3 chiều song song hóa hoàn chỉnh trên nền tảng tính toán hiệu năng cao (HPC/GPU).
  • Tích hợp lực tách phân tử (Disjoining Pressure) vào điều kiện biên đường tiếp xúc ba pha để mô tả hiện tượng bám dính ở thang đo nano.
  • Khảo sát ảnh hưởng của trường nhiệt độ biến thiên tuần hoàn theo thời gian thực ($AC$ thermal fields) và bề mặt có cấu trúc siêu kỵ nước phân cấp (hierarchical re-entrant microstructures).
  • Thực hiện các đo đạc thực nghiệm quang học vi dòng (Micro-PIV và Shadowgraphy) trên hệ giọt kép để đối sánh định lượng trường vận tốc thực tế.

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

Luận án của NCS Nguyễn Đình Hòe thể hiện năng lực học thuật xuất sắc và năng suất công bố khoa học vượt trội, tạo ra tầm ảnh hưởng rõ nét trên bình diện quốc tế và quốc gia:

  • Năng suất học thuật đỉnh cao: Các kết quả cốt lõi của luận án đã được công bố trong 04 bài báo khoa học quốc tế uy tín (ISI/Scopus Q1/Q2) với tư cách là tác giả chính (First Author), 01 bài báo quốc tế đồng tác giả, và 01 công trình báo cáo tại Hội nghị Cơ học toàn quốc. Đây là minh chứng thuyết phục nhất cho tính mới, độ tin cậy và giá trị khoa học đạt chuẩn quốc tế của công trình.
  • Tác động tới ngành công nghiệp bán dẫn và MEMS: Cung cấp mô hình tính toán nền tảng cho việc thiết kế các kênh vi lỏng tích hợp trong thiết bị giải trình tự gen thế hệ mới và các hệ vi cơ điện tử (MEMS) ứng dụng mao dẫn nhiệt.
  • Tiềm năng trích dẫn học thuật: Với việc giải quyết bài toán giọt đa lớp phi đẳng nhiệt - một chủ đề đang thu hút sự quan tâm lớn của cộng đồng vật lý chất lưu (Physics of Fluids, Journal of Fluid Mechanics), công trình được dự báo sẽ đạt chỉ số trích dẫn cao trong các nghiên cứu về hệ nhũ tương và công nghệ vi lỏng.

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

  • Nghiên cứu sinh và Giới học thuật (Doctoral & Academic Researchers): Tiếp cận một khung lý thuyết và thuật toán số hoàn chỉnh (Front-Tracking + Dynamic Contact Angle) để tiếp tục phát triển các bài toán dòng chảy đa pha phức hợp.
  • Kỹ sư thiết kế hệ thống Vi lưu (Microfluidic & Bio-MEMS Engineers): Sở hữu bộ thông số thiết kế hình học kênh dẫn (độ thắt, góc côn) và biểu đồ pha vận hành tối ưu nhằm kiểm soát hạt lỏng mà không làm phá vỡ cấu trúc nhũ tương.
  • Chuyên gia R&D Dược phẩm và Công nghệ Thực phẩm: Nắm bắt cơ chế biến dạng và phá vỡ của giọt kép để hoàn thiện công nghệ đóng gói vi nang (microencapsulation) chứa hoạt chất sinh học, tăng thời hạn bảo quản và độ chính xác khi phóng thích thuốc.
  • Các nhà phát triển công nghệ tản nhiệt (Thermal Engineers): Ứng dụng hiện tượng nảy và dịch chuyển nhiệt của giọt lỏng để phát triển các giải pháp tản nhiệt không dùng năng lượng phụ trợ cho thiết bị điện tử siêu nhỏ.

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?

Đóng góp độc đáo nhất là việc mở rộng lý thuyết bôi trơn mao dẫn Bretherton - Schwartz và lý thuyết góc tiếp xúc động Kistler - Muradoglu từ hệ hạt đơn lớp sang hệ hạt lỏng đa lớp (compound droplets) dưới tác động của hiệu ứng mao dẫn nhiệt Marangoni. Luận án đã chứng minh hạt lõi bên trong làm tái phân bố phi đối xứng trường áp suất và cấu trúc xoáy nội tại, tạo ra các quy luật chuyển động và biến dạng hoàn toàn mới lạ so với giọt đơn pha truyền thống.

2. Sự đổi mới trong phương pháp nghiên cứu thể hiện như thế nào khi so sánh với các nghiên cứu trước?

So với nghiên cứu của Sui (2014) sử dụng điều kiện biên trượt Navier đơn giản hay nghiên cứu của Borthakur et al. (2018) chỉ dừng lại ở bài toán đẳng nhiệt trong ống thẳng, luận án đã:

  • Tự phát triển mã nguồn FORTRAN tích hợp thuật toán Front-Tracking xử lý đồng thời nhiều mặt phân cách di động Lagrange với lưới Euler cố định so le.
  • Kết hợp mô hình góc tiếp xúc động phi tuyến Muradoglu - Tasoglu với hàm Hoffman thực nghiệm trên lưới so le giải bằng SOR, xử lý chính xác góc tiếp xúc từ $60^\circ$ đến $150^\circ$ với sai số trạng thái tĩnh $\le 1.3%$.

3. Phát hiện bất ngờ nhất (counter-intuitive) có dữ liệu minh chứng là gì?

Phát hiện bất ngờ nhất là hiện tượng đảo ngược hướng chuyển động của giọt ngoài trên bề mặt dính ướt ($\theta_e < 90^\circ$). Trong khi lý thuyết kinh điển cho rằng giọt ưa nước luôn trôi về phía nhiệt độ thấp, luận án chứng minh rằng việc tăng tỷ số bán kính hạt lõi $R_{io}$ hoặc điều chỉnh tỷ số độ nhớt $\mu_{mo}$ sẽ kích hoạt sự chiếm ưu thế của xoáy nội tại thứ cấp, đẩy toàn bộ cấu trúc giọt đa lớp di chuyển ngược hướng về phía nguồn nhiệt cao hơn.

4. Quy trình tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Có. Luận án cung cấp đầy đủ hệ phương trình toán học chi tiết, các hệ số rời rạc hóa sai phân hữu hạn trên lưới so le (Phương trình 2.34 đến 2.44), thuật toán nội suy lực căng bề mặt Dirac Delta, điều kiện ổn định bước thời gian CFL, cũng như toàn bộ sơ đồ thuật toán tổng quát (Hình 2.12), cho phép tái lập mã nguồn mô phỏng với độ tin cậy tuyệt đối.

5. Định hướng nghiên cứu 10 năm tới được phác thảo ra sao?

Tập trung vào chuyển đổi mô hình sang $3D$ song song hóa quy mô lớn, tích hợp lực tương tác phân tử bề mặt (phương trình màng mỏng) để nghiên cứu tương tác giọt đa lớp trên các bề mặt siêu kỵ nước có cấu trúc vi gai nano (hierarchical nanostructured surfaces) phục vụ công nghệ y sinh cá thể hóa và siêu máy tính.

Kết luận

Luận án tiến sĩ của NCS Nguyễn Đình Hòe là một công trình khoa học mẫu mực, xuất sắc toàn diện cả về chiều sâu học thuật, tính tiên phong lý thuyết và độ chuẩn xác phương pháp luận trong lĩnh vực Cơ học chất lỏng tính toán.

Tóm lược 5 đóng góp đột phá cốt lõi:

  1. Xây dựng thành công bộ chương trình mô phỏng số chuyên dụng (FORTRAN in-house solver) ứng dụng phương pháp theo dấu biên (Front-Tracking) kết hợp mô hình góc tiếp xúc động tiên tiến, giải quyết trọn vẹn bài toán hạt chất lỏng đa lớp tương tác với bề mặt rắn dưới tác động của mao dẫn nhiệt.
  2. Làm sáng tỏ động lực học biến dạng và quy luật hình thành độ lõm bề mặt sau của giọt đa lớp khi di chuyển qua kênh vi lỏng thắt hẹp hình sin, xác lập các ngưỡng tham số không thứ nguyên ($Ca, R_{io}, \sigma_{io}$) dẫn đến sự phá vỡ cấu trúc giọt ngoài.
  3. Phát hiện và giải thích bản chất vật lý của hiện tượng nảy và cân bằng lơ lửng ổn định ($H_s$) của hạt đa lớp trên bề mặt không dính ướt dưới tác động của gradient nhiệt độ phương thẳng đứng.
  4. Khám phá cơ chế phân cực xoáy Marangoni nội tại và sự đảo chiều chuyển động của giọt đa lớp trên bề mặt dính ướt, chứng minh vai trò quyết định của kích thước hạt lõi trong việc điều hướng dòng chảy vi lỏng.
  5. Thiết lập hệ thống biểu đồ pha trạng thái động lực học hoàn chỉnh, cung cấp cơ sở dữ liệu và công cụ dự báo chuẩn xác phục vụ thiết kế các vi mạch Lab-on-a-Chip và công nghệ vận chuyển thuốc trúng đích.

Công trình không chỉ mở ra 3 hướng nghiên cứu mới (vi lưu phi đẳng nhiệt 3D, thao tác nhiệt vi hạt sinh học, màng bôi trơn nano đa pha) mà còn khẳng định bước tiến vượt bậc của nền khoa học cơ học tính toán Việt Nam trên bản đồ học thuật quốc tế.