Tổng quan luận án

Đề tài "Nghiên cứu chuyển động của hạt lưu chất dưới tác động của mao dẫn nhiệt" do nghiên cứu sinh Nguyễn Tuấn Vinh thực hiện, thuộc chuyên ngành Cơ học chất lỏng và chất khí (mã số: 9440108), bảo vệ 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 vào năm 2022. Đề tài được thực hiện dưới sự hướng dẫn khoa học của TS. Nguyễn Hồng Phan và PGS. TS. Vũ Văn Trường, đồng thời là một phần của đề tài mã số 107.03-2019.307 (107.307) do Quỹ Phát triển khoa học và công nghệ Quốc gia (NAFOSTED) tài trợ.

Tính cấp thiết của luận án xuất phát từ nhu cầu kiểm soát chuyển động của các cấu trúc vi lưu trong công nghệ y sinh và các ứng dụng không gian:

  • Trong kỹ thuật truyền dẫn thuốc hướng đích và điều trị xơ vữa động mạch vành, việc bao bọc dược chất trong cấu trúc lưu chất đa lớp (hạt lưu chất đa thành phần) giúp cải thiện độ hòa tan, kéo dài thời gian lưu thông và tập trung thuốc đúng vị trí bệnh lý. Việc điều khiển chuyển động của các hạt này thông qua chênh lệch nhiệt độ bề mặt mở ra phương thức dẫn thuốc không xâm lấn.
  • Trong môi trường vi trọng lực (không gian), khi lực nổi và trọng lực bị triệt tiêu, lực mao dẫn nhiệt (hiệu ứng Marangoni nhiệt) trở thành lực dẫn hướng chủ đạo, ứng dụng trong việc tách bong bóng khí khi sản xuất vật liệu thủy tinh, kiểm soát hệ thống làm lạnh nhân tạo và bình nhiên liệu tên lửa lỏng.
  • Trong các thiết bị vi lỏng (microfluidics), kích thước hình học ở quy mô micromét khiến lực mao dẫn nhiệt chi phối đáng kể động lực học dòng chảy, cho phép thao tác định tuyến, bẫy hoặc phân loại hạt lưu chất mà không cần hệ thống bơm cơ học phức tạp.

Khoảng trống nghiên cứu được tác giả xác định là việc thiếu hụt các mô hình và công cụ số hóa có khả năng mô tả chính xác tương tác thủy động lực học, biến dạng bề mặt, độ lệch tâm và quá trình phân tách của hạt/sợi lưu chất đa thành phần trong cả không gian mở lẫn môi trường bị giới hạn hình học (kênh thu hẹp) dưới gradient nhiệt độ.

Mục tiêu tổng quát của luận án là nghiên cứu chuyển động của hạt lưu chất đa thành phần dưới tác động của mao dẫn nhiệt bằng phương pháp theo dấu biên trên nền kỹ thuật sai phân hữu hạn.

Nội dung nghiên cứu chính bao gồm:

  1. Nghiên cứu động lực học biến dạng và tách hạt lưu chất trong môi trường đẳng nhiệt.
  2. Nghiên cứu động lực học mao dẫn nhiệt của hạt lưu chất đa thành phần trong không gian mở.
  3. Nghiên cứu động lực học mao dẫn nhiệt của hạt lưu chất đa thành phần khi chuyển động qua ống dẫn thu hẹp.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Sợi lưu chất đa thành phần và hạt lưu chất đa thành phần (dạng nhân - vỏ lỏng) có kích thước từ vài micromét đến vài milimét.
  • Phạm vi nghiên cứu: Tập trung vào phương pháp tính toán mô phỏng số thông qua giải hệ phương trình Navier-Stokes và phương trình bảo toàn năng lượng cho dòng chất lỏng Newton không nén được, giả thiết các pha lưu chất không trộn lẫn, dòng chảy tầng có số Reynolds nhỏ, bài toán được thiết lập ở dạng 2 chiều đối xứng trục trong hệ tọa độ trụ.

Tổng quan tài liệu và vị trí của luận án

