Tổng quan về luận án

Công nghệ vi lưu (microfluidics) và hệ thống vi tổng phân tích ($\mu\text{TAS}$) phát triển từ thập niên 1970 đã mở ra cuộc cách mạng trong y sinh học, chẩn đoán phân tử và tổng hợp vật liệu tiên tiến. Trong các hệ vi lưu kỹ thuật số (digital microfluidics), thao tác điều khiển hạt chất lỏng mà không cần hệ thống bơm van cơ học phức tạp là một thách thức cốt lõi. Việc khai thác lực mao dẫn nhiệt (thermocapillary force / Marangoni effect)—sinh ra từ gradient sức căng bề mặt dọc theo mặt phân cách do chênh lệch nhiệt độ—được xem là giải pháp điều khiển không tiếp xúc đầy tiềm năng. Tuy nhiên, phần lớn các công trình kinh điển chỉ tập trung vào hạt chất lỏng đơn lớp. Trong các ứng dụng dẫn truyền thuốc hướng đích, vi đóng gói thực phẩm và tổng hợp vi hạt, nhũ tương kép hay hạt chất lỏng đa lớp (compound/multi-layer droplets) mới là cấu trúc mang vật chất cốt lõi.

Khoảng trống nghiên cứu then chốt (specific research gap) nằm ở chỗ: các quy luật thủy động lực học của giọt đơn lớp không thể ngoại suy trực tiếp cho giọt đa lớp do sự hiện diện của mặt phân cách thứ cấp bên trong làm biến đổi hoàn toàn trường ứng suất cắt nội tại, gradient áp suất cục bộ và các cấu trúc xoáy Marangoni. Luận án tiến sĩ ngành Cơ học chất lỏng và chất khí (mã số: 9 44 01 08) của Nghiên cứu sinh Nguyễn Đình Hòe, 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 tại Học viện Khoa học và Công nghệ (Viện Hàn lâm Khoa học và Công nghệ Việt Nam), mang tiêu đề "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" đã tiên phong giải quyết khoảng trống học thuật này.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học của luận án được xác định rõ ràng:

  1. Câu hỏi nghiên cứu 1 (RQ1): Sự tương tác hình học và biến dạng cơ học của hạt chất lỏng đa lớp diễn ra như thế nào khi đi qua kênh vi lưu có biên dạng thắt lại (tương tác cản trở) trong điều kiện đẳng nhiệt?
  2. Câu hỏi nghiên cứu 2 (RQ2): Dưới tác động của gradient nhiệt độ theo phương thẳng đứng, động lực học va chạm của hạt chất lỏng đa lớp trên bề mặt rắn kị nước (không dính ướt) có tạo ra trạng thái cân bằng cơ nhiệt mới hay không?
  3. Câu hỏi nghiên cứu 3 (RQ3): Cơ chế nào điều khiển hướng dịch chuyển và vận tốc của hạt chất lỏng đa lớp trên bề mặt rắn dính ướt khi chịu tác động của gradient nhiệt độ theo phương ngang?
  4. Giả thuyết khoa học (H1): Cấu trúc giọt bao bọc bên trong (kích thước, tỷ số độ nhớt, sức căng bề mặt) đóng vai trò điều biến năng lượng then chốt, có khả năng làm đảo chiều chuyển động vĩ mô của toàn bộ giọt đa lớp so với quy luật truyền thống của giọt đơn lớp.

Khung lý thuyết của luận án xây dựng dựa trên sự liên kết chặt chẽ giữa Cơ học chất lưu phi ổn định (Unsteady Hydrodynamics), Nhiệt động lực học bề mặt (Surface Thermodynamics), phương trình bảo toàn động lượng Navier-Stokes kết hợp phương trình bảo toàn năng lượng và mô hình góc tiếp xúc động Kistler - Muradoglu - Tasoglu. Phạm vi nghiên cứu bao quát các chế độ dòng chảy không thứ nguyên với số Reynolds ($Re$), số Weber ($We$), số Mao dẫn ($Ca$), số Marangoni ($Ma$), số Bond ($Bo$), khảo sát góc tiếp xúc tĩnh từ $60^\circ$ đến $150^\circ$, mang lại đóng góp đột phá với 04 bài báo công bố trên các tạp chí quốc tế ISI/Scopus uy tín và 01 công bố đồng tác giả.

