Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Trong công nghệ y sinh, kỹ thuật truyền dẫn thuốc giữ vai trò thiết yếu nhằm nâng cao hiệu lực điều trị. Sự phát triển của công nghệ vi lỏng cho phép tạo ra các hạt cấu trúc kích thước từ vài micromét đến vài milimét chứa dược chất, hỗ trợ đưa thuốc đến chính xác vị trí bệnh lý (như điều trị xơ vữa động mạch vành) với liều lượng và thời gian được kiểm soát. Một cơ chế điều khiển dòng chảy không cần tiếp xúc cơ học là ứng dụng hiện tượng mao dẫn nhiệt (sự biến thiên của sức căng bề mặt theo gradient nhiệt độ) trên các hạt lưu chất đa thành phần (dạng thuốc lưu chất được bao bọc bởi màng lưu chất bên ngoài).

Trong môi trường không gian (vi trọng lực), lực nổi do trọng trường triệt tiêu, lực mao dẫn nhiệt trở thành yếu tố chi phối, được ứng dụng để tách bọt khí trong sản xuất vật liệu (như thủy tinh), kiểm soát bọt khí trong hệ thống nhiên liệu tên lửa lỏng và hệ thống làm lạnh nhân tạo. Tương tự, trong các thiết bị vi lỏng có kích thước hình học nhỏ, vai trò của trọng lực có thể bỏ qua và lực mao dẫn nhiệt là công cụ định hướng chuyển động của hạt lưu chất.

Khoảng trống nghiên cứu được luận án xác định bao gồm:

  1. Phần lớn các công trình trên thế giới chỉ tập trung vào hạt lưu chất một thành phần (hạt đơn) hoặc tương tác giữa các hạt đơn.
  2. Các nghiên cứu lý thuyết ban đầu về hạt lưu chất đa thành phần thường bỏ qua sự khác biệt về khối lượng riêng giữa các pha lưu chất và chưa xem xét tương tác động lực học chi tiết giữa nhân bên trong và màng bao bên ngoài.
  3. Về mặt thực nghiệm và mô phỏng số, còn thiếu các phân tích định lượng toàn diện về 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, đặc biệt là trong các kênh dẫn có hình học thu hẹp phức tạp (như kênh hình sin).
  4. Tại Việt Nam, các nghiên cứu mao dẫn nhiệt trước đây chỉ giới hạn ở hạt lưu chất đơn pha; chưa có công trình nào mô phỏng hoặc thực nghiệm trên hạt lưu chất đa thành phần.

Mục tiêu và nội dung nghiên cứu

  • Mục tiêu nghiên cứu: 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 (Front-Tracking) trên nền kỹ thuật sai phân hữu hạn.
  • Nội dung nghiên cứu:
    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.
    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 trong ố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 có kích thước hình học từ vài micromét ($\mu\text{m}$) đến vài milimét ($\text{mm}$).
  • Phạm vi nghiên cứu:
    • Mô phỏng số thông qua việc 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ảy lưu chất Newton, không nén được và không hòa trộn giữa các pha.
    • Chế độ dòng chảy giới hạn trong miền chảy tầng với số Reynolds nhỏ.
    • Cấu hình hình học khảo sát gồm môi trường không biên/kênh phẳng và ống dẫn thu hẹp đối xứng trục dạng hình sin.

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

Các hướng nghiên cứu trước đây

