Методика оценки возможности применения солнечной энергии для обеспечения длительного барражирования высотного дирижабля над заданной географической точкой

Научный Вестник МГТУ ГА_Том 21, № 04, 2018

Ovil Aviation High Technologies Vol. 21, No. 04, 2018

УДК 629.733

DOI: 10.26467/2079-0619-2018-21-4-96-109

МЕТОДИКА ОЦЕНКИ ВОЗМОЖНОСТИ ПРИМЕНЕНИЯ СОЛНЕЧНОЙ ЭНЕРГИИ ДЛЯ ОБЕСПЕЧЕНИЯ ДЛИТЕЛЬНОГО БАРРАЖИРОВАНИЯ ВЫСОТНОГО ДИРИЖАБЛЯ НАД ЗАДАННОЙ ГЕОГРАФИЧЕСКОЙ ТОЧКОЙ

И.Е. КОВАЛЕВ1, В.И. МАВРИЦКИЙ1, А.В. РЕДЬКИН1, И.Н. КАЧАРАВА1

1 ФГУП «Центральный аэрогидродинамический институт имени профессора Н.Е. Жуковского», г. Жуковский, Россия

Решение задачи оптимизации параметров беспилотного высотного дирижабля и его энергосистемы с целью получения максимальной продолжительности барражирования над заданной географической точкой требует эффективного сравнительного анализа применимости солнечной инсоляции в качестве основного источника энергии для силовой установки при различных высотных, широтных, сезонных и ветровых условиях. Наличие большого количества входных параметров при проектировании такого аппарата приводит к выводу о необходимости применения универсальных удельных энергетических характеристик высотного дирижабля и его энергоустановки для упрощения анализа и оценки эффективности технических решений. Основным фактором, определяющим потребляемую мощность силовой установки, является скорость ветрового потока над заданной точкой барражирования. Известно также, что наиболее доступным источником для длительного энергообеспечения является солнечная энергия. В большинстве случаев анализ решений сводится к оценке соотношения воздействия на высотный дирижабль двух природных явлений, солнечной энергии и энергии ветрового потока, для чего и была создана упрощенная модель высотного дирижабля, которая также может быть применена и к самолету на солнечной энергии. В работе в качестве основного показателя выбрано отнесенное к единице массы летательного аппарата значение мощности располагаемой и потребной. Располагаемая мощность определяется для любой географической точки с параметрами широта, высота, сезон (день года). Потребная мощность определяется высотой полета, скоростью и направлением ветра. Сравнивая значения удельных характеристик в пространственно-временных координатах, можно провести эффективный анализ избытка или недостатка располагаемой солнечной энергии для заданной географической зоны барражирования, для различных высот и дней года. Применяя отнесенные к взлетной массе летательного аппарата показатели, мы имеем возможность также выделить и оценить необходимую размерность летательного аппарата для реализации длительного барражирования в заданной точке пространства.

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

ВВЕДЕНИЕ

Техническая задача создания беспилотной аэростатической платформы (высотного дирижабля) в настоящее время является актуальной как реальная низкостоимостная альтернатива геостационарным спутниковым платформам, а также беспилотным носителям специального оборудования самолетного типа. В нашей стране реализация данной идеи может быть особенно эффективной для решения государственных задач по охране границ, обеспечению связи в труднодоступных районах, мониторингу атмосферы, поверхности Земли и прилегающих морей. Воздушное пространство Северного Ледовитого океана может контролироваться в основном авиационными комплексами большого радиуса действия или кораблями, способными преодолевать ледовые поля. В такой ситуации создание платформы, способной нести полезную целевую нагрузку на средних и больших высотах, может стать весьма экономичным и эффективным вариантом решения данных задач [1].

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

Vol. 21, No. 04, 2018

Ovil Aviation High Technologies

тельного аппарата (ЛА) в целом для обеспечения максимальной продолжительности полета над заданной точкой или в заданном районе. В идеале - полет должен продолжаться несколько месяцев и даже лет, аналогично современным спутниковым системам. Наиболее реальным источником энергии, способным обеспечить такой продолжительный полет, является солнечная энергия. Этот источник является возобновляемым в течение всего года на поверхности Земли от Южного полярного круга до Северного. Именно на ее использовании построена энергетика большинства реализуемых проектов высотных платформ [2], включая и самолетные варианты.

