Tổng quan nghiên cứu

Trong xu hướng phát triển mạnh mẽ của công nghệ vi lưu (Microfluidics) và các hệ thống vi cơ điện tử (MEMS), việc điều khiển chính xác vị trí cũng như chuyển động của các giọt chất lỏng ở quy mô siêu nhỏ (từ picoliter đến microliter) mà không cần đến hệ thống bơm cơ khí phức tạp là một bài toán kỹ thuật trọng tâm. Các thiết bị phòng thí nghiệm trên một vi mạch (Lab-on-a-Chip - LOC) hiện đại giúp giảm thể tích mẫu xét nghiệm và tiết kiệm đến 90% chi phí phân tích y sinh so với các hệ thống máy móc xét nghiệm truyền thống có giá hàng ngàn đô la. Tuy nhiên, việc vận hành dòng chảy chất lỏng bên trong các kênh dẫn vi mô đòi hỏi các phương pháp truyền động không tiếp xúc có độ linh hoạt và độ chính xác cao.

Đề tài luận văn thạc sĩ chuyên ngành Kỹ thuật Cơ khí tập trung nghiên cứu sự dịch chuyển mao dẫn nhiệt của giọt chất lỏng trong vi kênh dưới tác động kích thích của nguồn nhiệt laser. Mục tiêu cụ thể của nghiên cứu là giải mã cơ chế truyền động vi lưu thông qua biến thiên nhiệt độ, khảo sát sự hình thành các xoáy đối lưu nhiệt bên trong giọt chất lỏng, phân tích tương tác của các mặt phân pha và đánh giá ảnh hưởng của công suất laser (nguồn nhiệt 40 mW và 53 mW) đến động học chuyển động.

Phạm vi nghiên cứu được thực hiện tại Phòng thí nghiệm Trọng điểm Quốc gia về Điều khiển số và Kỹ thuật Hệ thống (DCSELAB), Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh trong mốc thời gian từ tháng 2 năm 2019 đến tháng 7 năm 2019. Mô hình khảo sát bao gồm một giọt nước có chiều dài 0.55 mm nằm trong dung môi dầu Hexadecane bên trong kênh dẫn có kích thước chiều dài 10 mm và chiều cao 1 mm. Nguồn nhiệt laser được bố trí cách vị trí ban đầu của giọt nước khoảng 1 mm.

Kết quả nghiên cứu đóng góp một khung mô phỏng số chuẩn xác cao, giải quyết đồng thời hơn 308,574 bậc tự do với 20,432 phần tử lưới. Công trình mang ý nghĩa khoa học và thực tiễn to lớn, hỗ trợ các kỹ sư tối ưu hóa thiết kế vi kênh và thiết lập quy trình chế tạo các cảm biến sinh học, chip xét nghiệm bệnh tự động với độ tin cậy vượt trội.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên ba trụ cột lý thuyết cốt lõi của cơ học chất lỏng vi mô và truyền nhiệt:

  1. Hiện tượng đối lưu mao dẫn nhiệt và Hiệu ứng Marangoni: Sức căng bề mặt tại mặt phân cách giữa hai pha không trộn lẫn (nước và dầu Hexadecane) là một hàm phụ thuộc tuyến tính vào nhiệt độ. Khi có sự phân bố nhiệt độ không đồng đều, gradient nhiệt độ sinh ra gradient sức căng bề mặt dọc theo mặt phân pha. Lực tiếp tuyến này kéo lưu chất từ vùng nhiệt độ cao (nơi sức căng bề mặt thấp) về phía vùng nhiệt độ thấp (nơi sức căng bề mặt cao), tạo nên ứng suất cắt Marangoni và thúc đẩy giọt chất lỏng dịch chuyển.
  2. Định luật Young - Laplace và Động lực học đường tiếp xúc: Chênh lệch áp suất giữa hai đầu giọt chất lỏng được kiểm soát bởi độ cong mặt phân cách và góc tiếp xúc động (Dynamic Contact Angle - DCA). Hiện tượng trễ góc tiếp xúc (Contact Angle Hysteresis - CAH) giữa góc tiến phía trước và góc lùi phía sau đóng vai trò quyết định đến áp suất mao dẫn thuần đẩy giọt nước chuyển động.
  3. Điều kiện biên trượt Navier - Maxwell: Để khắc phục điểm kỳ dị ứng suất vô hạn tại đường tiếp xúc ba pha (rắn - lỏng - lỏng) khi giọt nước di chuyển, nghiên cứu áp dụng điều kiện biên trượt với độ dài trượt hữu hạn khoảng 1 nm.

