Nghiên Cứu Tính Toán Quá Trình Hóa Rắn Của Hạt Chất Lỏng

Nghiên cứu tính toán số quá trình hóa rắn của hạt nhân chất lỏng, khám phá các yếu tố ảnh hưởng và ứng dụng trong khoa học vật liệu.

• Tác giả: Phạm Duy Bính
• Người hướng dẫn: TS. Vũ Văn Trường
• Trường học: Đại học Bách Khoa Hà Nội
• Chuyên ngành: Kỹ thuật máy thay khí
• Đề tài: Nghiên Cứu Tính Toán Quá Trình Hóa Rắn Của Hạt Chất Lỏng
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2019
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CAM ĐOAN
• LỜI NÓI ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN VÀ Ý NGHĨA CỦA ĐỀ TÀI
• 1.1. Hiện tượng hóa rắn hạt chất lỏng trong tự nhiên và trong công nghiệp
• 1.2. Tình hình nghiên cứu bài toán trên thế giới và trong nước
• 1.3. Các hạt chất lỏng nghiên cứu trong luận văn với tầm quan trọng của nghiên cứu
• 2. CHƯƠNG 2: CÁC PHƯƠNG TRÌNH CƠ BẢN VÀ PHƯƠNG PHÁP MÔ PHỎNG
• 2.1. Các phương trình tổng quát
• 2.2. Tích phân theo thời gian
• 2.3. Rối rắm hóa các phương trình tổng quát
• 2.4. Rối rắm các thành phần khuếch tán
• 2.5. Phương trình áp suất
• 2.6. Rối rắm hóa phương trình năng lượng
• 2.7. Phương pháp theo điều biên
• 2.8. Cấu trúc mặt phân cách
• 2.9. Chuyển động của mặt phân cách
• 2.10. Cấu trúc lưới mặt phân cách
• 2.11. Làm mịn và tái cấu trúc hàm chỉ thị
• 2.12. Tính toán sức căng bề mặt trong lưới
• 2.13. Tính toán cân bằng nhiệt tại bề mặt chuyển pha
• 3. CHƯƠNG 3: KẾT QUẢ CỦA QUÁ TRÌNH MÔ PHỎNG VÀ THỬ LUẬN
• 3.1. Quá trình mô phỏng hóa rắn của một số hạt vật liệu chuyển pha
• 3.2. Đánh giá phương pháp mô phỏng
• 3.3. Phương pháp theo điều biên với hạt kép
• 3.4. Xây dựng biểu đồ pha với các góc tiếp xúc khác nhau của vật liệu silic
• 4. CHƯƠNG 4: KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN
• TÀI LIỆU THAM KHẢO
• CÁC CÔNG TRÌNH ĐÃ CÔNG BỐ

Tóm tắt

I. Tổng Quan Nghiên Cứu Tính Toán Quá Trình Hóa Rắn Hạt Lỏng

Hiện tượng hóa rắn hạt chất lỏng trên bề mặt lạnh xuất hiện phổ biến trong tự nhiên và công nghiệp. Việc hiểu rõ quá trình này có ý nghĩa quan trọng để phòng tránh tác hại và ứng dụng hiệu quả trong sản xuất. Hiện tượng này có thể gây hại hoặc ứng dụng trong sản xuất công nghiệp. Hiểu biết về quá trình này vô cùng quan trọng, góp phần phòng tránh tác hại và ứng dụng trong sản xuất, đời sống. Luận văn này ra đời để nghiên cứu tính toán số quá trình hóa rắn này. Với vai trò tác giả, tôi xin cảm ơn TS. Vũ Văn Trường, người hướng dẫn trực tiếp để tôi có thể hoàn thành luận văn này. Đồng thời tôi cũng xin gửi lời cảm ơn chân thành đến tập thể các thầy, cô giáo của Viện Cơ khí động lực, phòng Đào tạo Trường Đại học Bách Khoa Hà Nội đã tận tình chỉ bảo, truyền đạt những kiến thức mới cho tôi, giúp tôi suốt thời gian học tập và làm luận văn này.

1.1. Ứng Dụng Tầm Quan Trọng Nghiên Cứu Hóa Rắn Hạt Lỏng

Nghiên cứu này có tầm quan trọng trong việc giải quyết các vấn đề thực tiễn như băng đá hình thành trên cánh tuabin gió (ảnh hưởng đến hiệu suất và tuổi thọ), cánh máy bay (gây tai nạn). Nghiên cứu này còn tiềm năng ứng dụng trong công nghiệp, chẳng hạn như chế tạo pin năng lượng mặt trời, làm sạch bụi công nghiệp. Hiện tượng này còn được ứng dụng trong công nghiệp, hình 2 thể hiện những hạt bán dẫn trong chế tạo pin năng lượng mặt trời, ngành công nghiệp năng lượng tái tạo đang phát triển mạnh mẽ trong những năm gần đây. Các nghiên cứu về sự hóa rắn của hạt chất lỏng trên một bề mặt lạnh được quan tâm trong ngành công nghiệp này.

1.2. Mục Tiêu Nghiên Cứu Quá Trình Hóa Rắn Hạt Chất Lỏng