Luận án đã tổng quan hệ thống các công trình nghiên cứu trên thế giới và tại Việt Nam theo ba nhánh tiếp cận:

  1. Nghiên cứu về sự hình thành và phân tách của sợi/tia lưu chất:
  • Về sợi đơn lớp: Stone và cộng sự (1986, 1999) đã xác định cơ chế co lại và hiện tượng thắt đứt đầu sợi ("end-pinching"). Schulkes (1994) mở rộng nghiên cứu sợi lưu chất trong môi trường thụ động bằng phương pháp phần tử hữu hạn. Notz và Basaran (2004) thiết lập biểu đồ pha xác định ranh giới giữa chế độ co không phân tách và chế độ phân tách phụ thuộc vào tỷ số hình dạng ban đầu và số Ohnesorge ($Oh$). Castrejón-Pita và cộng sự (2012), Hoepffner và Paré (2013), Anthony và cộng sự (2019) tiếp tục hoàn thiện các cơ chế cuộn xoáy nội tại và sơ đồ chuyển pha hình thái của sợi nhớt.
  • Về tia/sợi đa thành phần: Hertz và Hermanrud (1984) mở đầu việc phân tách tia đa lớp. Utada và cộng sự (2005), Pan và cộng sự (2014) thực nghiệm tạo hạt nhũ tương kép (W/O/W) trong kênh vi lưu. Vũ và cộng sự (2013, 2020), Bhagat và cộng sự (2018), Zhou và cộng sự (2014) mô phỏng quá trình tạo giọt đa nhân bằng các phương pháp trường pha và theo dấu biên.
  1. Nghiên cứu lý thuyết và thực nghiệm về mao dẫn nhiệt:
  • Nghiên cứu lý thuyết: Young, Goldstein và Block (1959) đặt nền móng giải tích cho chuyển động mao dẫn nhiệt của bóng khí đơn trong dòng chảy Stokes. Subramanian (1981, 1983), Balasubramaniam và Chai (1987), Haj-Hariri và cộng sự (1990) phân tích nghiệm tiệm cận và định lý đối ứng Lorentz cho hạt đơn khi có truyền nhiệt đối lưu và ảnh hưởng quán tính. Meyyapan, Subramanian và Acrivos (1981, 1983), Anderson (1985) nghiên cứu tương tác giữa hai bóng khí. Morton và cộng sự (1990), Jadhav và Ghosh (2013) nghiên cứu tác động của mao dẫn nhiệt lên hạt đa thành phần lệch tâm trong dòng chảy Poiseuille.
  • Nghiên cứu thực nghiệm: Selva và cộng sự (2010), Won và cộng sự (2011) thực nghiệm chuyển động mao dẫn nhiệt trong chip vi lỏng và điều khiển hạt bằng xung bức xạ laser.
  1. Phương pháp mô phỏng số cho dòng mao dẫn nhiệt:
  • Phương pháp Level-Set (LS): Haj-Hariri và cộng sự (1997), Zhao và cộng sự (2011), Brady và cộng sự (2011), Balcára và cộng sự (2013) giải phương trình Navier-Stokes bắt biên ẩn, khảo sát ảnh hưởng của số Marangoni ($Ma$) và số mao dẫn ($Ca$).
  • Phương pháp thể tích lưu chất (VOF): Samareh và cộng sự (2011), Tripathi và Sahu (2018), Capobianchi và cộng sự (2018) mô phỏng bóng khí chịu đồng thời lực nổi và lực mao dẫn nhiệt trong không gian 3 chiều.
  • Phương pháp mạng Boltzmann (LBM): Liu và cộng sự (2011, 2013), Xie và cộng sự (2014), Mitchell và cộng sự (2018) phân tích quá trình hợp nhất giọt và vi kênh nhiệt.
  • Phương pháp theo dấu biên (Front-Tracking - FT): Tryggvason và cộng sự (2001), Nas và Tryggvason (1993, 2003), Yin và cộng sự (2008, 2011), Lu và cộng sự (2009), Kalichetty và cộng sự (2019) giải hệ phương trình nhiều pha, khảo sát tương tác đa hạt và ảnh hưởng của vật cản.
  • Tại Việt Nam: Nhóm nghiên cứu của PGS. TS. Vũ Văn Trường tại Trường Đại học Phenikaa và Viện Cơ học đã công bố nhiều công trình về động lực học hạt đa lớp trong kênh hẹp và dòng cắt (Nguyễn và cs., 2017, 2020; Vũ và cs., 2019).

Vị trí của luận án: Luận án kế thừa kỹ thuật theo dấu biên của Tryggvason nhưng tập trung phát triển sâu vào bài toán chuyển động mao dẫn nhiệt của hạt lưu chất đa thành phần (compound droplet) hai lớp chất lỏng, khảo sát đồng thời cả hai hiện tượng: quá trình co rút phân tách đẳng nhiệt của sợi đa lớp và động lực học chuyển động - biến dạng của hạt đa lớp trong kênh dẫn thu hẹp hình sin.

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