Tác giả đã hệ thống hóa các công trình nghiên cứu trong và ngoài nước theo các nhóm nội dung:

  1. Sự hình thành hạt và sự co lại của sợi/tia lưu chất:

    • Sợi đơn lớp: Stone và cộng sự (1986) thực nghiệm sự co lại của sợi lưu chất do sức căng bề mặt và phát hiện cơ chế tách đầu sợi ("end-pinching"). Schulkes (1996) dùng phương pháp phần tử hữu hạn mở rộng nghiên cứu này. Notz và Basaran (2004) thiết lập biểu đồ pha dựa trên tỷ số hình dạng ban đầu và số Ohnesorge ($Oh$). Castrejon-Pita và cộng sự (2012), Hoepffner và Paré (2013), Anthony và cộng sự (2019) khảo sát cấu trúc xoáy nội tại và điều kiện phân tách. Các nghiên cứu liên quan khác gồm Cuellar và cộng sự, Contò và cộng sự, Driessen và cộng sự, Dziedzic và cộng sự, Pierson và cộng sự, Wang và cộng sự, Yang và cộng sự.
    • Tia và sợi đa thành phần: Hertz và Hermanrud (1983) nghiên cứu sự phân tách tia đa thành phần. Evangelio và cộng sự (2006), Yu và cộng sự (2009), Utada và cộng sự (2005), Pan và cộng sự (2014) thực nghiệm tạo hạt đa thành phần nước-trong-dầu-trong-nước. Vũ và cộng sự (2013), Bhagat và cộng sự nghiên cứu tia đa thành phần có truyền nhiệt và đông đặc. Các mô phỏng số về tia đa thành phần được thực hiện bởi Zhou và cộng sự (phương pháp trường pha), Vũ và cộng sự (phương pháp theo dấu biên), Liu và cộng sự, Wang và cộng sự, Zhang và cộng sự (phương pháp VOF).

  2. Chuyển động mao dẫn nhiệt của hạt lưu chất trên thế giới:

    • Lý thuyết: Young và cộng sự (1959) đặt nền móng giải tích cho dòng Stokes. Subramanian (1981, 1983) mở rộng cho đối lưu nhiệt với số Marangoni ($Ma$) nhỏ. Balasubramaniam và Chai (1987) tìm nghiệm chính xác cho vận tốc hạt đơn khi đối lưu không đáng kể. Haj-Hariri và cộng sự (1997) dùng định lý đối ứng Lorentz khảo sát ảnh hưởng quán tính. Meyyapan và cộng sự (1981, 1983), Anderson (1985), Acrivos và cộng sự (1990) phân tích chuyển động tương tác giữa hai hoặc nhiều bọt khí. Morton và cộng sự nghiên cứu hạt đa thành phần trong môi trường đồng nhất nhiệt. Jadhav và Ghosh nghiên cứu độ biến dạng và độ lệch tâm của hạt đa thành phần trong dòng Poiseuille.
    • Thực nghiệm: Young và cộng sự (1959) đo vận tốc bọt khí giữa hai đe kim loại có kiểm soát nhiệt độ. Selva và cộng sự (2010) điều khiển giọt trong kênh vi lỏng phân nhánh bằng gradient nhiệt. Won và cộng sự (2020) nghiên cứu lực mao dẫn nhiệt từ nguồn nhiệt laser tác động lên chip vi lỏng.
    • Mô phỏng số:
      • Phương pháp Level-Set (LS): Haj-Hariri và cộng sự, Zhao và cộng sự, Brady và cộng sự, Balcara và cộng sự.
      • Phương pháp thể tích lưu chất (VOF): Samareh và cộng sự, Tripathi và Sahu (kết hợp mao dẫn nhiệt và lực nổi), Capobianchi và cộng sự (kết hợp VOF-LS khảo sát hiệu ứng thành kênh và tắc nghẽn).
      • Phương pháp Lattice Boltzmann (LB): Liu và cộng sự (mô phỏng đốt nóng laser và tương tác hai hạt), Xie và cộng sự, Mitchell và cộng sự.
      • Phương pháp theo dấu biên (FT): Tryggvason và cộng sự (1992, 2001), Nas và Tryggvason (1993, 2003), Yin và cộng sự, Lu và cộng sự (2017 - khảo sát kênh có vật cản), Kalichetty và cộng sự (2019 - sử dụng OpenFOAM).

  3. Tình hình nghiên cứu tại Việt Nam:

    • Nhóm Cơ học lưu chất ứng dụng (Trường Đại học Phenikaa) do PGS. Vũ Văn Trường đứng đầu: Nguyễn và cộng sự khảo sát hạt đa thành phần qua ống hẹp; Vũ và cộng sự khảo sát sự phân tách nhân trong qua kênh hình nón; Nguyễn và cộng sự nghiên cứu kênh hình học phức tạp; các công trình khác nghiên cứu va chạm, tương tác động lực học và biến dạng trong dòng chảy cắt.
    • Trường Đại học Bách khoa TP. Hồ Chí Minh: Lê, Nguyễn và Lê sử dụng phần tử hữu hạn kết hợp Level-Set mô phỏng mao dẫn nhiệt hạt đơn do nguồn nhiệt laser trong vi kênh và góc tiếp xúc động.