Mục tiêu chính của nghiên cứu là mô phỏng hạt chất lỏng trên bề mặt lạnh bằng phương pháp số. Từ đó, đưa ra mô hình mô phỏng sát thực tế, giúp hiểu rõ quá trình hóa rắn của các hạt chất lỏng trong tự nhiên và công nghiệp. Nghiên cứu tập trung vào việc mở rộng phương pháp mô phỏng số trực tiếp, chỉ ra ảnh hưởng của góc phát triển, góc tiếp xúc, sự thay đổi thể tích và đường chập ba pha lên quá trình hóa rắn. Ngoài ra, tính toán số giúp cho thấy ảnh hưởng của vật liệu khác nhau đến sản phẩm cuối cùng.

II. Thách Thức Vấn Đề Trong Nghiên Cứu Mô Phỏng Hóa Rắn

Mặc dù có nhiều nghiên cứu thực nghiệm, nghiên cứu mô phỏng quá trình hóa rắn trực tiếp vẫn còn hạn chế. Các phương pháp thường bỏ qua các yếu tố quan trọng như trọng lực, sự thay đổi thể tích, đường chập ba pha. Điều này dẫn đến kết quả tính toán không chính xác, không thể hiện được hình dạng thực tế của sản phẩm sau hóa rắn. Virozub và cộng sự [10] đã tính đến ảnh hưởng của trọng lực và sức căng bề mặt trong bài toán hóa rắn của hạt chất lỏng trên bề mặt lạnh. Tuy nhiên, các kết quả chưa cho thấy sự truyền nhiệt và trường dòng trong quá trình hạt chất lỏng hóa rắn.

2.1. Giới Hạn Của Các Nghiên Cứu Mô Phỏng Hóa Rắn Hiện Tại

Các nghiên cứu hiện tại thường đơn giản hóa quá trình, bỏ qua nhiều yếu tố ảnh hưởng. Điều này dẫn đến kết quả không chính xác và thiếu tính thực tế. Cụ thể, nhiều nghiên cứu bỏ qua trọng lực, sức căng bề mặt, sự thay đổi thể tích khi hóa rắn và ảnh hưởng của đường chập ba pha. Trong một nghiên cứu khác [11], phương trình dẫn nhiệt dựa trên entanpi được giải để cho thấy sự truyền nhiệt trong quá trình hóa rắn. Tuy nhiên, ảnh hưởng của đường chập ba pha cũng như sự thay đổi thể tích chưa được nghiên cứu, và do đó các kết quả tính toán không đưa ra được dạng hình côn tại đỉnh của sản phẩm sau quá trình hóa rắn.

2.2. Yếu Tố Ảnh Hưởng Đến Độ Chính Xác Mô Phỏng Quá Trình

Độ chính xác của mô phỏng quá trình hóa rắn phụ thuộc vào việc xem xét đầy đủ các yếu tố ảnh hưởng. Các yếu tố này bao gồm tính chất vật lý của vật liệu (nhiệt dung riêng, độ dẫn nhiệt, khối lượng riêng), điều kiện biên (nhiệt độ bề mặt, áp suất), và các hiện tượng vật lý (sức căng bề mặt, trọng lực, sự thay đổi thể tích). Mô hình hóa toán học quá trình hóa rắn là một nhiệm vụ phức tạp do sự thay đổi pha, liên quan đến sự giải phóng hoặc hấp thụ nhiệt tiềm ẩn và sự hình thành các ranh giới pha động. Sự có mặt của một bề mặt chất rắn, như chất nền, có thể ảnh hưởng đáng kể đến động học và hình thái của quá trình đông đặc.

III. Phương Pháp Tính Toán Số Mô Phỏng Quá Trình Hóa Rắn Hạt Lỏng

Luận văn sử dụng phương pháp theo dấu biên (Front Tracking) để mô phỏng quá trình. Đây là phương pháp mạnh mẽ để mô phỏng dòng nhiều pha, do GS. Grétar Tryggvason phát triển. Bằng cách tổng hợp kết quả từ các bài báo thực nghiệm của các tác giả trên thế giới, tác giả so sánh kết quả mô phỏng với các kết quả đó để đánh giá tính ưu việt của phương pháp Front Tracking. Cụ thể, phương pháp tập trung giải quyết các bài toán về truyền nhiệt, trường dòng, sự thay đổi thể tích và đường chập ba pha.

3.1. Giới Thiệu Phương Pháp Theo Dấu Biên Front Tracking

Phương pháp Front Tracking là một phương pháp mô phỏng quá trình hóa rắn hiệu quả, đặc biệt trong việc theo dõi và mô tả sự di chuyển của các bề mặt phân pha. Phương pháp này sử dụng một lưới rời rạc để biểu diễn miền tính toán và một tập hợp các điểm Lagrangian để theo dõi vị trí và hình dạng của bề mặt phân pha. Các phương trình bảo toàn (khối lượng, động lượng, năng lượng) được giải trên lưới Euler, trong khi vị trí của các điểm Lagrangian được cập nhật theo thời gian bằng cách sử dụng vận tốc nội suy từ trường vận tốc trên lưới Euler.

3.2. Ưu Điểm Nổi Bật Của Phương Pháp Front Tracking

Phương pháp Front Tracking có nhiều ưu điểm so với các phương pháp mô phỏng khác, bao gồm khả năng mô tả chính xác hình dạng và sự di chuyển của bề mặt phân pha, khả năng xử lý các bài toán có độ biến dạng lớn của bề mặt phân pha, và khả năng mô phỏng các hiện tượng vật lý phức tạp như sức căng bề mặt, mao dẫn và sự thay đổi pha. Phương pháp theo dấu biên còn được chứng minh là hiệu quả trong việc mô phỏng các vấn đề đa vật lý, nơi có sự tương tác giữa nhiều hiện tượng vật lý, ví dụ như truyền nhiệt, truyền khối lượng và động lực học chất lỏng.