Luận án sử dụng hệ phương trình cơ bản của cơ học chất lỏng liên tục mô tả dòng chảy Newton không nén được, bao gồm:

  • Phương trình bảo toàn khối lượng (phương trình liên tục): $$\nabla \cdot \mathbf{u} = 0$$
  • Phương trình bảo toàn động lượng (Navier-Stokes cho dòng nhiều pha có bổ sung lực căng bề mặt): $$\frac{\partial (\rho \mathbf{u})}{\partial t} + \nabla \cdot (\rho \mathbf{u} \mathbf{u}) = -\nabla p + \nabla \cdot \left[ \mu \left( \nabla \mathbf{u} + \nabla \mathbf{u}^T \right) \right] + \rho \mathbf{g} + \mathbf{F}_\sigma$$
  • Phương trình bảo toàn năng lượng nhiệt: $$\frac{\partial (\rho c_p T)}{\partial t} + \nabla \cdot (\rho c_p \mathbf{u} T) = \nabla \cdot (k \nabla T)$$
  • Phương trình trạng thái của sức căng bề mặt theo nhiệt độ: $$\sigma(T) = \sigma_0 - \gamma_T (T - T_0)$$ trong đó $\gamma_T = -\frac{\partial \sigma}{\partial T}$ là hệ số suy giảm sức căng bề mặt theo nhiệt độ.

Lực căng bề mặt trên mặt phân cách được biểu diễn thông qua tích phân mặt kết hợp với hàm Dirac delta: $$\mathbf{F}_\sigma = \int_f \left( \sigma \kappa \mathbf{n}_f + \nabla_s \sigma \right) \delta(\mathbf{x} - \mathbf{x}_f) , dS$$

Phương pháp nghiên cứu thực tế và công cụ tính toán:

  • Phương pháp theo dấu biên (Front-Tracking Method - FT): Phối hợp hai hệ thống lưới gồm lưới Euler cố định dạng so le (Staggered MAC Grid) để giải các trường vận tốc, áp suất, nhiệt độ và lưới Lagrange di động (tập hợp các điểm nút liên kết phi cấu trúc 1D) để nắm bắt mặt phân cách giữa các pha chất lỏng.
  • Kỹ thuật số: Sử dụng phương pháp sai phân hữu hạn (FDM) bậc hai trong không gian; tích phân thời gian bằng lược đồ dự đoán - hiệu chỉnh (predictor-corrector); phương pháp lặp SOR (Successive Overrelaxation) để giải phương trình Poisson áp suất.
  • Thuật toán tái cấu trúc biên: Tự động thêm hoặc bớt các điểm Lagrange để duy trì khoảng cách phần tử biên đồng đều khi giao diện bị biến dạng hoặc kéo dãn mạnh.
  • Công cụ triển khai: Chương trình tính toán được lập trình hoàn toàn bằng ngôn ngữ Fortran.

Hệ thống tham số không thứ nguyên được sử dụng trong mô hình:

Tham số không thứ nguyên Công thức xác định Ý nghĩa vật lý
Số Reynolds ($Re$) $Re = \frac{\rho_o U_c L}{\mu_o}$ Tỷ số giữa lực quán tính và lực nhớt
Số Marangoni ($Ma$) $Ma = \frac{\rho_o c_{p,o} U_c L}{k_o}$ Tỷ số giữa vận chuyển đối lưu nhiệt Marangoni và khuếch tán nhiệt
Số mao dẫn ($Ca$) $Ca = \frac{\mu_o U_c}{\sigma_o}$ Tỷ số giữa lực cản nhớt và lực căng bề mặt
Số Ohnesorge ($Oh$) $Oh = \frac{\mu_o}{\sqrt{\rho_o \sigma_o R_o}}$ Tương quan lực nhớt với lực quán tính và sức căng bề mặt
Số Weber ($We$) $We = \frac{\rho_o L U_c^2}{\sigma_o}$ Tỷ số giữa lực quán tính và lực căng bề mặt
Tỷ số bán kính ($R_{io}$) $R_{io} = \frac{R_i}{R_o}$ Tỷ số kích thước giữa lõi trong và hạt ngoài
Tỷ số độ nhớt ($\mu_{mo}, \mu_{io}$) $\mu_{mo} = \frac{\mu_m}{\mu_o}, \mu_{io} = \frac{\mu_i}{\mu_o}$ Tỷ số độ nhớt động lực học giữa các pha lưu chất
Tỷ số co sợi ($Ar_o, Ar_i$) $Ar_o = \frac{L_o}{2R_o}, Ar_i = \frac{L_i}{2R_i}$ Tỷ số giữa nửa chiều dài và bán kính sợi ban đầu

Bảng so sánh phương pháp số (tổng hợp từ luận án):

