Luận á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

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, mô phỏng quá trình nhiệt thủy động lực học chuyên sâu và chính xác.

• Trường học: Học viện Khoa học và Công nghệ
• Chuyên ngành: Cơ học chất lỏng và chất khí
• Đề 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
• Loại tài liệu: Luận án tiến sĩ
• Năm xuất bản: 2023
• Địa điểm: Hà Nội

Tóm tắt

I. Tổng quan về nghiên cứu động lực học hạt lưu chất đa thành phần

Nghiên cứu động lực học hạt lưu chất đa thành phần có truyền nhiệt chuyển pha đóng vai trò quan trọng trong nhiều lĩnh vực công nghiệp như hóa chất, năng lượng và y sinh. Hạt lưu chất đa thành phần hình thành khi các dòng lưu chất khác nhau tương tác, tạo ra cấu trúc đa lớp hoặc đa nhân. Quá trình này liên quan đến cơ chế tách hạt, co thắt dòng chảy và tương tác giữa các pha. Các nghiên cứu thực nghiệm và lý thuyết đã chứng minh khả năng kiểm soát kích thước, hình dạng hạt thông qua thiết kế kênh dòng chảy. Ứng dụng thực tiễn bao gồm sản xuất vật liệu nano, dược phẩm và công nghệ vi lưu. Sự phát triển của phương pháp mô phỏng số đã nâng cao hiểu biết về cơ chế hình thành hạt, đặc biệt trong điều kiện truyền nhiệt và chuyển pha.

1.1. Cơ chế hình thành hạt lưu chất đa thành phần

Hạt lưu chất đa thành phần được tạo ra thông qua các phương pháp vật lý như sử dụng ống co thắt, vòi phun đồng trục hoặc kênh giao nhau. Khi lưu chất đi qua vùng co thắt, áp suất giảm đột ngột gây tách dòng thành các hạt rời rạc. Cơ chế này phụ thuộc vào vận tốc dòng chảy, độ nhớt và sức căng bề mặt. Các nghiên cứu cho thấy kích thước hạt có thể điều chỉnh bằng cách thay đổi lưu lượng hoặc thiết kế hình học của kênh. Hạt đa nhân hình thành khi dòng lưu chất chính bị cắt bởi dòng phụ, tạo ra cấu trúc lớp hoặc lõi-vỏ.

1.2. Ứng dụng trong công nghiệp và khoa học

Ứng dụng chính của hạt lưu chất đa thành phần bao gồm sản xuất vật liệu composite, thuốc nano và hệ thống dẫn thuốc. Trong y sinh, hạt đa lớp được sử dụng để vận chuyển thuốc đích đến tế bào ung thư. Công nghệ vi lưu tận dụng hạt đa thành phần để kiểm soát chính xác phản ứng hóa học. Trong năng lượng, hạt lưu chất đa thành phần cải thiện hiệu suất truyền nhiệt trong các hệ thống nhiệt điện. Các nghiên cứu gần đây tập trung vào tối ưu hóa quá trình hình thành hạt nhằm nâng cao hiệu suất sản xuất.

II. Phân tích các thách thức trong mô hình hóa hạt lưu chất đa pha

Mô hình hóa động lực học hạt lưu chất đa thành phần gặp nhiều thách thức do sự phức tạp của tương tác giữa các pha và điều kiện biên. Các vấn đề chính bao gồm biến dạng bề mặt tự do, truyền nhiệt không ổn định và chuyển pha phức tạp. Phương trình Navier-Stokes kết hợp với điều kiện biên di động đòi hỏi giải pháp số phức tạp. Sự không ổn định của dòng chảy gây khó khăn trong việc dự đoán kích thước và hình dạng hạt. Ngoài ra, tính chất vật liệu thay đổi theo nhiệt độ và áp suất, ảnh hưởng đến độ chính xác của mô hình. Các nghiên cứu gần đây tập trung vào phát triển phương pháp theo dấu biên và lưới thích ứng để cải thiện độ chính xác.

2.1. Ảnh hưởng của truyền nhiệt đến quá trình hình thành hạt