Literature Review và Positioning

Lịch sử nghiên cứu thủy động lực học bề mặt phân cách ghi nhận các đóng góp nền tảng từ Young và cộng sự (1959) khi thiết lập công thức lý thuyết giải tích về vận tốc di chuyển ổn định của bọt khí/hạt chất lỏng dưới tác động của gradient nhiệt độ. Tiếp sau đó, Bretherton (1961) và Schwartz và cộng sự (1986) công bố các phương trình xác định bề dày màng chất lỏng bôi trơn giới hạn khi bọt khí di chuyển trong ống mao dẫn thẳng: $W = 1{,}29(3\mu U / \sigma)^{2/3}$ và $W = \frac{0{,}643}{\kappa}(3Ca)^{2/3}$. Đối với tương tác trực tiếp trên bề mặt rắn, Ford và Nadim (1994) áp dụng lý thuyết bôi trơn (lubrication theory) và điều kiện biên trượt để thiết lập quan hệ vận tốc dịch chuyển của giọt đơn lớp nhưng bị giới hạn nghiêm ngặt ở góc tiếp xúc tĩnh nhỏ ($\theta_e < 45^\circ$).

Các tranh luận khoa học sâu sắc xuất hiện khi các nghiên cứu thực nghiệm và mô phỏng số mở rộng sang dải góc tiếp xúc lớn và hình học phức tạp. Dai và cộng sự (2008, 2010) phát hiện quy luật giọt dầu hỏa dãn ra rồi co lại khi chuyển động từ nguồn nhiệt cao sang nhiệt thấp trên bề mặt phẳng và bề mặt có rãnh vi mô. Ngược lại, trong mô phỏng số, Sui (2014) khi khảo sát góc tiếp xúc tĩnh lên đến $150^\circ$ đã chỉ ra trạng thái đứng yên tại $\theta_e = 100^\circ - 105^\circ$, nhưng khi tiếp tục tăng góc tiếp xúc thì hạt lại đổi hướng di chuyển về phía nhiệt độ cao hơn. Fath và Bothe (2015) tiếp tục mở rộng tranh biện khi chứng minh rằng các xoáy đối lưu nhiệt nội tại và sự mất cân bằng lực căng bề mặt tại đường tiếp xúc ba pha là nguyên nhân quyết định hướng chuyển động. Ở hướng tương tác cản trở, Olbricht và Kung (1992), Nath và cộng sự (2017) quan sát thấy hiện tượng lõm bề mặt sau và vỡ hạt ở số mao dẫn lớn ($Ca \sim 1$). Borthakur và cộng sự (2018) và Vu và cộng sự (2020) bắt đầu tiếp cận hạt hai lớp nhưng chỉ dừng lại ở kênh dẫn thẳng hoặc hình côn đẳng nhiệt.

Định vị của luận án này tạo ra bước tiến vượt bậc bằng cách lấp đầy khoảng trống chưa từng được khảo sát trên thế giới: tổng hợp đồng thời tính chất đa lớp của giọt chất lỏng với các điều kiện nhiệt độ không đồng nhất, trên cả ba dạng bài toán: hình học thắt cản trở, va chạm trên bề mặt không dính ướt và di chuyển nhiệt mao dẫn trên bề mặt dính ướt. So với nghiên cứu của Sui (2014) hay Borthakur và cộng sự (2018), luận án đã giải quyết trọn vẹn sự tương tác phi tuyến giữa mặt phân cách kép và bề mặt chất rắn trong trường nhiệt độ biến thiên.

