Модель парникового эффекта

Анализируется состояние потепления Земли, предлагается механизм неустойчивости системы “Земля – атмосфера” к воздействию выбросов углекислого газа. Предлагаются модели, на основании которых находится область устойчивости. Анализируются возможные направления борьбы с потеплением Земли. Обосновывается...

Повний опис

Збережено в:
Бібліографічні деталі
Опубліковано в: :Промышленная теплотехника
Дата:2007
Автори: Присняков, В.Ф., Присняков, К.В.
Формат: Стаття
Мова:Російська
Опубліковано: Інститут технічної теплофізики НАН України 2007
Теми:
Онлайн доступ:https://nasplib.isofts.kiev.ua/handle/123456789/61306
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Модель парникового эффекта / В.Ф. Присняков, К.В. Присняков // Промышленная теплотехника. — 2007. — Т. 29, № 6. — С. 72-82. — Бібліогр.: 9 назв. — рос.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
_version_ 1860065700332175360
author Присняков, В.Ф.
Присняков, К.В.
author_facet Присняков, В.Ф.
Присняков, К.В.
citation_txt Модель парникового эффекта / В.Ф. Присняков, К.В. Присняков // Промышленная теплотехника. — 2007. — Т. 29, № 6. — С. 72-82. — Бібліогр.: 9 назв. — рос.
collection DSpace DC
container_title Промышленная теплотехника
description Анализируется состояние потепления Земли, предлагается механизм неустойчивости системы “Земля – атмосфера” к воздействию выбросов углекислого газа. Предлагаются модели, на основании которых находится область устойчивости. Анализируются возможные направления борьбы с потеплением Земли. Обосновывается реальность охлаждения Земли космическими средствами. Аналізується стан потепління Землі, пропонується механізм нестійкості системи «Земля – атмосфера» до впливу викидів вуглекислого газу. Пропонуються найпростіші моделі, на підставі яких в площині робочих параметрів знаходиться область стійкості. Аналізуються можливі напрямки боротьби з потеплінням Землі. Обґрунтовується реальність охолодження Землі космічними засобами. We analyze the state of warning of the Earth and propose a mechanism of the instability of the “Earth-atmosphere” system to the action of carbon dioxide emission. We propose models that make it possible to find the domain of stability. Different directions of the struggle against warming of the Earth, are analyzed, and the reality of cooling of the Earth by space means is grounded.
first_indexed 2025-12-07T17:07:24Z
format Article
fulltext Система “Земля – атмосфера” как динамическая система Обычно рассматриваются модели парниково; го эффекта в общем квазистатическом подходе, который позволяет компьютерным моделирова; нием рассчитать тренды изменения Земного тем; пературного режима. Но система “Земля – атмо; сфера” – динамическая система, которая требует рассмотрения уравнений динамики с последую; 72 ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ Аналізується стан потепління Землі, пропонується механізм нестійкості сис) теми «Земля – атмосфера» до впливу викидів вуглекислого газу. Пропонують) ся найпростіші моделі, на підставі яких в площині робочих параметрів знаходить) ся область стійкості. Аналізуються мож) ливі напрямки боротьби з потеплінням Землі. Обґрунтовується реальність охо) лодження Землі космічними засобами. Анализируется состояние потепле) ния Земли, предлагается механизм не) устойчивости системы “Земля – атмо) сфера” к воздействию выбросов углекислого газа. Предлагаются моде) ли, на основании которых находится об) ласть устойчивости. Анализируются возможные направления борьбы с по) теплением Земли. Обосновывается ре) альность охлаждения Земли космичес) кими средствами. We analyze the state of warning of the Earth and propose a mechanism of the instability of the “Earth)atmosphere” sys) tem to the action of carbon dioxide emis) sion. We propose models that make it pos) sible to find the domain of stability. Different directions of the struggle against warming of the Earth, are analyzed, and the reality of cooling of the Earth by space means is grounded. ББК 65.8 ПРИСНЯКОВ В.Ф., ПРИСНЯКОВ К.В. Институт геотехнической механики НАН Украины МОДЕЛЬ ПАРНИКОВОГО ЭФФЕКТА А – антропогенное энергетическое воздействие на атмосферу; А = 27,8 ·109 кВт; E – темп поступления энергии из ядра Земли; f – константа скорости разветвления цепи; j – плотность тока пара с поверхности испарения, в воздухе; m, mСО2 , mw – масса атмосферы (равная 5,15 ·1018кг), углекислого газа, водяных паров; mngh – масса непарниковых газов в атмосфере; m′w 0 и m′w – темп испарения воды, соответствующий максимуму излучения и текущий (m′w = 1,65·1010 кг/сек); m′СО2 – секундный приход углекислого газа; pngh и pgh парциальные давления непарниковых и парниковых газов; pH2O – парциальное давление водяных паров; Q – приход лучистой энергии от Солнца (Q = 10,5·109 Дж/м2год); dQ/dτ – величина тепловых потоков от нагретой поверхности Земли, поглощаемых парниковыми газами, равная Smq; Qs – темп поглощения солнечной энергии системой “атмосфера; поверхность Земли”; S – мидель Земли; Sm – площадь поверхности длинноволнового излучения Земли; Sw – площадь водных бассейнов и увлажненных поверхностей; Tа – средняя по объему температура; Tav – средняя температура атмосферы; Тt – приповерхностная температура; V – объем атмосферы, V = 7,65·109 км3; β – коэффициент массоотдачи; εст – cтепень черноты поверхности Земли; (при оледенении учитывает сезонное изменение отражательной способности); ζ – коэффициент отражения солнечной радиации от внешней поверхности атмосферы; τw – запаздывание в поступлении водяных паров в атмосферу; τα – период изменения наклона оси Земли в прецессии, нутации или при изменении эксцентриситета. Индекс: ( )0 – фиксированные параметры в окрестности разложения функции в ряд Тейлора. щим исследованием ее на устойчивость. Поэтому мы рассмотрим усредненную для всей Земли мо; дель с учетом запаздывающих эффектов. Особен; ностью тепломассобменных процессов в атмосфе; ре Земли как в газовом объеме, ограниченного с одной стороны поверхностью, большая часть ко; торой находится на линии насыщения, есть то, что ее температура вследствие постоянного испа; рения воды зависит только от давления. Кроме того, интенсивность испарения носит периоди; ческий характер вследствие смены дня и ночи. На такую закономерность накладывается испа; рение от конвекции воздуха, обтекающего вод; ные или мокрые участки. В нестационарных ус; ловиях меняется во времени как давление в атмосфере, так и температура пове;рхности Tt, а значит и приход испарившейся воды в атмосферу m′w. Известно, что в атмосфере в нормальных условиях содержится mw0 = (1,3...1,5)·1016 кг водяных паров, за год испаряется m′w= 5,2·1017 кг/год воды, т.е. в среднем секундный приход водяных паров в атмо; сферу равняется m′w=1,65·1010 кг/с. Отсюда следует, что среднеинтегральное время пребывания в атмосфере определенного ко; личества воды от момента испарения до кон; денсации (с выделением скрытой теплоты ис; парения) или время чистого запаздывания равняется τпр= 8,5 ·105 с ≈ 240 ч = 10 суток. В действительности реализация такой идеаль; ной картины маловероятна. Кроме того, на про; цесс конденсации водяных паров и выпадение осадков накладываются такие физические про; цессы как прогрев атмосферы с постоянной вре; мени 2 месяца, перемешивание и образование собственно центров конденсации по достижении относительной влажностью критического значе; ния. Темп процесса конденсации определенного количества воды происходит значительно быст; рее, чем накопление водяных паров за счет испа; рения. Процесс выпадения дождя или снега рас; тянут во времени. Очевидно, что ключевым моментом этого процесса является испарение воды с поверхности морей и океанов, с влажных участков суши, подъем их в атмосферу, переме; шивание и конденсация водяных паров с после; дующим выпадением в виде дождя или снега и с одновременным выделением скрытой теплоты фазового перехода. Это выделение теплоты в свою очередь вызывает дополнительное испаре; ние воды. При описании нестационарных про; цессов с подобными явлениями реальная зависи; мость доли испарившейся воды в атмосфере в период ее накопления до выпадении осадков за; меняется скачкообразным процессом, в котором каждая порция поступивших в атмосферу паров не конденсируется сразу, не переходит в жидкую или твердую (снег) фазу в течение времени τк, а мгновенно превращается в конденсат через неко; торое время запаздывания τк (см. рис. 1). Дейст; вительные кривые конденсации ϕ(τк) (на рис. 1 кривые 1 и 2) характеризуются сложным и посте; пенным процессом конденсации, закономер; ность которого трудно описать аналитически с учетом всех деталей физических процессов, про; текающих в атмосфере. От момента испарения воды происходит сложный процесс турбулентно; го перемешивания и прогрева воздуха атмосферы. С помощью указанной приближенной замены неразрешимая задача описания процесса кон; денсации паров воды в атмосфере через количе; ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 73 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ Рис. 1. Схематическое представление кривой образования конденсата воды (кривой конденсации паров в атмосфере): 1 – плохо перемешанные водяные пары, 2 – хорошо перемешанные водяные пары, 3 – ступенчатая линия скачкообразного процесса образования дождя или снега как представление кривой конденсации паров воды чистым временным запаздыванием. ственные характеристики его промежуточных этапов намного упрощается, поскольку единст; венной величиной, которая необходима для опи; сания процесса во времени будет только время τк от момента испарения воды до мгновенного пре; вращения паров в жидкую или твердую фазу. Это время будем называть временем преобразования или временем запаздывания конденсации паров воды. Очевидно, что трудности описания кривой образования конденсата воды переносятся в за; дачу определения времени фк. Время прогрева атмосферы равно τh = 5,2·106 с. Оценка времени перемешивания в представлении времени пре; бывания паровой фазы в атмосфере Земли как емкости дает величину ~ 104 с. Поэтому, если рас; сматривать наиболее неблагоприятный случай, то мы приходим к τк = 6,1·106 с. Это оценка этой величины из общих соображений. Таким обра; зом, наличие времени преобразования означает, что процессы фазовых превращений с выделе; нием скрытой теплоты парообразования в атмо; сфере запаздывают по сравнению с процессами испарения. Так, при изменении темпа испаре; ния вследствие возросшего количества паров воды и соответственно выделение скрытой теп; лоты испарения во время конденсации паров принимает значение, соответствующее изме; ненному приходу паров в атмосферу не сразу, а по истечении времени τк. Поэтому секундный расход пара в конденсат при образовании туч и выпадении осадков в момент времени τ будет равен m′k (τ)= m′w(τ–τк). При этих предположениях рассмотрим меха; низм неустойчивости течения парникового эф; фекта в атмосфере к возмущениям, вызванным антропогенной деятельностью с выбросом угле; кислого газа. Мы будем рассматривать колебания параметров системы “атмосфера – поверхность Земли”, связанные с тепловыми процессами не; стационарного поглощения длинноволнового излучения Земли и с выделением скрытой тепло; ты фазового перехода при выпадении осадков на Землю. В этом случае изменение давления в ат; мосфере влияет на время запаздывания конден; сации (и может считаться величиной второго по; рядка малости). Рассматривая эту систему как ав; токолебательный контур, выделим в нем колеба; тельное звено, источник энергии и обратную связь, обеспечивающую подпитку энергией ко; лебательного звена. Как известно, обратная связь дает колебательной системе энергию в нужной фазе, обеспечивая компенсацию неизбежных по; терь энергии колебаний. Энергия раскачки сис; темы выделяется в процессе конденсации паров воды и за счет поглощения части длинноволно; вого излучения парами води, углекислым газом и другими многоатомными молекулами. Обратная связь, управляющая потоком энергии, осуществ; ляется через время запаздывания. Это имеет мес; то в силу зависимости количества поступающих в атмосферу водяных паров от температуры по; верхности Земли (в свою очередь определяемой указанным поглощением и конденсацией (кос; венно)), и с другой стороны, в силу уменьшения величины парциального давления водяных паров после выпадения осадков. Таким образом, нали; чие запаздывания конденсации паров воды и их выпадение на землю из атмосферы, приводит к появлению положительной обратной связи меж; ду изменением парциального давления паров во; ды и нагревом атмосферы. Если рассматривать всю атмосферу в сосредоточенном виде (т. е. име; ющую одинаковые параметры по всему объему), то процесс повышения ее температуры и сопут; ствующие как следствия природные катаклизмы имеют чисто нестационарный характер и должны описываться общим уравнениями динамики ат; мосферы в обычных производных. Такое уравне; ние дает возможность находить изменение ее температуры во времени в зависимости от антро; погенного загрязнения СО2, величина которого характеризуется массовыми секундными выбро; сами m′СО2 или парциальным давлением рСО2 рав; ным рСО2 = mСО2 RСО2 Tа/V. (1) В начальный момент масса углекислого газа в атмосфере равняется mСО2 = 2,6·1015 кг, (отметим различные ее значения согласно разным источ; никам: mСО2 = 2,6·1015...5·1013 кг и даже 2,4·1017 кг или 2,6·1013 кг – естественно, это не влияет на ме; ханизм парникового эффекта, от этой величины 74 ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ могут зависеть только расчетные границы устой; чивости). Еще большее, чем в атмосфере, количество СО2 (98%) содержится в водоемах, где масса это; го газа достигает 130·1015 кг и он служит главным источником (на экваторе) и резервуаром – по; глотителем атмосферного СО2. Массовый при; ход СО2 по разным источникам равняется m′СО2 = 3,9·1015...2,6·1013 кг/год. По данным [6] в 2000 году эта величина равна 6,5·1012 кг/год = = 2,06·105 кг/сек (отметим, что некоторые авто; ры ее принимают равной 2...2,3·1014 кг/год = =1,24·108 кг/сек или (6,4...7,4)·106 кг/с – мы приводим все эти данные, чтобы читатель обра; тил внимание на необходимость их обоснова; ния). В литературе также указывается, что в сред; нем в последние годы парциальное давление углекислого газа растет 2 ppm/год. В 2004 г. его изменение составляло 2,6 ppm/год. Не исключе; но, что эти изменения вызваны вариациями его выделения из океана. Аналогично парциальное давление водяных паров рН2О = mw RН2О Tа/Vа, (2) где в начальный момент mw0 = (1,3...1,5)·1016 кг. Уравнение динамики атмосферы будет опре; деляться не только поступлением углекислого газа, но и поступлением паров воды, приход которых превышает равновесное из;за дополнительного ис; парения вследствие возрастающего поглощения длинноволнового излучения ее парами и СО2. Запишем уравнение баланса массы атмосферы Земли dm/dτ= m′w(τ)+ m′СО2 (τ)– m′k(τ) = m′w(τ)+ + m′СО2 (τ) – m′w(τ –τк), (3) где масса атмосферы m = mngh + mw +mСО2 = 4,5 ·1018 кг); масса углекислого газа mСО2 принимается равной 2,6·1015 кг, масса паров воды в атмосфере при обычных условиях mw = (1,3...1,5)·1016 кг, во влаж; ном воздухе содержится водяного пара у земной поверхности от 3...4% в тропиках до 2·10–5% в Антарк; тиде; секундный приход водяных паров в атмосферу будем принимать равным m′w=1,65·1010 кг/с), се; кундный приход углекислого газа – m′СО2 =1,24 ·108 кг/с. Величина массовой скорости испарения воды с поверхности Земли (приход пара в атмосферу) равняется [3] m′w = Sw j = Sw j0(ps–pH2O)/p=χ(ps–pH2O). (4) Здесь χ = Sw j0/p = Swβ – коэффициент, связан; ный с молекулярным весом, удельным объемом, молекулярной теплотой испарения, давлением насыщенных паров воды и температурой, при которой происходит испарение, находится по некоторой реперной точке, например, при давлении р = 105 Па на стационарном режиме (для воды при 15 oС χ = ~ 9,17·106 м·с), плот; ность потока пара с поверхности испарения j = 1,3·10–4 кг/м2·с), в воздухе при 20 оС j ~ 0,0002 кг/м2·с [3], плотность потока пара с поверхности в случае свободного испарения (в вакуум), для воды j0 = 0,1013 кг/м2·с. Масса атмосферы в виде смеси различных га; зов связана с другими ее параметрами усреднен; ным уравнением состояния m = pV/Ra Tа, (5) где давление в атмосфере равняется p = pngh + pH2O + pСO2 = pngh+ pgh. (6) Парциальное давление парниковых газов pgh= pH2O + pСO2 . Так как парциальное давление непарниковых газов pngh = const, то производная dm/dτ находит; ся следующим образом dm/dτ = [V/ Ra Tа,] {dpH2O /dτ + dpСO2 /dτ} – – pV/ (Ra Tа,)2 d(Ra Tа )/ dτ (7) Производная d(Ra Tа)/dτ может быть преобра; зована следующим образом: d(Ra Tа)/ dτ =[∂(Ra Tа)/∂pgh] dpgh/dτ. Отсюда получаем с учетом (3) и (7) {V/RaTа–pV/(Ra Tа)2[∂(Ra Tа)/∂pgh]}dpgh/dτ= =χ(ps –pH2O) + m′СO2 (τ) – χ{ps(τ–τк) – pH2O(τ–τк))}). (8) ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 75 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ Рассмотрим модель первого приближения, когда за время запаздывания конденсации давление углекислого газа меняется значи; тельно слабее, т.е. при dpH2O / dτ >> dpСО2 /dτ и [∂(Ra Tа)/∂pgh]dpgh/dτ = [∂(Ra Tа)/∂pH2O] dpH2O /dτ. Тогда уравнение динамики атмосферы запишет; ся в таком виде [V/Ra Tа,] – pV/(Ra Tа)2[∂(Ra Tа)/∂pH2O] dpH2O / / dτ = χ(ps–pH2O) + m′СО2 (τ) – χ{ps – pH2O (τ–τк)}. (9) В качестве нормирующих величин для парци; ального давления паров воды и темпа испарения введем значение давления в атмосфере в некото; рый стационарный начальный период p и испа; рения m′w 0 = χ (p0 s – p0 H2O). Умножим все части (9) на дробь p/m′w 0: pV/(m′w 0 Ra Tа)]{1–p/(Ra Tа)[∂(Ra Tа)/∂pH2O]} dpH2O/dτ = (p/m′w 0)χ(ps–pH2O) + pm′СО2 (τ) /m′w 0 – – (p/m′w 0)χ{ps(τ–τк) – pH2O (τ–τк))}. (10) Введем обозначения: pH2O = pH2O/p, θ = pV/m′w Ra Tа – постоянная времени атмосферы по водяному па; ру (равная 5,5·108 с), K = 1 – p/(Ra Tа)[∂(Ra Tа)/∂pH2O] – коэффициент теплового усиления парникового эффекта по водяному пару. 0 = (p0 s – p0 H2O)/p – относительная разность между давлением на; сыщенных паров при температуре поверхнос; ти испарения p0 s и парциальным давлением паров воды pH2O, в начальный момент при 15 оС 0 = 0,017). После некоторых преобразований получаем θKd p H2O/dτ = ( –p H2O)/ 0+ m′СО2 /m′w0 – ( –p H2O(τ–τк) )/ 0. (11) В парниковом эффекте ключевым считается избыток углекислого газа (чтобы не загромож; дать выкладки, здесь мы не будем учитывать ме; тан, который также провоцирует этот эффект, ибо учесть его не представляет особых труднос; тей). В нашем анализе основной переменной мы выбрали относительное парциальное давление водяных паров p H2O. Введем в рассмотрение константу скорости разветвления цепи f = pH2O/рСО2 (по аналогии с теорией цепного воспламенения Н. Семенова) и выразим через нее отношение m′СО2 /m′w 0. Величи; на загрязнения атмосферы углекислым газом, выражаемая его массовым секундным приходом m′СО2 (τ) = aτ, определяется не физическими зако; нами, а человеческой деятельностью (мы рассмо; трим наипростейший случай линейной зависи; мости). Поэтому она представляет исходное известное эмпирическое значение в некоторый рассматриваемый момент времени. В отличие от паров воды, которые конденсируются и перио; дически выпадают в виде осадков на Землю, уг; лекислый газ накапливается в атмосфере и ито; говая величина его парциального давления определяется интегралом mСО2 = aτ dτ= aΔτ2/2= m′СО2 Δτ/2. Значение Δτ представляет эмпирическое время, на интервале которого рассматривается устойчивость атмосферы к воздействию угле; кислого газа. Это время примем равным одно; му году, т.е. Δτ = 3,15·107 с. Такое приближение нивелируется переходом к парциальному дав; лению углекислого газа, которое определяется именно суммарным, накопленным за все вре; мя газом. Поэтому найдем с учетом (1) отно; шение m′СО2 /m′w 0=2mСО2 /(Δτm′w 0) = A pH2O/f 0, (12) A= 2θ 0/0,68Δτ, здесь учтено отличие газовых постоянных углекислого газа и воздуха RСО2 /Rа = 0,68, а также, что для реперной темпе; ратуры 15 оС 0 = 0,016. Мы принимаем случай 0 >0, так как величина относительной разни; цы между давлением насыщенных паров и пар; циальным давлением паров воды по условиям за; дачи положительная. Для случая ледяного или снежного покрова поверхности требуется допол; нительный анализ. Здесь есть определенное при; ближение относительно величины массы угле; кислого газа, накопленной к рассматриваемому времени, и переменности количества испарен; ной воды в связи с ростом ее температуры, но это мало сказывается на конечном результате, учи; тывая малость А по сравнению с другими величи; нами. Тогда для 15 оС получаем А = 0,82. pΔ pΔ pΔ pΔ 0 Δτ ∫ pΔs p pΔs p pΔ pΔ 76 ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ С учетом (12) уравнение (11) перепишем сле; дующим образом: θK 0 dpH2O/dτ = ( ;pH2O) +ApH2O/f – [ (τ–τк) – pH2O(τ–τк)] (13) Область устойчивости системы “Земля – атмосфера” Запишем (13) в отклонениях после разложе; ния функции в ряд Тейлора в окрестности стаци; онарной точки с ограничением величинами пер; вого порядка малости. Тогда с учетом обозначения θK 0 = G получаем G dδpH2O/dτ = δ –δpH2O +АδpH2O/f – –[δ (τ–τк) – δpH2O(τ–τк)]. (14) Вариацию давления насыщения при темпера; туре поверхности испарения вследствие потепле; ния климата за счет парникового эффекта δ опре; делим, учитывая, что ps = ps{Tt {рСО2 (pH2O)}}, в виде: δ =(∂ps/∂Tt) (∂Tt/∂рСО2 ) exp (∂рСО2 /∂pH2O) δpH2O = ={(∂Tt/∂рСО2 )exp/(∂Tt /∂ps)}(∂рСО2 /∂pH2O) δpH2O. Производная (∂ps/∂Tt) находится при помощи уравнения Клапейрона;Клазиуса (∂ps/∂Tt) =L/TtΔv (Δv – разница между удельным объемом пара и жидкости). Производная (∂Tt/∂рСО2 )exp находится по экспериментальной зависимос; ти потепления климата в функции от увеличе; ния содержания углекислого газа в атмосфе; ре(см. Табл.). Производная (∂рСО2 /∂pH2O) может быть выра; жена через константу скорости разветвления це; пи (∂рСО2 /∂pH2O) = 1/f. Таким образом, получаем для вариации давления насыщения при темпера; туре поверхности испарения зависимость δ =[(L/TtΔv)(∂Tt/∂рСО2 )exp/f]δpH2O= =B/f δpH2O , (15) где B = (L/TtΔv)(∂Tt/∂рСО2 )exp= (∂Tt/∂рСО2 )exp/ /(TtΔv/L) – параметр тепловой чувствительности атмосферы к воздействию углекислого газа. Подставляем (15) в (14) и после преобразова; ний получаем Gf/(f–A–B)) dδpH2O/dτ+δpH2O – – (f–B)/(f–A–B) δpH2O(τ – τк)]=0. (16) Дальнейшее исследование на устойчивость этого уравнения хорошо известно ([4]). Условие устойчивости можно найти методом D;разбие; ния, отыскав решение уравнения (16) в виде сум; мы функций типа δpH2O(τ)=Сexp(sτ). Соответст; вующее характеристическое уравнение имеет вид Gf/(f–A–B)s+1–(f–B)/(f–A–B)exp(–sτк)=0. Определяющими параметрами устойчивости рассматриваемой системы являются величины G, B, τк и f, от величины которых зависят численные значения корней s характеристического уравне; ния. При s = iω получаем уравнение границы ус; тойчивости системы “атмосфера – поверхность Земли”: i [Gf/(f–А–B)] ω +1 –(f–B)/(f–B–А))cos(ωτк)+ +i(f–B)/(f–B–А)sin(ωτк)=0. sp sp sp s p s p pΔ s p s ppΔ ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 77 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ Та б л и ц а . /* влажность 70% которое преобразуется в формулу τк/θ=fK 0 {kπ – arctg[2A(f–B–A/2)]1/2/(f–A–B)}/ /(f–B–A/2)]1/2, (17) где k – бесконечная последовательность целых и положительных чисел (k = 1,2,3,…). Расчет по (17) в плоскости параметров τк/θ= τк(B, f) представлен на рис. 2. Уточним константу скорости разветвления це; пи f, используя следующий прием. Примем, что на границе перехода к неустойчивости атмосфе; ры к парниковому нагреву (ориентировочно, 1960 год) при внешнем возмущении в виде угле; кислого газа константа скорости разветвления цепи f =1 (что очевидно, так как предполагает от; сутствие дополнительного испарения воды вследствие отсутствия роста углекислого газа сверх стационарной величины). Величина B для этого года, имеет порядок 0,1. Расчет для этого случая определяет границу τк/θ= 0,0285, что дает время запаздывания τк=1,5·107 с. Хорошо известными методами нетрудно опре; делить, что область «справа;внизу» от границы есть область неустойчивости. Граница устойчиво; сти, рассчитанная при k = 2, лежит слева от кри; вой, соответствующей k = 1, т.е. случай k = 1 включает все другие. Соотношение (18) дает возможность анализи; ровать влияние приведенных переменных на ус; тойчивость атмосферы Земли к парниковому на; греву за счет увеличения поступления углекислого газа, а значит также и ее устойчивость к другим возмущениям в виде поступления метана. Видно, что в такой простейшей модели устойчивость оп; ределяется тремя параметрами: отношением вре; мен τк/θ, параметром тепловой чувствительности атмосферы к воздействию углекислого газа B и константой скорости разветвления цепи f. Кроме этого, граница устойчивости зависит от относи; тельной разности между давлением насыщенных паров при температуре поверхности испарения и парциальным давлением паров воды в атмосфере 0, от коэффициент теплового усиления по во; дяному пару K. Экспериментальные данные показывают ка; чественное согласование полученных результа; тов с наблюдениями. Так, потепление климата сопровождается увеличением частоты выпадения осадков, т. е. уменьшением времени запаздыва; ния конденсации паров воды τк, что согласуется с предсказаниями теоретической модели. Умень; шение параметра тепловой чувствительности ат; мосферы к воздействию углекислого газа приво; дит к увеличению области устойчивости системы “Земля; атмосфера” к выбросу углекислого газа. Из выражения для В видно, что чем меньше x – степень сухости пара (в термодинамической тер; минологии) или паросодержание атмосферы (в терминологии двухфазных жидкостей), тем боль; ше Δv = v′′–v′ = (v–v′)/x и тем больше область ус; тойчивости. И наоборот, повышение интеграль; ной влажости приводит к уменьшению запаса устойчивости атмосферы к парниковому эффекту. Полученные зависимости позволяют оцени; вать количественные методы управления парни; ковым эффектом (которые в общем;то очевид; ны) изменением внутренних параметров, в частности изменением скорости испарения с по; верхности воды или влажных участков суши. Уменьшение m′w приводит к росту постоянной времени атмосферы Земли θ, а этот рост увели; чивает область устойчивости земной атмосферы. Модель позволяет оценивать техническую эф; фективность способов борьбы с потеплением за счет уменьшения испаряемости воды. Указан; ный эффект усиливается падением парциального давления паров воды, которое приводит, с другой pΔ pΔ 78 ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ Рис. 2. Область устойчивости при f=2,3,4,5,10 (на рисунке кривые, идущие снизу вверх), θ= 5,5·108с, K=1, k=1, относительный перепад давления на поверхности испарения 0= 0,016). pΔ стороны, к увеличению разности и тем са; мым к увеличению области устойчивости. Воз; можно, больше приемлемым может стать путь увеличения времени запаздывания конденсации паров τк, для чего необходимо выявлять факторы, которые задерживают дожди или снегопады и приводят к большему времени запаздывания. Это первые очевидные качественные выводы, которые должны сопровождаться анализом реа; лизуемости и стоимости. Динамика системы “Земля – атмосфера” Полученные результаты дают возможность бо; лее широко подойти к исследованию динамики системы “Земля – атмосфера”, путем опреде; ления ее передаточных функций и амплитудно; фазовых и частотных, действительных и мнимых характеристик. Так, передаточная функция меж; ду возмущениями массы выбрасываемого угле; кислого газа и парциального давления водяного пара W(s)=δpH2O(s)/δmCO2 (s)= 0/[sθK 0–(B/f –1) {1 – exp (–sτк)}]. Отсюда, принимая s = ± iω, несложно найти выражения для АФЧХ, действительных и мни; мых характеристик системы. Представленное выражение динамической ха; рактеристики «Земля;атмосфера» получено в простейшем случае, когда учитывается баланс массы. Более общим и корректным будет случай с учетом баланса энергии в атмосфере, ограни; ченной поверхностью Земли с одной стороны и границей с открытым Космосом – с другой для случая нестационарных процессов в открытой системе. Изменение во времени теплового состо; яния атмосферы определяется балансовым дина; мичес;ким уравнением в виде уравнения первого за;кона термодинамики для открытой системы: d(em)/dτ = Rатм + dQ/dτ + H+–Hk +E, где количество теплоты, выделяющейся при кон; денсации паров воды при выпадении осадков (дождя или снега) с расходом m′k(τ)= m′w(τ–τк)= Sw j(τ–τк) = Sw j0(ps(τ –τк) – – pH2O(τ –τк))/p = χ (ps –pH2O), H+ = Sw j0(ps(τ–τк) – pH2O(τ–τк))/p)L. Величина поглощаемой Землей и ее атмосфе; рой энергии в довольно длительном времени мо; жет считаться постоянной величиной. Так как на поглощение коротковолнового излучения изме; нения в количестве углекислого газа и водяных паров влияют слабо, то в этом смысле величина темпа поглощения солнечной энергии системой “Земля;воздух” QS может считаться постоянной величиной. В долгосрочной перспективе поступ; ление солнечной энергии действительно нерав; номерное и поэтому необходимо учитывать пе; риодическое (в течение десятков тысяч лет) изменение положения Земли в космическом пространстве. Поэтому будем для долговремен; ной перспективы считать, что получаемая Зем; лей солнечная энергия QS= QS(α), где α – некото; рый угол наклона Земли к перпендикуляру к линии между Землей и Солнцем. Темп поглощения коротковолнового излуче; ния атмосферой Rатм также можно считать в пер; вом приближении величиной постоянной. Тогда на динамическом режиме при переменной вели; чине значения приповерхностной температуры имеем cv Тt (RCO2 /RH2O) (∂рH2O/∂рCO2 +1) dmCO2 / dτ + +mcv ξ dТt/dτ – {cv ξ pRCO2 /{(Ra)2]} [∂(Ra) /∂(рH2O + + рCO2 )] (∂рH2O/∂рCO2 + 1) dmCO2 / dτ = Rатм + A + + Sm σ0 [(1+εст)/2] εg(Tt 4 –Tg 4)+ Sw j0(ps(τ–τк) – – pH2O(τ–τк)/p)L +Е. Отклонения в тепловом балансе от точки ста; ционарного режима будут равны: dδQатм/dτ = δQs +δIs+ δА+ δЕ δQs = Qs0 (∂S/∂α) δα(τ–τα), δIs = Sm δq +δL+δP. Газы обладают возможностью испускать и по; глощать энергию, но для различных газов эта особенность различна. В нашем случае нагретая до температуры Tt поверхность излучает энер; гию, которая проходит через объем атмосферы с pΔpΔ pΔ ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 79 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ увеличенным содержанием СО2 и Н2О (и СН4) (слой с большими их парциальными давлениями рСО2 , рН2О и рСН4 ) и больше поглощается в срав; нении с установившимся балансом в атмосфере. Тепловой поток от поверхности Земли, который поглощается увеличением количеством парни; ковых газов, приближенно определяется хорошо известной формулой q= σ0 [(1+εст)/2] εg(Tt 4 –Tg 4), εg = εСО2 + εН2О – εСО2 εН2О – суммарная поглоща; тельная способность смеси водяного и углекис; лого газов (принимаемая равной приведенной степени их черноты), которая определяется по таблицам в зависимости от величины парциаль; ного давления и длины пробега луча l = 3,6V/S; иногда учитывают, что εСО2 ~(рСО2 l)1/3, εН2О~(pl)2/3. Необходимо иметь в виду, что большой вклад в потепление косвенно оказывает метан. Если в этом случае ограничиваться перекрытием длин волн излучения только водяных паров и углекис; лого газа, то суммарная поглощательная способ; ность смеси этих трех газов будет равна: εg= εСО2 + εН2О – εСО2 εН2О + εCН4 . Представляемая модель позволяет учитывать поглощение и аэрозолей и др. газов. Мы ограни; чимся аналитическим рассмотрением только смеси водяных паров и углекислого газа. Поэто; му запишем: δq= σ0/2 {εg(Tt 4 –Tg 4)}0 δεст + σ0 [(1+εст0)/2] (Tt 4 –Tg 4)0 δεg + 2 σ0 (1+εст0) εg0 Tt 3δTt . Входящие в это уравнение отклонения степе; ни черноты различных газов определяются фор; мулами: δεg=δεСО2 +δεН2О – εН2О δεСО2 – εСО2 δεН2О = = (1– εН2О)0 δεСО2 + (1– εСО2 )0 δεН2О ; δεСО2 =[∂εСО2 /∂рСО2 ]0 (RСО2 Tav /V)δmСО2 ; δεН2О = [∂εН2О/∂рН2О]0 δрН2О(τ –τw); δεg = δεСО2 +δεН2О – εН2О δ εСО2 – εСО2 [∂εСО2 /∂рСО2 ]0 (RСО2 Tav /V)δmСО2 = (1–εН2О)0 [∂εСО2 /∂рСО2 ]0 (RСО2 Tav /V)δmСО2 + [1–εСО2 ]0 [∂εН2О/δрН2О] δрН2О(τ –τw). После подстановки этих выражений имеем δq= σ0 /2 {εg(Tt 4 –Tg 4)}0 δεст + σ0 [(1 + εст0)/2] (Tt 4 –Tg 4)0 {(1–εН2О)0 [∂εСО2 /∂рСО2 ]0 (RСО2 Tav /V) δmСО2 + [1–εСО2 ]0 [∂εН2О/∂рН2О] δрН2О(τ –τw)}+2 σ0 [(1+εст0)] εg0 Tt 3δTt . Если принимать обычную формулу для опреде; ления турбулентной теплоотдачи от поверхности Земли к атмосфере P = α(Tt–Tg), то легко записать значение отклонения ее величины вследствие уве; личения температуры δP = αδTt. Темп поступле; ния в атмосферу скрытой теплоты при испарении с водной поверхности и с увлажненных участков суши определяется скоростью испарения по (4) [3] L′= Sw jL = Sw Lj0(ps –pН2О(τ –τw))/p. Отсюда получаем отклонение темпа поступле; ния скрытой теплоты при испарении δL′= (Sw Lj/p0 ){(∂ps/∂Tt)δTt – δpН2О(τ –τw)}, где (∂ps/∂Ts) = L/[Ts(v′′ – v′)]. Записанные формулы дают возможность за; писать уравнение динамики системы “Земля;ат; мосфера”, для определения изменения при;по; верхностной температуры Земли: mcv ξ dδТt/dτ – {2 Sm σ0 [(1+εст0)] εg0 Tt 3+(Sw L j0/p)(∂ps/∂Tt) +α}δTt = {– cv Тt (RСО2 /RН2О) (∂рН2О/∂рСО2 +1) + cv ξ pRСО2 / {(Ra)2]} [∂(Ra) /∂(рН2О + рСО2 )] (∂рН2О/∂рСО2 +1)}dδmСО2 /dτ + +{Sm σ0 [(1+εст0)/2] (Tt 4 –Tg 4)0 {(1–εН2О)0 [∂εСО2 /∂рСО2 ]0 (RСО2 Tav /V) + [1–εСО2 ]0 [∂εН2О/∂рН2О]} – –(Sw L j0/p) {∂рН2О/∂рСО2 } (RСО2 Tav /V)} δmСО2 + +Qs0 (∂S/∂α) δα(τ–τα) +δA+ σ0 /2 {εg(Tt 4–Tg 4)}0 δεст +δЕ. Таким образом, мы видим, что если в первой простейшей модели, полученной на основании ба; ланса массы, выходным сигналом есть парциаль; ное давление водяных паров, реагирующих на входной сигнал в виде выбросов углекислого газа, то в модели, основанной на более полном описа; нии явления парникового эффекта выходом есть непосредственно интересуемая величина припо; верхностной температуры Тt, которая зависит от возмущений выбросов углекислого газа δmСО2 и его темпа dδmСО2 /dτ, а также от возмущений поло; 80 ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ жения Земли δα(τ–τα), которое некоторые авторы считают причиной потепления, от выбросов теп; лоты за счет антропогенной деятельности человека δA, от изменения степени черноты поверхности Земли за счет таяния ледников, помутнения воды и т.п. причины, которые также могут влиять на парниковый эффект и могут стать механизмом борьбы с ним (при определении количественного воздействия), от поступления теплоты из ядра пла; неты. Отсюда найти передаточные функции и ди; намические характеристики по каждому воздейст; вию системы “Земля;атмосфера” особого труда не представляет. Космические системы против глобального потепления Так как энергопотребление в мире на душу на; селения растет, то самый очевидный способ борьбы с потеплением – ее охлаждение путем уменьшения тепловых антропологических вы; бросов за счет приостановления роста численно; сти населения, что естественно не находит под; держки в наиболее вероятных кандидатов на это – в слаборазвитых стран. Но с другой стороны, ко; личество получаемой человеком дополнительной энергии вследствие прироста населения Земли (член δA) можно компенсировать уменьшением теплоты, приходящей от Солнца. Французский климатолог Э. Бар среди возможных геоинже; нерных методов рукотворного изменения клима; та в случае резкого и неожиданного ухудшения климатической ситуации (а это, возможно, при; дется уже через пару десятилетий) выделяет наи; более перспективное, по его мнению, эффектив; ное и менее всего опасное средство охлаждения Земли – зеркало между Землей и Солнцем, уменьшающее освещение Земли. В развитие вы; сказанных Э. Баром идей мы проанализировали вариант размещения на некоторой орбите кос; мических экранов (“космического зонта”), со; ставленных из уже разработанных сейчас и прохо; дящих экспериментальную проверку “солнечных парусов”, которые рассматриваются как основ; ной элемент перспективных космических двига; телей. Такой парус представляет собой огромное полотно из сверхтонкого материала (к примеру, алюминированного пластика), которое отражает падающий на него свет и под давлением фотонов создает тягу. В настоящее время разработаны ва; рианты пленок, конструктивное оформление, системы развертывания на орбите. Имеются и очень легкие с малой тягой высокоэкономичные электроракетные двигатели, которые, с одной стороны, будут стабилизировать огромных раз; меров космический “зонт”, а с другой, поддер; живать его на орбите. Таким образом, для созда; ния космических экранов уже имеется большой научный и конструкторский задел. Оценка основных параметров космического эк; рана (или нескольких экранов) показала следую; щее. В качестве базового мы выбрали вариант эк; рана, способного уменьшить количество падающей солнечной энергии на величину, кото; рая соответствует увеличению тепловой нагрузки на биосферу вследствие ежегодного современного прироста населения в 78 млн человек. При потреб; лении 2500 Вт на душу населения это дает величи; ну в 195 млрд Вт/год, что требует экран площадью 212 км2 (диаметр 16 км). При поверхностной плот; ности пленки 1 г/м2 получаем массу экрана 212 т. Такой парус за счет солнечного давления создает тягу по направлению к Земле в 90 Н, что требует приблизительно 100 ЭРД с тягой около 10 Н для ее компенсации. При фокусировании и определен; ном размещении компенсаторных двигателей это количество может быть сведено до нескольких. Основной проблемой, возникающей вследствие масштабных факторов отличия космического экрана от солнечного паруса (размеры больше на несколько порядков) является прочность пленки экрана и ее масса. Очевидно, что должен быть найден оптимум этих взаимоисключающих друг друга характеристик. При этом существенное облегчение для создания эк; рана в сравнении с солнечным парусом есть возмож; ность его размещения вдали от Солнца, в условиях менее напряженного теплового режима. В настоящее время некоторые изготовленные солнечные паруса имеют поверхностную плот; ность менее 10 г/м2. Прочность солнечного пару; са может быть увеличена за счет включения в па; рус металлических ребер, что уменьшит повреждения от столкновений с микрометеори; тами. Сейчас можно говорить о возможности со; здании тонкого экрана, состоящего из несколь; ких слоев крепкого, жаропрочного композитного ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 81 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ волокна с отражающим покрытием, который бу; дет иметь поверхностную плотность менее 1 г/м2. В будущем следует ожидать создания солнечных парусов с поверхностной плотностью 0,1 г/м2. Уже реально рассматривается создание инфраст; руктуры производства пленок для солнечных па; русов непосредственно в космосе. Используя ва; куумное напыление, можно изготовлять в космосе большие металлические листы толщи; ной 20...30 нм. При этом возможно получить по; верхностную плотность экрана до 0,05 г/м2. В бу; дущем нанотехнология может позволить создавать перфорированные или ячеистые плен; ки. Если перфорации в ней будут существенно меньше длины волны падающего на него света, малая масса и высокая отражательная способ; ность могут, дополнив друг друга, значительно по; высить производительность солнечных экранов. Что касается электро;ракетных двигателей, то для них наступает время широкой эксплуатации. Они имеют значительно большую удельную тягу, дающую возможность работать длительное время с небольшим запасом рабочего тела при неболь; шой потребной тяге, характерной именно для ЭРД. Сейчас мы имеем наш разработанный 25 лет назад ЭРД типа ТХД и российский СПД;290 с характеристиками достаточно близкими к тре; буемым в нашей задаче. Теперь относительно стоимости запуска всей си; стемы на орбиту. В перспективе космические само; леты будут выводить на орбиту объекты за 1;2 тыс. долларов за кг. Конечно, эта стоимость определя; ется орбитой, но как раз ЭРД могут также выпол; нить роль дешевых маршевых двигателей, кото; рые будут переводить космический экран с ближайшей орбиты на необходимую. Поэтому по расчетам стоимость солнечного экрана, компен; сирующего увеличение производства тепловой энергии на Земле за счет роста народонаселения, будет около 3 млрд долл. Эта величины сущест; венно меньше потерь в стихийных бедствиях или затрат на покупку квот по углекислому газу. Бо; лее того, предложенная программа позволяет де; лить экраны на отдельные части разных разме; ров, которые запускают разные страны, в частности которые дают наибольший вклад в прирост население и потребляют повышенное количество энергии. Если выделить такие “кос; мические” страны и континенты, то можно раз; делить экран например на 10 частей – для каждо; го субъекта – участника программы. Тогда для от; дельной страны будет необходимо обеспечить вывод на орбиту около 20 т, затратив на это при; близительно 300 млн долл. Достоинством предла; гаемой системы экранирования Земли является, с одной стороны, малая стоимость экрана, двигате; лей и всей установки, а с другой возможность вы; вода ее из эксплуатации в любое время без по; следствий для Земли, чего нельзя сказать о других вариантах вмешательства в биосферу Земли, ко; торые могут нанести своей неконтролируемостью больший вред, нежели ожидаемую пользу. ЛИТЕРАТУРА 1. Prisniakov V.F., Prisniakov K.V. Space Systems Against Global Warming//58th Intern. Astr. Cong. Paper. IAC;07;E5.3. Hyderabad, India, – September 24;28, 2007. 2. Присняков В.Ф. Динамика жидкостных ра; кетных двигательных установок. – М.: Машино; строение, 1983. – 248 с. 3. Ирисов А.С. Испаряемость топлив для пор; шневых двигателей. ГНТИ нефтяной и горно; топливной литературы. – М.: – 1955. 4. L. Hansen, M. Sato, R. Ruedy et all. Climate Simulations for 1880;2003 with GISS modelе. Draft of 13 Oct. 2006. 37p. 5. L. Hansen, M. Sato, R. Ruedy et all. Dangerous human;made interference with climate: A GISS model steady. Draft of 13 Oct. 2006. 6. Ахмад К.К., Анисимов О., Найджел Арнелл и др. Изменение климата, 2001 г. Доклад МГЭИК. Женева, 2001. – 60 с. 7. Шаманов Н.П. Энергетика в 21 веке. Науч; ный отчет НИЧ Снт;П ГМТУ. 8. Володин Е.М., Дианский Н.А. Воспроизве; дение изменений климата в 19;22 ст. // ИАН, Физика атмосферы и океана. – №3. – 2006. 9. Терез Э.И. Устойчивое развитие и проблемы изменения климата Земли. http://www/poteplenie.ru/ doc/terez.htm. 10. Смирнов Б.М. Углекислый газ в атмосфере Земли // УФН. – Т. 126, вып. 3. – 1978. – С. 527–530. Получено 18.09.2007 г. 82 ISSN 0204�3602. Пром. теплотехника, 2007, т. 29, № 6 ЭКОЛОГИЯ ТЕПЛОЭНЕРГЕТИЧЕСКИХ ОБЪЕКТОВ
id nasplib_isofts_kiev_ua-123456789-61306
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
issn 0204-3602
language Russian
last_indexed 2025-12-07T17:07:24Z
publishDate 2007
publisher Інститут технічної теплофізики НАН України
record_format dspace
spelling Присняков, В.Ф.
Присняков, К.В.
2014-04-30T07:08:45Z
2014-04-30T07:08:45Z
2007
Модель парникового эффекта / В.Ф. Присняков, К.В. Присняков // Промышленная теплотехника. — 2007. — Т. 29, № 6. — С. 72-82. — Бібліогр.: 9 назв. — рос.
0204-3602
https://nasplib.isofts.kiev.ua/handle/123456789/61306
65.8
Анализируется состояние потепления Земли, предлагается механизм неустойчивости системы “Земля – атмосфера” к воздействию выбросов углекислого газа. Предлагаются модели, на основании которых находится область устойчивости. Анализируются возможные направления борьбы с потеплением Земли. Обосновывается реальность охлаждения Земли космическими средствами.
Аналізується стан потепління Землі, пропонується механізм нестійкості системи «Земля – атмосфера» до впливу викидів вуглекислого газу. Пропонуються найпростіші моделі, на підставі яких в площині робочих параметрів знаходиться область стійкості. Аналізуються можливі напрямки боротьби з потеплінням Землі. Обґрунтовується реальність охолодження Землі космічними засобами.
We analyze the state of warning of the Earth and propose a mechanism of the instability of the “Earth-atmosphere” system to the action of carbon dioxide emission. We propose models that make it possible to find the domain of stability. Different directions of the struggle against warming of the Earth, are analyzed, and the reality of cooling of the Earth by space means is grounded.
ru
Інститут технічної теплофізики НАН України
Промышленная теплотехника
Экология теплоэнергетических объектов
Модель парникового эффекта
A model of the hotbed effect
Article
published earlier
spellingShingle Модель парникового эффекта
Присняков, В.Ф.
Присняков, К.В.
Экология теплоэнергетических объектов
title Модель парникового эффекта
title_alt A model of the hotbed effect
title_full Модель парникового эффекта
title_fullStr Модель парникового эффекта
title_full_unstemmed Модель парникового эффекта
title_short Модель парникового эффекта
title_sort модель парникового эффекта
topic Экология теплоэнергетических объектов
topic_facet Экология теплоэнергетических объектов
url https://nasplib.isofts.kiev.ua/handle/123456789/61306
work_keys_str_mv AT prisnâkovvf modelʹparnikovogoéffekta
AT prisnâkovkv modelʹparnikovogoéffekta
AT prisnâkovvf amodelofthehotbedeffect
AT prisnâkovkv amodelofthehotbedeffect