拉格朗日乘数空间在数学上过于丰富

计算科学 有限元 数值分析 稳定
2021-12-24 10:17:24

背景:

拉格朗日乘子法已被应用于许多领域,例如接触问题、材料界面、相变、刚性约束或沿界面滑动。
众所周知,拉格朗日乘数空间的错误选择或设计将在拉格朗日乘数上产生振荡结果(不稳定问题)。大量文献已经说明了这一观察结果,并且已经进行了一些修改或改进以消除通常由 inf-sup 条件偏差引起的振荡。

问题:

在阅读有关 XFEM 的文献时,我遇到了以下以红色突出显示的论点,这是非常数学的。如何解释或理解空间局部过于丰富而导致违反 inf-sup 条件?感谢您的任何贡献。

拉格朗日乘数

2个回答

您必须求解的鞍点矩阵采用以下形式:

[ABTB],
其中是无约束矩阵,是观察浸入界面沿线节点处解值的矩阵。AB

鞍点矩阵的条件部分由的最小奇异值控制。的最小奇异值越小,条件越差。有关详细信息,我强烈推荐以下论文:BB

克伦德尔、沃尔夫冈、瓦莱里娅·西蒙奇尼和沃尔特·祖勒纳。“鞍点问题的稳定性估计和结构光谱特性。” 数值数学 124.1 (2013): 183-213。https://arxiv.org/pdf/1202.3330.pdf

现在让我们详细看一下浸没在网格中的界面的矩阵,如下图所示(黑色为网格,红色为浸没界面):B

边界锁定

矩阵采用以下形式:B

B=[1/21/21/21/21/21/2]

但是,如果将乘以一个高度振荡的向量,则结果很小,因为正负贡献往往会抵消。例如:B

[1111]large norm B=[1/2000]small norm

因此,具有较小的奇异值,对应于振荡界面向量。因此,鞍点矩阵是病态的,并且网格越细,情况就越差。B

每个拉格朗日乘数对应一个约束。因此,如果拉格朗日乘数的空间太大,那么在不显着限制可用于满足问题物理性质的未知数的数量的情况下,您将无法再同时满足太多约束。这就是锁定中发生的情况:每个约束都会将未知数的数量减少一个,并且最终会导致未知数太少而无法在物理上准确。