Đó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 và thách thức các lý thuyết cơ bản về mao dẫn nhiệt truyền thống:

  1. Mở rộng lý thuyết nhiệt mao dẫn Marangoni kinh điển: Luận án chứng minh rằng định luật chuyển động truyền thống—cho rằng giọt chất lỏng luôn di chuyển về phía nhiệt độ thấp hơn khi góc dính ướt nhỏ và về phía nhiệt độ cao hơn khi góc dính ướt lớn—bị chi phối mạnh mẽ bởi hạt nhân chất lỏng bên trong. Tỷ số bán kính lõi $R_{io}$ và tỷ số độ nhớt $\mu_{mo}$ có thể triệt tiêu hoặc đảo ngược hoàn toàn lực đẩy Marangoni vĩ mô.
  2. Bổ sung lý thuyết động lực học đường tiếp xúc ba pha (Dynamic Contact Line): Bằng việc tích hợp mô hình góc tiếp xúc động Kistler cải tiến bởi Muradoglu và Tasoglu (2005), luận án giải thích cơ chế hình thành vi màng bôi trơn và sự cân bằng động giữa lực mao dẫn nhiệt cục bộ với lực cản nhớt phân tán.
  3. Xác lập mô hình trạng thái cân bằng lơ lửng (Thermo-levitation state): Thiết lập mệnh đề lý thuyết chứng minh sự cân bằng bền vững giữa lực đẩy mao dẫn nhiệt hướng lên và trọng lực/lực nảy thủy động, tạo ra chiều cao ổn định $H_s$ của giọt đa lớp phía trên bề mặt không dính ướt.
       [Gradient nhiệt độ: ∇T] 
         [Hành vi chuyển động & Biến dạng vĩ mô]
         - Chuyển tiếp độ lõm trong kênh hẹp
         - Nảy & Ổn định lơ lửng trên mặt kị nước
         - Đảo chiều quỹ đạo trên mặt dính ướt

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 lý thuyết: Cơ học môi trường liên tục (Continuum Mechanics), Nhiệt động lực học mặt phân cách (Interfacial Thermodynamics) và Động học vi lưu đa pha (Multiphase Micro-hydrodynamics). Cách tiếp cận mới sử dụng hàm chỉ thị Heaviside trơn kết hợp phương pháp theo dấu biên Front Tracking (FT) cho phép theo dõi chính xác từng milimét biến dạng của cả mặt phân cách trong ($f_i$) và mặt phân cách ngoài ($f_o$).

Điều kiện biên điều kiện xác định (boundary conditions) được kiểm soát nghiêm ngặt: dòng chảy tầng không nén được ($Re < 100$), dao động nhiệt không làm thay đổi tỷ trọng và độ nhớt thể tích ($\text{D}\rho/\text{D}t = 0, \text{D}\mu/\text{D}t = 0$), hàm phụ thuộc sức căng bề mặt tuyến tính $\sigma = \sigma_0 - \sigma_T(T-T_0)$, và loại trừ bức xạ nhiệt.

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 duy lý (Positivism / Computational Realism) thông qua phương pháp mô phỏng số thủy động lực học tính toán (CFD) độ chính xác cao. Thiết kế mô phỏng đa tầng (multi-level design) bao gồm:

  • Cấp độ vĩ mô: Trường dòng chảy áp suất - vận tốc và trường nhiệt độ trên lưới Euler toàn cục.
  • Cấp độ vi mô: Sự biến dạng hình học, độ cong và lực căng bề mặt trên lưới Lagrange cục bộ di động của mặt phân cách.
  • Cấp độ nano/micro: Động lực học đường tiếp xúc ba pha tại tường rắn thông qua mô hình giải tích góc tiếp xúc động.

Kích thước miền tính toán được số hóa chính xác trên hệ tọa độ Descartes và tọa độ trụ 2D với lưới so le (Staggered Grid / MAC Grid).

Quy trình nghiên cứu rigorous

