关于comsol模拟的一些问答

1 1 近来用近来用 COMSOLCOMSOL 计算光子晶体光纤的模场分布 可是不知道计算光子晶体光纤的模场分布 可是不知道 PMLPML 的参数如何设置 以及边的参数如何设置 以及边 界条件怎么设置 计算出来的结果不对界条件怎么设置 计算出来的结果不对 实验室老板催得急 算不出来特别郁闷 不想读实验室老板催得急 算不出来特别郁闷 不想读 的心思都有了 的心思都有了 请用过的人帮帮忙吧 请用过的人帮帮忙吧 我也是用 comsol 算光纤的 关于 pml 层的设定问题 如果不考虑损耗的话 pml 层可以 不设 你可以试一试就知道了 pml 对模场分布基本没有影响 2 COMSOLCOMSOL MultiphysicsMultiphysics 如何模拟带隙光子晶体光纤 如何模拟带隙光子晶体光纤 要用 COMSOL Multiphysics 模拟带隙光子晶体光纤 也就是要加入 kz 可以用如下方法 1 用平面波模式 将模型边界条件改为电场 输入一个表达式的名字 例如 E1 2 定义该边界表达式 E1 菜单 选项 表达式 边界表达式 选择不同的边界 分别写入 该边界上电场 E1 的表达式 这样就能加入 kz 将所需的周期性边界方程写入 COMSOL Multiphysics 3 3 如何准确求光子晶体光纤的限制损耗即有效折射率的虚部如何准确求光子晶体光纤的限制损耗即有效折射率的虚部 我在模拟 PCF 时 为了求其限制损耗即有效折射率的虚部 在 PCF 结构的外面加了 PML 但是在加了 PML 后 却发现光束不能约束在纤芯中了 不知道哪里出了问题 还望各位高手给予指点 谢谢 8v p w q T V1W d attach 219885 attach T o 1O B0j P 加了 PML 后的结果如下 attach 219886 attach beautycatcher 发表于 2009 10 21 07 31 我也是初学 也在做一些光子晶体的方法 目前还不懂帮你顶顶 大家多多讨论 caoer 发表于 2009 10 21 11 17 有限元做光子 这个挺有新意 不过要注意是否适用 mahui 发表于 2009 11 5 09 59 能说一下有限元做光子为什么不合适吗 不过用 FDTD 做光子的还蛮多的 Feit 发表于 2009 11 5 12 22 PML 的几何不对 应该是加个六边形的 PML 才对吧 fangany 发表于 2009 11 8 13 29 纤芯比外面的小 当然有可能找到外面的那个模式 多找几个模式或者将外面的区域减小应该就可以了 shanyrain 发表于 2009 11 8 20 35 加个圆形的就可以了 PML 要考虑模型的对称性 比如这个模型可以只计算 1 4 或者 1 6 xwx000000 发表于 2009 11 13 22 31 楼主具体交流下怎么划分格点的 我算光子晶体光纤的模式 伪模很多阿 比如设 neff 1 5 附近寻找 设 200 个 它就给找出 200 个 neff 出来 n 在在 comsol 如何输入积分和二重积分符号 如何输入积分和二重积分符号 我要在选项 表达式里输入表达式 在表达式里输入积分符号和二重积分符号 这两个符号都怎么输入 Lynn An 发表于 2009 12 28 09 17 COMSOL 没有这种运算符 实现这个功能的是 积分耦合变量 如果你把这个积分耦合变量定义在三维实体上 那么就是三重体积分的含义 定义在二维域上 就是二重 面积分 定义在一维线上 就是单重的线积分 如果你要做的不是对空间的积分 而是对时间的积分 那么 就需要你把这个时间的积分整理成一个普通 的常微分方程 在 COMSOL 中可以添加 ODE 的约束 在在 comsol 里如何操作运行有效模场面积里如何操作运行有效模场面积 在模拟光子晶体光纤的时候想知道有效模场面积和波长的关系 我现在知道有效模场面积的计算公式 但 不知道在 COMSOL 里如何操作 运行 如何计算有效模场面积 先谢了 Lynn An 发表于 2009 12 28 09 19 后处理菜单中 有个域积分计算 你在那里面对你有效面积公式的分子和分母那两个积分分别进行计算 然后手工去做个除法 GYZ 发表于 2010 1 14 20 49 我还是不明白 能不能再详细一点 谢谢 Lynn An 发表于 2010 1 22 11 23 to GYZ 首先 你知道那个 Aeff 的公式不 如果不知道 可以去 SCI 数据库查一下 K Saitoh 做 PCF 的论文 02 03 年附近的 可以找到这个公式 我论坛回复没法用公式编辑器 所以只能让你自己去查 一查了 查一下也好 你写论文的时候总要引用这个公式的出处的 然后 这个公式是一个除法 分子和分母均是对电场分布的二重积分 当然积分域并不一样 积分变量也 不完全一样 虽然都是电场分布 但是分子分母有一个平方 一个是四次方 总之这个二重域积分的功能 COMSOL 里面内置了 可以简单的分别计算出来 joasm 发表于 2010 4 22 18 06 i s 本帖最后由 joasm 于 2010 4 22 18 11 编辑 i 附上 Aeff 公式 我用范例里的 color red step index fiber color 计算 Aeff 但结果有点不太正确 主任帮我看看 我的过程 Options integration coupling variable subdomain variable 直接对 Subdomain2 core 积分 name expression Int normE rfwv 2 分母 Int1 normE rfwv 4 分子 更新计算后 从 data display global 观看结果 但结果似乎不太正确 在在 comso 里 如何画正三角形和正六变形里 如何画正三角形和正六变形 求求大家帮忙 指点一下如何画正三角形和正六边形 非常着急 由于不会画 下面的论文都做不下去了 再次谢谢大家 glassfiber 发表于 2009 12 1 13 22 i s 本帖最后由 glassfiber 于 2009 12 1 13 24 编辑 i 有个笨办法 画一条线 然后复制一下 再把原来的线旋转 60 度 粘贴原来的线出来 这样就有两根线 了 再把两个末端连起来 三个线同时选中后强制生成实体 就是一个正三角形 或者 画一个旋转 60 度的正方形 一个不旋转的正方形 然后里面就有一个正三角形的形状了 用线连 起来就是一个正三角形 kobe010525 发表于 2010 1 7 10 29 可以用 matlab 先编一个程序 导入就行了吧 yiri 发表于 2010 4 11 12 00 能说说如何用 matlab 编 在导入呢 还是全部在脚本中实现 Lynn An 发表于 2010 4 12 10 15 就是使用 comsol with matlab 这个时候会先打开 matlab 然后打开 comsol comsol 里面选 2D 然后确定就可以了 这个时候 comsol 进入工作界面 最小化它就可以了 在 matlab 界面写程序 画正三角形什么的都有现成的命令 这个随便 google 就知道了 有很多方法 画好了以后 使用一下 geomcoerce 命令 把曲线几何强制成实体 然后在 comsol 界面写 选 文件 导入 导入几何对象 这个时候你就能在 comsol 中看见 matlab 工作 空间中的几何对象 导入就可以了 Lynn An 发表于 2010 4 12 10 17 其实说实在的 我觉得 2 楼说的画线再复制旋转的方法更容易 正三角形而已 不至于用脚本 呵呵 因为基模的损耗一般最小 对应于虚部也最小 因为基模的损耗一般最小 对应于虚部也最小 求光纤传播模式和损耗的问题求光纤传播模式和损耗的问题 1 那个那个 in plane wave 是面内波 不是光学里面说的平面波 是面内波 不是光学里面说的平面波 plain wave comsol 里面的垂直里面的垂直 波也可以是平面波 面内波和垂直波 只是描述的光的入射方向 你计算的模式的话 画光纤截面多好 波也可以是平面波 面内波和垂直波 只是描述的光的入射方向 你计算的模式的话 画光纤截面多好 因为模式分布是一个稳态光场分布 你弄个面内波做一段无意义的传输干什么呢 唉 被因为模式分布是一个稳态光场分布 你弄个面内波做一段无意义的传输干什么呢 唉 被 fdtd 毒害太深毒害太深 了 了 2 求解模场使用的是特征频率分析这种求解器 你要告诉他一个大概的寻找范围 比如你做个光纤仿真 求解模场使用的是特征频率分析这种求解器 你要告诉他一个大概的寻找范围 比如你做个光纤仿真 芯的折射率芯的折射率 1 6 包层折射率 包层折射率 1 5 你就可以告诉软件从你就可以告诉软件从 1 6 附近开始找 因为模式的折射率一定在附近开始找 因为模式的折射率一定在 n1 和和 n2 之间 你要是告诉软件从折射率之间 你要是告诉软件从折射率 1 附近找 他肯定算不出来 附近找 他肯定算不出来 3 步骤就是画光纤截面 用垂直波模场分析这种应用模式 指定波长 指定各求解域的材料属性 外围 步骤就是画光纤截面 用垂直波模场分析这种应用模式 指定波长 指定各求解域的材料属性 外围 边界用吸收边界 散射边界或者边界用吸收边界 散射边界或者 PML 内部边界不用管它 默认连续 然后告诉他在哪个折射率附近 内部边界不用管它 默认连续 然后告诉他在哪个折射率附近 找 然后就计算 在后处理中看功率分布找到基膜 然后记下软件给出的有效折射率 是个复数 虚部表找 然后就计算 在后处理中看功率分布找到基膜 然后记下软件给出的有效折射率 是个复数 虚部表 示了损耗 有公式的 翻翻论坛上的帖子 示了损耗 有公式的 翻翻论坛上的帖子 散射边界和散射边界和 PML 都属于吸收边界 都属于吸收边界 区别就是 散射边界是个一阶的吸收边界 而区别就是 散射边界是个一阶的吸收边界 而 PML 是二阶的 是二阶的 换句话说 如果波矢精确已知 散射边界就是完美吸收边界 没有反射 但是如果波矢无法给出 那散射换句话说 如果波矢精确已知 散射边界就是完美吸收边界 没有反射 但是如果波矢无法给出 那散射 边界就不可能完美吸收了 会有反射存在 边界就不可能完美吸收了 会有反射存在 PML 要高级一点 那可以不要求用户指定波矢 反正波来了就吸收 衰减掉 只要要高级一点 那可以不要求用户指定波矢 反正波来了就吸收 衰减掉 只要 PML 层够厚 层够厚 i m A y s 但是大多数都是伪模 光的能量都不在中心 而是在各种各样乱七八糟的区域 200 里面可能有一两个 是限制在芯区的传输模 v Q e y B n D 如果光子晶体柱半径再小的话 可能一个传输模也找不到 wd19830203 发表于 2009 11 17 20 43 原因就是没找到模式 具体应该参照第一幅图的折射率设定 然后都找几个模 应该可以找到正确的模式 我目前还在做 希望和大家多交流 beautycatcher 发表于 2009 11 24 08 32 楼主 最后问题出在哪啊 给个总结贝 xwx000000 发表于 2009 11 30 13 42 quote 原因就是没找到模式 具体应该参照第一幅图的折射率设定 然后都找几个模 应该可以找到正 确的模式 我目前还在做 希望和大家多交流 size 2 color 999999 wd19830203 发表于 2009 11 17 20 43 color url 你根据成品的 neff 在这个值附近把求解的本征值个数设多一点看看 yeyong0 发表于 2007 12 30 10 38 回复回复 3 shanyrain 的帖子的帖子 谢谢 有空向你请教 我来增加点本征值个数试试 yeyong0 发表于 2008 1 3 15 20 回复回复 3 shanyrain 的帖子的帖子 感谢番茄 按你的指导我已经找到模了 图很漂亮 有问题再向你询问 谢谢 yeyong0 发表于 2008 1 6 11 03 COMSOL 仿真仿真 3D 的光子晶体光纤时网格不能划分 是什么原因啊 的光子晶体光纤时网格不能划分 是什么原因啊 COMSOL 仿真 3D 的光子晶体光纤时网格不能划分 是什么原因啊 弹出的对话框是这样的 f J f Zd 9Y J Z O S 这个网格我粗化了也不行 mahui 发表于 2008 1 6 14 18 你的模型不是很规则吧 手动划分实施 shanyrain 发表于 2008 1 7 10 32 有没有试试从 2D 来拉伸 建立 3D COMSOLCOMSOL MultiphysicsMultiphysi