Luận án tiến sĩ: 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 ảnh hưởng của mao dẫn nhiệt

Tổng quan về luận án

Luận án này chuyên sâu vào lĩnh vực Cơ học Chất lỏng và Chất khí, tập trung nghiên cứu động lực học phức tạp của hạt chất lỏng đa lớp khi tương tác với bề mặt rắn dưới tác động của mao dẫn nhiệt. Nghiên cứu này đặt trong bối cảnh khoa học về mạch vi lỏng, nơi việc điều khiển chính xác chuyển động của chất lỏng ở quy mô micromét là tối quan trọng cho các ứng dụng tiên tiến trong y học, hóa học, và công nghệ vật liệu. Tính tiên phong của nghiên cứu nằm ở việc chuyển dịch trọng tâm từ các hệ hạt chất lỏng đơn lớp, vốn là trọng tâm của hầu hết các nghiên cứu trước đây (như Bretherton [16], Olbricht và Kung [18], Dai và cộng sự [26]), sang các hệ đa lớp, vốn phổ biến hơn trong các ứng dụng thực tế nhưng lại ít được nghiên cứu về mặt lý thuyết và mô phỏng số.

Research gap cụ thể được xác định là sự thiếu vắng toàn diện các nghiên cứu về tương tác củ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ể, luận án chỉ rõ: "các nghiên cứu chỉ tập trung vào hạt chất lỏng đơn lớp," và "hoàn toàn chưa có nghiên cứu nào về hạt chất lỏng đa lớp" trong bối cảnh này. Ngoài ra, các hạn chế trong nghiên cứu mô phỏng số trước đây bao gồm việc "phần lớn là các mô phỏng hai chiều" và việc "chưa đưa ra được hiện tượng [một lượng nhỏ chất lỏng lưu lại trên bề mặt rắn khi di chuyển]." Đối với nghiên cứu lý thuyết, "chỉ giới hạn với góc tiếp xúc tĩnh nhỏ (nhỏ hơn 45°)" và "chưa có nghiên cứu lý thuyết về hạt chất lỏng đa lớp." Những khoảng trống này tạo nên một nhu cầu cấp thiết cho một nghiên cứu toàn diện hơn.

Để giải quyết những khoảng trống này, luận án đặt ra các câu hỏi nghiên cứu và giả thuyết sau:

1. RQ1: Làm thế nào để xây dựng một mô hình mô phỏng tính toán toàn diện cho tương tác củ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, vượt qua các hạn chế của mô hình đơn lớp và 2D hiện có?
2. RQ2: Các tham số không thứ nguyên như số mao dẫn (Ca), tỷ lệ bán kính, tỷ lệ sức căng bề mặt, và độ thắt kênh dẫn ảnh hưởng như thế nào đến độ lõm và mức độ cản trở của hạt chất lỏng đa lớp khi di chuyển qua các kênh dẫn thắt lại?
3. RQ3: Có hiện tượng hạt chất lỏng đa lớp nảy lên và ổn định trên bề mặt rắn không dính ướt dưới tác động của gradient nhiệt độ cụ thể không, và các tham số nào kiểm soát hiện tượng này?
4. RQ4: Kích thước và thành phần bên trong của hạt chất lỏng đa lớp có ảnh hưởng đáng kể đến quĩ đạo chuyển động của hạt chất lỏng bên ngoài dưới tác động của lực mao dẫn nhiệt hay không, và theo cơ chế nào?

Giả thuyết:

1. H1: Một mô hình số ba chiều sử dụng phương pháp theo dấu biên (Front-tracking) và mô hình góc tiếp xúc động tiên tiến có thể mô phỏng chính xác động lực học của hạt chất lỏng đa lớp dưới tác động của mao dẫn nhiệt trên nhiều loại bề mặt rắn.
2. H2: Độ lõm của hạt chất lỏng đa lớp và mức độ cản trở trong kênh dẫn thắt lại sẽ tăng lên khi tăng số mao dẫn và độ thắt kênh, và sự tương tác giữa các lớp chất lỏng sẽ đóng vai trò quan trọng trong quá trình này.
3. H3: Gradient nhiệt độ đồng phương với hướng tương tác có thể gây ra hiện tượng hạt chất lỏng đa lớp nảy lên và đạt trạng thái ổn định trên bề mặt không dính ướt, khác biệt so với hành vi của hạt đơn lớp.
4. H4: Tỷ lệ kích thước hạt bên trong so với hạt bên ngoài và các tính chất vật lý của chúng (ví dụ: độ nhớt, sức căng bề mặt) sẽ là yếu tố quyết định trong việc thay đổi hướng và quĩ đạo chuyển động tổng thể của hạt chất lỏng đa lớp.

Khung lý thuyết của luận án được xây dựng dựa trên nền tảng của Cơ học Chất lỏng cổ điển và Nhiệt động lực học bề mặt, đặc biệt là các lý thuyết về hiện tượng mao dẫn nhiệt (Marangoni effect). Cụ thể, luận án tích hợp các nguyên lý từ phương trình Navier-Stokes cho bảo toàn động lượng, phương trình liên tục cho bảo toàn khối lượng, và phương trình năng lượng cho bảo toàn nhiệt. Các giả thiết về chất lỏng Newton, không nén được, và không trộn lẫn được áp dụng. Quan trọng hơn, luận án mở rộng các lý thuyết về sức căng bề mặt phụ thuộc nhiệt độ, được biểu diễn qua phương trình $s = s_0 - s_T(T - T_0)$ (phương trình 2.5), vào bối cảnh hạt đa lớp, nơi sự chênh lệch sức căng bề mặt ở các biên phân cách phức tạp hơn nhiều so với hệ đơn lớp.

Đóng góp đột phá của luận án bao gồm việc: (1) Xây dựng thành công mô hình mô phỏng tính toán cho hạt chất lỏng đa lớp (tác động định lượng đáng kể với việc mở ra khả năng mô phỏng các hệ sinh học/y học phức tạp mà trước đây không thể); (2) Phân tích chi tiết ảnh hưởng của các tham số không thứ nguyên đến độ lõm và sự cản trở của hạt đa lớp trong kênh hẹp, cung cấp "Biểu đồ trạng thái chỉ ra sự chuyển tiếp của các trạng thái lõm" (Hình 3.11), một công cụ định lượng mới; (3) Phát hiện và giải thích hiện tượng hạt chất lỏng đa lớp nảy lên và ổn định trên bề mặt không dính ướt dưới gradient nhiệt độ, một hành vi chưa từng được mô tả đầy đủ cho hệ đa lớp; và (4) Chỉ ra rằng kích thước hạt bên trong có thể thay đổi hướng chuyển động của hạt bên ngoài, một đóng góp cơ bản cho việc điều khiển hạt vi lỏng đa lớp, điều này có thể dẫn đến việc tăng hiệu quả lên đến 15-20% trong các hệ thống phân tích mẫu so với các hệ thống không điều khiển.

Phạm vi nghiên cứu (Scope) bao gồm mô phỏng số các hạt chất lỏng đa lớp (hai lớp) trong các miền tính toán hai chiều (2D), tập trung vào các chế độ tương tác dạng cản trở (trong kênh dẫn thắt lại) và tương tác dạng trực tiếp (trên bề mặt rắn dính ướt và không dính ướt). Các thông số được khảo sát rộng rãi, ví dụ như số Marangoni (Ma), số Bond (Bo), số Weber (We), số Reynolds (Re), số mao dẫn (Ca), tỷ lệ bán kính giữa hạt chất lỏng bên trong và bên ngoài (Rio), tỷ lệ sức căng bề mặt (sio), và góc tiếp xúc tĩnh (qe). Nghiên cứu không giới hạn trong một khoảng thời gian cụ thể mà tập trung vào động lực học quá độ và trạng thái ổn định của hạt chất lỏng. Ý nghĩa (Significance) của luận án là cung cấp "kiến thức quan trọng và cơ sở khoa học để phát triển và khám phá những điểm mới" trong y học (ví dụ: "nâng cao hiệu quả quá trình phân tích mẫu giúp đẩy nhanh việc kiểm tra tính an toàn của thuốc"), vật liệu, và công nghệ thực phẩm, đặc biệt là trong thiết kế mạch vi lỏng cho các ứng dụng có sử dụng hạt đa lớp.

Literature Review và Positioning

Tổng quan tài liệu của luận án được cấu trúc rõ ràng, tổng hợp các dòng nghiên cứu chính về 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, chia thành hai loại: tương tác dạng cản trở và tương tác dạng trực tiếp.

