Tổng quan luận án

Luận án tiến sĩ chuyên ngành Cơ học chất lỏng và chất khí (Mã số: 9.08) với đề tài "Nghiên cứu tính toán động lực học hạt lưu chất đa thành phần có truyền nhiệt, chuyển pha" được thực hiện 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. Nghiên cứu tập trung vào việc mô phỏng số quá trình truyền nhiệt và chuyển pha (cụ thể là quá trình hóa rắn) của hạt lưu chất rỗng – một dạng hạt lưu chất đa thành phần chứa nhân khí bên trong.

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

Hiện tượng chuyển pha của các hạt lưu chất xuất hiện phổ biến trong tự nhiên và nhiều quy trình công nghiệp:

  • Chế tạo pin năng lượng mặt trời: Trong công nghệ chế tạo pin mặt trời từ các hạt bán dẫn hình cầu Silic, vùng có chức năng biến đổi quang điện chủ yếu nằm ở lớp mỏng bề mặt ngoài của hạt. Việc sử dụng hạt bán dẫn đa thành phần rỗng thay thế cho hạt đặc giúp giảm khối lượng vật liệu bán dẫn đắt tiền và hạ giá thành sản xuất mà vẫn duy trì hiệu suất chuyển đổi năng lượng.
  • Hiện tượng đóng băng trong kỹ thuật hàng không và năng lượng gió: Khi các hạt nước trong không khí (có thể lẫn bọt khí) va chạm và bám dính vào bề mặt cánh máy bay hoặc cánh tuabin gió có nhiệt độ thấp hơn nhiệt độ đóng băng, chúng sẽ bị hóa rắn. Sự hình thành lớp băng này làm biến dạng hình học của biên dạng khí động, làm suy giảm nghiêm trọng lực nâng, tăng lực cản, giảm hiệu suất vận hành và là nguyên nhân trực tiếp dẫn tới nhiều sự cố tai nạn hàng không.

Mặc dù các nghiên cứu lý thuyết, thực nghiệm và mô phỏng số về sự hóa rắn của hạt lưu chất đơn thành phần đã được thực hiện khá phong phú, các nghiên cứu chuyên sâu về động lực học truyền nhiệt và chuyển pha của hạt lưu chất đa thành phần (đặc biệt là hạt rỗng có chứa bóng khí) vẫn còn rất hạn chế do tính chất phi tuyến phức tạp tại các biên phân cách pha và sự xuất hiện của các điểm chập ba pha bên trong lẫn bên ngoài hạt. Do đó, việc xây dựng mô hình tính toán số để khảo sát chi tiết cơ chế vật lý của quá trình này là yêu cầu cấp thiết.

Mục tiêu nghiên cứu

  1. Mục tiêu tổng quát: Nghiên cứu quá trình truyền nhiệt và chuyển pha (quá trình hóa rắn) của hạt lưu chất rỗng bằng phương pháp theo dấu biên (Front-Tracking Method).
  2. Mục tiêu cụ thể:
    • Xây dựng mô hình tính toán mô phỏng số cho bài toán hạt lưu chất rỗng hóa rắn trên bề mặt phẳng lạnh và hạt lưu chất rỗng lơ lửng chịu tác động của dòng đối lưu cưỡng bức.
    • Khảo sát có hệ thống ảnh hưởng của các thông số không thứ nguyên và thông số hình học đến động học biên hóa rắn, hình dạng hình học cuối cùng và thời gian hóa rắn hoàn toàn của hạt lưu chất rỗng.

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

  • Đối tượng nghiên cứu: Hạt lưu chất rỗng (hạt đa thành phần đơn nhân khí) có kích thước nhỏ từ vài trăm micromet đến vài milimet.
  • Phạm vi nghiên cứu:
    • Giả thiết lưu chất là chất lỏng và chất khí không nén được, không trộn lẫn và tuân theo định luật lưu chất Newton.
    • Sử dụng hệ phương trình Navier-Stokes kết hợp phương trình bảo toàn năng lượng trong hệ tọa độ trụ đối xứng ($Orz$).
    • Các đặc tính nhiệt vật lý trong từng pha (khối lượng riêng, độ nhớt, nhiệt dung riêng, hệ số dẫn nhiệt) được xem là hằng số.

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