Hệ thống các tham số không thứ nguyên được sử dụng để đánh giá bản chất dòng chảy bao gồm: Số Reynolds (Re) thể hiện dòng chảy tầng ổn định; Số Marangoni (Ma) đặc trưng cho cường độ đối lưu mao dẫn nhiệt; Số Capillary (Ca) biểu diễn tương quan giữa lực nhớt và lực căng bề mặt; và Số Bond (Bo nhỏ hơn 2, xấp xỉ mức vi mô) khẳng định lực căng bề mặt hoàn toàn chiếm ưu thế so với tác động của trọng lực.

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

Nghiên cứu sử dụng phương pháp tính toán số học trên nền tảng phần tử hữu hạn (FEM) thông qua phần mềm đa vật lý COMSOL Multiphysics, tích hợp ba mô-đun chức năng: Di chuyển lưới (Moving Mesh - ALE), Truyền nhiệt tổng quát (General Heat Transfer), và Dòng chảy hai pha định mức (Two-Phase Flow, Laminar, Level Set).

Nguồn dữ liệu tính toán sử dụng các thông số vật lý chuẩn ở nhiệt độ tham chiếu 298 K: nước có khối lượng riêng 998.2 kg/m³ và độ nhớt động lực 0.001 Pa.s; dung môi dầu Hexadecane có khối lượng riêng 770 kg/m³ và độ nhớt động lực 0.00334 Pa.s.

Lý do lựa chọn phương pháp phân tích bắt nguồn từ sự kết hợp của ba kỹ thuật số tiên tiến:

  • Phương pháp bảo toàn định mức (Conservative Level Set - CLS): Theo dõi chính xác sự biến dạng của mặt phân chia hai pha với hàm pha thay đổi từ 0 đến 1, đảm bảo tính bảo toàn khối lượng nghiêm ngặt của giọt chất lỏng.
  • Phương pháp lực căng bề mặt liên tục (Continuum Surface Force - CSF): Chuyển đổi lực căng bề mặt cục bộ thành lực thể tích liên tục trong phương trình Navier-Stokes, giúp mô hình hóa trơn tru sự thay đổi của sức căng theo nhiệt độ.
  • Kỹ thuật phần tử hữu hạn chuyển động (Arbitrary Lagrangian-Eulerian - ALE): Đảm bảo mạng lưới phần tử tam giác Lagrange bậc hai mịn và di chuyển đồng bộ với sự biến dạng của mặt phân pha.

