Tổng quan về luận án

Sự phát triển vượt bậc của công nghệ vi lưu (microfluidics) và cơ học lưu chất y sinh đang mở ra những bước tiến mang tính cách mạng trong các phương pháp điều trị trúng đích, đặc biệt là trong can thiệp xơ vữa động mạch vành và các liệu pháp phân phối thuốc nano. Trong bối cảnh vật lý vi mô và môi trường vi trọng lực (microgravity), khi quy mô không gian giảm xuống từ vài micromet đến vài milimet, tỷ số diện tích trên thể tích tăng vọt khiến trọng lực và lực nổi tự nhiên trở nên thứ yếu, nhường vị trí chi phối cho các hiện tượng mặt ngoài, đặc biệt là lực căng bề mặt và gradient nhiệt độ. Nghiên cứu của luận án được thực hiện trong khuôn khổ đề tài "Nghiên cứu tính toán động lực học hạt lưu chất đa lớp với truyền nhiệt và chuyển pha" (mã số 107.307 do Quỹ Phát triển Khoa học và Công nghệ Quốc gia - NAFOSTED tài trợ), tập trung làm sáng tỏ bản chất thủy nhiệt động lực học của hạt lưu chất đa thành phần (compound droplets) dưới tác động của hiện tượng mao dẫn nhiệt (thermocapillary migration / Marangoni effect).

Khoảng trống nghiên cứu (research gap) cốt lõi trong y văn quốc tế nằm ở chỗ: phần lớn các công trình kinh điển từ Young et al. (1959), Subramanian (1981, 1983) đến Haj-Hariri et al. (1997) chỉ tiếp cận hạt lưu chất đơn thành phần (single droplets/bubbles) trong môi trường chất lỏng vô hạn hoặc dòng chảy Stokes tuyến tính. Các mô hình lý thuyết hiếm hoi về hệ đa pha phức hợp thường áp dụng giả thiết đơn giản hóa triệt để, bỏ qua bước nhảy khối lượng riêng giữa các pha, không xem xét sự tương tác thủy động phi tuyến giữa hạt nhân bên trong (inner core) và vỏ bao bọc bên ngoài (outer shell), đồng thời chưa khảo sát hành vi biến dạng của hạt khi dịch chuyển qua các kênh dẫn vi lưu bị thu hẹp cục bộ (constricted microchannels).

Để giải quyết khoảng trống học thuật này, luận án xác lập ba câu hỏi nghiên cứu (Research Questions) và ba giả thuyết khoa học (Hypotheses) tương ứng:

  • RQ1: Cơ chế biến dạng và bản đồ phân tách pha (pinch-off regime map) của sợi lưu chất đa thành phần trong môi trường đẳng nhiệt chịu sự chi phối như thế nào bởi các tham số không thứ nguyên hình học và động học?
  • RQ2: Sự hiện diện của gradient nhiệt độ tương tác với cấu trúc đa lớp tạo ra cơ chế truyền động mao dẫn nhiệt ra sao đối với trường vận tốc, trường áp suất và độ lệch tâm (eccentricity) giữa nhân trong và vỏ ngoài?
  • RQ3: Mức độ thu hẹp hình học của kênh dẫn (dạng hàm sin) cản trở và làm thay đổi động lực học biến dạng, thời gian dịch chuyển của hạt lưu chất đa thành phần như thế nào dưới tác động đồng thời của số Marangoni ($Ma$) và số mao dẫn ($Ca$)?
  • H1: Quá trình co lại của sợi lưu chất đa thành phần sẽ xuất hiện ít nhất ba chế độ phân tách đặc thù (phân tách nhân trong, phân tách vỏ ngoài và phân tách hỗn hợp) tùy thuộc vào tương quan giữa số Ohnesorge ($Oh$) và tỷ số hình dạng ban đầu ($A_{ro}, A_{ri}$).
  • H2: Sự hiện diện của hạt lưu chất bên trong tạo ra gradient nhiệt độ cục bộ bất đối xứng, làm biến đổi phân bố ứng suất Marangoni trên bề mặt ngoài và gia tốc vận tốc dịch chuyển của toàn bộ cấu trúc hạt so với trường hợp hạt đơn.
  • H3: Tương tác thủy động lực học thành kênh trong ống dẫn thu hẹp gây ra hiện tượng trễ vận tốc phi tuyến và gia tăng độ lệch tâm cực đại, đòi hỏi một ngưỡng số Marangoni tới hạn để khắc phục lực cản nhớt cục bộ.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp chặt chẽ giữa hệ phương trình Navier-Stokes cho dòng đa pha không nén được, phương trình bảo toàn năng lượng với số hạng đối lưu - dẫn nhiệt, và phương trình trạng thái suy giảm tuyến tính của sức căng bề mặt theo nhiệt độ: $$\sigma(T) = \sigma_0 - \sigma_T(T - T_0)$$ Đóng góp đột phá của luận án là việc xây dựng thành công chương trình mô phỏng số tự phát triển bằng ngôn ngữ Fortran dựa trên phương pháp theo dấu biên (Front-Tracking Method - FT) trên lưới so le Marker-and-Cell (MAC). Nghiên cứu thực hiện khảo sát toàn diện không gian tham số: số Reynolds $Re \in [1, 50]$, số Marangoni $Ma \in [5, 100]$, số Ohnesorge $Oh \in [0.01, 1.0]$, số mao dẫn $Ca \in [0.02, 0.08]$, tỷ số bán kính $R_{io} \in [0.2, 0.8]$, và tỷ số độ nhớt $\mu_{io}, \mu_{mo} \in [0.1, 3.2]$ trên các hệ lưới tinh chỉnh từ $192 \times 384$ đến $512 \times 2048$ phần tử. Kết quả mang ý nghĩa nền tảng cho việc thiết kế các chip vi lưu điều khiển hướng đích hạt bọc thuốc và tối ưu hóa quy trình chế tạo vật liệu tiên tiến trong công nghệ vũ trụ.