Luận án đã tổng hợp và phân loại các công trình nghiên cứu liên quan đến hạt đa thành phần và bài toán chuyển pha theo các nhóm nội dung chính:

Ứng dụng của hạt lưu chất đa thành phần

  • Thực phẩm và mỹ phẩm: Pulatsü và cộng sự [21] nghiên cứu cấu trúc hạt nước-dầu-nước; Buyukkestelli và cộng sự [3] ứng dụng mô hình này để tăng cảm nhận độ ngọt lên 75% với cùng hàm lượng sucrose; Lee và cộng sự [22], Stasse và cộng sự [23], Paulo và Santos [24] ứng dụng hạt dầu-nước-dầu để bao bọc retinol, chất béo và giải phóng axit caffeic.
  • Y dược: McClements [1], Iqbal và cộng sự [25], Nabi-Meibodi và cộng sự [26], Peres và cộng sự [27] khảo sát việc đóng gói, vận chuyển thuốc ưa nước, baclofene và lipid rắn.
  • Năng lượng và kỹ thuật: Biancardo và cộng sự [5] nghiên cứu hạt bán dẫn Silic trong pin mặt trời với hiệu suất chuyển đổi đạt 13,5% dưới bức xạ $100\text{ mW/cm}^2$; các nghiên cứu về đóng băng cánh máy bay và cánh tuabin gió của [10], [11], [12], [13], [14]; ứng dụng lọc bụi [7] và vật liệu hấp thụ âm thanh [6].

Sự hình thành hạt lưu chất đa thành phần

  • Cơ chế tách hạt từ dòng tia qua các vòi phun đồng trục hoặc kênh vi lưu cắt nhau được chứng minh qua các công trình của Chen và cộng sự [16], Ji và cộng sự [17], Vũ và cộng sự [18], Bhagat và cộng sự [19].
  • Cơ chế co lại của sợi lưu chất đa thành phần dưới tác động của sức căng bề mặt và số Ohnesorge ($Oh$) được mô tả bởi Hồ và cộng sự [20].

Các công trình nghiên cứu về hóa rắn hạt lưu chất

  • Nghiên cứu lý thuyết: Snoijer và Brunet [30] phát triển mô hình giải tích cho hạt nước đơn, xác định tỉ số khối lượng riêng đặc biệt $\rho_c = 3/4$ và giải thích sự hình thành đỉnh chóp nhọn; Zhang và cộng sự [31] tính toán chiều cao biên hóa rắn lý thuyết là 2,16 so với thực nghiệm là 2,0; Tembely và Dolatabadi [32] xây dựng mô hình một chiều xét đến độ võng của biên hóa rắn; Anderson và cộng sự [33] khảo sát động học góc phát triển $\phi_i$ dựa trên các vận tốc trượt $V_v, V_s$ tại điểm chập ba pha.
  • Nghiên cứu thực nghiệm: Huang và cộng sự [34] khảo sát góc ướt trong dải $76,0^\circ - 154,9^\circ$; Pan và cộng sự [35] làm thực nghiệm trên mặt nghiêng $30^\circ$ (góc ướt $77^\circ - 145^\circ$); Zhang và cộng sự [36], Ju và cộng sự [37] nghiên cứu trên bề mặt cầu; Satunkin [38] đo đạc góc phát triển của Silic, Gecmani và InSb; Bhagat và cộng sự [19] thực hiện hóa rắn hạt Gecmani rỗng và quan sát thấy chóp nhọn ở mặt ngoài nhưng mặt trong bóng khí nhẵn phẳng; Nadler và cộng sự [43] (hợp kim Al-2%Ca), Kendall và cộng sự [44].
  • Nghiên cứu mô phỏng số: Shetabivash và cộng sự [45] (phương pháp Level-Set chia biên âm/dương); Vahab và cộng sự [46] (Level-Set kết hợp Moment of Fluid); Sun và cộng sự [47] (Lattice Boltzmann); Zhou và cộng sự [48] (Volume of Fluid cho tia kim loại); Vũ và cộng sự [49], [50] (Front-Tracking cho hạt đơn kéo dãn bởi trọng trường); Hua và cộng sự [51], Vũ và cộng sự [52], [54], [55], Chen và cộng sự [53], Phạm và cộng sự [56] mô phỏng biến dạng, tách hạt và tương tác va chạm của hạt đa thành phần trong dòng chảy cắt.

