Tổng quan nghiên cứu

Trong kỹ thuật hàng không hiện đại, các mẫu máy bay chiến đấu và phương tiện bay không người lái thường xuyên sử dụng cấu hình cánh tam giác (cánh Delta) nhằm đạt hiệu suất cơ động tối ưu ở dải vận tốc siêu âm. Tuy nhiên, mọi chuyến bay đều bắt buộc phải trải qua giai đoạn cất cánh và hạ cánh ở dải tốc độ thấp, thông thường dao động từ 180 km/h đến 260 km/h. Tại dải vận tốc này, khi máy bay hoạt động ở góc tấn lớn từ 15 độ đến 30 độ, dòng chảy xung quanh cánh trở nên phi tuyến và mất ổn định nghiêm trọng, làm gia tăng nguy cơ mất điều khiển lên tới 35% so với chế độ bay bằng.

Vấn đề nghiên cứu cốt lõi của đề tài tập trung vào việc phân tích đặc tính khí động học của cánh Delta ở tốc độ thấp, đặc biệt là cơ chế hình thành, phát triển và tan rã của cặp cuộn xoáy trên mặt lưng cánh khi có sự biến đổi về hình dạng hình học. Mục tiêu cụ thể của luận văn là làm sáng tỏ ảnh hưởng của góc quét mép trước và tỷ số dạng cánh đến trường phân bố áp suất, trường vận tốc và các hệ số lực nâng, lực cản; từ đó xây dựng cơ sở khoa học phục vụ thiết kế tối ưu hóa cánh.

Nghiên cứu được triển khai trong thời gian 6 tháng tại Trường Cơ khí thuộc Đại học Bách Khoa Hà Nội, kết hợp chặt chẽ giữa tính toán mô phỏng số CFD trên phần mềm ANSYS Fluent và đo đạc kiểm chứng thực nghiệm trong ống khí động dưới âm. Kết quả nghiên cứu có ý nghĩa thực tiễn to lớn trong việc nâng cao hệ số lực nâng cực đại thêm khoảng 15% đến 20%, hỗ trợ các kỹ sư kiểm soát hiện tượng vỡ xoáy và nâng cao độ an toàn cho phương tiện bay trong giai đoạn hạ cánh.

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 dựa trên nền tảng của hệ phương trình Navier-Stokes trung bình Reynolds (RANS) mô tả dòng lưu chất nhớt chịu nén và không nén được, kết hợp với các mô hình đóng kín trường rối hiện đại như mô hình Shear Stress Transport k-omega (SST k-omega) và mô hình Spalart-Allmaras. Bên cạnh đó, luận văn ứng dụng lý thuyết tương đương lực hút mép trước của Polhamus, chia lực nâng tổng cộng của cánh Delta thành hai thành phần cơ bản: lực nâng thế tuyến tính gắn liền với dòng dính thông thường và lực nâng phi tuyến sinh ra từ độ xoáy của cặp xoáy mép trước.

Khung lý thuyết của đề tài xoay quanh các khái niệm then chốt trong khí động học phi tuyến:

  1. Lực nâng xoáy (Vortex Lift): Lực nâng phát sinh do vùng áp suất cực thấp hình thành dưới lõi của hai cuộn xoáy bám dọc mép trước trên mặt lưng cánh.
  2. Hiện tượng vỡ xoáy (Vortex Breakdown): Sự suy sụp đột ngột của cấu trúc dòng xoáy khi gặp građien áp suất nghịch, dẫn đến hiện tượng phân ly dòng và mất lực nâng cục bộ.
  3. Góc quét mép trước (Leading-Edge Sweep Angle): Góc hợp bởi mép trước của cánh và trục ngang thân máy bay, thường nằm trong khoảng từ 50 độ đến 70 độ.
  4. Tỷ số dạng cánh (Aspect Ratio) và tỷ số thon: Các thông số hình học quyết định trực tiếp đến cường độ cuộn xoáy và hiệu suất khí động toàn phần.
  5. Số Reynolds (Re): Đại lượng không thứ nguyên đặc trưng cho tỷ số giữa lực quán tính và lực nhớt trong môi trường khí quyển, đạt giá trị khoảng 2,5 x 10^5 trong các kịch bản thử nghiệm.

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

Quy trình nghiên cứu sử dụng phương pháp mô phỏng số động lực học dòng chảy (CFD) trên nền tảng phần mềm chuyên dụng ANSYS Fluent, được kiểm chứng đối soát thông qua hệ thống thực nghiệm ống khí động dưới âm.

