隐式DAE Mass Matrix Python

Question

我有一个问题M * y＆＃39; = f（y）将在Python中解决，其中M是质量矩阵，y＆＃39;导数和y是矢量，使得y1，y2等指r中的不同点。 有没有人在Python中使用类似问题的质量矩阵？

问题是r和z方向的2D问题。 r方向被离散化以将问题减少到1D问题。质量矩阵是对角线矩阵，对角线上有1和0。

Answer 1

由于质量矩阵是单数，因此这是一个差分代数方程＆＃34;。您可以找到DAE的现成求解器，例如SUNDIALS库中的IDA求解器。 SUNDIALS在scikit.odes包中有python绑定。