与线性弹性的瞬态 PDE 相比,稳态方程的条件较差。如果没有短暂的术语,它们是否天生就病态?
运算符矩阵的条件数似乎相当大(至少) 对于我已经解决的所有问题。这对于线性弹性来说是相当正常的吗?对于其他一些计算科学问题,我通常不会看到这么大的条件数。
我没有尝试对线性弹性进行 FEM 离散化;FEM 会产生比 FVM/FD 更好的调节系统吗?
与线性弹性的瞬态 PDE 相比,稳态方程的条件较差。如果没有短暂的术语,它们是否天生就病态?
运算符矩阵的条件数似乎相当大(至少) 对于我已经解决的所有问题。这对于线性弹性来说是相当正常的吗?对于其他一些计算科学问题,我通常不会看到这么大的条件数。
我没有尝试对线性弹性进行 FEM 离散化;FEM 会产生比 FVM/FD 更好的调节系统吗?
应用于二阶微分算子的任何方法(有限元、有限体积、有限差分)的刚度矩阵的条件数总是增长为其中指数等于微分算子的阶数。因此,如果你只是让网格足够细,你可以让条件数任意大——事实上,是一个非常小的值,可能是由于您的单元格或网格点数量相对较少。
您没有观察到此问题的原因是因为对于时间相关(动态)问题,需要反转的矩阵是以下任一形式(一次导数的抛物线问题)或(具有两个时间导数的双曲问题)。质量矩阵有条件数因此您需要反转的总矩阵的条件数为要么在上述两种情况下。换句话说,坏矩阵乘以一个小数使问题的条件更好。