Cỡ mẫu nghiên cứu bao gồm 5 cấu hình cánh Delta có sự biến thiên góc quét mép trước từ 50 độ đến 70 độ và thay đổi độ vát mép cắt. Phương pháp chọn mẫu là chọn mẫu kỹ thuật có chủ đích theo dải thông số khí động học chuẩn quốc tế nhằm cô lập biến số hình dạng. Lý do lựa chọn mô hình mô phỏng CFD kết hợp mô hình rối SST k-omega là khả năng nắm bắt chính xác hiện tượng phân tách dòng khí động phức tạp và sự phát triển của xoáy cuộn ba chiều sát bề mặt mà các công thức giải tích đơn giản không thể giải quyết.

Miền tính toán 3D được số hóa bằng hệ thống lưới đa khối với số lượng phần tử dao động từ 2,2 triệu đến 3,8 triệu phần tử. Lớp biên sát bề mặt cánh được phân chia tinh vi với chỉ số khoảng cách không thứ nguyên y+ duy trì dưới 1,0 nhằm nắm bắt trọn vẹn lớp lót nhớt. Thực nghiệm kiểm chứng được tiến hành tại ống khí động hở với vận tốc dòng thổi từ 15 m/s đến 30 m/s, sử dụng hệ thống đo áp suất vi sai đa kênh gắn dọc các mặt cắt dây cung x/c = 0,5 và x/c = 0,7. Toàn bộ tiến trình nghiên cứu kéo dài 6 tháng, phân chia thành 3 giai đoạn: 1 tháng nghiên cứu tổng quan lý thuyết, 3 tháng xây dựng mô hình số và mô phỏng chuyên sâu, 2 tháng tiến hành đo thực nghiệm và hiệu chỉnh kết quả.

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

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

Nghiên cứu đã ghi nhận và định lượng các quy luật khí động học quan trọng của dòng chảy qua cánh Delta ở vận tốc thấp:

Thứ nhất, cặp cuộn xoáy mép trước đóng vai trò quyết định trong việc tạo lực nâng phi tuyến, đóng góp tới 42,5% tổng hệ số lực nâng của cánh khi góc tấn đạt mốc 20 độ. Khi góc tấn tăng dần từ 0 độ đến 20 độ, hệ số lực nâng tăng tuyến tính mạnh mẽ từ 0,08 lên mức cực đại 1,15.

Thứ hai, hiện tượng vỡ xoáy bắt đầu xuất hiện tại góc tấn 25 độ, khởi phát từ vị trí dây cung tương đối x/c = 0,65 và nhanh chóng lan ngược về phía mép sau của cánh. Sự tan rã cấu trúc xoáy này làm sụt giảm 22,4% đỉnh áp suất hút tại mặt lưng cánh, trực tiếp gây ra hiện tượng rung lắc khí động (wing rock).

Thứ ba, việc điều chỉnh hình dạng cánh từ dạng tam giác đơn thuần góc quét 50 độ sang dạng cánh có góc quét 65 độ giúp trì hoãn hiện tượng vỡ xoáy, dời điểm khởi phát vỡ xoáy về phía sau thêm 18% chiều dài dây cung. Điều này cho phép mở rộng dải góc tấn hoạt động an toàn từ 22 độ lên tới 28 độ trước khi xảy ra hiện tượng thất tốc hoàn toàn.

Thứ tư, kết quả so sánh đối chiếu giữa mô phỏng số và thực nghiệm đo đạc áp suất tại mặt cắt x/c = 0,5 cho thấy độ tin cậy vượt trội của mô hình tính toán, với sai số trung bình chỉ dao động trong biên độ từ 3,6% đến 5,8%, hoàn toàn nằm trong giới hạn cho phép của các tiêu chuẩn khí động học công nghiệp.

Thảo luận kết quả

Nguyên nhân cốt lõi của việc duy trì lực nâng vượt trội ở góc tấn cao là do sự hình thành gradient áp suất âm cực lớn tại tâm cuộn xoáy, tạo ra lực hút mạnh ép chặt dòng khí lên bề mặt cánh. Tuy nhiên, khi góc tấn vượt qua ngưỡng 25 độ, gradient áp suất nghịch dọc theo trục dòng xoáy vượt quá năng lượng động lực của lõi xoáy, dẫn đến việc lõi xoáy bị phình to và vỡ tung.

Khi so sánh với các nghiên cứu kinh điển về cánh tam giác trong nước và quốc tế, kết quả của luận văn hoàn toàn tương đồng về mặt xu hướng khí động học nhưng cung cấp bức tranh dữ liệu chi tiết hơn về trường vận tốc tức thời ở dải vận tốc thấp dưới 30 m/s.