Truyền nhiệt ảnh hưởng đáng kể đến sự hình thành và ổn định hạt lưu chất đa thành phần. Gradient nhiệt độ gây ra sự thay đổi sức căng bề mặt, dẫn đến biến dạng hạt. Trong quá trình đông đặc, nhiệt độ đóng vai trò quan trọng trong việc kiểm soát tốc độ chuyển pha. Các nghiên cứu cho thấy sự không ổn định nhiệt có thể gây ra hiện tượng tách lớp không mong muốn. Phương pháp mô phỏng cần tích hợp đầy đủ các hiệu ứng nhiệt để dự đoán chính xác cấu trúc hạt.

2.2. Thách thức trong mô phỏng chuyển pha

Chuyển pha từ lỏng sang rắn hoặc khí tạo ra biên phân cách phức tạp, đòi hỏi xử lý đặc biệt trong mô phỏng. Sự không liên tục của các tính chất vật liệu tại biên chuyển pha gây khó khăn cho giải pháp số. Phương pháp theo dấu biên được sử dụng rộng rãi để theo dõi chuyển động của biên. Tuy nhiên, việc cập nhật biên theo thời gian đòi hỏi tính toán lặp lại, tăng chi phí tính toán. Các nghiên cứu gần đây tập trung vào phát triển thuật toán hiệu quả để giảm thời gian tính toán.

III. Phương pháp tính toán động lực học hạt lưu chất đa thành phần

Phương pháp tính toán động lực học hạt lưu chất đa thành phần kết hợp phương trình Navier-Stokes, truyền nhiệt và mô hình chuyển pha. Phương pháp theo dấu biên được sử dụng để theo dõi chuyển động của biên phân cách giữa các pha. Lưới thích ứng được áp dụng để tăng độ phân giải tại vùng biên. Các điều kiện biên bao gồm vận tốc, nhiệt độ và cân bằng năng lượng tại biên chuyển pha. Phương pháp phần tử hữu hạn kết hợp với kỹ thuật theo dấu biên cho phép mô phỏng chính xác quá trình hình thành hạt. Các nghiên cứu gần đây tập trung vào tối ưu hóa thuật toán để cải thiện hiệu suất tính toán.

3.1. Mô hình toán học và điều kiện biên

Mô hình toán học bao gồm phương trình bảo toàn khối lượng, động lượng và năng lượng. Điều kiện biên bao gồm vận tốc không trượt tại thành ống, nhiệt độ không đổi tại thành và cân bằng năng lượng tại biên chuyển pha. Phương pháp theo dấu biên sử dụng các điểm biên di động để theo dõi chuyển động của biên phân cách. Vận tốc biên được xác định từ điều kiện cân bằng năng lượng hoặc nội suy từ lưới tính toán. Phương pháp này cho phép mô phỏng chính xác sự tương tác giữa các pha.

3.2. Giải pháp số và tối ưu hóa thuật toán

Giải pháp số sử dụng phương pháp phần tử hữu hạn kết hợp với lưới thích ứng. Thuật toán giải quyết các phương trình phi tuyến bằng phương pháp lặp Newton-Raphson. Tối ưu hóa thuật toán bao gồm sử dụng kỹ thuật chia sẻ tải tính toán và song song hóa. Các nghiên cứu gần đây tập trung vào phát triển thuật toán hiệu quả để giảm thời gian tính toán. Phương pháp này cho phép mô phỏng các quá trình phức tạp trong thời gian hợp lý.

IV. Kết luận và triển vọng ứng dụng trong thực tế

Nghiên cứu đã cung cấp cái nhìn toàn diện về độ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ác phương pháp tính toán tiên tiến đã được phát triển để mô phỏng chính xác quá trình hình thành hạt. Kết quả nghiên cứu góp phần nâng cao hiểu biết về cơ chế tương tác giữa các pha. Ứng dụng thực tế bao gồm sản xuất vật liệu nano, dược phẩm và hệ thống nhiệt điện. Các nghiên cứu tiếp theo sẽ tập trung vào tối ưu hóa quá trình hình thành hạt và mở rộng ứng dụng trong các lĩnh vực công nghiệp khác.

4.1. Đánh giá kết quả nghiên cứu và đóng góp khoa học

