Tổng quan nghiên cứu

Trong kỹ thuật hàng không hiện đại, nhiệt độ dòng khí cháy đi vào tuabin khí thường vượt quá 1400 K đến 1750 K, cao hơn đáng kể so với giới hạn chịu nhiệt cơ học của các hợp kim siêu bền (chỉ chịu được khoảng 1100 K đến 1150 K). Để ngăn chặn nguy cơ hư hỏng cấu trúc và kéo dài tuổi thọ cánh tuabin, công nghệ làm mát bằng màng khí qua các tấm nhiều lỗ là giải pháp kỹ thuật sống còn. Vấn đề cốt lõi đặt ra là làm thế nào để tối ưu hóa đồng thời các thông số hình học và động học nhằm đạt hiệu quả bảo vệ nhiệt cao nhất nhưng chỉ tiêu tốn lượng khí trích làm mát tối thiểu từ máy nén.

Nghiên cứu tập trung giải quyết bài toán tương tác phức hợp giữa dòng khí nóng chính và các tia khí làm mát phản lực phun qua bề mặt thành buồng cháy và lá tuabin máy bay. Mục tiêu cụ thể là khảo sát chi tiết sự biến thiên của trường vận tốc, trường nhiệt độ và hiệu suất làm mát khi thay đổi góc nghiêng lỗ, tỉ số phun khối lượng và cách bố trí các hàng lỗ, đồng thời đánh giá ứng xử nhiệt trên cấu trúc đa lớp gồm lớp phủ gốm Ceramic cách nhiệt và vật liệu thép nền.

Phạm vi nghiên cứu được thực hiện thông qua mô hình hóa số động lực học dòng chảy tại Trường Đại học Bách khoa Hà Nội, thiết lập trên các dải thông số vận hành mô phỏng điều kiện thực tế của động cơ phản lực. Nghiên cứu mang lại ý nghĩa kinh tế và kỹ thuật vượt bậc: việc tối ưu hóa cấu hình màng khí giúp giảm nhiệt độ bề mặt kim loại từ 120 K đến 150 K, nâng cao tuổi thọ làm việc của cánh tuabin thêm khoảng 25% đến 30% và tiết kiệm tới 12% lưu lượng khí nén phục vụ làm mát.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng kết hợp của ba trụ cột lý thuyết nhiệt - thủy khí động lực học chuyên sâu:

Thứ nhất là Thuyết màng khí bảo vệ và truyền nhiệt đối lưu trong lớp biên. Khi dòng khí làm mát có nhiệt độ thấp được phun qua các lỗ nghiêng trên bề mặt, một lớp màng chất lưu mỏng được hình thành để ngăn cách dòng khí cháy nhiệt độ cao tiếp xúc trực tiếp với bề mặt kim loại. Hiệu quả làm mát đoạn nhiệt được xác định thông qua chênh lệch nhiệt độ không thứ nguyên giữa dòng khí chính, dòng phun và bề mặt thành.

Thứ hai là Lý thuyết mô hình hóa dòng chảy rối Navier-Stokes trung bình thời gian. Trong môi trường tuabin với xoáy cuộn mạnh, các mô hình rối hai phương trình thuộc họ k-epsilon đóng vai trò quyết định trong việc đóng kín hệ phương trình vi phân mô tả chuyển động và năng lượng của dòng khí.

Thứ ba là Lý thuyết truyền nhiệt liên hợp đa lớp và ứng xử nhiệt cấu trúc. Quá trình truyền nhiệt không chỉ dừng lại ở đối lưu bề mặt mà còn kết hợp dẫn nhiệt qua lớp rào chắn nhiệt Ceramic có độ dẫn nhiệt thấp và truyền nhiệt vào bản thép nền chịu lực.

Năm khái niệm then chốt được lượng hóa bao gồm: Tỉ số phun khối lượng M biểu thị tương quan mật độ và vận tốc giữa dòng phun và dòng chính; Hiệu suất làm mát màng ký hiệu eta; Khoảng cách không thứ nguyên x trên d dọc theo chiều dòng chảy; Bước lỗ tương đối P trên d; và Chiều dày tương đối của lớp gốm phủ cách nhiệt.

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

Nghiên cứu ứng dụng phương pháp mô phỏng số động lực học dòng chảy (CFD) thông qua phần mềm ANSYS Fluent, kết hợp thuật toán phân tích truyền nhiệt liên hợp.

Cỡ mẫu mô phỏng được thiết kế toàn diện với 18 kịch bản kết hợp đa biến hình học và nhiệt động lực học. Miền tính toán 3 chiều được rời rạc hóa với mật độ lưới từ 1,2 triệu đến 2,8 triệu phần tử lục giác có cấu trúc chất lượng cao. Phương pháp chọn mẫu lưới áp dụng kỹ thuật làm mịn cục bộ thích nghi cao độ tại vùng miệng lỗ phun và sát vách, khống chế giá trị y-cộng xấp xỉ 1 nhằm nắm bắt trọn vẹn hiện tượng truyền nhiệt trong lớp biên nhớt.

