Tổng quan nghiên cứu

Trong kỷ nguyên công nghệ không gian hiện đại, các dòng vệ tinh nhỏ như MicroSat và CubeSat đang trở thành xu hướng trọng điểm nhờ chi phí chế tạo thấp, thời gian phát triển ngắn và hiệu quả kinh tế vượt bậc. Tuy nhiên, kích thước nhỏ gọn cùng lớp vỏ mỏng khiến các vệ tinh này rất nhạy cảm với sự biến đổi nhiệt khắc nghiệt trên quỹ đạo tầm thấp (LEO). Khi hoạt động ở độ cao 600 km với chu kỳ quỹ đạo 5800 giây, vệ tinh liên tục chuyển tiếp giữa vùng chiếu sáng bởi bức xạ Mặt Trời có cường độ khoảng 1360 W/m² và vùng bóng tối của Trái Đất. Sự chênh lệch nhiệt độ lớn này đe dọa trực tiếp đến tuổi thọ của các linh kiện điện tử bên trong và các tấm pin năng lượng mặt trời bên ngoài.

Vấn đề nghiên cứu cốt lõi được đặt ra là làm thế nào để dự đoán chính xác, nhanh chóng ứng xử nhiệt phi tuyến của vệ tinh mà không phụ thuộc hoàn toàn vào các phần mềm mô phỏng số nặng nề và tốn kém thời gian. Mục tiêu cụ thể của luận văn là xây dựng mô hình nhiệt rời rạc hai nút (nút trong đại diện cho khối thiết bị điện tử, nút ngoài đại diện cho vỏ vệ tinh và tấm pin), thiết lập hệ phương trình vi phân cân bằng nhiệt phi tuyến, đồng thời phát triển phương pháp tuyến tính hóa tương đương cùng phương pháp cân bằng điều hòa để tìm nghiệm giải tích.

Nghiên cứu được triển khai trên mô hình vệ tinh giả định chuyển động quanh Trái Đất ở độ cao 600 km với tỷ số thời gian chiếu sáng đạt 63% chu kỳ. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm thiểu sai số tính toán nhiệt xuống dưới 1,5% so với phương pháp giải số Runge-Kutta bậc 4, đồng thời rút ngắn hơn 40% thời gian tính toán trong giai đoạn thiết kế sơ bộ hệ thống kiểm soát nhiệt vệ tinh.

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 dựa trên nền tảng vững chắc của nhiệt động lực học và lý thuyết dao động phi tuyến trong cơ học vật thể rắn. Hai lý thuyết vật lý nền tảng được áp dụng xuyên suốt là định luật Stefan-Boltzmann về bức xạ nhiệt và định luật Fourier về dẫn nhiệt. Theo định luật Stefan-Boltzmann, công suất bức xạ nhiệt của một bề mặt tỷ lệ thuận với lũy thừa bậc bốn của nhiệt độ tuyệt đối theo công thức $P = \varepsilon \sigma A T^4$, trong đó hằng số Stefan-Boltzmann có giá trị xấp xỉ $5,67 \times 10^{-8} \text{ W/m}^2\text{K}^4$. Bên cạnh đó, định luật Fourier xác định mật độ dòng nhiệt dẫn truyền qua bề mặt tỷ lệ với gradient nhiệt độ và hệ số dẫn nhiệt của vật liệu.

Mô hình nghiên cứu kế thừa và phát triển lý thuyết cân bằng nhiệt bức xạ trong không gian, tích hợp thêm lý thuyết tuyến tính hóa tương đương của GS.TSKH Nguyễn Đông Anh để giải quyết số hạng bức xạ phi tuyến bậc bốn. Các khái niệm chính được định nghĩa chặt chẽ bao gồm: nút đẳng nhiệt (vùng thể tích đồng nhất về nhiệt độ), hệ số hiển thị bức xạ (tỷ lệ năng lượng bức xạ rời khỏi một bề mặt bị chặn bởi bề mặt khác), suất phản chiếu Albedo của Trái Đất (trung bình đạt giá trị khoảng 0,31 đến 0,33), và hệ số tắt dần không thứ nguyên của hệ dao động nhiệt.

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

Để kiểm chứng mô hình toán học, nghiên cứu sử dụng bộ dữ liệu tham số tiêu chuẩn của một vệ tinh nhỏ hình khối với kích thước diện tích vỏ ngoài là 3,14 m² và diện tích khối thiết bị bên trong là 1,13 m². Cỡ mẫu khảo sát gồm 12 thông số vật lý then chốt, bao gồm nhiệt dung nút ngoài 30.000 J/K, nhiệt dung nút trong 20.000 J/K, công suất tiêu tán nhiệt nội bộ 50 W, hệ số dẫn nhiệt liên kết 10 W/K, và hệ số phát xạ hồng ngoại đạt 0,83. Phương pháp chọn mẫu tham số dựa trên các tiêu chuẩn thiết kế thực tế của các dòng vệ tinh quỹ đạo thấp đang hoạt động nhằm đảm bảo tính đại diện cao nhất.