IV. Phân Tích Ảnh Hưởng Các Yếu Tố Đến Quá Trình Hóa Rắn Số

Nghiên cứu này tập trung vào phân tích ảnh hưởng của các yếu tố như góc phát triển, góc tiếp xúc và sự thay đổi thể tích đến quá trình hóa rắn. Sử dụng phương pháp tính toán số, luận văn chỉ ra vai trò của từng yếu tố và ảnh hưởng của chúng đến hình dạng, cấu trúc của sản phẩm sau hóa rắn. Bên cạnh đó, luận văn còn so sánh kết quả mô phỏng với các thí nghiệm thực tế để kiểm chứng độ chính xác của phương pháp.

4.1. Vai Trò Của Góc Phát Triển Trong Quá Trình Hóa Rắn

Góc phát triển là một yếu tố quan trọng ảnh hưởng đến hình dạng và cấu trúc của sản phẩm sau quá trình hóa rắn. Góc phát triển được định nghĩa là góc giữa bề mặt phân pha rắn-lỏng và chất nền. Góc phát triển phụ thuộc vào tính chất vật lý của vật liệu và điều kiện biên. Các nghiên cứu cho thấy rằng góc phát triển có thể ảnh hưởng đến tốc độ hóa rắn, hình dạng tinh thể và sự phân bố tạp chất.

4.2. Ảnh Hưởng Của Góc Tiếp Xúc Và Sự Thay Đổi Thể Tích

Góc tiếp xúc và sự thay đổi thể tích cũng là những yếu tố quan trọng cần xem xét trong mô phỏng quá trình hóa rắn. Góc tiếp xúc là góc giữa bề mặt phân pha lỏng-khí và chất nền. Sự thay đổi thể tích xảy ra do sự khác biệt về khối lượng riêng giữa pha lỏng và pha rắn. Các yếu tố này có thể ảnh hưởng đến hình dạng của giọt chất lỏng, sự phân bố nhiệt độ và tốc độ hóa rắn. Phân tích ứng suất trong quá trình hóa rắn cũng quan trọng để dự đoán các khuyết tật tiềm ẩn hoặc biến dạng trong sản phẩm cuối cùng.

V. Ứng Dụng Mô Phỏng Số Trong Nghiên Cứu Hóa Rắn Vật Liệu

Luận văn trình bày kết quả mô phỏng quá trình hóa rắn của hạt chất lỏng khác nhau (nước, silic, gecmani) trên bề mặt lạnh. Kết quả mô phỏng được so sánh với các thí nghiệm thực tế để đánh giá độ chính xác của phương pháp. Nghiên cứu này đóng góp vào việc hiểu rõ hơn về quá trình hóa rắn và cung cấp công cụ hiệu quả để thiết kế, tối ưu hóa quá trình sản xuất vật liệu.

5.1. Mô Phỏng Hóa Rắn Nước So Sánh Với Thực Nghiệm

Kết quả mô phỏng quá trình hóa rắn nước cho thấy sự phù hợp với các thí nghiệm thực tế về hình dạng, tốc độ hóa rắn và sự phân bố nhiệt độ. Các kết quả cũng cho thấy ảnh hưởng của góc phát triển đến hình dạng của sản phẩm sau hóa rắn. Qua phân tích và so sánh, ta có thể đánh giá và điều chỉnh các tham số mô phỏng để đạt độ chính xác cao.

5.2. Mô Phỏng Silic Gecmani Vật Liệu Bán Dẫn Tiềm Năng

Mô phỏng quá trình hóa rắn silic và gecmani cho thấy sự khác biệt về hình dạng và cấu trúc do sự khác biệt về tính chất vật lý. Các kết quả này có ý nghĩa quan trọng trong việc thiết kế, tối ưu hóa quá trình sản xuất vật liệu bán dẫn. Sự phát triển của các kỹ thuật mô phỏng mới giúp hiểu rõ hơn về mối quan hệ giữa các thông số xử lý và đặc tính vật liệu.

VI. Kết Luận Hướng Phát Triển Nghiên Cứu Quá Trình Hóa Rắn

Nghiên cứu đã thành công trong việc mô phỏng quá trình hóa rắn hạt chất lỏng trên bề mặt lạnh bằng phương pháp Front Tracking. Kết quả mô phỏng cho thấy sự phù hợp với các thí nghiệm thực tế. Nghiên cứu này mở ra hướng phát triển mới trong việc thiết kế, tối ưu hóa quá trình sản xuất vật liệu và ứng dụng trong nhiều lĩnh vực khác nhau. Tuy nhiên, cần tiếp tục nghiên cứu để hoàn thiện phương pháp, xem xét thêm các yếu tố ảnh hưởng và áp dụng vào các bài toán thực tế phức tạp hơn.

6.1. Tóm Tắt Kết Quả Nghiên Cứu Đóng Góp Mới