Khoảng trống luận án lựa chọn giải quyết

Luận án giải quyết vấn đề mô phỏng chi tiết động lực học chuyển động và biến dạng của hạt lưu chất đa thành phần (gồm pha ngoài và nhân bên trong) dưới tác động của gradient nhiệt độ và lực mao dẫn nhiệt trong môi trường không biên cũng như trong ống dẫn thu hẹp hình sin, có xét đầy đủ sự khác biệt về độ nhớt, khối lượng riêng, sức căng bề mặt và tương tác giữa các mặt phân cách.

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

Khung lý thuyết và các phương trình chi phối

Toàn bộ miền dòng chảy nhiều pha được biểu diễn bằng một hệ phương trình bảo toàn duy nhất với các đặc tính vật chất thay đổi theo vị trí:

  1. Phương trình Navier-Stokes dạng bảo toàn: $$\frac{\partial (\rho \mathbf{u})}{\partial t} + \nabla \cdot (\rho \mathbf{u} \mathbf{u}) = -\nabla p + \nabla \cdot \left[ \mu (\nabla \mathbf{u} + \nabla \mathbf{u}^T) \right] + \int_f \sigma \kappa \mathbf{n}_f \delta(\mathbf{x} - \mathbf{x}_f) , dS + \mathbf{f}$$ Trong đó $\rho$ là khối lượng riêng, $\mathbf{u}$ là trường vận tốc, $p$ là áp suất, $\mu$ là độ nhớt động lực học, $\sigma$ là hệ số sức căng bề mặt, $\mathbf{n}_f$ là véc tơ pháp tuyến tại mặt phân cách $f$, $\delta(\mathbf{x} - \mathbf{x}_f)$ là hàm Dirac delta biểu thị lực căng bề mặt tập trung tại biên, và $\mathbf{f}$ là số hạng áp đặt điều kiện biên rắn.

  2. Phương trình bảo toàn năng lượng: $$\rho c_p \left( \frac{\partial T}{\partial t} + \mathbf{u} \cdot \nabla T \right) = \nabla \cdot (k \nabla T) + c_p h$$ Với $T$ là nhiệt độ, $k$ là hệ số dẫn nhiệt, $c_p$ là nhiệt dung riêng, và $h$ là số hạng năng lượng biên.

  3. Phương trình bảo toàn khối lượng (phương trình liên tục cho lưu chất không nén được): $$\nabla \cdot \mathbf{u} = 0$$

  4. Phương trình trạng thái đặc tính vật chất: $$\frac{D\rho}{Dt} = 0; \quad \frac{D\mu}{Dt} = 0; \quad \frac{Dk}{Dt} = 0; \quad \frac{Dc_p}{Dt} = 0$$

  5. Mô hình sức căng bề mặt phụ thuộc nhiệt độ: Sức căng bề mặt giảm tuyến tính theo nhiệt độ: $$\sigma(T) = \sigma_0 - \gamma_T (T - T_0)$$ Trong đó $\gamma_T = -\frac{d\sigma}{dT} = \text{const} > 0$, $\sigma_0$ là sức căng bề mặt tại nhiệt độ tham chiếu $T_0$.

