如何使用comsolMultiphysics®为为各向建模建模建模建模

2017年年12月4日

1669年晴朗夜晚,伊拉斯mus bartholinus教授教授把一块放在长凳长凳上冰岛方解石解石晶体晶体。突然发现这个光学折射,是是一一从晶体射出时分裂成成两两个个平行平行引起的。。我们我们将将在在本篇博客中

了解各向材料

伊拉斯mus bartholinus到穿过直射称为寻常光线,另寻常寻常寻常寻常晶体时发生弯曲弯曲的的的,是光束光束弯曲弯曲弯曲弯曲弯曲弯曲弯曲弯曲弯曲弯曲弯曲各向,例如上述方和实验的晶体。。

光线穿过向异性材料。。
穿过各异性的寻常和非常光线。

在环境,当当电磁波束向异性介质材料传播传播时,将电偶极子的电荷这种导致向异性介质内出现应场应场应场应场应场,其中其中应场应场种种波波的的不同不同((寻常波波波波波波波波波

在发生偏振偏振轨迹传播。

在硅波导引入各向异性

之前的一博客文章cmos制造与兼容性兼容性兼容性兼容性兼容性兼容性兼容性兼容性兼容性的如何光子集成芯片芯片中广泛广泛。具有具有各各向同性特性的硅硅材料材料材料于于由于具有和偏振的效应光学特性,各特性特性特性异性异性材料。。

制造波导中,硅光子学中中中各各向有光学光学光学光学效应效应效应效应效应效应((n〜1.44))(n〜1.44)(n〜1.44)(n〜1)(n〜1)()在内实现非常的光学传感器。。

光学传播模式

为了对向介质定性,研究分析,研究人员了了光能如何在波导内传播(((())。。平面,我们e^{x} _ {p,q}e^{y} _ {p,q}(2),定义,其中,其中Xy表示偏振,,p表示Xy坐标中最大值。

想象:你你在幅幅e^{x} _ {2,1}““风景画”(“上”),“风”(“风”沿()沿±x方向方向从从从从-X到 +X方向时,你会两个峰。当从从-y方向朝+y方向移动,可以可以观察两峰值。。

EX11平面波导模式分析。
EY11平面平面的。。。
EX12平面波导模式分析。
EY12平面波导模式。。。
EX21平面波导模式分析。
EY21平面波导模式。。。

平面波导模式。顶行,从从::e^{x} _ {1,1}e^{y} _ {1,1}。中间行,从左::e^{x} _ {1,2}e^{y} _ {1,2}。底行,从左::e^{x} _ {2,1}
e^{y} _ {2,1}
箭头图表示;(蓝色外流(高高

在comsolMultiphysics®软件软件各向结构结构结构

在使用光束之前,了解之前了解模式的/包层/包层包层包层包层尺寸内持续存在非常非常重要重要的的comsolMultiphysics®软件)进行,有助于,有助于分别和分析模式和。。。。

引入对角各异性

对各进行模态分析都个个复数值,而复数值,而复数值复数值向向异性异性异性异性材料材料的这里,,张量指的是个个个个个个,它矩阵矩阵矩阵\ Epsilonxx,,,,\ Epsilon是的,,,,\ EpsilonZ Z)和和角线((\ Epsilonxy,,,,\ Epsilonxz,,,,\ EpsilonYX,,,,\ EpsilonYz,,,,\ EpsilonZX,,,,\ Epsilonzy)项,如下所。

\ epsilon = \ begin
{bmatrix}
\ epsilon_ {xx}&\ epsilon _ {xy}&\ epsilon _ {xz} \\
\ epsilon _ {yx}&\ epsilon _ {yy}&\ epsilon _ {yz} \\
\ epsilon _ {zx}&\ epsilon _ {zy}&\ epsilon _ {zz}
\ end {bmatrix}

,,对于来说,你可以系系系,在坐标,在坐标系中个坐的,相应,相应,介电常数中对角线元素三个个值被材料的的主介。

:单轴主要两:单轴单轴单轴双轴晶体。在选择适当坐标系系((其中其中只有常数常数张量张量是元素元素是元素的的的的的的)单轴晶体仅考虑,即,即\ Epsilonxx=\ Epsilon是的=(((no2,,,,\ EpsilonZ Z=(((ne2,其中none分别是折射率非常折射率,当,当\ epsilon_
{xx} \ neq \ epsilon_ {yy} \ neq \ epsilon_ {zz}
时,它它称为双轴晶体

为了将纳入的角度角度,我们可以角度扩展扩展这篇硅光子学设计文章文章中隐埋导示例。我们可以对方对方形芯形芯长度分别分别分别分别分别分别分别分别分别分别分别分别分别分别分别分别分别分别的的的的的的二维横((()1.55 [um]。

具有各异性芯的的注释模型。
三维隐埋肋光导示意图,其中对二维横执行分析强度图箭头图分别代表代表电场的模式极化。
光波导的向异性芯特写视图。
标出沿X(((((((())和主轴)的的肋波。

在传统材料情况,我们我们假设即即即即C(()沿着沿着X(((),并并将角常数常数常数\ Epsilon是的\ EpsilonZ Z项(与C((())视为视为折射率(2〜2.31)。。C轴的\ Epsilonxx分量元素认为非常折射率((〜1.47992(〜2.19)(3)(3)。参考文献。。,非此外项被认为是零((示示示示)2)。。模式是上面所所所示的的的种模式请注意注意:“”nxx-n是的”称为,其中,其中nxx=\ sqrt {\ epsilon_ {xx}}}n是的=\ sqrt {\ epsilon_ {yy}}

\ epsilon =
\ begin {bmatrix}
2.19&0&0 \\
0&2.31&0 \\
0&0&2.31
\ end {bmatrix}

具有对元素的介电常数张量。

色散曲线

通过模式,我们可以地理解。。此外此外此外此外此外此外此外此外。。。。进行进行定量定量也也方便。通常色色散曲线表示有效折射率

对角各向异性

0.5 um到4um范围内对长度进行参数化模式模式分析分析分析分析分析分析分析分析分析分析分析分析模式模式模式出出各向的的色散曲线曲线曲线曲线曲线曲线曲线各向异性(即\ Epsilonxx= 2.19,\ Epsilon是的=\ EpsilonZ Z= 2.31,(((所有都))。。。将结果与与与等等的结果结果结果((进行进行进行进行进行进行进行。。。

横向各异性芯的色散。。
横向各异性芯的色散。。

(xy平面平面平面)

当光轴即,,C(轴)位于xy平面并与与X轴成\ theta角,对,对分量\ Epsilonxx,,,,\ Epsilon是的,,,,\ EpsilonZ Z和非对分量\ Epsilonxy\ EpsilonYz非零,而其余为。可以使用如下所的矩阵矩阵[r]来来计算计算计算全全相对相对介电张量张量张量C轴。由于C轴沿着主X轴,因此\ Epsilonxx是非常折射率(〜2.19),而\ Epsilon是的\ EpsilonZ Z是寻常((〜2.31)。。非角线元素元素\ Epsilonxy\ EpsilonYz是根据乘法导出,如下的,如下所。。

光轴位于xy平面平面光波导芯的。。。
c轴xy平面,与x与成成\ theta角。

\ epsilon = [r]。[\ epsilon_ {m}]。[r] ^t =
\ begin {bmatrix}
cos(\ theta)&-sin(\ theta)&0 \\
sin(\ theta)&cos(\ theta)&0 \\
0&0&1
\ end {bmatrix}
\ begin {bmatrix}
\ epsilon_ {xx}
&0&0 \\
0&\ epsilon_ {yy}&0 \\
0&0&\ epsilon_ {zz}
\ end {bmatrix}
\ begin {bmatrix}
cos(\ theta)&sin(\ theta)&0 \\
-sin(\ theta)&cos(\ theta)&0 \\
0&0&1
\ end {bmatrix}
=
\ begin {bmatrix}
(\ epsilon_ {xx})cos^2(\ theta) +(\ epsilon_ {yy}
)sin^2(\ theta)&(\ epsilon_ {xx})sin(\ theta)cos(\ theta) - (\ epsilon_ {yy})sin(\ theta)cos(\ theta)cos(\ theta)
(\ epsilon_ {xx}
)sin(\ theta)cos(\ theta) - (\ epsilon_ {yy})sin(\ theta)cos(\ theta)&(\ epsilon_ {yy}
)cos^2(\ theta) +(\ epsilon_ {xx})sin^2(\ theta)&0 \\
0&0&\ epsilon_ {zz}
\ end {bmatrix}

相对介电常数张量ε与旋转一处理,,在xy平面内将将C轴旋转角度\ theta

最后,对具有向和各同性的进行模式分析,其中分析分析,其中其中X轴轴轴度度度度度度度度度度度度度度角度度度度度度度度度度度度度度度度度度度度度度度度度度度度。从图图图图中中我们我们可以可以观察观察到到面面内范围内和包层长度执行参数化,同时扫描扫描扫描扫描扫描\ theta视为45°,绘制色散曲线所述,结果结果类似于各向的。。。

θ为0度度的。。。
θ为15度时模式。。。
θ为30度度模式。。。
θ为45度时模式。。。

θ= 0°(,θ= 15°()θ= 30°(θ)θ= 45°(,,,,,),包括(右右右右右右右分析分析分析分析分析分析分析分析模式模式芯内的磁场线。

((Yz())

最后,当,当光轴即,C(轴)位于yz平面平面与与y轴成\ phi角的向异性时,对对分量分量分量\ Epsilonxx\ Epsilon是的\ EpsilonZ Z和非对分量\ EpsilonYz\ Epsilonzy非零,而分量零。可以使用所示旋转旋转矩阵[r]来来来来计算计算相对相对介电常数常数C轴。由于C轴沿着主y轴,因此\ Epsilon是的是非常折射率(〜2.19),而\ Epsilonxx,\ EpsilonZ Z是寻常((〜2.31)。。非角线元素元素\ EpsilonYz\ Epsilonzy是根据乘法导出,如下的,如下所。。

光轴Yz平面平面光波导芯的。。。
C轴位于Yz平面,与x轴成\ phi角。

\ epsilon = [r]。[\ epsilon_ {m}]。[r] ^t =
\ begin {bmatrix}
1&0&0 \\
0&cos(\ phi)&-sin(\ phi)\\
0&sin(\ phi)&cos(\ phi)
\ end {bmatrix}
\ begin {bmatrix}
\ epsilon_ {xx}
&0&0 \\
0&\ epsilon_ {yy}&0 \\
0&0&\ epsilon_ {zz}
\ end {bmatrix}
\ begin {bmatrix}
1&0&0 \\
0&cos(\ phi)&sin(\ phi)\\
0&-sin(\ phi)&cos(\ phi)
\ end {bmatrix}
=
\ begin {bmatrix}
\ epsilon_ {xx}&0&0 \\
0&(\ epsilon_ {yy}
)cos^2(\ phi) +(\ epsilon_ {zz})sin^2(\ phi)&(\ epsilon_ {yy})sin(\ phi)cos(\ phi) - (\ phi) - (\ epsilon_
)sin(\ phi)cos(\ phi)\\
0&(\ epsilon_ {yy})sin(\ phi)cos(\ phi) - (\ epsilon_ {zz})sin(\ phi)cos(\ phi)&(\ phi)&(\ epsilon_ {zz})phi) +(\ epsilon_ {yy}
)sin^2(\ phi)
\ end {bmatrix}

相对介电常数张量ε通过旋转处理,表示在yz平面平面旋转\ phi

0.5 um到4um范围范围内芯的长度参数化来执行执行模式模式模式\ phi= 45°(,,C轴Yz平面,与,与y(45°角)(3)3)

纵向各异性芯的色散。。
纵向各异性芯的色散。。

关于各异性材料建模思考思考

(comsol多物理学)(((((和各来它们的色散关系。通过分析分析分析分析,我们类型分析分析分析分析分析分析分析分析分析分析分析类型类型在优化优化优化优化优化优化优化材料材料材料材料和和

阶跃折射率光纤这个教程帮助你入门,该你入门对三维光纤二维横截面的。

后续操作

如需这些,请请单击下面按钮进入进入进入进入案例下载下载页面页面页面页面页面页面

硅光子学系列博客更新列表

参考文献

  1. E. Hecht,光学,皮尔逊。
  2. E.A.J.Marcatili,“集成光学元件的介电矩形波导和定向耦合器”,贝尔系统。技术。J.,卷。48,第2071–2102页,1969年。
  3. M. Koshiba,K。Hayata和M. Suzuki,“具有任意张量介电常数的各向异性波导的有限元溶液”,”轻波技术杂志,卷。4,不。2,第121–126页,1986年。


评论(0)

留言
登录|注册
正在... ...
浏览comsol博客