Quy trình mô phỏng sử dụng phương pháp theo dấu biên Front Tracking (Unverdi & Tryggvason, 1992) kết hợp giải thuật tách biến Fractional Step Method. Các điểm nút Lagrange trên biên phân cách được liên tục tái cấu trúc (thêm/xóa điểm) để duy trì mật độ tối ưu từ 2 đến 4 điểm trên mỗi ô lưới Euler cố định, đảm bảo tính liên tục của hàm delta Dirac $\delta(\mathbf{x} - \mathbf{x}_f)$ và triệt tiêu dao động giả.

       [Bước t = n: Vận tốc u^n, Áp suất p^n, Nhiệt độ T^n, Vị trí biên x_f^n]
             [Tái cấu trúc biên Lagrange: Thêm/Xóa điểm (2-4 điểm/ô)]
             [Tính phân bố hàm chỉ thị Heaviside I qua phương trình Poisson]
             [Cập nhật thuộc tính vật lý: ρ(I), μ(I), k(I), cp(I)]
             [Tính toán lực căng bề mặt Marangoni: δf_s = (σt)_2 - (σt)_1]
             [Giải phương trình Năng lượng tường minh: Tính T^(n+1)]
       [Giải phương trình Navier-Stokes qua Fractional Step: Dự báo u*]
          [Giải phương trình Poisson áp suất bằng phương pháp lặp SOR]
       [Cập nhật vận tốc u^(n+1) & Mô hình góc tiếp xúc động θ_D]
             [Di chuyển biên Lagrange x_f^(n+1) = x_f^n + Δt.v_f]

Độ tin cậy và tính hợp lệ của phương pháp số được chứng thực qua tam giác hóa phương pháp (Methodological Triangulation) với 5 chuẩn kiểm nghiệm quốc tế độc lập:

  • Kiểm chứng kênh dẫn thẳng: Trùng khớp hình thái hạt hai lớp của Borthakur và cộng sự (2018).
  • Kiểm chứng biến dạng va chạm: Tỷ lệ co dãn hình học $D_h$ tương thích hoàn toàn dữ liệu thực nghiệm của Qu và Hu (2012).
  • Kiểm chứng vận tốc mao dẫn nhiệt: Tiệm cận chính xác nghiệm giải tích của Young và cộng sự (1959) ở $Ma = 0{,}1; 1{,}0; 10{,}0$ và kết quả số của Kalichetty và cộng sự (2018).
  • Kiểm nghiệm trạng thái tĩnh: Khảo sát các góc $\theta_e = 60^\circ, 90^\circ, 120^\circ, 150^\circ$, sai số tương đối cực đại chỉ dừng lại ở mức $1{,}3%$.
  • Kiểm nghiệm động học dính ướt: Đường cong bán kính dính ướt $R_w(t)$ bám sát thực nghiệm giọt nước trên sáp paraffin ($\theta_e = 93^\circ$) của Šikalo và cộng sự (2005) ở điều kiện $Re = 36$, $We = 93$, $Eo = 0{,}285$.

Data và phân tích

Hệ thống phương trình rời rạc được giải trên nền tảng mã nguồn FORTRAN tự lập trình:

  • Thành phần đối lưu và khuếch tán được rời rạc hóa bằng sai phân hữu hạn (Finite Difference Method) với lược đồ sai phân trung tâm và sai phân tiến bậc 2.
  • Phương trình Poisson áp suất được giải lặp bằng phương pháp Thư giãn liên tiếp quá mức (Successive Over-Relaxation - SOR) với hệ số thư giãn tối ưu nhằm đảm bảo tính hội tụ nghiệm tuyệt đối.
  • Bước thời gian $\Delta t$ được kiểm soát nghiêm ngặt theo tiêu chuẩn ổn định Courant-Friedrichs-Lewy (CFL) và giới hạn độ nhớt khuếch tán: $\frac{\mu \Delta t}{\rho h^2} \le \frac{1}{4}$.

Mô hình góc tiếp xúc động tuân theo hàm thực nghiệm Hoffman: $$\theta_D = f_{Hoff}\left(Ca + f_{Hoff}^{-1}(\theta_e)\right)$$ với: $$f_{Hoff}(x) = \arccos\left{1 - 2\tanh\left[5{,}16\left(\frac{x}{1 + 1{,}31x^{0{,}99}}\right)^{0{,}706}\right]\right}$$

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

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