Quy trình phân tích dữ liệu được thực hiện thông qua 4 phương pháp tính toán độc lập: phương pháp giải số Runge-Kutta bậc 4 (đóng vai trò làm chuẩn đối sánh), phương pháp giải tích tuyến tính hóa của Grande, phương pháp tuyến tính hóa tương đương, và phương pháp cân bằng điều hòa kết hợp khai triển chuỗi Fourier. Lý do lựa chọn kết hợp các phương pháp này là nhằm tìm kiếm một lời giải giải tích tường minh, giúp các kỹ sư nắm bắt bản chất vật lý của quá trình truyền nhiệt mà không cần lặp số phức tạp. Toàn bộ quá trình mô phỏng được thực hiện trên khung thời gian 10 chu kỳ quỹ đạo liên tiếp (tương đương 58.000 giây) để khảo sát toàn diện quá trình chuyển tiếp và trạng thái bình ổn nhiệt.

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

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

Phân tích số trị và giải tích từ luận văn đã đem lại 4 phát hiện quan trọng mang tính đột phá:

Thứ nhất, hệ nhiệt của vệ tinh nhanh chóng đạt trạng thái bình ổn nhiệt động học sau khoảng 5 đến 7 chu kỳ quỹ đạo (tương đương 29.000 đến 40.600 giây). Biên độ dao động nhiệt của nút ngoài diễn ra rất mạnh với mức chênh lệch nhiệt độ ngày - đêm lên tới khoảng 28 K đến 34 K. Trong khi đó, nhiệt độ nút trong biến thiên rất êm với biên độ dao động được triệt tiêu hơn 65%, chỉ dao động trong khoảng 4 K đến 7 K nhờ khả năng tích nhiệt của nhiệt dung 20.000 J/K.

Thứ hai, phương pháp tuyến tính hóa tương đương thể hiện sự chính xác vượt trội khi so sánh với phương pháp giải số Runge-Kutta bậc 4. Độ sai lệch cực đại về nhiệt độ giữa hai phương pháp này luôn được khống chế ở mức dưới 1,2%, tốt hơn đáng kể so với phương pháp tuyến tính hóa cổ điển của Grande (vốn có sai số khoảng 3,5% đến 4,8% tại các điểm cực trị).

Thứ ba, nghiên cứu xác định được tỷ số nhiệt dung giữa nút ngoài và nút trong ở mức 1,5 cùng hệ số liên kết nhiệt 10 W/K là ngưỡng tối ưu để giữ cho thiết bị bên trong luôn duy trì nhiệt độ trong dải an toàn từ 285 K đến 298 K (tương đương 12 °C đến 25 °C).

Thứ tư, tồn tại một độ trễ pha nhiệt học rõ rệt giữa dao động của nguồn bức xạ Mặt Trời đầu vào và phản ứng nhiệt của nút trong. Độ lệch pha đo được xấp xỉ 15 độ đến 18 độ trên vòng quỹ đạo, phản ánh quán tính nhiệt của kết cấu bên trong vệ tinh.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự khác biệt về biên độ nhiệt giữa hai nút là do nút ngoài tiếp xúc trực tiếp với các dòng nhiệt bức xạ biến thiên tuần hoàn theo dạng sóng vuông từ Mặt Trời và sóng điều hòa từ bức xạ Albedo Trái Đất. Ngược lại, nút trong được bao bọc kín và đóng vai trò như một bộ lọc thông thấp nhiệt học, chỉ tiếp nhận nhiệt gián tiếp qua dẫn nhiệt và bức xạ bên trong. Khi so sánh với các nghiên cứu trước đây về mô hình một nút đẳng nhiệt, mô hình hai nút trong luận văn này đã phản ánh chân thực hơn hiện tượng chênh lệch nhiệt độ giữa vỏ và lõi vệ tinh.

Để biểu diễn trực quan các kết quả này, toàn bộ quá trình biến thiên nhiệt độ có thể được trình bày thông qua biểu đồ không gian pha hai chiều thể hiện quỹ đạo khép kín (vòng giới hạn nhiệt động học giữa nhiệt độ nút ngoài và nút trong). Đồng thời, một bảng so sánh chi tiết nhiệt độ cực đại, cực tiểu và giá trị trung bình trên 10 chu kỳ quỹ đạo giữa bốn phương pháp tính toán giúp các nhà nghiên cứu dễ dàng nhận diện độ hội tụ và độ chính xác của từng giải thuật. Kết quả này cung cấp một công cụ phân tích giải tích mạnh mẽ cho giai đoạn thiết kế sơ bộ cấu trúc vệ tinh.

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

