Как рассчитать планетарную равновесную температуру в бинарных системах

Я работаю над программой, которая генерирует базовые характеристики планет земной группы в бинарных звездных системах. Я не силен в такой математике, поэтому у меня проблемы с вычислением температуры планетарного равновесия.

Формула, на которую я сейчас ссылаюсь, была опубликована в теме почти 4 года назад: Как рассчитать ожидаемую температуру поверхности планеты .

Формула

Это кажется достаточно точным, поскольку (если мне не изменяет память) его использование дает хороший ответ для Земли и Марса. Однако это относится к системе с одним излучающим телом, а не к системе с двумя излучающими телами (двойные звезды).

Мой инстинкт состоит в том, чтобы относиться к двум звездам как к одной большой звезде и добавлять мощность ( ж / м 2 ) планета получает от различных звезд. Таким образом, общее количество ватт = ватт звезда 1 + ватт звезда 2. Мне это кажется разумным, по крайней мере, для целей расчета планетарных температур.

Я уже проделал работу по расчету ватт, получаемых планетами из различных источников. Я использовал это уравнение, подставляя значения для обеих звезд.

Ф ( л / С 2 )

  • L = светимость (в долях солнечной), звезды 1 и 2 (соответственно): 0.0159&0.7758
  • S = большая полуось (в а.е.), звезды 1 и 2 (соответственно): 0.1228&3.6000
  • F = Средний солнечный поток Земли ( ж / м 2 ) 1362

Таким образом, приведенный выше расчет дает средний солнечный поток для примера планеты 1344 ж / м 2 от звезды 1 и 92 ж / м 2 от Star 2. Сложив их вместе, я получаю среднюю глобальную мощность 1436 ж / м 2 .

Однако я не уверен, как изменить исходное уравнение, чтобы использовать приведенный выше расчет солнечного потока (в идеале без ссылки на исходные переменные светимости и большой полуоси). Как правильно изменить формулу?

Любая помощь будет очень высоко ценится!

Ответы (1)

Мощность, излучаемая планетой, если предположить, что она имеет одинаковую температуру по всей своей поверхности:

п р а д "=" ϵ о Т е д 4 4 π р п л 2

Где ϵ это коэффициент излучения (чтобы соответствовать формуле в вопросе, установите это значение равным 1), о постоянная Стефана -Больцмана , Т е д - равновесная температура, а р п л это радиус планеты.

Мощность, поглощаемая планетой от звезд, если предположить, что она находится достаточно далеко от каждой звезды, что ее освещение незначительно за пределами 90 ° от подзвездной точки:

п а б с "=" я ф я ( 1 А я ) π р п л 2

Где ф я это поток от я ая звезда и А я - планетарное альбедо по отношению к я й звезды (может быть иначе, потому что отражательная способность будет зависеть от длины волны, а температура звезды влияет на то, какие длины волн излучаются).

Итак, приравняйте излучаемую и поглощаемую мощность и переставьте:

Т е д "=" [ я ф я ( 1 А я ) 4 ϵ о ] 1 / 4

Ограничение здесь состоит в том, что потоки будут меняться во времени из-за изменения расстояния планеты от звезд (то же самое относится к планете на эксцентрической орбите вокруг одной звезды), и нетривиально выяснить, как усредните это. Настоящей планете понадобится время, чтобы нагреться и остыть, и за это время изменятся расстояния и потоки.

Так что относитесь к числам, которые вы получаете из этого расчета «равновесной температуры», с щепоткой соли.


В соответствии с просьбой, вот несколько примеров расчетов. Во-первых, проверка здравомыслия: Земля. Принимая 1361 Вт/м 2 за значение солнечной постоянной и используя альбедо 0,3 и коэффициент излучения 1, температура получается как

Т е д "=" [ 1361   Вт м 2 × ( 1 0,3 ) 4 × 1 × ( 5.6704 × 10 8   Вт м 2 К 4 ) ] 1 / 4 255   К

Немного прохладно, но не так уж и плохо: парниковый эффект делает реальную Землю теплее, что можно учесть с помощью коэффициента излучения. ϵ < 1 .

Теперь для примера в вашем вопросе с потоками 1344 Вт/м 2 и 92 Вт/м 2 . Давайте также предположим, что в то время как первая звезда похожа на Солнце (поэтому я оставлю значение 0,3 для альбедо), вторая звезда холоднее Солнца. Эта звезда излучает больше света в красном цвете, где планета имеет меньшую отражательную способность. Я буду использовать альбедо 0,25 для этой звезды.

Расчет выполняется следующим образом (отбрасывая единицы измерения на шаге расчета для экономии места):

Т е д "=" [ 1344 × ( 1 0,3 ) + 92 × ( 1 0,25 ) 4 × 1 × ( 5.6704 × 10 8 ) ] 1 / 4 258   К

Надеюсь, это поможет!

Спасибо за этот ответ. У меня возникли некоторые проблемы с пониманием того, как применить уравнение. Не могли бы вы привести пример, подставив числа для двух планет? Если это слишком много, не беспокойтесь!
@n_bandit - конечно, не беспокойся.
Спасибо! Это отличный ответ. :) Один оставшийся вопрос касается степени 0,25. Изменится ли этот показатель при любых мыслимых условиях? Например, имеет ли она какое-либо отношение к скорости вращения планеты?
@n_bandit - степень 1/4 обусловлена ​​​​законом Стефана-Больцмана для излучения черного тела. Это включено в первое уравнение для мощности излучения планеты, которое содержит член σT⁴. Изменение этого показателя означало бы, что имеет место совершенно другой физический процесс (а также потребовало бы некоторых изменений в задействованных физических константах, чтобы единицы работали правильно).