Luận án mang lại 4 phát hiện khoa học đột phá có giá trị thực nghiệm và mô phỏng cao:

  1. Phát hiện chuyển tiếp trạng thái lõm trong kênh thắt hẹp: Khi hạt đa lớp di chuyển qua vùng thắt hình sin, bề mặt phía sau của giọt ngoài xuất hiện hiện tượng lõm sâu. Mức độ lõm phụ thuộc phi tuyến vào số mao dẫn $Ca$ và tỷ số bán kính $R_{io}$. Luận án đã xây dựng biểu đồ trạng thái (Regime Diagram) chỉ rõ ranh giới chuyển tiếp từ trạng thái không lõm sang lõm đơn và lõm kép.
  2. Phát hiện hiện tượng nảy và ổn định lơ lửng (Rebound & Hovering): Khi va chạm vào bề mặt rắn không dính ướt dưới tác động của trường nhiệt độ thẳng đứng ($T_{hot}$ phía trên, $T_{cold}$ tại mặt đáy), lực mao dẫn nhiệt hướng lên hỗ trợ thắng lực cản nhớt, làm giọt đa lớp nảy ngược lên đạt chiều cao cực đại $H_m$ và dao động tắt dần để duy trì trạng thái lơ lửng ổn định tại khoảng cách $H_s$ so với bề mặt rắn.
  3. Phát hiện đảo chiều quỹ đạo di chuyển do hạt nhân bên trong: Trên bề mặt dính ướt có gradient nhiệt ngang, sự chênh lệch cường độ của hai cấu trúc xoáy đối lưu nhiệt bên trong giọt ngoài quyết định hướng chuyển động. Khi tăng tỷ số bán kính hạt trong $R_{io}$, hạt nhân bên trong cản trở dòng tuần hoàn nội tại, làm thay đổi phân bố áp suất và có thể đảo ngược hoàn toàn hướng di chuyển từ phía lạnh sang phía nóng hoặc ngược lại.
  4. Bằng chứng dữ liệu trích xuất trực tiếp: Dữ liệu mô phỏng chứng minh rằng khi $\theta_e < 90^\circ$, hướng di chuyển phụ thuộc vào tương quan lực xoáy, trong khi ở góc dính ướt lớn ($\theta_e > 90^\circ$), ứng suất nhiệt mao dẫn tại chân giọt chiếm ưu thế tuyệt đối.

Implications đa chiều

  • Về mặt lý thuyết: Mở rộng cơ học chất lưu đa pha từ hệ thống một mặt phân cách sang hệ thống đa mặt phân cách ghép cặp trường nhiệt, đặt nền móng lý thuyết mới cho nhiệt vi lưu (Thermal Microfluidics).
  • Về mặt phương pháp luận: Thuật toán Front Tracking cải tiến có khả năng mô hình hóa sự biến dạng cực hạn của nhũ tương kép mà không bị mất khối lượng hay khuếch tán số bề mặt như các phương pháp VOF hoặc Level-Set thông thường.
  • Về mặt ứng dụng thực tiễn: Mở đường cho việc thiết kế các kênh vi lưu thông minh trong phân tách tế bào, vi bọc thuốc chống ung thư tự động điều hướng bằng vi cảm biến nhiệt mà không làm phá vỡ màng sinh học.
  • Về mặt chính sách và công nghệ: Cung cấp cơ sở khoa học để các viện nghiên cứu và doanh nghiệp công nghệ cao tại Việt Nam làm chủ công nghệ thiết kế chip Lab-on-a-Chip và Organs-on-a-Chip thế hệ mới.

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học, luận án thẳng thắn chỉ ra các giới hạn nội tại:

  1. Giới hạn chiều không gian mô hình: Mô phỏng hiện tại được thực hiện chủ yếu trên không gian 2D và 2D đối xứng trục, chưa phản ánh đầy đủ các bất đối xứng 3D phức tạp khi giọt bị vặn xoắn trong không gian ba chiều.
  2. Giả định màng chất lỏng dư: Trong thực nghiệm thực tế luôn tồn tại một lượng nhỏ màng chất lỏng cực mỏng lưu lại trên bề mặt rắn khi giọt di chuyển, tuy nhiên mô hình số hiện tại chưa tái hiện hoàn toàn lớp màng vi mô này do giới hạn độ phân giải lưới.
  3. Bỏ qua biến thiên đặc tính nhiệt: Giả định khối lượng riêng $\rho$, độ nhớt $\mu$ và độ dẫn nhiệt $k$ là hằng số theo nhiệt độ (chỉ xét biến thiên của sức căng bề mặt $\sigma(T)$) có thể gây sai số nhỏ khi gradient nhiệt độ thực tế quá lớn.