Dựa trên các kết quả nghiên cứu đã đạt được, luận văn đưa ra 4 nhóm giải pháp kỹ thuật cụ thể nhằm tối ưu hóa quá trình thiết kế và kiểm soát nhiệt cho vệ tinh nhỏ:

Thứ nhất, tích hợp phương pháp tuyến tính hóa tương đương vào phần mềm thiết kế sơ bộ hệ thống nhiệt. Đội ngũ kỹ sư thiết kế nhiệt thuộc các viện nghiên cứu hàng không vũ trụ nên áp dụng giải thuật này trong 3 tháng đầu của giai đoạn định hình dự án nhằm rút ngắn 45% thời gian tính toán mô phỏng so với việc sử dụng ngay các mô hình phần tử hữu hạn nặng nề.

Thứ hai, tối ưu hóa tỷ lệ diện tích bề mặt phủ pin quang điện và lớp sơn tản nhiệt trên vỏ ngoài vệ tinh. Nhóm phát triển vật liệu kết cấu cần duy trì tỷ lệ 85% diện tích pin mặt trời (hệ số hấp thụ 0,82) kết hợp 15% diện tích sơn trắng quang học (hệ số hấp thụ 0,88 và hệ số phát xạ 0,83) trong vòng 6 tháng thử nghiệm mẫu để đảm bảo nhiệt độ cân bằng trung bình của vệ tinh ổn định quanh mức 293 K (20 °C).

Thứ ba, thiết kế và lắp đặt hệ thống cách nhiệt đa lớp (Multi-Layer Insulation - MLI) kết hợp cầu truyền nhiệt có kiểm soát giữa vỏ và khối thiết bị. Đội ngũ kỹ sư cơ điện tử cần hoàn thiện cấu trúc dẫn nhiệt với hệ số truyền dẫn mục tiêu đạt 10 W/K trong thời gian 4 tháng trước khi tích hợp toàn hệ thống, nhằm hạn chế biên độ dao động nhiệt của khối xử lý trung tâm xuống dưới 5 K.

Thứ tư, mở rộng phát triển mô hình tính toán nhiệt từ hệ hai nút lên hệ nhiều nút không gian ba chiều (từ 30 đến 50 nút). Nhóm nghiên cứu thuộc các trường đại học chuyên ngành kỹ thuật cần tiếp tục hợp tác triển khai dự án nghiên cứu mở rộng trong 12 tháng tiếp theo để xây dựng thư viện tính toán nhiệt giải tích hoàn chỉnh cho các dòng vệ tinh khối lập phương CubeSat từ 1U đến 6U.

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

Luận văn thạc sĩ này là tài liệu học thuật và kỹ thuật giá trị cao cho 4 nhóm đối tượng chuyên môn:

Nhóm thứ nhất là các kỹ sư thiết kế hệ thống và kỹ sư nhiệt hàng không vũ trụ tại các trung tâm vũ trụ và doanh nghiệp công nghệ không gian. Nhóm này có thể sử dụng trực tiếp khung tính toán giải tích của mô hình hai nút để đánh giá nhanh phân bố nhiệt và lựa chọn vật liệu kiểm soát nhiệt thụ động cho các dự án vệ tinh nhỏ.

Nhóm thứ hai là các nhà khoa học, giảng viên và nghiên cứu viên chuyên ngành Cơ học vật thể rắn, Cơ học công trình và Vật lý kỹ thuật. Luận văn cung cấp phương pháp luận mẫu mực trong việc áp dụng kỹ thuật tuyến tính hóa tương đương vào việc giải quyết các hệ phương trình vi phân phi tuyến bậc cao có số hạng bức xạ nhiệt.

Nhóm thứ ba là học viên cao học và sinh viên các ngành Kỹ thuật Không gian, Kỹ thuật Hàng không, Cơ điện tử và Toán ứng dụng. Đây là tài liệu tham khảo bổ ích giúp người học nắm vững phương pháp mô hình hóa bài toán truyền nhiệt phức hợp (dẫn nhiệt kết hợp bức xạ) trong môi trường chân không ngoài không gian.

Nhóm thứ tư là các nhà quản lý dự án công nghệ cao và chuyên gia thẩm định kỹ thuật vệ tinh. Tài liệu cung cấp cơ sở lý thuyết chuẩn xác để đánh giá tính khả thi và độ tin cậy của các giải pháp kiểm soát nhiệt trước khi tiến hành các thử nghiệm nhiệt chân không (TVAC) tốn kém.

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

Tại sao việc phân tích nhiệt lại đặc biệt quan trọng đối với vệ tinh nhỏ ở quỹ đạo thấp?