Khoảng trống nghiên cứu luận án giải quyết

Trong khi các nghiên cứu số về hạt đơn hoặc tương tác hạt đa thành phần đẳng nhiệt đã phát triển mạnh, chưa có công trình mô phỏng số nào khảo sát chi tiết quá trình truyền nhiệt kết hợp chuyển pha hóa rắn của hạt lưu chất đa thành phần rỗng. Luận án đã lấp đầy khoảng trống này bằng việc xây dựng mô hình tính toán số chuyên biệt bằng phương pháp theo dấu biên.

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

Hệ phương trình toán học điều khiển

Quá trình truyền nhiệt và động lực học dòng chảy nhiều pha có chuyển pha được mô tả thông qua hệ phương trình Navier-Stokes tổng quát, phương trình liên tục có số hạng nguồn chuyển pha và phương trình bảo toàn năng lượng:

  1. Phương trình Navier-Stokes dạng liên tục một trường (One-field formulation): $$\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 \sigma \kappa \mathbf{n}_f \delta(\mathbf{x} - \mathbf{x}_f) dS + \mathbf{f} + \rho \mathbf{g}$$

Trong đó: $\mathbf{u} = (u, v)$ là véc-tơ vận tốc, $p$ là áp suất, $\rho$ là khối lượng riêng, $\mu$ là độ nhớt động lực, $\sigma$ là hệ số sức căng bề mặt, $\kappa$ là độ cong trung bình của bề mặt phân cách, $\mathbf{n}_f$ là véc-tơ pháp tuyến đơn vị tại biên phân cách $\mathbf{x}_f$, $\delta$ là hàm Dirac delta, $\mathbf{g}$ là gia tốc trọng trường, và $\mathbf{f}$ là lực thể tích áp đặt điều kiện biên không trượt tại pha rắn ($\mathbf{u} = 0$).

Trong hệ tọa độ trụ đối xứng ($Orz$), thành phần vận tốc theo phương hướng tâm $r$ và phương trục $z$ được mở rộng chi tiết với các số hạng khuếch tán và đối lưu tương ứng.

  1. Phương trình bảo toàn năng lượng: $$\frac{\partial (\rho C_p T)}{\partial t} + \nabla \cdot (\rho C_p T \mathbf{u}) = \nabla \cdot (k \nabla T) + \int q \delta(\mathbf{x} - \mathbf{x}_f) dS$$

Với $T$ là nhiệt độ, $C_p$ là nhiệt dung riêng đẳng áp, $k$ là hệ số dẫn nhiệt. Thông lượng nhiệt $q$ tại biên chuyển pha được xác định bởi định luật Fourier: $$q = k_s \left. \frac{\partial T}{\partial n} \right|_s - k_l \left. \frac{\partial T}{\partial n} \right|_l$$

  1. Phương trình liên tục khi có biến thiên khối lượng riêng do chuyển pha: $$\nabla \cdot \mathbf{u} = \frac{1}{L_h} \left( \frac{1}{\rho_s} - \frac{1}{\rho_l} \right) q \delta(\mathbf{x} - \mathbf{x}_f)$$ Trong đó $L_h$ là nhiệt ẩn chuyển pha, chỉ số $s, l, g$ lần lượt ký hiệu cho pha rắn, pha lỏng và pha khí.