Lý do lựa chọn mô hình dòng rối Realizable k-epsilon kết hợp hàm truyền nhiệt sát vách nâng cao là vì mô hình này vượt trội trong việc dự đoán chính xác sự phát triển của cặp xoáy ngược chiều hình quả thận và hiện tượng khuếch tán nhiệt trong dòng phản lực nghiêng. Phương pháp số giúp tiết kiệm hơn 85% chi phí nghiên cứu so với thử nghiệm buồng đốt thực tế, đồng thời kiểm soát sai số dưới ngưỡng 4,8% so với các dữ liệu thực nghiệm kinh điển. Toàn bộ chu trình tính toán và kiểm chứng độc lập lưới được hoàn tất trong thời gian 12 tháng.

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

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

Quá trình mô phỏng số đã làm sáng tỏ các quy luật biến thiên nhiệt - động lực học với bốn phát hiện then chốt sau:

Thứ nhất, ảnh hưởng của góc nghiêng lỗ phun thể hiện rõ nét qua hiệu suất làm mát cực đại. Lỗ nghiêng góc 35 độ duy trì tia làm mát bám dính bề mặt tốt nhất, đạt hiệu suất cực đại 0,68 tại vị trí x bằng 1d ngay sau miệng lỗ. Khi góc nghiêng tăng lên 45 độ, hiệu suất giảm 19,5%; và khi góc nghiêng đạt 60 độ, hiệu suất sụt giảm mạnh tới 36,2% do tia phun bị tách dòng bắn sâu vào dòng khí nóng chính.

Thứ hai, tác động của tỉ số phun M quyết định hiện tượng bám dính hoặc tách dòng. Tại tỉ số phun tối ưu M bằng 0,5 đến 1,0, màng khí trải rộng ổn định dọc theo bề mặt. Tuy nhiên, khi tỉ số phun tăng lên mức M bằng 1,5, động năng tia phun quá lớn gây ra hiện tượng bật khỏi vách, khiến hiệu suất làm mát cục bộ tại vị trí x bằng 3d sụt giảm 22,4% trước khi màng khí tái bám dính ở vị trí x từ 10d đến 20d.

Thứ ba, cấu hình bố trí số hàng lỗ mang lại hiệu ứng tích lũy vượt trội. Cấu hình 3 hàng lỗ bố trí so le giúp tăng diện tích bao phủ của màng lạnh lên 58%, duy trì hiệu suất làm mát trên 0,45 tại khoảng cách xa x bằng 20d, cao hơn 41,8% so với trường hợp chỉ sử dụng 1 hàng lỗ đơn lẻ.

Thứ tư, vai trò của lớp phủ gốm Ceramic thể hiện sự ưu việt trong việc giảm ứng suất nhiệt. Khi tỉ lệ chiều dày lớp phủ Ceramic trên tổng chiều dày tấm đạt mức 0,25, nhiệt độ truyền tới lớp thép nền giảm mạnh từ 125 K đến 145 K, giúp triệt tiêu hơn 27% ứng suất nhiệt cục bộ phát sinh trong kết cấu kim loại.

+-------------------------------------------------------------------------+
|                  SO SÁNH HIỆU SUẤT LÀM MÁT THEO GÓC NGHIÊNG             |
|                                                                         |
|  Hiệu suất (η)                                                          |
|   1.0 |                                                                 |
|   0.8 |    [35°]                                                        |
|   0.6 |     *---\                                                       |
|   0.4 |          \---\   [45°]                                          |
|   0.2 |               *----\   [60°]                                    |
|   0.0 +---------------------*--------------------> Khoảng cách (x/d)    |
|       x=1d       x=3d       x=5d       x=10d     x=20d                  |
+-------------------------------------------------------------------------+

Thảo luận kết quả

Cơ chế vật lý sâu xa của các phát hiện trên bắt nguồn từ cấu trúc cặp xoáy ngược chiều (Kidney Vortices) sinh ra do sự tương tác giữa tia phun phản lực và dòng khí chính chuyển động ngang. Ở góc nghiêng 35 độ, thành phần vận tốc tiếp tuyến chiếm ưu thế giúp xoáy cuộn yếu hơn, giữ lớp biên màng lạnh áp sát vách. Ngược lại, ở góc 60 độ, thành phần vận tốc pháp tuyến lớn đẩy tia phun vượt qua lớp biên, cuốn dòng khí nóng 1500 K tràn xuống quét trực tiếp vào bề mặt thành.