Literature Review và Positioning

Nghiên cứu về động lực học mặt phân cách và chuyển động mao dẫn nhiệt trải qua nhiều giai đoạn phát triển với các trường phái lý thuyết và thực nghiệm đa dạng. Về động lực học co và phân tách sợi lưu chất đơn, công trình thực nghiệm tiên phong của Stone et al. (1999) đã đặt nền móng khi phát hiện cơ chế đứt gãy ở đầu sợi (end-pinching mechanism). Tiếp nối hướng đi này, Schulkes (1996) và đặc biệt là Notz & Basaran (2004) đã sử dụng phương pháp phần tử hữu hạn nâng cao để xây dựng biểu đồ pha hoàn chỉnh, phân định ranh giới giữa chế độ co không phân tách và phân tách tạo hạt vệ tinh (satellite droplets) dựa trên tỷ số hình dạng ban đầu $A_r$ và số Ohnesorge $Oh$. Gần đây hơn, Castrejon-Pita et al. (2012) và Anthony et al. (2019) đã mở rộng nghiên cứu sang chất lỏng có độ nhớt cao, trong khi Hoepffner & Paré (2013) làm rõ vai trò của các cấu trúc xoáy bên trong giọt trong việc ức chế sự đứt gãy sợi.

Đối với hệ lưu chất đa thành phần, nghiên cứu thực nghiệm bắt đầu từ sự phân tách tia chất lỏng đồng trục của Hertz & Hermanrud (1983), tiếp nối bởi Evangelio et al. (2006) và Yu et al. (2006) trong các cơ cấu hội tụ dòng (flow-focusing devices). Utada et al. (2005) và Pan et al. (2014) đã thiết lập các chế độ tạo giọt đa lớp dạng nước-trong-dầu-trong-nước (W/O/W). Trên phương diện tính toán, các công trình của Vũ et al. (2010, 2013), Zhou et al. (2014) và Zhang et al. (2017) đã mô phỏng quá trình tạo hạt đa thành phần bằng phương pháp trường pha (Phase-Field) và thể tích lưu chất (VOF). Tuy nhiên, hầu hết các tác giả chỉ tập trung vào cơ chế phân tách cơ học dưới tác động của trường vận tốc ngoài, bỏ qua tương tác nhiệt động.

Về chuyển động mao dẫn nhiệt của hạt lưu chất, nghiên cứu lý thuyết bắt đầu từ công trình kinh điển của Young et al. (1959) (mô hình YGB) xác định vận tốc cuối của giọt hình cầu trong giới hạn dòng chảy Stokes ($Re \to 0, Ma \to 0$). Subramanian (1981, 1983) và Balasubramaniam & Chai (1987) đã bổ sung các hiệu ứng truyền nhiệt đối lưu ở số Marangoni hữu hạn bằng kỹ thuật mở rộng tiệm cận. Haj-Hariri et al. (1997) dùng định lý đối ứng Lorentz để đánh giá ảnh hưởng của quán tính yếu. Trong không gian nhiều hạt, Meyyapan & Subramanian (1984), Anderson (1985) và Acrivos et al. (1990) chứng minh rằng khi $Re=Ma=0$, các giọt khí có kích thước bằng nhau không tương tác với nhau. Tuy nhiên, khi $Re, Ma > 0$, các mô phỏng của Nas & Tryggvason (1993, 2003), Yin et al. (2008, 2012) và Kalichetty et al. (2019) đã chỉ ra hiện tượng tương tác thủy nhiệt mạnh mẽ: giọt ở vùng lạnh di chuyển nhanh hơn và có xu hướng bắt kịp giọt ở vùng nóng do sự biến dạng của lớp biên nhiệt.

