A hybrid $C^{0}$-interior penalty method for the nematic Helmholtz--Korteweg equation

技術要約：ネマティック・ヘルムホルツ・コートウェーク方程式に対する $C^0$ ハイブリッド・イントリア・ペナルティ法

問題設定
本研究は、ネマティック液晶などのネマティック・コートウェーク流体における定常音響波の伝搬をモデル化する、4階偏微分方程式であるネマティック・ヘルムホルツ・コートウェーク方程式の数値解法に取り組むものである。支配方程式は、古典的なヘルムホルツ方程式に2つの追加の4階項を加えたものであり、そのうちの1つはネマティック・ディレクター場 $\mathbf{n}$ に対して異方性を持つ。問題は、混合音響軟（sound-soft）およびインピーダンス境界条件を持つ、有界リプシッツ領域 $\Omega \subset \mathbb{R}^d$ ($d=2,3$) 上で定式化される。

解析においては、主に以下の3つの課題に直面する。

1. 高階性と異方性： 4階の異方性項の存在により、標準的なヘルムホルツ方程式やバイハーモニック問題と比較して、解析が複雑化する。
2. 正則性と安定性： 従来の解析は等方的なケース（$\beta=0$）の周りの摂動引数に依存しており、異方性パラメータ $\beta$ が等方的な問題の安定定数に対して小さくなければならないという制約があった。この定数は共振周波数付近で退化するため、異方的な問題がウェル・ポーズド（適切に解が存在）である場合であっても、人工的に $\beta \to 0$ を強いることになっていた。
3. 離散化の複雑さ： 4階問題のための標準的な適合有限要素法（$C^1$ 適合要素）は、高い多項式次数（例：2次元では Argyris 要素で $p \geq 5$、3次元では $p \geq 9$）を必要とし、曲面領域や3次元における実装が極めて困難である。

手法
本論文では、**$C^0$ ハイブリッド・イントリア・ペナルティ（$C^0$-HIP）**離散化を提案する。この手法は、体積未知数には $H^1$ 適合有限要素空間を用い、$H^2$ 適合性の欠如を扱うために不連続ガラーキン（DG）手法を利用する。

• 連続解析： 著者らは、異方性の影響を定量化するために**コードス条件（Cordes condition）**を採用することで、連続解析を洗練させた。彼らは、問題が改善された $\beta$ の小ささの条件の下で、インデックス0のフレドホルム性を満たすことを確立した。決定的なのは、この条件が等方的な問題の安定定数に依存しないことであり、これにより、より広範な物理的に関連のあるパラメータ範囲が可能となった。
• 離散化： 本手法は、体積未知数の法線微分を表す追加のファセット（facet）未知数を導入する。これらのファセット変数は体積未知数をデカップリングし、計算コストを削減するための**静的縮退（static condensation）**を可能にする。離散定式化では、離散ラプラシアンを定義するためにリフティング演算子を用い、要素界面における法線微分のジャンプを制御するための安定化項を含む。
• 理論的枠組み： 解析は、強コアクティブな作用素の非適合離散化に関する抽象的な枠組みに基づいている。著者らは、連続問題の「コンパクト性へのコアクティブ性」の構造を離散レベルでも再現する、離散**ガルーディング不等式（Gårding inequality）**を証明している。この構造を、漸近的一貫性と離散コンパクト性論理と組み合わせることで、安定性と収束性が保証される。

主な貢献

1. 改善されたウェル・ポーズド解析： 本論文は、ネマティック・ヘルムホルツ・コートウェーク方程式が、以前確立されていたものよりも大幅に広い範囲の異方性パラメータ $\beta$ に対してウェル・ポーズドであることを証明している。コードス条件を利用することで、導出された $\beta$ の小ささの条件は、等方的な問題の共振周波数において退化しない。
2. $C^0$ ハイブリッド・イントリア・ペナルティ法： 著者らは、$C^1$ 適合要素（Argyris や HCT 要素など）の実装の難しさを回避しつつ、完全な不連続ガラーキン法よりも少ない自由度を必要とする $C^0$-HIP 離散化を導入した。
3. 安定性と収束性の証明： 本手法は、異方性が十分に小さい場合、空間次元に依存せず、任意の多項式次数 $p \geq 2$ に対して安定であることが証明されている。論文では、最小限の正則性仮定（$u \in H^2(\Omega)$）の下での連続解への収束を確立し、追加の正則性の下で最適次数の収束率を導出している。
4. 3次元シミュレーション： 本研究は、この特定の式を用いた、ネマティック場による音響波伝搬への異方性の影響に関する初の3次元数値シミュレーションを提示している。

結果

• 理論的結果： 離散双線形形式は離散ガルーディング不等式を満たし、漸近的安定性を保証する。離散解の列は、メッシュサイズ $h \to 0$ の下で連続解に収束する。
• 数値実験：
  • 収束性： 単位円上の数値テストにより、多項式次数 $p \in \{2, 3, 4, 5\}$ に対して予測された収束率 $O(h^{p-1})$ が確認された。
  • 安定性： 実験により、2次元では $\beta > -\alpha$、3次元では $\beta \in (-\alpha, 3\alpha)$ において安定であることが検証され、これはコードス条件から導かれた理論的境界と一致している。
  • 異方性： 2次元および3次元のシミュレーションは、波の伝搬における異方的な性質を示している。結果は、波の速度と減衰がネマティック・ディレクター $\mathbf{n}$ の方位に依存して変化し、等方的なケースと比較して明確な散乱パターンを生じさせることを示している。
  • 効率性： 自由度と非ゼロ要素の比較により、内部バブルを静的縮退によって排除することにより、ハイブリダイゼーション戦略が高次多項式の場合において標準的な $C^0$-イントリア・ペナルティ法と比較して計算コストを大幅に削減することを実証している。

意義
本論文の主要な意義は、ネマティック・ヘルムホルツ・コートウェーク方程式を解くための堅牢かつ柔軟な数値フレームワークを提供することにある。これはネマトアコースティクス（nematoacoustics）のモデリングに不可みである。異方性の安定性を等方的な安定定数から切り離すことにより、本研究は異方性パラメータに対する人工的な制限を取り除き、より強いネマティック効果のシミュレーションを可能にした。さらに、提案された $C^0$-HIP 法は、$C^1$ 適合法の実際的な障壁を克服し、3次元および曲面領域における4階の異方性問題の高次シミュレーションを可能にする。数値例は、複雑な異方性波現象を捉えるための本手法の能力を示す概念実証として機能している。