Chuẩn hóa không thứ nguyên

Hệ phương trình được chuyển về dạng không thứ nguyên thông qua các đại lượng tham chiếu ($L, U_t, t_t = L/U_t, \sigma_o, T_t$), dẫn đến các chuẩn số chi phối:

Ký hiệu Tên gọi chuẩn số / Thông số Biểu thức định nghĩa Ý nghĩa vật lý
$Re$ Số Reynolds $\frac{\rho_l L U_t}{\mu_l}$ Tỉ số giữa lực quán tính và lực nhớt
$We$ Số Weber $\frac{\rho_l L U_t^2}{\sigma_o}$ Tỉ số giữa lực quán tính và lực căng bề mặt
$Pr$ Số Prandtl $\frac{C_{pl} \mu_l}{k_l}$ Tỉ số giữa khuếch tán động lượng và khuếch tán nhiệt
$St$ Số Stefan $\frac{C_{pl} (T_m - T_c)}{L_h}$ Tỉ số giữa nhiệt cảm và nhiệt ẩn chuyển pha
$Bo$ Số Bond $\frac{\rho_l g L^2}{\sigma_o}$ Tỉ số giữa trọng lực và lực căng bề mặt
$Oh$ Số Ohnesorge $\frac{\mu_l}{\sqrt{\rho_l \sigma_o L}}$ Mối liên hệ lực nhớt với lực quán tính và căng bề mặt
$\rho_{sl}, \rho_{gl}$ Tỉ số khối lượng riêng $\frac{\rho_s}{\rho_l}, \frac{\rho_g}{\rho_l}$ Sai khác mật độ giữa các pha
$R_{io}$ Tỉ số bán kính $\frac{R_i}{R_o}$ Tỉ lệ kích thước nhân khí so với hạt ngoài
$\theta$ Nhiệt độ không thứ nguyên $\frac{T - T_c}{T_m - T_c}$ Phân bố trường nhiệt độ chuẩn hóa

Phương pháp số và kỹ thuật tính toán

  • Môi trường lập trình: Toàn bộ chương trình tính toán được lập trình bằng ngôn ngữ Fortran.
  • Lưới tính toán: Phương pháp sai phân hữu hạn trên lưới Euler so le cố định (Marker-and-Cell - MAC). Vận tốc được bố trí tại các mặt ô lưới, áp suất và nhiệt độ được đặt tại tâm ô lưới.
  • Phương pháp theo dấu biên (Front-Tracking Method):
    • Các biên phân cách (lỏng-khí, rắn-lỏng, rắn-khí) được rời rạc hóa thành các chuỗi điểm liên kết di động $\mathbf{x}_f$.
    • Tọa độ điểm biên được cập nhật theo thời gian: $\mathbf{x}_f^{n+1} = \mathbf{x}_f^n + \mathbf{V}_f \Delta t$. Vận tốc pháp tuyến $\mathbf{V}_f$ được nội suy từ lưới Euler hoặc xác định từ điều kiện cân bằng nhiệt tại mặt chuyển pha.
    • Sử dụng hai hàm chỉ thị thể tích $I_1, I_2$ (hoặc $I_s, I_l$) phân bố từ 0 đến 1 để phân định trạng thái 3 pha và nội suy liên tục các thuộc tính vật lý ($\rho, \mu, C_p, k$) trên toàn miền tính toán.

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ơ sở lý thuyết về sự hình thành hạt lưu chất đa thành phần (qua vòi phun đồng trục, kênh giao cắt và sự co lại của sợi lưu chất), phân tích các ứng dụng thực tế trong công nghệ pin mặt trời, bảo vệ cánh tuabin/máy bay, y sinh và thực phẩm. Đồng thời, chương này so sánh ưu nhược điểm của các phương pháp số (Level-Set, VOF, Lattice Boltzmann, Front-Tracking), từ đó khẳng định tính ưu việt của phương pháp Front-Tracking trong việc xác định chính xác lực căng bề mặt và quản lý biên chuyển pha phức tạp.