TIẾN TRÌNH NGHIÊN CỨU LÝ THUYẾT & MÔ PHỎNG MAO DẪN NHIỆT:

Trong y văn tồn tại hai luồng quan điểm và tranh luận khoa học sâu sắc:

  1. Tranh luận về ảnh hưởng của quán tính và số Marangoni cao: Một trường phái (dẫn đầu bởi Zhao et al., 2001 và Brady et al., 2005 bằng phương pháp Level-Set) cho rằng ở số Marangoni lớn, sự bao bọc của các đường đẳng nhiệt quanh mặt phân cách làm suy giảm gradient nhiệt tiếp tuyến, dẫn đến việc vận tốc dịch chuyển luôn đạt trạng thái ổn định và giảm dần khi $Ma$ tăng. Ngược lại, trường phái thứ hai (Yin et al., 2008 và Liu et al., 2013) quan sát thấy hiện tượng dao động vận tốc dạng sóng (oscillatory velocity) do sự tích tụ và giải phóng năng lượng nhiệt không ổn định phía sau đuôi hạt.
  2. Tranh luận về tương tác hạt - vách trong không gian hẹp: Lu et al. (2014) cho rằng số mao dẫn $Ca$ không có ảnh hưởng đáng kể đến thời gian di chuyển nếu khe hẹp lớn hơn đường kính hạt, trong khi Capobianchi et al. (2017) khẳng định hiệu ứng tắc nghẽn (confinement blockage) luôn tạo ra lớp biên thủy động làm biến dạng phi đối xứng nghiêm trọng màng chất lỏng bao quanh hạt.

So sánh với hai nghiên cứu quốc tế tiêu biểu:

  • So với công trình của Morton et al. (2012) (nghiên cứu giải tích tiệm cận hạt đa lớp trong không gian vô hạn, bỏ qua quán tính và giả định hình dạng giọt cầu cố định), luận án đã vượt bậc khi giải quyết trọn vẹn phương trình Navier-Stokes phi tuyến đầy đủ kết hợp truyền nhiệt đối lưu, mô phỏng chính xác sự biến dạng hình học động và độ lệch tâm của nhân trong.
  • So với nghiên cứu của Jadhav & Ghosh (2013) (khảo sát hạt đa thành phần trong dòng Poiseuille đẳng nhiệt), luận án mở rộng sang trường đa vật lý có tương tác mao dẫn nhiệt, đồng thời khảo sát hình học thành kênh thu hẹp dạng sin thay vì kênh dẫn thẳng đồng nhất.

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

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

Luận án đã mở rộng căn bản lý thuyết mao dẫn nhiệt kinh điển của Young, Goldstein & Block (1959) và lý thuyết mất ổn định mao dẫn Rayleigh-Plateau sang hệ thống đa pha, đa lớp phân cấp (multi-layered nested interfaces).

Mô hình lý thuyết của luận án xác lập hệ thống 4 mệnh đề khoa học (Theoretical Propositions):

  • Proposition 1: Trong quá trình co sợi đa thành phần, sự tương tác nhớt giữa ba pha $(\mu_1, \mu_2, \mu_3)$ và sức căng bề mặt tại hai mặt phân cách $(\sigma_i, \sigma_o)$ thiết lập một trường sóng áp suất mao dẫn phân cấp, quyết định vị trí thắt nút (necking) xảy ra độc lập hoặc đồng thời tại biên trong và biên ngoài.
  • Proposition 2: Vận tốc dịch chuyển mao dẫn nhiệt của hạt đa thành phần là một hàm hợp của tỷ số bán kính $R_{io}$, tỷ số dẫn nhiệt $k_{io}, k_{mo}$, và tỷ số nhớt $\mu_{io}, \mu_{mo}$. Sự hiện diện của nhân trong làm gia tăng độ dốc nhiệt độ hiệu dụng tại bề mặt ngoài, qua đó tạo ra lực đẩy Marangoni vượt trội so với hạt đơn có cùng thể tích tương đương.
  • Proposition 3: Độ lệch tâm của nhân trong không duy trì vị trí đồng tâm tĩnh mà tiến hóa theo một hàm động học phụ thuộc vào số Reynolds ($Re$) và số Marangoni ($Ma$); tại $Re \ge 30$, lực quán tính tạo ra chuyển dịch vị trí nhân về phía cực ấm của vỏ ngoài.
  • Proposition 4: Trong ống dẫn thu hẹp, độ thắt hình học $d/R_c$ kích hoạt sự biến đổi ứng suất nhớt thành vách cực đại, gây ra hiện tượng dịch chuyển pha giữa đỉnh vận tốc và đỉnh biến dạng hình học.
