Analytical and Gegenbauer spectral solutions of damped wave equations with singular potentials in ℝn

技術要約：$\mathbb{R}^n$ における減衰波方程式の解析的およびゲンツェンバッハ・スペクトル解法

問題定式化
本研究は、長方形領域 $\Omega = [0, \pi]^n$ で定義される多次元減衰波方程式の解の存在、一意性、および数値近似を調査するものである。支配方程式は、形式 $|x|^{-\alpha}$ を持つ特異なソース項を含む、混合双曲型・放物型演算子である：
$$\frac{\partial^2 u}{\partial t^2} + \gamma \frac{\partial u}{\partial t} - \Delta u + \frac{u}{|x|^\alpha} = f(x, t)$$
本問題は、コーシー初期条件 ($u(x,0)=\phi(x)$, $u_t(x,0)=\psi(x)$) および非斉次ノイマン境界条件 ($\partial u / \partial \nu = g_j$) の下で定式化されている。解析は、特異性パラメータが $0 < \alpha < n$ かつ減衰係数が $\gamma > 0$ である領域に焦点を当てている。特異ポテンシャルの存在は、可積分性と正則性に関して重大な解析的課題をもたらすが、同時に、この混合的な性質を持つ演算子は、波の伝播と拡散的な散逸の両方を同時に扱うことを要求する。

手法
本論文は、解析的手法と数値的手法の二角的なアプローチを採用している：

1. 解析的枠組み（固有関数展開）：

   • 解は、過渡的な斉次部分 ($u_h$) と定常状態の特異部分 ($u_p$) に分解される。
   • 特異解は、時間に関するフーリエ級数を用いて構成され、これにより時間依存の偏微分方程式は、各時間調和 $k$ に対する結合楕円方程式系へと還元される。
   • 複素数値空間関数 $U_k(x) = A_k(x) - iB_k(x)$ を導入することで、系は複素ヘルムホルツ型方程式 $(L - \lambda_k)U_k = F_k$ （ここで $L = -\Delta + |x|^{-\alpha}, \lambda_k = k^2 - i\gamma k$）へと簡約される。
   • 解は、ノイマン・ラプラシアンの固有関数 $\Phi_m(x) = \prod \cos(m_j x_j)$ を用いた二重級数として明示的に表現される。係数は、内部ソースと境界フラックスの両方を考慮した上で、直交基底への投影によって決定される。

2. 数値的枠組み（ゲンツェンバッハ・スペクトル・コロケーション法）：

   • パラメータ $\lambda=1$（これは最適なノード分布のためにチェビシェフ多項式へと還元される）を用いたゲンツェンバッハ多項式による高次スペクトル法を開発した。
   • 領域の幾何学的形状を扱うために、テンソル積格子上のチェビシェフ・ガウス・ロバット（CGL）ノードを用いた空間離散化を行う。
   • 特異ポテンシャル $|x|^{-\alpha}$ は、漸近的挙動を維持しつつ数値的なオーバーフローを防ぐため、原点付近において $V_\epsilon(x) = (|x|^2 + \epsilon^2)^{-\alpha/2}$ を用いて正則化される。
   • 時間積分には、硬い（stiff）拡散項を暗黙的に扱い、双曲型項を陽的に扱うことで、制限的なCFL条件なしに安定性を確保する、二次のインプリシット・エクスプリシット（IMES）スキームを用いる。

主要な貢献および理論的結果
本論文は、解の厳密な数学的性質を確立している：

• 係数の収束： 展開係数が漸近的に $O(|m|^{-n-\alpha})$ で減衰することを証明しており、$0 < \alpha < n$ である限り、級数がウェルディファインド（定義可能）であることを保証している。
• 安定性と非共鳴： 減衰係数 $\gamma > 0$ が極めて重要であることを示している。これは、複素スペクトル分母 $|\mu_m - \lambda_k| \geq \gamma k > 0$ を保証し、ゼロ除算を防ぎ、純粋な共鳴現象を排除する。
• 正則性とソボレフ空間：
  • 解 $u$ は $L^2(\Omega)$ において有界かつ適切に定まる（well-posed）。
  • 空間勾配 $\nabla u$ が平方可積分（すなわち $u \in H^1(\Omega)$）であるのは、より厳格な条件 $0 < \alpha < n-2$ の下でのみである。$n-2 \leq \alpha < n$ の場合、勾配は原点において特異性を示す。
• 大域的安定性： 解の $L^2$ ノルムが $e^{-\gamma t/2}$ として指数関数的に減衰することを示す明示的なノルム推定を導出し、大域的安定性とエネルギー散逸を裏付けている。
• 誤差解析： ゲンツェンバッハ・コロケーション法の理論的な誤差境界を確立した。空間離散化誤差は $O(N^{-(n-\alpha-1/2)})$ の代数的レートで収束するが、これは特異性によって誘発される正則性の喪失を反映している。弱い特異性の場合、準スペクトル収束が観察される。

数値実験の結果
1次元、2次元、および3次元における包括的な数値実験を実施した。

• 検証： シミュレーションは理論的な収束率を検証している。弱い特異性（$\alpha = 0.3$）の場合、本手法は高次の収束（レート $> 3.5$）を達成する。臨界閾値に近い強い特異性（$\alpha \approx n - 0.5$）の場合、収束率は理論的な代数的限界へと緩やかに低下する。
• エネルギー挙動： 数値プロットは、減衰パラメータ $\gamma$ と一致してエネルギーが時間とともに指数関数的に減衰することを確認している。
• 次元のスケーリング： 結果は、本手法が次元 $n=1, 2, 3$ において頑健かつ安定であることを示しており、計算コストは $O(N^n)$ でスケーリングする。

意義
本論文は、標準的な存在定理では構成的な解を提供することが困難な、特異ポテンシャル下での $n$ 次元混合双曲型・放物型演算器に対する明示的な解析的表現を提供することで、既存の文献における空白を埋めるものである。本研究は、減衰係数が特異性の存在下における共鳴に対する基本的な安定化因子として機能することを実証している。さらに、統合された解析的・数値的アプローチは、局所的な特異性を持つ非均質媒体における波の伝播と拡散を特徴付けるための堅牢なモデルを提供する。著者らは、これらの知見が熱弾性学、流体構造相互作用、および散逸媒体における電磁波伝搬などの物理分野に直接適用可能であると位置づけている。