Phương pháp tính toán số

  • Phương pháp theo dấu biên (Front-Tracking - FT): Sử dụng lưới Eulerian cố định để giải các phương trình Navier-Stokes và năng lượng; sử dụng lưới Lagrangian di động gồm các phần tử liên kết (danh sách liên kết chứa con trỏ) để theo dõi chính xác vị trí mặt phân cách, tái cấu trúc biên và chuyển đổi thuộc tính bề mặt/lực căng bề mặt lên lưới cố định.
  • Lưới tính toán: Lưới so le MAC (Marker-and-Cell) tiêu chuẩn 2D/đối xứng trục. Áp suất được đặt tại tâm thể tích khống chế; các thành phần vận tốc $u$ và $v$ được đặt tại các biên tương ứng của ô lưới.
  • Lược đồ tích phân:
    • Tích phân thời gian bằng phương pháp chiếu (tách biến) dạng dự đoán - hiệu chỉnh (Predictor-Corrector).
    • Khử phân kỳ vận tốc dẫn đến phương trình Poisson cho áp suất rời rạc.
    • Giải phương trình áp suất bằng phương pháp lặp SOR (Successive Over Relaxation) với hệ số thư giãn $\beta \in [1{,}2; 1{,}5]$.
    • Điều kiện biên trên tường rắn áp dụng thông qua các điểm ảo (ghost cells), ví dụ: $u_{\text{tường}} = \frac{1}{2}(u_{i,1} + u_{i,2}) \Rightarrow u_{i,1} = 2u_{\text{tường}} - u_{i,2}$.
  • Ngôn ngữ lập trình: Toàn bộ chương trình tính toán mô phỏng được tác giả tự xây dựng bằng ngôn ngữ lập trình Fortran.

Bảng so sánh các phương pháp mô phỏng mặt phân cách

Bảng 1: So sánh các đặc tính của phương pháp FT, LS và VOF (theo phân tích trong luận án)

Tính năng và đặc điểm Phương pháp Theo dấu biên (FT) Phương pháp Level-Set (LS) Phương pháp Thể tích lưu chất (VOF)
Bản chất Dịch chuyển tất cả các điểm đánh dấu Lagrangian thuộc biên Dựa vào trường hàm khoảng cách ẩn Dựa vào phương trình vận chuyển phân số thể tích (hoặc nồng độ)
Độ chính xác tái cấu trúc biên Độ chính xác cao (2 hoặc 4 điểm được chọn trên mỗi ô lưới) Chính xác (phương pháp hình học) Không dễ dàng (tính toán biên tuyến tính đơn giản)
Tính toán sức căng bề mặt Dễ dàng Dễ dàng Không dễ dàng
Mở rộng bài toán 3D Tương đối phức tạp/không dễ dàng Không dễ dàng Tương đối đơn giản

Các tham số không thứ nguyên sử dụng trong nghiên cứu

Bảng 2: Danh mục các tham số không thứ nguyên chính

Ký hiệu Tên tham số Định nghĩa / Ý nghĩa vật lý
$Oh$ Số Ohnesorge Tỷ số giữa lực nhớt với lực quán tính và lực căng bề mặt: $Oh = \frac{\mu}{\sqrt{\rho \sigma R}}$
$Re$ Số Reynolds Tỷ số giữa lực quán tính và lực nhớt trong dòng chảy
$Ma$ Số Marangoni Tỷ số giữa vận tốc do dòng Marangoni và vận tốc khuếch tán nhiệt
$Ca$ Số mao dẫn Tỷ số giữa lực cản nhớt và lực căng bề mặt: $Ca = \frac{\mu U}{\sigma}$
$Ar_o, Ar_i$ Tỷ số hình dạng ban đầu Tỷ số giữa chiều dài và đường kính của sợi lưu chất ngoài và trong
$R_{io}$ Tỷ số bán kính Tỷ số bán kính giữa hạt/sợi bên trong và bên ngoài ($R_i / R_o$)
$\sigma_{io}$ Tỷ số sức căng bề mặt Tỷ số sức căng bề mặt giữa biên trong và biên ngoài ($\sigma_i / \sigma_o$)
$\mu_{io}, \mu_{mo}$ Tỷ số độ nhớt Tỷ số độ nhớt giữa các pha lưu chất tương ứng
$d/R_c$ Tỷ số co thắt Tỷ số độ thắt lớn nhất của vùng thu hẹp hình sin so với bán kính ống

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