Chương 2: Nghiên cứu và xây dựng chương trình tính toán cho bài toán có truyền nhiệt, chuyển pha

Nội dung chương tập trung xây dựng thuật toán số giải hệ phương trình Navier-Stokes và phương trình năng lượng trên lưới so le MAC.

Kiểm chứng chương trình tính toán: Để đánh giá độ tin cậy của thuật toán, tác giả thực hiện đối chứng kết quả mô phỏng với dữ liệu thực nghiệm đã công bố:

  1. Hạt nước có góc ướt nhỏ hơn $90^\circ$: Mô phỏng hạt nước với góc ướt $\theta_0 = 78^\circ$ và $\theta_0 = 85^\circ$. Kết quả hình thái biên hóa rắn và sự biến thiên thể tích hạt theo thời gian cho thấy sự trùng khớp với kết quả thực nghiệm của Zhang và cộng sự [69].
  2. Hạt nước có góc ướt lớn hơn $90^\circ$: Mô phỏng hạt nước kỵ nước với $\theta_0 = 124^\circ$ và siêu kỵ nước với $\theta_0 = 155^\circ$. Kết quả thể hiện chính xác động học của trường nhiệt độ $\theta$, trường vận tốc không thứ nguyên $U_t$ và quá trình co kéo hình thành đỉnh chóp nhọn.

Chương 3: Nghiên cứu quá trình truyền nhiệt, chuyển pha của hạt lưu chất rỗng trên một bề mặt lạnh

Chương này khảo sát quá trình đóng băng của hạt lưu chất rỗng đối xứng đặt trên bề mặt phẳng lạnh duy trì ở nhiệt độ $T_c$.

Nghiên cứu độ hội tụ lưới: Thực hiện đánh giá sai số giữa các độ phân giải lưới $128 \times 256$ và $256 \times 512$, chứng minh sai số trung bình nằm trong giới hạn cho phép và đảm bảo tính độc lập lưới.

Kết quả phân tích ảnh hưởng của các thông số:

  • Số Bond ($Bo = 0,18$ và $Bo = 1,78$): Số $Bo$ lớn làm tăng tác động của trọng lực, khiến hạt bị bẹp hơn, làm tăng diện tích tiếp xúc với bề mặt lạnh, dẫn đến chiều cao hạt sau hóa rắn ($H_{out}, H_{in}$) giảm và thời gian hóa rắn thay đổi.
  • Số Prandtl ($Pr = 0,032$ và $Pr = 0,316$): $Pr$ điều khiển tỉ lệ giữa khuếch tán động lượng và nhiệt, tác động trực tiếp đến gradien nhiệt độ trong pha lỏng và hình thái mặt phân cách.
  • Số Stefan ($St = 0,032$ và $St = 0,316$): Khi $St$ tăng, nhiệt cảm chiếm ưu thế so với nhiệt ẩn, làm tăng tốc độ dịch chuyển của biên hóa rắn $V_n$ và rút ngắn đáng kể tổng thời gian hóa rắn $\tau_s$.
  • Tỉ số khối lượng riêng ($\rho_{sl} = 0,8$ đến $1,1$): Khi $\rho_{sl} < 1$ (hiện tượng dãn nở thể tích như nước đóng băng), đỉnh chóp nhọn xuất hiện rõ rệt tại đỉnh hạt ngoài do chất lỏng bị đẩy lên trong giai đoạn cuối. Ngược lại, khi $\rho_{sl} > 1$, đỉnh hạt bị lõm xuống.
  • Tỉ số bán kính ($R_{io} = 0,3$ đến $0,7$): Bóng khí bên trong càng lớn ($R_{io}$ tăng) thì thể tích chất lỏng cần chuyển pha càng ít, thời gian hóa rắn kết thúc càng nhanh. Quá trình hóa rắn xuất hiện một điểm gấp khúc đặc trưng trên đường cong chiều cao biên hóa rắn ($H_s$) tại thời điểm mặt đông đá bao phủ hoàn toàn hạt bóng khí.
  • Ảnh hưởng của các góc tiếp xúc:
    • Khảo sát góc ướt trong $\theta_{0i}$ ($50^\circ - 120^\circ$), góc ướt ngoài $\theta_{0o}$ ($60^\circ - 130^\circ$) và trường hợp đồng thời $\theta_0 = 180^\circ - \theta_{0i} = \theta_{0o}$ ($60^\circ - 130^\circ$). Góc ướt quyết định diện tích bám dính ban đầu và độ cong mặt ngoài/mặt trong.
    • Góc phát triển $\phi_{gr}$ ($0^\circ - 25^\circ$): Chi phối trực tiếp góc nghiêng tiếp tuyến tại điểm chập ba pha, làm biến đổi mạnh mẽ độ tăng chiều cao đỉnh ngoài ($H_d - h_{0o}$) và bóng khí ($h_i - h_{0i}$).