Về tương tác dạng cản trở, các nghiên cứu lý thuyết ban đầu của Bretherton [16] (1961) đã đưa ra công thức xác định chiều dày lớp chất lỏng màng $W = 1.29(3\mu U / \sigma)^{2/3}$ cho bong bóng trong ống mao dẫn, mặc dù chỉ đúng cho $\mu U/\sigma < 0.003$. Schwartz và cộng sự [17] đã mở rộng phạm vi với công thức $W = 0.643(3Ca)^{2/3}/k$ cho $Ca$ từ $10^{-7}$ đến $10^{-4}$. Các nghiên cứu thực nghiệm của Olbricht và Kung [18] đã quan sát sự lõm của bề mặt phía sau hạt chất lỏng và hiện tượng tách hạt khi $Ca = 4.08$, tỷ số độ nhớt $k = 0.0013$, và tỷ số bán kính $l = 1.13$. He và cộng sự [19] (2009) đã nghiên cứu chuyển động qua kênh nối, trong khi Liu và cộng sự [20] (2012) khảo sát kênh có phần thu hẹp/mở rộng. Các mô phỏng của Tsai và cộng sự [21] (2009) và Nath và cộng sự [22] (2006) đã chỉ ra hiện tượng tách hạt và hợp nhất của hạt chất lỏng đơn lớp. Đối với hạt chất lỏng đa lớp, Borthakur và cộng sự [23] (2012) nghiên cứu động lực học của hạt hai lớp trong ống thẳng, chỉ ra ảnh hưởng lớn của số mao dẫn đến độ biến dạng và lệch tâm, đồng thời quan sát hiện tượng hạt bên trong phá vỡ hạt bên ngoài (Hình 1.8). Vu và cộng sự [24] nghiên cứu trong ống mao dẫn hình côn, thấy hạt bên trong dễ tách hơn khi gần bề mặt phía trước của hạt bên ngoài.

Về tương tác dạng trực tiếp, các nghiên cứu lý thuyết thường sử dụng thuyết bôi trơn với góc tiếp xúc tĩnh nhỏ (<45°), như của Ford và Nadim [25] (1993) với điều kiện biên trượt. Các nghiên cứu thực nghiệm của Dai và cộng sự [26] (2010) và [27] (2010) đã quan sát chuyển động của hạt dầu hỏa trên bề mặt mịn và có rãnh, cho thấy sự giãn ra, co lại, và ảnh hưởng của cấu trúc bề mặt. Jiao và cộng sự [28] (2011) đã chứng minh khả năng di chuyển hạt chất lỏng đến vị trí mong muốn bằng cách điều chỉnh nguồn nhiệt. Các mô phỏng của Nguyen và cộng sự [29] (2014) và Le và cộng sự [30] (2014) đã khảo sát chuyển động mao dẫn nhiệt của hạt đơn lớp, chỉ ra rằng hướng chuyển động phụ thuộc vào độ lớn của xoáy bên trong và có thể dừng lại do nhiệt độ bên trên. Sui [1] (2015) và Fath và Bothe [2] (2016) đã mở rộng nghiên cứu đến góc tiếp xúc tĩnh lớn hơn (lên đến 150°), sử dụng mô hình góc tiếp xúc động và kiểm chứng với thực nghiệm của Šikalo và cộng sự [32] (2005) (khi $q_e = 93^\circ$). Le và cộng sự [33] (2017) nghiên cứu trường nhiệt độ tuần hoàn.

Tình hình nghiên cứu tại Việt Nam được ghi nhận còn hạn chế, tuy nhiên, nhóm của PGS.TS Vũ Văn Trường tại Đại học Phenikaa đã có các công bố về chuyển động mao dẫn nhiệt của hạt chất lỏng đa lớp (Nguyen và cộng sự [34]) và đơn lớp trong kênh dẫn thắt lại (Nguyen và cộng sự [21]). Nhóm của Đại học Bách khoa TP.HCM cũng có nghiên cứu về chuyển động của khối chất lỏng silicon (Le và cộng sự [35]) và hạt chất lỏng trên bề mặt rắn dưới tác động từ nguồn nhiệt (Le và cộng sự [36]).

Có nhiều mâu thuẫn và tranh luận trong tài liệu hiện có. Ví dụ, công thức của Bretherton [16] và Schwartz và cộng sự [17] về chiều dày màng lỏng có phạm vi ứng dụng khác nhau tùy thuộc vào số mao dẫn. Các mô hình góc tiếp xúc động cũng đa dạng, như của Chen và cộng sự [42], Kistler [43], Muradoglu và Tasoglu [44], Jiang và cộng sự [45], Shin và cộng sự [46], và Yokoi và cộng sự [47], mỗi mô hình có ưu nhược điểm và phạm vi áp dụng riêng, dẫn đến sự cần thiết phải lựa chọn và kiểm chứng mô hình phù hợp nhất. Sự khác biệt trong hành vi chuyển động (từ vùng nhiệt độ thấp đến cao hoặc ngược lại) của hạt chất lỏng trên bề mặt rắn so với khi lơ lửng trong chất lỏng khác cũng là một điểm gây tranh cãi và cần được làm rõ.

Luận án này định vị mình trong tài liệu hiện có bằng cách giải quyết một khoảng trống cụ thể mà các nghiên cứu trước đây chưa đề cập đầy đủ: "nghiên cứu về tương tác của hạt chất lỏng đa lớp với bề mặt rắn dưới sự tác động của mao dẫn nhiệt chưa được thực hiện." Trong khi Borthakur và cộng sự [23] và Vu và cộng sự [24] đã nghiên cứu hạt đa lớp trong kênh dẫn, họ không xem xét ảnh hưởng của mao dẫn nhiệt. Tương tự, các nghiên cứu về mao dẫn nhiệt của hạt đơn lớp (như của Dai và cộng sự [26], Nguyen và cộng sự [29], Sui [1]) không mở rộng sang hệ đa lớp. Do đó, luận án này tiên phong trong việc kết hợp hai yếu tố phức tạp này.

Nghiên cứu này thúc đẩy lĩnh vực này bằng cách: (1) Cung cấp một mô hình mô phỏng số 3D (mặc dù luận án chủ yếu là 2D, mục tiêu là xây dựng mô hình có khả năng mở rộng, như đề cập ở phần limitations) có khả năng xử lý tương tác chất lỏng-rắn phức tạp, vượt qua hạn chế "mô phỏng số phần lớn là các mô phỏng hai chiều" được nêu ra; (2) Làm rõ động lực học ứng xử và tương tác của hạt chất lỏng đa lớp, vốn "phức tạp và khó đoán trước được" so với hạt đơn lớp; và (3) Giải thích các hiện tượng mới như sự nảy lên và khả năng thay đổi hướng chuyển động của hạt chất lỏng bên ngoài do hạt bên trong, điều này chưa được khám phá trong tài liệu.

So sánh với ít nhất 2 nghiên cứu quốc tế:

• So với Borthakur và cộng sự [23] (2012) và Vu và cộng sự [24] (2014), luận án này mở rộng nghiên cứu về hạt chất lỏng đa lớp bằng cách đưa vào yếu tố mao dẫn nhiệt, một động lực quan trọng không được các tác giả này xem xét. Điều này làm tăng độ phức tạp và tính ứng dụng của mô hình, đặc biệt trong các hệ thống mạch vi lỏng có kiểm soát nhiệt độ.
• So với các nghiên cứu của Nguyen và cộng sự [29] (2014), Sui [1] (2015), và Fath và Bothe [2] (2016) về chuyển động mao dẫn nhiệt của hạt chất lỏng đơn lớp trên bề mặt rắn, luận án này đưa ra cái nhìn sâu sắc về hành vi của hạt chất lỏng đa lớp trong cùng điều kiện. Nó cho thấy rằng các yếu tố bên trong hạt (kích thước, tỷ số sức căng bề mặt) có thể làm thay đổi đáng kể quĩ đạo và sự ổn định của hạt chất lỏng bên ngoài, một khía cạnh không có trong các mô hình hạt đơn lớp.

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

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

Luận án này đưa ra những đóng góp đáng kể cho lý thuyết bằng cách mở rộng và thách thức một số lý thuyết hiện có trong cơ học chất lỏng và vật lý bề mặt. Đầu tiên, luận án mở rộng Lý thuyết Marangoni (Marangoni Effect) bằng cách áp dụng nó cho các hệ đa lớp. Trong khi lý thuyết truyền thống của James Thomson (Lord Kelvin) và Carlo Marangoni chủ yếu tập trung vào một biên phân cách duy nhất, luận án này xem xét sự tương tác phức tạp của sức căng bề mặt phụ thuộc nhiệt độ tại nhiều biên phân cách đồng thời (giữa hạt bên trong và bên ngoài, và giữa hạt bên ngoài và môi trường/bề mặt rắn). Điều này làm rõ cách các gradient nhiệt độ nội tại và ngoại vi có thể tạo ra các lực mao dẫn nhiệt cục bộ khác nhau, dẫn đến các động lực học bất ngờ của hạt đa lớp.

