Вычислительный эксперимент с линзой Френеля Текст научной статьи по специальности «Компьютерные и информационные науки»

Н.Л. Казанский

ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ С ЛИНЗОЙ ФРЕНЕЛЯ

Для создания элементов плоской оптики (ЭПО) характерны [1-4]:

- расчет необходимой функции пропускания ЭПО в рамках геометрооптической модели;

- дискретизация ЭПО при расчете полученной функции пропускания на ЭВМ;

- сложный технологический процесс изготовления, обусловливающий квантование фазовой составляющей функции пропускания ЭПО;

#

- необходимость моделирования функционирования ЭПО в рамках дифракционной модели.

Для оценки качества работы ЭПО на этапе проектирования актуальным является вопрос исследования и компенсации дифракционных эффектов, связанных с технологией расчета и изготовления ЭПО• Такие факторы, как большой объем расчетных данных, значительное число вариантов значений параметров исследуемого световог< поля, наличие набора критериев качества, требования естественности восприятия результатов расчета выводят решаемую задачу из рамок расчетных, позволяя охарактеризовать ее как вычислительный эксперимент [3, 4]. Специфика вычислительного эксперимента как диалога с исследователем диктует требования к повышению

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

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

1. Дифракционный расчет при аксиальной симметрии

На аксиально-симметричную оптическую систему (рис. 1) с функцией амплитуд-

но-фазового пропускания Т(г) падает световой пучок с комплексной амплитудой

Е(г)=А(г)ехр(хф(г)), где А(г) и ф(г) - соответственно амплитуда и фаза падающей

волны. С учетом технологии расчета и изготовления функция амплитудно-фазового

пропускания ЭПО с аксиальной симметрией имеет вид [4, 5]:

N Г-Г 1

T(r)= 2 t • rect(-g-E + , r<a, (1)

p=i Р р

где N - число колец разрешения ЭПО? а - радиус ЭПО (rN=a) ; ôp - ширина, а -

внешний радиус р-го кольца с амплитудно-фазовым пропусканием t .

Р

Считая, что ô=minô »X, воспользуемся для оценки комплексной амплитуды

р Р /-

w (р, z) в точке наблюдения (р, z) , где р=\/х2+у2 (см. рис. 1), интегралом Кирхгофа. Комплексная амплитуда светового поля в точке наблюдения согласно (1) представляется в виде суперпозиции световых полей, созданных отдельными коль-

= 1

цами разрешения (г , г ), р=1,N

p-i р

N

w(p, z)= Е t • А (г - 6/2) • [w (р, z)-w (р, z)], (2)

р Р Р Р Р~1

где w (р, z) - комплексная амплитуда поля от круглого отверстия радиуса г . При

Р Ьг

условии 2a/f«l комплексную амплитуду светового поля w^(p, z) вблизи базовой точки (р , z) можно вычислить следующим образом [б, 4].

Е(г)

(Р, Z)

Рис. 1. Оптическая схема работы плоского оптического

элемента

кг? Р3_Ро

V5' 2) = 1-ь % ' ехр11к<1"ог + тгт»!-

Р ог ог

ехр(- зр) + 102(чр, 5р)] при |др/Бр|<1,

Б2 Д-Ч

-±ехр(1^£ )-ехр(- зр)-1У0(Чр, вр)] при |др/зр1>1,'

Р

1 1 кг • р

