2025-11-23T13:40:16.956535

Time-harmonic scattering of plane waves from an infinite periodically inhomogeneous medium

Hu, Rathsfeld, Zhang et al.

We propose a new radiation condition for an infinite inhomogeneous two-dimensional medium which is periodic in the vertical direction and remains invariant in the horizontal direction. The classical Rayleigh-expansion radiation condition does not apply to our case, because this would require the medium to be inhomogeneous in a half plane. We utilize the Floquet theory to derive upward/downward wave modes and define radiation conditions by expansions w.r.t. these modes. The downward radiation conditions leads to a downward Dirichlet-to-Neumann map which can be used to truncate the infinite inhomogeneous domain in the vertical direction. So we prove mapping properties of the upward/downward Dirichlet-to-Neumann maps based on the asymptotic behavior of high-order wave modes. Finally, we verify the strong ellipticity of the sesquilinear form corresponding to the new scattering problem and show the unique solvability for all wavenumbers with the exception of a countable set of numbers bounded below by a small positive constant.

academic

無限周期不均質媒質からの平面波の時間調和散乱

基本情報

論文ID: 2510.14070
タイトル: Time-harmonic scattering of plane waves from an infinite periodically inhomogeneous medium
著者: Guanghui Hu, Andreas Rathsfeld, Jiayi Zhang, Ruming Zhang
分類: math.AP（数学解析 - 偏微分方程式の解析）
発表日: 2025年10月17日（プレプリント）
論文リンク: https://arxiv.org/abs/2510.14070

要約

本論文は、鉛直方向に周期的で水平方向に不変の無限不均質二次元媒質に対して、新しい放射条件を提案している。古典的なRayleigh展開放射条件は、媒質が半平面内で不均質であることを要求するため、本論文の状況には適用できない。著者らはFloquet理論を用いて上行波/下行波モードを導出し、これらのモードの展開により放射条件を定義している。下行放射条件は下行Dirichlet-to-Neumann写像をもたらし、無限不均質領域を鉛直方向に切断するために使用できる。高次波モードの漸近挙動に基づいて、上行/下行Dirichlet-to-Neumann写像の写像性質を証明し、新しい散乱問題に対応する双線形形式の強楕円性を検証し、可算個の下に有界な波数集合を除くすべての波数に対して一意可解性を証明している。

研究背景と動機

問題記述

本論文が解決しようとする中核的な問題は、鉛直方向に周期的で水平方向に不変の無限不均質媒質において、時間調和平面波散乱問題に対して適切な放射条件をいかに確立し、対応する境界値問題の適切性を証明するかである。

重要性

実用的価値: 不均質媒質はフォトニック結晶、周期波導などの実用的応用において重要な役割を果たしている
理論的意義: 古典的なRayleigh展開放射条件の適用範囲を拡張している
数値計算: 効率的な数値スキーム設計と収束性解析の理論的基礎を提供している

既存方法の限界

古典的Rayleigh展開: 均質媒質内の周期構造散乱にのみ適用可能
既存の不均質媒質研究: 主に均質背景媒質または限定吸収原理（LAP）の使用に限定されている
適用範囲の制限: 無限周期不均質媒質に対する明示的な放射条件が欠けている

研究動機

Floquet理論とHill方程式の解析を通じて、無限周期不均質媒質に対する新しい放射条件を確立し、LAP パラメータの使用を回避し、より強い可解性結果を提供する。

中核的貢献

新しい放射条件の提案: Floquet理論に基づいて、鉛直方向に周期的な無限不均質媒質に対する上行および下行放射条件を確立した
DtN写像の構築: 上行および下行Dirichlet-to-Neumann写像を確立し、非有界領域の有界化処理を実現した
写像性質の証明: 高次波モードの漸近挙動に基づいて、DtN写像のSobolev空間における連続性を証明した
適切性理論の確立: 変分形式の強楕円性を検証し、可算個の例外波数を除く一意可解性を証明した
明示的波モードの提供: Hill方程式を解くことにより明示的波モードを得、LAP方法の使用を回避した

方法の詳細

問題設定

二次元領域における時間調和Helmholtz方程式を考える： $\Delta u(x_1, x_2) + k^2\tilde{q}(x_1, x_2)u(x_1, x_2) = 0$

ここで屈折率関数は： $\tilde{q}(x_1, x_2) = \begin{cases} 1, & \text{in } \Omega_d^+ \\ q_0(x_1, x_2), & \text{in } \Omega \\ q(x_2), & \text{in } \Omega_b^- \end{cases}$

目標は下半空間 $\Omega_b^-$ における $2\pi$ -周期関数 $q(x_2)$ に対して適切な放射条件を確立することである。

中核的理論枠組み

1. Floquet理論の応用