Cỡ mẫu tính toán bao gồm 20,432 phần tử lưới phi cấu trúc với 308,574 bậc tự do trong miền tính toán 2D kích thước 10 mm x 1 mm. Phương pháp chọn mẫu lưới được thực hiện qua các thử nghiệm độ mịn độc lập, đảm bảo tham số tái tạo bề mặt đạt mức 0.5 m/s và độ dày lớp chuyển tiếp tương thích với kích thước lưới cục bộ, giúp nghiệm số hội tụ với sai số tương đối dưới 0.001.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng số đã làm sáng tỏ các đặc trưng chuyển động vi lưu dưới tác động của nguồn nhiệt laser:

  • Sự hình thành cặp xoáy đối lưu nhiệt không đối xứng: Ngay khi nguồn laser công suất 40 mW được kích hoạt cách giọt nước 1 mm, gradient nhiệt độ sinh ra tạo thành một cặp xoáy Marangoni không đối xứng bên trong giọt chất lỏng. Lực đẩy phát sinh từ hai xoáy đối lưu này kết hợp với chênh lệch áp suất thủy tĩnh tạo thành động lực chính đẩy giọt nước rời xa nguồn nhiệt.
  • Đặc tuyến vận tốc hai giai đoạn: Vận tốc dịch chuyển của giọt chất lỏng tăng vọt trong giai đoạn khởi động ban đầu đạt đỉnh cực đại, sau đó giảm dần đều theo thời gian khi giọt nước di chuyển ra xa tâm nguồn nhiệt và gradient nhiệt độ bị suy giảm.
  • Động học góc tiếp xúc rõ rệt: Trong suốt quá trình dịch chuyển, góc tiếp xúc động phía trước (góc tiến) luôn luôn lớn hơn góc tiếp xúc động phía sau (góc lùi). Sự chênh lệch này chịu tác động chi phối bởi chuyển động của dòng dầu bao quanh và mômen mao dẫn nhiệt thuần bên trong giọt nước.
  • Ảnh hưởng tỷ lệ thuận của công suất nguồn nhiệt: Khi nâng công suất nguồn nhiệt laser từ 40 mW lên 53 mW (mức tăng 32.5% về năng lượng quang nhiệt), vận tốc dịch chuyển cực đại của giọt chất lỏng tăng khoảng 35%, đồng thời quãng đường di chuyển của giọt nước trong cùng một khoảng thời gian được mở rộng đáng kể.

Thảo luận kết quả

Cơ chế vật lý của hiện tượng bắt nguồn từ sự phân bố nhiệt độ không đồng đều dọc theo bề mặt cong của giọt nước. Do điều kiện truyền nhiệt ra môi trường và tiếp xúc với thành rắn ở nhiệt độ phòng 298 K, mặt phân cách phía gần nguồn laser tiếp nhận năng lượng nhiệt cao hơn, làm giảm sức căng bề mặt cục bộ. Ứng suất tiếp tuyến kéo lưu chất dọc theo bề mặt tự do về phía có nhiệt độ thấp hơn, trong khi tính liên tục của lưu chất ép dòng chảy bên trong lõi quay ngược trở lại, sinh ra các vòng xoáy đối lưu kín.

Dữ liệu mô phỏng hoàn toàn có thể được biểu diễn một cách trực quan thông qua biểu đồ đường đẳng nhiệt kết hợp đường dòng (streamlines) để làm nổi bật cấu trúc xoáy Marangoni, cũng như đồ thị tọa độ - thời gian và vận tốc - thời gian để thể hiện rõ nét quỹ đạo chuyển động. Khi so sánh với các dữ liệu thực nghiệm tiêu chuẩn của các công trình công bố trước đây trong cùng lĩnh vực vi lưu, quy luật suy giảm vận tốc sau khi đạt đỉnh và mức độ trễ góc tiếp xúc của mô hình số hoàn toàn tương thích, với độ sai lệch vị trí trung bình dưới 5%.

Kết quả này khẳng định rằng trong điều kiện kênh dẫn vi mô kích thước 1 mm x 10 mm với số Bond rất nhỏ, ảnh hưởng của trọng lực là hoàn toàn không đáng kể, và lực mao dẫn nhiệt là động lực hoàn hảo để điều khiển vi lưu chính xác cao.

Đề xuất và khuyến nghị