KHUNG KHÁI NIỆM TƯƠNG TÁC ĐA LỚP & ĐA VẬT LÝ:

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

Khung phân tích của nghiên cứu là sự tổng hợp đồng thời của 4 lý thuyết vật lý cốt lõi:

  1. Lý thuyết thủy động lực học Navier-Stokes: Mô tả dòng chảy tầng, nhớt, không nén được của môi trường nhiều pha.
  2. Lý thuyết truyền nhiệt đối lưu - khuếch tán: Mô tả sự vận chuyển năng lượng nhiệt qua các môi trường có hệ số dẫn nhiệt $k$ và nhiệt dung riêng $c_p$ biến thiên theo từng pha.
  3. Lý thuyết nhiệt động lực học bề mặt Marangoni: Thiết lập mối liên hệ trực tiếp giữa gradient nhiệt độ tiếp tuyến và ứng suất cắt bề mặt: $$\mathbf{\tau} \cdot \mathbf{n}_f = \nabla_s \sigma = -\sigma_T \nabla_s T$$
  4. Lý thuyết mất ổn định bề mặt mao dẫn: Phân tích sự cân bằng giữa lực căng bề mặt co rút và lực cản nhớt/quán tính thông qua các nhóm tham số không thứ nguyên $Oh, Ca, Re, Ma$.

Điều kiện biên (Boundary Conditions) và giới hạn bài toán được thiết lập chuẩn xác: miền tính toán hai chiều trục đối xứng $(r, z)$, tường trên và tường dưới áp đặt điều kiện không trượt $\mathbf{u} = 0$ với nhiệt độ cố định $T_{nóng}$ và $T_{lạnh}$; các mặt biên bên áp dụng điều kiện đối xứng hoặc tuần hoàn phù hợp với từng bài toán cụ thể.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan bản thể luận thực chứng (Positivist Paradigm) và phương pháp luận mô hình hóa số tất định (Deterministic Computational Modeling). Để khắc phục các hạn chế cố hữu của các phương pháp nắm bắt biên truyền thống, luận án áp dụng phương pháp theo dấu biên (Front-Tracking Method - FT) trên lưới thể tích hữu hạn so le MAC tiêu chuẩn.

Lý do lựa chọn phương pháp FT thay vì Level-Set (LS) hay Volume of Fluid (VOF) bắt nguồn từ bản chất cơ học của bài toán mao dẫn nhiệt:

  • Phương pháp VOF tuy bảo toàn khối lượng hoàn hảo nhưng việc tính toán độ cong và đạo hàm bậc cao của sức căng bề mặt phụ thuộc nhiệt độ thường sinh ra dòng chảy ký sinh (parasitic currents) lớn.
  • Phương pháp LS giải quyết tốt thay đổi topo nhưng hàm khoảng cách dấu (signed distance function) có xu hướng làm tiêu tán độ sắc nét của gradient nhiệt tại mặt phân cách.
  • Phương pháp FT kết hợp ưu điểm của việc sử dụng một hệ phương trình liên tục duy nhất trên lưới Eulerian cố định cho toàn miền với một lưới Lagrangian di động bao gồm các điểm đánh dấu liên kết (marker points) để theo dõi tường minh vị trí mặt phân cách. Điều này cho phép tính toán chính xác lực căng bề mặt và chuyển đổi mượt mà các bước nhảy vật tính chất lỏng thông qua hàm phân phối delta xấp xỉ.

Quy trình nghiên cứu rigorous

