删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

二维近零折射率声学材料的负向Schoch位移

本站小编 Free考研考试/2021-12-29

摘要:对二维声学超材料与常规材料界面处的负向Schoch位移进行了研究. 研究表明, 在界面处增加合适厚度的覆盖层, 当声波在某一频率附近从常规材料向近零折射率材料传播时, 反射声波相对于入射声波会在界面处发生负向Schoch位移, 此时, 超材料的有效体积模量倒数值趋近于零从而使它成为一种折射率近零的声学材料; 同时, 超材料有效阻抗的极大值和反射系数的极大值都在这一特殊频率处且反射系数虚部相位在对应频率处有π rad的相位突变; 研究还发现, 发生负向Schoch位移的频率位于MK方向第一带隙中且靠近上边界的通带频率. 常规材料界面处的Schoch位移通常是正向的且大小可忽略不计的, 本文利用近零折射率声学材料实现了负向的Schoch位移, 为设计出基于界面声波的声学器件提供了一种新的理论参考.
关键词: 近零折射率材料/
有效参数/
Schoch位移

English Abstract


--> --> -->
当一束声波入射到两种介质界面并发生反射时, 实际的反射波束相对于几何预测的反射波束有一定的正向(或负向)偏移, 这个正向(或负向)的偏移称为Schoch位移, 这种现象称为Schoch效应[1]. Schoch效应最初是由Schoch于1950年通过纹影摄影技术观察超声波在界面上的反射波束时, 偶然发现的一种非镜面反射现象[1,2], 而后Bertoni和Tamir[3]发展了一种统一的近似理论来解释这一现象. 自2000年以来, 科学家开始聚焦于可调的周期性结构材料的Schoch效应. 2003年, Declercq等[4]基于非均匀波衍射理论利用一维周期性突起的黄铜-水结构探究超声波的反射 时, 获得了与实验结果一致的仿真数据. 2004年, Declercq等[5]将非均匀波衍射理论拓展到负向Schoch位移的声学衍射, 研究结果表明Schoch效应与界面处类斯通利波的产生密切相关. 2009年, Herbison等[6]采用脉冲波代替时域谐波的新方法对Schoch位移进行定量测量. 2014年, Declercq[7]实验研究了产生Schoch位移时超声波的能量转换. 2017年, 陈宗旺等[8]通过在常规材料与声子晶体的界面处引入缺陷(覆盖层), 通过调节覆盖层的厚度实现了巨大的正向Schoch位移.
超材料因其具有常规材料不具备的超常物理性质而备受关注[9-11], 声波的负向Schoch位移就是其中的一种. 在未来声学器件趋于小型化的前提下Schoch位移是一个必须考虑的问题. 近年来对于近零折射率材料的研究方法无论是声学还是光学方面都已经相当成熟, Schoch效应是受到古斯-汉欣效应的启发发现的, 目前对于古斯-汉欣效应的研究无论是在理论上还是高精度的实验器件上都已取得显著成果[12-16]. 文献[12]基于古斯-汉欣位移提出了一种双通道窄带滤波器. 文献[13]理论推导并仿真验证了线性偏振光从空气倾斜入射到近零折射率光学材料的情况. 文献[14]研究了不同的入射波波长以及温度对近零折射率光学材料的古斯-汉欣位移的影响. 文献[15]基于在近零折射率光学材料表面覆盖可调电压的单层石墨烯, 实现了古斯-汉欣位移大小可调的机制. 文献[16]通过在常规材料与光学超材料之间加一厚度可调的覆盖层, 以光束的入射角和覆盖层厚度为变量, 在厚度选取合适的条件下实现了超胞尺寸数倍的巨大负向古斯汉欣位移. 但是目前在特性更为复杂的近零折射率声学材料中并没有见到关于声波Schoch效应的相关报道, 为了进一步探究声学的Schoch效应, 基于文献[16]中对光学超材料的研究方法和实验结论, 本文设计了一种模型, 通过调控散射体橡胶圆柱半径以及覆盖层水的厚度, 有效实现了不同参数下伴随着不同物理特性的近零折射率声学超材料的负向Schoch位移. 在此基础上又探究了超材料中Schoch效应与其带结构之间的联系.
光波(电磁波)在光学超材料中传播时, 以介质的介电常数和磁导率表征光波与光学超材料的相互作用. 而声波在声学超材料中传播时, 以周期性复合介质的声学参数(即质量密度和体积模量)表征声波与声学超材料的相互作用, 但不同的是对于声学超材料而言, 还可以用剪切模量等参数进行调控, 这就决定了声学超材料具有更加复杂、更加有趣的物理特性. 光学超材料中的古斯-汉欣位移通常发生在兆赫兹及以上的频率区间, 其位移尺寸对现代光学器件精度的影响可以忽略. 相对的声学超材料的Schoch位移通常发生在低频领域, 对于工作频率处于低频领域的中小型声学器件精度的影响必不可忽略. 而本文声波的入射频率选取范围为400—1000 Hz, 负向Schoch位移的尺寸最大达到了晶格常数的8倍, 因此对近零折射率声学超材料的研究将有更丰富的内涵.
本文使用的模型系统由常规材料与声学超材料组成, 其具体结构如图1(a)所示, 从上到下由水银(mercury)、水(water)和声学超材料依次排列组成. 声学超材料由二维的橡胶圆柱散射体以三角晶格的形式周期性排列在水中构成. 图1(b)为声学超材料的单个原胞, 其中橡胶圆柱散射体半径为R, 间距为a, 基体为水. 本文采用有限元软件COMSOL Multiphysics模拟声波的Schoch效应, 数值模拟的材料参数为: 散射体橡胶的密度ρ = 1.3 × 103 kg/m3, 拉姆常数λ = 2.1 × 108 N/m2, μ = 5.2 × 107 N/m2; 水的密度ρ0 = 1.0 × 103 kg/m3, 波速c0 = 1490 m/s; 水银的密度ρ1 = 13.6 × 103 kg/m3, 声速c1 = 1400 m/s. 为了更易激发出声表面波和漏波, 在超材料与水银溶液的界面处加一材料为水的覆盖层, 覆盖层的厚度选定为${h_0} = 0.3 a$. 以高斯声波(半腰宽$w=5 a$)从水银溶液中倾斜45°方向入射到声学超材料中, 以入射声波的频率f为变量进行仿真测试.
图 1 (a) 模型结构, 从上到下分别为水银、水层和声学超材料; (b) 超材料的单个晶胞, 其中圆柱散射体为橡胶, 基体为水
Figure1. (a) Model structure, which consists of mercury, water layer and acoustic metamaterial from top to bottom; (b) single crystal cell of metamaterial, where the cylindrical scatterer is rubber and the matrix is water.