Chương 1. Tổng quan

  • Giới thiệu các khái niệm cơ bản: lưu chất Newton không nén được, sợi lưu chất, hạt lưu chất đa thành phần, hiện tượng mao dẫn và hiện tượng mao dẫn nhiệt do biến thiên gradient nhiệt độ.
  • Trình bày tổng quan cơ chế hình thành hạt từ sự co lại và phân tách của sợi/tia lưu chất đơn và đa thành phần, cơ chế tách đầu sợi ("end-pinching") và sự tạo thành hạt vệ tinh.
  • Tổng hợp tình hình nghiên cứu lý thuyết, thực nghiệm và mô phỏng số (LS, VOF, LB, FT) về chuyển động mao dẫn nhiệt trên thế giới và ở Việt Nam.
  • So sánh ưu nhược điểm của các phương pháp bắt biên và xác lập tính cấp thiết của việc phát triển mô hình theo dấu biên bằng ngôn ngữ Fortran cho bài toán hạt đa thành phầ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

  • Trình bày chi tiết hệ phương trình bảo toàn Navier-Stokes và phương trình năng lượng kết hợp với mô hình sức căng bề mặt phụ thuộc nhiệt độ.
  • Xây dựng thuật toán giải số: kỹ thuật rời rạc sai phân hữu hạn trên lưới so le MAC, thuật toán tích phân thời gian tách biến dự đoán - hiệu chỉnh, dẫn xuất phương trình Poisson cho áp suất và thuật toán lặp SOR.
  • Trình bày cấu trúc dữ liệu mảng và danh sách liên kết cho các phần tử biên Lagrangian, phương pháp nội suy chuyển đổi thuộc tính giữa lưới biên và lưới cố định.
  • Kiểm chứng phương pháp:
    • So sánh kết quả mô phỏng sự co lại của sợi lưu chất đơn lớp với kết quả tính toán của Notz và Basaran (2004) tại hai chế độ: $Oh = 1{,}0; Ar_o = 15$ và $Oh = 0{,}01; Ar_o = 15$, cho thấy sự trùng khớp chính xác về hình dạng co lại theo thời gian.
    • Kiểm chứng chuyển động của hạt lưu chất đơn dưới tác động của gradient nhiệt độ.

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

  • Khảo sát sự co lại và phân tách của sợi lưu chất đa thành phần trong môi trường lưu chất đẳng nhiệt không biên.
  • Kiểm tra sự hội tụ của lưới tính toán thông qua biến thiên theo thời gian của tham số biến dạng bên trong ($T_i$) và bên ngoài ($T_o$).
  • Khảo sát ảnh hưởng của các tham số không thứ nguyên:
    • Ảnh hưởng của số Ohnesorge ($Oh$): Khảo sát các giá trị $Oh = 0{,}01; 0{,}063; 1{,}0$. Khi $Oh$ nhỏ, lực quán tính chiếm ưu thế gây ra dao động mạnh và phân tách; khi $Oh$ lớn, lực nhớt dập tắt dao động, sợi co lại thành hình cầu mà không phân tách.
    • Ảnh hưởng của tỷ số hình dạng ngoài ($Ar_o$) và trong ($Ar_i$): Khảo sát $Ar_o = 7{,}5$ và $15$; $Ar_i = 5, 15, 20, 30$. Tỷ số hình dạng lớn thúc đẩy hiện tượng tách đầu sợi.
    • Ảnh hưởng của tỷ số bán kính ($R_{io}$): Khảo sát các giá trị $R_{io} = 0{,}2; 0{,}5; 0{,}8$.
    • Ảnh hưởng của tỷ số sức căng bề mặt ($\sigma_{io}$): Khảo sát $\sigma_{io} = 0{,}4$ và $3{,}2$.
    • Ảnh hưởng của tỷ số độ nhớt ($\mu_{32}$ và $\mu_{12}$): Khảo sát từ $0{,}2$ đến $3{,}2$.
  • Xây dựng sơ đồ trạng thái phân tách: Luận án thiết lập các biểu đồ pha chuyển tiếp giữa các chế độ:
    • Chế độ không phân tách ("non").
    • Phân tách sợi bên trong ("i-").
    • Phân tách sợi bên ngoài ("o-").
    • Phân tách hỗn hợp cả trong và ngoài ("i-/o-").

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

  • Mô hình hóa hạt lưu chất đa thành phần đặt trong miền lưu chất có gradient nhiệt độ tuyến tính giữa tường nóng ($T_{\text{nóng}}$) và tường lạnh ($T_{\text{lạnh}}$).
  • Nghiên cứu hội tụ lưới với các độ phân giải $192 \times 384$ và $256 \times 512$.
  • Phân tích động lực học dịch chuyển từ vùng lạnh sang vùng nóng:
    • Ảnh hưởng của số Reynolds ($Re = 1, 10, 30, 50$): Đánh giá vận tốc dịch chuyển trục của hạt ngoài ($U_o$) và hạt trong ($U_i$), sự phân bố sức căng bề mặt $\sigma_i$ dọc theo góc cực $\theta \in [0^\circ, 180^\circ]$ tại vị trí $Z_{ci} = 1{,}7R_o$, và độ lệch tâm tương đối.
    • Ảnh hưởng của số Marangoni ($Ma = 20, 60, 100$): $Ma$ tăng làm biến dạng trường nhiệt độ xung quanh hạt, làm thay đổi vận tốc dịch chuyển cực đại và độ biến thiên sức căng bề mặt $\sigma_o$ tại $Z_{co} = 3R_o$.
    • Ảnh hưởng của tỷ số bán kính ($R_{io} = 0{,}4; 0{,}5; 0{,}6$): Ảnh hưởng đến tương tác thủy động lực học và độ lệch tâm giữa nhân trong và vỏ ngoài.
    • Ảnh hưởng của tỷ số độ nhớt ($\mu_{io} = 0{,}2; 0{,}8; 3{,}2$ và $\mu_{mo} = 0{,}1; 0{,}4; 1{,}6$): Xác định sự cản trở chuyển động và phân bố vận tốc nội tại.

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

  • Thiết lập mô hình ống đối xứng trục có phần thắt hình sin, tọa độ cổ ống $(0, z_c)$ với bán kính $R_n$.
  • Kiểm tra hội tụ lưới bằng hai hệ lưới $256 \times 1024$ và $512 \times 2048$.
  • Khảo sát trường áp suất không thứ nguyên $p_n = p / (0{,}5\rho_o U_c^2)$ và nhiệt độ không thứ nguyên $T_n = (T - T_{\text{cold}}) / (T_{\text{hot}} - T_{\text{cold}})$:
    • Ảnh hưởng của số Marangoni ($Ma = 5, 20, 80$): Xác định vận tốc dịch chuyển trục, thời gian đến (arrival time) và độ lệch tâm lớn nhất.
    • Ảnh hưởng của số mao dẫn ($Ca = 0{,}02; 0{,}04; 0{,}08$): $Ca$ chi phối mức độ biến dạng của màng ngoài khi đi qua vùng cổ thắt.
    • Ảnh hưởng của tỷ số bán kính ($R_{io} = 0{,}3; 0{,}7; 0{,}8$): Tác động đến khả năng chui qua khe hẹp và nguy cơ va chạm giữa nhân và vỏ.
    • Ảnh hưởng của tỷ số co thắt ($d/R_c = 0; 0{,}4; 0{,}6$): Đo lường lực cản hình học và sự thay đổi độ lệch tâm cực đại khi hạt đi qua tiết diện thu hẹp.

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