Cấu trúc thuật toán mô phỏng được lập trình trong Fortran bao gồm các module tương hỗ chặt chẽ:

  1. Mô tả mặt phân cách: Biên giữa các pha lưu chất được rời rạc hóa thành chuỗi các phần tử đường liên kết bởi các điểm nút Lagrangian. Cấu trúc dữ liệu sử dụng danh sách liên kết kép (doubly linked list) cho phép chèn, xóa và tái cấu trúc biên (re-meshing/re-structuring) linh hoạt khi sợi hoặc hạt bị biến dạng kéo dài.
  2. Chuyển đổi đặc tính vật lý: Khối lượng riêng $\rho$, độ nhớt $\mu$, hệ số dẫn nhiệt $k$, và nhiệt dung riêng $c_p$ được tái cấu trúc trên toàn bộ lưới Eulerian thông qua hàm chỉ thị pha $I(\mathbf{x})$, ví dụ: $$\rho(\mathbf{x}) = \rho_o + (\rho_m - \rho_o) I_o(\mathbf{x}) + (\rho_i - \rho_m) I_i(\mathbf{x})$$
  3. Phân phối lực căng bề mặt: Lực căng bề mặt trên mỗi phần tử biên Lagrangian $\delta \mathbf{F}s$ được chuyển giao sang lưới Eulerian cố định thông qua hàm delta Dirac rời rạc: $$\mathbf{F}{sv}(\mathbf{x}) = \sum_l \mathbf{F}_{s,l} D_h(\mathbf{x} - \mathbf{x}_l) \Delta s_l$$
  4. Giải phương trình Navier-Stokes bằng phương pháp tách bước (Projection Method):
    • Bước dự đoán (Predictor step): Bỏ qua gradient áp suất, tính toán trường vận tốc trung gian $\mathbf{u}^$: $$\frac{\rho^{n+1} \mathbf{u}^ - \rho^n \mathbf{u}^n}{\Delta t} = - \nabla_h \cdot (\rho^n \mathbf{u}^n \mathbf{u}^n) + \nabla_h \cdot \left[ \mu^n (\nabla_h \mathbf{u}^n + \nabla_h^T \mathbf{u}^n) \right] + \mathbf{F}_{sv}^n$$
    • Bước giải phương trình áp suất Poisson: $$\nabla_h \cdot \left( \frac{1}{\rho^{n+1}} \nabla_h p \right) = \frac{1}{\Delta t} \nabla_h \cdot \mathbf{u}^*$$ Phương trình áp suất được giải lặp bằng phương pháp thư giãn liên tiếp SOR (Successive Over-Relaxation) với hệ số gia tốc $\beta \in [1.2, 1.5]$ nhằm đảm bảo tốc độ hội tụ nhanh và triệt tiêu sai số phân kỳ vận tốc $\nabla \cdot \mathbf{u}^{n+1} = 0$.
    • Bước hiệu chỉnh (Corrector step): $$\mathbf{u}^{n+1} = \mathbf{u}^* - \frac{\Delta t}{\rho^{n+1}} \nabla_h p$$
  5. Giải phương trình năng lượng và cập nhật vị trí biên: Trường nhiệt độ $T^{n+1}$ được cập nhật tường minh, sau đó các điểm biên Lagrangian được di chuyển theo vận tốc nội suy từ lưới Eulerian: $$\mathbf{x}_f^{n+1} = \mathbf{x}_f^n + \mathbf{V}_n \Delta t$$
SƠ ĐỒ KHỐI THUẬT TOÁN FRONT-TRACKING GIẢI COUPLED NS-ENERGY:

Data và phân tích

Tính xác thực và độ tin cậy khoa học của chương trình mô phỏng được kiểm chứng qua 2 bước nghiêm ngặt:

  1. Kiểm chứng độc lập (Benchmark Validation): So sánh hình dạng co lại của sợi lưu chất đơn lớp tại các mốc thời gian không thứ nguyên khác nhau với kết quả giải tích và phần tử hữu hạn của Notz & Basaran (2004) tại $Oh = 1.0; 0.01$ và $A_{ro} = 15$. Độ trùng khớp gần như tuyệt đối giữa hai đường biểu diễn đã khẳng định tính chính xác của thuật toán.
  2. Kiểm tra tính độc lập của lưới (Mesh Convergence Study):
    • Trong bài toán chuyển động hạt đa thành phần trong không gian mở: khảo sát trên hai hệ lưới $192 \times 384$ và $256 \times 512$. Sự biến thiên của vận tốc vỏ ngoài $U_o$, vận tốc nhân trong $U_i$, cũng như tọa độ tâm $Z_{co}, Z_{ci}$ giữa hai cấp lưới có độ sai lệch dưới $1.2%$.
    • Trong bài toán hạt di chuyển qua ống thu hẹp: kiểm chứng trên lưới $256 \times 1024$ và lưới siêu mịn $512 \times 2048$. Biên dạng hình học của hạt tại vùng cổ thắt hoàn toàn đồng nhất.

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

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

Quá trình phân tích dữ liệu mô phỏng số đã mang lại 4 phát hiện mang tính đột phá về cơ học lưu chất:

BẢN ĐỒ CÁC CHẾ ĐỘ PHÂN TÁCH SỢI LƯU CHẤT ĐA THÀNH PHẦN:
Tỷ số hình dạng ngoài ($A_{ro}$)
  0.001       0.01        0.1        1.0

  1. Thiết lập bản đồ chế độ phân tách sợi lưu chất đa thành phần trong môi trường đẳng nhiệt: Lần đầu tiên, một hệ thống 6 sơ đồ trạng thái pha toàn diện $(A_{ro} \text{ vs } Oh, A_{ro} \text{ vs } R_{io}, A_{ro} \text{ vs } A_{ri}, A_{ro} \text{ vs } \mu_{32}, A_{ro} \text{ vs } \sigma_{io}, A_{ro} \text{ vs } \mu_{12})$ được xây dựng chi tiết. Luận án chỉ ra sự tồn tại của 4 chế độ động lực học rõ rệt:

    • Chế độ không phân tách ("non"): Xảy ra khi $Oh$ lớn ($Oh \ge 1.0$) hoặc $A_{ro}$ nhỏ ($A_{ro} \le 7.5$), lực nhớt dập tắt các dao động mao dẫn, sợi co cụm êm dịu thành hạt đa thành phần hình cầu.
    • Chế độ phân tách sợi bên trong ("i-"): Xảy ra khi $R_{io}$ nhỏ hoặc $\sigma_{io}$ lớn, màng trong bị thắt nút trước tạo thành chuỗi hạt nhân đơn lơ lửng trong vỏ ngoài.
    • Chế độ phân tách sợi bên ngoài ("o-"): Xuất hiện khi độ nhớt vỏ ngoài thấp, hình thành hạt vệ tinh ngoài bao bọc sợi trong.
    • Chế độ phân tách hỗn hợp ("i-/o-"): Đứt gãy phức hợp đồng thời ở cả hai mặt phân cách ở tỷ số hình dạng cao ($A_{ro} = 15 - 30, Oh = 0.01$).

  2. Cơ chế gia tốc vận tốc và biến thiên độ lệch tâm dưới tác động mao dẫn nhiệt: Khi hạt lưu chất đa thành phần chuyển động từ vùng lạnh sang vùng nóng, lực kéo Marangoni sinh ra do gradient sức căng bề mặt $\partial \sigma / \partial T < 0$ đẩy lớp chất lỏng sát biên về phía sau, tạo phản lực đẩy hạt về phía trước. Dữ liệu chứng minh rằng vận tốc cực đại của vỏ ngoài $U_{o,\max}$ tăng tuyến tính theo số Marangoni ($Ma = 20 \to 100$) và tăng phi tuyến theo số Reynolds ($Re = 1 \to 50$). Đặc biệt, sự hiện diện của hạt bên trong làm tăng vận tốc di chuyển tổng thể so với hạt đơn do làm tăng độ méo cục bộ của các đường đẳng nhiệt. Độ lệch tâm trục $e = (Z_{ci} - Z_{co})/R_o$ ban đầu tăng nhanh, sau đó đạt trạng thái tiệm cận ổn định; giá trị độ lệch tâm tỉ lệ thuận với tỷ số bán kính $R_{io}$ và giảm khi tỷ số độ nhớt $\mu_{io}$ tăng.

  3. Hiện tượng trễ thủy động và biến dạng tới hạn trong ống dẫn thu hẹp dạng sin: Khi hạt đi qua cổ thắt ống vi lưu, lực cản nhớt từ thành ống tăng vọt làm suy giảm vận tốc trục. Tuy nhiên, nếu số Marangoni đủ lớn ($Ma = 80$), lực đẩy mao dẫn nhiệt thắng được lực cản hình học, rút ngắn đáng kể "thời gian đến" (arrival time) của hạt tại hạ lưu. Số mao dẫn $Ca$ đóng vai trò quyết định mức độ biến dạng: tại $Ca = 0.08$, hạt bị kéo dài thành dạng elip thon nhọn dọc theo trục đối xứng, trong khi tại $Ca = 0.02$, hạt duy trì hình dạng gần như hình cầu nhưng chịu sự suy giảm vận tốc mạnh hơn khi tiếp cận tâm cổ thắt ($z = z_c$).

  4. Tác động của mức độ thu hẹp thành ống ($d/R_c$): Khảo sát tỷ số co thắt từ $d/R_c = 0$ (ống thẳng) đến $d/R_c = 0.6$ (thu hẹp $60%$) chứng minh rằng độ lệch tâm cực đại giữa nhân và vỏ đạt đỉnh ngay trước khi hạt đi qua tiết diện hẹp nhất, sau đó hồi phục nhanh chóng nhờ lực hồi phục mao dẫn.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình toán lý hoàn chỉnh cho dòng chảy đa pha tương tác nhiệt động, lấp đầy khoảng trống về động lực học đa lớp phân cấp trong cơ học vi lưu.
  • Về mặt phương pháp luận: Khẳng định tính ưu việt và độ tin cậy của mã nguồn Fortran tự phát triển dựa trên thuật toán Front-Tracking kết hợp giải Poisson SOR, cung cấp một công cụ mô phỏng số mạnh mẽ có thể chuyển giao cho các bài toán dòng chảy có chuyển pha (solidification/melting) và dòng khí-lỏng phức tạp.
  • Về mặt ứng dụng thực tiễn: Cung cấp cơ sở khoa học chính xác để thiết kế các hệ thống vi lưu y sinh (Bio-MEMS / Lab-on-a-chip). Bằng cách điều khiển nguồn nhiệt cục bộ (như chiếu xạ vi laser), các kỹ sư y sinh có thể kiểm soát chính xác quỹ đạo chuyển động, phân tách liều lượng, và kích hoạt giải phóng thuốc trúng đích tại các vị trí mảng xơ vữa động mạch mà không cần sử dụng các bộ phận truyền động cơ học phức tạp.
  • Trong công nghệ không gian: Kết quả mô phỏng là tiền đề để tính toán hệ thống thu gom bong bóng khí trong sản xuất thủy tinh quang học tinh khiết ở môi trường không trọng lực và kiểm soát bọt khí trong hệ thống nhiên liệu tên lửa lỏng.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu (Boundary Conditions & Limitations):

  1. Mô hình số hiện tại được xây dựng trên giả định đối xứng trục hai chiều (2D axisymmetric), chưa nắm bắt được các cấu trúc bất ổn định ba chiều (3D non-axisymmetric instabilities) xuất hiện ở số Reynolds hoặc Marangoni rất cao.
  2. Giả thiết các pha lưu chất đều là chất lỏng thuần Newton và không nén được; chưa xét đến đặc tính phi Newton (viscoelastic / shear-thinning) vốn rất phổ biến ở dịch cơ thể và máu người.
  3. Nghiên cứu tập trung vào gradient nhiệt độ thuần túy, chưa tích hợp hiện tượng truyền khối (mass transfer) và phản ứng chuyển pha hóa - lý tại mặt phân cách.