Специфика такой системы заключается в зависимости от инсоляционных условий (интенсивности солнечного излучения), зависимых от времени года (дня года), широты района и высоты полета [1]. Критической границей для создания круглогодичной системы энергоснабжения аппарата является широта, близкая к широте полярного круга в период минимальной продолжительности дня (22 декабря). Солнечная энергия в дневное время преобразуется в электрическую посредством батареи солнечных элементов (БСЭ), расположенных на оболочке дирижабля либо на других специальных элементах конструкции, обеспечивает потребности силовой установки (маршевых электродвигателей), бортовых систем и целевой полезной нагрузки, а также запасается в системе аккумулирования для энергообеспечения аппарата в ночное время [2, 3].

МЕТОДЫ И МЕТОДОЛОГИЯ ИССЛЕДОВАНИЙ

Для оценки соотношения воздействия двух природных явлений, солнечной энергии и энергии ветрового потока, используем упрощенную модель высотного дирижабля (ВД), которая также может быть применена и к самолету на солнечной энергии. Представим себе оболочку дирижабля в виде единой элементарной плоской поверхности, равной ей по площади и создающей эквивалентное аэродинамическое сопротивление с учетом дополнительного сопротивления формы и других элементов конструкции (гондол, оперения и др.).

Одновременно на этой же поверхности могут быть размещены солнечные элементы, соединенные в батареи, которые преобразуют солнечную энергию в электрическую. По соображениям рационального распределения эффективной полезной площади батареи солнечных элементов, практически она не должна превышать 50 % от площади поверхности оболочки. Зона их размещения определяется суточной траекторией движения солнца относительно оболочки аппарата и диапазоном изменения курсового угла 0.

В случае реализации задачи зависания над заданной точкой направление полета ВД будет противоположным направлению ветрового потока, а скорость его полета будет совпадать по величине со скоростью ветра относительно Земли. Основные геометрические параметры, характеризующие положение оболочки (корпуса) ВД относительно Земли и Солнца при фиксированном в течение суток направлении ветра, показаны на рис. 1.

Для любой плоской поверхности при фиксированном в течение суток курсовом угле 0 мы можем определить значение поступающей суммарной суточной инсоляционной энергии Есут. Данная величина зависит только от широты, высоты и дня года [3]. Реальная поверхность БСЭ расположена на оболочке с двойной кривизной поверхности. В связи с этим определение значения поступающей Есут должно выполняться с учетом соотношения между площадью БСЭ £сэ, расположенной на криволинейной поверхности, и площадью проекции 5прсэ этой поверхности на эквивалентную поверхность.

Представление высотного дирижабля в виде эквивалентной поверхности близко к истинной физической природе процесса, так как основную долю аэродинамического сопротивления дирижабля составляет сопротивление трения, которое может быть рассчитано как сопро-

Ovil Aviation High Technologies

Vol. 21, No. 04, 2018

тивление плоской пластины заданной площади соответствующей площади оболочки £об. Число Рейнольдса Яв определяется в зависимости от геометрических параметров оболочки.

Одновременно данная эквивалентная поверхность, ориентированная относительно Солнца в соответствии с курсовым углом 0, определяемым по направлению ветрового потока, является носителем солнечных элементов, воспринимающих солнечную энергию, которая посредством энергоустановки преобразует ее в поступательную энергию движения дирижабля.

Рис. 1. Положение в пространстве оболочки ВД, эквивалентной поверхности и ее проекции на нормаль к направлению на Солнце в случае, когда угол атаки а Ф 0 и курсовой угол 9 Ф 0 Fig. 1. The position in the space of the envelope high altitude airship envelope, equivalent surface and its projection on the normal to the direction of the Sun in the case where the angle of attack а Ф 0 and the course angle 9 Ф 0

Соотношение мощности или энергии, потребляемой и поступающей на эквивалентную поверхность, является показателем баланса эффекта воздействия на дирижабль двух природных факторов (солнечной энергии и энергии ветрового потока), определяющих реализуемость длительного барражирования над заданной географической точкой (широта, долгота, высота) в конкретный день года. Дополнительным важнейшим параметром исследования будет являться собственно размерность эквивалентной поверхности или площадь оболочки £об, приводимая к прямой зависимости от взлетной массы аппарата давзл. Теперь рассмотрим подробнее алгоритм расчета относительных значений потребной и располагаемой мощности для обобщенной модели высотного дирижабля.

