Python : 常微分方程式

odeint, BDFの利用について Stiff な方程式なのか、うまく解けなかったときのメモです。 連続槽型反応器の濃度変化をグラフ化したく、常微分方程式を数値的に解くことを考えました。 問題反応器サイズ　1L のCSTRに、反応速度式 $${-r_A = k C_A^{1.558}}$$となる原料を 100mmol/Lの濃度で、流量 10, 3, 1.2, 0.5 L/hrの各流量変化させて供給した。各流量切り替え前には、タンク内の反応成分は押し出して、ゼロの状態から

Feynman Lectures on Physicsを読破したい。 https://www.feynmanlectures.caltech.edu/

OpenFOAM : Cavity flow チュートリアルをしっかりやってみます。その１

以下のチュートリアルケースをしっかりとトレースしてみました。 このチュートリアルでは、以下の形状で、上面以外は、壁に囲まれた空間です。上部壁は、x 方向に 1m/sで動いている。以下の図は、サイトから引用しています 初期は層流状態からスタートし、均一なメッシュで、icoFoam（非定常層流解析）ソルバーで、層流、等温、非圧縮流れとして解析します。 準備設定ファイル、計算結果を保存するフォルダを作成、その中にチュートリアルケースをコピーします。以下のコマンドを実行します。