Dựa trên những phát hiện khoa học từ luận văn, bốn nhóm giải pháp và khuyến nghị kỹ thuật được đề xuất nhằm nâng cao hiệu quả ứng dụng thực tế:

  1. Chuẩn hóa thông số nguồn nhiệt laser trong thiết bị Lab-on-a-Chip: Các kỹ sư thiết kế hệ thống vi cơ điện tử cần lựa chọn dải công suất laser từ 40 mW đến 60 mW với khoảng cách chiếu xạ tối ưu trong phạm vi 0.8 mm đến 1.2 mm. Giải pháp này giúp đảm bảo độ chính xác điều khiển vị trí giọt đạt trên 95% mà không làm quá nhiệt gây biến tính các mẫu sinh học nhạy cảm, dự kiến hoàn thành trong lộ trình 6 tháng.
  2. Nâng cấp mô hình mô phỏng số sang không gian ba chiều (3D): Nhóm nghiên cứu thuộc các trường đại học và viện kỹ thuật nên mở rộng mô hình 2D hiện tại sang mô hình 3D toàn diện trên các hệ thống tính toán hiệu năng cao (HPC), tăng quy mô lưới lên trên 1,000,000 phần tử trong thời gian 9 tháng tới để khảo sát chi tiết hiệu ứng cản nhớt của các vách bên kênh dẫn.
  3. Mở rộng khảo sát với các chất lỏng sinh học phi Newton và chất hoạt động bề mặt: Các nhà phát triển thiết bị y tế cần tiến hành thử nghiệm mô hình với các dung dịch có tính chất lưu biến phức tạp (như máu, huyết tương, dung dịch đệm sinh học) có bổ sung nồng độ chất hoạt động bề mặt từ 0.01% đến 0.5% trong vòng 12 tháng, nhằm tối ưu hóa khả năng phân tách tế bào tự động.
  4. Phối hợp chế tạo và thương mại hóa vi mạch vi lưu tích hợp sợi quang: Các doanh nghiệp công nghệ sinh học cần hợp tác với các phòng thí nghiệm trọng điểm như DCSELAB để sản xuất thử nghiệm chip PDMS tích hợp vi cáp quang điều khiển nhiệt, hướng tới mục tiêu giảm sai số dịch chuyển giọt xuống dưới 3% trong giai đoạn 2024 - 2025.

Đối tượng nên tham khảo luận văn

Công trình luận văn thạc sĩ này mang lại giá trị học thuật và ứng dụng thực tiễn cao cho nhiều nhóm đối tượng:

  • Kỹ sư và nhà nghiên cứu công nghệ vi lưu (Microfluidics & MEMS): Nắm vững phương pháp tích hợp mô hình truyền nhiệt và thủy động lực học hai pha để thiết kế các vi kênh dẫn không van, cơ cấu trộn mẫu và cơ cấu định lượng mẫu tự động không tiếp xúc.
  • Chuyên gia nghiên cứu và phát triển thiết bị chẩn đoán y sinh (Lab-on-a-Chip): Ứng dụng nguyên lý mao dẫn nhiệt laser công suất thấp (40 mW - 53 mW) để thiết kế các bộ xét nghiệm nhanh bệnh truyền nhiễm, máy khuếch đại ADN/ARN và thiết bị sàng lọc tế bào ung thư sớm lưu thông trong máu.
  • Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Cơ khí, Cơ kỹ thuật và Vật lý kỹ thuật: Sử dụng luận văn như tài liệu tham khảo mẫu mực về kỹ thuật mô phỏng phần tử hữu hạn nâng cao, khai thác hiệu quả sự phối hợp giữa phương pháp bảo toàn định mức (CLS), phần tử chuyển động (ALE) và lực căng bề mặt liên tục (CSF) trên phần mềm COMSOL.
  • Kỹ sư phát triển hệ thống quang cơ điện tử (Optomechatronics): Khai thác các dữ liệu tương tác nhiệt quang để tích hợp các module laser bán dẫn kích thước nhỏ vào hệ thống điều khiển tự động hóa các thao tác vi mô trên chip.

Câu hỏi thường gặp

Hiện tượng mao dẫn nhiệt (Marangoni) hoạt động như thế nào trong kênh vi lưu?

Khi nguồn nhiệt laser 40 mW chiếu vào kênh dẫn, nhiệt độ phân bố không đều tạo ra sự chênh lệch sức căng bề mặt dọc theo mặt phân cách giữa nước và dầu. Lưu chất di chuyển từ vùng nhiệt độ cao về vùng nhiệt độ thấp hơn, hình thành một cặp dòng xoáy đối lưu đẩy giọt nước dài 0.55 mm tiến về phía trước một cách tự động.