当入射声波的约化频率为F = fa/c0 = 0.456 = FS (伴有最大负向Schoch位移的约化频率记为FS), 反射声波具有明显的负向Schoch位移和较好的高斯波形, 如图2(a)所示, 其中黑色实线箭头代表声波真实的传播方向, 虚线箭头代表几何预测的方向, 在虚线箭头右侧还有一微弱的反射波束, 所以图中共有两反射波束. 反射过程中能流较多的一部分以表面波和漏波的形式分别在水银与水的界面和覆盖层中传播, 较少的一部分用于波束的正常反射. 左边的反射波束位置明显偏离了声波的入射点, 向左偏离了入射波束中心位置$8 a$, 此反射波束就是由表面波和漏波主导的负向反射[17]. 同时, 在两反射波束之间出现了空白区域, 这是由反射波束和出现在水层和超材料界面附近的漏波共同作用形成的[18,19]. 如图3(a)所示, 当有负向Schoch位移出现时, 在水层和超材料界面附近会形成大量的后向漏波. 改变F = 0.503作为对照, 没有Schoch位移出现其声压场图如图3(b), 此时在水层和超材料界面附近完全没有漏波出现. 研究还发现在这一结构模型下, 约化频率F仅在0.447—0.466范围内有负向Schoch位移出现, 如图2(b)所示, 并且在这一范围内负向Schoch位移随着F的增大先增大后减小为零, 在F = 0.456处位移取得最大负值$8 a$.
图 2 (a) 具有明显负向Schoch位移的声压场图($R = $$ 0.224 a$, F = fa/c0 = 0.456 = FS, a = 1 m), 其中纵轴和横轴分别表示本文结构的高度y和宽度x, 右侧的颜色条对应的物理量为总声压场强Pt; (b) 在FS = 0.456附近, 负向Schoch位移随约化频率F的变化
Figure2. (a) Acoustic pressure field with significant negative Schoch displacement ($R = 0.224 a$, F = fa/c0 = 0.456 = FS, a = 1 m); (b) relation of Schoch displacement to reduced frequency F near FS = 0.456.