где Ь = др=*г»(±- - в ' ип и Уп - функции Ломме-

°* * * о о г ог

ля [б], - точка фокусировки (2а/Ь0«1) падающего на ЭПО пучка Е(г).

Описанный способ вычисления дифракционного интеграла базируется на быстр методе расчета функций Ломмеля, предложенном в [2].

2. Программный комплекс для вычислительного эксперимента

Программный комплекс для вычислительного эксперимента создан на базе пакета прикладных программ обработки изображений и цифровой голографии ППП [7] и включает реализованные на языках ФОРТРАН и ПЛ/1 алгоритмы предлагаемых методов вычисления дифракционного интеграла, а также модуль связи ППП с комплексом графических программ ГРАФОР [8].

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

Выходными данными являются массивы значений интенсивности и фазы светового поля, записываемые в базу данных ППП.

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

3, Анализ возможностей вычислительного эксперимента

Созданные на основе описанных методов программные средства позволяют легко воспроизвести результаты работы [б] для случая дифракции сходящейся сферическо волны на круглом отверстии, выявив при этом асимметрию распределения относительно фокальной плоскости [9, 10], и наглядно представить результат в виде трехмерного изображения (рис. 2). На рис. 2 по вертикальной оси отображена интенсивность, а в горизонтальной плоскости выбраны безразмерные координаты

и ), которые изменяются в пределах: 0<р<4тх# ^Влйё^Ьк.

О 2

В то же время использование при расчетах метода, основанного на Функциях Ломмеля, не позволяет выявить эффект фокального сдвига максимума интенсивности по направлению к линзе при малых числах Фпенеля Г9-121.

ПОДЕ ОТ ЛИНЗЫ

Рис. 2. Пространственное распределение интенсивности в

фокальной области линзы (0<р<4тх, -8и<2<8п)

Исследованию точности используемого приближения Френеля посвящены работы [13, 14], в которых показано, что при достаточно больших £-числах (£/2а>12) и к-Ь »1 результаты для модуля амплитуды отличаются от точного значения не бо-

О 2

лее чем на 2%. Для обеспечения более точных расчетов необходимо включение в программное обеспечение вычислительного эксперимента методов, основанных на не посредственном вычислении интеграла Кирхгофа [9, 10]. Однако такие расчеты тре буют существенного увеличения машинного времени и применения более мощной вычислительной техники.

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

Заключение

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

шах

р/а

р/а

Рис. 3. Дифракция плоской волны на круглом отверстии -графики распределений интенсивности для различных чисел Френеля Е

а

- Зтх

-2тх

-ti

Рис. 4. Форма фокального пятна линзы Френеля

Литература

1. Сисакян И.Н.,Сойфер В.А. // Физические основы и прикладные вопросы голографии. Л.: ЛИЯФ, 1984. С, 142.

2. Васин А.Г., Голуб М.А., Данилов В.А., Казанский Н.Л., К а р п е е в С.В, Сисакян И.Н., С о й -ф е р В.А., Уваров Г.В. М.: Препринт ФИАН СССР, 1983, » 304.

3. С о й ф е р В.А. // Современные проблемы физики и ее приложений. Ч. I. М.: ВИНИТИ, 1987. С. 5.

4. Голуб М.А., Казанский Н.Л., Сисакян И.Н., С о й ф е р В.А. Автометрия, 1987, № 5.

5. Бобров С.Т., Г р е й с у х Г.И., Туркевич Ю.Г. Оптика дифракционных элементов и систем. Л.: Машиностроение, 1986. С. 200.

6. Борн М., Вольф Э. Основы оптики. М.: Наука, 1973. С. 397.

7. Пакет прикладных программ обработки изображений и цифровой голографии. Государственный фонд алгоритмов и программ. Регистрационный № П004582 от 1.09.1980.

8. Баяковский Ю.М., Галактионов A.B., Михайлова Т.П. Графор. Графическое расширение Фортрана. М.: Наука, 1985.

9. К и н б е р Б.Е., Новоселов C.B. РЭ, ¿0, 1469 (1985)

10. L i Y., Wolf Е. JOSA, AI, 801 (1984).

11. Дементьев A.C., Домаркене Д.П. Опт. и спектр 56, 869 (1984).

12. L i Y-, H o r n g J.S. Optik, 127 (1985).

13. S o u t h w e 1 1 W.H. JOSA, 71, 7 (1981).

14. G o o d man D. JOSA, 72, 1560 (1982).

15. K a t h u r i a Y.P. Optica Applicata, JL4, 509

16. B o t t e m a M. JOSA, 59, 1632 (1969) .