Потребная мощность, необходимая для преодоления аэродинамического сопротивления ВД, в значительной степени определяется сопротивлением трения поверхности материала, которое в свою очередь зависит от числа Рейнольдса Яв. Сопротивление давления, зависящее от формы корпуса, является постоянным компонентом, не зависящим от Яв. Для оценки

Vol. 21, No. 04, 2018

Civil Aviation High Technologies

суммарного коэффициента сопротивления корпуса (оболочки) Схо6 для угла атаки а = 0° существует эмпирическая формула, показавшая достоверные результаты при сравнительных расчетах для реальных дирижаблей, указанная в работах [4] и [5]:

_ *^об

гп

хоб и2/3 об

или примем коэффициент формы Кф = и -=-

Уоб* 1 Х

Схоб —Кф-Cf

1+15 G)1^ er

L

где Д Ь - соответственно диаметр и длина оболочки (корпуса) дирижабля, а отношение Я = — является удлинением оболочки [4, 5]; С^ - коэффициент трения плоской пластины площадью 1 м , который может быть определен по формуле в зависимости от числа Рейнольдса Яв [5, 6]:

Cf = 0,0307- ,

где число Рейнольдса для дирижабля

Re = EILJLI:

где ц - динамическая вязкость воздуха.

В скобках выражения (1) первый компонент эквивалентен силе сопротивления трения, второй характеризует увеличение трения из-за возрастания скорости воздушного потока при обтекании формы, а третий - сопротивление, возникающее при отрыве потока [4].

Используя уравнение для определения аэростатической подъемной силы Ф [7]

9 • /ГгеП ^об ^ВЫП

где Ф= д • твзл - в случае, если масса ВД уравновешена аэростатической подъемной силой; g - ускорение свободного падения; ^вып - коэффициент выполнения оболочки; /НеН - удельная подъемная сила гелия на высоте И;

выражая объем оболочки иоб через объем цилиндра, в который вписана форма оболочки, с учетом показателя полноты объема оболочки (корпуса) у, являющегося характеристикой конкретно выбранной формы [6], выразим длину Ь, как характерный размер для определения числа Яв, через иоб и твзл:

£ _ 3 |4-Я2-Уоб = 3 I 4-Я2-Швзл

ф-П ^ ^-П-Квып-Гпеп ' К '

Дополнительное сопротивление, создаваемое хвостовым оперением, и сопротивление других элементов аппарата (гондола, расчалки, мотогондолы), согласно [5, 6], в среднем сум-

Civil Aviation High Technologies

Vol. 21, No. 04, 2018

марно оценим в 37 % от сопротивления оболочки (корпуса) Хоб, учтем это, подставив коэффициент 1,37 в формулу (1). В результате подстановки (2), (3), (4) в (1) получим уравнение для полного значения аэродинамического сопротивления ВД Хвд:

Хвд = 0,042

- ¥ 13 / I-N

КФ'РЬ7' % • / Швзл ^ /1 • 3 Î402

2 \KBMnfheh) ' \ u J тр'П'

Мощность для силовой установки или энергосистемы (в случае использования солнечной энергии) #потрэс, определенная в [4, 5] при скорости vкр:

N

^кр • ^вд

потрэс

Vb ' ^ред • ^эп • ^эд

где - КПД воздушного винта, ^ред - КПД редуктора, ^эп - коэффициент, учитывающий потери в силовых кабелях, проводах и соединениях между БСЭ, аккумуляторной батареей (АКБ) и маршевыми ЭД, ^эд - КПД электродвигателя, включая преобразователь (контроллер). Использованные при расчетах значения коэффициентов приведены в табл. 1 аналогично [3].

Таблица 1 Table 1

Расчетные значения КПД для элементов электрической силовой установки Estimated efficiency values for the elements of an electric power plant

Коэффициент VB ^ред ^эп ^эд Vb '^ред '^эп '^эд

Значение 0,75 0,97 0,95 0,94 0,65

В результате из (5) и (6) потребная мощность #поТрЭс, отнесенная к взлетной массе ВД:

1 + +

N,

потрэс _

ТПвЧЛ

0,042 • Кф • ph7 • VkP • Швзл 21

¿'Vb • ^ред • ^эп • ^эд • (Хвып • fheh) 21

Располагаемая мощность системы энергообеспечения с использованием солнечной энергии зависит от внешнего фактора - поступающей на единицу поверхности оболочки ВД солнечной инсоляции, а также от внутреннего фактора - удельных весовых характеристик всей системы и КПД собственно солнечных элементов [3]. Принципиальная схема системы энергообеспечения выглядит так, как показано на рис. 2.

Рис. 2. Принципиальная схема, основные компоненты и параметры энергосистемы на солнечной энергии Fig. 2. Structure diagram, main components and parameters of solar power system

Vol. 21, No. 04, 2018

Civil Aviation High Technologies

Массы БСЭ и системы аккумулирования являются определяющими для всей системы. Алгоритм их расчета подробно описан в [3]. Задача расчета - определить соотношение между массами этих элементов так, чтобы количество вырабатываемой днем солнечными элементами энергии точно соответствовало сумме дневной потребляемой энергии и энергии, необходимой для заряда АКБ. Именно при таком условии будет обеспечена равномерная среднесуточная мощность #сс. Пример типового суточного цикла распределения энергии, построенный по результатам расчета поступающей на поверхность БСЭ инсоляции, показан на рис. 3.

Рис. 3. Суточный цикл работы энергосистемы на солнечных элементах Fig. 3. The daily cycle of the power system operation by solar cells

На основании геометрической интерпретации рис. 3 определяется баланс поступающей энергии Есутсэ, которая уходит на питание постоянных потребителей Еднпотр (под отрезком АВ) днем и на заряд системы аккумулирования (АКБ) Езаракк:

Есутсэ Еднпотр + Езаракк-

При этом в процессе заряда-разряда АКБ часть энергии переходит в тепловую энергию, являясь потерями, определяемыми КПД цикла аккумулирования щакк:

Еакк Езаракк • Щакк-

Соотношение между Еднпотр и Еакк определяется с помощью геометрической пропорции заштрихованных трапеций (рис. 3), сторонами которых являются продолжительности дня ночи ^оч с учетом времени перехода ^ер:

24 - t„„ + t,

дн ^ *-пер

Еднпотр ^дн ^пер

Научный Вестник МГТУ ГА_Том 21, № 04, 2018

Civil Aviation High Technologies Vol. 21, No. 04, 2018

Привязка к весовым параметрам для БСЭ определяется выражениями:

Есутсэ ^сэ • ./сутсэ • Щсэ

и

тбсэ ^сэ • 7сэ • ^ком,

где усэ - удельный вес квадратного метра солнечных элементов; кком - коэффициент, учитывающий увеличение веса за счет дополнительных элементов; тбсэ - масса БСЭ; усутсэ - удельная суммарная суточная инсоляция, приходящаяся на квадратный метр БСЭ. Для массы системы аккумулирования такк:

такк • куст:>

Уакк

где уакк - удельный вес системы аккумулирования по отношению к запасенной энергии, куст -коэффициент, учитывающий дополнительный вес установочных и соединительных элементов аккумуляторов.

Суммарная полезная солнечная энергия ХЕп0л и среднесуточная мощность системы Ысс связаны уравнениями:

^Епол Еакк + Еднпотр Nсс • 24.

Используя уравнения (8), (9), (10), (11), (12), (13), (14), получаем расчетное выражение для Бсэ и

S = _:_^__

сэ Усутсэ • ^сэ • куст , 5 ^ ^

7 1 \ Усэ ' кком

+ • Уакк

^77 _ _™-Ж ' (1 + ^цикл)_

¿^ол = Куст + Усэ • кком . i 1 + к Л ,

Уакк Усутсэ • ^сэ ^акк циклу

где тэс - масса энергосистемы на солнечной энергии:

тэс = тбсэ + такк.

Используя статистическую формулу для массы конструкции шкон дирижаблей мягкой конструкции согласно [5]:

т = 3 53 • т 0,809

кон '"взл

и принимая массу полезной нагрузки тпн = 0, получим выражение для массы энергосистемы через взлетную массу:

ШэС= твзл - 3,53 • твзл0,809

и, используя (16), получим для :

Том 21, № 04, 2018_Научный Вестник МГТУ ГА

Vol. 21, No. 04, 2018 Civil Aviation High Technologies

Wcc _ 1

твзл 24 KZ£L + Усэ • кком • f i + к \ '

Уакк Усутсэ ' ^сэ ^акк циклу

Полученная формула (18) позволяет выполнить расчет, но при наличии известного значения поступающей среднесуточной инсоляции усутсэ на единицу площади БСЭ (Втч/м2сут). Точный расчет усутсэ может быть выполнен для каждого конкретного случая (широта, высота, день года), а также для заданной геометрии оболочки дирижабля, расположения и площади БСЭ £сэ с помощью разработанной автором в САПР 3D-модели Земли и ВД (рис. 4), [3].

Однако в таком случае мы не получим уравнения для в общем виде. Эффективной

^взл

площадью, воспринимающей солнечную инсоляцию, будет являться проекция поверхности БСЭ на нормальную к направлению на солнце плоскость в данный момент времени. Мы же имеем фиксированные расчетные значения среднесуточной инсоляции только для плоских пластин, зафиксированных в течение суток под заданным углом относительно солнца.

Рис. 4. 3D-модель Земли и дирижабля для расчета поступающей солнечной энергии на поверхность БСЭ Fig. 4. 3D model of the Earth and airship to estimate the incoming solar energy on the surface of the solar cells battery

Поэтому «краеугольным» камнем данного исследования является задача о нахождении соотношения между площадью БСЭ £сэ и площадью ее проекции 5прсэ. При этом зависимость должна выражаться через фиксированные параметры и основной параметр давзл. Задача носит геометрический характер, основные параметры и условия для ее решения показаны на рис. 5.

Угол ф назовем углом раствора сектора, на котором расположены солнечные элементы. В случае если БСЭ расположены по всей длине оболочки X = 0, Х2 = Ь), то методом интегрирования, согласно [8], мы получим достаточно простые уравнения:

5сэ - ф

и

5прсэ sin<p

^ =

5об

Civil Aviation High Technologies

Vol. 21, No. 04, 2018

Рис. 5. Геометрические параметры для определения зависимости площади проекции и полной площади БСЭ Fig. 5. Geometric parameters for determining the dependence of the projection area and the total area of the solar cells

battery

Отметим, что случаи неполного сектора (X1Ф 0, X2 ф L), а также других более сложных форм были также рассмотрены авторами и решены с помощью определенных интегралов, но являются частными случаями, которые могут быть описаны в отдельной работе. Определение такой формы является отдельной оптимизационной задачей.

Выражая площадь поверхности оболочки £об через взлетную массу давзл, коэффициент формы поверхности Кф, как это сделано ранее для уравнений (1), (4), получим для ф

2/

т— 5сэ • П .^вып • fheh\ '3 ^ Кф К Шв,л /

Подставляя (20) в (18) и выражая Sсэ, получим

2/:

asín

п • ^пр (К,

Кф

вып • fheh\ /3\. £ф Швзл J ) П \Кв

г/ъ

:•fheh'

(22).

Уравнение (22) отражает взаимосвязь между площадью БСЭ и площадью проекции £прсэ в зависимости от основных летно-технических и геометрических характеристик дирижабля: взлетной массы, высоты полета, коэффициента формы, степени выполнения оболочки на расчетной высоте полета к. Однако данное уравнение не позволяет определить значения массовых характеристик всех элементов системы. Используя уравнение (17) для компонентов взлетной массы аппарата, а также уравнение соотношения усутш и усутсэ

„/сутсэ

_ ^прсэ

„/сутэщ

получим уравнение, в котором исключены «криволинейные» параметры у'сутсэ,

^прсэ • 7сутэп • ^сэ • куст I  I Л • ^прсэ ÍК,

I ÚiS í n

(чакк + ^цикл) • Уакк

Кф

— швзл + 3,53 • тв

2/

\ / 3 ^ Кф • Усэ • ^ком

0,809

Значения у^™ (Вт^ч/м ), н, ¿пер могут быть определены для плоской пластины при заданном направлении движения ВД, сориентированной перпендикулярно плоскости движения