Tiêu chí Front-Tracking (FT) Level-Set (LS) Volume of Fluid (VOF)
Bản chất phương pháp Dịch chuyển điểm đánh dấu Lagrangian Dựa vào hàm khoảng cách ẩn Dựa vào phương trình thể tích phân đoạn
Độ chính xác giao diện Rất cao (xác định tường minh biên) Cao Trung bình (phụ thuộc tái tạo biên)
Tái cấu trúc biên Cần thuật toán thêm/bớt điểm Tự động cập nhật qua hàm mức Tự động nhưng biên dạng tuyến tính đơn giản
Tính lực căng bề mặt Rất chính xác và trực tiếp Dễ dàng Khó khăn, dễ sinh vận tốc giả
Bảo toàn thể tích Bảo toàn tốt Kém hơn, dễ bị tiêu tán Bảo toàn khối lượng tuyệt đối

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 hệ thống hóa các định nghĩa và khái niệm cơ bản về lưu chất, sợi lưu chất đa thành phần, hạt lưu chất đa thành phần (compound drop), hiện tượng mao dẫn thông thường và hiện tượng mao dẫn nhiệt (dòng Marangoni). Khi tồn tại gradient nhiệt độ dọc theo bề mặt hạt lưu chất, sức căng bề mặt thay đổi không đều dẫn đến xuất hiện ứng suất cắt tiếp tuyến kéo lưu chất xung quanh chuyển động ngược chiều gradient nhiệt, đẩy hạt lưu chất di chuyển về phía vùng có nhiệt độ cao hơn (vùng lạnh sang vùng nóng).

Chương này trình bày chi tiết bức tranh tổng quan về các phương thức hình thành hạt lưu chất trong tự nhiên và công nghệ vi lỏng, phân tích các nghiên cứu lý thuyết, thực nghiệm và mô phỏng số trên thế giới (sử dụng các phương pháp LS, VOF, LBM, FT) và tình hình nghiên cứu tại Việt Nam. Tác giả chỉ ra tính ưu việt của phương pháp FT trong việc tính toán sức căng bề mặt biến thiên theo nhiệt độ mà không bị ảnh hưởng bởi lỗi làm nhòe biên.

Chương 2: Xây dựng chương trình mô phỏng chuyển động của hạt lưu chất dưới tác động của mao dẫn nhiệt

Chương 2 trình bày chi tiết thuật toán xây dựng chương trình số hóa:

  • Thiết lập hệ phương trình Navier-Stokes và phương trình năng lượng ở dạng biến số thực và dạng không thứ nguyên trên lưới so le MAC.
  • Quy trình tích phân thời gian bằng lược đồ dự đoán - hiệu chỉnh: Vận tốc dự đoán $\mathbf{u}^*$ được tính từ số hạng đối lưu, khuếch tán và lực căng bề mặt; sau đó hiệu chỉnh áp suất qua phương trình Poisson giải bằng phương pháp lặp SOR.
  • Phương pháp theo dấu biên: Cấu trúc biên 1D gồm các đỉnh nối nhau bằng các đoạn thẳng; kỹ thuật nội suy vận tốc từ lưới Euler sang các đỉnh Lagrange; kỹ thuật phân phối lực căng bề mặt Lagrange lên lưới Euler thông qua hàm delta rời rạc của Peskin.
  • Thuật toán tái cấu trúc biên (Re-meshing): Tự động chèn điểm nút mới khi khoảng cách giữa hai đỉnh vượt quá $1.5 \Delta x$ và gộp điểm nút khi khoảng cách nhỏ hơn $0.5 \Delta x$.

Kiểm chứng chương trình mô phỏng:

  • Bài toán 1: Sự co lại của sợi lưu chất đơn lớp không biên. Kết quả mô phỏng hình thái biên dạng co rút ở các mốc thời gian với $Oh = 1.0, Oh = 0.01$ và $Ar_o = 15$ hoàn toàn trùng khớp với kết quả giải tích và số của Notz và Basaran (2004).
  • Bài toán 2: Chuyển động mao dẫn nhiệt của hạt lưu chất đơn trong gradient nhiệt độ đồng nhất. Vận tốc ổn định tiệm cận chính xác với nghiệm giải tích của Young, Goldstein và Block (1959).

Chương 3: Nghiên cứu động lực học biến dạng và tách hạt lưu chất trong môi trường đẳng nhiệt

Chương 3 tập trung giải quyết câu hỏi về cơ chế hình thành hạt lưu chất đa thành phần thông qua nghiên cứu động lực học co lại của sợi lưu chất hai lớp (sợi trong bán kính $R_i$ nằm đồng trục trong sợi ngoài bán kính $R_o$) trong môi trường chất lỏng vô hạn khi chưa có tác động của nhiệt độ.

