ro en

Электронный журнал №2(43)2019

"ПРОБЛЕМЫ РЕГИОНАЛЬНОЙ ЭНЕРГЕТИКИ"

СОДЕРЖАНИЕ

0	ТЕПЛОЭНЕРГЕТИКА
	Аннотация:
1	Оптимальные режимы работы волнового компрессора микротурбины Авторы: Левихин А.А.1, Волобуев И.А.2 1 - Балтийский государственный технический университет «ВОЕНМЕХ» им. Д.Ф. Устинова, Санкт-Петербург, Российская Федерация, 2 - Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики (Университет ИТМО), Санкт-Петербург, Российская Федерация
	Аннотация: Работа посвящена проблеме повышения эффективности сжатия воздуха в компрессоре микротурбинной установки. Одним из современных направлений исследований является применение сверхзвукового волнового компрессора, в котором сжатие воздуха происходит не в лопаточной машине, в ударно-волновой структуре, подобно тому, как это происходит в сверхзвуковом воздухозаборнике летательных аппаратов. Цель работы - найти оптимальные режимы работы волнового компрессора с наименьшими потерями полного давления. Поставленная цель достигается тем, что изучаются режимы работы волнового компрессора при сжатии воздуха в системе двух ударных волн одного направления и в центрированной изоэнтропической волне. Аналитически исследованы области существования ударно-волновых структур с отраженным разрывом: волной разрежения, ударной волной и слабым разрывом. Слабый отраженный разрыв соответствует наибольшей эффективности сжатия. Получены следующие наиболее важные результаты. Ударно-волновой компрессор с двумя ударными волнами одного направления может работать на оптимальном режиме в диапазоне чисел Маха на входе от 2,089 до 2,46. При этом может быть получена степень сжатия в компрессоре 5,76. Для получения более высокой степени сжатия в одной ступени необходимо переходить к трем и более ударным волнам, что конструктивно сложно. Центрированная волна сжатия может обеспечить степень сжатия в одной ступени до 10-12. Число Маха на входе в колесо компрессора должно составлять от 3,2 до 3,48. При более высоких степенях сжатия образуются нерегулярные ударно-волновые структуры. Численные эксперименты показали, что теоретическая эффективность волнового компрессора выше, чем у традиционного центробежного и осевого, при скорости набегающего потока на входе в компрессор, соответствующей числу Маха, больше 2,2. Значимость результатов заключается в том, что определены границы оптимальных режимов работы при использовании двух типов ударно-волновых структур.
	Ключевые слова: микротурбина, волновой компрессор, ударная волна, изоэнтропическая волна сжатия, догоняющие ударные волны, ударные волны одного направления DOI: https://doi.org/10.5281/zenodo.3367056
2	Численное исследование теплообмена в канале с продольным градиентом давления Авторы: Цынаева А.А. Самарский государственный технический университет Самара, Российская Федерация
	Аннотация: Работа посвящена исследованию теплообмена в каналах с продольным градиентом давления. Целью работы является разработка методов повышения интенсивности теплообмена в каналах с продольным градиентом давления. Поставленная цель была достигнута за счёт разработки оригинальной конструкции поверхностного интенсификатора, выполненного в виде подковообразной лунки, а также проведением исследования теплообмена в канале с продольным градиентом давления и лунками различной геометрии и расположения. Было установлено, что численное решение с использованием k-epsilon модели турбулентности показывает недостаточно стабильные результаты, в том числе, значения локальных коэффициентов теплоотдачи могут до 2...2.5 раз превышать их действительные значения для одиночной лунки. При этом использование k-w-sst модели турбулентности показывает более стабильные результаты. Наиболее значимыми результатами, полученными в работе являются следующие: выявлено, что применение лунок, их геометрия и расположение оказывают значительное влияние на интенсивность теплообмена в канале с продольным градиентом давления; разработана конструкция оригинальных подковообразных лунок, применение которых позволяет интенсифицировать теплообмен на 13 ... 47% в зависимости от числа Рейнольдса при их расположении в три ряда в коридорном порядке по сравнению с полусферическими лунками в шахматном порядке; при сравнении эффективности применения разработанных лунок с полусферическими лунками, установленными в три ряда в коридорном порядке, рост теплоотдачи составил 21...51% при изменении числа Рейнольдса от 3000 до 15000. Исследования выполнялись методом численного решения уравнения сохранения массы (уравнения неразрывности), уравнения сохранения количества движения (в форме Навье-Стокса), уравнения сохранения энергии. Численные исследования проведены в достаточно широком диапазоне изменения чисел Рейнольдса Re=3000...70000. Исследование проводилось с использованием численных методов, адекватность которых была подтверждена сопоставлением с опытными данными.
	Ключевые слова: системы охлаждения, лопатки турбин, теплообмен, градиент давления, моделирование, канал, лунки. DOI: https://doi.org/10.5281/zenodo.3367076