Chương 4: Nghiên cứu quá trình truyền nhiệt, chuyển pha của hạt lưu chất rỗng chịu ảnh hưởng của đối lưu cưỡng bức

Chương 4 mở rộng mô hình cho bài toán hạt lưu chất rỗng lơ lửng trong không gian tự do chịu tác động của một dòng khí lạnh cưỡng bức thổi từ đáy lên với vận tốc đầu vào $U_{in}$ và nhiệt độ $T_0$.

Nghiên cứu độ hội tụ lưới: So sánh sai số giữa các cấu hình lưới khác nhau với lưới chuẩn phân giải cao $384 \times 1536$.

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

  • Số Reynolds ($Re$): Tăng $Re$ làm gia tăng cường độ đối lưu cưỡng bức xung quanh hạt, màng nhiệt bị ép mỏng ở phía đón dòng (cực dưới), làm tăng tốc độ hóa rắn ở phần đáy và gây bất đối xứng mạnh về hình thái.
  • Tỉ số hình dạng trong ($A_{ri}$) và ngoài ($A_{ro}$): Sự phân bố áp suất và lực ma sát của dòng ngoài làm hạt bị biến dạng cơ học trước và trong khi đông đặc. Luận án định lượng sự thay đổi của $A_{ri}, A_{ro}$ và thời gian hóa rắn theo dải biến thiên của $Re$.
  • Ảnh hưởng của các thông số khác:
    • Khảo sát ảnh hưởng của số Stefan ($St$), tỉ số khối lượng riêng ($\rho_{sl}$), kích thước nhân hóa rắn ban đầu ($r_0/R$), độ lệch tâm ban đầu của bóng khí bên trong ($\varepsilon_0$), tỉ số bán kính ($R_{io}$) và góc phát triển ($\phi_{gr}$).
    • Độ lệch tâm $\varepsilon_0$ phá vỡ tính đối xứng truyền nhiệt nội tại, dẫn đến việc đóng băng không đồng đều xung quanh chu vi bóng khí.

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

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

  1. Thiết lập mô hình số chuyên biệt: Xây dựng thành công chương trình tính toán mô phỏng 2D đối xứng trục dựa trên phương pháp theo dấu biên (Front-Tracking) kết hợp lưới so le Euler MAC để giải bài toán động lực học 3 pha (rắn - lỏng - khí) có truyền nhiệt và chuyển pha hóa rắn.
  2. Xử lý hình học đa biên và điểm chập ba pha kép: Mô hình hóa chính xác sự tồn tại đồng thời của hai điểm chập ba pha (tại mặt ngoài hạt lưu chất và tại mặt phân cách bóng khí bên trong), điều mà các mô hình giải tích và phương pháp mặt mức thông thường gặp khó khăn khi mô tả.
  3. Làm rõ quy luật động lực học chuyển pha: Phát hiện và giải thích bản chất vật lý của điểm gấp khúc trên đường cong chiều cao biên hóa rắn $H_s(t)$ – thời điểm đánh dấu sự kết thúc quá trình hóa rắn bao quanh hạt bóng khí bên trong.

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

  1. Cung cấp bộ cơ sở dữ liệu định lượng về sự phụ thuộc của thời gian hóa rắn và hình thái chóp nhọn vào các chuẩn số ($Bo, Pr, St, Re, \rho_{sl}, R_{io}, \phi_{gr}$).
  2. Cung cấp công cụ dự báo và điều khiển hình dạng hạt phục vụ tối ưu hóa quy trình chế tạo hạt bán dẫn rỗng trong sản xuất pin mặt trời màng mỏng.
  3. Cung cấp các thông số thủy nhiệt hỗ trợ nghiên cứu cơ chế bám dính và phát triển các giải pháp phá băng thụ động/chủ động trên cánh máy bay và cánh tuabin gió.

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