Chương trình nghiên cứu tương lai (Future Research Agenda) 5-10 năm bao gồm:

  • Mở rộng thuật toán Front-Tracking sang không gian 3 chiều đầy đủ (3D DNS) trên nền tảng tính toán song song hiệu năng cao (MPI/OpenMP).
  • Tích hợp mô hình lưu biến phi Newton (Carreau-Yasuda hoặc Oldroyd-B) để mô phỏng chính xác sự lưu thông của hạt thuốc trong hệ tuần hoàn mạch máu thực tế.
  • Khảo sát các cơ chế kích hoạt hỗn hợp: kết hợp mao dẫn nhiệt với điện di (electrophoresis) hoặc sóng âm vi lưu (acoustophoresis) nhằm tối ưu hóa khả năng dẫn hướng vi hạt.

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

Luận án tạo ra tác động học thuật sâu sắc với khả năng trích dẫn cao trong các tạp chí quốc tế đầu ngành thuộc danh mục ISI/Scopus (Q1) như Journal of Fluid Mechanics, Physics of Fluids, và International Journal of Heat and Mass Transfer. Công trình cung cấp hệ thống dữ liệu chuẩn (benchmark data) có độ chính xác cao cho cộng đồng cơ học tính toán toàn cầu.

Về mặt kinh tế - xã hội, nghiên cứu góp phần trực tiếp vào chiến lược phát triển công nghệ cao quốc gia trong lĩnh vực sản xuất thiết bị y tế thông minh và công nghiệp vật liệu bán dẫn tiên tiến.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận mã nguồn mô phỏng tiên tiến, phương pháp luận nghiên cứu thực chứng và các khoảng trống học thuật về cơ học hạt đa pha.
  • Giảng viên & Giáo sư chuyên ngành Cơ học lưu chất / Truyền nhiệt: Tài liệu tham khảo chuẩn mực trong đào tạo sau đại học về dòng chảy đa pha và phương pháp số.
  • Kỹ sư R&D công nghệ vi lưu & Dược phẩm sinh học: Ứng dụng các biểu đồ tham số để chế tạo chip định lượng thuốc và kiểm soát độ bền vi nang bọc dược chất.
  • Cơ quan quản lý khoa học & Quỹ NAFOSTED: Minh chứng rõ nét cho hiệu quả đầu tư nghiên cứu cơ bản định hướng ứng dụng thực tiễn chất lượng cao tại Việt Nam.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp độc đáo nhất là việc mở rộng lý thuyết mao dẫn nhiệt kinh điển YGB (Young et al., 1959) và lý thuyết mất ổn định co sợi Notz & Basaran (2004) sang hệ hạt và sợi đa thành phần phân cấp. Luận án đã chứng minh cơ chế tương tác nhiệt - thủy động lực học phi tuyến giữa nhân trong và vỏ ngoài, xác lập quy luật gia tốc vận tốc và biến thiên độ lệch tâm động mà các lý thuyết đơn pha trước đây không thể giải thích.