Vol. 21, No. 04, 2018

Civil Aviation High Technologies

солнца относительно оболочки аппарата, с помощью разработанной автором в САПР 3D-модели Земли и ВД (рис. 4). Уравнение (24) может быть решено математическими либо численными методами в программе Excel с помощью функционала «подбор параметра», в результате чего получаем значение £прсэ, затем, используя (23), - значения усутсэ, 5сэ. Остальные параметры и, собственно, величину соотношения —— можно определить, используя (18).

^взл

Для проведения предварительных расчетных исследований вариантов ВД были использованы на основании данных [2, 5, 9] следующие значения коэффициентов и удельных характеристик системы, приведенные в табл. 2.

Таблица 2 Table 2

Расчетные значения КПД и удельные характеристики элементов системы энергоснабжения ВД на солнечной энергии Estimated values of efficiency and specific characteristics of the power supply systems elements in airship on solar energy

Коэффициент Уакк Усэ, ^акк Г)сэ кком куст

Единицы изм. Вт-ч/кг кг/м2 - -

Значение 256 0,8 0,8 0,22 1,15 1,15

Значения куст, уакк и ^акк приняты по наилучшим показателям для серийных образцов литий-ионных аккумуляторов корпорации Кокат1.

Общая последовательность нахождения основных параметров при решении задачи нахождения удельной располагаемой энергии может быть записана следующим образом:

Уакк:> Усэ, Цак^ кком:> куст:> ^-вып

кф , fheh, mвзл, ^иад ,/сутэп > $прсэ > $сэ > .

^взл

РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЙ

Расчетные исследования и анализ эффективности применения солнечной энергии для длительного энергоснабжения высотного дирижабля был выполнен для шести характерных точек территории РФ и близлежащих стран: 40° с. ш. 80° в. д.; 50° с. ш. 80° в. д.; 60° с. ш. 80° в. д.; 40° с. ш. 140° в. д.; 50° с. ш. 140° в. д.; 60° с. ш. 140° в. д. получены, ориентируясь на характерные сезонные состояния атмосферы, как описано в [10], для расчетов использованы ветровые данные (у^ и 0) согласно стандарту2 в зимнем (январь) и летнем (июль) сезонах, характеристики воздуха р\,, л для стандартной атмосферы30 на высотах h.

Был рассмотрен размерный ряд высотных дирижаблей с твзл = 3000, 10000, 30000 и 100000 кг, расчет выполнялся для трех значений высоты барражирования ^ = 15, 17 и 20 км.

Результаты исследований были получены в виде расчетных таблиц с параметрами и в виде диаграмм со значениями удельных располагаемой и потребной мощностей в Вт/кг. Пример результатов расчета в виде диаграмм для 60° с. ш. 80° в. д., высот полета ^ = 15, 17, 20 км, зимнего сезона (январь) показан на рис. 6, а, летнего сезона на рис. 6, б. Левые столбцы соответствуют удельной потребной мощности, правые - располагаемой.

1 Технические характеристики литий-ионных полимерных аккумуляторов [Электронный ресурс] // Корпорация Kokam ltd. Режим доступа: http://kokam.com/cell/

2 ГОСТ 24728-81. Ветер. Пространственное и временное распределение характеристик: государственный стандарт СССР. М.: Изд-во стандартов, 1981.

3 ГОСТ 4401-81. Атмосфера стандартная. Параметры: межгосударственный стандарт. М.: Изд-во стандартов, 1981.

Civil Aviation High Technologies

Vol. 21, No. 04, 2018

Рис. 6. Диаграммы удельной потребной и располагаемой мощностей энергоустановки ВД на 60° с. ш. 80° в. д.

для сезонов: а - летнего

С использованием данной методики может быть рассчитан любой массив данных для ВД с охватом всей территории Земного шара. Заданная сетка ветровых данных и их точность будет определяющим фактором при выполнении расчетных исследований по возможности реализации длительного барражирования беспилотного ВД.

Vol. 21, No. 04, 2018

Civil Aviation High Technologies

ЗАКЛЮЧЕНИЕ