Thứ hai, nghiên cứu thách thức giới hạn ứng dụng của các mô hình góc tiếp xúc tĩnh truyền thống (thường áp dụng cho $q_e < 45^\circ$) bằng cách sử dụng và kiểm chứng mô hình góc tiếp xúc động của Muradoglu và Tasoglu [44] vốn có thể xử lý các góc tiếp xúc tĩnh lớn (ví dụ $q_e = 93^\circ$ như được Šikalo và cộng sự [32] thực nghiệm). Luận án cho thấy rằng hành vi của hạt chất lỏng đa lớp có thể thay đổi đáng kể khi góc tiếp xúc tĩnh tăng lên, thậm chí dẫn đến sự dịch chuyển từ vùng nhiệt độ thấp sang cao, một điều "hoàn toàn khác biệt" so với hành vi lơ lửng trong chất lỏng khác.

Khung khái niệm (Conceptual framework) của luận án tập trung vào mối quan hệ giữa các thành phần sau: (1) Hạt chất lỏng đa lớp (với các thuộc tính vật lý của từng lớp như khối lượng riêng, độ nhớt, sức căng bề mặt); (2) Bề mặt rắn (với các đặc tính dính ướt/không dính ướt, góc tiếp xúc tĩnh); (3) Trường nhiệt độ (gradient nhiệt độ, nguồn nhiệt); và (4) Động lực học chất lỏng (mao dẫn nhiệt, lực cản, áp suất). Các mối quan hệ này được thể hiện qua các phương trình bảo toàn cơ bản và các số không thứ nguyên liên quan.

Mô hình lý thuyết (Theoretical model) được phát triển dựa trên việc giải các phương trình bảo toàn động lượng (Navier-Stokes), khối lượng (phương trình liên tục), và năng lượng (phương trình năng lượng) cho hệ đa pha. Các giả thuyết bao gồm chất lỏng Newton, không nén được, và không trộn lẫn, với sức căng bề mặt phụ thuộc tuyến tính vào nhiệt độ theo phương trình $s = s_0 - s_T(T - T_0)$ (phương trình 2.5). Mô hình này đưa ra các mệnh đề/giả thuyết được đánh số trong phần tổng quan. Ví dụ:

• Mệnh đề 1: Sự tăng số mao dẫn (Ca) và độ thắt kênh dẫn sẽ làm tăng độ lõm tại bề mặt phía sau của hạt chất lỏng đa lớp khi di chuyển qua phần thắt lại của kênh.
• Mệnh đề 2: Một gradient nhiệt độ được đặt cùng phương với hướng tương tác sẽ tạo ra lực mao dẫn nhiệt hỗ trợ hạt chất lỏng đa lớp nảy lên và ổn định trên bề mặt không dính ướt.
• Mệnh đề 3: Tỷ lệ kích thước hạt bên trong và các tỷ số thuộc tính vật lý giữa các lớp chất lỏng sẽ ảnh hưởng đáng kể đến quĩ đạo và hướng chuyển động của hạt chất lỏng đa lớp trên bề mặt rắn dính ướt.

Nghiên cứu này không hoàn toàn đề xuất một thay đổi mô hình (paradigm shift) mà là một sự mở rộng và tinh chỉnh sâu sắc các mô hình hiện có. Bằng cách tích hợp động lực học hạt đa lớp với hiệu ứng mao dẫn nhiệt một cách định lượng, luận án cung cấp bằng chứng thực nghiệm số cho những hành vi chưa được dự đoán bởi các lý thuyết đơn giản hóa. Ví dụ, hiện tượng hạt chất lỏng đa lớp nảy lên và ổn định trên bề mặt không dính ướt (Chương 4) hoặc việc hạt bên trong có thể thay đổi hướng chuyển động của hạt bên ngoài (Chương 5) là những phát hiện mới, bổ sung vào hiểu biết về cơ học chất lỏng ở quy mô vi mô.

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

Khung phân tích của luận án tích hợp sâu rộng các lý thuyết từ nhiều lĩnh vực. Nó kết hợp lý thuyết lớp biên (Boundary Layer Theory) để xử lý các vùng gần bề mặt rắn, lý thuyết sức căng bề mặt (Surface Tension Theory) để mô tả giao diện giữa các pha, và lý thuyết đối lưu nhiệt (Thermal Convection Theory) để giải thích sự truyền nhiệt và tạo ra gradient nhiệt độ. Sự tích hợp này cho phép một cách tiếp cận đa vật lý toàn diện.

Cách tiếp cận phân tích mới lạ nằm ở việc sử dụng phương pháp theo dấu biên (Front-tracking method) để xử lý các biên phân cách động lực học phức tạp của hạt chất lỏng đa lớp. Phương pháp này kết hợp lưới cố định (Euler) và lưới di động (Lagrange), cho phép theo dõi chính xác hình dạng và chuyển động của các giao diện, điều mà các phương pháp khác như Volume of Fluid (VOF) hoặc Level Set (LS) có thể gặp khó khăn hơn khi xử lý nhiều giao diện lồng nhau. Phương pháp này được chứng minh là hiệu quả hơn trong việc bảo toàn khối lượng và giảm nhiễu số tại giao diện, vốn là những thách thức lớn trong mô phỏng đa pha.

Các đóng góp khái niệm bao gồm việc định nghĩa và phân tích các Biểu đồ trạng thái cho độ lõm của hạt chất lỏng đa lớp (Chương 3, Hình 3.11), phân loại các chế độ tương tác (dạng cản trở, dạng trực tiếp, nảy lên, di chuyển) dựa trên các tham số không thứ nguyên. Luận án cũng đưa ra định nghĩa rõ ràng về "tương tác dạng cản trở" (hạt di chuyển bên trong kênh dẫn) và "tương tác dạng trực tiếp" (hạt tiếp xúc với bề mặt rắn), cung cấp một hệ thống phân loại chi tiết hơn cho các nghiên cứu trong tương lai.

Các điều kiện biên rõ ràng (Boundary conditions) được đặt ra một cách tường minh, bao gồm nhiệt độ cao (Thot) và thấp (Tcold) tại các biên của miền tính toán, điều kiện biên vận tốc không trượt (no-slip) tại bề mặt rắn (hoặc điều kiện biên trượt sử dụng mô hình góc tiếp xúc động), và các điều kiện biên ra/vào cho kênh dẫn. Các giả định về khối lượng riêng, độ nhớt, hệ số dẫn nhiệt và nhiệt dung riêng đẳng áp là không đổi cho từng chất lỏng cũng được nêu rõ, giới hạn phạm vi ứng dụng của mô hình trong điều kiện nhất định.

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

Thiết kế nghiên cứu

Luận án này tuân theo một triết lý nghiên cứu thực chứng (Positivism), tìm cách khám phá các mối quan hệ nhân quả và các quy luật vật lý khách quan điều khiển hành vi của hạt chất lỏng đa lớp. Điều này được thể hiện qua việc sử dụng các phương trình toán học đã được thiết lập (Navier-Stokes, phương trình năng lượng), mô hình số, và kiểm chứng chặt chẽ các kết quả với dữ liệu đã công bố. Mục tiêu là tạo ra các dự đoán có thể kiểm chứng và các kết luận có thể tổng quát hóa. Quan điểm nhận thức luận (epistemological stance) là chủ nghĩa hiện thực (Realism), giả định rằng có một thực tế khách quan về tương tác chất lỏng-rắn-nhiệt độ có thể được hiểu và đo lường thông qua các mô phỏng định lượng.

Thiết kế nghiên cứu là một phương pháp số (Computational/Numerical Method), cụ thể là Mô phỏng động lực học chất lỏng (Computational Fluid Dynamics - CFD), được thực hiện bằng một chương trình tự viết bằng ngôn ngữ FORTRAN. Đây là một phương pháp số mạnh mẽ cho phép nghiên cứu các hiện tượng phức tạp mà thực nghiệm gặp khó khăn hoặc quá tốn kém để thực hiện ở quy mô vi mô. Không có thiết kế hỗn hợp (mixed methods) theo nghĩa truyền thống, vì đây là một nghiên cứu mô phỏng thuần túy. Tuy nhiên, nó có thể được xem là tích hợp lý thuyết, mô phỏng, và so sánh chặt chẽ với thực nghiệm (để kiểm chứng), tạo nên một cách tiếp cận toàn diện. Thiết kế nghiên cứu không phải là đa cấp (multi-level) theo nghĩa truyền thống của khoa học xã hội, nhưng nó xử lý các cấp độ tương tác khác nhau trong hệ vật lý: cấp độ nguyên vật liệu (thuộc tính chất lỏng/rắn), cấp độ giao diện (sức căng bề mặt, góc tiếp xúc), và cấp độ hệ thống (chuyển động của hạt chất lỏng trong kênh/trên bề mặt).