Đóng góp về mặt lý luận và khoa học

  1. Xây dựng và hoàn thiện chương trình mô phỏng số dựa trên phương pháp theo dấu biên (Front-Tracking) kết hợp kỹ thuật sai phân hữu hạn trên lưới so le MAC bằng ngôn ngữ Fortran, giải quyết bài toán dòng lưu chất đa pha, đa thành phần có truyền nhiệt và mao dẫn nhiệt.
  2. Phân tích định lượng và làm rõ cơ chế biến dạng, phân tách sợi lưu chất đa thành phần trong môi trường đẳng nhiệt, xây dựng các biểu đồ pha phân loại chế độ phân tách (không phân tách, phân tách trong, phân tách ngoài, phân tách hỗn hợp) theo các tham số không thứ nguyên ($Oh, Ar_o, Ar_i, R_{io}, \sigma_{io}, \mu$).
  3. Làm sáng tỏ quy luật chuyển động mao dẫn nhiệt của hạt đa thành phần trong kênh thẳng và ống thu hẹp hình sin, mô tả chi tiết sự biến thiên của độ lệch tâm giữa nhân trong và vỏ ngoài theo các số $Re, Ma, Ca$ và tỷ số hình học.

Đóng góp về mặt thực tiễn

  1. Cung cấp dữ liệu tính toán và cơ sở khoa học hỗ trợ việc thiết kế, kiểm soát chuyển động của các hạt vi lưu chứa dược chất trong các thiết bị vi lỏng phục vụ y sinh (truyền dẫn thuốc điều trị xơ vữa động mạch).
  2. Đề xuất thông số tối ưu cho quy trình công nghệ chế tạo vật liệu lưu chất đa lớp và điều khiển bọt khí/hạt lưu chất trong môi trường vi trọng lực (công nghệ không gian, hệ thống nhiên liệu tên lửa lỏng).

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