Nghiên cứu đã giải quyết thành công các thách thức trong mô phỏng động lực học hạt lưu chất đa thành phần. Phương pháp theo dấu biên kết hợp lưới thích ứng đã được chứng minh là hiệu quả trong việc mô phỏng quá trình hình thành hạt. Các kết quả thu được phù hợp với dữ liệu thực nghiệm, xác nhận độ tin cậy của phương pháp. Đóng góp chính bao gồm phát triển mô hình toán học và thuật toán tính toán tiên tiến.

4.2. Hướng phát triển và ứng dụng trong tương lai

Các nghiên cứu tiếp theo sẽ tập trung vào tối ưu hóa quá trình hình thành hạt bằng cách điều chỉnh các thông số vận hành. Ứng dụng trong sản xuất vật liệu nano và dược phẩm sẽ được mở rộng. Ngoài ra, tích hợp trí tuệ nhân tạo vào quá trình mô phỏng sẽ giúp tối ưu hóa thiết kế hệ thống. Các nghiên cứu tương lai cũng sẽ tập trung vào ứng dụng trong lĩnh vực năng lượng tái tạo và y sinh.

Trích đoạn nội dung tài liệu

MỞ ĐẦU Tính cấp thiết của đề tài Các lưu chất (chất lỏng hoặc chất khí) khi nhúng vào lưu chất khác mà bề mặt các lưu chất này có xuất hiện sức căng bề mặt sẽ tạo thành các hạt lưu chất. Các hạt lưu chất này có thể xảy ra quá trình hóa rắn trong tự nhiên hoặc trong các quá trình công nghiệp. Ngày nay, quá trình chuyển pha hạt lưu chất như vậy được nhiều nhà nghiên cứu quan tâm và được ứng dụng rộng rãi [1]–[7]. Trong đó, hạt lưu chất đa thành phần (sau đây gọi tắt là hạt đa thành phần) có tiềm năng ứng dụng rất lớn trong việc chế tạo pin năng lượng mặt trời.

Có rất nhiều nhà nghiên cứu quan tâm đến nguồn năng lượng có thể tái tạo này vì chúng gần như vô tận và ít gây ra những tác hại tiêu cực đến môi trường như đối với nguồn nhiên liệu hóa thạch. Biancardo và cộng sự [5] đã nghiên cứu các hạt bán dẫn để biến đổi năng lượng ánh sáng mặt trời thành năng lượng điện. Các hạt Silic cầu được đặt trong một tấm lưới. Theo đó, giá trị sử dụng của hạt bán dẫn chỉ chủ yếu là một lớp mỏng bên ngoài hạt.

Do đó, việc sử dụng hạt rỗng (hay hạt lưu chất đa thành phần hóa rắn) sẽ giúp tiết kiệm chất bán dẫn hơn là làm cho tấm lưới bán dẫn được nhẹ hơn mà hiệu quả chuyển đổi năng lượng là không đổi. Một số các công trình khác có nghiên cứu việc chế tạo pin mặt trời có thể tham khảo tại [4], [8], [9]. Việc nghiên cứu sự chuyển pha của các hạt đa thành phần bán dẫn như vậy có thể giúp việc chế tạo những tấm pin mặt trời đại trà, rẻ hơn và rộng rãi hơn. Việc nghiên cứu chuyển pha (trong luận án này là hóa rắn) hạt lưu chất đa thành phần còn giúp tìm ra giải pháp phá băng đá trên bề mặt cánh máy bay [10] và tuabin gió [11], [12] nhằm nâng cao tuổi thọ và hiệu suất máy.

Theo đó, hạt nước trong không khí khi tiếp xúc với bề mặt cánh của máy bay hoặc tuabin gió sẽ bị hóa rắn nếu nhiệt độ tại bề mặt cánh nhỏ hơn nhiệt độ hóa rắn của nước. Các hạt nước như vậy có thể lẫn các bóng khí bên trong và chúng cũng bị hóa rắn khi bám vào các bề mặt cánh. Sự hóa rắn của hạt nước trên bề mặt cánh có thể là một nguyên nhân nghiêm trọng giảm hiệu suất và tuổi thọ của máy. Nguy hiểm hơn khi nó là một trong những nguyên nhân trong các vụ tai nạn hàng không nghiêm trọng do ảnh hưởng tới khí động lực học của dòng khí khi qua cánh máy bay [10], [13], [14].