图 3 (a) $F = 0.456$, 水层和超材料界面附近出现大量后向漏波; (b)$F = 0.{{503}}$, 水层和超材料界面附近完全没有漏波
Figure3. (a) $F = 0.456$, there are a large number of backward leaky Rayleigh waves appear near the interface between the water layer and the metamaterial; (b)$F = $$ 0.{{503}}$, there is no leaky Rayleigh wave near the water layer and the metamaterial interface.

为了探究伴随着这种负向Schoch位移出现时超材料的物理特性, 参照周期性复合流固材料的有效参数[20,21], 计算并获取了超材料不同的物理参数随约化频率F的变化关系.
图4(a)显示当负向Schoch位移发生时, 超材料的reffZeff极大值出现在同一约化频率$F = $$ 0.456$(标记为A点), 这与图2(a)发生负向Schoch位移时FS的值相同; 图4(b)显示neff, ρeff, 1/κeffF的变化关系, 图中添加两条平行坐标轴的虚线作为辅助线, 其交点记为B (0.456, 0), B点的横坐标取值与A点处的F相同, 图像右下角内嵌图形为B点附近的放大. 以B作为参考点, 在F = 0.456时, ρeff小于零, 1/κeff趋向于零, neff的值小于且接近零, 它们之间满足关系[22] $\mathop n\nolimits_{{\rm{eff}}}^2 =\! 1/{\kappa _{{\rm{eff}}}} \cdot\! {\rho _{{\rm{eff}}}}$, 此时超材料为近零折射率声学材料.
图 4 声学超材料不同的物理参数特性 ($R = 0.224 a$, FS = 0.456) (a) 相对阻抗(Zeff)和反射系数(reff)随F的变化; (b) 有效折射率(neff)、有效质量密度(ρeff)、有效体积模量的倒数(1/κeff)随F的变化, 图像右下角内嵌图形为B点附近的放大
Figure4. Different physical parameters of acoustic metamaterials ($R = 0.224 a$, FS = 0.456): (a) Relationship between relative impedance (Zeff) and reflection coefficient(reff) with F; (b) relationship of effective refractive index (neff)、effective mass density (ρeff) and inverse of effective volume modulus (1/κeff) to F, the embedded figure in the lower right corner of the image is an enlargement near point B.

古斯-汉欣位移主要是由反射波束的相位突变导致[23], Schoch位移是否也与此有关?
图5$R = 0.224 a$, FS = 0.456条件下发生Schoch位移时超材料的另外两个物理特性. 图5(a)表示反射系数的虚部相位发生突变, 突变前后的相位差值为π rad. 图5(b)为是偶极子激发形成的本征模的声压场图的三分之一, 根据超材料形成的原理[24]知, 偶极子表征超材料有效质量密度为负, 与图4(b)ρeff < 0相符.
图 5 (a) 反射系数虚部相位在$F = 0.456$处发生突变($R = $$ 0.224 a$, FS = 0.456); (b) 晶格的本征模声压场图
Figure5. (a) Phase of the imaginary part of the reflection coefficient mutates at $F = 0.456$ ($R = 0.224 a$, FS = 0.456); (b) eigen mode acoustic pressure field diagram of the lattice.