2. Đột phá phương pháp luận so với các nghiên cứu quốc tế trước đây? Trả lời: So với các nghiên cứu dùng phương pháp Level-Set của Haj-Hariri et al. (1997) hay VOF của Samareh et al. (2014), luận án đã phát triển thành công thuật toán Front-Tracking trên nền lưới MAC so le với bộ giải áp suất Poisson SOR ($\beta = 1.2 - 1.5$) viết bằng Fortran. Phương pháp này tính toán chính xác lực căng bề mặt Marangoni trực tiếp từ các phần tử Lagrangian mà không làm nhòe mặt phân cách, đồng thời tiết kiệm tài nguyên tính toán hơn so với mô hình 3D LBM của Liu et al. (2013).

3. Phát hiện bất ngờ nhất (counter-intuitive finding) có dữ liệu minh chứng? Trả lời: Phát hiện bất ngờ nhất là sự hiện diện của hạt nhân bên trong không làm cản trở hay làm chậm chuyển động do ma sát nội tại, mà ngược lại, làm tăng vận tốc di chuyển tổng thể của hạt đa thành phần so với hạt đơn có cùng kích thước hình học khi đặt trong cùng gradient nhiệt độ. Điều này được minh chứng qua dữ liệu phân bố nhiệt: sự chênh lệch hệ số dẫn nhiệt $k_{io}$ và nhiệt dung riêng $c_{io}$ làm gia tăng cục bộ độ dốc nhiệt độ tiếp tuyến tại mặt ngoài, tạo ra lực đẩy Marangoni mạnh hơn.

4. Quy trình lặp lại nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không? Trả lời: Luận án cung cấp tường minh toàn bộ hệ phương trình vi phân toán tử, quy trình rời rạc hóa thời gian/không gian (bước dự đoán - hiệu chỉnh), cấu trúc dữ liệu danh sách liên kết kép của biên Lagrangian, phương pháp nội suy hàm delta $D_h$, và các điều kiện biên vật lý. Mọi nhà nghiên cứu độc lập đều có thể tái lập mã nguồn và mô phỏng lại chính xác các kết quả đã công bố.

5. Bản kế hoạch nghiên cứu 10 năm được định hình như thế nào? Trả lời: Kế hoạch 10 năm tập trung vào 3 trụ cột: (1) Nâng cấp mã nguồn lên 3D DNS đa vi xử lý GPU; (2) Tích hợp mô hình sinh lý học dòng máu và giải phóng dược chất vi nang; (3) Chế tạo thử nghiệm chip vi lưu thực tế tại phòng thí nghiệm để kiểm chứng chéo với kết quả mô phỏng số.

Kết luận

Luận án tiến sĩ đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 5 đóng góp học thuật cốt lõi:

  1. Xây dựng thành công bộ chương trình mô phỏng số hoàn chỉnh 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 lưới so le MAC và kỹ thuật sai phân hữu hạn, giải quyết trọn vẹn bài toán ghép đôi Navier-Stokes và bảo toàn năng lượng đa pha.
  2. Thiết lập hệ thống bản đồ chế độ phân tách sợi đa thành phần đầu tiên trong môi trường đẳng nhiệt, phân định chính xác ranh giới của 4 chế độ động lực học ("non", "i-", "o-", "i-/o-") theo các số không thứ nguyên $Oh, A_{ro}, A_{ri}, R_{io}, \sigma_{io}, \mu_{32}, \mu_{12}$.
  3. Làm sáng tỏ cơ chế mao dẫn nhiệt của hạt lưu chất đa thành phần, chứng minh định lượng ảnh hưởng của các số $Re, Ma, R_{io}, \mu_{io}$ đến quy luật gia tốc vận tốc hạt và độ lệch tâm động của nhân trong.
  4. Khám phá quy luật tương tác hạt - thành vách trong ống dẫn thu hẹp dạng sin, xác định mối quan hệ tương hỗ giữa số Marangoni $Ma$, số mao dẫn $Ca$ và mức độ co thắt $d/R_c$ lên sự biến dạng hình học và thời gian dịch chuyển của hạt.
  5. Mở ra hướng ứng dụng đột phá trong công nghệ dẫn truyền thuốc vi lưu trúng đích điều trị tim mạch và kiểm soát bọt khí trong công nghệ không gian vi trọng lực.

Nghiên cứu khẳng định vị thế tiên phong của các nhà khoa học Việt Nam trong lĩnh vực cơ học lưu chất tính toán tiên tiến, tạo nền tảng vững chắc cho các phát triển công nghệ y sinh và vật liệu mới trong tương lai.