Tại sao cần kết hợp đồng thời ba phương pháp CLS, ALE và CSF trong mô phỏng số?

Phương pháp CLS đảm bảo bảo toàn khối lượng của giọt chất lỏng qua hàm pha, CSF chuyển hóa lực căng bề mặt thành lực khối tương thích với phương trình Navier-Stokes, và ALE cho phép lưới phần tử biến dạng liên tục theo biên phân pha. Sự phối hợp này giúp kiểm soát chính xác 308,574 bậc tự do mà không làm méo lưới tính toán.

Góc tiếp xúc động (DCA) ảnh hưởng như thế nào đến sự chuyển động của giọt nước?

Chênh lệch giữa góc tiếp xúc phía trước và phía sau sinh ra sự bất đối xứng về bán kính cong của hai đầu giọt nước. Theo phương trình Young-Laplace, sự chênh lệch này tạo ra một độ sụt áp suất mao dẫn dọc theo thân giọt, trực tiếp cộng hưởng cùng lực Marangoni để quyết định gia tốc và vận tốc di chuyển của giọt chất lỏng.

Việc tăng công suất nguồn nhiệt từ 40 mW lên 53 mW mang lại lợi ích gì?

Việc tăng 32.5% công suất nhiệt laser tạo ra gradient nhiệt độ dốc hơn dọc theo chiều dài kênh dẫn 10 mm. Điều này giúp tăng cường độ của cặp xoáy Marangoni, nâng vận tốc cực đại của giọt nước lên thêm khoảng 35%, giúp rút ngắn thời gian vận chuyển mẫu phẩm trong các quy trình xét nghiệm sinh hóa.

Điều kiện biên trượt Navier có vai trò quan trọng ra sao trong mô hình?

Nếu áp dụng điều kiện ranh giới không trượt truyền thống, tại đường tiếp xúc ba pha rắn - lỏng - lỏng sẽ phát sinh điểm kỳ dị toán học với ứng suất cắt tiến tới vô hạn, khiến mô phỏng số bị phân kỳ. Việc thiết lập độ dài trượt hữu hạn 1 nm loại bỏ điểm kỳ dị này và phản ánh chính xác động học trượt ở cấp độ phân tử.

Kết luận

  • Xây dựng thành công mô hình mô phỏng số hai chiều giải quyết toàn diện bài toán dịch chuyển mao dẫn nhiệt của giọt chất lỏng trong vi kênh dưới ảnh hưởng của nguồn nhiệt laser.
  • Chứng minh cơ chế truyền động vi lưu xuất phát từ sự kết hợp tương hỗ giữa cặp xoáy đối lưu nhiệt Marangoni không đối xứng và chênh lệch áp suất mao dẫn do biến thiên góc tiếp xúc động.
  • Thiết lập giải pháp tính toán ổn định và tin cậy cao trên nền tảng phần tử hữu hạn với 20,432 phần tử và 308,574 bậc tự do, kết hợp hoàn hảo các thuật toán CLS, ALE và CSF.
  • Xác lập định lượng tác động của công suất laser, khẳng định nguồn nhiệt 53 mW giúp gia tăng khoảng 35% vận tốc dịch chuyển so với nguồn nhiệt 40 mW, mở ra khả năng điều khiển vi lưu linh hoạt bằng quang nhiệt.
  • Đặt nền tảng lý thuyết và số trị quan trọng phục vụ chế tạo các thiết bị Lab-on-a-Chip và cảm biến y sinh thế hệ mới không cần hệ thống bơm cơ học.

Trong giai đoạn 2024 - 2026, các nhóm nghiên cứu cần nhanh chóng thúc đẩy các thử nghiệm 3D và triển khai thực nghiệm trên vi mạch chip thực tế. Hãy áp dụng ngay các giải pháp thiết kế và bộ thông số mô phỏng từ công trình này để tối ưu hóa hiệu năng cho các sản phẩm vi lưu và thiết bị vi cơ điện tử y sinh của bạn.