Дифракция в решёточных оптических фильтрах с угловым селективным светопропусканием
Автор: Закируллин Рустам Сабирович
Журнал: Компьютерная оптика @computer-optics
Рубрика: Дифракционная оптика, оптические технологии
Статья в выпуске: 3 т.44, 2020 года.
Бесплатный доступ
Исследовано влияние дифракции на коэффициент направленного светопропускания решёточных оптических фильтров, предназначенных для смарт-окон. Фильтр имеет угловое селективное светопропускание за счёт двух тонкоплёночных решёток, образованных «непропускающими» (поглощающими, отражающими или рассеивающими) параллельными полосами на поверхностях окна с одинарным или двойным остеклением. Приведены методы расчёта оптимального угла наклона двух решёток фильтра, их взаимного расположения на разных поверхностях окна и ширин полос обеих решёток для минимизации светопропускания в конкретную дату и время суток при заданных значениях широты и долготы здания, азимута ориентации окна. Разработан метод расчёта дифракции на входной и выходной решётках фильтра при сложной траектории движения Солнца относительно окна и соответствующем изменении угла падения солнечных лучей на решётки. Для оценки уменьшения коэффициента светопропускания фильтра из-за влияния дифракции в расчётную формулу введён дополнительный сомножитель. Проведены расчёты геометрических параметров фильтров и дифракции. Значения дифракционного сомножителя составляют 0,9999992 и 0,9999998 для окон с одинарным или двойным остеклением соответственно, т. е. влиянием дифракции на светопропускание решёточного фильтра можно пренебречь.
Дифракция, дифрагированный свет, решётка, фильтры, тонкие плёнки, решёточный оптический фильтр, угловое селективное светопропускание, коэффициент направленного светопропускания
Короткий адрес: https://sciup.org/140249998
IDR: 140249998 | DOI: 10.18287/2412-6179-CO-654
Текст научной статьи Дифракция в решёточных оптических фильтрах с угловым селективным светопропусканием
Угловая толерантность оптических фильтров, т. е. обеспечение постоянства характеристик пропускания или отражения при изменении углов падения лучей, является одним из важнейших направлений их совершенствования. В фильтре с резонансной решёткой эту задачу предложено решать путём введения второй решётки, параллельной первой [1]. В таких же фильтрах при наклонном падении коллимированного пучка лучей угловую независимость предложено увеличивать путём применения многомодового плоского волновода и би-атомной решётки [2]. Описаны фильтры с низкой угловой чувствительностью на основе одномерного фотонного кристалла [3], а также на основе стека интегрированных двумерных фотонных кристаллов [4]. Разработан фильтр с оптическими метаповерхностями на гибкой подложке, обеспечивающий независимость отклика от поляризации и угла падения световых лучей [5].
В обзоре [6] описаны субволновые дифракционные решётки с четырьмя секторами с линиями решёток, ориентированными под разными углами, и показана зависимость интенсивности отражённого света от угла наклона секторальных линий. При исследовании резонансных фотонно-кристаллических структур с дифракционной решёткой [7] получены зависимости модуля коэффициента отражения от угла падения для различных сред. В статье [8] выведены зависимости угловой и спектральной селективности от углов падения опорного пучка при получении отражательных и пропускающих цветных голограмм.
Известны оптические фильтры, использующие угловую зависимость спектра пропускания. Описаны подобные резонансные фильтры с кремниевыми наномембранами на стеклянных подложках [9] и жидкокристаллические фильтры с резкими резонансными отражательными характеристиками в угловой зависимости [10]. В работах [11, 12] исследованы тонкие плёнки с наклонными столбчатыми микроструктурами, имеющие угловую селективность коэффициента пропускания в узких угловых и спектральных диапазонах.
Дифракция накладывает ограничения не только на микро- и наноструктуры, рассмотренные выше, но и, например, на муаровые решётки [13], период которых может достигать миллиметровых значений. При наклонном падении пучка параллельных световых лучей на такие решётки их период уменьшается [14], и это позволяет исследовать дифракцию рентгеновских лучей на «обычных» дифракционных решётках.
При увеличении угла падения параллельного пучка лучей дифракционные ограничения могут возникнуть и в решётках с миллиметровыми периодами и более. Работа посвящена исследованию влияния дифракции на коэффициент направленного пропускания решёточных оптических фильтров с угловым селективным светопропусканием, имеющих периоды решёток до десятков миллиметров и предназначенных для смарт-окон с одинарным и двойным остеклением.
1. Решёточный оптический фильтр
Угловая селективность фильтра обусловлена пропусканием пучка параллельных световых лучей через две поверхностные тонкослойные решётки [15–19]. Схемы фильтров для одинарного и двойного остекления представлены на рис. 1. Входная и выходная решётки формируются нанесёнными на поверхность стекла на постоянной основе поглощающими, отражающими или рассеивающими параллельными тонкослойными полосами, между которыми находятся направленно пропускающие полосы. Доля проходящего светового потока зависит от угла падения Θ , ширин полос с 1 ... с 4 , расстояния между решётками s и относительного расположения решёток, оцениваемого характеристическим углом фильтра Θ c (луч проходит через центры полос обеих решёток). В отличие от горизонтальных или вертикальных жалюзи, фильтр может быть адаптирован к траектории Cолнца относительно окна с учётом азимута его ориентации из-за возможности расположения полос под любым углом наклона. Способ регулирования направленного све-топропускания запатентован автором (патенты РФ № 2509324 и № 2677069).