Dữ liệu nghiên cứu được trình bày trực quan thông qua đồ thị đường 2 chiều biểu diễn sự suy giảm hiệu suất theo khoảng cách không thứ nguyên từ 1d đến 20d, song song với các bảng đối chiếu thông số nhiệt động học. Bên cạnh đó, các biểu đồ đẳng nhiệt 3 chiều đa sắc màu thể hiện rõ rệt thang nhiệt độ giảm dần từ vùng tâm tia phun ra hai biên. Kết quả này hoàn toàn tương thích với các nghiên cứu thực nghiệm trong ngành khí động hàng không, chứng minh tính tin cậy tuyệt đối của phương pháp mô hình hóa số đã chọn.

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

Dựa trên các kết quả tính toán chi tiết, bốn giải pháp kỹ thuật cụ thể được đề xuất nhằm ứng dụng trực tiếp vào quy trình thiết kế và chế tạo động cơ máy bay:

Thứ nhất, chuẩn hóa góc gia công lỗ làm mát nghiêng từ 30 độ đến 35 độ trên bề mặt cánh tuabin và thành buồng đốt. Các nhà máy chế tạo động cơ hàng không cần triển khai ứng dụng công nghệ khoan laser hoặc gia công tia lửa điện đa trục trong thời gian 6 đến 9 tháng, hướng tới mục tiêu duy trì hiệu suất màng khí đạt tối thiểu 0,55 trên 80% diện tích bề mặt chịu nhiệt.

Thứ hai, thiết lập hệ thống kiểm soát tỉ số phun biến thiên M từ 0,8 đến 1,0 tùy theo chế độ hoạt động của động cơ (chế độ cất cánh, bay hành trình hoặc giảm tốc). Đội ngũ kỹ sư điều khiển động lực học thực hiện trong lộ trình 12 tháng, với chỉ tiêu giảm 10% đến 15% lượng khí làm mát tiêu hao mà vẫn đảm bảo an toàn nhiệt động lực học.

Thứ ba, áp dụng lớp phủ rào cản nhiệt Ceramic đa lớp với tỉ lệ chiều dày tương đối đạt từ 0,20 đến 0,30 trên nền hợp kim thép chịu nhiệt. Doanh nghiệp vật liệu hàng không phối hợp cùng viện nghiên cứu triển khai trong 18 tháng, nhằm hạ nhiệt độ bề mặt kim loại nền thêm ít nhất 130 K, kéo dài niên hạn sử dụng của chi tiết máy lên 35%.

Thứ tư, áp dụng cấu hình màng khí tối ưu với 3 hàng lỗ bố trí so le, bước lỗ tương đối P trên d bằng 3,0 và khoảng cách giữa các hàng bằng 5d. Bộ phận thiết kế khí động tuabin hoàn thiện bộ thông số hình học trong 6 tháng, đảm bảo triệt tiêu các điểm tập trung nhiệt độ cục bộ dọc theo chiều dài dòng chảy.

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

Nội dung và kết quả nghiên cứu trong luận văn mang lại giá trị thực tiễn to lớn cho bốn nhóm đối tượng chuyên ngành:

Nhóm thứ nhất là các kỹ sư thiết kế khí động và truyền nhiệt trong ngành hàng không vũ trụ và năng lượng tuabin khí. Tài liệu cung cấp cơ sở dữ liệu định lượng chính xác về góc nghiêng và tỉ số phun, phục vụ trực tiếp cho công tác tối ưu hóa hệ thống làm mát cánh tầng đầu tuabin.

Nhóm thứ hai là các nhà nghiên cứu, giảng viên chuyên ngành Kỹ thuật Động cơ, Máy và Thiết bị Thủy khí. Luận văn đóng vai trò là tài liệu tham khảo chuẩn mực về phương pháp luận mô hình hóa số CFD, thiết lập lưới phân tích truyền nhiệt liên hợp và phương pháp xử lý mô hình rối trong dòng chảy phức tạp.

Nhóm thứ ba là học viên cao học và sinh viên đại học ngành Kỹ thuật Cơ khí Động lực, Hàng không và Năng lượng Nhiệt. Công trình là ví dụ điển hình về việc ứng dụng phần mềm ANSYS Fluent giải quyết bài toán nhiệt - thủy khí thực tế, hỗ trợ đắc lực cho việc thực hiện đồ án tốt nghiệp và luận văn chuyên sâu.

Nhóm thứ tư là các chuyên gia nghiên cứu phát triển vật liệu chịu nhiệt cao cấp và công nghệ màng phủ cách nhiệt Ceramic. Luận văn cung cấp bức tranh rõ nét về tương tác truyền nhiệt giữa pha khí và kết cấu vật liệu đa tầng, giúp định hướng nâng cao độ bền cho lớp phủ TBC trong môi trường nhiệt độ khắc nghiệt.

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

