暗黙のDAE質量行列Python

Question

Pythonで解くことになるM * y '= f（y）問題があります。ここで、Mは質量行列、y'は微分、yはベクトルです.y1 、y2等はrの異なる点を指す。 Pythonで同様の問題に質量行列を使用した人はいますか？

問題は、rおよびz方向の2D問題です。 r方向は離散化され、問題を1次元問題に還元する。質量行列は、対角線上に1と0の対角行列です。

Answer 1

質量行列は特異であるため、これは「微分代数方程式」です。 SUNDIALSライブラリのIDAソルバなど、DAEの既製のソルバを見つけることができます。 SUNDIALSには、scikit.odesパッケージにpythonバインディングがあります。