为了探究图4图5中的物理参数特性与负向Schoch位移的联系, 改变散射体R, 寻找反射系数峰值对应的F, 当图4图5中的物理参数特性同时得到满足时, 验证该约化频率下是否出现负向Schoch位移. 经过多次的仿真发现, 改变散射体R, 物理参数特性与$R = 0.2{{24}}a$时的结果类似. 发生负向Schoch位移时散射体半径RFS(fa/c0)的变化关系以图6(a)表示, 整体来看约化频率FS与散射体橡胶半径R近似呈现线性关系, 随着R的增大FS逐渐减小.
图 6 (a) 发生Schoch位移时FSR的变化; (b) 负向Schoch位移随覆盖层厚度${h_0}$的变化($R = 0.224 a$, FS = 0.456); (c) 负向Schoch位移随入射声波半腰宽wi的变化; (d) 反射声波半腰宽wr随覆盖层厚度的变化
Figure6. (a) Relationship of FS to R when the Schoch displacement happens; (b) relationship between negative Schoch displacement and overburden thickness ${h_0}$; (c) relationship between negative Schoch displacement and the half-waist width of the incident acoustic wave; (d) relationship between the half-width of the reflected acoustic wave and the thickness of the covering layer.

常规材料与光学超材料界面处合适的覆盖层厚度能极大激发表面波和漏波[16,17,25], 在常规材料与声学超材料界面处增加覆盖层(水层)存在类似现象. 经过计算发现, 界面处覆盖层厚度对负向Schoch十分重要, 如图6(b)所示. 只有当覆盖层厚度处于某一区间时才有不同程度的负向位移出现. 当${h_0}$处于0.2a—0.4a时, 负向位移值随着${h_0}$的增大先增大后减小, 在${h_0} = 0.3a$处取得最大负值${{8}}a$. 这是因为随着${h_0}$的增大, 界面处产生表面波和漏波的条件[16]${h_0}$处于${{0}}{{.3}}a$附近近似得到满足, 进而在界面处诱导出巨大的负向位移. 但当${h_0}$较大或较小时, 不满足界面处激发表面波和漏波的条件, 反射声束的位移可忽略. 为了衡量负向Schoch位移的大小, 进一步探究了声波半腰宽w与负向Schoch位移的关系如图6(c)6(d)所示. 随着入射声波半腰宽wi的增大, 负向Schoch位移的变化情况如图6(c)所示. 图6(d)为反射声波半腰宽wr随覆盖层厚度${h_0}$的变化关系, 随着覆盖层厚度的增大反射声波半腰宽wr会先增大后减小. 对比图6(b)图6(d), 在入射声波半腰宽wi${{5}}a$的条件下, 当反射声波半腰宽wr最大时负向Schoch位移达到最大值.
在带隙边缘附近, 光子晶体对光束的调控与有效介电常数接近零的材料类似, 因此带隙附近的超材料通常具有许多新颖的光学特性, 这些特性统称为超折射光学[26], 其主要现象为入射光束的分裂或放大, 或负的古斯-汉欣位移[27]. 消逝场的存在会影响波的传播, 最明显的情况就是波在带隙中的传播[28], 不仅古斯-汉欣效应与带隙有着密切的关系[29], Schoch效应与超材料的带隙同样有着联系. 在全内反射情况下, 不同的介质界面发生Schoch效应时通常伴随有表面波或漏波的产生, 消逝波存在导致了反射波束相位发生突变, 使得反射波束与入射波束之间有一个较大的相位差. 研究发现, 发生负向Schoch位移的约化频率总是位于第一带隙上边界附近.
由前面知$R = 0.2{{24}}a$, 发生负向Schoch位移时FS = 0.456, 其带结构如图7(a)所示, 红色实线是添加的纵坐标为F = FS = 0.456的辅助线. 从图上看, 第一带隙上边界与红色实线有两个交点, 左右两交点分为位于ΓMMK方向, 记右边交点为C点. 因为声波入射角为45°, 导致在ΓM方向激发较复杂的模式[30], 故仅需探究C点附近的情况就能说明问题. 多次的仿真发现, 若入射声波的约化频率在C点及以下0.002个单位的约化频率范围内, 声压场图显示有负向Schoch位移的出现, C点以上则没有. 负向的古斯-汉欣位移需要在后向漏波结构(或负折射率材料)中激发[31,32], 同样可认为负向的Schoch位移也需要在后向漏波结构(或有效折射率为负的声学超材料)中激发. 对于带隙下边界, 由图7(a)知在MK方向, 带隙下边界的F值在0.309附近波动. 为了分析带隙下边界的情况, 图7(b)给出了F = 0.309时超材料的有效参数, 图中作有两条平行于坐标轴的黑色辅助虚线(neff = 0和F = 0.309), 其交点为B, B点的有效参数即为声波以F = 0.309入射时超材料的有效参数. 从图7(b)中的B点知, 此时1/κeff不再趋向于零, neff大于零, 因此推测声波以带隙下边界频率F = 0.309入射时并不能激发负向的Schoch位移. 图7(c)所示的声压场图验证了本文的推测, 当F = 0.309时, 水层和超材料界面附近并不能激发后向漏波, 没有负向Schoch位移出现. 所以不再探究带隙下边界的情况, R为其他值时情况与此类似.
图 7 (a) 超材料的带结构($R = 0.2{{24}}a$); (b) neff, ρeff, 1/κeffF的变化; (c) 声波以带隙下边界频率($F{{ = 0}}{{.309}}$)入射时超材料的声压场图
Figure7. (a) Band structures of metamaterial ($R = 0.2{{24}}a$); (b) relationship of neff, ρeff and 1/κeff to F; (c) acoustic pressure field diagram of metamaterial when acoustic waves incident at the lower band-gap boundary frequency.