Выполненные в соответствии с разработанной методикой расчеты для шести характерных географических точек РФ и близлежащих государств с использованием разработанной методики позволяют сделать вывод о возможности и достаточной эффективности применения солнечной энергии в качестве основного источника для ВД длительного барражирования в диапазоне широт 40-60° при расчетной взлетной массе аппарата > 30 тонн. Дирижабли меньшей размерности могут быть использованы в течение летнего сезона (возможно, с марта по октябрь), либо на них можно использовать ДВС с химическим топливом при ограниченной продолжительности полета.

Определена рациональная высота барражирования. Для барражирования на широтах 60-70° с. ш. она находится в диапазоне 15-17 км, для более южных широт высота минимума потребной удельной мощности возрастает до 18-20 км.

Разработанная методика может быть использована для расчета и оптимизации параметров конкретного ВД с определенной зоной барражирования под заданную полезную нагрузку, величина которой тпн суммируется и учитывается в весовом уравнении (17). Равенство потребной и располагаемой мощности (красный и синий столбик диаграмм) будет соответствовать решению уравнения весового проектирования высотного дирижабля для заданной географической точки и высоты. Также может быть оптимизирована траектория полета ВД в условиях изменения направления и скорости ветра.

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

СПИСОК ЛИТЕРАТУРЫ

1. Маврицкий В.И., Редькин А.В. Перспективы и возможности создания высотных аэростатических платформ - носителей телекоммуникационного и другого целевого оборудования // Труды ЦАГИ. 2009. Вып. 2682. С. 3-9.

2. Stavros P. Androulakakis. Status and Plans of High Altitude Airship (HAA) Program. AIAA Lighter-Than-Air Systems Technology (LTA) Conference 25-28 March 2013, Daytona Beach, Florida.

3. Редькин А.В. Оптимизация параметров комбинированной системы энергоснабжения беспилотного высотного дирижабля для длительного барражирования в северных широтах // Техника воздушного флота. 2015. № 3-4(720-721). С. 58-65.

4. Грумондз В.Т., Семенников Н.В., Яковлевский О.В. Аэромеханика дирижабля. М.: Наука, 2017. 424 с.

5. Grant E.C., Leland M.N. Fundamentals of Aircraft and Airship Design. AIAA education series. 2013.

6. Фомина Н.Н. Атлас форм корпусов дирижаблей // Труды ЦАГИ. 1935. Вып. 238.

7. Семенов В.А. Механика свободного аэростата. М.: ВВИА им. проф. Жуковского,

1959.

8. Фитенгольц Г.М. Основы математического анализа. М.: Наука, 1968.

9. Schmidt D.K., Stevens J., Roney J. Near-space station-keeping performance of a large high-altitude notional airship: preprints. Atmospheric Flight Mechanics Conference and Exhibit. Keystone, CO: Amer. Inst. Of Aeronautics and Astronautics, 2006.

10. Бугаева И.В. Велопауза и режим летней стратосферной циркуляции. Л.: Гидромет-издат, 1969. 98 с.

Civil Aviation High Technologies

Vol. 21, No. 04, 2018

СВЕДЕНИЯ ОБ АВТОРАХ

Ковалев Игорь Евгеньевич, доктор технических наук, профессор, заместитель генерального директора ФГУП «ЦАГИ» - начальник комплекса управления научными проектами, igor.kovalev@tsagi.ru.

Маврицкий Владимир Иванович, кандидат технических наук, главный научный сотрудник научно-исследовательского отделения № 10 ФГУП «ЦАГИ», mavritsky@tsagi.ru.

Редькин Андрей Владимирович, аспирант, ведущий инженер научно-исследовательского отделения № 10 ФГУП «ЦАГИ», redy68@rambler.ru.

Качарава Ираклий Нугзарович, начальник сектора научно-исследовательского отделения № 11 ФГУП «ЦАГИ», Iraklykacharava@rambler.ru.

METHODOLOGY OF ESTIMATING THE POSSIBLE APPLICATION OF SOLAR ENERGY TO PROVIDE CONTINUOUS AIR STANDING PATROL OF HIGH ALTITUDE AIRSHIP OVER THE GIVEN GEOGRAPHICAL POINT

Igor E. Kovalev1, Vladimir I. Mavritsky1, Andrey V. Redkin1, Irakly N. Kacharava1

1Central Aerohydrodynamic Institute named after prof. N.E. Zhukovsky