Các kết quả và quy luật chính:

  • Kiểm tra tính hội tụ của lưới: Thực hiện trên ba mật độ lưới khác nhau ($64 \times 128$, $128 \times 256$, $256 \times 512$), khẳng định tính độc lập của lưới đối với tham số biến dạng sợi trong $T_i$ và sợi ngoài $T_o$.
  • Nhận diện các chế độ biến dạng:
    1. Chế độ không phân tách (Non-breakup): Sợi co rút êm dịu thành một hạt lưu chất đa thành phần hình cầu duy nhất.
    2. Chế độ phân tách sợi bên trong (Inner pinch-off, ký hiệu "i-"): Lớp sợi trong bị thắt đứt thành nhiều hạt nhân nhỏ trong khi vỏ ngoài vẫn nguyên vẹn.
    3. Chế độ phân tách sợi bên ngoài (Outer pinch-off, ký hiệu "o-"): Vỏ ngoài bị thắt đứt trước ở đầu sợi hoặc thân sợi.
    4. Chế độ phân tách hỗn hợp (Mixed pinch-off, ký hiệu "i-/o-"): Cả sợi trong và vỏ ngoài đồng thời bị đứt gãy phức tạp, sinh ra các hạt vệ tinh (satellite droplets).
  • Ảnh hưởng của các thông số vật lý và hình học:
    • Tăng số Ohnesorge ($Oh$ từ 0.01 lên 1.0): Lực nhớt tăng làm triệt tiêu dao động mao dẫn, chuyển trạng thái từ phân tách mạnh sang co êm dịu không phân tách.
    • Tỷ số hình dạng ngoài $Ar_o$ và trong $Ar_i$: Sợi càng dài ($Ar_o \ge 15$) càng thúc đẩy hiện tượng tách đầu sợi ("end-pinching").
    • Tỷ số bán kính $R_{io}$: Khi $R_{io}$ tăng (lõi trong dày lên), tương tác giữa hai mặt phân cách tăng mạnh, thúc đẩy sự phân tách hỗn hợp.
    • Tỷ số sức căng bề mặt $\sigma_{io}$ và tỷ số độ nhớt ($\mu_{32}, \mu_{12}$): Chi phối trực tiếp tần số dao động co rút và thời gian phân tách.
  • Chương này đã xây dựng hệ thống sơ đồ pha (phase diagrams) biểu diễn ranh giới phân tách giữa $Ar_o$ với $Oh, R_{io}, Ar_i, \sigma_{io}, \mu_{32}, \mu_{12}$.

Chương 4: Nghiên cứu động lực học mao dẫn nhiệt của hạt lưu chất đa thành phần

Chương 4 khảo sát sự dịch chuyển của hạt lưu chất đa thành phần (gồm một hạt nhân bên trong và một lớp vỏ bên ngoài) đặt trong môi trường chất lỏng có gradient nhiệt độ không đổi dọc theo trục đối xứng $z$.

Các kết quả định lượng chính:

  • Khảo sát hội tụ lưới: Kiểm chứng với các kích thước lưới $192 \times 384$ và $256 \times 512$, sai lệch vị trí tâm và vận tốc đạt dưới 1%.
  • Động lực học dịch chuyển: Do chênh lệch sức căng bề mặt, hạt di chuyển từ cực lạnh sang cực nóng. Sau giai đoạn gia tốc ban đầu, vận tốc hạt ngoài $U_o$ và hạt trong $U_i$ đạt trạng thái ổn định và chuyển động cùng nhau.
  • Ảnh hưởng của số Reynolds ($Re = 1, 10, 30, 50$): Khi $Re$ tăng, lực quán tính tăng làm vận tốc dịch chuyển cực đại tăng lên, hạt dịch chuyển nhanh hơn dọc trục; độ lệch tâm giữa hạt trong và hạt ngoài biến thiên rõ nét hơn ở giai đoạn đầu trước khi ổn định.
  • Ảnh hưởng của số Marangoni ($Ma = 20, 60, 100$): Tăng $Ma$ đồng nghĩa với việc truyền nhiệt đối lưu chiếm ưu thế hơn khuếch tán nhiệt. Hiệu ứng quét nhiệt làm biến dạng các đường đẳng nhiệt xung quanh hạt, làm giảm gradient nhiệt độ cục bộ tại bề mặt, dẫn đến vận tốc dịch chuyển cực đại của hạt giảm nhẹ và độ lệch tâm cực đại giảm.
  • Ảnh hưởng của tỷ số bán kính ($R_{io} = 0.4, 0.5, 0.6$): Khi tăng kích thước hạt trong ($R_{io}$ lớn), lực cản thủy động tăng lên làm vận tốc dịch chuyển ổn định của toàn bộ cấu trúc giảm nhẹ.
  • Ảnh hưởng của tỷ số độ nhớt ($\mu_{mo} = 0.1, 0.4, 1.6$): Tăng độ nhớt của lớp lưu chất giữa làm gia tăng lực ma sát nội tại, làm giảm tốc độ dịch chuyển của hạt đa thành phần.