上面的结论表明, 负向Schoch位移发生的约化频率位于第一带隙之中且靠近上边的通带; 不同系统的约化频率范围稍有不同.
本文主要研究了二维近零折射率声学材料发生负向Schoch位移时的物理参数特性以及模型结构参数对位移的影响. 通过在常规材料与声学超材料界面处增加材料为水的覆盖层, 选取合适的厚度参数, 实现了巨大的负向Schoch位移. 计算结果表明, 合适的覆盖层厚度能够诱导界面处反射波束的部分能流, 使其进入覆盖层并且转化为向后传播的漏波, 进而在低频领域获得波长数量级(米尺寸级别)的负向Schoch位移. 大的负向Schoch位移的发生会伴随着以下几个物理特性, 超材料的有效体积模量倒数值趋近于零从而使它成为一种折射率近零的声学材料; 超材料有效阻抗的极大值和反射系数的极大值几乎在同一约化频率处; 反射系数虚部相位在对应约化频率处有π rad的相位突变; 发生负向Schoch位移的约化频率位于MK方向第一带隙中且靠近上边界的通带频率. 常规材料界面处声波的Schoch位移通常是正向的且大小可忽略的, 本文获得了负向的Schoch位移, 不但丰富了声波Schoch现象方面的研究成果, 而且为设计基于界面声波的声学器件提供理论参考.
相关话题/材料 声学 结构 物理 光学

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 纳米结构及浸润性对液滴润湿行为的影响
    摘要:液滴在纳米结构表面的润湿模式研究(Dewetting,Cassie,PartialWenzel及Wenzel)对强化冷凝、表面自清洁、油水分离等具有重要意义,现有文献主要研究了液滴在微柱阵列纳米结构表面的润湿行为.本文采用分子动力学模拟,研究了纳米结构倾角及表面浸润性对氩液滴在铂固体壁面上润湿 ...
    本站小编 Free考研考试 2021-12-29
  • 高效硫硒化锑薄膜太阳电池中的渐变能隙结构
    摘要:硫硒化锑(Sb2(S,Se)3)薄膜太阳电池因其原材料丰富、制备方法简单、性能稳定等优势近年来得到了快速发展.本文基于Sb2(S,Se)3吸光层能隙可调的特点,应用wx-AMPS软件对具有渐变能隙Sb2(S,Se)3太阳电池进行建模仿真和结构设计,并与50%Se含量的恒定能隙Sb2(S,Se) ...
    本站小编 Free考研考试 2021-12-29
  • 基于微结构气体探测器对单能和连续谱快中子的模拟解谱
    摘要:本工作探索了基于TripleGEM探测器对快中子能谱的测量,利用MCNPX和Geant4软件分别模拟了两种在TripleGEM阴极耦合由多层聚乙烯组成的堆栈式中子转化质子的转化模型,研究对象包含了5种单能中子源和Am-Be连续谱中子源.模拟得到了探测系统对160条单能中子的响应函数和上述源的反 ...
    本站小编 Free考研考试 2021-12-29
  • 碳纤维增强复合材料褶皱缺陷的超声成像
    摘要:碳纤维增强复合材料在制造和使用过程中会产生褶皱,褶皱缺陷的存在使得纤维层起伏不平,超声波传播至纤维界面时不同声束方向的声能量存在差异.针对此问题,提出一种基于声束指向性函数校正的全聚焦阵列成像方法:制备了含褶皱缺陷的碳纤维增强复合材料试样,建立超声阵列数据采集实验系统,捕获全矩阵数据;考虑超声 ...
    本站小编 Free考研考试 2021-12-29
  • 二维材料<i>X</i>Te<sub>2</sub> (<i>X</i> = Pd, Pt)热电性能的第一性原理计算
    摘要:利用密度泛函理论结合玻尔兹曼输运方程,预测了二维层状热电材料XTe2(X=Pd,Pt)的热电性质.两种材料都具有较低的热导率,材料的晶格热导率随温度的升高而降低,且表现出各向异性.而电子热导率随温度的升高而升高.在较低温时,晶格热导率对总热导率的贡献占据主导地位.较高的载流子迁移率、电导率及塞 ...
    本站小编 Free考研考试 2021-12-29
  • 准二维湿颗粒体系融化过程中的结构与缺陷
    摘要:研究颗粒体系中的结构与缺陷对于研究固-液融化的物理机制具有重要的价值.本文实验研究了垂直振动下单层湿颗粒在固-液融化过程中的结构与缺陷.根据实验及理论分析构建了湿颗粒体系的接触模型,量化了准二维湿颗粒体系融化过程中颗粒的结构变化.然后以颗粒为点建立Voronoi图对颗粒体系的“相”转变进行研究 ...
    本站小编 Free考研考试 2021-12-29
  • 纳米光子学结构对GaInAsSb p-n结红外光电性能的调控
    摘要:GaInAsSb在红外光电领域具有重要应用价值,但是窄带隙材料较高的本征载流子浓度和俄歇复合系数使其室温暗电流密度较高,需要进行制冷才能获得满足应用要求的光电性能.本文利用表面宽带隙半导体纳米柱阵列和背面高反射率金属对GaInAsSbp-n结有源区进行双面光调控,将光限制在较薄的有源区进行吸收 ...
    本站小编 Free考研考试 2021-12-29
  • 周期性分流微通道的结构设计及散热性能
    摘要:微通道散热器在集成电路中具有重要应用,但目前传统的长直微通道散热过程导致温度不均匀,散热效率较低.本文设计了一种周期性分流微结构并与传统微通道进行集成,实现了一种高效率的周期性分流微通道散热器.基于以上周期性分流微通道,系统研究了单根微通道内微结构数目、微结构的排布方式及结构参数对其散热性能的 ...
    本站小编 Free考研考试 2021-12-29
  • 含钨难熔高熵合金的制备、结构与性能
    摘要:高熵合金(high-entropyalloys,HEAs)作为一种新型多主元合金,原子排列有序、化学无序,具有高熵、晶格畸变、缓慢扩散、“鸡尾酒”等四大效应,表现出优异的组合性能,有望作为新型高温结构材料、耐磨性材料、抗辐照材料应用于航空航天、矿山机械、核聚变反应堆等领域.本文介绍了目前含钨H ...
    本站小编 Free考研考试 2021-12-29
  • 超电容储能电极材料的密度泛函理论研究
    摘要:由于量子限域效应和态密度的限制,石墨烯、硅烯等二维材料的量子电容在费米能级附近趋近于零.基于密度泛函理论的第一性原理研究发现,掺杂和吸附使石墨烯等二维电极材料的电子结构得以有效的调制,它促进狄拉克点附近局域电子态的形成和/或费米能级的移动,从而使量子电容得到了提高.比较Ti(Au,Ag,Cu, ...
    本站小编 Free考研考试 2021-12-29