Фізика теплозахисного щита при поверненні — абляція, плазма, виживання

1. Гіперзвуковий режим і числа Маха

Режими течії класифікують за числом Маха M = v/a, де a — місцева швидкість звуку. Для атмосфери Землі на висоті a ≈ 295 м/с на 30 км, ≈ 340 м/с на рівні моря.

Режим	Діапазон Маха	Типовий приклад
Дозвуковий	M < 0.8	Крейсерський політ лайнера 0,82
Трансзвуковий	0.8–1.2	Зона звукового бар'єра
Надзвуковий	1.2–5	«Конкорд», SR-71
Гіперзвуковий	5–25	«Спейс Шаттл», повернення капсули
Орбітальне повернення	M ≈ 25–28	7,9 км/с із НОО; 11 км/с від Місяця

Понад Мах 5 домінують кілька нових ефектів:

Ефекти реального газу — O₂ і N₂ дисоціюють, потім іонізуються; показник адіабати γ падає з 1,4 до 1,2
Високотемпературна хімія — утворення NO, рекомбінація, каталітичні реакції на стінці
В'язкісна взаємодія — товсті примежові шари взаємодіють з ударними хвилями
Радіаційне нагрівання — за v > 10 км/с гарячий газ випромінює як абсолютно чорне тіло, іноді перевищуючи конвективне нагрівання

Орбітальна швидкість МКС: 7,66 км/с = 27 576 км/год = Мах 23 на висоті 400 км. Кінетична енергія 1 кг на цій швидкості становить ½·(7660)² = 29,4 МДж — еквівалент 7 кг тротилу. Усе це має бути перетворене на тепло, звук і випромінювання під час повернення.

2. Головна ударна хвиля та температура в точці гальмування

На гіперзвукових швидкостях перед апаратом утворюється відірвана головна ударна хвиля. Через прямий стрибок ущільнення повна ентальпія (ентальпія гальмування) зберігається. Для ідеального газу:

h₀ = h + v²/2 T₀ = T\infty + v²/(2\cdotcₚ) (температура гальмування) За v = 7900 м/с, T\infty = 240 K (30 км): T₀ = 240 + 7900²/(2\cdot1005) \approx 240 + 31 100 \approx 31 340 K Реальний газ T_stag \approx 8 000-11 000 K (енергія йде на дисоціацію, а не на температуру)

Величезна «втрата» температури (31 000 K теоретично проти 11 000 K реально) — насправді добра новина: енергія, яка нагріла б газ, натомість витрачається на розрив молекулярних зв'язків — дисоціюючи O₂ за ~5 еВ/молекулу та N₂ за ~9,8 еВ/молекулу. Це називають термохімічним стоком енергії.

Холодний набіжний потікСтиснутий ударною хвилеюГальмуванняПонад температуру поверхні Сонця

Точка гальмування — точка на носі апарата, де швидкість потоку дорівнює нулю — отримує найбільший тепловий потік. Навіть частка відсотка цієї енергії, що досягає конструкції, була б катастрофічною без теплозахисного щита.

3. Формула нагрівання в точці гальмування Чапмена

Детра, Кемп і Ріддел (1957) та Чапмен (1958) вивели напівемпіричні формули для конвективного теплового потоку в точці гальмування. Наближення Чапмена (холодна стінка) в одиницях СІ:

q̇ = C \cdot (ρ/ρ_sl)^0.5 \cdot (v/v_c)^3 / R_n^0.5 де: q̇ = тепловий потік [W/cm²] ρ = місцева густина набіжного потоку [kg/m³] ρ_sl = густина на рівні моря = 1.225 kg/m³ v = швидкість апарата [m/s] v_c = колова орбітальна швидкість \approx 7905 m/s R_n = радіус носа [m] C \approx 18 470 W\cdots³/cm²\cdotm\cdotkg^0.5 (холодна стінка)

Ключові висновки з цієї формули:

Залежність v³ — подвоєння швидкості збільшує нагрівання у 8 разів. Повернення від Місяця (11 км/с) у 2,7 раза гарячіше за повернення з НОО (7,9 км/с)
Залежність ρ^0.5 — нагрівання сягає піку в середній атмосфері (~50–70 км), а не на межі входу (~120 км)
Залежність R_n^−0.5 — більший радіус носа різко зменшує піковий тепловий потік. Apollo: R_n = 4,69 м → q̇_max ≈ 480 W/cm²

Пікове нагрівання Apollo 11: ~480 W/cm² упродовж ~40 секунд. Це 4,8 МВт/м². Побутова електродуховка дає близько 2 кВт на 0,1 м² = 20 кВт/м² — у 240 разів менш інтенсивно.

Радіаційне нагрівання

На надорбітальних швидкостях (v > 10 км/с) гарячий газ ударного шару випромінює як абсолютно чорне тіло. Для швидкості повернення від Місяця ~11 км/с радіаційне нагрівання може перевищувати конвективне. Наближення Таубера–Саттона для радіаційного потоку:

q̇_rad \propto ρ^1.5 \cdot v^8.5 \cdot R_n^1.0 Зауважте: залежність від R_n додатна для випромінювання, але від'ємна для конвекції. Оптимальний радіус носа мінімізує суму q̇_conv + q̇_rad.

4. Парадокс затупленого тіла

Г. Джуліан Аллен у NACA (1951) довів, що затуплений ніс переживає повернення краще за загострений. Це було контрінтуїтивно — гостріші об'єкти створюють слабші ударні хвилі в надзвуковому польоті (менший опір). Але затуплені тіла мають три переваги для повернення:

Відірвана ударна хвиля

Сильний прямий стрибок ущільнення стоїть далеко перед затупленим носом. Більшість енергії осідає в ударній хвилі та випромінюється геть від апарата, а не в нього.

Великий R_n зменшує q̇

Формула Чапмена: q̇ ∝ 1/√R_n. Радіус Apollo 4,7 м дає у 3 рази менший піковий тепловий потік, ніж ніс радіусом 0,5 м.

Високий опір для гальмування

Більший опір = швидше гальмування = менше часу на високій швидкості. Це різко зменшує сумарне теплове навантаження (проінтегроване за часом).

Стійке балансування

Широкі форми капсул (Apollo, Orion, «Союз») аеродинамічно стійкі під кутами повернення, що зменшує потребу в активному керуванні.

Сумарне теплове навантаження проти пікового теплового потоку:
Крутий вхід (−15°) максимізує опір і швидко гальмує → високий пік q̇, але коротка тривалість → менше сумарне теплове навантаження.
Пологий вхід (−5°) має нижчий пік q̇, але довшу тривалість → таке саме або більше сумарне теплове навантаження. Конструктивний оптимум залежить від питомої теплоємності теплозахисту.

5. Абляційні теплозахисні щити — як вони працюють

Абляційний теплозахист навмисно жертвує собою, щоб захистити конструкцію. Процес абляції охоплює три одночасні механізми:

Фазовий перехід і піроліз

У міру проникнення тепла в поверхню аблятора органічні в'яжучі піролізують (термічно розкладаються) за 300–500°C. Ця ендотермічна реакція поглинає тепло: ~1–3 МДж/кг. Гази піролізу витікають назовні, забезпечуючи ефект транспіраційного охолодження (тонкий прохолодний шар між гарячою ударною хвилею та гарячою поверхнею).

Поверхневі реакції

Решта обвугленого шару (вуглецева матриця) реагує з газами ударного шару. За температур понад ~3000°C вуглець абелює через окиснення (C + O → CO, де є кисень) і сублімацію (C → C₂, C₃ тощо). Кожна реакція відносить від поверхні значну приховану енергію.

Випромінювання

Розжарена обвуглена поверхня випромінює енергію назад за законом T⁴ (член чорного тіла). За 3000 K випромінювальна здатність = σT⁴ = 5.67e-8 × 3000⁴ ≈ 4,6 МВт/м² — порівнянна з вхідним потоком під час пікового нагрівання.

Тепловий баланс на поверхні аблятора (ṁ = швидкість втрати маси абляції на одиницю площі): q̇_conv + q̇_rad_in = ṁ\cdotH_eff + ε\cdotσ\cdotT_w⁴ + q̇_cond H_eff = ефективна ентальпія абляції [J/kg] \approx H_pyrolysis + H_surface_rxn + H_transpiration (~10-30 МДж/кг) q̇_cond = тепло, що проводиться в конструкцію («витік», який ми маємо мінімізувати)

Ефективність теплозахисту PICA: марсіанський посадковий модуль Phoenix використовував PICA (фенольно-просочений вуглецевий аблятор), розроблений NASA Ames. Ефективна ентальпія абляції ≈ 40 000–50 000 кДж/кг. Густина 0,27 г/см³. Загальна втрачена маса під час входу, спуску й посадки на Марсі: ~30 кг зі щита діаметром 1,5 м.

6. Матеріали теплозахисту: плитки «Шаттла» проти PICA проти AVCOAT

Матеріал	Тип	T_max (°C)	Густина (g/cm³)	Використання
AVCOAT	Абляційний, стільниковий	~2800	0.51	Командний модуль Apollo, SLS Orion
PICA	Абляційний, низької густини	~1650	0.27	Stardust, Phoenix, Dragon
PICA-X	Абляційний (варіант SpaceX)	>1650	~0.25	SpaceX Dragon 2
HRSI (LI-900)	Багаторазові плитки RSI	1260	0.144	Нижня поверхня «Спейс Шаттла»
FRSI	Багаторазове, гнучке покривало	371	—	Верхня поверхня «Шаттла» (низький нагрів)
Армований вуглець-вуглець	Конструкційний, багаторазовий	1650+	1.5	Носовий обтічник, передні кромки «Шаттла»
Starship TUFI	Керамічні плитки (багаторазові)	~1650	—	SpaceX Starship

Архітектура теплозахисту «Спейс Шаттла»

«Спейс Шаттл» застосовував зонований підхід на своїй поверхні, кожна зона відповідала піковій температурі, спрогнозованій аеротермічним аналізом:

Носовий обтічник і передні кромки крил (~1650°C) — RCC (армований вуглець-вуглець)
Навітряна нижня поверхня (~650–1260°C) — чорні плитки HRSI LI-900 з кремнеземної піни (~20 000 плиток)
Боки фюзеляжу та підвітряна поверхня (~400°C) — білі плитки LRSI (низькотемпературний RSI)
Верхня частина крил і верхні поверхні (<371°C) — гнучкі покривала FRSI з номексу

Катастрофа «Колумбії» (STS-107, 2003): удар шматка піноізоляції під час запуску пробив отвір 15–25 см у передній кромці RCC лівого крила. Під час повернення 1 лютого 2003 року плазма за ~1500°C проникла в конструкцію крила. За 7 хвилин конструкційна температура крила перевищила 1600°C — понад точку плавлення RCC. Крило зруйнувалося на Маху 18,4, на висоті 63 км.

7. Коридор повернення та балістичний коефіцієнт

Коридор повернення

Траєкторія повернення має пройти крізь вузький коридор, визначений:

Верхньою межею (надто полого): кут > −3° → недостатнє гальмування → відскік від атмосфери назад на орбіту
Нижньою межею (надто круто): кут < −8° (для Apollo) → надто швидке гальмування → перевантаження > 10–15g → конструкційний або фізіологічний збій

Коридор Apollo: кут нахилу траєкторії −5,5° ± 1,5° на межі входу 120 км. Це вузьке вікно (загалом 3°) потребувало точної навігації, яка була б неможливою без бортового комп'ютера наведення.

Балістичний коефіцієнт

Балістичний коефіцієнт β визначає, наскільки глибоко в атмосферу проникає апарат до значного гальмування:

β = m / (C_D \cdot A) [kg/m²] Високий β \to щільний, малий \to швидко пробиває атмосферу \to крута крива гальмування \to високе пікове перевантаження Низький β \to великий, легкий \to гальмує високо в розрідженій атмосфері \to нижче пікове перевантаження, нижча Tmax Командний модуль Apollo: β \approx 390 kg/m² «Союз ТМА»: β \approx 300 kg/m² SpaceX Dragon: β \approx 450 kg/m² Mars Pathfinder: β \approx 63 kg/m² (розріджене марсіанське повітря)

Проінтегроване теплове навантаження

Піковий тепловий потік визначає вимогу до температури поверхні теплозахисту. Проінтегроване теплове навантаження (J/cm²) визначає, скільки матеріалу вам потрібно:

Q_total = \int q̇(t) dt [J/cm²] Потрібна товщина аблятора \approx Q_total / (ρ_ablator \cdot H_eff) Командний модуль Apollo: Q_total \approx 30 000 J/cm² \to товщина AVCOAT 5-7 см

8. Радіомовчання та плазмова оболонка

На піку нагрівання плазма ударного шару достатньо щільна, щоб послаблювати або повністю блокувати радіозв'язок. Це радіомовчання під час повернення.

Фізика плазмової оболонки

Температури ударного шару іонізують повітря: N₂ + енергія → N₂⁺ + e⁻; O₂ → O⁺ + e⁻. Густина електронів n_e навколо капсули Apollo сягала піку ~10¹⁶ м⁻³. Радіохвиля може поширюватися крізь плазму, лише якщо її частота перевищує плазмову частоту:

f_p = (1/2π) \cdot \sqrt(n_e \cdot e² / (ε₀ \cdot m_e)) n_e = 10¹⁶ м⁻³ \to f_p \approx 900 MHz (діапазон УВЧ блокується) n_e = 10¹⁸ м⁻³ \to f_p \approx 9 GHz (X-діапазон блокується — GPS, радар теж блокуються)

Тривалість і способи пом'якшення

Місія	Початок радіомовчання (висота)	Тривалість
Apollo	~90 км (Мах 18)	~4 хвилини
«Союз»	~85 км	~3 хвилини
«Спейс Шаттл»	~73 км (Мах 12)	~16 хвилин (довша траєкторія)
Starship	Не визначено (вхід «черевом» подовжує радіомовчання)	~20–25 хв (оцінка)

Стратегії пом'якшення охоплюють:

Електрофільні упорскування: розпилення речовини, що захоплює вільні електрони → локально зменшує n_e, щоб відновити зв'язок
Наддовгі хвилі (VLF): сигнали нижче 100 кГц проникають крізь плазмові оболонки — використовувалися експериментально для каналів ядерного командування
Ретрансляція через супутникову геометрію: передача перпендикулярно до траєкторії польоту, де плазма тонша
Магнітогідродинамічне вікно: прикладання локального магнітного поля зменшує густину електронів у зазорі антени

Експеримент SPRITE (1960-ті): NASA випробовувала упорскування йонів калію під час повернень Gemini, частково відновлюючи зв'язок у S-діапазоні. Сучасні аблятори (PICA) мають нижчий вміст лужних йонів, ніж старіші аблятори, що зменшує природне самопом'якшення — виклик для майбутніх пілотованих апаратів.