Hạn chế của nghiên cứu

  • Mô hình tính toán giới hạn trong không gian 2D và đối xứng trục, chưa mở rộng ra không gian 3D đầy đủ cho các dạng hình học bất đối xứng.
  • Phạm vi dòng chảy chỉ tập trung vào dòng chảy tầng với số Reynolds nhỏ.
  • Chưa xét đến các hiệu ứng chuyển pha nhiệt phức tạp (như đông đặc hoàn toàn hoặc bay hơi) và biến thiên phi tuyến của sức căng bề mặt ở dải nhiệt độ rộng.

Hướng phát triển tiếp

  • Mở rộng chương trình tính toán sang không gian 3 chiều (3D) để mô tả các hiện tượng biến dạng bất đối xứng phức tạp.
  • Nghiên cứu bổ sung các bài toán có truyền nhiệt kèm chuyển pha (đông đặc, nóng chảy) cho các hạt lưu chất kim loại hoặc bán dẫn.
  • Khảo sát các cấu hình kênh vi lỏng có hình học phân nhánh hoặc có chướng ngại vật phức tạp hơn ở dải thông số dòng chảy rộng hơn.

Giá trị tham khảo

  • Đối tượng hưởng lợi:
    • 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 và chất khí, Kỹ thuật Nhiệt, Kỹ thuật Cơ khí, Kỹ thuật Vi lỏng và Toán tin ứng dụng.
    • Kỹ sư và nhà nghiên cứu công nghệ: Trong lĩnh vực thiết kế chip vi lỏng (Lab-on-a-Chip), công nghệ dẫn thuốc y sinh học, và kỹ thuật hàng không - vũ trụ (quản lý chất lỏng trong vi trọng lực).
  • Phần nội dung có giá trị tham khảo cao:
    • Thuật toán số và cấu trúc mã nguồn theo dấu biên (Front-Tracking) kết hợp bộ giải Poisson áp suất SOR trong Chương 2.
    • Biểu đồ pha phân tách của sợi lưu chất đa pha trong Chương 3.
    • Phương pháp phân tích độ lệch tâm và tương tác mao dẫn nhiệt nhiều pha trong Chương 4 và Chương 5.

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

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