Hạn chế của đề tài

  • Mô hình tính toán hiện tại bị giới hạn trong hệ tọa độ hai chiều đối xứng trục ($Orz$), chưa mô phỏng được các dao động không đối xứng ba chiều (3D) hoàn chỉnh.
  • Giả thiết các pha là lưu chất Newton thuần túy với các hệ số nhiệt vật lý không phụ thuộc nhiệt độ.
  • Do giới hạn tài nguyên tính toán, số lượng tham số khảo sát mới tập trung vào các dải thông số đại diện chính.

Hướng phát triển trong tương lai

  • Mở rộng thuật toán Front-Tracking cho bài toán không gian 3 chiều đầy đủ nhằm mô tả hiện tượng lệch trục và biến dạng bất đối xứng phức tạp.
  • Phát triển mô hình cho các bài toán đa nhân khí/lỏng (nhiều hạt bóng khí phân tán bên trong) và lưu chất phi Newton.
  • Nghiên cứu quá trình va chạm động học của hạt lưu chất rỗng lên bề mặt kèm theo chuyển pha tức thời.

Giá trị tham khảo

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

  • Nghiên cứu sinh và giảng viên: Chuyên ngành Cơ học chất lỏng, Thủy khí động lực học tính toán (CFD), Kỹ thuật Nhiệt - Năng lượng, và Vật lý tính toán; đặc biệt là phương pháp số giải bài toán nhiều pha có biên chuyển pha (Stefan problem).
  • Kỹ sư và nhà khoa học vật liệu: Nghiên cứu công nghệ chế tạo hạt cầu rỗng kích thước micro/nano ứng dụng trong pin mặt trời, vật liệu hấp thụ âm thanh hoặc bao bọc vi nang trong dược - thực phẩm.
  • Kỹ sư hàng không và năng lượng gió: Nghiên cứu mô hình toán học và cơ chế vật lý của hiện tượng đóng băng khí quyển trên các bề mặt chuyển động.

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

1. Tại sao việc mô phỏng hạt bán dẫn rỗng lại có ý nghĩa quyết định trong công nghệ pin mặt trời?

Theo công trình của Biancardo và cộng sự [5] được trích dẫn trong luận án, đối với các hạt bán dẫn hình cầu Silic, vùng tham gia chuyển đổi quang điện thành điện năng chủ yếu tập trung ở lớp mỏng ngoài bề mặt. Do đó, việc chế tạo các hạt bán dẫn rỗng (chứa nhân khí bên trong) giúp tiết kiệm đáng kể lượng vật liệu bán dẫn đắt tiền và làm nhẹ tấm pin mà hiệu suất chuyển đổi vẫn được duy trì. Mô phỏng số giúp xác định các thông số nhiệt động lực học để kiểm soát độ dày vỏ và độ đồng đều hình học của hạt trong quá trình hóa rắn.