Nghiên cứu đã xây dựng và kiểm chứng phương pháp mô phỏng quá trình hóa rắn hạt chất lỏng bằng Front Tracking. Kết quả cho thấy ảnh hưởng của góc phát triển, góc tiếp xúc và sự thay đổi thể tích đến hình dạng và cấu trúc của sản phẩm sau hóa rắn. Nghiên cứu này cung cấp công cụ hiệu quả để nghiên cứu, thiết kế và tối ưu hóa quá trình sản xuất vật liệu.

6.2. Hướng Nghiên Cứu Tiếp Theo Về Quá Trình Hóa Rắn

Hướng nghiên cứu tiếp theo cần tập trung vào việc hoàn thiện phương pháp Front Tracking, xem xét thêm các yếu tố ảnh hưởng (đối lưu nhiệt, bức xạ nhiệt, tương tác pha) và áp dụng vào các bài toán thực tế phức tạp hơn (hóa rắn trong điều kiện chân không, hóa rắn với sự có mặt của tạp chất). Nghiên cứu cần chú trọng đến việc tích hợp các mô hình hóa học và vật lý để mô tả đầy đủ quá trình hóa rắn không cân bằng.

Nội dung chính

Tổng quan nghiên cứu

Hiện tượng hóa rắn của hạt chất lỏng trên bề mặt lạnh là một hiện tượng phổ biến trong tự nhiên và công nghiệp, gây ảnh hưởng lớn đến hiệu suất và tuổi thọ thiết bị. Ví dụ điển hình là băng đá hình thành trên cánh tuabin gió hoặc hiện tượng đóng băng trên cánh máy bay, làm giảm hiệu quả hoạt động và tiềm ẩn nguy cơ tai nạn nghiêm trọng. Theo báo cáo ngành, hiện tượng này cũng được ứng dụng trong sản xuất pin năng lượng mặt trời và làm sạch bụi công nghiệp, góp phần bảo vệ môi trường và tiết kiệm nguyên liệu. Mục tiêu nghiên cứu của luận văn là mô phỏng số quá trình hóa rắn của hạt chất lỏng trên bề mặt lạnh, tập trung vào các vật liệu silic, gec-ma-ni và nước, nhằm hiểu rõ cơ chế và ảnh hưởng của các thông số như góc tiếp xúc, góc phát triển và đặc tính vật liệu đến quá trình hóa rắn. Nghiên cứu được thực hiện tại Hà Nội trong giai đoạn 2017-2019, sử dụng phương pháp mô phỏng số theo phương pháp theo dõi biên (front-tracking). Kết quả nghiên cứu có ý nghĩa quan trọng trong việc dự đoán hình dạng hạt sau khi hóa rắn, từ đó đề xuất giải pháp ngăn ngừa hiện tượng đóng băng không mong muốn và ứng dụng trong sản xuất công nghiệp.

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

Khung lý thuyết áp dụng

Luận văn dựa trên hai khung lý thuyết chính: phương trình Navier-Stokes mô tả chuyển động chất lỏng và phương trình năng lượng liên quan đến truyền nhiệt và chuyển pha. Phương pháp mô phỏng số theo dõi biên (front-tracking) được áp dụng để mô hình hóa sự chuyển động của mặt phân cách giữa các pha rắn, lỏng và khí. Các khái niệm chính bao gồm:

• Góc tiếp xúc (ϕ₀): Góc giữa hạt chất lỏng và bề mặt lạnh, ảnh hưởng đến thời gian và hình dạng hóa rắn.
• Góc phát triển (θ): Góc sinh ra khi hạt chất lỏng hóa rắn, đặc trưng cho sự phát triển của pha rắn.
• Sức căng bề mặt (σ): Lực tác động tại mặt phân cách giữa các pha, ảnh hưởng đến hình dạng hạt sau hóa rắn.
• Thông số không thứ nguyên: Bao gồm số Prandtl (Pr), số Stefan (St), số Bond (Bo), số Weber (We), và số Ohnesorge (Oh), dùng để mô tả các đặc tính truyền nhiệt, chuyển pha và động lực học của hệ.

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

Nguồn dữ liệu chính là các kết quả thực nghiệm từ các nghiên cứu trước về hóa rắn của hạt nước, silic và gec-ma-ni trên bề mặt lạnh. Phương pháp phân tích sử dụng mô phỏng số dựa trên phương trình Navier-Stokes kết hợp với phương trình năng lượng, áp dụng phương pháp theo dõi biên để mô hình hóa mặt phân cách giữa các pha. Cỡ mẫu mô phỏng là lưới 192×192 điểm, được lựa chọn dựa trên phân tích hội tụ và độ chính xác của kết quả. Thời gian nghiên cứu kéo dài từ 2017 đến 2019, với các bước tính toán chi tiết bao gồm giải phương trình động lực học, truyền nhiệt và cập nhật vị trí mặt phân cách theo thời gian. Phương pháp lưới so le được sử dụng để tăng độ ổn định và chính xác của mô phỏng. Các điều kiện biên được thiết lập phù hợp với thực tế, bao gồm điều kiện không trượt và nhiệt độ cố định trên bề mặt lạnh.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

1. Ảnh hưởng của góc tiếp xúc đến thời gian hóa rắn:
   Mô phỏng cho thấy với hạt nước có thể tích khoảng 56,34 μL, khi góc tiếp xúc tăng từ 76° lên 155°, thời gian hóa rắn kéo dài đáng kể. Cụ thể, thời gian hóa rắn với góc 76° nhanh hơn khoảng 30% so với góc 155°. Điều này được giải thích do diện tích tiếp xúc giữa hạt và bề mặt lạnh giảm khi góc tiếp xúc tăng, làm giảm tốc độ truyền nhiệt.

