Содержание к диссертации
Введение
1. Состояние исследований о влиянии струи электрореактивного двигателя на радиосистемы космического аппарата 18
2. Дифракция электромагнитных волн на двух сферах 36
2.1. Сферические функции, основные характеристики электромагнитного поля 37
2.2.Метод преобразований амплитуд волн 46
2.2.1. Смещение начала координат 46
2.2.2. Вращение системы координат 52
2.3.Строгое решение задачи дифракции волн на двух сферах 55
2.3.1. Поле произвольно ориентированного диполя 55
2.3.2. Возбуждение двух сфер турникетной антенной 65
2.3.3. Дифракция на сфере, расположенной над проводящей плоскостью 67
2 4 Метод переотражений при решении задач дифракции на двух телах 71
2.5.Дифракция волн на неоднородном шаре без центральной симметрии 74
2.6.Метод зеркальных изображений в сферической системе координат 77
2.7.Метод расчета антенной решетки на сфере 78
3. Дифракция волн на N телах сложной формы 78
3.1. Дифракция волн на двух телах: диске и сфере 78
3.2.Дифракции волы на N телах 89
3.3.Дифракция волн на плазменном образовании произвольного размера 98
4. Плоскослоистое приближение при заданной точности расчета 109
4.1 . Метод самосогласованных конечных разностей 114
4.2.Принцип предельного перехода, особенность в нуле є 140
4.3 .Исследования полей вблизи нуля є 147
4.4.Взаимодействие волны с холодной плазмой 156
4.5.Обратная задача рассеяния волн 169
Приложения 1,2 177
5. Проектирование антенн космического аппарата 180
5 Л . Моделирование антенных систем 183
5.2.Метод расчета дифракции волн на двух телах сложной формы 185
5.3.Расчетные модели и результаты 189
Приложение, Особенность решения задачи дифракции волн на бесконечно тонком идеально проводящем диске 198
3аключение 210
Литература 213
- Состояние исследований о влиянии струи электрореактивного двигателя на радиосистемы космического аппарата
- Сферические функции, основные характеристики электромагнитного поля
- Дифракция волн на двух телах: диске и сфере
- Метод самосогласованных конечных разностей
- Моделирование антенных систем
Введение к работе
Актуальность проблемы. Развитие космической техники требует широкого применения математического моделирования при создании космических систем, что обусловлено их существенным усложнением, удорожанием и необходимостью сокращения времени проектирования. Создание антенных систем в СВЧ диапазоне радиоволн на космическом аппарате (КА) связано с жесткими ограничениями по месту расположения, весовым, габаритным и другим характеристикам, что сопряжено со значительными техническими трудностями. Положение усугубляется при размещении антенн вблизи острых выступающих частей или при наличии около КА неоднородной холодной плазмы с критической концентрацией. Такая плазма может быть, например, от струи электрореактивного двигателя (ЭРД) или появляться при возрастании активности Солнца. В обоих случаях при анализе работы антенн надо учитывать особенности решения задач дифракции электромагнитных волн на острой кромке и распространения волн в неоднородной среде. Все это требуетразвития строгой теории дифракции волн на нескольких телах.
Особенности решения задачи взаимодействия электромагнитного излучения с холодной плазмой при резком изменении свойств среды приводят к усложнению проблемы электромагнитной совместимости радиосистем в присутствии вблизи КА холодной плазмы, концентрация которой достигает критической. Роль математического моделирования в этом случае возрастает из-за невозможности экспериментальной наземной отработки всех аспектов нелинейного взаимодействия холодной плазмы около КА с СВЧ полем. В результате появляется необходимость моделирования также космических радиофизических экспериментов. Возросшие требования к безопасности космических полетов подчеркивают значимость предполагаемых исследований.
Теория дифракции электромагнитных волн на двух телах систематически изложена Ивановым Е.А. Метод решения для двух шаров в сферических координатах применили при многочисленных расчетных исследованиях Bruning J.H. и Lo Y.T. Однако решение и расчет задач дифракции волн на сложных препятствиях, которыми можно моделировать КА, антенны с диском, неоднородную плазму, представляют существенные трудности. Аналитические исследования дифракции волн на диэлектрическом шаре с шаровым включением проводились Uzunoglu N.K. с помощью векторной теоремы сложения. Анализ особенности решения такой задачи при резком изменении диэлектрической проницаемости можно провести в плоскослоистом приближении, развитом в работах Бреховских Л.М., Каценеленбаума Б.З., которое позволяет сократить число независимых переменных задачи.
Плоскослоистое приближение актуально в связи с воздействием мощных сверхкоротких лазерных и СВЧ импульсов на мишени, из-за которых возникает облако холодной плазмы с закритической концентрацией. В этом направлении требуется развитие работ Гинзбурга В.Л. и Ландау Л.Д. Сразу отметим, что рассматриваемая задача сводится к решению уравнения Шредингера, имеющего приложения в различных областях физики.
Основная цель исследований. Работа направлена на создание машинных методов проектирования антенн КА с учетом электромагнитной совместимости; на исследования взаимодействия электромагнитных волн с холодной плазмой вблизи критической концентрации. Предполагается развитие теории дифракции волн в неоднородных, в том числе локально неоднородных, средах.
Научная новизна работы состоит в следующем: - создан метод решения и расчета задачи дифракции волн на нескольких телах вращения сложной формы размером до нескольких длин волн; решены задачи дифракции волн на двух телах типа диск с полусферами и неоднородный шар без центральной симметрии; выявлена аномалия в решении одномерного волнового уравнения; теоретически предсказана качественная зависимость решения задачи нормального падения плоской волны на холодную плоскослоистую плазму вблизи критической концентрации от малых параметров задачи; разработано плоскослоистое приближение при заданной точности задачи взаимодействия электромагнитных волн с холодной плазмой; - решена задача дифракции волн на плазменном образовании произвольного размера с учетом особенности в нуле диэлектрической проницаемости; - разработаны математические методы моделирования антенн КА с учетом взаимодействия радиоволн с холодной плазмой;
Практическая значимость работы подтверждается тем, что разработанные методы использованы при создании антенн КА РКК «Энергия», прошедших летные испытания. Развитые методы могут быть использованы при проектировании антенных систем КА с учетом влияния локализованного вблизи КА источника холодной плазмы. Численное моделирование с использованием развитых методов позволяет значительно улучшить качество проектов при сокращении времени и стоимости разработки.
Положения, выносимые на защиту:
Метод преобразования амплитуд волн при смещении и вращении сферической системы координат. Метод позволяет решать и проводить вычисления по итерационной схеме задач дифракции волн на двух телах сложной формы характерного размера до нескольких длин волн.
Метод решения задачи о распространении волн в нелинейной слабо поглощающей среде в плоскослоистом приближении с заданной точностью. Метод сокращает число независимых переменных задачи. Вывод количественных критериев коротковолнового и длинноволнового приближений, критерия появления поверхностной волны.
Выявлена аномалия в решении одномерного волнового уравнения при резком изменении свойств среды.
Метод предварительного проектирования антенн космических аппаратов при наличии неоднородного плазменного образования и взаимодействия антенн между собой и корпусом космического аппарата.
Задача предварительного проектирования антенн КА состоит в определении типа и параметров антенной системы по заданной диаграмме направленности (ДЫ), которая, в свою очередь, определяется из программы полета КА. Предварительное проектирование предполагает не просто расчет параметров антенн, удовлетворяющих заданным требованиям. Предполагается моделирование работы радиосистемы в процессе выполнения программы полета КА с учетом электромагнитной совместимости. Работоспособность радиосистемы КА обеспечивается учетом влияния других СВЧ систем, учетом влияния радиофизической обстановки вблизи КА, особенно во время проведения радиофизических экспериментов, связанной с наличием плазменных образований. Проектирование предполагает, в том числе, анализ влияния мощного СВЧ излучения передатчика на устойчивость плазменной струи ЭРД.
Обратная задача определения параметров антенн нами решается методом подбора - прежде всего, методом простого моделирования, использующего программные средства решения прямой задачи расчета ДН антенн (по заданным синусоидальным токам на вибраторах). Существенным является присутствие на линии радиосвязи с КА или вблизи КА плазменных областей (искусственного или естественного происхождения) с критической концентрацией, соответствующей нулю диэлектрической проницаемости в модели сплошной среды (для рабочей длины радиоволны).
При предварительном проектировании антенн КА необходимо, прежде всего, проведение качественного анализа явлений, влияющих на радиофизические характеристики антенн. Из основных явлений выделим дифракцию электромагнитных волн на антеннах, корпусе КА с его острыми выступающими частями, плазменном образовании и нелинейное взаимодействие поля передающих антенн с неоднородной холодной плазмой, например, струи ЭРД. Особое внимание надо обратить на проблему электромагнитной совместимости. Изучение перечисленных выше явлений приводит к необходимости решения задачи дифракции на двух телах. На этапе предварительного проектирования можно ограничиться простыми моделями, описываемыми полностью или частично сферическими координатами. Такие модели позволяют выявить основные закономерности и учесть особенности решения краевых задач. При этом точность расчетов не нужна выше (7-10) % пополю.
Данная монография посвящена разработке методов предварительного проектирования антенных систем КА на основе математического моделирования, предполагающего использование строгой теории дифракции на двух телах. При этом существенное внимание уделяется особенностям решения краевых задач.
Задача дифракции волн на диске со сферой является модельной задачей при расчете турникетной антенны с диском (радиуса до длины волны) около КА. На этапе предварительного проектирования антенн, когда не требуется высокая точность расчетов, рассматривается осесимметричная система двух тел. Космический аппарат довольно сложной формы, с точки зрения выявления общих закономерностей, можно аппроксимировать идеально проводящим диском, сопряженным с двумя полусферами. Наличие кромки диска важно при предварительном анализе влияния, например, солнечных батарей. Плазменное образование вблизи КА моделируется неоднородным шаром, образованным вложенными друг в друга без центральной симметрии однородными шарами. Такая модель холодной плазмы вблизи критической концентрации позволяет учесть трансформацию волн ортогональных поляризаций и нелинейные явления, из-за чего СВЧ поле передатчика может попасть в приемную антенну.
Машинные методы проектирования СВЧ систем находят все большее применение в науке и технике. Так, например, стоящие проблемы, принципы создания и перспективы применения объемных интегральных схем СВЧ и КВЧ изложены в книге [1]. Тем более это касается - космической техники, стоимость создания которой возрастает. Математическое моделирование работы радиосистем космического аппарата (КА) позволяет существенно улучшить качество проектов, осуществить системный подход к созданию космического комплекса. Важной проблемой является электромагнитная совместимость радиосистем, когда на борту КА находится источник неоднородной низкотемпературной плазмы. В работах этого направления автор имел поддержку академика РАН Б,Е. Чертока, профессоров Б.З. Каценеленбаума, Е.И. Нефедова, Ю.А. Ильина, к.т.н. Н.А. Яблочкина.
В настоящее время главные усилия в теории дифракции на двух или нескольких телах (когда эффективные размеры тел и расстояния между ними порядка длины волны) направлены на развитие численных методов. Вычислительные методы в электродинамике обобщены при блестящем изложении в книге под редакцией Р. Митры [2]. Примером применения интегральных методов, частично совмещенных с аналитическими методами, являются работы Е.Н. Васильева [3]. Аналитические исследования решения краевых задач с острыми краями проведены В.Ф. Апельциным и А.Г. Кюркчаном [4]. Среди наиболее известных методов решения задачи рассеяния волн группой тел отметим такие, как метод специальной ортогонализации Я.Н. Фельда, метод вспомогательных токов [5], метод дискретных источников [6], метод диаграммных уравнений [7], метод В. Тверского [8]. Аналитический метод разделения переменных систематически изложен Е.А. Ивановым [9].
В предлагаемой работе применяется метод собственных функций. Этот метод в сферических координатах позволяет обобщить хорошо известные решения для одиночных препятствий на решения для нескольких тел. Кроме этого, сферические функции хорошо исследованы, что делает вычисления не слишком громоздкими. Осевая симметрия задачи за счет выделения одной азимутальной гармоники упрощает решение и его анализ. Разделение переменных в сферической системе координат с использованием векторной теоремы сложения позволяет расширить класс решаемых задач дифракции электромагнитных волн от уединенного шара до произвольной совокупности непересекающихся шаров [10 - 12]. Однако в ряде работ, основанных на использовании скалярных потенциалов Дебая и теоремы сложения, например [9, 13], принцип суперпозиций полей (представленньїх разложением по собственным функциям в локальных системах координат со смещенными центрами) ошибочно переносился на суперпозицию скалярных потенциалов - векторная задача заменялась скалярной. Нами в [14, 15] уточнено это решение и дан метод решения задачи дифракции на непересекающихся шарах. Получены преобразования амплитуд волн при переходе между двумя локальными сферическими системами координат, полученными смещением начала сферической системы координат. Преобразования показывают, что возбуждение одного типа волн в первой локальной сферической системе координат при переходе смещением ко второй приводит к появлению двух типов волн.
Аналитический метод решения задачи дифракции волн на двух телах дает возможность разработать наглядные модели, позволяющие изучать физику взаимодействия двух отражателей. Сначала для поля вблизи каждого тела выделены неизвестные - амплитуды волн, обусловленные влиянием соседнего отражателя, а затем показана взаимосвязь полей около этих отражателей через преобразования неизвестных амплитуд волн при переходе между двумя локальными системами координат. Из граничных условий на поверхностях отражателей в локальных координатах получена неполная бесконечная система линейных алгебраических уравнений, которая дополняется преобразованиями неизвестных амплитуд волн до полной системы. Доказана разрешимость этой бесконечной системы уравнений методом усечения.
Важным свойством преобразований амплитуд волн, если известно поле вблизи одного из взаимодействующих тел (при наличии второго тела), является то, что эти преобразования допускают многократное применение (при этом амплитуды волн не изменяются). Поэтому приближенное решение задачи дифракции на двух телах может быть уточнено многократным применением преобразований амплитуд волн. В результате разработан итерационный метод расчета при введении малого параметра [16, 17], когда, например, при многократном применении преобразований амплитуд волн постепенно увеличивается размер одного из препятствий от минимального (для которого итерационный процесс сходится) до заданного. Развитый таким образом метод преобразований амплитуд волн позволяет проводить вычисления для тел сложной формы размеров до нескольких длин волн.
Из литературных источников известны многочисленные расчеты дифракции на двух шарах [11]. Известны расчеты для ансамбля шаров, произвольно расположенных в пространстве [18], причем используется итерационный метод расчета, близкий к нашему методу, но опубликованный значительно позже. При этом нам удалось использовать в расчетах тела существенно более сложной формы.
Разработанный метод преобразований амплитуд волн использован для решения и расчета задач дифракции волн на двух шарах; на неоднородном диэлектрическом шаре, образованном из вложенных друг в друга однородных шаров без центральной симметрии (шары не пересекаются). Хотя ранее были известны решения этих задач (см. библиографию), на них продемонстрирован предложенный метод решения, Кроме этого даны решения новых задач дифракции, которые значительно сложнее и представляют интерес в практике моделирования антенн. Так решена задача дифращии волн на двух телах типа идеально проводящий диск с полусферами и неоднородный шар, образованный вложенными друг в друга (при осевой симметрии задачи) однородными диэлектрическими шарами. Получено решение задачи дифракции на системе N тел. Дано решение задачи дифракции волн на неоднородном диэлектрическом теле произвольного размера, когда решение по методу геометрической оптики "сшивается" со строгим решением для неоднородного шара вблизи нуля диэлектрической проницаемости (где необходимо строгое решение) с помощью принципа Гюйгенса.
Выделена особенность решения задачи дифракции волн на неоднородном шаре, образованном из вложенных один в другой однородных шаров без центральной симметрии, в точке «почти касания» шаров, когда свойства среды резко изменяются. Проведены исследования этой особенности в плоскослоистом приближении, позволяющем улучшить сходимость решения.
Исследованы особенности решения задачи дифракции на краю диска. Выделены особенности решения при осесимметричном возбуждении диска в квазиоптической области и при скользящем падении вертикально поляризованной плоской волны на диск. При предварительном проектировании рассматривается возбуждение диска радиуса меньше половины длины электромагнитной волны, поэтому указанную особенность решения при осесимметричном возбуждении диска можно не учитывать при расчетах. Численные исследования в [11] показали, что при расчетах дифракции волн на двух шарах хорошо «работает» лучевая оптика. Поскольку расстояние между диском и сферой, по крайней мере, больше радиуса диска, то это дает нам основание рассматривать диск в модели двух тел на уровне строгости метода частичных областей (при «сшивании» полей в частичных областях по методу наименьших квадратов), что затем было подтверждено расчетами.
Показано практическое применение теории дифракции волн на нескольких телах при проектировании антенн. Создан метод предварительного проектирования антенн КА на основе решения задачи дифракции волн на двух телах сложной формы. Метод позволяет проводить вычисления при размерах тел сложной формы до нескольких длин волн. Электродинамическая модель двух тел позволяет учесть взаимное влияние антенн между собой, с корпусом КА или его выступающими частями, в присутствии плазменной струи ЭРД. Модель позволяет исследовать явления, связанные с электромагнитной совместимостью радиосистем при наличии на борту КА источника плазмы.
Важным этапом предварительного проектирования являются исследования особенностей поля в точке «почти касания» шаров и на кромке диска, а также условий, в которых эти особенности сказываются.
Рассмотрено решение стационарной задачи взаимодействия плоской электромагнитной волны с неоднородной плазмой в плоскослоистом приближении - частного случая решения фундаментального уравнения Шредиыгера. Предложенный метод самосогласованных конечных разностей предполагает разбиение неоднородного слоя диэлектрика на однородные подслои переменной толщины, согласованной с изменением функции s(z) на слое. Автомодельное решение позволило сократить количество независимых параметров слоя. Разработанные количественные критерии применимости коротковолнового и длинноволнового приближений, критерий появления поверхностной волны позволяют проводить вычисления с заданной точностью.
В результате теоретически предсказана качественная зависимость решения задачи взаимодействия электромагнитной волны с холодной плазмой в плоскослоистом приближении вблизи критической концентрации от малых параметров задачи. Получено решение для плотности потока энергии, которое хорошо сходится и этим привлекательно для практики.
В окрестности нуля є, когда неприменимо плоскослоистое приближение, можно использовать модель плазмы в виде неоднородного шара, образованного из вложенных друг в друга однородных шаров без центральной симметрии, что приводит к трансформации полей, двух ортогональных поляризаций. Локальное применение плоскослоистого приближения при исследовании дифракции волн на таком образовании позволяет улучшить сходимость решения и одновременно учесть слабую нелинейность задачи.
Трансформация электромагнитных волн в плазменные волны может вызвать в струе плазменный резонанс. Исследования возникающих нестационарных колебаний позволят изучать специфические явления около КА и оценить их влияние, прежде всего, на радиосвязь, устойчивость и другие параметры струи ЭРД вблизи критической концентрации плазмы. Таким образом, появляется нелинейная, нестационарная задача возбуждения в неоднородной плазме плазменных волн. Плазменные волны, в свою очередь, могут привести опять к электромагнитным колебаниям или могут изменить форму струи, в результате оказать влияние на тягу двигателя.
Состояние исследований о влиянии струи электрореактивного двигателя на радиосистемы космического аппарата
Развитие космической техники требует широкого использования электрореактивных двигателей (ЭРД), как двигателей высокой удельной тяги. Однако из-за выявленной аномалии решения волнового уравнения в среде без поглощения, плазменная струя ЭРД может оказывать существенное влияние на работу радиотехнических систем космического аппарата (КА) в СВЧ диапазоне. Специфическое взаимодействие электромагнитных волн с холодной плазмой струи для радиоволн заданной частоты происходит вблизи нуля диэлектрической проницаемости (є) в модели сплошной среды. Это касается КА различного назначения от связных ИСЗ до аппаратов астрофизических исследований. В результате усугубляется проблема электромагнитной совместимости (ЭМС) антенн КА в присутствии плазменной струи ЭРД из-за нелинейного ее взаимодействия с электромагнитной волной и собственного СВЧ излучения ЭРД. Возможен сбой в работе радиосистем КА, вероятность сбоя резко возрастает при максимальной тяге ЭРД.
Математическое моделирование взаимодействия электромагнитных волн с неоднородной плазмой приводит к необходимости строгого решения задачи дифракции электромагнитных волн на двух телах. Плазменное образование ЭРД моделируется диэлектрическими шарами, вложенными один в другой со смещенными центрами. В этом случае в точке «почти касания» шаров (при резком изменении диэлектрической проницаемости) возможно нелинейное взаимодействие электромагнитной волны с плазмой.
Сферические функции, основные характеристики электромагнитного поля
Задача дифракции электромагнитных волн на двух телах, поверхности которых совпадают с координатными поверхностями сферической системы координат, решается с применением векторной теоремві сложения [11-12].
В настоящей книге рассматриваются преобразования амплитуд волн при переходе от одной локальной сферической системы координат к другой, начала которых разнесены, а полярные оси совпадают. Рассматриваются также преобразования амплитуд волн при вращении системы координат.
Задача дифракции на нескольких телах решается методом преобразований амплитуд волн следующим образом. С каждым телом связывается локальная система координат, электромагнитное поле в которой представляется амплитудами волн сторонних токов и неизвестными амплитудами волн, определяемыми присутствием других тел. Из граничных условий на поверхности каждого тела в локальной системе координат образуется неполная система бесконечных линейных уравнений относительно неизвестных амплитуд волн. Преобразования неизвестных амплитуд волн при переходе от одной локальной системы координат к каждой другой, дополняют эту систему уравнений до полной. Система может быть решена методом усечения. Решение ее будет строгим решением задачи дифракции на нескольких телах.
Сферические функции, основные характеристики электромагнитного поля
Рассматривается электродинамическая задача возбуждения произвольными сторонними токами системы из двух идеально проводящих тел: диска, сопряженного с двумя полусферами, и сферы. Решения для произвольной модуляции т (по ср), в силу осевой симметрии системы тел, независимы.
Задача возбуждения диска, сопряженного с двумя полусферами, в сферической системе координат сводится к решению бесконечной системы линейных алгебраических уравнений относительно амплитуд волн [85].
Дифракция волн на двух телах: диске и сфере
В результате особенное значение при исследовании полей вблизи точки є = 0 приобретает постановка задачи. Так если, например, в [32, 86] исследовано решение для полубесконечного слоя (є 1) с линейной зависимостью e(kz), то нами в разделе 4.1 рассматривается решение для линейного слоя при є 0, а затем принимается sN - 0.
В разделе 4.2 - 4.4 уже исследуется физическая задача о распространении плоской электромагнитной волны в плоском слое произвольной толщины с линейной зависимостью s(kz). Слой проходит через точку є = 0 при направлении падения волны, которое близко к нормали. Предполагается сверхмалое поглощение в среде (критерий дан). Слой в окрестности точки s = 0 заменяется нелинейным слоем с зависимостью s(kz), для которой имеется точное решение в элементарных функциях. Предельным переходом выделяется основное решение, соответствующее случаю нормального падения волны на полубесконечный слой без поглощения с линейной зависимостью e(kz). Показано, что нуль диэлектрической проницаемости (є) в неоднородной плоскослоистой среде без поглощения является математической критической точкой. Вблизи нуля є решение физической задачи качественно зависит от малых параметров, прежде всего, затухания (є"), отклонения угла падения волны (8) от нормального, кривизны (1/г) поверхности s=const, а также изменения вида функции с(г, 0, ф).
Моделирование антенных систем
Подготовительным этапом к проектированию является расчет линии радиосвязи "КА - Земля", "КА - КА", в результате которого определяется требуемая ДН антенны. После этого задача проектирования антенн сводится к выбору типа антенной системы, ее геометрии и положения на КА, определению принципа функционирования системы. Трудности проектирования антенной системы связаны с необходимостью моделирования антенн вместе с КА при наличии возмущающих поле антенны объектов, таких, например, как выступающие элементы конструкций, острые кромки, плазменные образования. При этом КА может находиться в неоднородной атмосфере планеты - вариант посадки КА. Существенным фактором является наличие на линии радиосвязи плазменных областей с критической для каждой длины волны излучения концентрацией, при которой диэлектрическая проницаемость плазмы обращается в нуль. Наличие приведенных выше особенностей требует отказа от моделей, основанных на геометрической оптике и связанных с лучевой геометрией распространения излучения.
Выбор физико-математической модели сильно зависит от наличия препятствий в так называемой резонансной области, когда характерный размер препятствия L приблизительно равен длине волны излучения X, и особенно при нахождении этого препятствия в ближней зоне антенны.
Математическое моделирование применяется на предварительном этапе проектирования антенн КА при использовании строгих решений задач дифракции на двух телах, когда требуемая точность расчета характеристик по полю составляет 7-10 %.
