Belgorod, Belgorod, Russian Federation
Belgorod, Russian Federation
V stat'e priveden metod ocenivaniya osobennostey shemnyh resheniy raspredelennnoy vneshney molniezaschity vysokoenergeticheskih ob'ektov. V osnove metoda lezhit opredele-nie udel'nogo raspredeleniya liniy funkciy potoka (LFP) pri resheniyah kraevyh zadach dlya sistem molniezaschity (SMZ) ob'ektov, nahodyaschihsya v elektrostaticheskom pole grozovogo oblaka. Na osnove matematicheskogo modelirovaniya i vychislitel'nyh eksperimentov etot metod aprobirovan pri issledovaniyah shemnyh resheniy dlya raspredelennoy SMZ vysoko-energeticheskih ob'ektov. Issledovalis' shemy vneshnih SMZ kak v vide sostavnyh molnie-priemnikov tipa «metall-izolyator-metall» (MIM), tak i shemy s raspredelennymi me-tallicheskimi molniepriemnikami. Na osnove analiza rezul'tatov vychislitel'nyh ekspe-rimentov dlya kraevyh zadach pokazano, chto sistemy MIM yavlyayutsya bolee effektivnymi v sravnenii s metallicheskimi molniepriemnikami (rezul'taty otlichayutsya v desyatki raz).
elektrostaticheskoe pole, molniepriemniki, molniezaschita, matematicheskoe modelirova-nie, vychislitel'nyy eksperiment
Введение. Принято считать, что основу электрической природы молнии установили исследователи в 18 веке, а именно, американский Бенджамин Франклин (изобретатель «громоотвода», 1775 г.) и российские Михаил Ломоносов и Георг Рихман (погиб при экспериментах по изучению молнии, 1753 г.). Однако и в настоящее время, как заметил Константин Богданов, доктор биологических наук, кандидат физико-математических наук, относительно физики молнии и молниезащиты имеется больше вопросов, чем ответов.
Например, что касается «громоотвода» Франклина в виде «металлического заостренного стержня», находящегося над объектом защиты (ОЗ), т.е. строением, то в [1] (работе Э.М. Базеляна и Ю.П. Райзера) относительно Останкинской телебашни (высотой 540 м) отмечается следующее: «Молниеотвод не в состоянии защитить себя, – звучит анекдотически, но такова реальность!» Авторы работы [1] отмечают систематически проявляющийся эффект, связанный с ударами молний не только в вершину телебашни, а и в ее конструкцию, причем на уровнях до 200 м.
Если представить, что высота бетонной части телебашни с условной крышей порядка 385 м и эту часть можно условно считать ОЗ. Тогда выше ОЗ идет в основном металлическая антенная часть (155 м) с флагштоком на вершине. Согласно [1] антенную часть телебашни можно считать молниеприемником внешней СМЗ, которая поражается молнией от вершины до нижней части, причем молниеприемник находится над ОЗ, т.е. является по Франклину «громоотводом» для строения. В этом случае Э.М. Базеляном и Ю.П. Райзером фактически зафиксирован эффект возможности попадания молнии в конструкцию молниеприемника, а не только в его вершину. Учитывая, что это не единичный факт, указывающий на существующие противоречия в рамках теоретических предположений об ударах молнии только в заостренные вершины молниеприемников с максимальной напряженностью поля на их остриях, однако на практике наблюдаются и другие случаи [1]. В целом продолжается создание как теоретической основы физических процессов возникновения и развития молнии, так и математических моделей для исследования различных аспектов молнии, в том числе и внешних СМЗ объектов.
Краткий обзор математических моделей и методик с учетом особенностей СМЗ для энергообъектов, в том числе, зданий и сооружений. К высокоэнергетическим объектам относятся ракеты-носители, находящиеся на стартовых площадках (столах) космодромов. Удары молний для них представляют большую угрозу, например, в [2] молния миновала центральный стержневой молниеприемник СМЗ, расположенный на высоте значительно превосходящей шаттл на пусковом устройстве, а ударила в стартовый стол (типа как в [1]). В другом случае удар молнии в ракету-носитель Falcon 9 (2016 г.) компании SpaceX с телекоммуникационным спутником привел к взрыву и это существенно затормозило продвижение проекта Internet.org, в рамках которого планировалось подключение к интернету более миллиарда людей, удаленных от цифровых сетей [3]. При этом взрыв ракеты-носителя, например, типа VLS-1 на старте космодрома Алкантара, Бразилия, (2003 г.), мгновенно превратил ракету в фрагменты, а пусковую систему в массу искореженного металла [4]. Эта авария приостановила развитие бразильской космической программы на многие годы [5].
Для исследований физики молнии и СМЗ следует отметить один из важных обзоров математических моделей, методов и методик [6] с учетом последних разработок в этой области. В обзорах [7, 8] по СМЗ для ракетно-космической техники представлены не только математические модели, но и особенности СМЗ, их характеристики, из которых следует выделить эффективные внешние СМЗ составного типа на базе МИМ.
В [9] отмечается, что для изучения нелинейных процессов, связанных с исследованиями ударов молний в объекты, разработаны и доступны математические модели, представленные в работах Ericsson, Dellera и Garbagnati, Рizk, Becerra и Cooray и др., в которых определяются лидерные каналы разрядов. В [9] также обобщены особенности нелинейных моделей и приведены примеры расчетов, выполненных авторами Dellera и Garbagnati для линий электропередач (ЛЭП), в которых траектории каналов разряда, представлены условно сходящимися линиями на ЛЭП с учетом рельефа местности. Для исследования реального атмосферного электрического разряда, т.е. разветвления молнии возможно применение, например, математических моделей на основе фрактальных подходов [10].
Вместе с тем на практике для определения зон защиты молниеотводов используются весьма упрощенные модели, а именно, электрогеометрические методы (по защитному углу), а также Rolling Sphere Model, т.е. на базе катящейся сферы. Эти подходы в основе многих стандартов, например, США [11] и других стран (в том числе, для РФ, но с некоторыми особенностями).
Упрощенные электрогеометрические методы, другие математические модели и методики, а также эмпирические зависимости рассмотрены, например, в [12, 13] с указанием их недостатков и ограничений. Например, отмечается, что в целом их невозможно применять на практике для строений, как распределенных на относительно большой площади земли, так и для высотных зданий и сооружений.
Анализ источников, например, [1, 6–13] и др. показал, что отсутствуют математические модели и результаты расчетов об определении особенностей внешних СМЗ на базе молниеприемников типа МИМ как для высокоэнергетических объектов, так и для других типов объектов. Вместе с тем при отсутствии в открытой печати данных математического или физического моделирования и др., следует заметить, что идет развитие этих внешних СМЗ с МИМ, причем их применение можно обнаружить на космодромах в США, Франции, Индии, на совместном российско-французском космодроме в Куру (Французская Гвиана) и др.
Цель работы связана с представлением метода оценивания схемых решений для распределенной внешней СМЗ (с применением молниеприемников типа МИМ или металлических) для высокоэнергетических объектов. За основу принимается подход [14, 15] (Потапенко А.Н., Кумара У. и др.), связанный с математическим моделированием и вычислительными экспериментами при решениях краевых задач для исследований распределенных СМЗ с различными типами молниеприемников.
Особенности предлагаемого метода на основе определения удельного распределения ЛФП при исследованиях СМЗ объекта в условиях грозового облака. Согласно современным представлениям [16] математическое моделирование и вычислительный эксперимент, позволяют исследовать объекты в условиях «виртуальной реальности». В целом это необходимо для познавательных исследований в рамках современных информационных технологий, а именно, для получения новых знаний, их накоплению (архивированию), а также для применения в других областях знаний, например, для процессов, совершенно различных по физической природе, но описываемых одинаковыми по форме уравнениями в частных производных.
Следует отметить, что для распределенных внешних СМЗ (в виде составных молниеприемников типа МИМ и обычных металлических) с объектами в рамках применения предложенного метода определялось удельное распределение ЛФП при решениях краевых задач. На этой основе оценивались различные схемные решения для внешней СМЗ. Расчет электростатических полей относительно грозового облака выполнялся на основе эллиптического уравнения относительно электрической функции потока μ (аналогичного уравнению Лапласа относительно потенциала ψ). За основу решения уравнения в частных производных был принят метод конечных разностей с учетом стандартных приемов и с использованием соответствующих процедур. При этом для дискретной модели были разработаны разностные схемы, алгоритмы и программы для решения краевых задач с учетом граничных условий.
Главная особенность этого метода связана с определением удельного распределения ЛФП при решении краевой задачи. Это решение представим на упрощенном примере для двумерной области D(x, y). Поле распределения μ определяется из уравнения вида:
Δ μ = 0. (1)
Вместе с тем необходимо заметить, что граничные условия для вычисления поля μ из уравнения (1) и соответственно поля ψ из уравнения Лапласа существенно различаются, т.е. при расчетах функций μ и ψ отсутствует идентичность граничных условий. На базе поля функции потока μ определяются линии семейств
μi(х, у) = const с заданным шагом, причем каждая из этих линий μi(х, у) = const для потенциальных полей начинается на одной из проводящих границ области D(x, y) и заканчивается на другой проводящей границе D(x, y) либо наоборот.
При определении ЛФП в виде μi(х, у) = const записываются граничные условия исследуемых схем моделей СМЗ с объектами в D(x, y). В общем случае для расчета электростатического поля относительно μ между облаком и поверхностью земли считаем, что на грозовом облаке, т.е. на верхней границе Gi в D(x, y), задается условие вида:
∂μ/∂n = 0, (x, y) Î Gi (2)
Тогда на поверхности земли, т.е. на нижней границе Gj в D(x, y) также записывается граничное условие аналогичное (2), т.е. в виде:
∂μ/∂n = 0, (x, y) Î Gj. (3)
В этом случае считается, что границы Gi и Gj являются проводящими в D(x, y), а также считается, что даже при наличии грозового облака в исследуемой D(x, y) отсутствуют как коронные разряды, так и распределенные электрические заряды.
Затем на одной из границ, например, на Gk (перпендикулярной границам Gi и Gj в D(x, y)) записывается условие вида:
μ(x, y) = μ01, (x, y) Î Gk . (4)
Здесь μ01 = const для потенциальных полей и задается в относительных величинах.
Тогда на противоположной границе по отношению к Gk, т.е. на Gn условие следующее:
μ(x, y) = 0, (x, y) Î Gn . (5)
Границы Gk и Gn вводятся условно для ограничения области D(x, y), т.е. для сведения внешней краевой задачи к внутренней, например, с учетом подхода [17, 18]. Этот метод не применим к внешним краевым задачам.
Вместе с тем при определении электростатического поля из уравнения Лапласа относительно ψ на проводящих границах Gi и Gj задаются граничные условия как в [19]:
ψ (x, y) = ψк0, ψ (x, y) = 0 Î Gi и Gj , (6)
Также вместо условий (4) и (5) относительно Gk и Gn для уравнения Лапласа по ψ задаются граничные условия также как в [19]
∂ψ/∂n = 0, (x, y) Î Gk и Gn. (7)
Граничное условие вида (7) для ψ в электрических полях – это условие не протекания тока (для электростатических полей название иное). Тогда из уравнения Лапласа при соответствующих граничных условиях определяются эквипотенциальные линии в виде семейства ψi(х, у) = const. Важно заметить, что линии семейств
μi(х, у) = const и ψi(х, у) = const, (8)
в потенциальных полях ортогональны. Таким образом, для краевых задач, решаемых с использованием математического моделирования и вычислительного эксперимента граничные условия для нахождения функций μ и ψ должны различаться.
Для определения компонент напряженности поля Ex и Ey относительно функций ψ и µ записываются зависимости вида:
Ex = – ∂ψ/∂x = – ∂µ/∂у; (9)
Ey = – ∂ψ/∂y = ∂µ/∂х, (10)
причем Ei = grad ψ = (ψi+1 - ψi) / (ri+1 - ri) = Δ ψ / Δ r (r – соответствующая координата по нормали, например, к ψi(х, у) = const). При этом Ei = grad µ практически не применяется.
Дополнительная особенность этого метода. Для решения краевых задач с учетом в исследуемой, например, области с(x, y) реальных соотношений высот между грозовым облаком, поверхностью земли и объектами Ωi возникает необходимость выполнить масштабирование. Известно, что в современных программных продуктах автоматизирован процесс наложения сетки на исследуемые Ωi. Также предусмотрена возможность изменения шага сетки, например, в краевых областях (где напряженность электрического поля высокая) с применением значительно меньших тетраэдрических элементов (например, для метода конечных элементов).
Учитывая, что в исходной задаче имеется возможность определения линий семейства вида (8), тогда за основу принимается подход [20], при котором они применяются как, например, показано на рис. 1, а для масштабирования.
|
а |
б |
Рис. 1. Пример моделирования с использованием подхода [20]; определяются линии семейства
μi(х, у) = const и ψi(х, у) = const (а); выделяется новая область abcd (б)
В исследуемой области (см. рис. 1, а) выделяется некоторая область abcd (обозначена пунктиром), подлежащая масштабированию (см. рис. 1, б, для более детального изучения объекта) на основе применения ортогональной криволинейной сетки. Однако в исходной задаче применяется модифицированная методика по отношению к [20] связанная с тем, что в решаемых краевых задачах применяется универсальная ортогональная прямоугольная сетка для выделения прямоугольной области (аналогичной abcd) для масштабирования.
Отличительной чертой в модифицированной методике с учетом перехода от одной области к другой является следующее. Необходимо соблюдать условие, при котором верхняя линия из семейства ψi(х, у) = const в исследуемой области
с(x, y) должна быть прямой. Это связано с необходимостью выделения некоторой части электростатического поля, где исключено влияние исследуемого объекта на верхнюю границу ψk0(х, у) = const, а также учитывается степень масштабирования. В этой связи вместо исходной с(x, y) появляется новая область С(x, y). Принципиальное отличие заключается в том, что вместо реальной верхней границы, например, в виде облака в с(x, y) появляется новая С(x, y) с верхней границей в виде некоторой прямой линии ψk0(х, у) = const из семейства эквипотенциальных линий. Учитывая, что в потенциальных полях важны не их абсолютные значения, а соотношения этих функций, поэтому на границах С(x, y) задаются в относительных величинах значения требуемых функций.
Таким образом, принципиальная особенность этого метода связана с необходимостью использования удельного распределения ЛФП в рамках решения краевой задачи для оценки схемных решений относительно распределенной молниезащиты высокоэнергетического объекта. Оценка схемных решений выполняется на основе сравнения результатов расчета удельного распределения ЛФП как для внешних СМЗ, так и для защищаемых объектов. Установлено, что чем больше величина удельного распределения ЛФП на внешних СМЗ и чем меньше она на объекте защиты, тем выше эффективность исследуемого схемного решения.
Основные примеры для 3-х вариантов краевых задач.
Вариант 1. В [14] представлены данные расчетов удельного распределения ЛФП для объекта защиты и одного из распределенных металлических СМЗ (см. рис. 2, а и рис. 2, б) Обозначения в вариантах следующие: 1) χ1i – зависимость удельного распределения ЛФП для объекта защиты, как функция χ1i = f (l*1i, H*2k); 2) χ3i – зависимость удельного распределения ЛФП для молниеприемника СМЗ, как функция χ3i = f (l*1i, H*2k).
Эти зависимости представлены в относительных величинах при изменениях дистанций l*1i между объектом защиты и молниеприемником. Высота объекта защиты H*3 = Hi / H3 = 1, высота молниеприемника СМЗ H*21 = Hi / H3 = 2, дистанция l1i* = li / H3, причем H3 – исходная высота объекта защиты принята за основу для нормирования; знак(*) – указывает на нормированные параметры или величины.
Анализ показывает (см. рис. 2), что зависимость удельного распределения ЛФП для объекта защиты χ1i меньше чем χ3i в представленном дипазоне изменения l*1i. Однако при увеличении дистанции l1i* более чем l*1i = 4 величины χ1i и χ3i стремятся к некоторому условному пересечению за счет нелинейностей. Также из анализа зависимостей на рис. 2 следует, что, например, при минимальном значении χ1min = 0,03 (l*1i = 0,8) наблюдается, что значение χ3i = 0,67, т.е. χ3i > χ1i, причем более чем в 20 раз при l*1i = 0,8. Однако при максимальном значении χ3max = 0,735 (l*1i = 1,6) наблюдается χ1i = 0,09, т.е. также удельное распределения ЛФП для СМЗ имеет зависимость χ3i > χ1i, но в пределах до 8 раз (за счет нелинейности). При дальнейшем увеличении дистанции l*1i более чем l*1i = 1,6 величина χ3i начинает стремительно уменьшаться. Эти данные необходимо учитывать при проектировании систем молниезащит.
|
а |
б |
Рис. 2. Изменения удельного распределения ЛФП (χ) при H*21 = const и
l*1i = var для χ1 – объекта защиты (a) и для χ3 – молниеприемника СМЗ (б)
Вариант 2. В случае, если увеличить высоту молниеприемника СМЗ до величины H*22 = 3 [14]. (Например, при условии H3 = 50 м исследуемая модель будет в пределах высот H*22 = Hi / H3 = 3, что характерно для двух основных молниеприемников (H22 =150 м) на новом космодроме Восточный). Тогда значение χ3max увеличится до χ3i = 0,87 (l*1i = 1,6), а значение для объекта зашиты уменьшится до χ1i = 0,045, т.е. при H*22=3 соблюдается условие χ3i > χ1i, но уже в пределах до 19 раз. Сравнение этих данных с результатом при H*21 = 2 показывает, что увеличение высот молниеприемников СМЗ естественно приводит к увеличению зависимости удельного распределения ЛФП χ3i и соответственно к уменьшению зависимости удельного распределения ЛФП χ1i, т. е. в целом приводит к повышению эффективности исследуемой схемы.
Вариант 3. Вместе с тем, если исследовать внешнюю СМЗ в виде составных молниеприемников типа МИМ [15]. Тогда в условиях как для варианта 1, но при этом высота каждого элемента составного молниеприемника следующая. Высота нижней части (решетчатчатая металлическая башня) H*023 = 1,6; высота средней части (изолятор цилиндрического типа) H*022 = 0,3 и высота верхней части (металлический заземленный стержнь) H*021 = 0,1, т.е. в целом высота H*21 = H*023 + H*022 + H*021 = 2 (как в варианте 1). В этом варианте значение χ3max увеличится до χ3max = 0,875 (l*1i = 1,6), а величина χ1i, для объекта защиты уменьшится до χ1i = 0,035, т.е. соблюдается условие χ3i > χ1i, но уже в пределах до 25 раз. Сравнение данных варианта 3 с результатами при H*21 = 2 (вариант 1) показывает, что в исследуемых условиях это приводит к увеличению зависимости удельных распределений ЛФП для СМЗ и уменьшению χ1i для объекта, причем без увеличения высоты молниеприемника (для варианта 3 отношение χ3max / χ1i стремиться к 25, а для варианта 1 отношение χ3max / χ1i к 8 раз), т.е. в целом повышается эффективность исследуемой схемы. Следует заметить, что вариант 3 в сравнении с вариантом 2 имеет практически одинаковые результаты по эффективности, но вариант 3 позволит дать значительную экономию затрат, например, для космодрома Восточный как по металлоемкости для каждой конструкции внешней СМЗ, так и по расходам на фундаменты для каждой СМЗ и др.
Выводы. Предлагаемый метод, основанный на решениях краевых задач, позволит дополнить существующие электрогеометрические методы в случаях, где они не дают однозначную оценку, а для других случаев, например, позволит получить единственное решение в виде оценки исследуемых схем внешней СМЗ с молниеприемниками типа МИМ для высокоэнергетических объектов в условиях грозового облака. Выделенные основные особенности предлагаемого метода по определению электростатического поля относительно электрической функции потока и удельного распределения ЛФП для внешних СМЗ с объектом позволят существенно упростить процесс освоения этого метода. На примерах показаны особенности применения методики анализа при оценках схемных решений распределеннной молниезащиты объекта в условиях грозового облака.
Развитие этого метода даст возможность расширить его применение для объектов энергетики (ЛЭП, подстанций, ветроэнергетических установок и др.), для зданий и сооружений, для телекоммуникационных систем и др.
*Исследование выполнено при финансовой поддержке РФФИ в рамках международного научного проекта №17-58-45155.
1. Bazelyan E.M., Rayzer Yu.P. Fizika molnii i molniezaschity. M.: Fizmatlit, 2001. 320 s.
2. NASA Facts AC 321/867-2468, Lightning and Space Program, FS-1998-08-16-KSC, Au-gust 1998.
3. Raketu Falcon 9 vzorvala sharovaya mol-niya? – URL: http://paranormal-news.ru/news/raketu_falcon_9_vzorvala_sharovaja_molnija/2016-09-27-12679 (data obrasche-niya 01.03.2017).
4. Lantratov K. Brazil'skaya raketa vzo-rvalas' na starte // Kommersant' (25 avgusta 2003), s. 11.
5. 50 let brazil'skogo osvoeniya kosmosa – URL: http://genby.livejournal.com/420452.html (data obrascheniya 01.03.2017).
6. Rakov V.A., Rachidi F. Overview of Re-cent Progress in Lightning Research and Light-ning Protection, IEEE Transactions on Electro-magnetic Compatibility. 2009. Vol. 51. №3. Pp. 428-442.
7. Kumar U. Lightning Protection of Satellite Launch Pads, in Lightning Protection, edited by V. Cooray, IEE Power and Energy Series 58, United Kingdom: IET, 2010. R. 789–819.
8. Kumar U., Nagabhushana G.R. Analysis of Lightning Protection System for Indian Satel-lite Launch Pad, National Interagency coordina-tion group Lightning Conference, Orlando, Flori-da, USA, 2000. 12th April.
9. Malgesini R., Fenomeni di fulminazione diretta ed indiretta di linee elettriche: model-lazione e sperimentazione in laboratorio, Rap-porto RSE 09004921 – URL: http://www.rseweb.it/documenti.page?RSE_originalURI=/documenti/documento/2822&RSE_manipulate-Path=yes&country=ita#sthash.CenPpSEj.dpuf (data obrascheniya 01.03.2017).
10. Yazhong Xu, Mingli Chen, Striking dis-tance calculation for flat ground and lightning rod by a 3D self-organized Leader Propagation Model, 2012 International Conference on Light-ning Protection, Vienna, Austria, 2012.
11. NFPA 780: Standard for the Installation of Lightning Protection Systems, 2004. 108 p.
12. D’Alessandro F.A., Gumley J.R., ‘Col-lection Volume Method’ for the placement of air terminals for the protection of structures against lightning, // Electrostatics. 2001. № 50. R. 279–302.
13. Preventing Direct Lightning Strikes, Rev. B - Carpenter, et al. - March 2014. – URL: http://www.lightningprotection.com/ (data ob-rascheniya 01.09.2017).
14. Potapenko A.N., Kumar U., Potapenko T.A., Shtifanov A.I. Issledovanie raspre-delennoy molniezaschity startovyh ploschadok kosmodromov na osnove vychislitel'nogo eks-perimenta // Fundamental'nye issledovaniya. 2015. № 11-4. S. 727–734.
15. Shtifanov A.I., Potapenko A.N., Po-tapenko T.A., Kumar U. Vychislitel'nye eks-perimenty s vozmozhnost'yu opredeleniya os-novnyh osobennostey sostavnyh molniepri-emnikov // Fundamental'nye issledovaniya. 2015. № 11-6. S. 824–831.
16. Plohotnikov K.E. Matematicheskoe modelirovanie i vychislitel'nyy eksperi-ment: metodologiya i praktika // Intellektu-al'nye sistemy. Teoriya i prilozheniya. 2009. T. 13. № 1-4. S. 5–32.
17. Potapenko T.A., Kanunnikova E.A., Potapenko A.N. Research of 3-Dexterior bound-ary problems related to electric fields in atmos-phere by inversion method // Electric Power Sys-tems Research. 2014. Vol. 113. P. 10–14.
18. Potapenko T.A., Kanunnikova E.A., Potapenko A.N. Research of 3-D Exterior Boundary Problems Related to Electric Fields in Atmosphere by Inversion Method // 2012 Inter-national Conference on Lightning Protection (ICLP), Vienna, Austria. 2012.
19. Demirchyan K.S., Neyman L.R., Korov-kin N.V., Chechurin V.L. Teoreticheskie osnovy elektrotehniki: V 3- t. Uchebnik dlya vuzov. Tom 3. 4-e izd. SPb.: Piter, 2003. 377 s.
20. Ryazanov G.A. Opyty i modelirovanie pri izuchenii elektromagnitnogo polya. M.: Nauka, 1966. 208 s.
