流体解析のあれこれ

流体解析をこよなく愛する有限体積法の信者によるブログ

常微分方程式

常微分方程式の数値解法(6)

SciPyのodeintを用いて以下のOrdinary Differential Equations (ODE)を解き,文献[1]の数値解を再現する. ここで,,,,であり,初期値は,,とした. [1] Field, R. J. and Noyes, R. M., Oscillations in chemical systems. IV. Limit cycle behavior in…

常微分方程式の数値解法(5)

SciPyのodeintを用いて以下のOrdinary Differential Equations (ODE)を解き,文献[1,2]のChemical Akzo Nobel problemの数値解を再現する. [1] Amat, S. et al., On a Variational Method for Stiff Differential Equations Arising from CHemistry Kinetic…

常微分方程式の数値解法(4)

SciPyのodeintを用いて以下のOrdinary Differential Equations (ODE)を解き,文献[1]のProblem HIRESの数値解を再現する. [1] Amat, S. et al., On a Variational Method for Stiff Differential Equations Arising from Chemistry Kinetics, Mathematics, …

常微分方程式の数値解法(3)

SciPyのodeintを用いて以下の厳密解のあるOrdinary Differential Equations (ODE)を解き,文献[1]の数値解を再現する. なお,文献[1]のの右辺第3項のはの間違いである.これら常微分方程式の厳密解は以下に示すとおりである. [1] Thohura, S. and Rahman, …

常微分方程式の数値解法(2)

SciPyのodeintを用いて以下のOrdinary Differential Equations (ODE)を解き,文献[1,2]のRobertson problemの数値解を再現する. ここで,初期値,,である. [1] Thohura, S. and Rahman, A., Numerical Approach for Solving Stiff Differential Equations…

常微分方程式の数値解法(1)

SciPyのodeintを用いて以下のOrdinary Differential Equations (ODE)を解き,文献[1]の数値解を再現する. ここで,,初期値,である. [1] Petzold, L., Automatic selection of methods for solving stiff and nonstiff systems of ordinary differential …