趣味でスペクトル法_#003_2次元周期流 - Qiita

概要 文献12を参考に，2次元非圧縮性流体の数値解析プログラムをFortranで作成した．2次元周期流を仮定し，渦度輸送方程式をスペクトル法で解いた．離散フーリエ変換にはIntel MKL3を用いた． $1024^2$自由度の解析．Paraview4で可視化した． 記事内へのリンク 支配方程式の導出 初期条件の与え方 2次元DFT 非線形項の評価（変換法） 時間積分法 計算例 ソースコード 更新履歴 2019.12.01 : 公開 支配方程式の導出 2次元の非圧縮性流体を仮定し，スペクトル法で解く支配方程式を導出する．156 渦度輸送方程式の導出 基礎方程式は次に示す連続の式（質量保存則）とNavier-Stokes方程式（運動量保存則）である．56 連続の式 Navier-Stokes方程式 $$ \begin{align} \frac{\partial \mathbf{u}}{\par