[論文レビュー] Clusters of Bloch waves in three-dimensional periodic media

Yuri A. Godin, B. Vaĭnberg|arXiv (Cornell University)|Jul 17, 2022

Advanced Mathematical Modeling in Engineering参考文献 30被引用数 9

ひとこと要約

本稿は、漸近解析とディリクレ・トゥ・ノイマン作用素を用いて、任意の形状の小さな不純物を有する3次元周期的媒体における音響波の分散関係を厳密に導出する。波数ベクトルが特定の例外的な値をとる場合、異なる方向または周波数の波がクラスタを形成し、分散関係が一意でなくなる現象を同定する。これは、ブラッグ反射に類似した新しい非一意性現象を示している。

ABSTRACT

We consider acoustic wave propagation through a periodic array of the inclusions of arbitrary shape. The inclusion size is much smaller than the array period while the wavelength is fixed. We derive and rigorously justify the dispersion relation for general frequencies and show that there are exceptional frequencies for which the solution is a cluster of waves propagating in different directions with different frequencies so that the dispersion relation cannot be defined uniquely. Examples are provided for the spherical inclusions.

研究の動機と目的

• 3次元周期的媒体に小さな不純物を有する音響波の伝播に関する分散関係を厳密に導出すること。
• 標準的な分散関係が波のクラスタリングにより失敗する例外的な波数ベクトルを同定・特徴付けること。
• 摂動理論とディリクレ・トゥ・ノイマン作用素を用いて、漸近近似の数学的裏付けを提供すること。
• 既存の漸近的手法を一般の不純物形状および境界条件（ノイマンおよびディリクレ境界を含む）に拡張すること。
• 特定の波数ベクトルにおいて、非一意な解が存在することを示し、ブラッグ反射に類似した現象を解明すること。

提案手法

• 不純物境界での伝達条件を満たす周期的単位胞における波動問題を定式化する。
• 各不純物を囲む球面上でのディリクレ・トゥ・ノイマン（DtN）作用素を用いて、特異的摂動問題を正則問題に還元する。
• 不純物サイズパラメータ 'a' のべき級数展開を内側および外側のDtN作用素に適用する。
• 摂動理論を用いてDtN作用素の差を解析し、有効な分散関係を導出する。
• 球ベッセル関数および層電位展開を用いて漸近係数を計算する。
• 単純立方格子における球形不純物について、明示的計算により結果を検証する。

実験結果

リサーチクエスチョン

• RQ1 3次元周期的媒体におけるブロッホ波の分散関係が、どのような条件下で一意に定義されなくなるか？
• RQ2 任意の形状の小さな不純物が、周期的媒体における波動伝播およびスペクトル構造にどのように影響を与えるか？
• RQ3 例外的な波数ベクトルにおける異なる方向または周波数の波のクラスタが出現する背後にある数学的メカニズムは何か？
• RQ4 特異的摂動周期系において、波動解の漸近近似をどのように厳密に裏付けることができるか？
• RQ5 不純物の幾何学的形状と材料対比が、例外点の存在を決定づける役割を果たすか？

主な発見

• 本稿は、一般の周波数に対して有効な、3次元周期的媒体に小さな不純物を有する場合の分散関係に対する厳密な漸近展開を確立した。
• |k| = |k − m| を満たす複数のブロッホモードが共存する波数ベクトル k が例外的であることが同定され、これにより分散関係が非一意となる。
• 例外点において、解は異なる方向に進むが異なる空間周波数を持つ波のクラスタを形成する。
• 球形不純物の場合、分散関係への一次補正には、|k|² および |k| に比例する項を含む行列 M が関与し、これにより不純物の幾何学的形状と対比が反映される。
• 非一意性は、元の非摂動問題における退化性に起因し、|k| = |k − m| を満たす複数の逆格子ベクトル m が存在するためである。
• 解析により、漸近近似誤差が O(εa³ + a⁴) であることが確認され、ここで ε = (k₊ − |k|)/|k| であり、a が十分に小さい場合の収束性が保証される。