Tại sao góc nghiêng lỗ 35 độ lại mang lại hiệu quả làm mát cao hơn góc 45 độ và 60 độ?

Góc nghiêng 35 độ giúp tia khí phun ra có thành phần vận tốc tiếp tuyến lớn, giữ dòng chảy luôn bám sát mặt phẳng vách. Ở góc 45 độ và 60 độ, thành phần vận tốc pháp tuyến mạnh đẩy tia khí bật sâu vào dòng khí nóng, làm suy giảm hiệu suất bảo vệ nhiệt từ 19,5% đến 36,2% tại vị trí ngay sau lỗ phun.

Tỉ số phun M bằng bao nhiêu là tối ưu nhất cho hệ thống làm mát màng khí?

Dải tỉ số phun tối ưu dao động từ M bằng 0,5 đến 1,0. Trong khoảng này, động năng của tia khí vừa đủ để dàn trải đều trên bề mặt mà không gây xoáy cuộn mạnh. Khi M vượt quá 1,2, tia phun bị tách dòng cục bộ, làm nhiệt độ bề mặt tăng lên tại vùng x từ 1d đến 5d.

Lớp phủ gốm Ceramic có tác động cụ thể như thế nào đến độ bền nhiệt của kim loại nền?

Lớp phủ Ceramic đóng vai trò là vật liệu cản nhiệt có hệ số dẫn nhiệt rất thấp. Khi đạt tỉ lệ chiều dày 0,25 so với tổng chiều dày, lớp gốm này hấp thụ và làm suy giảm nhiệt lượng truyền vào tấm thép, giúp hạ nhiệt độ lớp nền khoảng 135 K và giảm 27,4% ứng suất nhiệt phát sinh.

Mô hình dòng rối nào được chứng minh là phù hợp nhất cho bài toán màng khí làm mát?

Mô hình Realizable k-epsilon kết hợp hàm truyền nhiệt sát vách nâng cao được chứng minh là chính xác nhất. Mô hình này mô phỏng chân thực sự khuếch tán của cặp xoáy ngược chiều và trường phân bố vận tốc 3 chiều với độ lệch dưới 4,8% so với số liệu thực nghiệm.

Làm thế nào để áp dụng kết quả từ mô hình tấm phẳng vào cánh tuabin thực tế có độ cong?

Dữ liệu từ mô hình tấm phẳng cung cấp các định luật cơ bản về tỉ số phun và góc nghiêng tối ưu. Khi đưa vào cánh tuabin thực tế có độ cong và gradient áp suất, các kỹ sư chỉ cần hiệu chỉnh thông số góc nghiêng cục bộ dựa trên bán kính cong và áp suất dòng khí cục bộ theo các biểu thức đại số đã xây dựng.

Kết luận

Luận văn đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 kết luận và đóng góp học thuật nổi bật:

  • Xác lập chính xác quy luật khí động học màng khí, khẳng định góc nghiêng 35 độ cho hiệu suất bảo vệ nhiệt cao nhất đạt 0,68 tại vị trí sát lỗ.
  • Lượng hóa giới hạn tới hạn của tỉ số phun khối lượng trong dải M từ 0,5 đến 1,0 nhằm ngăn chặn hiện tượng tách dòng tia lạnh.
  • Chứng minh cấu hình 3 hàng lỗ so le giúp tăng cường diện tích bảo vệ nhiệt thêm 58% và nâng hiệu suất màng khí ở vùng hạ lưu lên 41,8%.
  • Đánh giá toàn diện cơ chế truyền nhiệt liên hợp, khẳng định lớp phủ Ceramic giúp hạ nhiệt độ kim loại nền 135 K và giảm hơn 27% ứng suất nhiệt phá hủy.
  • Xây dựng thành công hệ thống biểu thức đại số đặc trưng liên hệ giữa hiệu suất làm mát, vận tốc dòng chảy và chiều dày màng bảo vệ theo khoảng cách không thứ nguyên.

Kế hoạch phát triển tiếp theo trong 12 đến 24 tháng tới là mở rộng mô hình tính toán trên bề mặt cong 3 chiều cánh tuabin có xét đến lực ly tâm và trường xoáy quay thực tế. Các viện nghiên cứu, trường đại học và doanh nghiệp hàng không quan tâm đến mô hình tính toán và cơ sở dữ liệu số chuyên sâu hãy liên hệ Bộ môn Máy và Thiết bị Thủy khí - Trường Đại học Bách khoa Hà Nội để thúc đẩy hợp tác nghiên cứu và chuyển giao công nghệ làm mát tiên tiến.