Tổng quan luận án
Bài toán phân tích và kiểm soát nhiệt đóng vai trò then chốt trong quy trình thiết kế, chế tạo và vận hành vệ tinh nhân tạo nhằm bảo đảm các thiết bị bên trong (thiết bị điện tử, quang học, pin ắc quy) hoạt động ổn định trong miền nhiệt độ giới hạn kỹ thuật. Trong môi trường không gian khắc nghiệt, vệ tinh chịu tác động đồng thời của các tải nhiệt môi trường có tính chu kỳ và phi tuyến. Đối với vệ tinh nhỏ (khối lượng dưới 500 kg, đặc biệt là nhóm dưới 100 kg hoạt động trên quỹ đạo thấp quanh Trái đất), việc phân tích nhiệt thường dựa trên các công cụ mô phỏng số đóng gói trong các phần mềm thương mại. Tuy nhiên, phương pháp số đòi hỏi thời gian tính toán kéo dài và tiêu tốn nhiều tài nguyên máy tính, gây bất lợi trong giai đoạn tiền thiết kế khi các tham số hình học, vật liệu liên tục bị thay đổi.
Ngược lại, các công cụ giải tích có ưu thế cung cấp nghiệm nhanh và làm rõ bản chất vật lý của bài toán lại rất hạn chế do sự hiện diện của số hạng bức xạ nhiệt phi tuyến bậc bốn ($T^4$). Các hướng tiếp cận giải tích trước đây thường phải dựa trên giả thiết dao động nhiệt quanh giá trị trung bình là nhỏ, điều này làm giảm độ tin cậy khi biên độ dao động nhiệt lớn (như sự chuyển tiếp giữa vùng chiếu sáng và vùng bóng tối của Trái đất). Xuất phát từ khoảng trống nghiên cứu này, luận án tiến sĩ kỹ thuật với đề tài "Nghiên cứu đáp ứng nhiệt của vệ tinh nhỏ trên quỹ đạo thấp chịu tác dụng của môi trường nhiệt vũ trụ" của tác giả Phạm Ngọc Chung (người hướng dẫn: GS.TSKH Nguyễn Đông Anh và PGS.TS Đinh Văn Mạnh) đã đề xuất việc mở rộng và áp dụng phương pháp tuyến tính hóa tương đương theo tiêu chuẩn đối ngẫu để giải bài toán cân bằng nhiệt phi tuyến của vệ tinh.
Mục tiêu nghiên cứu của luận án được xác định gồm 3 nội dung chính:
- Xây dựng các mô hình nhiệt một nút, hai nút và nhiều nút với các mô hình tải nhiệt khác nhau tác động lên vệ tinh nhỏ trên quỹ đạo thấp quanh Trái đất.
- Tìm được nghiệm dưới dạng giải tích của các phương trình cân bằng nhiệt của vệ tinh bằng phương pháp tuyến tính hóa tương đương theo tiêu chuẩn đối ngẫu.
- Nghiên cứu và phân tích được một số ứng xử định tính của nhiệt độ vệ tinh trong các mô hình nhiệt.
Đối tượng và phạm vi nghiên cứu của luận án tập trung vào:
- Đối tượng nghiên cứu: Đáp ứng nhiệt và trường nhiệt độ của vệ tinh cỡ nhỏ (tập trung khảo sát nhóm vệ tinh có khối lượng dưới 100 kg) chuyển động trên quỹ đạo thấp quanh Trái đất.
- Phạm vi không gian và mô hình: Quỹ đạo thấp quanh Trái đất (LEO - Low Earth Orbit) ở độ cao từ 300 km đến 1000 km, độ nghiêng quỹ đạo tròn từ $97^\circ$ đến $99^\circ$ (quỹ đạo đồng bộ mặt trời); mô hình nhiệt phân tích giới hạn ở các cấu trúc một nút, hai nút, sáu nút (thân vệ tinh hình hộp chữ nhật) và tám nút (thân hình hộp chữ nhật gắn cánh pin mặt trời).
- Phạm vi thời gian: Khảo sát diễn tiến nhiệt độ theo thời gian chuyển động trên quỹ đạo (chu kỳ $P_{orb}$) ở chế độ dừng, trạng thái bình ổn nhiệt và dao động nhiệt tuần hoàn theo chu kỳ quỹ đạo.
Tổng quan tài liệu và vị trí của luận án
Luận án đã tổng hợp và phân loại các hướng nghiên cứu chính trong lĩnh vực phân tích nhiệt vệ tinh trên thế giới:
- Hướng nghiên cứu bài toán mạng nhiệt và phương pháp số: Phương pháp tham số phân bổ (lưới nút nhiệt) được sử dụng để xây dựng các phần mềm chuyên dụng như ESATAN-TMS, SINDA, THERMICA. Krishnaprakas (1995) đã khảo sát các thuật toán số và chỉ ra phương pháp Crank-Nicolson cho nghiệm ổn định không điều kiện đối với hệ mạng nhiệt lớn. Milman và Petrick nghiên cứu định tính về tính duy nhất nghiệm và sự hội tụ của phương trình cân bằng nhiệt chứa số hạng bức xạ bậc bốn. Papalexandris nghiên cứu thuật toán điều khiển phản hồi hệ nhiệt. Appel và đồng nghiệp xây dựng thuật toán xác định ma trận dẫn nhiệt. Liu và cộng sự (2010s) kết hợp phương pháp hệ số biến thiên với phương pháp chiếu tia Monte-Carlo để hiệu chỉnh hệ số hấp thụ mặt trời và tính toán các hệ số truyền nhiệt bức xạ.
- Hướng nghiên cứu rút gọn mô hình (Model Reduction): Nhằm giảm chi phí tính toán cho các mô hình nhiệt bậc cao, Deiml và đồng nghiệp phát triển thuật toán rút gọn mô hình và tối ưu hóa cấu hình thiết kế; Rico và cộng sự đề xuất kỹ thuật xử lý ma trận thưa và đối xứng để giảm số phần tử mạng nhiệt; Gorlani và Rossi áp dụng phương pháp rút gọn lặp tự động cho thiết bị IRES N2 của vệ tinh GALILEO AVIONICA Firenze; Jacques và đồng nghiệp phát triển thuật toán kết đám phần tử hữu hạn cho vệ tinh Meteosat; Frey và đồng nghiệp xây dựng quy trình tự động hóa mô hình nhiệt hai giai đoạn.
- Hướng nghiên cứu tương tác nhiệt - kết cấu (Cơ - Nhiệt): Li và đồng nghiệp phân tích nhiệt của cánh vệ tinh có cấu trúc lõi tổ ong; Li và Yan nghiên cứu dao động nhiệt kết hợp của cánh vệ tinh LEO; Azadi và cộng sự khảo sát biến dạng nhiệt và độ võng của cánh vệ tinh dạng tấm ngàm chịu tải quỹ đạo tròn; Liu và cộng sự sử dụng thuật toán hiện phân tích tương tác kết cấu - nhiệt bức xạ.
- Hướng nghiên cứu yếu tố ngẫu nhiên và bất định: Akita và đồng nghiệp ứng dụng bộ lọc Kalman và phương pháp ước lượng thích nghi để xác định trở kháng nhiệt tiếp xúc và xử lý tham số bất định phi tuyến trong phương trình cân bằng nhiệt, kế thừa các hướng tiếp cận của Stumpel và Chalmers, Annan và cộng sự, Ueno và cộng sự.
- Hướng nghiên cứu giải tích cho bài toán nhiệt: Oshima và đồng nghiệp tìm nghiệm giải tích cho mô hình một nút và hai nút dưới tải hằng số; Arduini và cộng sự dùng kỹ thuật tuyến tính hóa cho bài toán ngược điều khiển nhiệt; Gadalla và Wahba khảo sát trường nhiệt giải tích và sai phân hữu hạn cho vệ tinh hình trụ quay; Grande và đồng nghiệp (1998) dùng khai triển Fourier và tuyến tính hóa phân tích dao động nhiệt mô hình hai nút; Gaite và cộng sự (2000s) phân tích vòng giới hạn và phương pháp nhiễu cho mô hình một nút, hai nút và mở rộng cho nhiều nút.
- Phương pháp tuyến tính hóa tương đương: Được khởi xướng bởi Booton (1953), Kazakov (1954), Caughey (1959, 1963) và phát triển từ nền tảng của Krylov - Bogoliubov (1937), tổng kết bởi Roberts và Spanos (1990, 2003), Socha (2008). Gần đây, Nguyễn Đông Anh và đồng nghiệp (2009, 2011, 2013) đề xuất tiêu chuẩn đối ngẫu cho hệ động lực phi tuyến ngẫu nhiên và tiền định.
Khoảng trống nghiên cứu được luận án tập trung giải quyết là: Khắc phục hạn chế của các phương pháp tuyến tính hóa truyền thống và giả thiết dao động nhỏ của Grande trong phân tích nhiệt giải tích; đưa tiêu chuẩn đối ngẫu của phương pháp tuyến tính hóa tương đương vào việc giải phương trình vi phân cân bằng nhiệt phi tuyến bậc bốn của vệ tinh nhỏ, cung cấp công cụ giải tích xác thực, có thời gian tính toán nhanh và độ chính xác tiệm cận với nghiệm số Runge-Kutta bậc 4.
Cơ sở lý thuyết và phương pháp nghiên cứu
Luận án xây dựng mô hình dựa trên nguyên lý cân bằng nhiệt động lực học và phương pháp tham số phân bổ (mô hình mạng nhiệt). Vệ tinh được rời rạc hóa thành các nút nhiệt đẳng nhiệt. Mỗi nút $i$ có nhiệt dung $C_i = \rho_i V_i C_{p,i}$ (trong đó $\rho_i$ là mật độ khối lượng, $V_i$ là thể tích, $C_{p,i}$ là nhiệt dung riêng), nhiệt độ $T_i(t)$ và nhiệt hao tán nội bộ $Q_{dis,i}$.
Các hình thức truyền nhiệt trong mô hình gồm:
- Dẫn nhiệt giữa nút $i$ và nút $j$: $Q_{c,ij} = k_{ij}(T_j - T_i)$, với $k_{ij} = \lambda_{mat} A / L$ là hệ số dẫn nhiệt.
- Bức xạ nhiệt giữa nút $i$ và nút $j$: $Q_{r,ij} = r_{ij}(T_j^4 - T_i^4) = \sigma A_i F_{ij}(\varepsilon_i T_j^4 - \varepsilon_j T_i^4)$, với $\sigma$ là hằng số Stefan-Boltzmann, $F_{ij}$ là hệ số quan sát hình học, $\varepsilon$ là hệ số phát xạ bề mặt.
Phương trình cân bằng nhiệt tổng quát cho nút $i$ trong hệ $n$ nút:
$$C_i \frac{dT_i}{dt} = Q_{external,i} + Q_{dis,i} + \sum_{j=1}^n k_{ij}(T_j - T_i) + \sum_{j=0}^n r_{ij}(T_j^4 - T_i^4)$$
Trong đó, nút $0$ biểu diễn môi trường không gian vũ trụ; $Q_{external,i} = Q_{sol,i} + Q_{alb,i} + Q_{e,i}$ là tổng các tải nhiệt môi trường không gian tại quỹ đạo thấp:
- Bức xạ mặt trời trực tiếp: $Q_{sol} = G_s A_{sp} \alpha_s f_s(t)$, với hằng số mặt trời $G_s = 1366.1 \text{ W/m}^2$, $\alpha_s$ là hệ số hấp thụ, $f_s(t)$ là hàm sóng vuông thể hiện chu kỳ chiếu sáng ($\rho = P_{il}/P_{orb}$).
- Bức xạ albedo của Trái đất: $Q_{alb} = a_e G_s A_{sc} F_{se} \alpha_s f_a(t)$, với hệ số phản xạ Trái đất $a_e \approx 0.3$ ($0.24 \div 0.42$), $F_{se}$ là hệ số quan sát Trái đất, $f_a(t)$ là hàm điều hòa mô tả góc thiên đỉnh.
- Bức xạ hồng ngoại Trái đất: $Q_e = \varepsilon A_{sc} F_{se} \sigma T_e^4$, với nhiệt độ vật thể đen tương đương của Trái đất $T_e \approx 255 \text{ K}$ (tương ứng dòng nhiệt trung bình $230 \text{ W/m}^2$, biến thiên $150 \div 350 \text{ W/m}^2$).
Về phương pháp xử lý toán học:
- Phương pháp tuyến tính hóa tương đương theo tiêu chuẩn đối ngẫu: Phương trình phi tuyến không thứ nguyên $d\theta/d\tau + f(\theta) = \xi(\tau)$ (với $f(\theta) = \theta^4$) được thay thế bởi phương trình tuyến tính $d\theta/d\tau + a + b\theta = \xi(\tau)$. Các hệ số $a, b$ và tham số tỷ lệ $\lambda'$ được xác định từ phiếm hàm cực tiểu hóa sai số đối ngẫu:
$$J = (1 - \lambda)\langle [f(\theta) - a - b\theta]^2 \rangle + \lambda \langle [a + b\theta - \lambda' f(\theta)]^2 \rangle \to \min_{a, b, \lambda'}$$
với $\lambda = 1/2$ đại diện cho sự bình đẳng giữa bước thay thế thuận và thay thế đối ngẫu.
- Phương pháp Newton-Raphson: Sử dụng để giải hệ phương trình đại số phi tuyến xác định các hệ số tuyến tính hóa $a$ và $b$.
- Phương pháp Runge-Kutta bậc 4 (RK4): Giải số hệ phương trình vi phân cân bằng nhiệt phi tuyến, đóng vai trò nghiệm chuẩn để đánh giá sai số của phương pháp giải tích và tính toán trường nhiệt cho các mô hình nhiều nút.
Nội dung chính theo từng chương
Chương 1: Tổng quan về bài toán phân tích nhiệt vệ tinh
Chương 1 hệ thống hóa các kiến thức cơ sở về phân loại vệ tinh theo khối lượng (pico, nano, micro, mini, trung bình, lớn) và theo quỹ đạo (LEO, MEO, HEO). Tác giả trình bày chức năng của 7 phân hệ chính trên vệ tinh (điều khiển tư thế, lực đẩy, điều khiển nhiệt, thông tin liên lạc với các băng tần S, C, X, Ku, xử lý dữ liệu, năng lượng, kết cấu) và quy trình 3 nhiệm vụ điều khiển nhiệt (phân tích nhiệt, thiết kế nhiệt thụ động/chủ động, thử nghiệm).
Chương này trình bày chi tiết miền nhiệt độ hoạt động tiêu chuẩn của các thiết bị vệ tinh:
| Thiết bị / Thành phần |
Nhiệt độ nhỏ nhất $T_{min}$ ($^\circ\text{C}$) |
Nhiệt độ lớn nhất $T_{max}$ ($^\circ\text{C}$) |
| Bộ cảm biến hồng ngoại |
-223 |
-173 |
| Cánh vệ tinh (tấm pin mặt trời) |
-100 |
+120 |
| Vỏ ngoài vệ tinh |
-100 |
+100 |
| Ăng ten đĩa |
-65 |
+95 |
| Thiết bị điện tử |
-10 |
+50 |
| Ắc quy |
0 |
+20 |
| Bể chứa Hydrazine |
+10 |
+50 |
Chương 1 cũng thiết lập toán học cho ba loại nút nhiệt (nút khuếch tán $C_i \neq 0$, nút số học $C_i = 0$, nút biên $T = \text{const}$), công thức truyền nhiệt dẫn nhiệt, bức xạ nhiệt, tích phân hệ số góc quan sát $F_{ij}$, định nghĩa góc quỹ đạo $\beta$, quy trình phân tích nhiệt 5 bước và tổng quan các công trình nghiên cứu trong nước và quốc tế.
Chương 2: Phân tích đáp ứng nhiệt của vệ tinh nhỏ trên quỹ đạo thấp dựa trên mô hình nhiệt một nút
Chương 2 thiết lập mô hình nhiệt một nút đẳng nhiệt cho vệ tinh nhỏ chịu tác động của bức xạ mặt trời, albedo và bức xạ hồng ngoại. Phương trình vi phân cân bằng nhiệt được chuẩn hóa về dạng không thứ nguyên:
$$\frac{d\theta}{d\tau} + \theta^4 = \mu_1 f_s(\tau) + \mu_2 f_a(\tau) + \mu_3$$
Khai triển Fourier giữ lại xấp xỉ bậc nhất đưa vế phải về dạng kích động $P + H\cos\tau$. Tác giả áp dụng phương pháp tuyến tính hóa tương đương theo tiêu chuẩn đối ngẫu ($\lambda = 1/2$) để thay thế số hạng phi tuyến $\theta^4$ bằng $a + b\theta$. Hệ phương trình phi tuyến xác định $a, b$ được giải bằng thuật toán lặp Newton-Raphson (khảo sát kèm miền hút hội tụ).
Nghiệm giải tích xấp xỉ của đáp ứng nhiệt không thứ nguyên có dạng:
$$\theta(\tau) = R + A\cos\tau + B\sin\tau$$
Các kết quả tính toán giải tích theo tiêu chuẩn đối ngẫu (DC) được so sánh định lượng với phương pháp tuyến tính hóa thông thường (CL), phương pháp Grande (G) và nghiệm số Runge-Kutta bậc 4 (RK). Kết quả chỉ ra rằng phương pháp đối ngẫu cho độ chính xác cao về nhiệt độ trung bình và biên độ dao động nhiệt khi khảo sát sự thay đổi của nhiệt dung $C$, hệ số hấp thụ mặt trời $\alpha_s$, và hệ số albedo $a_e$.
Chương 3: Phân tích đáp ứng nhiệt của vệ tinh nhỏ trên quỹ đạo thấp dựa trên mô hình nhiệt hai nút
Chương 3 mở rộng phương pháp tuyến tính hóa tương đương theo tiêu chuẩn đối ngẫu cho mô hình hai nút nhiệt: nút 1 đại diện cho vỏ ngoài vệ tinh và nút 2 đại diện cho khối thiết bị bên trong, có tính đến tương tác dẫn nhiệt ($k_{12}$) và bức xạ nhiệt nội bộ ($r_{12}$).
Nội dung chương tập trung làm rõ các đặc trưng động lực nhiệt:
- Diễn tiến nhiệt độ và không gian pha 3D: Khảo sát quỹ đạo pha và sự hội tụ của nghiệm nhiệt độ từ các điều kiện đầu khác nhau về một vòng giới hạn ổn định (limit cycle), chứng minh tính không nhạy cảm với điều kiện đầu khi thời gian chuyển động đủ lớn.
- Phân tích sai số và thời gian tính toán: Đánh giá sai số tuyệt đối của nghiệm giải tích so với nghiệm Runge-Kutta 4. Kết quả khẳng định phương pháp đối ngẫu DC đạt sai số nhỏ hơn đáng kể so với phương pháp tuyến tính hóa thông thường CL và phương pháp Grande G, đồng thời rút ngắn thời gian tính toán qua số chu kỳ quỹ đạo so với giải số thuần túy.
- Khảo sát tham số: Phân tích quy luật biến thiên của nhiệt độ trung bình ($\bar{\theta}_1, \bar{\theta}_2$) và biên độ dao động nhiệt ($\tilde{\theta}_1, \tilde{\theta}_2$) theo sự thay đổi của nhiệt dung nút trong $C_2$.
Chương 4: Tính toán đáp ứng nhiệt cho vệ tinh nhỏ trên quỹ đạo thấp sử dụng mô hình nhiệt nhiều nút
Chương 4 áp dụng phương pháp số Runge-Kutta bậc 4 để tính toán chi tiết trường nhiệt độ cho ba cấu hình kết cấu vệ tinh trong các kịch bản quỹ đạo khắc nghiệt:
- Mô hình nhiệt hai nút cho cánh vệ tinh: Nghiên cứu cánh gồm mặt trước (nút 1) và mặt sau (nút 2) trong hai kịch bản quỹ đạo và tư thế:
- Kịch bản 1: Cánh hướng vuông góc với bức xạ mặt trời; xác định dòng nhiệt $q_{sol,1}, q_{alb,1}, q_{IR,1}, q_{alb,2}, q_{IR,2}$.
- Kịch bản 2: Vệ tinh thay đổi tư thế quay; xác định phân bố nhiệt độ trên mặt trước và mặt sau cánh.
- Mô hình nhiệt sáu nút cho vệ tinh hình hộp chữ nhật: Rời rạc hóa 6 mặt vỏ vệ tinh thành 6 nút nhiệt ($+X, -X, +Y, -Y, +Z, -Z$), thiết lập hệ phương trình vi phân tương tác dẫn nhiệt và bức xạ. Tính toán trường nhiệt độ trong hai kịch bản:
- Kịch bản lạnh (Cold Case - CC): Góc $\beta = 0^\circ$, tải nhiệt môi trường cực tiểu.
- Kịch bản nóng (Hot Case - HC): Tải nhiệt môi trường cực đại. Xác định nhiệt độ cực đại ($T_{max}$) và cực tiểu ($T_{min}$) của từng nút, đối chiếu với ngưỡng an toàn.
- Mô hình nhiệt tám nút cho vệ tinh hình hộp gắn cánh: Bổ sung 2 nút cho cánh năng lượng mặt trời (mặt trước nút 7, mặt sau nút 8) liên kết với 6 nút thân vệ tinh (nút 1 đến 6). Khảo sát ba kịch bản: Cold Case (CC), Hot Case cho thân (HC1) và Hot Case cho cánh (HC2). Khảo sát ảnh hưởng của hệ số hấp thụ mặt trời $\alpha_F$, hệ số phát xạ mặt trước $\varepsilon_F$ và mặt sau $\varepsilon_R$ lên diễn tiến nhiệt độ của nút cánh và nút thân. Kết quả mô phỏng số xác nhận tất cả các nút nhiệt đều nằm trong ngưỡng giới hạn cho phép của thiết bị.
Kết quả và những đóng góp mới
Luận án đã đạt được các kết quả khoa học và đóng góp mới cụ thể sau:
- Về mặt lý luận và phương pháp:
- Lần đầu tiên phát triển và mở rộng thành công phương pháp tuyến tính hóa tương đương theo tiêu chuẩn đối ngẫu cho bài toán phân tích nhiệt phi tuyến của vệ tinh nhỏ trên quỹ đạo thấp quanh Trái đất.
- Thiết lập quy trình giải tích xác định nghiệm xấp xỉ dạng giải tích cho phương trình cân bằng nhiệt chứa số hạng bức xạ phi tuyến bậc bốn ($T^4$) trong mô hình nhiệt một nút và hai nút mà không cần phụ thuộc vào giả thiết dao động nhiệt nhỏ như các nghiên cứu trước đây.
- Chứng minh được tính hiệu quả, độ chính xác cao và sự ổn định của tiêu chuẩn đối ngẫu so với phương pháp tuyến tính hóa thông thường và phương pháp xấp xỉ Grande khi so sánh với nghiệm số Runge-Kutta bậc 4.
- Về mặt thực tiễn và ứng dụng:
- Xây dựng hoàn chỉnh mô hình toán học và bộ thuật toán tính toán tải nhiệt môi trường không gian (mặt trời, albedo, hồng ngoại Trái đất) cho các cấu hình vệ tinh nhỏ: mô hình cánh hai nút, mô hình thân vệ tinh hình hộp sáu nút và mô hình vệ tinh có gắn cánh tám nút.
- Cung cấp công cụ tính toán nhanh phân bố nhiệt độ và biên độ dao động nhiệt của các bộ phận vệ tinh trong các kịch bản quỹ đạo nóng nhất (Hot Case) và lạnh nhất (Cold Case), phục vụ trực tiếp cho việc lựa chọn vật liệu, lớp phủ bề mặt và thiết kế nhiệt thụ động trong giai đoạn tiền khả thi và thiết kế sơ bộ vệ tinh.
Hạn chế và hướng nghiên cứu tiếp
Dựa trên phạm vi và nội dung nghiên cứu được trình bày trong văn bản, luận án có các hạn chế và định hướng phát triển tiếp theo như sau:
- Hạn chế: Nghiên cứu giải tích của phương pháp tuyến tính hóa đối ngẫu mới dừng lại ở mô hình nhiệt một nút và hai nút; đối với các mô hình nhiều nút (sáu nút, tám nút), việc tìm nghiệm giải tích gặp khó khăn do số lượng phương trình đại số phi tuyến lớn nên tác giả phải chuyển sang sử dụng phương pháp số Runge-Kutta; các thông số vật liệu và hình học được xem xét trên các dạng kết cấu hình học cơ bản (hình hộp, tấm phẳng) ở quỹ đạo tròn đồng bộ mặt trời, chưa xét đến các quỹ đạo elip dẹt hoặc biến dạng cơ học phức tạp.
- Hướng nghiên cứu tiếp: Mở rộng phương pháp tuyến tính hóa tương đương theo tiêu chuẩn đối ngẫu để tìm nghiệm giải tích cho hệ mạng nhiệt có số nút lớn hơn ($n \ge 3$); tích hợp tính toán cơ - nhiệt đồng thời giữa dao động kết cấu tấm pin mặt trời và phân bố nhiệt độ; đưa vào mô hình các tham số ngẫu nhiên của môi trường nhiệt vũ trụ và khảo sát đáp ứng nhiệt ngẫu nhiên của vệ tinh.
Giá trị tham khảo
Luận án có giá trị tham khảo chuyên môn đối với các nhóm đối tượng:
- Nghiên cứu sinh, học viên cao học và giảng viên ngành Cơ học, Kỹ thuật Hàng không Vũ trụ, Kỹ thuật Cơ điện tử: Cung cấp tài liệu tham khảo có hệ thống về phương pháp giải tích giải hệ động lực phi tuyến bậc bốn, phương pháp tuyến tính hóa đối ngẫu, và phương pháp luận thiết lập mô hình nhiệt tham số phân bổ cho vật thể bay trong không gian.
- Kỹ sư thiết kế và chuyên gia công nghệ vũ trụ tại các viện nghiên cứu, trung tâm vũ trụ: Phần thuật toán tính toán tải nhiệt quỹ đạo LEO (Chương 1), phương trình cân bằng nhiệt nhiều nút và dữ liệu kịch bản mô phỏng Cold Case/Hot Case cho vệ tinh hình hộp và cánh pin (Chương 4) là tài liệu kỹ thuật hữu ích để ứng dụng trong thiết kế sơ bộ hệ thống điều khiển nhiệt vệ tinh nhỏ.
Câu hỏi thường gặp
1. Tại sao bài toán phân tích nhiệt vệ tinh lại xuất hiện tính phi tuyến bậc bốn và gây khó khăn cho phương pháp giải tích?
Phương trình cân bằng nhiệt của vệ tinh chứa số hạng truyền nhiệt bằng bức xạ nhiệt giữa các nút với nhau và giữa nút với môi trường không gian bên ngoài. Theo định luật Stefan-Boltzmann, năng lượng bức xạ tỷ lệ thuận với lũy thừa bậc bốn của nhiệt độ tuyệt đối ($T^4$ hoặc $T_j^4 - T_i^4$). Số hạng phi tuyến bậc bốn này khiến việc tìm nghiệm giải tích trực tiếp trở nên phức tạp, đòi hỏi phải sử dụng các kỹ thuật xấp xỉ tuyến tính hóa.
2. Tiêu chuẩn đối ngẫu trong phương pháp tuyến tính hóa tương đương có điểm gì khác biệt so với tiêu chuẩn sai số bình phương trung bình thông thường?
Tiêu chuẩn thông thường chỉ thực hiện một bước thay thế hàm phi tuyến $f(\theta)$ bằng hàm tuyến tính $a + b\theta$ và cực tiểu hóa bình phương sai số. Tiêu chuẩn đối ngẫu kết hợp hai bước thay thế: thay thế hàm phi tuyến bằng hàm tuyến tính, sau đó thay thế hàm tuyến tính này bằng một hàm phi tuyến cùng lớp $\lambda' f(\theta)$. Phiếm hàm sai số kết hợp cả hai bước với trọng số $\lambda = 1/2$, giúp cân bằng vai trò của hai bước thay thế và giảm thiểu đáng kể sai số nghiệm đối với các hệ có tính phi tuyến mạnh.
3. Ba nguồn nhiệt chính từ môi trường vũ trụ tác động lên vệ tinh trên quỹ đạo thấp (LEO) là gì?
Ba nguồn nhiệt chính gồm: (1) Bức xạ mặt trời trực tiếp ($Q_{sol}$), với hằng số mặt trời $G_s = 1366.1 \text{ W/m}^2$; (2) Bức xạ albedo ($Q_{alb}$) do ánh sáng mặt trời phản chiếu từ bề mặt Trái đất và các tầng mây vào vệ tinh; (3) Bức xạ hồng ngoại Trái đất ($Q_e$), phát ra từ bề mặt Trái đất và khí quyển với nhiệt độ vật thể đen tương đương $T_e \approx 255 \text{ K}$.
4. Kịch bản Cold Case (CC) và Hot Case (HC) trong Chương 4 được thiết lập nhằm mục đích gì?
Kịch bản Cold Case (CC) và Hot Case (HC) là hai trường hợp biên nhiệt độ giới hạn trên quỹ đạo. Kịch bản CC ứng với tải nhiệt môi trường nhỏ nhất (thời gian trong bóng tối dài nhất, góc quỹ đạo $\beta = 0^\circ$) để kiểm tra xem nhiệt độ các nút có bị giảm xuống dưới ngưỡng $T_{min}$ cho phép hay không. Kịch bản HC ứng với tải nhiệt môi trường lớn nhất để kiểm tra xem nhiệt độ các nút có vượt quá ngưỡng $T_{max}$ cho phép hay không, từ đó bảo đảm an toàn nhiệt cho toàn bộ thiết bị.
Kết luận
Luận án tiến sĩ của tác giả Phạm Ngọc Chung đã giải quyết bài toán phân tích đáp ứng nhiệt của vệ tinh nhỏ trên quỹ đạo thấp quanh Trái đất thông qua việc áp dụng phương pháp tuyến tính hóa tương đương theo tiêu chuẩn đối ngẫu kết hợp với phương pháp số Runge-Kutta bậc 4. Công trình đã xây dựng thành công hệ phương trình cân bằng nhiệt cho các mô hình một nút, hai nút, sáu nút và tám nút, đồng thời thiết lập biểu thức nghiệm giải tích cho nhiệt độ trung bình và biên độ dao động nhiệt. Các kết quả tính toán mô phỏng trong các kịch bản Cold Case và Hot Case đã chứng minh tính tin cậy của phương pháp và đảm bảo nhiệt độ vận hành của các thiết bị nằm trong giới hạn thiết kế. Đây là công trình có đóng góp thiết thực cho cơ sở lý thuyết động lực học nhiệt phi tuyến và công nghệ thiết kế nhiệt vệ tinh.