2. Hình dạng hạt sau hóa rắn phụ thuộc vào góc phát triển:
   Với vật liệu silic và gec-ma-ni, góc phát triển lần lượt là 12° và 14°, mô phỏng cho thấy hạt sau hóa rắn có hình dạng chóp nón rõ rệt. Hình dạng này phù hợp với kết quả thực nghiệm, xác nhận tính chính xác của phương pháp mô phỏng. Đặc biệt, hạt silic có hình chóp nón nhọn hơn so với gec-ma-ni do góc phát triển nhỏ hơn.

3. Sự khác biệt về đặc tính vật liệu ảnh hưởng đến quá trình hóa rắn:
   Các thông số vật liệu như khối lượng riêng, nhiệt dung riêng, hệ số dẫn nhiệt và nhiệt ẩn chuyển pha được tính toán và đưa vào mô hình. Kết quả cho thấy hạt silic có tốc độ hóa rắn nhanh hơn so với gec-ma-ni do hệ số dẫn nhiệt cao hơn (31 W/mK so với 60 W/mK cho pha rắn). Tỷ lệ khối lượng riêng rắn/lỏng cũng ảnh hưởng đến sự giãn nở thể tích khi hóa rắn, từ đó ảnh hưởng đến hình dạng cuối cùng của hạt.

4. Mô hình mô phỏng số khớp với kết quả thực nghiệm:
   So sánh kết quả mô phỏng với các thí nghiệm của Huang và cộng sự (2012) về hạt nước và Satunkin về silic, gec-ma-ni cho thấy sự tương đồng cao về hình dạng và thời gian hóa rắn. Sai số giữa mô phỏng và thực nghiệm nằm trong khoảng 5-10%, chứng tỏ phương pháp mô phỏng số theo dõi biên là công cụ hiệu quả để nghiên cứu hiện tượng này.

Thảo luận kết quả

Nguyên nhân chính của các phát hiện trên là do sự thay đổi diện tích tiếp xúc và đặc tính truyền nhiệt giữa các vật liệu và góc tiếp xúc khác nhau. So với các nghiên cứu trước đây chỉ tập trung vào thực nghiệm, luận văn đã phát triển một mô hình mô phỏng số chi tiết, cho phép dự đoán chính xác hơn về quá trình hóa rắn và hình dạng hạt sau khi đóng băng. Kết quả mô phỏng có thể được trình bày qua biểu đồ thể hiện mối quan hệ giữa góc tiếp xúc và thời gian hóa rắn, cũng như bảng so sánh các thông số vật liệu và kết quả thực nghiệm. Ý nghĩa của nghiên cứu nằm ở việc cung cấp công cụ dự đoán và kiểm soát hiện tượng hóa rắn trong các ứng dụng công nghiệp, từ đó giảm thiểu thiệt hại và nâng cao hiệu quả sản xuất.

Đề xuất và khuyến nghị

1. Phát triển hệ thống kiểm soát góc tiếp xúc trên bề mặt lạnh:
   Thiết kế bề mặt có khả năng điều chỉnh góc tiếp xúc để giảm thời gian hóa rắn, hướng tới giảm thiểu hiện tượng đóng băng không mong muốn. Chủ thể thực hiện: các nhà sản xuất thiết bị công nghiệp; Thời gian: 1-2 năm.

2. Ứng dụng mô hình mô phỏng số trong thiết kế sản phẩm:
   Áp dụng mô hình mô phỏng để dự đoán và tối ưu hóa hình dạng hạt đóng băng trong sản xuất pin năng lượng mặt trời và các ngành công nghiệp liên quan. Chủ thể thực hiện: các viện nghiên cứu và doanh nghiệp công nghệ; Thời gian: 2 năm.

3. Nghiên cứu mở rộng các vật liệu khác và điều kiện môi trường:
   Mở rộng nghiên cứu sang các vật liệu khác và điều kiện nhiệt độ, áp suất đa dạng để tăng tính ứng dụng của mô hình. Chủ thể thực hiện: các nhóm nghiên cứu khoa học; Thời gian: 3 năm.

4. Phát triển phần mềm mô phỏng tích hợp:
   Xây dựng phần mềm mô phỏng thân thiện, tích hợp các tham số vật liệu và điều kiện thực tế, hỗ trợ các kỹ sư trong việc thiết kế và kiểm soát quá trình hóa rắn. Chủ thể thực hiện: các công ty phần mềm kỹ thuật; Thời gian: 1-2 năm.

Đối tượng nên tham khảo luận văn

1. Nhà nghiên cứu và giảng viên kỹ thuật cơ khí, nhiệt động lực học:
   Học hỏi phương pháp mô phỏng số và ứng dụng trong nghiên cứu chuyển pha và truyền nhiệt.

2. Kỹ sư thiết kế và vận hành thiết bị công nghiệp:
   Áp dụng kết quả nghiên cứu để cải tiến thiết kế bề mặt và quy trình vận hành nhằm giảm thiểu hiện tượng đóng băng.

3. Doanh nghiệp sản xuất pin năng lượng mặt trời và vật liệu bán dẫn:
   Tận dụng mô hình để tối ưu hóa quá trình sản xuất, nâng cao chất lượng sản phẩm.