Do đó, việc ứng dụng nghiên cứu sự truyền nhiệt, chuyển pha của hạt lưu chất nói 2 chung và hạt đa thành phần nói riêng góp phần tìm ra giải pháp, tránh các thiệt hại và nâng cao hiệu suất của máy. Ngoài ra, hạt lưu chất đa thành phần còn được sử dụng để loại bỏ các hạt nước bẩn [7], trong công nghệ thực phẩm [2], [3] và sản xuất vật liệu hấp thụ âm thanh [6]. Việc nghiên cứu hóa rắn hạt đa thành phần có truyền nhiệt, chuyển pha có tiềm năng ứng dụng rất lớn đặc biệt trong việc chế tạo pin năng lượng mặt trời và loại bỏ băng đá trên bề mặt cánh máy bay, tuabin. Việc hiểu thấu đáo các đặc tính vật lý sâu bên trong của động lực học truyền nhiệt, chuyển pha của hạt lưu chất đa thành phần (điều mà khó có thể đạt được bằng lý thuyết hay thực nghiệm) là cần thiết.

Đây cũng là mục đích chính cho sự ra đời của đề tài luận án tiến sĩ “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”. Mục tiêu của đề tài Mục tiêu tổng quát của đề tài: nghiên cứu quá trình truyền nhiệt, chuyển pha (ở đây là quá trình hóa rắn) của hạt lưu chất rỗng (một loại hạt lưu chất đa thành phần có nhân khí bên trong) bằng phương pháp theo dấu biên (front-tracking method). Mục tiêu cụ thể của đề tài: xây dựng mô hình tính toán hạt lưu chất rỗng hóa rắn và khảo sát ảnh hưởng của các thông số không thứ nguyên đến quá trình hóa rắn hạt lưu chất rỗng. Ý nghĩa khoa học và thực tiễn của đề tài Ý nghĩa khoa học của đề tài: thiết lập một chương trình tính toán mô phỏng số để mô tả các quá trình vật lý động lực học chuyển pha của hạt lưu chất rỗng.

Các thông số không thứ nguyên được xây dựng để đặc trưng cho các tính chất của lưu chất nghiên cứu và qua đó cho thấy sự ảnh hưởng của chúng đến quá trình chuyển pha của hạt. Ý nghĩa thực tiễn của đề tài: thông qua việc khảo sát ảnh hưởng của các thông số không thứ nguyên ảnh hưởng đến quá trình hóa rắn và hình dạng của hạt lưu chất rỗng, kết quả nghiên cứu cứu có thể ứng dụng trong quá trình chế tạo pin 3 mặt trời và giúp các nhà nghiên cứu ứng dụng trong việc đối phó với quá trình hóa rắn của hạt chất lỏng trên bề mặt cánh,… Đối tượng và phạm vi nghiên cứu của đề tài Đối tượng nghiên cứu của đề tài: các hạt lưu chất rỗng có kích thước nhỏ từ vài trăm micromet đến vài milimet. Chúng được hóa rắn khi tiếp xúc với bề mặt lạnh hoặc ở trong một môi trường lạnh. Phạm vi nghiên cứu của đề tài: tập trung vào các lưu chất không nén được, không trộn lẫn và là lưu chất Newton.

Phương trình Navier-Stokes cùng phương trình năng lượng được sử dụng. Các bài toán trong luận án này được giải quyết trong hệ tọa độ trụ. Các thông số không thứ nguyên được xem xét và chúng được khảo sát qua các dải giá trị khác nhau. Do hạn chế về mặt thời gian nên chỉ một số thông số không thứ nguyên và thông số chính được lựa chọn và khảo sát.

Phương pháp nghiên cứu Ngôn ngữ lập trình Fortran được sử dụng để xây dựng chương trình tính toán. Phương pháp sai phân hữu hạn được sử dụng để xấp xỉ các phương trình tính toán trên một nền lưới cố định (lưới Euler). Thêm nữa, do có sự xuất hiện của các lưu chất khác nhau nên xuất hiện các biên phân cách giữa các pha. Trong luận án này, phương pháp theo dấu biên (front-tracking method) được sử dụng để theo dõi sự di chuyển các biên phân cách của các pha đó.

