本篇論文使用Dirichlet to Neumann法來計算彈性波在由多殼層介質柱子與流體背景構成的聲子晶體之頻帶結構。傳統常用之平面波展開法 (plane wave expension) 在計算聲子晶體之頻帶時，若遇到固體與流體共存之系統，會由於固體區與流體區之波函數分量數量不同，而遭遇無法執行的困難。本論文使用之Dirichlet-to-Neumann法能夠輕易的藉由連接邊界條件來解決上述困難，並且大幅縮小計算資源。第二章首先介紹二維基本結構以及彈性波之基本理論，然後推導含有剪力波之彈性固體柱在在流體背景中之聲子晶體頻帶計算細節，以盡可能模擬真實情況。在第三章中，我們將原本的Dirichlet to Neumann法推廣至三維，並提出了可能的執行方法，以計算三維聲子晶體之頻帶。最後探討在二維結構中使用不同的殼層材料與厚度搭配對於聲波在週期介質中傳播情形之影響。

In this thesis the Dirichlet- to-Neumann map method is applied to the calculation of the band structure of the phononic crystal which consists of multi-shell rods and fluid background. Due to the inconsistency of the total number of dynamical wave variables in the solid and fluid regions of the phononic crystal, the conventional plane wave expansion method is not directly applicable. We thus turn to use the Dirichlet-to-Neumann map method to overcome such a problem. In chapter 2, we first introduce the typical structure of the phononic crystal and the basic theory of elastic waves, and then explain how to derive the wave fields in every region in a unit cell and give the calculation details for the band structures of the phononic crystals in the text and in the appendix. In order to mimic the real acoustic/elastic systems, in most calculations the shear waves are taken into account. In chapter 3, we then generalize the original Dirichlet-to-Neumann map method to 3D and propose a possible procedure for its implementation to the band structure calculations of 3D phononic crystals. Finally, we discuss the effects of the materials used in the phononic crystal and the thicknesses of the shells on the resultant band structures.

[1] C. Kittle, Introduction to Solid State Physics, 7th Ed., John Wiley & Sons (2001)
[2] David K. Cheng, Field and wave electromagnetic, 2th Ed., Addison-Wesley
(1989)
[3] Zhengyou Liu, Xixiang Zhang, Yiwei Mao, Y. Y. Zhu,Zhiyu Yang, C. T. Chan,
Ping Sheng, Science, 289, 1734 (2000)
[4] Ping Sheng, X. X. Zhang, Z. Liu, C. T. Chan, Physica B. Condensed Matter, 338,
201 (2003)
[5] H.H. Huang a, C.T. Sun , G.L. Huang, International Journal of Engineering
Science, 47, 610 (2009)
[6] Wu Fu-Gen, Liu Zheng-You, Liu You-Yan, Chinese Phys. Lett., 18, 785 (2001)
[7] I. E. Psarobas, A. Modinos, R. Sainidou, N. Stefanou , Phys. Rev. B, 65, 064307
(2002)
[8] Ph. Lambin, A. Khelif, J. O. Vasseur, L. Dobrzynski, B. Djafari-Rouhani, Phys.
Rev. E, 63, 066605 (2001)
[9] Suxia Yang, J. H. Page, Zhengyou Liu, M. L. Cowan1, C. T. Chan, Ping Sheng,
Phys. Rev. Lett., 88, 104301 (2002)
[10] Liu ZY, Chan CT, Sheng P, Phys. Rev. B, 65, 165116 (2002)
[11] Zalipaev VV, Movchan AB, Poulton CG, et al., Proceedings of the Royal
Society of London series A-Mathematical Physical and Engineering Sciences,
458, 1887 (2002)
[12] Xin Zhang, Zhengyou Liub, Youyan Liua, Fugen Wu, Phys. Lett. A, 313, 455
(2003)
[13] Jensen JS, Journal of Sound and Vibration, 266, 1053 (2003)
[14] Hou ZL, Fu XJ, Liu YY, Phys. Lett. A, 317, 127 (2003)
[15] Wu FG, Liu ZY, Liu YY, Journal of Physics D-Applied Physics, 35, 162 (2002)
103
[16] Fugen Wu, Zhilin Hou, Zhengyou Liu, Youyan Liu, Solid State
Communications, 123, 239 (2002)
[17] Li XL, Wu FG, Hu HF, et al., Journal of Physics D-Applied Physics, 36, 15
(2003)
[18] Wu FG, Hou ZL, Liu ZY, et al., Phys. Lett. A, 292, 198 (2001)
[19] Wu FG, Liu ZY, Liu YY, Chinese Physics Letters, 18, 785 (2001)
[20] Sigmund O, Jensen JS, Philosophical Transactions of the Royal Society of
London Series A-Mathematical Physical and Engineering Sciences, 361, 1001
(2003)
[21] Cervera F, Sanchis L, Sanchez-Perez J V, Martinez-Sala R, Rubio C, Meseguer F,
Lopez C, Caballero D, Sanchez-Dehesa, J Phys. Rev. Lett., 88, 023902 (2001)
[22] Kuo C-H, Ye Z, J. Phys. D: Appl. Phys., 37, 2155 (2004)
[23] Kushwaha M S, Halevi P, Martinez G, Dobrzynski L, Djafari-Rouhani B, Phys.
Rev. B, 49, 2313 (1994)
[24] I. E. Psarobas,,N. Stefanou ,A. Modinos, Phys. Rev. B, 62, 278 (2000)
[25] S. Guo and S. Albin et. al. “Simple plane wave implementation for photonic crystal
calculations” Opt. Express Vol. 11, pp. 167-175, 2003.
[26] Dennis M. Sullivan, Electromagnetic Simulation Using The FDTD Method,
Wiley-IEEE Press, New York, 2001.
[27] Kane S.Yee et. al. “Numerical Solution of Initial Boundary Value Problems
Involving Maxwell’s Equations in Isotropic Media” , IEEE. Trans. Antennas.
Propag. ,Vol.14, pp. 302-307, 1966.
[28] O. C. Zienkiewicz, Robert Leroy Taylor, J. Z. Zhu, The Finite Element Method: Its Basis and Fundamentals, Butterworth-Heinemann press, 2005.
[29] Rahman, B.M.A.; “A review on the characterization of photonic devices using the finite element method,” Electrotechnical Conference, 1996. MELECON ‘96., 8th
Mediterranean , vol.2, no., pp.705-708 vol.2, 13-16 May, 1996.
[30] Jianhua Yuan, Ya Yan Lu, J. Opt. Soc. Am. A, 23, 3217 (2006)
[31] 欒丕綱、陳啟昌，光子晶體：從蝴蝶翅膀到奈米光子學，五南出版社，台北市，民國九十九年。
[32] Arfken G B, Weber H J, Harris F, Mathematical Methods for Physicists, 6th Ed.,Academic Press (2005)
[33]Murray R. Spiegel ,Theory and problems of advanced mathematics for engineers & scientists. McGRAW-HILL (1983)
[34] Feng-Lian Li Yue-Sheng Wang and Chuanzeng Zhang, Phys. Scr. 84 (2011) 055402 (9pp)