3	Энергетический мониторинг инновационных энерготехнологий переработки растительного сырья Авторы: Бурдо О.Г.1, Гаврилов А.В.2, Кашкано М.В.1, Левтринская Ю.О.1, Сиротюк И.В.1, Пилипенко Е.А.1, Терзиев С.Г. 3 1 Одесская национальная академия пищевых технологий 2 Академия биоресурсов и природопользования «КФУ им. В.И. Вернадского» 3 Публичное акционерное общество «Енни Фудз» г. Одесса, Украина
	Аннотация: Анализ многочисленных научных работ показывает, что интенсивное развитие инновационных образцов техники опережает уровень развития методологических основ энергетического менеджмента. Отсутствуют чёткие сравнения энергоэффективности электротехнологий и теплотехнологий, поскольку используются различные виды энергии; отсутствуют объективные показатели эффективности использования энергии в различных технологиях обезвоживания сырья. В данном исследовании предложена методология, в основе которой положена гипотеза, что объективные результаты при сравнении эффективности использования энергии при переработке сырья можно получить на основе системного анализа всей цепи конверсии энергии от топлива до готового продукта. Цель работы – экспериментально доказать объективность данной гипотезы. Для достижения этой цели предлагается использовать показатель доли энергии топлива в готовом продукте и количество удаленной влаги при сжигании 1 кг топлива, который не зависит от колебания цен на энергоносители, которые могут колебаться и отличаться для разных стран. Выполнен анализ тепловых балансов сушильной и выпарной установок, показано что при одинаковых технических задачах неоптимизированная выпарка в разы эффективней оптимизированного процесса сушки. Приведены структурные модели конверсии энергии при комбинированных процессах получения концентрированных пищевых продуктов. Рассчитаны эффективность использования энергии топлива в традиционных схемах сушки, выпарки, криоконцентрирования. Наиболее важным результатом работы является сравнение этих параметров с данными для инновационных образцов техники, разработанных авторами. Значимость полученных результатов состоит в том, что разработанные авторами выпарные установки не уступают по эффективности традиционным, и позволяют получать концентраты до 90 °brix. Так, для топлива с нефтяным эквивалентом 40 МДж на 1 кг традиционные сушилки могут удалить не более 3 кг влаги, криоконцентраторы – 20 кг. Достигнуты результаты, соответственно, до 6 и до 100 кг влаги.
	Ключевые слова: энергетический менеджмент, пищевые производства, энергоэффективность, сушка, криоконцентрирование, микроволновое поле, обезвоживание, пищевые концентраты. DOI: https://doi.org/10.5281/zenodo.3367058