Cụ thể, các biên phân cách giữa các pha được đại diện bởi các điểm liên kết và các điểm liên kết này di chuyển trên nền lưới cố định. Điểm mới của luận án - Đã xây dựng được một mô hình tính toán của hạt lưu chất rỗng hóa rắn đặt trên một bề mặt lạnh và hạt lưu chất rỗng lơ lửng hóa rắn trong môi trường tự do dưới tác động của đối lưu cưỡng bức. Cụ thể, hạt bóng khí hạt lưu chất được thêm vào bên trong hạt lưu chất tạo thành mô hình một hạt lưu chất rỗng. Từ đó, xuất hiện thêm một điểm chập ba pha ở biên phía bên trong của hạt; - Đã phân tích ảnh hưởng của các thông số không thứ nguyên đến quá trình hóa rắn của hạt lưu chất rỗng trên bề mặt lạnh.

Qua đó, có thể điều chỉnh 4 được hình dạng và thời gian hóa rắn của hạt lưu chất rỗng trên bề mặt lạnh đó bằng cách điều chỉnh các thông số không thứ nguyên; - Đã phân tích ảnh hưởng của các thông số không thứ nguyên đến quá trình hóa rắn hạt lưu chất rỗng lơ lửng trong môi trường tự do dưới tác động của đối lưu cưỡng bức. Điều chỉnh các thông số không thứ nguyên này sẽ gây ra thay đổi về hình dạng và thời gian hóa rắn của hạt lưu chất rỗng lơ lửng trong môi trường tự do dưới tác động của đối lưu cưỡng bức. Cấu trúc của luận án Nội dung chính của luận án bao gồm: Chương 1: Tổng quan Trình bày tổng quan về động lực học hạt đa thành phần bao gồm: một số khái niệm cơ bản, sự hình thành, các ứng dụng, các công trình nghiên cứu và các phương pháp phổ biến dùng trong mô phỏng động lực học chuyển pha của của hạt đa thành phần. 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 Nghiên cứu xây dựng chương trình tính toán cho bài toán hạt lưu chất hóa rắn.

Các phương trình tổng quát Navier-Stokes và phương trình năng lượng được rời rạc hóa và đưa vào các mô đun tính toán. Sau khi xây dựng được chương trình tính toán, để kiểm chứng độ chính xác, kết quả mô phỏng được so sánh với các kết quả thực nghiệm của các nhà nghiên cứu. 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 giải quyết bài toán hạt lưu chất rỗng hóa rắn trên một bề mặt lạnh. Một loạt các thông số không thứ nguyên và các thông số hình học được xem xét trong chương này.

Qua đó, cho thấy được sự ảnh hưởng của chúng lên hình dạng và thời gian hóa rắn của hạt lưu chất rỗng. 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 5 Chương này xem xét bài toán hạt lưu chất rỗng lơ lửng hóa rắn trong môi trường lạnh dưới sự tác động của một dòng cưỡng bức. Các thông số không thứ nguyên chính và góc phát triển cũng được xem xét để khảo sát ảnh hưởng của chúng lên hình dạng và thời gian hóa rắn của hạt. Cuối cùng là phần kết luận đưa ra tóm tắt các kết quả đã đạt được và những vấn đề cần giải quyết trong tương lai.

Một số khái niệm, định nghĩa cơ bản Sau đây là một số khái niệm, định nghĩa cơ bản các kí hiệu, thuật ngữ được dùng trong luận án. - Hạt lưu chất là hạt được tạo thành bởi một hoặc nhiều lưu chất (ở đây là chất lỏng, chất khí) không trộn lẫn và được nhúng trong một lưu chất khác. Với sự có mặt của sức căng bề mặt trên bề mặt lưu chất, các hạt lưu chất được hình thành. - Hạt lưu chất đa thành phần là hạt lưu chất bên trong có chứa các hạt lưu chất con khác (xem Hình 1.

- Hạt lưu chất rỗng là hạt lưu chất đa thành phần mà có một hạt bóng khí bên trong. - Hiện tượng chuyển pha của hạt lưu chất (ở đây là hiện tượng hóa rắn) là hiện tượng hạt lưu chất chuyển từ pha lỏng sang pha rắn khi hạt lưu chất trong một môi trường lạnh hoặc hạt tiếp xúc với các tác nhân khác có nhiệt độ nhỏ hơn nhiệt độ hóa rắn của lưu chất.