Численное моделирование распределения электромагнитного поля в волноводе на основе фотонного кристалла
Выполнение численного моделирования двумерного волнового уравнения с использованием явной трехслойной схемы типа "крест". Распределение компоненты электромагнитного поля в волноводе при распространении по нему волны с различной несущей частотой.
Рубрика | Коммуникации, связь, цифровые приборы и радиоэлектроника |
Вид | статья |
Язык | русский |
Дата добавления | 04.11.2018 |
Размер файла | 454,1 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Московский государственный университет им. М.В. Ломоносова
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ РАСПРЕДЕЛЕНИЯ ЭЛЕКТРОМАГНИТНОГО ПОЛЯ В ВОЛНОВОДЕ НА ОСНОВЕ ФОТОННОГО КРИСТАЛЛА
А.С. Логгинов
А.Ш. Майоров
Методом численного моделирования рассчитаны распределения электромагнитных полей в волноводе на основе фотонного кристалла. Двумерное уравнение Даламбера, описывающее процесс распространения волны в интегрально-оптическом устройстве, решено методом конечных разностей с помощью явной численной схемы типа "крест". Показано, что волновод на основе фотонного кристалла обладает частотной селективностью. Это явление может найти применение при создании высокоскоростных систем связи со спектральным уплотнением каналов. Волноводы на основе фотонных кристаллов могут быть использованы в качестве базовых элементов нового поколения интегрально-оптических схем.
Стремительное развитие технологии изготовления устройств интегральной оптики и увеличение производительности персональных компьютеров открывают возможности для численного моделирования и экспериментальной реализации приборов, использующих фотонные кристаллы (ФК) [1]. ФК называют строго периодичную структуру, период которой сравним с длиной волны падающего на нее оптического излучения. При взаимодействии света с ФК возникает ряд кардинально новых явлений.
Многократная интерференция электромагнитных волн, дифрагировавших на каждой элементарной ячейке кристалла, приводит к образованию диапазонов частот (фотонных запрещенных зон), для которых распространение излучения по ФК невозможно. На этих частотах свет испытывает идеальное отражение от поверхности кристалла вследствие деструктивной интерференции, в то время как на других частотах он проходит через ФК. Это явление может быть использовано для решения широкого спектра физических задач, связанных с локализацией и эмиссией света. Материалы с фотонными запрещенными зонами позволяют создавать лазеры с микрорезонаторами, высокоэффективные волноводы, высокоскоростные оптические переключатели [2]. ФК обладают и многими другими интересными свойствами, связанными с дисперсией, оптической анизотропией и поляризационными характеристиками излучения. Такие свойства открывают возможность создания эффективных дисперсионных компенсаторов, преобразователей частоты, современных суперпризм, оптических поляризаторов, фильтров и т. д. [3-5]
Целью настоящей работы является численное моделирование распространения электромагнитной волны и расчет характеристик волновода на основе фотонного кристалла (рис.1). В нем создан волноведущий канал, ширина которого равна . Основная задача исследования -- изучение селективности подобного устройства. ФК, период которого обозначим символом , образован элементарными ячейками в форме квадрата со стороной .
Рис. 1, 2. Волновод на основе ФК.
Методы
Для определения характеристик изучаемого волновода для случая распространения ограничимся рассмотрением только одной компоненты поля --, для которой справедливо скалярное уравнение Даламбера [6]:
В работах [7] и [8], посвященных численному моделированию интегрально-оптических устройств, при вычислениях рассматривают уравнения Максвелла и для их решения используют алгоритм Кейна Йе. Непосредственное сравнение этой методики с результатами численного решения уравнения (1) показало, что его использование позволяет сократить время вычислений примерно на порядок. Этот выбор весьма важен, так как в рассматриваемой задаче объем вычислений значителен.
Явная численная схема типа "крест" [9] для двумерного волнового уравнения (1) может быть изображена в виде шаблона (рис.3). Сделаем некоторые пояснения. Пространство и компонента напряженности электрического поля дискретизированы так, что пространственная точка представлена как , а -- , где , и -- целые числа, и -- пространственные приращения, -- временное приращение. Пространственные и временные производные второго порядка в (1) заменяются так:
,
,
Если производные вычислять таким образом, тогда погрешность аппроксимации численной схемы порядка , и. Если учесть (2)-(4), то конечно-разностная форма (1) примет вид:
.
Рис. 3. Шаблон для явной трехслойной численной схемы типа “крест”.
Используя (5), найдем решение на шаге по времени. Оно вычисляется по (6), для этого необходимо знать распределение в двух предыдущих шагах времени.
где характеризует точку пространства, соответствующую текущим значениям , с диэлектрической проницаемостью и магнитной проницаемостью , а -- скорость света в вакууме.
Существуют определенные практические соображения при дискретизации пространства координат. Для того чтобы длина волны была представлена корректно, пространственные интервалы должны быть малы по сравнению с ней. При выборе временного шага дискретизации необходимо учитывать обобщенное условие устойчивости Куранта [9]:
Уравнение (6) необходимо дополнить граничными и начальными условиями. В задаче анализа электромагнитных полей в волноводе на основе ФК исследуемое пространство содержит поля, уходящие в бесконечность. Поскольку вычисления можно производить лишь в конечной области, рассматриваемое пространство необходимо ограничить, и на этих границах волны должны быть поглощены без возникновения существенных отражений [8]. электромагнитный поле волновод частота
Выбором начальных условий в систему вводится входной сигнал (гауссов профиль по оси , модулированный гармонической функцией с несущей частотой оптического сигнала). Математический образ входного сигнала можно представить в виде (8):
где - амплитуда входного сигнала, -- несущая частота, значения , и задают форму и пространственное положение функции Гаусса.
Результаты
Численное моделирование двумерного волнового уравнения (1) выполнено с использованием явной трехслойной схемы типа "крест" [9]. При расчете предполагались следующие параметры системы (рис.1, 2): показатель преломления материала ФК ; показатель преломления окружающей среды ; ширина волноведущего канала мкм; период ФК мкм, размер элементарной ячейки мкм. На рис. 4-11 показано распределение компоненты электромагнитного поля в устройстве (рис.1, ) при возбуждении волновода сигналом (8) различной несущей длины волны. ФК при определенных значениях длины волны является "зеркалами", обеспечивающими волноводный режим распространения света, отражающегося от них. при волноводном распространении содержит "боковые лепестки", образующиеся при интерференции света. Рис.4-12 демонстрируют спектральную селективность волновода на основе ФК. Это явление может найти применение при создании высокоскоростных систем связи со спектральным уплотнением каналов (WDM).
Рис. 4-11. Распределение компоненты электромагнитного поля в волноводе при распространении по нему электромагнитной волны с различной несущей частотой.
Рис. 12. Зависимость отношения амплитуды электромагнитного поля на выходе волновода на основе ФК к амплитуде на входе от длины волны входного сигнала.
Заключение
Численное моделирование позволяет определить распределение компоненты электромагнитного поля для предложенной конфигурации волновода на основе фотонного кристалла. Эти результаты являются основополагающими при проектировании (выборе требуемых параметров) и дальнейшей реализации устройств интегральной оптики. Компактность, высокое быстродействие и возможность интеграции ФК волноводов с ФК направленными ответвителями и резонаторами на одной подложке позволит использовать эти устройства в новом поколении интегральных схем.
Литература
[1] Joannopoulos J.D., Meade R.D., Winn J.N. Photonic Crystals: Molding the Flow of Light, Princeton: Princeton University Press, 1995. Pp. 45-50.
[2] Baba T. // IEEE J. Select. Topics Quantum Electron. 1997. Vol. 3. № 3. Pp. 808-830.
[3] Kosaka H., Kawashima T., Tomita A. et al. // J. Lightwave Technol. 1999. Vol. 17. № 11. Pp. 2032-2038.
[4] Steeland M.J., Osgood R.M. // J. Lightwave Technol. 2001. Vol. 19. № 4. Pp. 495-503.
[5] Imada M., Noda S., Chutinan A. et al // J. Lightwave Technol. 2002. Vol. 20. № 5. Pp. 873-878.
[6] Логгинов А.С., Майоров А.Ш., Кряжимский С.А. // Радиотехника. 2005. №1. Cтр. 24-29.
[7] Hagness S.C., Rafizadeh D., Ho S. T., and Taflove A. // J. Lightwave Technol. 1997. Vol. 15. № 11. Pp. 2154-2165.
[8] Chu S.T. and Chaudhuri S.K. // J. Lightwave Technol. 1989. Vol. 7. № 12. Pp. 2033-2038.
[9] Калиткин Н.Н. Численные методы. М.: Наука, 1978. 425-439 с.
Размещено на Allbest.ru
Подобные документы
Определение комплексных амплитуд составляющих вектора; диапазон частот. Расчет и построение графиков зависимостей поля от координат x, y, z. Вычисление среднего за период потока энергии через поперечное сечение волновода. Коэффициент затухания волны.
курсовая работа [831,3 K], добавлен 15.04.2014Структура электромагнитного поля основной волны. Распространение электромагнитных волн в полом прямоугольном металлическом волноводе. Резонансная частота колебаний. Влияние параметров реальных сред на процесс распространения электромагнитных волн.
лабораторная работа [710,2 K], добавлен 29.06.2012Основные понятия и классификация приборов для измерения напряженности электромагнитного поля и помех. Измерение напряженности электромагнитного поля. Метод эталонной антенны. Метод сравнения. Измерительные приемники и измерители напряженности поля.
реферат [31,8 K], добавлен 23.01.2009Уравнения Максвелла для анизотропной среды. Магнитная и электрическая проницаемость вещества. Представление решения системы уравнений в виде плоских волн. Анализ составляющих частей волновода. Уравнения непрерывности электрического и магнитного полей.
курсовая работа [218,7 K], добавлен 17.11.2010Исследование влияния электромагнитного поля на подземную антенну, расположенную на определенной глубине. Расчеты напряжения нагрузки проволочной антенны. Разработка программного продукта, позволяющего выполнять основные операции разработанного алгоритма.
дипломная работа [1,7 M], добавлен 07.06.2012Сигналы в системах (зондирующий, сигнал подсвета, запросный, собственное радиоизлучение объекта наблюдения, отраженный сигнал и т.п.). Электромагнитные поля. Поляризационная структура электромагнитного поля. Амплитудное равномерное распределение поля.
реферат [2,0 M], добавлен 14.12.2008Расчет компонентов и разработка вычислительного блока системы электромагнитного позиционирования. Обоснование выбора катушек индуктивности и изучение их влияния на стабильность системы. Измерение индукции электромагнитного поля при парной работе катушек.
дипломная работа [2,4 M], добавлен 16.07.2013Принципы работы детектора поля RD-14. Расположение закладного устройства в незаметном месте. Частота и мощность входного сигнала. Уровень и частота принимаемого сигнала. Интегральный метод измерения уровня электромагнитного поля в точке его расположения.
лабораторная работа [593,8 K], добавлен 15.03.2015Расчёт напряжённости электрического поля на входе радиоприёмного устройства при заданной мощности излучения. Определение скорости распространения и направления прихода электромагнитного поля. Изучение поляризационных характеристик и искажений сигнала.
курсовая работа [198,7 K], добавлен 23.12.2012Разработка и выбор функциональной схемы датчика электромагнитного расходомера. Формирование и исследование аналоговой, цифровой схемы. Расчет блока питания устройства. Порядок разработки алгоритма работы и программного обеспечения микроконтроллера.
курсовая работа [1,9 M], добавлен 19.08.2012