[论文解读] $C^0$ Discontinuous Galerkin Methods for a Kirchhoff Plate Contact Problem
本文针对Kirchhoff板接触问题中摩擦作用的四阶椭圆变分不等式，提出并分析了$C^0$间断伽辽金（DG）方法。作者将五种现有DG格式——C0 IP、NIPG、Wells-Dung、Baker-DG与LCDG——进行扩展，提供了严格的相容性与稳定性分析，并推导出在弱正则性假设下的最优阶误差估计，数值结果验证了能量范数与$H^1$范数下收敛率接近1的最优收敛性。
Discontinuous Galerkin (DG) methods are considered for solving a plate contact problem, which is a 4th-order elliptic variational inequality of second kind. Numerous $C^0$ DG schemes for the Kirchhoff plate bending problem are extended to the variational inequality. Properties of the DG methods, such as consistency and stability, are studied, and optimal order error estimates are derived. A numerical example is presented to show the performance of the DG methods; the numerical convergence orders confirm the theoretical prediction.
研究动机与目标
- 将$C^0$间断伽辽金方法扩展至求解源自Kirchhoff板摩擦接触问题的四阶椭圆变分不等式。
- 建立所提DG格式在该类变分不等式问题中的相容性与稳定性性质。
- 在弱正则性假设下，推导$C^0$ DG方法的最优阶误差估计。
- 通过$h$-型加密的基准问题，数值验证理论收敛速率。
提出的方法
- 通过引入数值迹与罚函数项，将五种现有$C^0$ DG格式——C0 IP、NIPG、Wells-Dung、Baker-DG与LCDG——适配于四阶变分不等式。
- 将离散问题表述为凸最小化问题：$\min_{u^h \in V^h} \frac{1}{2}a(u^h,u^h) + j_h(u^h) - (f,u^h)$，其中$j_h$表示接触势能。
- 采用原始-对偶定点算法求解离散优化问题，对偶步骤中应用近端算子处理$\ell_1$-范数项。
- 定义能量范数$\interleave \cdot \interleave$与$H^1(\Omega)$半范数，用于衡量误差收敛性。
- 采用一致且稳定的$C^0$ DG格式（LCDG），并给出可量化的稳定性条件，避免依赖过大的罚参数。
- 以最细网格（$h = 1/64$）的解作为参考解，计算较粗网格的误差，从而评估收敛阶。
实验结果
研究问题
- RQ1现有用于双调和方程的$C^0$ DG方法能否扩展至处理板接触问题中出现的四阶变分不等式？
- RQ2所扩展的$C^0$ DG格式在该变分不等式问题中具有怎样的相容性与稳定性性质？
- RQ3所提DG方法在能量范数与$H^1$范数下是否达到最优收敛速率？
- RQ4不同DG格式（如C0 IP、LCDG、Wells-Dung）在精度与收敛行为方面有何差异？
- RQ5原始-对偶定点算法能否高效求解含$\ell_1$-正则化的离散凸优化问题？
主要发现
- 所提$C^0$ DG方法在能量范数下实现了最优收敛阶，所有五种格式的数值收敛率均接近于1。
- 数值结果表明，$H^1(\Omega)$半范数下的收敛阶也接近于1，与理论预测一致。
- LCDG方法提供了稳定的格式，其罚参数可量化，避免了为保证稳定性而使用任意大的罚值。
- Wells-Dung与Baker-DG格式的性能与C0 IP及NIPG相当，且在不同罚参数$\eta$下表现出一致的收敛行为。
- 原始-对偶定点算法能高效求解离散问题，近端算子使$\ell_1$-正则化项可实现显式更新。
- 对所有方法，当$h$减半时，能量范数下的误差约减小为原来的二分之一，证实了一阶收敛性。
更好的研究，从现在开始
从阅读论文到最终审阅，大幅缩短您的研究时间。
无需绑定信用卡
本解读由 AI 生成，并经人工编辑审核。