Vệ tinh hoạt động ở quỹ đạo thấp (khoảng 600 km) phải chịu sự thay đổi nhiệt độ chu kỳ liên tục giữa vùng nhận bức xạ Mặt Trời 1360 W/m² và vùng bóng tối Trái Đất trong chu kỳ 5800 giây. Vỏ vệ tinh nhỏ rất mỏng nên nếu không phân tích nhiệt chính xác, linh kiện điện tử bên trong có thể bị quá nhiệt hoặc đóng băng, dẫn đến hỏng hóc toàn bộ sứ mệnh.

Mô hình nhiệt hai nút có đảm bảo độ tin cậy so với mô hình liên tục không?

Mô hình nhiệt hai nút hoàn toàn đảm bảo độ tin cậy trong giai đoạn thiết kế sơ bộ. Bằng cách gom cụm thành nút ngoài và nút trong đẳng nhiệt, mô hình cho kết quả sai lệch dưới 1,2% so với giải thuật số Runge-Kutta bậc 4, giúp nắm bắt chính xác xu hướng biến thiên nhiệt mà không đòi hỏi chi phí tính toán lớn.

Điểm khác biệt cốt lõi giữa phương pháp tuyến tính hóa tương đương và phương pháp Grande là gì?

Phương pháp Grande tuyến tính hóa số hạng bức xạ quanh nhiệt độ trung bình tĩnh, trong khi phương pháp tuyến tính hóa tương đương tối thiểu hóa sai số năng lượng trung bình của hệ phi tuyến trong cả chu kỳ dao động. Do đó, phương pháp tuyến tính hóa tương đương đạt độ chính xác cao hơn rõ rệt tại các vùng chuyển tiếp nhiệt độ cực trị.

Nhiệt độ của vệ tinh biến thiên như thế nào khi đi vào vùng bóng tối của Trái Đất?

Khi đi vào vùng bóng tối (chiếm khoảng 37% chu kỳ quỹ đạo), tải nhiệt bức xạ Mặt Trời và Albedo lập tức bằng 0. Nhiệt độ nút ngoài giảm nhanh chóng từ khoảng 305 K xuống 275 K, trong khi nhiệt độ nút trong giảm chậm hơn nhiều nhờ quán tính nhiệt và tiếp tục được sưởi ấm bởi công suất tiêu tán 50 W nội bộ.

Làm thế nào để kiểm soát nhiệt độ linh kiện bên trong vệ tinh một cách thụ động?

Kiểm soát nhiệt thụ động được thực hiện bằng cách lựa chọn vật liệu phủ bề mặt với hệ số hấp thụ mặt trời 0,67 và độ phát xạ hồng ngoại 0,83, kết hợp điều chỉnh hệ số dẫn nhiệt giữa các vách ngăn ở mức 10 W/K. Giải pháp này giúp duy trì nhiệt độ khối linh kiện ổn định quanh mức 20 °C mà không cần tiêu tốn điện năng.

Kết luận

Luận văn thạc sĩ đã giải quyết trọn vẹn và xuất sắc bài toán phân tích nhiệt cho vệ tinh nhỏ hoạt động trên quỹ đạo tầm thấp thông qua 5 kết quả cốt lõi:

  • Thiết lập thành công mô hình toán học hai nút rời rạc phản ánh đầy đủ các cơ chế truyền nhiệt dẫn nhiệt và bức xạ nhiệt trong môi trường không gian.
  • Ứng dụng thành công phương pháp tuyến tính hóa tương đương để giải hệ phương trình vi phân nhiệt phi tuyến, đạt độ sai lệch dưới 1,2% so với nghiệm số Runge-Kutta bậc 4.
  • Làm sáng tỏ quy luật hình thành vòng giới hạn nhiệt động học trong không gian pha và xác định độ trễ pha dao động nhiệt giữa vỏ ngoài và khối thiết bị bên trong.
  • Xác định được bộ thông số nhiệt dung và hệ số dẫn nhiệt tối ưu nhằm duy trì môi trường nhiệt độ an toàn từ 285 K đến 298 K cho các thiết bị điện tử.
  • Cung cấp một phương pháp luận giải tích nhanh chóng, tiết kiệm tài nguyên tính toán và có độ tin cậy cao cho các kỹ sư cơ học hàng không vũ trụ.

Đóng góp chính của công trình là tạo tiền đề lý thuyết vững chắc cho các chương trình phát triển vệ tinh nhỏ của Việt Nam trong giai đoạn 2026 - 2030, mở ra khả năng ứng dụng cho các hệ nhiệt nhiều nút phức tạp hơn. Các viện nghiên cứu, trường đại học và kỹ sư chuyên ngành hãy áp dụng ngay khung tính toán này vào các dự án thiết kế vệ tinh thực tế để tối ưu hóa thời gian và nâng cao độ an toàn cho các sứ mệnh không gian tương lai.