Chương 5: Nghiên cứu động lực học và biến dạng của hạt lưu chất đa thành phần trong ống dẫn thu hẹp dưới tác động của mao dẫn nhiệt

Chương 5 khảo sát hạt lưu chất đa thành phần chuyển động qua ống dẫn đối xứng trục có đoạn thu hẹp cục bộ dạng hình sin $r(z) = R_c - \frac{d}{2} \left[ 1 + \cos\left(\frac{\pi (z - z_c)}{L_c}\right) \right]$.

Các thông số cơ sở được cố định: $Re = 10$, $\rho_{io} = \mu_{io} = k_{io} = C_{io} = 1$, $\rho_{mo} = \mu_{mo} = k_{mo} = C_{mo} = 0.5$, $\sigma_{io} = 1$, $R_c / R_o = 2.5$. Mật độ lưới kiểm chứng đạt độ mịn $512 \times 2048$.

Các phát hiện chính:

  • Cơ chế vượt cổ thắt: Khi đến gần cổ thắt (constriction), hạt lưu chất chịu lực cản hình học lớn từ thành ống, vận tốc giảm mạnh. Vận tốc hạt biểu hiện ba đỉnh đặc trưng tương ứng với các giai đoạn: tăng tốc ban đầu ở thượng lưu, lọt vào khe hẹp và bứt tốc ở hạ lưu sau khi thoát khỏi cổ ống.
  • Ảnh hưởng của số Marangoni ($Ma = 5, 20, 80$):
    • Với $Ma = 5$ (khuếch tán nhiệt chi phối), lực quán tính mạnh hơn lực mao dẫn nhiệt khiến vận tốc cực đại xuất hiện ở đỉnh thứ ba (hạ lưu).
    • Với $Ma > 20$ (đối lưu nhiệt chi phối), vận tốc cực đại chuyển về đỉnh thứ nhất. Tăng $Ma$ làm giảm vận tốc cực đại và kéo dài thời gian đến ($t_{art}$) của hạt. Độ lệch tâm cực đại giảm khi $Ma$ tăng.
  • Ảnh hưởng của số mao dẫn ($Ca = 0.01, 0.02, 0.04, 0.05, 0.08$):
    • Khi $Ca \le 0.04$, sức căng bề mặt chi phối giúp hạt duy trì hình dạng cầu, tạo lực cản lớn khi vào khe hẹp, vận tốc tại cổ ống gần như bằng không.
    • Khi $Ca$ tăng lên 0.08, lực nhớt thắng sức căng bề mặt khiến hạt bị kéo dãn thuôn dài linh hoạt, dễ dàng trượt qua cổ thắt, rút ngắn đáng kể tổng thời gian chuyển tiếp ($t_{art}$) mặc dù độ lệch tâm cực đại tăng lên.
  • Ảnh hưởng của kích thước lõi trong ($R_{io} = 0.3 - 0.8$):
    • Độ lệch tâm cực đại thường đạt được khi hạt chuẩn bị rời cổ ống. Tuy nhiên, khi $R_{io} = 0.8$, độ lệch tâm cực đại chuyển dịch lên vùng thượng lưu ngay trước khi vào cổ ống do hạt trong bị nén chặt.
    • Thời gian đến ($t_{art}$) không biến thiên đơn điệu: $t_{art}$ tăng khi $R_{io}$ tăng từ 0.3 đến 0.6, nhưng sau đó giảm khi $R_{io}$ tăng từ 0.7 lên 0.8 do hiệu ứng tăng ứng suất kéo dãn mép trước của hạt.
  • Ảnh hưởng của độ thắt ống dẫn ($d/R_c = 0, 0.4, 0.6, 0.8$):
    • Trường hợp ống thẳng ($d/R_c = 0$), hạt chuyển động với vận tốc cao nhất và không bị cản trở.
    • Khi $d/R_c$ tăng từ 0 đến 0.6, lực cản tăng làm hạt bị kéo dãn mạnh, độ lệch tâm tăng và thời gian qua ống kéo dài.
    • Khi $d/R_c \ge 0.8$, độ thắt quá hẹp khiến lực cản hình học và áp suất ngược vượt quá lực kéo mao dẫn nhiệt; hạt lưu chất bị tắc nghẽn hoàn toàn và không thể đi qua cổ ống.

Kết quả và những đóng góp mới