4. Sinh viên cao học và nghiên cứu sinh chuyên ngành kỹ thuật máy thủy khí, kỹ thuật cơ khí động lực:
   Tham khảo phương pháp nghiên cứu, mô hình hóa và phân tích kết quả trong luận văn để phát triển đề tài nghiên cứu sâu hơn.

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

1. Phương pháp mô phỏng số theo dõi biên (front-tracking) là gì?
   Đây là phương pháp mô phỏng dòng đa pha, trong đó mặt phân cách giữa các pha được theo dõi trực tiếp bằng các điểm hoặc phần tử riêng biệt, giúp mô hình hóa chính xác sự chuyển động và biến dạng của mặt phân cách trong quá trình hóa rắn.

2. Góc tiếp xúc ảnh hưởng như thế nào đến quá trình hóa rắn?
   Góc tiếp xúc càng nhỏ thì diện tích tiếp xúc giữa hạt và bề mặt lạnh càng lớn, dẫn đến truyền nhiệt nhanh hơn và thời gian hóa rắn ngắn hơn. Ngược lại, góc lớn làm giảm diện tích tiếp xúc, kéo dài thời gian hóa rắn.

3. Tại sao chọn vật liệu silic và gec-ma-ni để nghiên cứu?
   Silic và gec-ma-ni là vật liệu bán dẫn phổ biến, có đặc tính vật lý khác biệt rõ rệt, giúp nghiên cứu sự ảnh hưởng của đặc tính vật liệu đến quá trình hóa rắn, đồng thời có nhiều ứng dụng công nghiệp quan trọng.

4. Lưới so le có ưu điểm gì trong mô phỏng?
   Lưới so le giúp tăng độ ổn định và chính xác của mô phỏng bằng cách lưu trữ áp suất và vận tốc tại các vị trí lệch nhau, giảm thiểu sai số số học và cải thiện khả năng hội tụ của phương pháp số.

5. Kết quả mô phỏng có thể ứng dụng thực tế như thế nào?
   Kết quả giúp dự đoán hình dạng và thời gian hóa rắn của hạt chất lỏng trên bề mặt lạnh, từ đó thiết kế bề mặt và quy trình vận hành để ngăn ngừa đóng băng, nâng cao hiệu quả và an toàn trong các ngành công nghiệp như hàng không, năng lượng tái tạo và sản xuất vật liệu.

Kết luận

• Luận văn đã phát triển thành công mô hình mô phỏng số theo dõi biên để nghiên cứu quá trình hóa rắn của hạt chất lỏng trên bề mặt lạnh với các vật liệu silic, gec-ma-ni và nước.
• Kết quả mô phỏng cho thấy góc tiếp xúc và góc phát triển là các yếu tố quan trọng ảnh hưởng đến thời gian và hình dạng hóa rắn.
• Mô hình mô phỏng số khớp với kết quả thực nghiệm, chứng minh tính chính xác và khả năng ứng dụng trong thực tế.
• Nghiên cứu góp phần nâng cao hiểu biết về hiện tượng hóa rắn, hỗ trợ phát triển các giải pháp kỹ thuật trong công nghiệp.
• Các bước tiếp theo bao gồm mở rộng nghiên cứu với nhiều vật liệu và điều kiện môi trường khác, đồng thời phát triển phần mềm mô phỏng tích hợp để ứng dụng rộng rãi hơn.

Hành động khuyến nghị: Các nhà nghiên cứu và kỹ sư nên áp dụng mô hình mô phỏng này để tối ưu hóa thiết kế và vận hành thiết bị, đồng thời tiếp tục phát triển nghiên cứu nhằm mở rộng ứng dụng trong các lĩnh vực công nghiệp khác.

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

B GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NI ---------------------------------------- PHẠM DUY BÍNH NGHIÊN CỨU TÍNH TOÁN SỐ QUÁ TRÌNH HÓA RN CỦA HẠT CHẤT LỎNG LUN VĂN THẠC SĨ KHOA HỌC Chuyên ngành: Kỹ thut Cơ khí Đng lực NGƯỜI HƯNG DẪN KHOA HỌC: TS. VŨ VĂN TRƯỜNG Hà Nội - 2019 17051113913231000000 Lun văn thc sỹ Chuyên ngành: Kỹ thut máy thy khí LI CAM ĐOAN Tôi ậ Phm Duy Bính - Xin cam đoan đơy lƠ nghiên cu ca riêng tôi di sự hng dn ca TS. Các s liu và kt qu trong lun văn lƠ trung thực, không phi là sao chép toƠn văn ca bt kỳ công trình nào khác. Hà Ni, 23 tháng 11 năm 2018 Tác gi Phm Duy Bính XÁC NHN CAăNGIăHNG DN KHOA HC Tôi, Vũ Văn Trng, xác nhn lun văn thc sỹ khoa hc ca tác gi Phm Duy Bính là công trình nghiên cu ca tác gi di sự hng dn ca tôi.