Рис. 1. Схемы оптических фильтров для окон: с одинарным (а) и двойным (б) остеклением
В работе [15] рассчитаны, построены и экспериментально подтверждены зависимости коэффициента светопропускания от углов падения (угловые характеристики светопропускания фильтров) в диапазоне от 0° до 60° через каждые 5° для 7 фильтров с разными ширинами полос и разными периодами решёток, в т.ч. с кратными периодами двух решёток. Для равномерного регулирования светопропускания по всей площади окна периоды двух решёток должны быть одинаковыми или кратными.
Метод проектирования и расчёта решёточных оптических фильтров с угловым селективным свето-пропусканием представлен в [16], где показано, что при изменении угла падения также меняется доля излучения, которая проходит через обе решётки. Там же приведён расчёт геометрических и оптических параметров фильтра для получения заданной угловой зависимости его светопропускания, т. е. для предварительной адаптации характеристики фильтра к углам падения, изменяющимся при перемещении источника света относительно фильтра.
Метод оптимизации угловой селективной фильтрации прямой солнечной радиации путём минимизации светопропускания окна в заданную дату и время суток с учётом азимута окна, географических координат здания, сезонного и суточного распределения интенсивности солнечной радиации описан в [17]. Оптимизация заключается в определении оптимального угла наклона решёток, оптимального характеристического угла и ширин полос решёток фильтра. Для разных соотношений ширин полос обеих решёток установлены границы диапазонов с постоянным максимальным, уменьшающимся, постоянным минимальным и увеличивающимся коэффициентом светопропускания, из которых состоят угловые характеристики светопропускания фильтра, для этих диапазонов получены уравнения для расчёта коэффициента светопропускания.
В [18] продемонстрирована эффективность окна с двойным остеклением со встроенным оптическим фильтром для обеспечения динамического контроля дневного освещения и солнечной энергии, проходящей в помещение, в течение всего года без вмешательства человека и использования устройств перераспределения дневного света. Для 15-го числа каждого месяца года проведено численное моделирование светопро-пускания смарт-окна с фильтром с рассчитанным минимальным пропусканием в день (15 июня) и время (11 ч. 30 мин.) максимальной солнечной радиации в г. Оренбурге и построены соответствующие угловые характеристики светопропускания для тех часов дня, когда солнечные лучи попадают в окно. Среднемесячные значения теоретического и реального (скорректированного с учётом отражения и поглощения) коэффициентов светопропускания фильтра являются минимальными в июне, а также в мае и июле, когда защита от солнечной радиации наиболее востребована. При необходимости изменять угловые характеристики светопропускания фильтра, например, для защиты от горизонтальных солнечных лучей при низкой высоте стояния Солнца зимой (когда лучи падают на сетчатку глаза практически под нулевым углом и вызывают неприятные ощущения), одна из решёток должна быть передвижной по вертикали для изменения характеристического угла фильтра.
Результаты расчёта оптимального угла наклона решёток, ширин полос и характеристического угла фильтра для минимизации направленного светопро-пускания окна в заданную дату и время суток с учётом азимута окна, географических координат здания, сезонного и суточного распределения интенсивности солнечной радиации экспериментально подтверждены на модели смарт-окна с двойным остеклением [19]. Экспериментальные данные хорошо согласуются с результатами расчётов. Несмотря на то, что смарт-окна с поглощающими, отражающими или рассеивающими полосами, нанесёнными на поверхность стекла на постоянной основе, рассмотренные в [15–18], оптимально регулируют светопропускание в светлое время суток в течение всего года, они не могут быть применены в обычных окнах из-за постоянного ограничения видимости через них. Их можно использовать лишь в световых проёмах, предназначенных для пропускания солнечного света в необходимом количестве в течение суток, например, в световых люках на крыше, в промышленных зданиях. «Непропускающие» полосы предпочтительно должны быть изготовлены из хромогенных материалов, которые активизируются (становятся непрозрачными) при достижении достаточно высокой температуры (термохромные материалы окрашиваются, термотропные материалы начинают рассеивать свет), при повышении интенсивности коротковолнового диапазона спектра солнечного излучения (фотохромные материалы), при принудительном изменении концентрации водорода внутри стеклопакета (газохромные окна). Многослойные электрохром-ные стёкла, а также стёкла с жидкими кристаллами и взвешенными частицами прозрачны при подаче слабого электрического напряжения, без напряжения они непрозрачны. Светопропускание всех этих видов современных смарт-стекол при неактивном состоянии незначительно меньше, чем у обычных стёкол. Таким образом, при низкой температуре, или низкой интенсивности солнечной радиации (зимой, а также в утреннее и вечернее время), или при отключённом напряжении смарт-окна со встроенными фильтрами будут прозрачными [19]. Современные смарт-стекла полностью покрывают всю площадь окна, в результате, когда окно переключается в непрозрачный режим, оно перестает пропускать не только прямое солнечное излучение, но также рассеянное и отражённое излучение. Однако для комфортного дневного освещения и инсоляции в помещении желательно защитить его только от прямых солнечных лучей, пропуская рассеянный свет неба и свет, отражённый от поверхности Земли и окружающих зданий [19]. Именно такую функцию и обеспечивают смарт-окна со встроенными фильтрами, при этом их угловые и временные характеристики с чередующимися диапазонами с максимальным, уменьшающимся, минимальным и увеличивающимся светопропусканием можно предварительно рассчитать, чтобы динамически адаптировать к положению Солнца, избегая максимумов и минимумов све-топропускания в нежелательное время дня.
2. Метод расчёта фильтра
Оптимальные геометрические параметры фильтра определяются для минимизации светопропускания окна в конкретную дату и время суток при заданных значениях широты и долготы здания, азимута ориентации окна, расстояния между решётками, показателя преломления и натурального коэффициента поглощения стекла. Угол наклона решёток фильтра определяется на основе метода [17], включающего: (1) вычисление высот и азимутов Солнца для выбранной даты (для дня максимальной солнечной радиации или середины самого жаркого периода года) и определение траектории движения Солнца относительно окна в течение светового дня; (2) линейную аппроксимацию траектории Солнца в выбранном интервале времени светового дня; (3) определение угла наклона решёток как угла между аппроксимирующей и горизонтальной линиями.
При сложной траектории движения Солнца для расчёта параметров фильтра определяется проекция угла 9 между лучом, прошедшим через входную решётку (на рис. 2 этот угол принадлежит плоскости падения е ), и перпендикуляром к плоскости окна на плоскость е', перпендикулярную к полосам решёток. Для одинарного остекления этот угол равен углу преломления 9 = arcsin(sin 0 / n ), для двойного - углу падения 9 = @ . Проекция угла преломления или падения определяется по формуле
9 ' = arctan ( {^ x 2 + у 2 cos [у- arctan ( x/y ) ] } Js ) , (1)
где x и у – координаты следа точки падения, оставляемого на выходной поверхности.