Các đóng góp mới về mặt khoa học và kỹ thuật của luận án:

  1. Về phương pháp tính toán: Đã xây dựng và hoàn thiện bộ chương trình mô phỏng số tự viết bằng ngôn ngữ Fortran dựa trên phương pháp theo dấu biên (Front-Tracking) kết hợp sai phân hữu hạn, giải quyết ổn định bài toán tương tác nhiệt - thủy động lực học của hệ lưu chất đa pha, đa thành phần có bề mặt tự do chịu tác động của lực mao dẫn nhiệt.
  2. Về cơ chế phân tách sợi đa lớp: Đã phân tích chi tiết cơ chế động lực học co rút của sợi lưu chất đa thành phần trong môi trường đẳng nhiệt; phát hiện và phân loại 4 chế độ phân tách (không phân tách, tách sợi trong, tách sợi ngoài, tách hỗn hợp); thiết lập hệ thống biểu đồ pha làm cơ sở cho việc kiểm soát hình thái tạo hạt trong công nghệ chế tạo vi hạt.
  3. Về động lực học mao dẫn nhiệt trong không gian mở: Xác định rõ quy luật định lượng của trường nhiệt độ, trường vận tốc, độ lệch tâm và tốc độ dịch chuyển của hạt lưu chất đa thành phần theo các thông số $Re, Ma, R_{io}, \mu_{mo}$, làm rõ tương tác thủy động giữa lõi trong và vỏ ngoài.
  4. Về động lực học hạt trong kênh thu hẹp: Làm sáng tỏ hành vi biến dạng phức tạp, sự dịch chuyển tâm và quy luật thay đổi thời gian hành trình ($t_{art}$) của hạt đa thành phần qua ống thắt hình sin dưới tác động của các tham số $Ma, Ca, R_{io}, d/R_c$; chỉ ra giới hạn hình học tới hạn ($d/R_c = 0.8$) dẫn đến hiện tượng tắc nghẽn dòng vi lưu.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án:

  • Mô hình tính toán mới chỉ giới hạn ở cấu trúc hạt đa thành phần dạng 1 lõi - 1 vỏ (đơn nhân đồng trục); chưa mô phỏng cấu trúc phức tạp chứa nhiều hạt nhân phân tán.
  • Không gian tính toán dừng lại ở mô hình 2 chiều đối xứng trục trong hệ tọa độ trụ, chưa phản ánh đầy đủ các chuyển động xoáy và biến dạng bất đối xứng trong không gian 3 chiều thực tế.
  • Các kết quả nghiên cứu hoàn toàn dựa trên mô phỏng số, chưa có điều kiện đối chứng trực tiếp với các thí nghiệm vi lưu hoặc thí nghiệm vi trọng lực tương ứng.

Hướng phát triển tiếp theo được tác giả đề xuất:

  1. Mở rộng mô hình cho sợi và hạt lưu chất chứa 2, 3 hoặc nhiều nhân phân bố ngẫu nhiên bên trong dưới tác động của mao dẫn nhiệt.
  2. Phát triển mã nguồn mô phỏng 3 chiều (3D) đầy đủ để nắm bắt các biến dạng không đối xứng trục và tương tác biên phức tạp.
  3. Thiết kế và thực hiện các thí nghiệm thực tế (trên chip vi lưu hoặc môi trường mô phỏng vi trọng lực) nhằm kiểm chứng và hiện thực hóa các ứng dụng trong y sinh học và chế tạo vật liệu không gian.

Giá trị tham khảo

Luận án là tài liệu tham khảo có giá trị chuyên môn cao cho các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Cơ học chất lỏng, Thủy khí động lực học, Kỹ thuật nhiệt và Toán ứng dụng: Cung cấp phương pháp luận và kỹ thuật chi tiết để xây dựng thuật toán Front-Tracking giải bài toán biên di động đa pha.
  • Các kỹ sư và nhà nghiên cứu trong lĩnh vực Công nghệ vi lỏng (Microfluidics) và Kỹ thuật y sinh: Tham khảo số liệu động lực học, độ biến dạng và ngưỡng hình học tới hạn để thiết kế các kênh vi lưu định tuyến thuốc, điều khiển giọt bằng laser hoặc nhiệt độ.
  • Các nhà khoa học công nghệ vật liệu: Ứng dụng các biểu đồ pha co rút sợi lỏng để kiểm soát kích thước hạt nhũ tương kép, phục vụ chế tạo vật liệu mới (hạt cầu rỗng kim loại, composite, kính quang học) trong môi trường vi trọng lực.

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

Hiện tượng mao dẫn nhiệt làm hạt lưu chất dịch chuyển theo cơ chế nào?