Các kt qu nêu trong lun văn lƠ trung thực vƠ đm bo đy đ tính khoa hc ca mt lun văn thc sỹ. HƠ ni ngƠy 03 tháng 05 năm 2019 Giáoăviênăhngădn TS.ăVǛăVĔNăTRNG Thc sỹ khoa hc Trang 1 Lun văn thc sỹ Chuyên ngành: Kỹ thut máy thy khí LI NÓI ĐU Hin tng hóa rắn ca ht cht lng trên mt b mặt lnh xut hin khá nhiu trong tự nhiên cũng nh trong công nghip. Hin tng này có th gây hi hoặc có th ng dng trong sn xut công nghip. Vic hiu bit v quá trình hóa rắn ca ht cht lng là vô cùng quan trng góp phn trong vic phòng tránh tác hi ca hin tng cũng nh ng dng ca hin tng này trong sn xut vƠ đi sng.

Vì tính cp thit cn mt đ tài nghiên cu v hin tng này lun văn “σghiên cứu tính toán số quá trình hóa r ắn của hạt chất lỏng” ra đi. Vi t cách lƠ tác gi ca lun văn nƠy, tôi xin gửi li cm n chân thành và sâu sắc đn TS. Vũ Văn Trng, ngi đư trực tip hng dn tôi đ tôi có th hoàn thành lun văn nƠy. Đng thi tôi cũng xin gửi li cm n chơn thƠnh đn tp th các thy, cô giáo ca Vin C khí đng lực, phòng Đào to Trng Đi hc Bách Khoa Hà Ni đư tn tình chỉ bo, truyn đt những kin thc mi cho tôi, giúp tôi sut thi gian hc tp và làm lun văn nƠy.

Cui cùng xin gửi li cm n đn gia đình, bn bè tôi đư giúp đỡ và ng h tôi trong sut thi gian tôi hc tp và làm lun văn nƠy. Hà Ni, ngƠy 23 tháng 11 năm 2018 Hc viên Thc sỹ khoa hc Trang 2 Lun văn thc sỹ Chuyên ngành: Kỹ thut máy thy khí MC LC M ĐU. Lý do chn đ tài. Mc đích nghiên cu ca lun văn.

Tóm tắt và những đóng góp mi trong lun văn. Phng pháp nghiên cu. 10 CHNG 1: TNG QUAN VĨ ụ NGHƾA CA Đ TÀI .1 Hin tng hóa rắn ht cht lng trong tự nhiên và trong công nghip 11 1.2 Tình hình nghiên cu bài toán trên th gii cũng nh trong nc .3 Các ht cht lng nghiên cu trong lun văn vƠ tm quan trng ca nghiên cu. 15 CHNG 2: CÁC PHNG TRÌNH C BN VĨ PHNG PHÁP MÔ PHNG .1 Các phng trình tng quát .2 Tích phân theo thi gian .4 Ri rc hóa các phng trình tng quát .2 Ri rc các thành phn khuch tán .3 Phng trình áp sut .5 Ri rc hóa phng trình năng lng .8 Phng pháp theo du biên .1 Cu trúc mặt phân cách .2 Chuyn đng ca mặt phân cách.3 Cu trúc li mặt phân cách .4 Làm mịn và tái cu trúc hàm chỉ thị.

34 Thc sỹ khoa hc Trang 3 Lun văn thc sỹ Chuyên ngành: Kỹ thut máy thy khí 2.5 Tính toán sc căng b mặt trong li .6 Tính toán cân bằng nhit ti b mặt chuyn pha. 35 CHNG 3: KT QU CA QUÁ TRÌNH MÔ PHNG VÀ THO LUN .1 Quá trình mô phng hóa rắn ca mt s ht vt liu chuyn pha .3 Gec-ma-ni .2 Đánh giá phng pháp mô phng .2 Silic và gec-ma-ni .3 Phng pháp theo du biên vi ht kép .3 Xây dựng biu đ pha vi các góc t khác nhau ca vt liu silic. 55 CHNG 4: KT LUN VĨ HNG PHÁT TRIN. 57 TÀI LIU THAM KHO .58 CÁC CÔNG TRÌNH Đĩ CÔNG B.

61 Thc sỹ khoa hc Trang 4 Lun văn thc sỹ Chuyên ngành: Kỹ thut máy thy khí DANH MC HÌNH VẼ Hình 1: Băng đá hình thƠnh trên cánh tuabin gió [1]. 8 Hình 2: Ht bán dn trong ch to pin năng lng mặt tri [2]. 1: Ht nc đóng băng trên mt b mặt lnh [5]. 2: Bên trái lƠ đng biên ht nc trc và sau hóa lng [6], bên phi là so sánh mô phng s ca silic vi kt qu mô phng ca satunkin [7].

3: Thí nghim hóa rắn ht nc vi các góc t khác nhau ca Huang và cng sự [5]. 1: B trí trong li so le [12]. 2: Li vn tc theo phng ngang [12]. 3: Li vn tc theo phng thẳng đng [12].

4: Li cu trúc thông thng. 5: Ký hiu đc sử dng cho li chuẩn so le [12]. 6: Ký hiu đc sử dng cho li chuẩn so le [12]. 7: KỦ pháp đc sử dng cho li MAC so le chuẩn [12].

8: Điu kin ranh gii trong li so le [13]. 9: Cu trúc mặt phân cách giữa các pha [13]. 11: Loi b mt thành phn [13]. 12: Di chuyn ca biên chuyn pha theo nhit đ.

1: Sự phát trin ca trng nhit đ theo thi gian không th nguyên vi ϕo=76o ca nc. 2: Sự phát trin ca trng nhit đ theo thi gian không th nguyên vi góc t ϕo=124o ca nc. 3: Sự phát trin ca trng nhit đ theo thi gian không th nguyên vi góc t ϕo=155o ca nc. 4: Sự phát trin ca khói lng riêng theo thi gian không th nguyên vi ϕo=33o ca silic.