Hướng phát triển nghiên cứu tương lai:

  • Nâng cấp chương trình tính toán lên không gian ba chiều (3D Front Tracking) trên hệ thống tính toán hiệu năng cao (HPC/GPU).
  • Tích hợp thêm các trường ngoại lực ghép cặp như điện trường (Electro-thermocapillarity) hoặc từ trường để điều khiển giọt đa chức năng.
  • Thực hiện các kiểm chứng thực nghiệm trực tiếp trên chip vi lưu chế tạo bằng công nghệ in 3D nano và quang khắc mềm (Soft-lithography).
  • Khảo sát nhũ tương đa nhân (multiple core droplets) với các phản ứng hóa sinh diễn ra đồng thời bên trong lõi.

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

Luận án tạo ra tác động sâu rộng trên nhiều bình diện:

  • Ảnh hưởng học thuật (Academic Impact): Đặt nền tảng nghiên cứu chuyên sâu về nhiệt mao dẫn hạt đa lớp tại Việt Nam; 04 công trình công bố trên các tạp chí quốc tế chuyên ngành Cơ học chất lỏng là minh chứng cho chất lượng học thuật đạt chuẩn thế giới, dự kiến thu hút trích dẫn cao trong cộng đồng vi lưu quốc tế.
  • Chuyển đổi công nghệ y dược: Cung cấp công cụ tính toán chính xác hỗ trợ các kỹ sư y sinh tối ưu hóa quy trình bao bọc và giải phóng hoạt chất thuốc hướng đích, rút ngắn thời gian và chi phí thử nghiệm tiền lâm sàng.
  • Phát triển công nghiệp bán dẫn và vật liệu: Hỗ trợ công nghệ in phun vi mạch chính xác cao, tự làm sạch bề mặt và làm mát vi mạch điện tử công suất lớn thông qua vi bọt mao dẫn nhiệt.

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

  • Nghiên cứu sinh và học viên cao học ngành Cơ học, Vật lý kỹ thuật: Tiếp cận phương pháp lập trình Front Tracking tiên tiến và mô hình toán học giải quyết bài toán đa pha ghép cặp nhiệt.
  • Các nhà nghiên cứu cấp cao (Senior Academics): Khai thác các phát hiện về cân bằng xoáy nhiệt Marangoni để phát triển các lý thuyết thủy động lực học phi tuyến.
  • Kỹ sư R&D trong công nghiệp Vi lưu và Y tế: Ứng dụng các biểu đồ trạng thái (Regime Diagrams) để thiết kế kích thước hình học kênh dẫn và chế độ nhiệt độ vận hành tối ưu cho thiết bị phân tích y sinh.
  • Các cơ quan hoạch định chiến lược khoa học công nghệ: Có thêm luận cứ thực tiễn để đầu tư phát triển phòng thí nghiệm trọng điểm về công nghệ vi lỏng và thiết bị y tế kỹ thuật cao.

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 Thuyết nhiệt mao dẫn Marangoni kinh điển sang hệ thống giọt đa lớp có tương tác thành rắn. Luận án đã chứng minh được rằng sự hiện diện của hạt nhân bên trong làm tái cấu trúc hoàn toàn hai vòng xoáy nhiệt đối lưu nội tại, tạo ra cơ chế đảo chiều chuyển động vĩ mô và thiết lập trạng thái cân bằng lơ lửng nhiệt mao dẫn (Thermo-levitation) mới trên bề mặt kị nước.

2. Tính đổi mới trong phương pháp nghiên cứu khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?