Khi có một gradient nhiệt độ đặt lên môi trường lưu chất, nhiệt độ dọc theo bề mặt ngoài của hạt lưu chất sẽ không đồng đều. Do sức căng bề mặt của hầu hết các chất lỏng giảm khi nhiệt độ tăng ($\gamma_T > 0$), vùng bề mặt tiếp xúc với nhiệt độ thấp sẽ có sức căng bề mặt lớn hơn vùng tiếp xúc với nhiệt độ cao. Chênh lệch sức căng bề mặt sinh ra ứng suất cắt tiếp tuyến (lực Marangoni), kéo lớp chất lỏng sát bề mặt chuyển động từ vùng ấm về vùng lạnh, phản lực thủy động này đẩy toàn bộ hạt lưu chất dịch chuyển về phía vùng có nhiệt độ cao hơn (từ lạnh sang nóng).

Phương pháp theo dấu biên (Front-Tracking) có ưu điểm gì so với Level-Set và VOF trong bài toán này?

Phương pháp Front-Tracking sử dụng một lưới Lagrange rời rạc rõ ràng để theo dõi chính xác vị trí của mặt phân cách. Nhờ đó, phương pháp này tính toán trực tiếp và chính xác độ cong bề mặt cũng như gradient tiếp tuyến của sức căng bề mặt khi có biến thiên nhiệt độ ($\nabla_s \sigma$), triệt tiêu các vận tốc giả (parasitic currents) và tránh được hiện tượng làm nhòe biên hoặc mất thể tích vốn là nhược điểm thường gặp ở phương pháp Level-Set và Volume of Fluid.

Vai trò của số Marangoni ($Ma$) đối với chuyển động của hạt qua ống thu hẹp là gì?

Số Marangoni thể hiện tương quan giữa vận chuyển nhiệt bằng đối lưu và khuếch tán. Khi $Ma$ nhỏ ($Ma = 5$), khuếch tán nhiệt chiếm ưu thế, trường nhiệt độ ít bị biến dạng bởi dòng chảy, vận tốc cực đại của hạt xuất hiện ở hạ lưu (đỉnh thứ ba). Khi $Ma$ lớn ($Ma = 80$), đối lưu nhiệt chiếm ưu thế, dòng chảy cuốn các đường đẳng nhiệt làm giảm gradient nhiệt độ cục bộ tại bề mặt hạt, dẫn đến giảm vận tốc chuyển dịch và kéo dài thời gian để hạt đi qua ống thu hẹp ($t_{art}$).

Vì sao số mao dẫn ($Ca$) lớn lại giúp hạt lưu chất vượt qua cổ thắt nhanh hơn?

Số mao dẫn $Ca$ biểu thị tỷ số giữa lực nhớt và lực căng bề mặt. Khi $Ca$ nhỏ ($Ca \le 0.04$), sức căng bề mặt lớn giữ cho hạt luôn có dạng hình cầu cứng, gây tắc nghẽn và lực cản rất lớn khi tiến vào cổ ống hẹp. Khi $Ca$ tăng ($Ca = 0.08$), lực căng bề mặt giảm tương đối so với lực nhớt, hạt dễ dàng bị kéo dãn và biến dạng thuôn dài để trượt qua khe hẹp, làm giảm đáng kể thời gian di chuyển qua cổ thắt.

Hiện tượng gì xảy ra khi tỷ số độ thắt ống dẫn $d/R_c \ge 0.8$?

Khi độ thắt $d/R_c \ge 0.8$, tiết diện của cổ ống bị thu hẹp quá mức so với kích thước ban đầu của hạt lưu chất. Lực cản thủy động và gradient áp suất ngược tại cổ thắt vượt qua lực kéo sinh ra từ hiệu ứng mao dẫn nhiệt. Kết quả mô phỏng cho thấy hạt lưu chất bị chặn lại ở vùng thượng lưu và không thể đi qua cổ ống.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Tuấn Vinh đã phát triển thành công công cụ mô phỏng số dựa trên phương pháp theo dấu biên (Front-Tracking) trên nền tảng Fortran để khảo sát toàn diện động lực học sợi và hạt lưu chất đa thành phần. Công trình đã làm sáng tỏ cơ chế co rút, phân tách đẳng nhiệt của sợi đa lớp cũng như quy luật chuyển động, biến dạng và độ lệch tâm của hạt đa lớp dưới tác động của mao dẫn nhiệt trong không gian mở và qua ống dẫn thu hẹp. Kết quả nghiên cứu cung cấp cơ sở dữ liệu khoa học chính xác và các biểu đồ tham số có ý nghĩa thực tiễn cao cho việc tối ưu hóa thiết bị vi lỏng, công nghệ truyền dẫn thuốc hướng đích và xử lý vật liệu trong môi trường vi trọng lực.