Рис. 2. Проекция 9' угла преломления (падения) 9 на плоскость, перпендикулярную к полосам решёток
На рис. 3 показана схема смарт-окна с наклонными решётками, адаптированными к траектории движения Солнца по отношению к окну, и координаты следов точки падения O 1 при характеристическом ( O с ) и произвольном ( O 2 ) углах падения.
Для выбранной расчётной даты координаты следа точки падения рассчитываются при одинарном и двойном остеклении, соответственно, по формулам x = 5 tan
( A - A o )
. f . 0 arcsin I sin— к n
У =
0 c и произвольном угле падения 0 (расстояние между точками O 2 и O c на рис. 1 или проекция этого расстояния на плоскость, перпендикулярную к полосам решёток, на рис. 3) рассчитывается по формуле
s tan
arcsin I sin—

. ( . 0 arcsin I sin— к n.
а= 7 x min + y min ■ cos

cos
x = 5 tan ( A - A 0 ) ,


f „0 arcsin I sin — к n.
, (2)
-7 x 2 + y 2 cos
t x Y- arctan— .
У _
tan ( 0-| A - A 1) cos ( A - A 0 )
где А – азимут Солнца, град; А 0 – азимут ориентации окна, град; n – показатель преломления стекла.

Рис. 3. Схема смарт-окна с наклонными решётками
Для расчёта характеристического угла по формуле (1) берутся значения x = x min и у = у min , рассчитываемые по формулам (2) при значениях азимута Солнца и угла падения, соответствующих времени максимальной солнечной радиации, когда направленное светопропускание окна должно иметь минимальное значение. Тогда характеристический угол равен
0 c = arctan (7 x mn + y mn x
X cos
Y
arctan
x min
"Ik
^^^^^^B
y min
s
-17
Расположение полос решёток на рис. 3 соответствует характеристическому углу падения – точка падения O 1 находится в середине пропускающей полосы входной решётки, а её след O c – в середине непропускающей полосы выходной решётки. Сдвиг А между следами входной решётки фильтра на плоскости выходной решётки при характеристическом угле
т r = т
f i - 0,5
sin 2
sin 2
0 - arcsin ( sin 0/ n )
0 + arcsin ( sin 0/ n )
tan 2
tan 2
Ширины полос с 1 , с 2 , с 3 и с 4 рассчитываются по методу, приведённому в [17] (первые две формулы для одинарного и двойного остекления соответственно)
25 sin 0 c 25 sin 0 av c 3 = I ■ I ■ ’
7n2 + sin2 0c nn2 + sin2 0av c 3 = 25 tan 0 c - 25 tan 0 av, max c4 = C 3--------, (5)
Tmax c 2 = C 3 -Tmin ( c 3 + c 4 ) , c 1 = c 3 + c4 - c 2 , где 0av - заданный средний угол (0av < 0c), соответствующий среднему значению коэффициента свето-пропускания (тav = 0,5 (Tmin+Tmax)); Tmin и Tmax - заданные минимальные и максимальные теоретические значения светопропускания фильтра.
Теоретическая угловая характеристика фильтра (зависимость коэффициента направленного свето-пропускания от угла падения солнечных лучей или соответствующего времени суток) состоит [15–17] из чередующихся диапазонов с постоянным максимальным, уменьшающимся, постоянным минимальным и увеличивающимся коэффициентом светопропуска-ния. В диапазонах с минимальным и максимальным коэффициентом направленного светопропускания она рассчитывается, соответственно, по формулам
c i - c 4 c 3
T ■ =-------. T______ =------- min ,max ci + c 2 ci + c 2
Диапазоны с уменьшающимся и увеличивающимся светопропусканием рассчитываются по формуле
= |A|- 0,5 c 2 + 0,5 c 3
c i + c 2
.
Реальный (скорректированный) коэффициент направленного светопропускания т r рассчитывается с учетом отражения по уравнениям Френеля и поглощения по закону Бугера–Ламберта
0 - arcsin ( sin 0/ n ) 0 + arcsin ( sin 0/ n )
q
exp
-
1 + sin 2 0
■a b\ ~;-----;—
V n 2 - sin 2 0
где а - натуральный коэффициент поглощения стекла, мм–1; q – показатель степени (равен 2, 4 и 6 для одинарного, двойного и тройного остекления соответственно); b – суммарная толщина слоёв остекления, мм.
3. Метод расчёта дифракции на двух решётках фильтра
Угол падения @ солнечных лучей на плоскость окна непрерывно изменяется и зависит от высоты стояния Солнца h и его азимута А. В случае вертикального окна угол падения можно вычислить по частному случаю первой теоремы косинусов для трёхгранного угла, когда двухгранный угол напротив искомого плоского угла @ равен 90° (между вертикальной плоскостью, проходящей через падающий луч, и горизонтальной азимутальной плоскостью) cos @ = cos h ■ cos(А - А о), где разность азимутов Солнца и окна А – А0 представляет собой азимут Солнца относительно окна при данном его положении на небе (в данный момент времени суток). Для выбранной расчётной даты высота и азимут Солнца вычисляются на онлайн-калькуляторе по заданным значениям широты и долготы здания (например, на сайте .
Как известно [14], при наклонном падении пучка параллельных лучей на решётку дифракционная картина наблюдается на плоскости, перпендикулярной к падающим лучам. В простейшем случае, когда угол падения пучка лучей на плоскость решётки изменяется, а направление пучка всегда остаётся перпендикулярным к щелям дифракционной решётки (рис. 1), период решётки d уменьшается согласно формуле d ' = d ■ cos @ [14] . Однако в рассматриваемом случае необходимо учитывать, что углы падения солнечных лучей одновременно изменяются как в азимутальной (горизонтальной), так и в вертикальной плоскости, проходящей через падающий луч. Кроме того, решётки расположены под углом у к горизонтали.
Периоды входной и выходной решёток фильтра равны d = с 1 + с 2 = с 3 + с 4 . Дифракционная картина после прохождения пучка параллельных лучей через входную решётку будет наблюдаться при уменьшении периода этой решётки по формуле d' = d ■ cos @' (рис. 4).
Ширины полос входной решётки уменьшатся пропорционально. При наклонном падении пучка параллельных лучей ширина пропускающей полосы и период решётки (после преобразований с учётом формулы (1)) будут, соответственно, составлять c1s
5 2 + COS 2

+ У 2 )
ds
с 1 =
d ' =
5 2 + COS 2

+ У 2 )
где координаты х и у при одинарном и двойном остеклении вычисляются по соответствующим формулам (2).

Рис. 4. Схема дифракции на входной решётке фильтра
На рис. 4 плоскость падения солнечного луча в данный момент времени (плоскость s на рис. 2) спроецирована на плоскость, перпендикулярную к полосам решёток (плоскость s' на рис. 2). Условно показаны направления максимумов интенсивности до третьего порядка и проекция ф 1' угла дифракции, соответствующего максимуму первого порядка.
Согласно формуле для дифракционных решёток [20], после прохождения солнечного света через входную решётку с количеством щелей, т.е. пропускающих полос, N , интенсивность света / ( ф ) при угле дифракции ф будет равна
/ ( ф ) = / о
| п C1 sin ф 1
. l X J sin
[ n c 1 sin ф 1
l X )
X
. [ Nпd sin ф 1 sin lx)
. ) пd' sin ф) sin l X )
, (10)
где I 0 – интенсивность света в направлении главного дифракционного максимума; X - длина волны.
При определении коэффициента направленного светопропускания по формулам (6) и (7) с последующей коррекцией по формуле (8) дифракция не учитывается. Понятие «направленное светопропускание» допустимо только при очень малых углах дифракции ф , которые могут наблюдаться при значительном превышении периода решётки над длиной волны. Расчёты показывают [15, 17], что периоды решёток оптических фильтров для окон превышают длины волн видимого излучения на четыре порядка при одинарном остеклении и на пять порядков при двойном остеклении. В соответствии с формулой d' = d ■ cos9' влияние дифракции усиливается при увеличении углов падения солнечных лучей на окно, особенно при приближении к «скользящим» лучам, т. е. к углам падения 90°. Однако при таких углах падения прошедшие через оконную конструкцию солнечные лучи попадают в зону в непосредственной близости к окну и необходимость регулирования светопропускания отсутствует.
Направления максимумов интенсивности определяются из условия [21] d' (sin ф m -sin 9 ) = m X , где m -порядок максимума. Это выражение справедливо для неограниченной решётки, однако в оконных конструкциях количество периодов достаточно большое и составляет порядка 50 при двойном остеклении и 375 при одинарном остеклении. Соответствующая углу дифракции ф m интенсивность вычисляется по формуле (10).
При значительном превышении периода решётки над длиной световой волны и, соответственно, при малых углах дифракции полная интенсивность после прохождения входной решётки [21] равна I прош1 = с 1 I пад / d . Эта интенсивность складывается из интенсивности света в направлении главного максимума и интенсивности дифрагированного света I прош1 = I 0 + I диф , где I 0 = ( с 1 / d )2 I пад . Тогда полная интенсивность прошедшего света [19] равна
двух краёв непропускающей полосы выходной решётки в пределах одного её периода. При больших расстояниях влиянием дифракции можно пренебречь. Интенсивность дифрагированного света уменьшится по сравнению с рассчитанной по формуле (12) на величину (2 × 780/ d )2 I пад . Тогда интенсивность света, дифрагированного на входной решётке и проходящего через выходную решётку, равна (при размерности длины в мм)
I диф1
C i ( d - C i ) - 2,4336 x i0 - 6 d 2 I пад .
Таким образом, из формул (11)–(13) следует, что интенсивность света, прошедшего через входную решётку и доходящего до выходной решётки с учётом дифракционных потерь, равна
I прош1

п 2 m = 1 m 2
. , п mci | r Sin 2 ---— 1 1 пад , d J
где второе слагаемое в скобках показывает долю дифрагированного света.
Интенсивность дифрагированного света связана с интенсивностью I пад падающего на входную решётку света соотношением
I диф1
C i ( d - C i ) d 2 I пад .
Максимальная интенсивность дифрагированного света по этому выражению достигается при с 1 = 0,5 d , т.е. когда ширины пропускающих и непропускающих полос решётки равны с 1 = с 2 .
При малых углах дифракции солнечное излучение, прошедшее через входную решётку, можно считать направленным (не рассеянным). Для оконных конструкций рассеянный дифрагированный свет допустим ввиду того, что при необходимости защиты от солнечного излучения невозможно одновременно обеспечивать хорошую видимость через окно.
Распространение параллельного пучка света, прошедшего через решётку, можно считать прямолинейным и подчиняющимся геометрической оптике на расстояниях от решётки r << D 2 / X , где D – диаметр пучка. Для реальных оконных конструкций и расстояний между входной и выходной решётками фильтра это неравенство полностью удовлетворяется.
Дифрагированный на входной решётке свет частично проходит через пропускающие полосы выходной решётки (рис. 4), остальная часть блокируется непропускающими полосами. Соответственно, уменьшается коэффициент светопропускания, вычисленный по формулам (6) и (7) с коррекцией по формуле (8). Оценим снижение коэффициента светопро-пускания из-за дифракции на расстояниях максимальной длины волны видимого излучения 780 нм от
прош1
C i2 + C i c 2 - 2,4336 x i0 - 6 d 2

Без учёта дифракционных потерь
прош1

+ C i C 2
d 2

Отношение выражений в скобках в формулах (14) и (15) показывает уменьшение коэффициента свето-пропускания фильтра из-за дифракции на входной решётке.
Проекция расстояния 780 нм от края непропускающей полосы выходной решётки на плоскость рис. 4 при расстоянии s между решётками фильтра равна
7,8 x io - 4 cos 9' = 5 [ tan ( 9' + ф' ) - tan 9'] . (i6)
Тогда проекция ф' угла дифракции ф на расстоянии 780 нм от края полосы определяется как
(9'
ф = arctanI 7,8xi0-4cos—+ tan9 1-9 .(i7)
I
Угол дифракции ф определяется из последней
формулы (2)
у - 7,8 x i0 - 4
откуда
ф = arctan
5 tan ( ® + ф-| A - A o I) cos ( A - Ao )
- cos ( A - A o )( у - 7,8 x io - 4 )
s
-0 + 1 A - A o |.
Из условия d' sin ф m =± m X определяются порядки максимумов интенсивности в пределах угла ф . По формуле (10) определяются соответствующие максимумы интенсивности и интенсивность в направлении угла ф .
Реальный коэффициент направленного светопро-пускания, рассчитываемый по формуле (10), равен отношению интенсивности, прошедшей через выходную решётку, к интенсивности, падающей на входную решётку т r = I прош/ 1 пад (без учёта потерь дифрагированного света на краях непропускающих полос выходной решётки). При учёте дифракции по формулам (14) и (15) в формуле (10) для расчёта коэффициента светопропускания добавляется соответствующий (дифракционный) сомножитель
Т r
Г
1 - 0,5
I
sin 2
sin 2
0- arcsin ( sin 0/ n ) 0 + arcsin ( sin 0/ n )
tan 2
+ tan2
x exp
q
0 - arcsin ( sin 0/ n )
0 + arcsin ( sin 0/ n )
-
, 1 + sin 2 0
■ a b.--------------
V n 2 - sin 2 0
2,4336 x 10 - 6
c 1 c 2
Дифракция на выходной решётке не влияет на значение коэффициента светопропускания, однако доля дифрагированного (рассеянного) света увеличится. По формуле (20) оцениваются дифракционные пределы углов падения и ширин пропускающих и непропускающих полос входной решётки, при которых влиянием дифракционного сомножителя можно пренебречь.
4. Результаты расчёта
Геометрические параметры фильтра определены по методу, приведённому в параграфе 2, для г. Оренбурга (51°47'00''N, 55°06'00''E, UTC+05:00) для 15.06.2018 (день максимальной солнечной радиации в г. Оренбурге) и 11 ч 30 мин (время максимальной солнечной радиации в г. Оренбурге) для окна с двойным остеклением при его азимуте 120 ° и расстоянии между решётками s =16 мм. Высота и азимут Солнца рассчитаны относительно 10 ч 49 мин, когда азимуты Солнца и окна равны ( A = A 0 ), через каждый час/полчаса с 6 ч 49 мин до 12 ч 19 мин. Линейная аппроксимация траектории Солнца проведена в интервале времени от 7 ч 49 мин до 12 ч 19 мин. Оптимальный угол наклона решёток фильтра составляет у = 42 ° .
По формуле (3) определён характеристический угол фильтра 0 c = 28,71 ° при x = x mm = 3,688 мм и у = у min = –15,1118 мм, рассчитанных по формулам (2) для момента времени максимальной солнечной радиации 11 ч 30 мин. Ширины полос решёток фильтра рассчитаны по формулам (5) при 0 av =19,807°, Т тШ = 0,45 и T max = 0,7: c 1 =6,4286; c 2 = 2,1429; c 3 = 6 и c 4 = 2,5714 мм.
По методу, приведённому в параграфе 3, оценено влияние дифракции на светопропускание оптического фильтра с определёнными выше параметрами. При углах падения солнечных лучей более 70° резко воз- растает коэффициент отражения, а также прошедшие лучи попадают в зону около окна, т.е. необходимости регулирования светопропускания нет. При принятых выше расчётных параметрах в 12 ч 49 мин угол падения составляет 0 = 69,9842°. При этом угле падения ширина пропускающей полосы входной решётки и период решётки по формулам (2) и (9) равны с 1 = 6,3193 и d' =8,4258 мм.
Полная интенсивность света, прошедшего через входную решётку, по формуле I прош1 = с 1 I пад / d равна I прош1 = 0,75 I пад . Она складывается из интенсивности света в направлении главного максимума и интенсивности дифрагированного света (по формулам (11) и (12)) I 0 = 0,5625 I пад и I диф =0,1875 I пад .
Угол дифракции на расстоянии 780 нм от края непропускающей полосы выходной решётки по формуле (19) равен ф = 0,0017°. Максимумы первого и второго порядков из условия d' sin ф m =± m X при X = 780 нм составляют ф 1 =0,0053° и ф 2 = 0,0106°. Т.е. в пределах угла ф нет максимумов интенсивности, кроме главного максимума. Вычисленные по формуле (10) при N =50 интенсивности, соответствующие углам дифракции ф , ф 1 и ф 2 , равны 1 ( ф ) = 0,5849 1 0 , I ( ф 1 ) = 0,0169 1 0 и 1 ( ф 2 ) = 0,0132.
Главный дифракционный максимум является резким (75% от полной интенсивности света, прошедшего через входную решётку) и узким – порядка длины волны на плоскости выходной решётки (из сравнения угла дифракции ф = 0,0017° и направления минимума между максимумами нулевого и первого порядков ф 1 /2 = 0,00265°).
Сдвиг между следами входной решётки фильтра на плоскости выходной решётки при характеристическом угле 0 c и угле падения 0 = 69,9842° по формуле (4) равен А =11,78 мм. Теоретический коэффициент светопропускания по формуле (7) равен т = 0,5973.
По формулам (8) и (20) рассчитан реальный коэффициент светопропускания без учёта и с учётом дифракции. В обоих случаях они равны т r = 0,3783. Коэффициенты практически не различаются, т.к. дифракционный сомножитель в формуле (20) составляет 0,9999998. При угле падения 0 = 69,9842° значение дифракционного сомножителя в формуле (20) было бы сравнимым с произведением второго и третьего сомножителей, оценивающих вклад отражения и поглощения, при значениях ширин полос в несколько нанометров, что невозможно в рассматриваемых оконных конструкциях. При принятых для расчёта значениях ширин полос при углах падения 85 °, 87 ° и 89 ° реальный коэффициент светопропускания без учёта и с учётом дифракции составляет, соответственно, 0,0848, 0,0365 и 0,0050, т.е. при приближении к углу падения 90° коэффициент стремится к нулю.
Для окна с одинарным остеклением при расстоянии между решётками (толщине стекла) s =4 мм и ширинах полос решёток фильтра c 1 =3; c 2 = 1; c 3 =2,5
и c 4 = 1,5 мм (фильтр с такими параметрами рассмотрен в [15]) значение дифракционного сомножителя в формуле (20) составляет 0,9999992 и влиянием дифракции также можно пренебречь.
Заключение
Двухрешёточные оптические фильтры для смарт-окон исследованы на предмет влияния дифракции на коэффициент их направленного светопропускания. Такие фильтры функционируют в условиях непрерывного изменения углов падения солнечных лучей при сложной траектории движения Солнца относительно окна. Для таких условий разработан метод расчёта дифракции солнечного света при последовательном прохождении его через входную и выходную решётки фильтра. Для оценки уменьшения коэффициента направленного светопропускания фильтра из-за влияния дифракции в расчётную формулу введён дополнительный сомножитель. Значения дифракционного сомножителя, полученные при расчётах для смарт-окон с одинарным или двойным остеклением при обычных для них эксплуатационных условиях, позволяют считать, что влиянием дифракции на све-топропускание решёточного фильтра можно пренебречь. При приближении угла падения солнечных лучей к 90° относительный период обеих решёток фильтра уменьшился бы, и влияние дифракции стало бы значительным, однако такие углы падения конструктивно невозможны из-за углубления оконных проёмов в стену здания. К тому же при углах падения более 70° резко возрастает коэффициент отражения, а также прошедшие лучи попадают в зону около окна, т.е. усиление влияния дифракции не будет препятствием для использования решёточных оптических фильтров в окнах.
Список литературы Дифракция в решёточных оптических фильтрах с угловым селективным светопропусканием
- Lemarchand, F. Increasing the angular tolerance of resonant grating filters with doubly periodic structures / F. Lemarchand, A. Sentenac, H. Giovannini // Optics Letters. - 1998. - Vol. 23. - Р. 1149-1151. - DOI: 10.1364/OL.23.001149
- Sentenac, A. Angular tolerant resonant grating filters under oblique incidence / A. Sentenac, A.-L. Fehrembach // Journal of the Optical Society of America A. - 2005. - Vol. 22. - Р. 475-480. - DOI: 10.1364/JOSAA.22.000475
- Chen, C.-H. Wide-field-of-view GaAs/AlxOy one-dimensional photonic crystal filter / C.-H. Chen, K. Tetz, W. Nakagawa, Y. Fainman // Applied Optics. - 2005. - Vol. 44. - Р. 1503-1511. - DOI: 10.1364/AO.44.001503
- Yang, F. Integrated 2D photonic crystal stack filter fabricated using nanoreplica molding / F. Yang, G. Yen, B.T. Cunningham // Optics Express. - 2010. - Vol. 18. - Р. 11846-11858. - DOI: 10.1364/OE.18.011846
- Di Falco, A. Optical metasurfaces with robust angular response on flexible substrates / A. Di Falco, Y. Zhao, A. Alu // Applied Physics Letters. - 2011. - Vol. 99. - 163110. - DOI: 10.1063/1.3655332
- Стафеев, С.С. Бинарные дифракционные решётки для управления поляризацией и фокусировкой лазерного света [обзор] / С.С. Стафеев, А.Г. Налимов, Л. О'Фаолейн, М.В. Котляр // Компьютерная оптика. - 2017. - Т. 41, № 3. - С. 299-314. -
- DOI: 10.18287/2412-6179-2017-41-3-299-314
- Кадомина, Е.А. Резонансные фотонно-кристаллические структуры с дифракционной решёткой для измерения показателя преломления среды / Е.А. Кадомина, Е.А. Безус, Л.Л. Досколович // Компьютерная оптика. - 2016. - Т. 40, № 2. - С. 164-172. -
- DOI: 10.18287/2412-6179-2016-40-2-164-172
- Лушников, Д.С. Расчёт и исследование угловой и спектральной селективностей для полнопараллаксных голограммных стереограмм / Д.С. Лушников, А.Ю. Жердев, В.В. Маркин, С.Б. Одиноков, А.В. Смирнов // Компьютерная оптика. - 2016. - Т. 40, № 6. - С. 802-809. -
- DOI: 10.18287/2412-6179-2016-40-6-802-809
- Chen, L. Polarization and angular dependent transmissions on transferred nanomembrane Fano filters / L. Chen, Z. Qiang, H. Yang, H. Pang, Z. Ma, W. Zhou // Optics Express. - 2009. - Vol. 17. - Р. 8396-8406. -
- DOI: 10.1364/OE.17.008396
- Yang, F. Optical characterization of liquid crystals by means of half-leaky guided modes / F. Yang, J.R. Sambles // Journal of the Optical Society of America B. - 1993. - Vol. 10. - P. 858-866. -
- DOI: 10.1364/JOSAB.10.000858
- Jahan, F. Investigation of angular selective optical properties of silver/titanium oxide cermet thin films / F. Jahan, G.B. Smith // Thin Solid Films. - 1998. - Vol. 333. - Р. 185-190. -
- DOI: 10.1016/S0040-6090(98)00844-X
- Palmer, S. Angular selective optical properties of thin films: Measurement of polar and azimuthal transmittance / S. Palmer, G.W. Mbise, G.A. Niklasson, C.G. Granqvist // Solar Energy Materials and Solar Cells. - 1996. - Vol. 44. - P. 397-403. -
- DOI: 10.1016/S0927-0248(95)00176-X
- Дюрелли, А. Анализ деформаций с использованием муара / А. Дюрелли, В. Паркс; пер. с англ. - М.: Мир, 1974. - 359 с.
- Ландсберг, Г.С. Элементарный учебник физики. Т. 3. Колебания и волны. Оптика. Атомная и ядерная физика / Г.С. Ландсберг. - 13-е изд. - М.: Физматлит, 2003. - 656 с.
- Закируллин, Р.С. Оптический фильтр с угловой селективностью светопропускания / Р.С. Закируллин // Оптический журнал. - 2013. - Т. 80, Вып. 8. - С. 16-24.
- Zakirullin, R.S. Creating optical filters with angular-selective light transmission / R.S. Zakirullin // Applied Optics. - 2015. - Vol. 54, Issue 21. - P. 6416-6419. -
- DOI: 10.1364/AO.54.006416
- Zakirullin, R.S. Optimized angular selective filtering of direct solar radiation / R.S. Zakirullin // Journal of the Optical Society of America A. - 2018. - Vol. 35, Issue 9. - P. 1592-1598. -
- DOI: 10.1364/JOSAA.35.001592
- Zakirullin, R.S. Optical filter for smart windows with angle-selective light transmission / R.S. Zakirullin // Journal of Optical Technology. - 2019. - Vol. 86, Issue 5. - P. 278-283. -
- DOI: 10.1364/JOT.86.000278
- Zakirullin, R.S. A smart window for angular selective filtering of direct solar radiation / R.S. Zakirullin // Journal of Solar Energy Engineering, Transactions of the ASME. - 2019. - Vol. 142, Issue 1. - 011001. -
- DOI: 10.1115/1.4044059
- Стюард, И.Г. Введение в фурье-оптику: Пер. с англ. / И.Г. Стюард. - 13-е изд. - М.: Мир, 1985. - 182 с.
- Сивухин, Д.В. Общий курс физики. Т. 4. Оптика / Д.В. Сивухин. - М.: Наука, 1980. - 752 с