Trong các báo cáo kỹ thuật, toàn bộ hệ thống dữ liệu này có thể được biểu diễn trực quan thông qua đồ thị phân bố hệ số áp suất Cp theo phương sải cánh tại các mặt cắt x/c = 0,3; 0,5; 0,7 và 0,9; kết hợp cùng biểu đồ vector vận tốc thể hiện vị trí mắt xoáy và bảng đối chiếu tỷ số lực nâng trên lực cản (L/D) giữa các phương án hình học. Những kết quả này chứng minh rằng việc biến đổi hình dạng mép trước cánh là giải pháp cơ học thụ động hiệu quả nhất để kiểm soát dòng xoáy mà không làm tăng khối lượng kết cấu máy bay.

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

Nhằm ứng dụng hiệu quả các kết quả nghiên cứu vào thực tiễn thiết kế phương tiện bay, các giải pháp và khuyến nghị kỹ thuật cụ thể bao gồm:

  1. Tối ưu hóa cấu hình hình học cánh tam giác kép (Double-Delta Wing): Nhóm kỹ sư thiết kế khí động học cần áp dụng cấu hình cánh với góc quét mép trước gốc cánh đạt 70 độ và góc quét đầu cánh đạt 50 độ. Giải pháp này giúp gia tăng góc tấn tới hạn thêm từ 3 độ đến 5 độ, nâng cao khả năng điều khiển phương tiện khi hạ cánh, triển khai trong lộ trình 12 tháng.

  2. Tích hợp công nghệ kiểm soát dòng chủ động (Active Flow Control): Viện nghiên cứu và phát triển công nghệ hàng không nên tích hợp các khe thổi khí vi mô dọc theo 40% chiều dài mép trước cánh. Mục tiêu là cung cấp thêm động năng cho lõi xoáy, giúp trì hoãn hiện tượng vỡ xoáy ở dải tốc độ dưới 40 m/s và tăng hệ số lực nâng thêm 15% trong vòng 18 tháng thử nghiệm.

  3. Chuẩn hóa quy trình tính toán khí động số CFD đa tầng: Các đơn vị nghiên cứu kỹ thuật cơ khí động lực cần chuẩn hóa quy trình phân chia lưới lớp biên với tiêu chuẩn kiểm soát y+ nhỏ hơn hoặc bằng 1,0 kết hợp mô hình SST k-omega. Quy trình này giúp giảm thiểu 25% thời gian tính toán mô phỏng mà vẫn bảo đảm sai số dưới 5%, thực hiện ngay trong vòng 6 tháng tới.

  4. Nâng cấp hệ thống đo lường quang học trong hầm gió thí nghiệm: Các trung tâm nghiên cứu cần trang bị thêm hệ thống đo vận tốc trường bằng hình ảnh hạt (Particle Image Velocimetry - PIV) nhằm trực quan hóa không gian ba chiều của cấu trúc xoáy với độ chính xác đạt trên 98%, hoàn thành trong lộ trình 24 tháng.

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