4	Расчетно-экспериментальные исследования концентрационных пределов горения попутных нефтяных газов Авторы: Бачев Н.Л., Матюнин О.О., Бульбович Р.В., Бетинская О.А., Шилова А.А. Пермский национальный исследовательский политехнический университет Пермь, Российская Федерация
	Аннотация: Вопрос об утилизации попутных нефтяных газов на малодебитных месторождениях с использованием газотурбинных энергоустановок в последние годы является весьма актуальным. Для разработки универсальных камер сгорания таких установок, которые бы обеспечивали устойчивое горение разнородных по составу и сильно забалластированных газов на разных месторождениях, требуется достоверная информация по пределам их горения в реальных условиях эксплуатации. Анализ составов попутных нефтяных газов различных месторождений Российской Федерации показал, что содержание негорючих компонентов (азот, углекислый газ, пары воды) может достигать 90 %, тяжелых углеводородов (пентан и выше) – до 10 %, сероводорода – до 6 %. Целью данной работы является теоретическое определение и экспериментальное подтверждение концентрационных пределов горения забалластированных попутных нефтяных газов, а также разработка рекомендаций по выбору коэффициента избытка воздуха в первичной зоне камеры на этапе проектирования. Поставленная цель достигается теоретическим определением коэффициента избытка воздуха на верхнем и нижнем пределах горения с использованием методов флегматизации и экспериментальным подтверждением расчетных результатов. Впервые экспериментальные результаты получены в процессе сжигания нефтяных газов, отобранных непосредственно на месторождениях, с содержанием балластирующих компонентов 26 % и 40 %. В процессе экспериментальных исследований определялись коэффициенты избытка воздуха на верхнем и нижнем пределах горения и продольные профили температуры в пламени по оси горелки. О положении фронта пламени в камере можно судить по показаниям термопар в зоне пламени и кадрам видеосъемки внутрикамерного процесса. Расчетно-экспериментальные исследования пределов горения забалластированных разнородных по составу нефтяных газов позволили выработать рекомендации по организации стабильного диффузионного или гомогенного горения в первичных зонах камер сгорания утилизационных микрогазотурбинных энергоустановок.
	Ключевые слова: попутные нефтяные газы, концентрационные пределы горения, коэффициент избытка воздуха, балластирующие компоненты, экспериментальная установка. DOI: https://doi.org/10.5281/zenodo.3367068
5	Влияние централизованного и индивидуального отопления многоквартирных зданий на экологию (на примере г. Кишинев, Республика Молдова) Авторы: Тыршу М.С., Лупу М.Л., Посторонкэ С. А. Институт энергетики Кишинев, Республика Молдова
	Аннотация: Целью настоящей работы является аргументация использования централизованной системы обеспечения тепловой энергией исходя из присущей ей безвредности по отношению к общественному здоровью и к окружающей среде, по сравнению с использованием индивидуального отопительного оборудования, смонтированного непосредственно в жилых помещениях муниципия Кишинев. Поставленная цель была достигнута за счет сравнения количества выбросов, образующиеся в результате работы этих двух типов систем отопления. Актуальность предмета исследования исходит из наличия существенного количества потребителей, отключенных от централизованной системы отопления, хотя в последние годы наблюдается резкое снижение такой тенденции. Продолжает пользоваться приоритетом индивидуальное отопление в недавно построенных жилых домах. Статистика, представленная специализированными учреждениями муниципия, говорит о высоком уровне инцидентов отравления угарными газами, в частности окисью углерода, за которым следуют отравления от утечки газа, ожоги, происходящие по причине нарушения правил эксплуатации домашнего газового оборудования, а также летальные исходы. Новыми результатами являются определение разности в значениях выбросов для обоих случаев. Важность полученных результатов состоит в том, что для генерирования тепловой энергии в централизованной системе отопления при использовании технологии когенерации, вредные выбросы в окружающую среду оказываются на 28% меньше, чем при использовании индивидуальных систем обеспечения теплом. Данная констатация позволяет аргументировать явный приоритет использования централизованной системы перед системами индивидуального отопления.
	Ключевые слова: отопление, централизованная система, оборудование, природный газ, общественное здоровье, вредные выбросы, когенерация. DOI: https://doi.org/10.5281/zenodo.3367070
6	Влияние эффективности градирни на показатели экономичности турбоустановки Авторы: Галацан М.П., Кравченко В.П., Киров В.С. Одесский национальный политехнический университет Одесса, Украина
	Аннотация: Анализ мирового электроэнергетического баланса показал важное место паротурбинных установок. Рост их эффективности путем повышения начальных параметров встречает трудности, связанные с со-зданием новых материалов. Повышение эффективности за счет снижения конечного давления связано с совершенствованием системы оборотного водоснабжения. Целью работы является определение влияния эффективности градирни на экономичность работы электростанции. Поставленная цель достигнута за счет разработки математической модели связанных друг с другом башенной градирни и турбоустановки с конденсатором и проведением расчетов для двух типов оросителя градирни: асбестоцементные листы и сетчатые элементы, выполненные из полиэтилена. Для расчета принят энергоблок Ровенской АЭС с турбоустановкой К-1000-5,8/50. Принимая среднемесячную температуру воздуха, в результате вари-антных расчетов определялась температура охлажденной в градирне воды. Далее в зависимости от нее с помощью математической модели конденсатора определяется температура конденсации пара, от ко-торой в свою очередь зависит расход пара в конденсатор, производительность турбоустановки и темпе-ратура охлаждающей воды на выходе конденсатора, то есть на входе в градирню. Таким образом, для каждого месяца были определены температуры охлажденной воды и показатели эффективности работы АЭС: КПД и количество произведенной электроэнергии. Новизна работы состоит в учете взаимодей-ствия градирни и турбоустановки. Значимость полученных результатов состоит в том, что с учетом взаимовлияния турбоустановки и градирни эффективность замены оросителя оказывается выше. В ре-зультате расчетов получено, что при замене оросителя средняя за полгода температура конденсации снизится на 2,34 °С. При этом производство электроэнергии за год увеличится на 41,72 ГВт∙ч. С учетом затрат на приобретение элементов оросителя и его монтаж срок окупаемости модернизации оросителя составит 2,5 года.
	Ключевые слова: система оборотного водоснабжения, математическая модель, градирня, турбоустановка, конденсатор. DOI: https://doi.org/10.5281/zenodo.3367064
