基于能带特性与变换介质的液体表面波调控研究.docVIP

　　基于能带特性与变换介质的液体表面波调控研究 第1章绪论 1.1研究背景与意义 波动是质点群联合起来表现出的周而复始的运动现象。水波广泛存在于自然界中，是最为直观的波动现象（如图1.1所示），并影响着人类生活的方方面面。海洋中的波动问题一直是研究的热点，根据其物理机制可按表1.1进行分类，巨大的波浪能量既可能造成灾害，也可以作为可再生能源被收集利用。一般在海边看到的波浪是由风力引起的，吹风区内称之为风浪，往往波峰尖削，在海面上的分布很不规律，波峰线短、周期小，当风大时常常出现破碎现象。风浪离开风吹的区域后称为涌浪，涌浪的波面比较平坦、光滑，波峰线长，周期、波长都比较大，在海上的传播比较规则[2]。地震津波（又称海喻）是由海底地震、火山爆发、海底滑坡或气象变化产生的破坏性海浪，其波速高达每小时700?800公里，波长可达数百公里，可以传播几千公里而能量损失很小。当海喻到迖海岸浅水地带时，波长减短而波高急剧增高，可形成高达数十米的水墙，其含有的巨大能量具有极大的破坏力上述波浪均发生在空气与水的界面上，且其回复力主要是重力，可统称为表面重力波，简称表面波。对水表面波的传播现象和调控机制进行研究，不仅有助于拓展和加深对波动问题的理论认识，而且在波浪灾害防护、波浪能源利用等方面也有着潜在的工程价值。由于海洋中的水波波长往往达到几十米甚至几百米，现场测试研究一般技术复杂、耗资较大。不过，基于表面波的方程满足标度不变性，可以通过研究小波长尺度的液体表面波讨论其传播现象、调控机理和调控效果。作为波动现象的一种，液体表面波与其他经典波（如电磁波与声波）在很多方面都具有相似性，因此可以借鉴这些波的调控理论和方法，来指导液体表面波调控研究。 . 1.2研究现状 液体表面波在周期结构中的传播行为是一个经典问题[87-88]。在早期的研究中，一些力学和工程领域的专家就采用周期结构来实现较强的多重散射，实现了局部方向的表面波带隙。他们的关注点在于周期结构对水波的阻挡耗散，利用带隙频段上对水波的强反射，可以对河口、海岸、海上平台等设施起到防护作用。他们的研究主要是关于水波在二维圆柱阵列中的绕射传播[89-91]以及周期圆柱阵列[92]的理论计算方面，只研究了柱体填充比较小的情况，仅发现了其中的方向带隙。此外，Chou研究了周期流体存在的表面波带隙的不稳定性[93]以及在周期交界面上传播的液体表面波的带隙现象随着以光子晶体为代表的超构材料研究的不断发展，一些物理学者也开展了液体表面波在周期结构中传播的研究，利用液体表面波来展示波动现象，研究或证实新的物理机制。1999年，Torres等人利用周期挖孔阵列，首次将Bloch波可视化(见图1.11)。在引入无序后，看到了两个畴之间的边界（domain 和5 cm时液体表面波频率与波长的关系曲线。从图中可以清晰地看到，液体表面波的频率与波长的对应关系为非线性。同时也可以注意到液体的深度与毛细长度对液体表面波波长及频率的影响。当波长较小时，液体的毛细长度对频率的影响较为明显，而当波长较大时，表面波频率主要与液体深度有关。自从MSE问世到现在，国内外众多专家学者围绕其做了大量的研究，不断有人针对更大的坡度情形，提出相应的波动方程，如Kirby提出的在缓变地形上考虑速度变化的振幅分量的护展型缓坡方程（EMSE) [126]，Chamberlain和Porter提出的只包含底坡作用高阶项而不包括无限小量的改进的缓坡方程（MMSE) Agnon提出的考虑水深偏导中所有线性项的增缓坡方程（AMSE) 等。针对变化底面深度上液体表面波中的布拉格散射问题，Ye于2003年也给出了 一个处理方程,在近似程度上与MSE样都是一阶近似，而且该方程在数学形式上与式（2.38)的类似性更为明显。 .. 2.2液体表面波在周期结构中传播时的能带和场的计算方法 传输矩阵法（transfer matrix method,简称TMM)是广泛应用于计算光子晶体与声子晶体[42，132_135]中一维多层体系的能带与传播特性的方法。其基本思想是从连续状态参数的基本方程入手，结合界面连续条件，得出单个周期的传递矩阵。通过引入周期边界条件求解得到相应的色散关系，即能带结构。通过有限个传输矩阵相乘可以得到有限周期结构的传播特性[42]。类似的，对于液体表面波，一般用传输矩阵法来处理一维台阶体系的色散关系和传输系数。对于一维体系，其对应波动方程简单，可以通过传递矩阵求得解析的能带结构。而在二维体系中方程比较复杂，必须借助数值法求解，平面波展开法（plane ethod,简称PWE)是在求解光子晶体与声子晶体的二维周期体系中最常用的方法之一。其基本思想是：利用结构的周期性，把结构参数按傅里汁级数展开，并与布洛赫定理结合，将波动方程在倒格矢空间以平面波叠加的形式