Hill方程式に対して： $u_n''(x_2) + (k^2q(x_2) - \alpha_n^2)u_n(x_2) = 0$

ここで $\alpha_n = \hat{\alpha} + \frac{2\pi n}{p}$ であり、Floquet理論を用いて特性乗数と特性指数を解析する。

2. 波モードの分類

特性パラメータ $\eta_n = w_{n,1}(2\pi) + w'_{n,2}(2\pi)$ の値に基づいて、波モードを5つのカテゴリに分類する：

ケース(a): $\eta_n > 2$ 、指数増減モード
ケース(b): $\eta_n < -2$ 、位相付き指数モード
ケース(c): $|\eta_n| < 2$ 、振動モード
ケース(d): $\eta_n = 2$ 、周期的または線形増長モード
ケース(e): $\eta_n = -2$ 、反周期モード

3. 放射条件の定義

上行放射条件(URC)： $u(x_1, x_2) = \sum_{n\in\mathbb{Z}} C_n^+ \cdot u_n^+(x_2) e^{i\alpha_n x_1} \text{ in } x_2 > a^+$

下行放射条件(DRC)： $u(x_1, x_2) = \sum_{n\in\mathbb{Z}} C_n^- \cdot u_n^-(x_2) e^{i\alpha_n x_1} \text{ in } x_2 < a^-$

DtN写像の構築

上行DtN写像

$T^+: f \mapsto \sum_{n\in\mathbb{Z}} \frac{u_n^{+'}(d)}{u_n^+(d)} f_n e^{i\alpha_n x_1}$

下行DtN写像

$T^-: f \mapsto \sum_{n\in\mathbb{Z}} \frac{-u_n^{-'}(b)}{u_n^-(b)} f_n e^{i\alpha_n x_1}$

技術的革新点

LAP方法の回避: Floquet理論を通じて明示的波モードを直接構築し、限定吸収原理を不要にした
漸近解析: $|n| \to \infty$ 時の波モードの漸近挙動を詳細に解析し、 $|\mu_{n,1}| = O(|n|)$ を証明した
強楕円性の証明: 変分形式を強制項とコンパクト項の和に分解することにより、強楕円性を確立した
小波数解析: $k \to 0$ 時の挙動を特別に解析し、小波数の場合の一意可解性を証明した

実験設定

理論検証枠組み

本論文は主に理論研究であり、以下の方法により方法の有効性を検証している：

特殊ケースの検証: $q(x_2) \equiv 1$ の場合、新しい放射条件が古典的Rayleigh展開に退化することを証明した
漸近挙動解析: 各種パラメータの漸近性質を詳細に解析した
写像性質の証明: Sobolev空間の枠組みにおいてDtN写像の連続性を証明した

主要な仮定条件

仮定A: 特定の退化ケースを排除し、DtN写像の良定義性を確保する
周期性: $q(x_2)$ は $2\pi$ 周期性を有する
実値性: $q(x_2) > 0$ は実値関数である

主要な理論結果

中核的定理

定理4.1（DtN写像の連続性）

DtN写像 $T^+: H_{\hat{\alpha}}^{1/2}(\mathbb{R}) \to H_{\hat{\alpha}}^{-1/2}(\mathbb{R})$ は連続であり、正定数 $C$ が存在して以下が成立する： $\|T^+f\|_{H_{\hat{\alpha}}^{-1/2}(\mathbb{R})} \leq C\|f\|_{H_{\hat{\alpha}}^{1/2}(\mathbb{R})}$

定理4.2（強楕円性）

変分形式 $a(\cdot,\cdot): H_{\hat{\alpha}}^1(C) \times H_{\hat{\alpha}}^1(C) \to \mathbb{C}$ は $H_{\hat{\alpha}}^1(C)$ 上で強楕円的である。

定理4.3（小波数可解性）

$q(x_2) > \sin^2\theta$ がすべての $b-2\pi < x_2 \leq b$ に対して成立すると仮定すれば、変分問題はすべての $k \in (0,k_0)$ と全ての入射方向に対して一意解を有する。

定理4.4（一般可解性）

固定入射角 $\theta \in (-\pi/2, \pi/2)$ に対して、変分問題は離散部分集合に属さないすべての波数 $k > 0$ に対して一意解を有し、この離散集合の唯一の集積点は無限遠である。

漸近解析結果

大 $|n|$ 挙動: $\eta_n = O(e^{2\pi|n|})$ 、 $|\mu_{n,1}| = O(|n|)$
小波数挙動: $k \to 0$ のとき、 $\mu_{n,1} > |n| - \frac{1}{2\pi}\ln 2 > 0$ （ $n \neq 0$ ）
収束性: 放射条件級数の点ごとの収束性を証明した