7	ЭЛЕКТРОЭНЕРГЕТИКА
	Аннотация:
8	Влияние технических параметров электрических аппаратов систем электроснабжения на их характеристики надежности Авторы: Грачева Е.И. 1, Ившин И.В. 1, Горлов А.Н. 2, Алимова А.Н. 1 1Казанский государственный энергетический университет Казань, Россия 2Юго-Западный Государственный Университет Курск, Россия
	Аннотация: Целью данной работы предлагается новый подход к исследованию эффективности функционирования коммутационных аппаратов в процессе их эксплуатации. Поставленная в работе цель достигнута в результате проведенных экспериментальных и теоретических исследований с помощью тепловизионного контроля. Наиболее важными результатами являются: разработка нового критерия оценки технического состояния аппаратов – коэффициента фактической кратности превышение сопротивления контактов по условиям эксплуатационных режимов и определение его численных для автоматических выключателей, контакторов, магнитных пускателей, рубильников и пакетных выключателей, которые позволяют оценивать их техническое состояние и эксплуатацию, не допускать развитие дефектов коммутационной аппаратуры и аварии, использовать разработанный критерий для вспомогательного регламента сроков про-водимых обследований коммутационной аппаратуры. Значимость полученных результатов состоит в получении по данным проведенных экспериментов значений эмпирических зависимостей начальных сопротивлений контактных соединений аппаратов от величины номинального тока. По данным тепловизионных обследований установлены виды повреждений аппаратов и получены зависимости сопротивлений контактных соединений дефектного оборудования от температуры контактных систем аппаратов при их различных режимах загрузки. Разработаны алгоритм и метод комплексной оценки эффективности функционирования коммутационных низковольтных аппаратов позволяющие установить зависимость между вероятностью времени безотказной работы и сопротивления контактов аппаратов с учетом количества циклов коммутаций. При этом сопротивления контактов рассматривается как случайная функция количества коммутационных циклов. По результатам проведенных экспериментальных исследований для достижения критического значения сопротивления аппарата принимается в среднем трехкратное увеличение начального сопротивления. Проведенные согласно разработанному алгоритму исследования с выявлением функциональной зависимости между параметрами сопротивлений контактов и вероятности времени безотказной работы аппаратов позволяют оценивать техническое состояние оборудования уточнять величину потерь электроэнергии и характеристик надежности систем электроснабжения.
	Ключевые слова: аппараты, безотказность, надежность, электроэнергия, силовые цепи, сопротивления контактных соединений. DOI: https://doi.org/10.5281/zenodo.3367062
