複素TDGL方程式をExcel VBAで計算してスパイラルパターンをGIFアニメにしてみる。

非線形偏微分方程式をExcel VBAで計算してGIFアニメにしてみるシリーズ。

まずは複素TDGL方程式（Complex Time Dependent Ginzburg Landau equation）

∂W/∂t = (1+iC0)W + (1+iC1)∇²W - (1+iC2)|W|²W

を計算する。C0=C1=0, C2=2として、128x128の格子、∇W=0 at 境界という条件で、初期条件はランダムとしてRe(W)を図示してGIFアニメにしてみた。

このパラメータだとらせん（スパイラル）パターンがでることが知られています。

初期条件その１

その２

なかなか気持ち悪い。。。

こういうのはBZ反応でも見られますよ。Excek VBAでこれも計算した結果がこちら。

https://sci.tea-nifty.com/blog/2011/11/bzexcel-2a42.html

これもGIFアニメにしてみようか（続く）。