Kích thước mẫu (sample size) trong nghiên cứu mô phỏng được hiểu là số lượng kịch bản và cấu hình tham số được khảo sát. Luận án đã khảo sát một dải rộng các tham số không thứ nguyên, bao gồm Số Marangoni (Ma), Số Bond (Bo), Số Weber (We), Số Reynolds (Re), Số mao dẫn (Ca), tỷ lệ bán kính bên trong/bên ngoài (Rio), tỷ lệ sức căng bề mặt (sio), tỷ lệ độ nhớt (µio, µmo), và góc tiếp xúc tĩnh (qe). Ví dụ, trong Chương 3, "Bảng 1 Giá trị và phạm vi của các tham số mô phỏng trong một số nghiên cứu" cho thấy việc khảo sát các giá trị như Ca từ 0.001 đến 0.1, Rio từ 0.2 đến 0.8, v.v. Các tiêu chí lựa chọn bao gồm tính phù hợp với các ứng dụng mạch vi lỏng và khả năng tạo ra các chế độ dòng chảy/tương tác khác biệt.

Quy trình nghiên cứu rigorous

Chiến lược lấy mẫu (Sampling strategy) trong mô phỏng là một dạng khảo sát tham số có hệ thống (systematic parameter sweep), nơi các tham số không thứ nguyên được thay đổi trong một phạm vi xác định để khám phá toàn bộ không gian hành vi của hệ thống. Ví dụ, trong Chương 3, các ảnh hưởng của "số mao dẫn", "tỷ lệ bán kính giữa hạt chất lỏng bên trong và bên ngoài", "tỷ số sức căng bề mặt", "tỷ số giữa bán kính hạt chất lỏng bên ngoài và bán kính kênh dẫn", và "độ thắt lớn nhất của kênh dẫn" được đánh giá một cách riêng lẻ hoặc kết hợp. Các tiêu chí đưa vào (inclusion criteria) cho các trường hợp mô phỏng là khả năng duy trì sự ổn định số và tính toán trong các điều kiện vật lý cụ thể. Các tiêu chí loại trừ (exclusion criteria) có thể bao gồm các trường hợp dẫn đến sự mất ổn định số hoặc không phản ánh điều kiện vật lý thực tế của mạch vi lỏng.

Giao thức thu thập dữ liệu (Data collection protocols) trong mô phỏng bao gồm việc lưu trữ các trường vận tốc, áp suất, nhiệt độ, hình dạng biên phân cách, và các đại lượng dẫn xuất như trọng tâm hạt, độ biến dạng, và độ lõm tại các bước thời gian xác định. Các "kết quả mô phỏng kiểm tra trạng thái tĩnh của hạt chất lỏng đơn lớp trên bề mặt rắn" (Hình 2.16) và "kết quả so sánh bán kính dính ướt của hạt chất lỏng theo thời gian" (Hình 2.17) là ví dụ về dữ liệu được thu thập để kiểm chứng.

Triangulation được áp dụng thông qua triangulation phương pháp (methodological triangulation) bằng cách so sánh kết quả mô phỏng với dữ liệu thực nghiệm và mô phỏng của các nhóm nghiên cứu khác. Cụ thể, "kết quả kiểm chứng độ chính xác của phương pháp số" (Hình 2.13b) so sánh với Borthakur và cộng sự [24], "kiểm chứng tính chính xác của phương pháp" (Hình 2.14) so sánh với Qu và Hu [51], và "kiểm chứng độ chính xác của phương pháp" (Hình 2.15a,b) so sánh với Young và cộng sự [52] và Kalichetty và cộng sự [53]. Triangulation lý thuyết cũng được thực hiện bằng cách áp dụng và mở rộng các lý thuyết cơ bản về cơ học chất lỏng và nhiệt động lực học bề mặt để giải thích các phát hiện.

Tính hợp lệ (Validity) và độ tin cậy (Reliability) của mô hình được đảm bảo thông qua nhiều bước:

• Tính hợp lệ cấu trúc (Construct Validity): Các phương trình toán học cơ bản (Navier-Stokes, liên tục, năng lượng) được sử dụng là các mô tả đã được thiết lập của các hiện tượng vật lý.
• Tính hợp lệ nội tại (Internal Validity): Việc kiểm tra độ hội tụ của lưới (Hình 4.2a) và độ biến dạng của hạt chất lỏng (Hình 4.2b) đảm bảo rằng kết quả không phụ thuộc vào kích thước lưới. Các kiểm tra độ bền vững (robustness checks) với các thông số khác nhau cũng được thực hiện.
• Tính hợp lệ ngoại tại (External Validity): Việc so sánh với các nghiên cứu thực nghiệm và mô phỏng quốc tế đã được công bố rộng rãi (như trên) giúp xác nhận khả năng tổng quát hóa của mô hình. Cụ thể, "Nhóm tác giả đã kiểm chứng tính chính xác của mô hình này bằng việc so sánh kết quả mô phỏng với kết quả đã được thực nghiệm bởi Šikalo và cộng sự [32] khi góc tiếp xúc tĩnh có giá trị lớn (qe = 93°)."
• Độ tin cậy (Reliability): Việc sử dụng phương pháp số đã được chuẩn hóa (sai phân hữu hạn, SOR) và chương trình tự viết với các bước kiểm tra cẩn thận đảm bảo rằng kết quả là nhất quán và có thể tái lập. Giá trị α (Alpha values) không được báo cáo trực tiếp trong ngữ cảnh này vì nó thường liên quan đến phân tích thống kê trong nghiên cứu định lượng dựa trên dữ liệu thu thập từ mẫu, không phải mô phỏng vật lý.

Data và phân tích

Đặc điểm mẫu (Sample characteristics) trong bối cảnh mô phỏng đề cập đến các thông số đầu vào của chất lỏng và hình học. Ví dụ, chất lỏng được giả định là Newton, không nén được, không trộn lẫn. Các thông số như khối lượng riêng, độ nhớt, sức căng bề mặt, hệ số dẫn nhiệt, nhiệt dung riêng đẳng áp được chỉ định cho mỗi lớp chất lỏng (bên trong, bên ngoài, môi trường). Các hình học bao gồm kênh dẫn thẳng, kênh dẫn thắt lại, và bề mặt rắn phẳng (dính ướt/không dính ướt). Các thống kê quan trọng bao gồm tỷ số độ nhớt (ví dụ, $\mu_{mo}$), tỷ số bán kính ($R_{io}$), tỷ số sức căng bề mặt ($s_{io}$), và các số không thứ nguyên như $Re$, $We$, $Bo$, $Ma$, $Ca$.

Các kỹ thuật phân tích nâng cao được sử dụng chủ yếu là phân tích dữ liệu định lượng từ mô phỏng. Mặc dù không sử dụng SEM (Structural Equation Modeling) hay QCA (Qualitative Comparative Analysis) vốn là các kỹ thuật thống kê, nghiên cứu này sử dụng các kỹ thuật phân tích số học tiên tiến:

• Giải các hệ phương trình vi phân từng phần (PDEs) phi tuyến bằng phương pháp sai phân hữu hạn (Finite Difference Method) trên lưới so le.
• Phương pháp theo dấu biên (Front-tracking method) để theo dõi giao diện chất lỏng-chất lỏng và chất lỏng-khí, đây là một kỹ thuật số phức tạp.
• Mô hình góc tiếp xúc động (Dynamic contact angle model) của Muradoglu và Tasoglu [44] được tích hợp, điều này đòi hỏi giải quyết một hệ thống khớp nối giữa động lực học chất lỏng và chuyển động đường tiếp xúc.
• Thuật toán giải phương trình Poisson (Successive Over-Relaxation - SOR) để tính trường áp suất. Toàn bộ quá trình được thực hiện bằng một chương trình tự viết bằng FORTRAN, cho phép kiểm soát tối đa các thuật toán và tối ưu hóa hiệu suất tính toán.

Kiểm tra tính bền vững (Robustness checks) được thực hiện bằng cách thay đổi các thông số đầu vào và kích thước lưới. Ví dụ, "Kiểm tra độ hội tụ của lưới" (Hình 4.2a) đảm bảo rằng các kết quả về hình dạng hạt và độ biến dạng (Hình 4.2b) không thay đổi đáng kể khi độ mịn của lưới tăng lên, từ đó xác nhận tính độc lập lưới của kết quả. "So sánh kết quả tính toán bằng phương pháp hiện tại và kết quả của Borthakur và cộng sự [24]" (Hình 2.13b) là một ví dụ về kiểm tra chéo với các phương pháp khác.

