Eigenを用いた2次元移流拡散方程式の数値計算 - Qiita

二次元空間での移流拡散方程式の数値計算です。拡散する物理量は濃度とか温度だと思ってください。 離散化手法として、時間離散化にクランク＝ニコルソン法、空間離散化に１次精度の風上差分法を用いています。 二次元での陰解法なので、TDMAで線順法（c.f.パタンカー）でもいいのですが、Eigenの疎行列ライブラリが使いたかったので一次方程式にして解きました。あんまりオブジェクト指向してませんがそこは適当です。 連立方程式のソルバーは適当に選んでます。このへんの知識をもう少し入れたいですね。 #include <iostream> #include <math.h> #include <fstream> #include <string> #include <Eigen/Sparse> #include <Eigen/Dense> using namespace std; using namespa

[lochtext]

数値計算しました