Nội dung và dữ liệu của luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cho 4 nhóm đối tượng chính:

  1. Kỹ sư thiết kế khí động học hàng không: Tiếp cận bộ thông số tối ưu về tỷ số dạng và góc quét cánh Delta để áp dụng vào các dự án phát triển máy bay tiêm kích phản lực, máy bay huấn luyện siêu âm và thiết bị bay không người lái (UAV) cơ động cao.

  2. Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Cơ khí Động lực: Sử dụng luận văn như tài liệu tham khảo tiêu chuẩn về phương pháp luận kết hợp mô phỏng số ANSYS Fluent với đo đạc thực nghiệm trong ống khí động hở, phục vụ các đề tài cơ học chất lưu phức tạp.

  3. Kỹ sư khí động lực học ô tô thể thao và xe đua chuyên nghiệp: Ứng dụng nguyên lý tạo lực nâng xoáy đảo ngược để thiết kế cánh gió trước và bộ khuếch tán gầm xe dạng tam giác, nhằm tối ưu hóa lực ép mặt đường (downforce) thêm khoảng 10% đến 18% khi xe vào cua ở tốc độ trên 150 km/h.

  4. Giảng viên và nhà nghiên cứu tại các trường đại học kỹ thuật: Khai thác làm học liệu giảng dạy chuyên sâu cho các học phần Khí động học ứng dụng, Động lực học chất lưu tính toán (CFD) và Kỹ thuật thí nghiệm hàng không.

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

  1. Tại sao cánh Delta lại hình thành hai cuộn xoáy lớn ở mặt trên khi bay ở tốc độ thấp? Khi dòng khí di chuyển qua mép trước sắc nhọn của cánh Delta ở một góc tấn nhất định, dòng chảy ở mặt dưới có áp suất cao hơn sẽ tràn qua mép trước lên mặt trên có áp suất thấp hơn. Sự chênh lệch áp suất này cuộn dòng khí lại thành hai cuộn xoáy mạnh mẽ chạy dọc theo mép cánh về phía đuôi, tạo ra vùng áp suất âm sâu giúp nâng đỡ máy bay.

  2. Hiện tượng vỡ xoáy (Vortex Breakdown) gây ra những nguy hiểm gì cho máy bay? Hiện tượng vỡ xoáy làm cấu trúc xoáy có tổ chức bị phân rã đột ngột thành dòng hỗn loạn, làm mất hoàn toàn vùng áp suất âm trên cánh. Hậu quả là máy bay bị sụt giảm hơn 20% lực nâng tức thời, gia tăng lực cản và xuất hiện hiện tượng rung lắc khí động bất đối xứng, dễ dẫn đến mất kiểm soát thăng bằng ở chế độ hạ cánh.

  3. Phương pháp mô phỏng ANSYS Fluent trong đề tài có độ tin cậy như thế nào so với thực nghiệm? Độ tin cậy của mô hình số được kiểm chứng qua sự trùng khớp cao với dữ liệu thực nghiệm hầm gió, với sai số đo phân bố áp suất chỉ từ 3,6% đến 5,8%. Việc sử dụng mô hình rối SST k-omega và lưới tính toán có y+ dưới 1,0 giúp mô phỏng chính xác sự biến đổi của dòng xoáy phức tạp.

  4. Vì sao việc nghiên cứu cánh Delta ở dải vận tốc thấp lại quan trọng dù máy bay vốn hoạt động ở tốc độ siêu âm? Mặc dù cánh Delta được thiết kế tối ưu cho chế độ bay siêu thanh, máy bay vẫn bắt buộc phải cất cánh và hạ cánh ở dải tốc độ thấp từ 180 km/h đến 260 km/h. Đây là giai đoạn máy bay phải bay ở góc tấn cao, dòng khí mất ổn định nhất nên việc nghiên cứu giúp đảm bảo an toàn tuyệt đối cho phi hành đoàn và phương tiện.

  5. Việc thay đổi góc quét mép trước cánh mang lại lợi ích cụ thể nào cho đặc tính khí động? Khi tăng góc quét mép trước từ 50 độ lên 65 độ, cường độ cuộn xoáy được duy trì ổn định hơn và điểm khởi phát hiện tượng vỡ xoáy bị đẩy lùi về phía sau thêm 18% chiều dài cánh. Điều này giúp nâng cao góc tấn thất tốc thêm 6 độ và duy trì lực nâng ổn định trong các thao tác cơ động ngặt nghèo.

Kết luận

  • Luận văn đã phân tích toàn diện cơ chế hình thành và biến đổi của cặp cuộn xoáy mép trước trên cánh Delta ở dải vận tốc thấp dưới tác động của sự thay đổi hình dạng cánh.
  • Xây dựng thành công mô hình mô phỏng khí động số CFD với độ chính xác cao, kiểm chứng sai số thực nghiệm duy trì ổn định dưới 5,8% tại các mặt cắt chuẩn.
  • Xác định chính xác ngưỡng góc tấn tới hạn 25 độ là điểm bùng phát hiện tượng vỡ xoáy làm suy giảm 22,4% đỉnh áp suất hút trên bề mặt cánh.
  • Đưa ra các phương án cải tiến hình học cánh giúp dời vị trí vỡ xoáy về phía sau thêm 18%, góp phần nâng cao hệ số lực nâng tổng thể thêm 15% đến 20%.
  • Mở ra nền tảng lý thuyết và số liệu thực nghiệm vững chắc để phục vụ việc thiết kế các mẫu máy bay quân sự thế hệ mới và UAV tốc độ cao.

Đóng góp lớn nhất của công trình là thiết lập hệ thống cơ sở dữ liệu khí động học chi tiết kết hợp giữa mô phỏng số và thực nghiệm hầm gió đối với cánh tam giác ở vận tốc thấp tại Việt Nam. Trong lộ trình 12 đến 24 tháng tới, các nhóm nghiên cứu nên tiếp tục mở rộng đề tài sang kiểm soát dòng chủ động và mô phỏng tương tác kết cấu - khí động học (FSI). Các kỹ sư và nhà nghiên cứu quan tâm có thể khai thác trực tiếp quy trình mô phỏng và bộ dữ liệu này để tối ưu hóa các thiết kế khí động học hàng không thực tế.