Zhukovsky, Russia

ABSTRACT

The task solution of optimizing the parameters for an unmanned high-altitude airship and its power system in order to obtain the maximum of air standing patrol duration over a given geographic point requires an effective comparative analysis of the applicability of solar insolation as the main energy source for a power plant under various altitude, latitudinal, seasonal and wind conditions. The availability of a large number of input parameters for designing such a vehicle leads to the conclusion that it is necessary to use the universal specific energy characteristics of a high-altitude airship and its power plant to simplify the analysis and evaluation of technical solutions effectiveness. The main factor determining the power consumption of the power plant is wind flow speed in a given air standing patrol area. It is also known that solar energy is the most affordable source for long-term energy supply. In most cases, the analysis of solutions is reduced to assessing the ratio between the impact of two natural phenomena, solar energy and wind energy, on the high-altitude airship, for this purpose a simplified model of an airship was designed, which can also be applied to a solar energy airplane. In this work, the value of the available and required power per unit of the aircraft mass is selected as the main performance criteria. The available power is determined for any geographical point with the parameters of latitude, altitude and season (day of the year). The required power is determined by the flight altitude, speed and wind direction. Comparing the values of specific characteristics in space-time coordinates, it is possible to perform an effective analysis of the excess or shortage of available solar energy at a given geographical air standing patrol zone, for different altitudes and days of the year. Using the parameters related to aircraft take-off mass, we are also able to define the required dimensionality of the aircraft to carry out continuous aircraft air standing patrol at a given space point.

Key words: high-altitude airship, power system, solar energy, wind flow, specific characteristics, relative required and available power, space-time coordinates.

REFERENCES

1. Mavritsky, V.I., Redkin, A.V.

Vol. 21, No. 04, 2018

Civil Aviation High Technologies

2. Stavros, P.A. (2013). Status and Plans of High Altitude Airship (HAA) Program. AIAA Lighter-Than-Air Systems Technology Conference 25-28 March 2013, Daytona Beach, Florida.

3. Redkin, A.V. (2015). Parametres optimisation of the combined power system for a high altitude long endurance unmanned airship in northern latitudes. Tekhnika vozdushnogo flota [Technics of air fleet], no. 720-721(3-4), pp. 58-65.

4. Grumondz, V.T., Semenchikov, N.V. and Yakovlevsky, O.V. (2017). Aeromechanika dir-igabllya [Aeromechanics of airship]. Moscow: Nauka.

5. Grant, E.C. and Leland, M.N. (2013). Fundamentals of Aircraft and Airship Design. AIAA education series.

6. Fomina, N.N.

7. Semenov, V.A. (1959). Mekhanika svobodnogo aehrostata [Mechanics of Free Aerostat]. Moscow: Zhukovsky Air Force Academy.

8. Fitengolts, G.M. (1968). Osnovy matematicheskogo analiza [Fundamentals of Mathematical Analysis, Moscow]. Moscow: Nauka.

9. Schmidt, D.K., Stevens, J. and Roney, J. (2006). Near-space station-keeping performance of a large high-altitude notional airship. Preprints. Atmospheric Flight Mechanics Conference and Exhibit. Keystone, CO: Amer. Inst. Of Aeronautics and Astronautics.

10. Bugaeva, I.V. (1969). Velopauza i rezhim letnej stratosfernoj cirkulyacii [Wind pause and summer stratospheric circulation regime]. Leningrad: Gidrometizdat, 98 p.

INFORMATION ABOUT THE AUTHORS

Igor E. Kovalev, Doctor of Technical Sciences, Professor, Deputy of General Director of FSUE "TsAGI", Head of Project Management Complex, kovalev@tsagi.ru.

Vladimir I. Mavritsky, Candidate of Technical Sciences, Chief Research Engineer, Research Department 10, FSUE "TsAGI", mavritsky@tsagi.ru.

Andrey V. Redkin, Postgraduate, Senior Engineer, Research Department 10, FSUE "TsAGI", redy68@rambler.ru.

Irakly N. Kacharava, Head of Sector, Research Department 11, FSUE "TsAGI", Iraklykacharava@rambler.ru.

Поступила в редакцию 27.04.2018 Received 27.04.2018

Принята в печать 17.07.2018 Accepted for publication 17.07.2018