9	Машинное обучение системы функционального диагностирования шахтной подъемной машины Авторы: Довбыш А. С., Зимовец В. И., Зубань Ю. А., Приходченко А. С. Сумской государственный университет Сумы, Украина
	Аннотация: Цель работы – повышение точности функционального диагностирования узлов шахтной подъемной машины путем синтеза и реализации метода информационно-экстремального машинного обучения с иерархической структурой данных. Задачи, которые ставились в работе: разработать категорийную модель; синтез на основе ее иерархического алгоритма машинного обучения для системы функционального диагностирования; оптимизация системы контрольных допусков. Известно, что при функциональном диагностировании возникает необходимость анализа большого количества диагностических признаков и классов распознавания, характеризующих не только возможные неисправности, но и промежуточные технические состояния узлов и агрегатов сложной машины. При этом предложенный алгоритм разрабатывается в рамках так называемой информационно-экстремальной интеллектуальной технологи анализа данных, основанной на максимизации информационной способности системы в процессе машинного обучения. Основная идея предлагаемого метода, как и в нейронных сетях, заключается в адаптации в процессе машинного обучения входного математического описания системы функционального диагностирования к максимальной достоверности диагностических решений. Но в отличии от нейронных сетей построенные в процессе информационно-экстремального машинного обучения решающие правила являются практически инвариантными к многомерности как пространства диагностических признаков, так и алфавита классов распознавания. Реализация предложенного метода информационно-экстремального машинного обучения осуществляется на примере функционального диагностирования многоканатной шахтной подъемной машины, являющейся сложным электромеханическим устройством с повышенными требованиями к безопасности эксплуатации. Исследовано влияние системы контрольных допусков на диагностические признаки. Наиболее значимыми результатами является повышение достоверности диагностических решений при использовании иерархического алгоритма машинного обучения системы функционального диагностирования по сравнению с линейным классификатором. Кроме того, построенные по оптимальным геометрическим параметрам гиперсферических контейнеров классов распознавания решающие правила позволяют принимать высокодостоверные диагностические решения в реальном темпе времени.
	Ключевые слова: информационно-экстремальная интеллектуальная технология, машинное обучение, информационный критерий, функциональное диагностирование, шахтная подъемная машина. DOI: https://doi.org/10.5281/zenodo.3367060
10	ВОЗОБНОВЛЯЕМЫЕ ИСТОЧНИКИ ЭНЕРГИИ
	Аннотация:
11	Синтез рациональной конструкции валкового пресса для производства композитного твердого топлива Авторы: Баюл К.В. Институт черной металлургии им. З.И. Некрасова НАН Украины Днепр, Украина
	Аннотация: Целью работы является синтез конструкции валкового пресса для производства альтернативных видов топлива из мелкофракционного энергетически ценного сырья методом брикетирования. Достижение поставленной в работе цели осуществляется путем формирования системного подхода к поиску рациональных конструктивных решений валковых прессов. Исследования проводились на основе анализа зарубежных и отечественных литературных источников, практическом опыте специалистов Института черной металлургии им. З.И. Некрасова НАН Украины в области теоретических исследований, разработки и внедрения прессового оборудования для брикетирования мелкофракционного сырья, применяемого в качестве альтернативных источников энергии. Установлено, что перспективным направлением развития альтернативной энергетики является производство композитных твердых видов топлива методом брикетирования из смеси энергетически ценных мелкофракционных материалов. Показано, что для производства композитных топливных брикетов перспективным является использование валковых прессов. Выявлена проблема, состоящая в том, что на сегодняшний день не существует самостоятельных или интегрированных в известные системы автоматизированного проектирования программных средств, позволяющих с учетом технологических особенностей, проводить разработку и моделирование всего жизненного цикла валковых брикетных прессов. Новизна работы состоит в том, что впервые для создания по модульному принципу рациональной конструкции валкового брикетного пресса для производства композитного твердого топлива применен системный подход, основанный на использовании структурно-параметрического синтеза и анализа. С использованием предложенного системного подхода разработана новая усовершенствованная модель валкового пресса для производства топливных брикетов, в том числе композитных, из энергетически ценных мелкофракционных сырьевых материалов, имеющих малую насыпную плотность (0,2…0,6г/см3). Развитие описанного в работе системного подхода позволит упростить разработку модификаций валковых прессов с учетом особенностей процесса брикетирования.
	Ключевые слова: альтернативные источники энергии, композитное твердое топливо, брикеты, валковый пресс, структурно-параметрический синтез. DOI: https://doi.org/10.5281/zenodo.3367048
