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

技术摘要：针对向列型 Helmholtz–Korteweg 方程的 $C^0$ 混合内罚函数法

问题陈述
本研究旨在解决向列型 Helmholtz–Korteweg 方程的数值解问题，该方程是一个四阶偏微分方程，用于模拟向列型 Korteweg 流体（如向列型液晶）中时谐声波的传播。该方程在有界 Lipschitz 区域 $\Omega \subset \mathbb{R}^d$ ($d=2,3$) 上进行建模，并采用混合声软（sound-soft）和阻抗边界条件。

分析面临三个主要挑战：

1. 高阶性与各向异性： 与标准的 Helmholtz 或双调和问题相比，存在四阶各向异性项使得分析变得更加复杂。
2. 正则性与稳定性： 先前的分析依赖于围绕各向同性情形（$\beta=0$）的摄动论证，这要求各向异性参数 $\beta$ 相对于各向同性问题的稳定性常数足够小。然而，该常数在共振频率附近会发生退化，从而人为地迫使 $\beta \to 0$，即使各向异性问题本身可能仍然是适定的。
3. 离散化复杂度： 用于解决四阶问题的标准共形有限元方法（$C^1$-conforming）需要高阶多项式次数（例如，在 2D 中需要 $p \geq 5$，在 3D 中需要 $p \geq 9$），且在处理弯曲区域或三维问题时实现难度极大。

方法论
本文提出了一种 $C^0$-混合内罚 ($C^0$-HIP) 离散化方法。该方法利用 $H^1$-共形有限元空间作为体变量（volume unknowns），并利用间断伽辽金（DG）技术来处理缺乏 $H^2$-共形性的问题。

• 连续分析： 作者通过采用 Cordes 条件 来量化各向异性效应，从而完善了连续分析。他们证明了在改进后的各向异性参数 $\beta$ 的小性条件下，该问题是 Fredholm 指数为零的。至关重要的是，该条件独立于各向同性情形的稳定性常数，从而允许更广泛的物理相关参数范围。
• 离散化： 该方法引入了代表体变量法向导数的附加面变量（facet unknowns）。这些面变量使体变量解耦，从而能够通过 静态凝聚（static condensation） 来降低计算成本。离散形式使用提升算子（lifting operator）来定义离散拉普拉斯算子，并包含一个稳定项以控制跨单元界面的法向导数跳跃。
• 理论框架： 分析依赖于一种针对算子的非共形离散化的抽象框架，该框架处理的是紧扰动强制算子。作者证明了一个离散 Gårding 不等式，它在离散层面重现了连续问题中“紧性修正下的强制性”结构。这种结构结合渐近一致性和离散紧性论证，确保了稳定性和收敛性。

核心贡献

1. 改进的适定性分析： 本文证明了向列型 Helmholtz–Koteweg 方程在比先前建立的更广泛的各向异性参数 $\beta$ 范围内是适定的。通过利用 Cordes 条件，推导出的 $\beta$ 小性条件不会在各向同性问题的共振频率处发生退化。
2. $C^0$-混合内罚法： 作者引入了一种 $C^0$-HIP 离散化方法，避免了实现 $C^1$-共形单元（如 Argyris 或 HCT 单元）的困难，同时所需的自由度比完全间断伽辽金方法更少。
3. 稳定性与收敛性证明： 只要各向异性足够小，该方法被证明对于任何多项式次数 $p \geq 2$ 都是稳定的，且与空间维度无关。作者在最小正则性假设（$u \in H^2(\Omega)$）下证明了其收敛至连续解，并在额外的正则性条件下推导出了最优阶收敛率。
4. 三维模拟： 本工作展示了首次使用该特定方程对向列场对声波传播之各向异性效应进行的三维数值模拟。

结果

• 理论结果： 离散双线性形式满足离散 Gårding 不等式，确保了渐近稳定性。当网格尺寸 $h \to 0$ 时，离散解序列收敛至连续解。
• 数值实验：
  • 收敛性： 在单位圆盘上的数值测试证实了对于多项式次数 $p \in \{2, 3, 4, 5\}$，其预测的收敛率为 $O(h^{p-1})$。
  • 稳定性： 实验验证了在 2D 中 $\beta > -\alpha$ 且在 3D 中 $\beta \in (-\alpha, 3\alpha)$ 时具有稳定性，这与基于 Cordes 条件推导出的理论界限一致。
  • 各向异性： 2D 和 3D 模拟展示了波传播的各向异性特征。结果表明，波速和衰减随向列型定向 $\mathbf{n}$ 的变化而变化，导致与各向同性情形截然不同的散射模式。
  • 效率： 通过对比自由度和非零项数量可以发现，由于通过静态凝聚消除了内部泡状函数（interior bubbles），与标准的 $C^0$-内罚法相比，混合策略显著降低了高阶多项式下的计算成本。

意义
本文声称其主要意义在于为求解向列型 Helmholtz–Korteweg 方程提供了一个稳健且灵活的数值框架，这对于模拟向列声学（nematoacoustics）至关重要。通过将稳定性分析与各向同性稳定性常数解耦，这项工作消除了对各向异性参数的人为限制，使得模拟更强的向列效应成为可能。此外，所提出的 $C^0$-HIP 方法克服了 $C^1$-共形方法的实际应用障碍，使得在三维及弯曲区域进行高阶四阶各向异性问题模拟变得可行。数值示例作为一种概念验证，展示了该方法捕捉复杂各向异性波现象的能力。