Phương pháp Front-Tracking biểu diễn mặt phân cách bằng các phần tử Lagrangian di động trên lưới cố định Eulerian, giúp xác định vị trí biên chính xác mà không bị khuếch tán số. Ưu điểm quan trọng nhất đối với bài toán mao dẫn nhiệt là việc tính toán độ cong và lực căng bề mặt biến thiên theo nhiệt độ được thực hiện trực tiếp và chính xác trên các phần tử biên, thuận lợi hơn so với VOF và LS.

2. Sự co lại của sợi lưu chất đa thành phần trong Chương 3 gồm những chế độ phân tách nào?

Quá trình co lại của sợi lưu chất đa thành phần được phân loại thành 4 chế độ:

  • Chế độ không phân tách ("non"): Sợi co lại thành hạt hình cầu duy nhất.
  • Chế độ phân tách bên trong ("i-"): Chỉ có sợi nhân bên trong bị đứt đoạn.
  • Chế độ phân tách bên ngoài ("o-"): Chỉ có lớp màng bên ngoài bị đứt đoạn.
  • Chế độ phân tách hỗn hợp ("i-/o-"): Cả sợi bên trong và màng bên ngoài đều bị phân tách, có thể tạo thành các hạt lưu chất vệ tinh.

3. Số Ohnesorge ($Oh$) ảnh hưởng như thế nào đến quá trình biến dạng của sợi lưu chất?

Số $Oh$ thể hiện tỷ số giữa lực nhớt với lực quán tính và lực căng bề mặt. Khi $Oh$ nhỏ ($Oh = 0{,}01$), lực quán tính chiếm ưu thế gây ra dao động biên mạnh và thúc đẩy hiện tượng phân tách tạo hạt vệ tinh. Khi $Oh$ lớn ($Oh = 1{,}0$), lực nhớt chiếm ưu thế làm suy giảm dao động, khiến sợi co lại êm dịu về dạng hình cầu ổn định mà không bị phân tách.

4. Cơ chế chuyển động của hạt lưu chất dưới tác động của mao dẫn nhiệt diễn ra theo hướng nào?

Hạt lưu chất dịch chuyển từ vùng nhiệt độ thấp (lạnh) sang vùng nhiệt độ cao (nóng). Nguyên nhân do sức căng bề mặt giảm khi nhiệt độ tăng ($\gamma_T > 0$); sự chênh lệch nhiệt độ dọc bề mặt hạt tạo ra gradient sức căng bề mặt, gây ra ứng suất trượt Marangoni đẩy chất lỏng xung quanh về phía vùng lạnh và đẩy hạt dịch chuyển ngược lại về phía vùng nóng.

5. Những yếu tố nào ảnh hưởng đến độ lệch tâm của hạt đa thành phần khi chuyển động qua ống thu hẹp trong Chương 5?

Độ lệch tâm giữa nhân bên trong và màng ngoài chịu ảnh hưởng trực tiếp bởi: số Marangoni ($Ma$), số mao dẫn ($Ca$), tỷ số bán kính ($R_{io}$), và tỷ số co thắt của kênh ($d/R_c$). Khi hạt đi qua vùng cổ thắt hình sin, sự biến dạng của màng ngoài và gradient áp suất làm thay đổi vị trí tương đối giữa tâm hạt trong ($Z_{ci}$) và tâm hạt ngoài ($Z_{co}$), dẫn đến biến thiên cực đại của độ lệch tâm.

Kết luận

Luận án đã xây dựng thành công chương trình mô phỏng số bằng ngôn ngữ Fortran dựa trên phương pháp theo dấu biên kết hợp sai phân hữu hạn để nghiên cứu động lực học của sợi và hạt lưu chất đa thành phần. Công trình đã làm rõ quy luật biến dạng, phân tách sợi đẳng nhiệt và cơ chế chuyển động mao dẫn nhiệt của hạt đa thành phần trong kênh dẫn thẳng và kênh thu hẹp hình sin. Các kết quả và biểu đồ pha không thứ nguyên thu được cung cấp cơ sở dữ liệu khoa học tin cậy, đóng góp thiết thực cho các ứng dụng điều khiển vi lưu trong công nghệ y sinh và xử lý chất lỏng trong môi trường vi trọng lực.