2016-06-28 9 views
0

Ich habe ein Problem M * y '= f (y), die in Python gelöst werden, wobei M die Massenmatrix, y' die Ableitung und y ein Vektor ist, so dass y1 , y2 usw. bezieht sich auf verschiedene Punkte in r. Hat jemand eine Massenmatrix für ein ähnliches Problem in Python verwendet?Implizite DAE-Massenmatrix Python

Das Problem ist ein 2D-Problem in r- und z-Richtung. Die r-Richtung wird diskretisiert, um das Problem auf ein 1D-Problem zu reduzieren. Die Massenmatrix ist eine Diagonalmatrix mit Einsen und Nullen auf der Diagonalen.

Antwort

0

Da die Massenmatrix singulär ist, ist dies eine "differential-algebraische Gleichung". Sie können Standard-Solver für DAEs finden, z. B. den IDA-Solver aus der SUNDIALS-Bibliothek. SUNDIALS hat Python-Bindungen im Paket scikit.odes.