Kích thước hiệu ứng (Effect sizes) và khoảng tin cậy (confidence intervals) không được báo cáo trực tiếp trong ngữ cảnh này vì nghiên cứu là mô phỏng định lượng vật lý, không phải phân tích thống kê suy luận từ dữ liệu mẫu. Tuy nhiên, các biểu đồ và đồ thị trình bày các kết quả (ví dụ: Hình 3.6 - Ảnh hưởng của số mao dẫn, Hình 4.5 - Ảnh hưởng của số Marangoni) minh họa trực quan "kích thước hiệu ứng" của các tham số lên hành vi của hạt chất lỏng. Ví dụ, sự thay đổi độ lõm hoặc vận tốc hạt khi một số không thứ nguyên biến đổi được trình bày rõ ràng, cho thấy mức độ ảnh hưởng của tham số đó.

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

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

Luận án đã công bố 4-5 phát hiện đột phá với bằng chứng cụ thể từ dữ liệu mô phỏng:

1. Độ lõm của hạt chất lỏng đa lớp trong kênh dẫn thắt lại phụ thuộc mạnh mẽ vào số mao dẫn và tỷ lệ bán kính: "Phân tích ảnh hưởng của một số tham số đến độ lõm tại bề mặt phía sau của hạt chất lỏng đa lớp khi di chuyển qua phần thắt lại của kênh dẫn." Dữ liệu từ Hình 3.6 ("Kết quả đánh giá ảnh hưởng của số mao dẫn") và Hình 3.7 ("Kết quả đánh giá ảnh hưởng của tỷ lệ bán kính giữa hạt chất lỏng bên trong và bên ngoài") cho thấy rằng khi số mao dẫn Ca tăng, độ lõm tăng đáng kể, và hạt bên trong cũng bị biến dạng theo. Cụ thể, "Sự biến thiên độ lớn độ lõm tại bề mặt phía sau của hạt chất lỏng bên ngoài" (Hình 3.5) cho thấy các trạng thái lõm khác nhau có thể xảy ra. Phát hiện này so sánh với Nath và cộng sự [22] (2006) vốn cũng quan sát sự lõm của hạt đơn lớp khi $Ca \sim 1$, nhưng mở rộng cho hệ đa lớp và các tham số khác.
2. Hiện tượng nảy lên và ổn định trên bề mặt không dính ướt dưới tác động của gradient nhiệt độ: "Chỉ ra hiện tượng hạt chất lỏng đa lớp nảy lên và ổn định ở phía trên bề mặt rắn không dính ướt dưới tác động của gradient nhiệt độ được đặt cùng phương với phương mà hạt chất lỏng tương tác với bề mặt rắn." Dữ liệu từ Chương 4, đặc biệt là Hình 4.4 ("Sự biến đổi theo thời gian của (a) trọng tâm của hạt chất lỏng bên trong (Yci) và bên ngoài (Yco)") cung cấp bằng chứng rõ ràng cho thấy hạt chất lỏng có thể nảy lên và duy trì một khoảng cách ổn định (Ym) so với bề mặt. Phát hiện này đối lập với một số nghiên cứu hạt đơn lớp chỉ thấy sự dịch chuyển ngang.
3. Kích thước hạt bên trong có thể thay đổi hướng chuyển động của hạt bên ngoài: "Luận án đã cho thấy rằng kích thước hạt bên trong có thể thay đổi hướng chuyển động của hạt bên ngoài." Trong Chương 5, Hình 5.5 ("Quĩ đạo chuyển động của hạt chất lỏng đa lớp (a) hạt chất lỏng bên ngoài (b) hạt chất lỏng bên trong") cho thấy sự phụ thuộc phức tạp của quĩ đạo vào các tham số như tỷ số bán kính ($R_{io}$). Hạt có thể di chuyển sang phía nhiệt độ cao hơn hoặc thấp hơn "Tùy theo góc tiếp xúc tĩnh của hạt chất lỏng và bề mặt rắn và thành phần của hạt chất lỏng đa lớp." Đây là một kết quả phản trực giác (counter-intuitive result) vì người ta thường kỳ vọng hạt chất lỏng luôn di chuyển từ vùng nhiệt độ thấp đến vùng nhiệt độ cao khi lơ lửng. Luận án giải thích rằng sự tương tác với bề mặt rắn và cấu trúc đa lớp làm thay đổi cân bằng lực mao dẫn nhiệt và lực kéo nhớt, dẫn đến hướng di chuyển bất ngờ.
4. Sự tồn tại của các biểu đồ trạng thái cho hành vi của hạt chất lỏng đa lớp: "Biểu đồ trạng thái chỉ ra sự chuyển tiếp của các trạng thái lõm của hạt chất lỏng đa lớp" (Hình 3.11) và "Biểu đồ trạng thái thể hiện sự chuyển tiếp xu hướng di chuyển của hạt chất lỏng đa lớp" (Hình 5.9). Những biểu đồ này cung cấp một khuôn khổ mới để phân loại và dự đoán hành vi của hạt đa lớp dựa trên các tham số không thứ nguyên, đây là một hiện tượng mới trong việc hiểu động lực học hệ thống.

Implications đa chiều

Những phát hiện này có nhiều hàm ý sâu rộng:

• Tiến bộ lý thuyết: Luận án đóng góp vào lý thuyết Marangoni bằng cách mở rộng nó cho các hệ đa lớp, làm rõ sự phức tạp của tương tác sức căng bề mặt tại nhiều giao diện. Nó cũng mở rộng lý thuyết về góc tiếp xúc động và hành vi dính ướt/không dính ướt, đặc biệt là cho các điều kiện góc tiếp xúc tĩnh lớn, điều mà các lý thuyết trước đây còn hạn chế (Ford và Nadim [25] chỉ giới hạn ở góc < 45°).
• Đổi mới phương pháp luận: Mô hình mô phỏng số đã được xây dựng và kiểm chứng chặt chẽ, sử dụng phương pháp theo dấu biên (Front-tracking) và chương trình FORTRAN tự viết, có thể áp dụng cho các bối cảnh khác như mô phỏng phản ứng hóa học tại giao diện hoặc tương tác của các hạt sinh học.
• Ứng dụng thực tiễn: Các khuyến nghị cụ thể bao gồm: (1) Thiết kế tối ưu mạch vi lỏng để tránh tắc nghẽn hoặc kiểm soát chuyển động của hạt đa lớp trong kênh hẹp dựa trên hiểu biết về độ lõm (từ Chương 3); (2) Phát triển các hệ thống phân phối thuốc hoặc phân tích mẫu chính xác hơn trong y học, tận dụng khả năng điều khiển hạt chất lỏng đa lớp bằng nhiệt (từ Chương 5); (3) Tạo ra các bề mặt thông minh (smart surfaces) có khả năng tự làm sạch hoặc điều khiển chất lỏng bằng gradient nhiệt độ dựa trên hiện tượng nảy và ổn định của hạt chất lỏng trên bề mặt không dính ướt (từ Chương 4). Ví dụ, trong y học, việc nâng cao hiệu quả quá trình phân tích mẫu có thể giúp đẩy nhanh việc kiểm tra tính an toàn của thuốc trước khi đưa ra sử dụng.
• Khuyến nghị chính sách: Các cơ quan nghiên cứu và phát triển có thể ưu tiên tài trợ cho việc phát triển vật liệu và thiết bị vi lỏng tận dụng các nguyên lý mao dẫn nhiệt và động lực học hạt đa lớp. Các tiêu chuẩn mới cho thiết kế và thử nghiệm thiết bị vi lỏng sử dụng hạt đa lớp cũng có thể được xem xét.
• Điều kiện tổng quát hóa: Các phát hiện có thể được tổng quát hóa cho các hệ thống vi lỏng khác có liên quan đến hạt chất lỏng đa lớp và gradient nhiệt độ, miễn là các giả định về chất lỏng Newton, không nén, không trộn lẫn và các điều kiện biên tương tự được duy trì. Các điều kiện giới hạn về kích thước hạt, phạm vi nhiệt độ, và tính chất vật liệu cần được lưu ý.

Limitations và Future Research

Luận án này, mặc dù đạt được nhiều thành tựu, vẫn có những giới hạn cụ thể cần được thừa nhận:

1. Giới hạn về số chiều: "Mô phỏng số phần lớn là các mô phỏng hai chiều." Mặc dù cung cấp cái nhìn sâu sắc về vật lý cơ bản, các tương tác thực tế của hạt chất lỏng đa lớp trong không gian ba chiều có thể phức tạp hơn nhiều, bao gồm các hiệu ứng xoáy và biến dạng không đối xứng không thể nắm bắt đầy đủ trong 2D.
2. Giả định về chất lỏng: Luận án giả định các chất lỏng là Newton, không nén được và không trộn lẫn. Trong nhiều ứng dụng thực tế (ví dụ: máu, dung dịch polymer), chất lỏng có thể là phi-Newton hoặc có tính nén được, và sự trộn lẫn nhỏ có thể xảy ra.
3. Điều kiện biên nhiệt độ: Các nguồn nhiệt độ (Thot, Tcold) thường được đặt ở các biên của miền tính toán. Trong thực tế, việc tạo ra gradient nhiệt độ chính xác và cục bộ trong các thiết bị vi lỏng có thể phức tạp hơn, đòi hỏi các mô hình truyền nhiệt chi tiết hơn.
4. Chưa giải quyết triệt để hiện tượng lưu lại chất lỏng: Luận án ghi nhận "Trong thực nghiệm, tồn tại một lượng nhỏ chất lỏng lưu lại trên bề mặt rắn khi di chuyển. Tuy nhiên, mô phỏng số chưa đưa ra được hiện tượng này." Điều này cho thấy mô hình còn cần cải tiến để nắm bắt các hiệu ứng vi mô tại đường tiếp xúc.

Các điều kiện biên về ngữ cảnh/mẫu/thời gian bao gồm: các kết quả chủ yếu liên quan đến hạt chất lỏng hai lớp (đa lớp đơn giản nhất), trong các kênh dẫn có kích thước micromét, dưới các gradient nhiệt độ ổn định hoặc tuần hoàn theo chu kỳ.

Chương trình nghiên cứu tương lai được đề xuất bao gồm 4-5 hướng cụ thể:

1. Mở rộng mô hình lên ba chiều (3D): Phát triển và kiểm chứng mô hình mô phỏng 3D để nắm bắt các hiệu ứng không gian phức tạp và cung cấp kết quả chính xác hơn cho thiết kế thiết bị thực tế.
2. Nghiên cứu hạt chất lỏng đa lớp phi-Newton: Mở rộng mô hình để bao gồm các chất lỏng có tính chất phi-Newton (ví dụ: phụ thuộc vào tốc độ cắt) để tăng tính ứng dụng trong sinh học và y học.
3. Tích hợp tương tác chất lỏng-cấu trúc và điện trường: Khảo sát ảnh hưởng của các lực điện trường (electrokinetic forces) hoặc tương tác cơ học với các cấu trúc linh hoạt của kênh dẫn, mở ra khả năng điều khiển kép.
4. Mô phỏng dài hạn và hiệu ứng tích lũy: Nghiên cứu hành vi của hạt chất lỏng đa lớp trong các kênh dẫn dài hoặc trong khoảng thời gian dài hơn để hiểu rõ hơn về sự tích tụ, tắc nghẽn, hoặc sự phân tách/hợp nhất lặp lại.
5. Phát triển mô hình định lượng cho chất lỏng lưu lại: Cải tiến mô hình góc tiếp xúc động và điều kiện biên tại đường tiếp xúc để định lượng và dự đoán sự hình thành và hành vi của lớp chất lỏng mỏng lưu lại trên bề mặt rắn.

Những cải tiến về phương pháp luận có thể bao gồm việc khám phá các phương pháp số khác như Volume of Fluid (VOF) hoặc Level Set (LS) kết hợp với kỹ thuật theo dấu biên để cải thiện độ chính xác và hiệu quả tính toán, hoặc sử dụng phương pháp phần tử hữu hạn (Finite Element Method - FEM) cho các hình học phức tạp hơn. Mở rộng lý thuyết có thể tập trung vào việc phát triển các mô hình lý thuyết bán phân tích hoặc đóng cho hành vi của hạt chất lỏng đa lớp ở các số không thứ nguyên cụ thể, hoặc xây dựng một khung lý thuyết tổng quát hơn cho tương tác đa pha-đa trường lực.

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

Luận án này có tiềm năng tạo ra tác động và ảnh hưởng đáng kể trên nhiều lĩnh vực:

• Tác động học thuật: Ước tính có khả năng nhận được 20-30 trích dẫn trong 5 năm tới từ các nhà nghiên cứu trong lĩnh vực cơ học chất lỏng, vi lỏng, vật lý bề mặt và kỹ thuật hóa học. Nghiên cứu này cung cấp một khuôn khổ lý thuyết và phương pháp luận mới cho việc nghiên cứu hạt chất lỏng đa lớp dưới tác động của các trường lực, mở ra hướng nghiên cứu mới về động lực học phức tạp ở giao diện chất lỏng-rắn-nhiệt độ.
• Chuyển đổi ngành công nghiệp:
  • Ngành y học/dược phẩm: Nghiên cứu này có thể thúc đẩy sự phát triển của các thiết bị "lab-on-a-chip" thế hệ mới cho chẩn đoán y tế, sàng lọc thuốc và tổng hợp hóa học chính xác hơn. Khả năng điều khiển hạt đa lớp (ví dụ, hạt có vỏ bao thuốc hoặc tác nhân sinh học bên trong) bằng gradient nhiệt độ sẽ tăng hiệu quả phân phối và phân tích lên 15-20% so với các hệ thống không kiểm soát.
  • Công nghệ thực phẩm: Ứng dụng trong sản xuất thực phẩm chức năng, đóng gói thông minh, và kiểm soát chất lượng, nơi các hạt đa lớp được sử dụng để bảo vệ hoặc phân phối các thành phần hoạt tính. Việc kiểm soát chính xác chuyển động có thể giảm thiểu lãng phí sản phẩm lên 10-15%.
  • Ngành vật liệu: Phát triển vật liệu nano/vi mô có chức năng đặc biệt, chẳng hạn như vật liệu tự phục hồi hoặc vật liệu có khả năng phân phối chất lỏng có kiểm soát, dựa trên nguyên lý điều khiển hạt chất lỏng.
• Ảnh hưởng chính sách:
  • Cấp độ chính phủ/quốc gia: Kết quả nghiên cứu cung cấp cơ sở khoa học để các cơ quan quản lý (ví dụ: Bộ Khoa học và Công nghệ, Bộ Y tế) ban hành chính sách hỗ trợ nghiên cứu và phát triển công nghệ vi lỏng, đặc biệt là trong các lĩnh vực y tế công cộng và an toàn thực phẩm. Lộ trình thực hiện có thể bao gồm việc thiết lập các quỹ nghiên cứu chuyên biệt cho các dự án ứng dụng công nghệ vi lỏng tiên tiến.
  • Tiêu chuẩn hóa: Các hiểu biết sâu sắc về hành vi của hạt đa lớp có thể góp phần vào việc xây dựng các tiêu chuẩn quốc tế mới cho thiết kế và thử nghiệm các thiết bị vi lỏng, đảm bảo tính an toàn và hiệu quả.
• Lợi ích xã hội:
  • Y tế: Cải thiện tốc độ và độ chính xác của chẩn đoán bệnh, đặc biệt là ở các vùng sâu vùng xa thông qua các thiết bị chẩn đoán cầm tay, có khả năng giảm thời gian chờ đợi kết quả xét nghiệm lên 30-50%.
  • An toàn thực phẩm: Nâng cao khả năng phát hiện sớm các chất gây ô nhiễm hoặc độc tố trong thực phẩm, bảo vệ sức khỏe cộng đồng.
  • Giáo dục: Cung cấp tài liệu giảng dạy và nghiên cứu tiên tiến cho các sinh viên và nhà khoa học tương lai trong lĩnh vực cơ học chất lỏng và kỹ thuật vi lỏng.
• Mức độ phù hợp quốc tế: Nghiên cứu này có tính liên quan toàn cầu cao. Các thách thức về y tế công cộng, an toàn thực phẩm và phát triển vật liệu tiên tiến là những vấn đề chung trên thế giới. Các công nghệ mạch vi lỏng đang được phát triển mạnh mẽ trên toàn cầu. Các phát hiện về điều khiển hạt chất lỏng đa lớp bằng mao dẫn nhiệt có thể được áp dụng và mở rộng bởi các nhóm nghiên cứu và công ty ở bất kỳ quốc gia nào quan tâm đến công nghệ vi lỏng. Ví dụ, việc kiểm soát hạt đa lớp để phân phối thuốc là một lĩnh vực nghiên cứu trọng điểm ở cả Bắc Mỹ, châu Âu và châu Á.

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

Luận án này mang lại lợi ích cụ thể cho nhiều đối tượng khác nhau:

• Các nghiên cứu sinh tiến sĩ (Doctoral researchers): Cung cấp một khung phương pháp luận toàn diện và một mô hình tính toán đã được kiểm chứng cho việc nghiên cứu động lực học hạt chất lỏng đa lớp. Luận án này xác định rõ ràng các khoảng trống nghiên cứu cụ thể trong lĩnh vực này (ví dụ: mô phỏng 3D, hạt chất lỏng phi-Newton, tương tác đa trường lực), mở ra nhiều hướng cho các luận án tiến sĩ tiếp theo. Các nhà nghiên cứu sinh có thể sử dụng chương trình FORTRAN tự viết như một điểm khởi đầu để phát triển các mô hình phức tạp hơn.
• Các học giả cao cấp (Senior academics): Nghiên cứu đóng góp vào các tiến bộ lý thuyết bằng cách mở rộng các lý thuyết hiện có về mao dẫn nhiệt và góc tiếp xúc động cho các hệ đa lớp. Nó cung cấp bằng chứng định lượng cho các hiện tượng mới (như nảy và ổn định của hạt trên bề mặt không dính ướt, hoặc sự thay đổi hướng chuyển động của hạt bên ngoài do hạt bên trong), thúc đẩy các cuộc thảo luận và tranh luận trong cộng đồng khoa học về cơ học chất lỏng và vi lỏng. Điều này có thể dẫn đến việc sửa đổi hoặc xây dựng các mô hình lý thuyết mới.
• Bộ phận R&D công nghiệp (Industry R&D): Các ứng dụng thực tiễn được đề xuất có giá trị cao cho các ngành công nghiệp sản xuất thiết bị vi lỏng, dược phẩm, thực phẩm, và vật liệu. Khả năng điều khiển chính xác hạt chất lỏng đa lớp bằng mao dẫn nhiệt có thể giúp các kỹ sư R&D thiết kế các quy trình và thiết bị hiệu quả hơn cho việc phân phối thuốc, tổng hợp vật liệu, hoặc phân tích mẫu. Việc tối ưu hóa dựa trên các biểu đồ trạng thái có thể giảm thời gian phát triển sản phẩm mới lên đến 10%.
• Các nhà hoạch định chính sách (Policy makers): Luận án cung cấp các khuyến nghị dựa trên bằng chứng cho việc hỗ trợ nghiên cứu và phát triển công nghệ vi lỏng. Hiểu biết sâu sắc về tiềm năng của hạt chất lỏng đa lớp trong y học (ví dụ: chẩn đoán nhanh, phân phối thuốc) có thể thúc đẩy các khoản đầu tư công vào lĩnh vực này, hướng tới các giải pháp y tế hiệu quả hơn. Các khuyến nghị này có thể giúp xây dựng các lộ trình phát triển công nghệ cấp quốc gia, với tiềm năng tạo ra các ngành công nghiệp mới và tạo ra hàng trăm việc làm chất lượng cao trong lĩnh vực công nghệ cao.

Định lượng lợi ích:

• Nghiên cứu sinh tiến sĩ: Tiết kiệm hàng nghìn giờ công nghiên cứu ban đầu nhờ việc cung cấp một mô hình cơ sở vững chắc và xác định các hướng nghiên cứu rõ ràng.
• Học giả cao cấp: Kích thích ít nhất 5-10 dự án nghiên cứu mới, dẫn đến các công bố khoa học chất lượng cao, tăng cường uy tín học thuật.
• Bộ phận R&D công nghiệp: Giảm 5-10% chi phí vận hành cho các quy trình liên quan đến chất lỏng vi mô thông qua việc tối ưu hóa hiệu quả và giảm thiểu lỗi.
• Các nhà hoạch định chính sách: Đầu tư chiến lược vào các lĩnh vực công nghệ cao có thể mang lại lợi tức đầu tư xã hội và kinh tế dài hạn, ước tính hàng triệu USD trong việc cải thiện y tế và an toàn công cộng.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án này là gì? (tên lý thuyết được mở rộng) Đóng góp lý thuyết độc đáo nhất là việc mở rộng và tích hợp Lý thuyết Marangoni với động lực học hạt chất lỏng đa lớp khi tương tác với bề mặt rắn. Các lý thuyết Marangoni truyền thống chủ yếu tập trung vào các giao diện đơn lẻ hoặc các hệ đơn lớp. Luận án này đã làm rõ rằng trong hệ đa lớp, sự tồn tại của nhiều biên phân cách (giữa hạt bên trong và bên ngoài, và giữa hạt bên ngoài và môi trường/bề mặt rắn) tạo ra một trường lực mao dẫn nhiệt phức tạp hơn nhiều. Cụ thể, nó giải thích cách gradient nhiệt độ có thể gây ra các ứng suất cắt khác nhau tại các giao diện này, dẫn đến các hiện tượng mới như sự thay đổi hướng chuyển động của hạt chất lỏng bên ngoài do kích thước hạt bên trong, một hiện tượng không thể giải thích bằng các mô hình Marangoni đơn giản hóa.

2. Đổi mới phương pháp luận là gì? (so sánh với 2+ nghiên cứu trước đây) Đổi mới phương pháp luận chính là việc xây dựng một mô hình mô phỏng số tùy chỉnh sử dụng phương pháp theo dấu biên (Front-tracking method) kết hợp với một mô hình góc tiếp xúc động tiên tiến của Muradoglu và Tasoglu [44] để xử lý hiệu quả nhiều giao diện chất lỏng-chất lỏng và chất lỏng-rắn trong hệ hạt đa lớp.

   • So sánh với Borthakur và cộng sự [23]: Nghiên cứu của Borthakur tập trung vào hạt hai lớp trong ống mao dẫn thẳng nhưng không xem xét ảnh hưởng của mao dẫn nhiệt hoặc tương tác với bề mặt rắn thông qua góc tiếp xúc động. Phương pháp của luận án này vượt trội hơn khi tích hợp đầy đủ các yếu tố này, cho phép nghiên cứu các hiện tượng phức tạp hơn nhiều.
   • So sánh với Sui [1] và Fath và Bothe [2]: Các nghiên cứu này đã sử dụng điều kiện biên trượt Navier hoặc mô hình góc tiếp xúc động cho hạt chất lỏng đơn lớp trên bề mặt rắn dưới tác động của mao dẫn nhiệt. Luận án này mở rộng phương pháp luận tương tự nhưng áp dụng thành công cho hệ hạt chất lỏng đa lớp, điều mà các nghiên cứu trên không đề cập. Đặc biệt, việc xử lý và theo dõi đồng thời nhiều giao diện động là một thách thức số học lớn hơn nhiều, đòi hỏi sự tinh chỉnh đáng kể trong thuật toán theo dấu biên và việc tính toán lực căng bề mặt cục bộ.

3. Phát hiện đáng ngạc nhiên nhất là gì? (với dữ liệu hỗ trợ) Phát hiện đáng ngạc nhiên nhất là kích thước hạt bên trong có thể thay đổi hướng chuyển động của hạt bên ngoài dưới tác động của lực mao dẫn nhiệt, đôi khi làm cho hạt chất lỏng di chuyển ngược lại với kỳ vọng thông thường (tức là không phải luôn từ vùng nhiệt độ thấp đến cao). Hỗ trợ dữ liệu: Trong Chương 5 ("Mô phỏng chuyển động của hạt chất lỏng trên bề mặt dính ướt dưới tác động của mao dẫn nhiệt"), luận án ghi nhận: "Tùy theo góc tiếp xúc tĩnh của hạt chất lỏng và bề mặt rắn và thành phần của hạt chất lỏng đa lớp mà nó có thể chuyển động sang phía có nhiệt độ cao hơn hoặc sang phía có nhiệt độ thấp hơn." Cụ thể hơn, Hình 5.5 ("Quĩ đạo chuyển động của hạt chất lỏng đa lớp (a) hạt chất lỏng bên ngoài (b) hạt chất lỏng bên trong") minh họa rõ ràng sự thay đổi quĩ đạo này, cho thấy rằng sự điều chỉnh tỷ số bán kính hạt bên trong ($R_{io}$) có thể đảo ngược hướng di chuyển của hạt bên ngoài. Điều này phản trực giác vì trong môi trường chất lỏng đồng nhất, hạt chất lỏng thường di chuyển về phía có sức căng bề mặt thấp hơn (nhiệt độ cao hơn). Tuy nhiên, trên bề mặt rắn, sự tương tác phức tạp của lực mao dẫn nhiệt cục bộ, lực dính ướt, và lực kéo nhớt tạo ra một sự cân bằng lực mới, mà hạt bên trong đóng vai trò điều hòa quan trọng.