12	Оценка выхода биогаза и электроэнергии при моносбраживании навоза крупного рогатого скота и при добавлении фуза в качестве косубстрата Авторы: Полищук Виктор Николаевич, Титова Людмила Леонидовна, Шворов Сергей Андреевич Национальный университет биоресурсов и природопользования Украины Киев, Украина Гунченко Юрий Александрович Одесский национальный университет имени И. И. Мечникова Одесса, Украина
	Аннотация: Целью работы является повышение выхода биогаза и выработки электроэнергии на биогазовых установках за счет использования стимулирующей добавки в виде фуза при переходе от периодической к квазинепрерывной загрузке метантенка. Для достижения поставленной цели решались следующие задачи: определялся выход биогаза из навоза крупного рогатого скота при различных температурных режимах с периодической загрузкой метантенка; оценивался выход биогаза при сбраживании навоза крупного рогатого скота с добавлением фуза; с помощью математической модели прогнозировался выход биогаза для квазинепрерывной загрузке метантенка. Исследования проводились на лабораторной биогазовой установке, состоящей из метантенка полезным объемом 30 л и газгольдера мокрого типа. Выход биогаза фиксировался по поднятию цилиндра-уровнемера мокрого газгольдера с помощью закрепленной на нем шкалы, отградуированной в сантиметрах. Биогаз сжигался на газовой плите. Теплота сгорания биогаза определялась по его элементному составу, который фиксировался газоанализатором. Наиболее важными результатами являются следующие: по данным экспериментальных исследований выхода биогаза при периодическом режиме загрузки с использованием данной модели обеспечивается прогнозирование выхода биогаза для квазинепрерывной загрузки метантенка. В результате экспериментальных исследований установлено, что при моносбраживании навоза максимальный выход биогаза наблюдается на 3 8 сутки брожения, после чего постепенно уменьшается. При добавлении к субстрату 1.3% фуза выход биогаза увеличивается и достигает максимального значения на 20 сутки, после чего плавно уменьшается. Значимость результатов исследований состоит в том, что использование фуза в качестве косубстрата позволит почти в три раза увеличить выход биогаза и выработку электроэнергии, снизить срок окупаемости биогазовой установки мощностью 4 МВт при использовании "зеленого" тарифа до 5.1 года.
	Ключевые слова: биогаз, субстрат, навоз крупного рогатого скота, фуз, сухое вещество, метантенк, биогазовая установка, метановое брожение. DOI: https://doi.org/10.5281/zenodo.3367054
13	Исследование системы молниезащиты ветроэнергетической установки Авторы: 1Потапенко А.Н., 2Кумар У., 3Штифанов А.И., 1Солдатенков А.С. 1Белгородский государственный технологический университет им. В.Г. Шухова, Белгород, Россия 2Индийский институт науки, Бангалор, Индия
	Аннотация: Работа посвящена исследованию влияния системы молниезащиты отдельно стоящего высотного сооружения на примере ветроэнергетической установки на возможность удара молнии в характерные зоны: стержень молниеприемника, защищаемый объект и поверхность земли. Целью работы является исследование атмосферных электрических полей для комплексного объекта: ветровая турбина, грозовое облако и поверхность земли, в двух вариантах – с применением системы молниезащиты и без нее. цель достигается за счет построения математической модели на основе уравнения Лапласа эллиптического типа и граничных условий для вычисления электрической функции потока. Проведены вычислительные эксперименты для установки в двух базовых конфигурациях – с применением системы молниезащиты стержневого типа и без неё, причем в конфигурации с применением системы молниезащиты варьировалась длина стержня молниеприемника. Наиболее существенным результатом работы является математическая модель с граничными условиями, которые определяют жесткие границы объекта исследования и константные значения на свободной и симметричной границах для вычисления электрической функции потока, методика нахождения удельного распределения линий функции потока вдоль границ исследуемой ветроэнергетической установки и результаты вычислительных экспериментов. В результате были найдены две зависимости в функции от длины стержня: первая - зависимость изменения количества эквипотенциальных линий функции потока, попадающих в стержень молниеприемника, и вторая - зависимость изменения количества эквипотенциальных линий эквипотенциальных линий функции потока, попадающих в ветровую турбину. Значимость полученных результатов состоит в определении эффективной длины молниеприемника для исследуемого объекта. Предложенная математическая модель позволяет анализировать атмосферное электростатическое поле вблизи остроконечных пиков и на поверхности земли рядом с защищаемым объектом, а также может быть применена для моделирования других объектов, использующих стержневой тип молниезащиты.
	Ключевые слова: математическое моделирование, молниезащита, ветроэнергетическая установка, вычислительный эксперимент, силовые линии поля, электростатическое поле. DOI: https://doi.org/10.5281/zenodo.3367052