5: Sự phát trin ca trng nhit đ theo thi gian không th nguyên vi góc t ϕo=33o ca gec-ma-ni. 6: So sánh giữa mô phng (bên phi) và thử nghim (trái [5]) ca git nc vi ϕo=76o. 7: So sánh giữa mô phng (bên phi) và thử nghim (trái [5]) ca git nc khi đóng băng vi ϕo=124o. 8: So sánh giữa mô phng (bên phi) và thử nghim (trái [5]) ca git nc khi đóng băng vi ϕo=155o.

49 Thc sỹ khoa hc Trang 5 Lun văn thc sỹ Chuyên ngành: Kỹ thut máy thy khí Hình 3. 9: So sánh giữa mô phng (bên phi) và thí nghim (trái [7]) vi ϕo=33o ca git silic kt tinh. 10: So sánh giữa mô phng (phi) và thí nghim (trái [7]) vi ϕo =33o ca git gec-ma-ni kt tinh. 11: Ht kép bin dng trong dòng chy trt.

12: Nghiên cu biên ht kép ti  = 9.1 bằng các li khác nhau. 13: So sánh sự bin dng ca ht kép trong dòng chy trt giữa các tính toán hin ti vi Hua và cng sự [21]. 14: Quá trình hóa rắn ht silic vi các góc t 33o, 45o, 60o trong mô phng s. 15: Quá trình hóa rắn ht silic vi các góc t 75o, 90o, 105o, 120o trong mô phng s.

16: Mi liên h giữa chiu cao hóa rắn vi góc t và sự thay đi chiu cao trc và sau khi hóa rắn vi góc t ca silic. 55 Thc sỹ khoa hc Trang 6 Lun văn thc sỹ Chuyên ngành: Kỹ thut máy thy khí DANH MC CÁC BNG Bng 3. 1: Các thông s đặc tính ca nc. 2: Mt s thông s không th nguyên quan trng ca nc.

3: Các thông s đặc tính ca silic. 4: Mt s thông s không th nguyên quan trng ca silic. 5: Các thông s đặc tính ca gec-ma-ni. 6: Mt s thông s không th nguyên quan trng ca gec-ma-ni.

46 Thc sỹ khoa hc Trang 7 Lun văn thc sỹ Chuyên ngành: Kỹ thut máy thy khí M ĐU 1. Lý do chnăđ tài Băng đá hình thƠnh trên các b mặt cánh máy cánh dn (nh hình 1 nc đóng băng trên cánh ca tuabin gió) gây nh hng ln đn kh năng hot đng ca máy [1]. Đơy lƠ nguyên nhơn dn đn tình trng máy hot đng rung lắc, công sut không đc nh tính toán ban đu, tui th ca máy suy gim. Ngoài ra hin tng này còn gây ra nhiu v tai nn thm khc trong ngành hàng không.

Do sự hóa rắn ca ht nc trên b mặt cánh máy bay mà dòng khí khi chy qua profile cánh không còn đc nh tính toán, lực nâng gim và lực cn tăng. Từ đó thy đc tm quan trng ca vic nghiên cu ht cht lng hóa rắn trên mt b mặt lnh đ tìm ra gii pháp đ ngăn chặn các hin tng này. Ngoài các hin tng tiêu cực ca hin tng hóa rắn ca ht cht lng  trên thì ta thy hin tng nƠy còn đc ng dng trong công nghip. Hình 2 th hin những ht bán dn trong ch to pin năng lng mặt tri [2].

Ngành công nghip năng lng tái to đang phát trin mnh m trong các năm gn đơy. Nguyên nhân là do nh hng tiêu cực ca vic sử dng quá mc ngun nhiên liu hóa thch. Các nghiên cu v sự hóa rắn ca ht cht lng trên mt b mặt lnh đc quan tâm trong ngành công nghip này. Hin tng này có th ng dng trong vic ch to các tm pin mặt tri thông qua hóa rắn các ht bán dn đ sn xut những tm pin mặt tri đi trà rẻ hn vƠ rng rưi hn.

Thêm mt ví d nữa là làm sch bi công nghip nh hin tng này. Các ht nc đc phun ra từ vòi phun qun ly ht bi công nghip và bám vào mt b mặt lnh đư đc chuẩn bị sẵn. Hóa rắn chúng và gom chúng li có th tái sử dng hoặc xử lỦ đ không gây ô nhim môi trng,ầ Hnh 1:ăBĕngăđáăhnhăthƠnhătrênăcánh tuabin gió [1]. Thc sỹ khoa hc Trang 8 Lun văn thc sỹ Chuyên ngành: Kỹ thut máy thy khí Hnh 2: Ht bán dn trong ch toăpinănĕngălng mặt tri [2].

Qua các tác hi ca hin tng trong tự nhiên cũng nh ng dng ca nó trong công nghip thì vic hiu rõ nó là rt cn thit. Vì vy lun văn “Nghiên cu mô phng sự hóa rắn ca cht lng trên mt b mặt lnh” ra đi nhằm đóng góp mt phn nƠo đó cho sự hiu bit ca bn đc và nêu ra gii pháp và ng dng hin tng đó. Đơy cũng lƠ c s đ phát trin các nghiên cu chuyên sâu sau này.