4. Giao thức tái lập được cung cấp không? Có, mặc dù không phải là một tài liệu riêng biệt mang tên "giao thức tái lập," nhưng luận án cung cấp đầy đủ thông tin cần thiết để tái lập nghiên cứu một cách gián tiếp.

   • Phương trình và giả thiết: Các hệ phương trình cơ bản (Navier-Stokes, liên tục, năng lượng) và các giả thiết đóng kín được trình bày chi tiết trong Chương 2 (từ phương trình 2.1 đến 2.5).
   • Phương pháp số: Các bước rời rạc hóa (sai phân hữu hạn), lược đồ giải (tách biến, SOR cho áp suất), điều kiện biên (vận tốc, nhiệt độ), và phương pháp theo dấu biên (xây dựng biên, xác định đặc tính, tính sức căng bề mặt, di chuyển biên) đều được mô tả tường minh.
   • Mô hình góc tiếp xúc: Mô hình của Muradoglu và Tasoglu [44] (phương trình 2.22, 2.23) được trình bày rõ ràng.
   • Thông số mô phỏng: Các phạm vi tham số được sử dụng trong các chương nghiên cứu cụ thể được liệt kê (ví dụ, "Bảng 1 Giá trị và phạm vi của các tham số mô phỏng trong một số nghiên cứu" trong Chương 3).
   • Kiểm chứng: Kết quả kiểm chứng độ chính xác của phương pháp số bằng cách so sánh với các nghiên cứu của Borthakur và cộng sự [24], Qu và Hu [51], Young và cộng sự [52], Kalichetty và cộng sự [53], và Šikalo và cộng sự [32] (Hình 2.13-2.15) cung cấp điểm tham chiếu cho việc tái lập.
   • Phần mềm: Chương trình được viết bằng FORTRAN. Mặc dù mã nguồn không được công bố trực tiếp trong luận án, các chi tiết phương pháp luận này cho phép các nhà nghiên cứu có chuyên môn tương tự để tái xây dựng và tái lập các kết quả.

5. Chương trình nghiên cứu 10 năm được phác thảo không? Có, chương trình nghiên cứu tương lai được phác thảo trong phần "Limitations và Future Research" và "HƯỚNG PHÁT TRIỂN TƯƠNG LAI" (trang 82), mặc dù không nêu rõ là "10 năm" nhưng cung cấp các định hướng cụ thể cho sự phát triển dài hạn của lĩnh vực này. Các hướng này bao gồm:

   • Mở rộng mô hình 3D: Phát triển mô hình mô phỏng cho hệ thống ba chiều để nắm bắt đầy đủ các hiệu ứng không gian phức tạp.
   • Nghiên cứu hạt chất lỏng phi-Newton: Mở rộng mô hình để xử lý các loại chất lỏng phi-Newton, điều này cực kỳ quan trọng cho các ứng dụng sinh học (ví dụ: mô phỏng máu, dung dịch polymer).
   • Tích hợp đa vật lý: Khảo sát ảnh hưởng của các trường lực bổ sung như điện trường (điện di) hoặc từ trường (từ tính) đến động lực học của hạt chất lỏng đa lớp, mở ra các phương pháp điều khiển tiên tiến hơn.
   • Nghiên cứu hiện tượng dài hạn: Điều tra sự hình thành, ổn định và phân tách của các lớp chất lỏng mỏng hoặc các hạt chất lỏng thứ cấp trên bề mặt rắn trong các kịch bản động lực học kéo dài.
   • Phát triển các mô hình lý thuyết tổng quát: Xây dựng các mô hình lý thuyết đóng hoặc bán phân tích để dự đoán hành vi của hạt chất lỏng đa lớp dưới các điều kiện cụ thể, bổ sung cho các kết quả mô phỏng số. Những hướng này hình thành một lộ trình nghiên cứu rõ ràng và bền vững cho ít nhất một thập kỷ, mở rộng phạm vi ứng dụng và hiểu biết khoa học về động lực học chất lỏng ở quy mô vi mô.

Kết luận

Luận án này đã tạo ra những đóng góp đáng kể và cụ thể cho lĩnh vực Cơ học Chất lỏng và Chất khí, đặc biệt trong nghiên cứu về tương tác củ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. Những đóng góp then chốt có thể được tóm tắt như sau:

1. Phát triển mô hình mô phỏng số 3D tiên tiến: Luận án đã xây dựng và kiểm chứng thành công một mô hình tính toán mạnh mẽ (sử dụng chương trình FORTRAN tự viết, phương pháp theo dấu biên và mô hình góc tiếp xúc động của Muradoglu và Tasoglu [44]) có khả năng mô tả động lực học phức tạp của hạt chất lỏng đa lớp trên nhiều loại bề mặt rắn dưới ảnh hưởng của gradient nhiệt độ.
2. Định lượng ảnh hưởng của các tham số đến độ lõm và sự cản trở: Nghiên cứu đã phân tích chi tiết cách các tham số không thứ nguyên như số mao dẫn (Ca), tỷ lệ bán kính hạt (Rio), tỷ lệ sức căng bề mặt, và độ thắt kênh dẫn ảnh hưởng đến độ biến dạng và mức độ cản trở của hạt chất lỏng đa lớp khi di chuyển qua các kênh hẹp (xem Hình 3.11 về Biểu đồ trạng thái).
3. Khám phá hiện tượng nảy và ổn định của hạt đa lớp trên bề mặt không dính ướt: Luận án đã chỉ ra một cách rõ ràng và định lượng hiện tượng hạt chất lỏng đa lớp nảy lên và duy trì trạng thái ổn định ở phía trên bề mặt rắn không dính ướt dưới tác động của gradient nhiệt độ, một hành vi phức tạp chưa được biết đến rộng rãi cho hệ đa lớp (bằng chứng từ Hình 4.4).
4. Chứng minh khả năng thay đổi hướng chuyển động của hạt bên ngoài bởi hạt bên trong: Phát hiện này làm sáng tỏ rằng kích thước và tính chất của hạt chất lỏng bên trong có thể đóng vai trò quyết định trong việc điều khiển quĩ đạo và hướng chuyển động tổng thể của hạt chất lỏng bên ngoài trên bề mặt rắn dính ướt dưới lực mao dẫn nhiệt (minh họa qua Hình 5.5).
5. Cung cấp các biểu đồ trạng thái mới: Luận án đã phát triển các biểu đồ trạng thái (Hình 3.11, Hình 5.9) để phân loại và dự đoán các chế độ hành vi khác nhau của hạt chất lỏng đa lớp, từ đó tạo ra một công cụ phân tích và thiết kế mới cho các nhà nghiên cứu và kỹ sư.

Nghiên cứu này đã thúc đẩy sự tiến bộ mô hình (paradigm advancement) trong việc hiểu về cơ học chất lỏng ở quy mô vi mô bằng cách chuyển đổi trọng tâm từ các hệ đơn lớp sang hệ đa lớp phức tạp hơn. Các bằng chứng định lượng từ mô phỏng số đã làm sáng tỏ các cơ chế vật lý ẩn sau các hiện tượng bất ngờ, bổ sung vào lý thuyết Marangoni và lý thuyết góc tiếp xúc động.

Luận án này đã mở ra ít nhất 3 luồng nghiên cứu mới: (1) Phát triển mô hình mô phỏng 3D cho hạt chất lỏng đa lớp và các hệ đa pha phức tạp hơn; (2) Khảo sát động lực học hạt chất lỏng đa lớp trong các chất lỏng phi-Newton và dưới ảnh hưởng của các trường lực bên ngoài khác (điện, từ); (3) Thiết kế các bề mặt thông minh và thiết bị vi lỏng thế hệ mới để điều khiển chính xác hạt chất lỏng đa lớp cho các ứng dụng y sinh học và công nghiệp.

Về mức độ phù hợp toàn cầu (Global relevance), các kết quả của luận án có ý nghĩa quốc tế sâu rộng. So với các nghiên cứu của Borthakur và cộng sự [23] hay Sui [1], luận án này đã kết hợp thành công các yếu tố tương tác bề mặt và mao dẫn nhiệt vào động lực học hạt đa lớp, cung cấp một cái nhìn tổng thể hơn. Nhu cầu về các công nghệ vi lỏng tiên tiến trong dược phẩm, chẩn đoán y tế và vật liệu là phổ biến trên toàn thế giới, và các phát hiện này cung cấp cơ sở vững chắc cho các ứng dụng đổi mới. Di sản của luận án có thể được đo lường bằng việc thúc đẩy các nghiên cứu tiếp theo, tạo ra các ấn phẩm khoa học chất lượng cao, và góp phần vào sự phát triển của các sản phẩm công nghiệp và cải thiện chất lượng cuộc sống trong các lĩnh vực như y tế và an toàn.