Trả lời: So với nghiên cứu của Sui (2014) sử dụng điều kiện biên trượt Navier chưa kiểm chứng ở góc tiếp xúc lớn, và nghiên cứu của Borthakur và cộng sự (2018) chỉ khảo sát giọt hai lớp trong kênh thẳng đẳng nhiệt, luận án đã:

  • Kết hợp thành công phương pháp theo dấu biên (Front Tracking) với mô hình góc tiếp xúc động Muradoglu - Tasoglu.
  • Giải quyết đồng thời bài toán biên dạng kênh thắt phi tuyến và tương tác tiếp xúc dính ướt dưới tác động của trường nhiệt độ không đồng nhất với sai số kiểm chứng dưới $1{,}3%$.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng số liệu chứng minh?

Trả lời: Phát hiện bất ngờ nhất là kích thước hạt bên trong có thể làm đảo ngược hoàn toàn hướng di chuyển của hạt bên ngoài trên bề mặt dính ướt. Khi tăng tỷ số bán kính $R_{io}$ vượt qua ngưỡng tới hạn, dòng xoáy Marangoni phụ bị biến dạng ép dòng xoáy chính đổi hướng, khiến giọt dịch chuyển về phía ngược lại với quy luật thông thường của giọt đơn lớp ở cùng góc dính ướt $\theta_e$.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không?

Trả lời: Hoàn toàn có. Luận án trình bày chi tiết từng bước toán học: từ hệ phương trình vi phân Navier-Stokes, phương trình năng lượng, cách tính hàm chỉ thị $I$, công thức rời rạc hóa lực căng bề mặt trên lưới so le MAC, điều kiện biên vận tốc ảo và nhiệt độ ảo, đến sơ đồ thuật toán tổng thể giải phương trình Poisson áp suất bằng phương pháp lặp SOR. Mã nguồn FORTRAN được cấu trúc hóa logic cho phép tái lập độc lập.

5. Lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Trả lời: Lộ trình 10 năm hướng tới:

  • Giai đoạn 1 (1-3 năm): Mở rộng hoàn chỉnh mã nguồn mô phỏng sang 3D song song hóa đa luồng trên GPU.
  • Giai đoạn 2 (4-6 năm): Chế tạo thực nghiệm chip vi lưu kiểm chứng hiện tượng nhiệt mao dẫn giọt đa lớp.
  • Giai đoạn 3 (7-10 năm): Ứng dụng thực tế vào sản xuất hệ thống thiết bị y tế Lab-on-a-Chip chẩn đoán sớm bệnh lý phân tử và dẫn truyền thuốc thông minh tại Việt Nam.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Đình Hòe đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra với các kết luận tổng kết:

  1. Xây dựng thành công mô hình tính toán số hoàn chỉnh: Thiết lập chương trình mô phỏng FORTRAN chính xác cao dựa trên phương pháp Front Tracking kết hợp mô hình góc tiếp xúc động Kistler - Muradoglu - Tasoglu cho bài toán giọt đa lớp chịu tác động của mao dẫn nhiệt.
  2. Làm sáng tỏ cơ chế thủy động học trong kênh dẫn thắt: Xác lập quy luật biến dạng và phân vùng biểu đồ trạng thái độ lõm bề mặt sau của giọt đa lớp qua phần thắt hình sin trong điều kiện đẳng nhiệt.
  3. Phát hiện và định lượng hiện tượng nảy lơ lửng: Chứng minh bằng số liệu hiện tượng hạt đa lớp nảy lên và đạt trạng thái cân bằng lơ lửng ổn định trên bề mặt không dính ướt dưới tác động của gradient nhiệt độ đứng.
  4. Khám phá quy luật đảo chiều quỹ đạo nhiệt mao dẫn: Chỉ ra vai trò quyết định của cấu trúc hạt nhân bên trong trong việc làm thay đổi phân bố xoáy nội tại và điều khiển hướng di chuyển của giọt trên bề mặt dính ướt.
  5. Giá trị công bố quốc tế vượt trội: Toàn bộ các kết quả cốt lõi đã được phản biện độc lập và công bố trên 04 bài báo tạp chí quốc tế ISI/Scopus uy tín, khẳng định vị thế tiên phong của nghiên cứu cơ học chất lưu vi mô tại Việt Nam trên trường quốc tế.