2. Ưu thế chính của phương pháp theo dấu biên (Front-Tracking) so với phương pháp Volume of Fluid (VOF) và Level-Set trong bài toán này là gì?

Phương pháp Front-Tracking sử dụng một chuỗi các điểm liên kết di động rời rạc trên nền lưới Euler cố định để biểu diễn biên phân cách. Kỹ thuật này cho phép:

  • Xác định trực tiếp và chính xác độ cong $\kappa$ và lực căng bề mặt $\sigma$, khắc phục nhược điểm tính toán lực căng bề mặt phức tạp và dễ sai số của VOF và Level-Set.
  • Theo dõi sắc nét vị trí mặt chuyển pha mà không bị khuếch tán số (numerical diffusion) hay phụ thuộc quá mức vào hàm Heaviside.
  • Dễ dàng định vị và xử lý điều kiện biên động tại các điểm chập ba pha giữa các pha rắn, lỏng và khí.

3. Nguyên nhân vật lý dẫn đến sự hình thành đỉnh chóp nhọn tại mặt ngoài của hạt sau khi hóa rắn là gì?

Sự hình thành chóp nhọn bắt nguồn từ hai yếu tố chính:

  • Sự sai khác khối lượng riêng ($\rho_{sl} < 1$): Khi nước đóng băng, thể tích pha rắn dãn nở lớn hơn pha lỏng. Trong quá trình hóa rắn từ đáy lên trên, chất lỏng chưa đông bị ép và dồn dần lên đỉnh hạt.
  • Động học góc tiếp xúc và góc phát triển ($\phi_{gr}$): Ở giai đoạn cuối cùng của quá trình đông đặc, góc tiếp xúc giữa mặt lỏng-khí và mặt phẳng ngang biến thiên đột ngột, kéo mặt phân cách co lại thành một đỉnh nhọn trước khi chuyển pha hoàn tất.

4. Điểm gấp khúc trên đường cong chiều cao biên hóa rắn $H_s(t)$ tại Chương 3 thể hiện hiện tượng vật lý gì?

Điểm gấp khúc trên đồ thị $H_s(t)$ biểu thị thời điểm biên hóa rắn di chuyển từ đáy lên đã hoàn thành việc bao bọc toàn bộ bề mặt của hạt bóng khí bên trong. Sau thời điểm này, diện tích truyền nhiệt và cấu trúc dòng nhiệt chuyển từ dạng truyền qua lớp vỏ quanh nhân khí sang dạng đông đặc khối chất lỏng phía trên đỉnh, làm thay đổi đột ngột tốc độ nâng cao của mặt phân cách pha.

Kết luận

Luận án đã giải quyết bài toán mô phỏng số động lực học truyền nhiệt và chuyển pha hóa rắn của hạt lưu chất rỗng đa thành phần bằng phương pháp theo dấu biên (Front-Tracking) trên nền lưới so le MAC. Công trình đã xây dựng thành công hệ mã nguồn Fortran tin cậy, được kiểm chứng chính xác qua các dữ liệu thực nghiệm về hạt nước trên bề mặt phẳng lạnh. Các kết quả mô phỏng đã làm rõ quy luật chi phối của các chuẩn số thủy nhiệt ($Bo, Pr, St, Re, \rho_{sl}, R_{io}, \phi_{gr}$) đến động học đông đặc, sự biến dạng bóng khí và cơ chế hình thành chóp nhọn. Công trình cung cấp cơ sở khoa học định lượng phục vụ các ứng dụng công nghệ trong chế tạo pin mặt trời màng mỏng và kiểm soát hiện tượng đóng băng trong kỹ thuật hàng không.