Excelで蔵本・シバシンスキー方程式を解いてみる。

Kuramoto-Sivashinsky方程式：

∂ｆ/∂ｔ + ｆ∂f/∂x　= -λ∂^2xf/∂x^2 - ν∂^4xf/∂x^4

を三次風上差分+4次精度の拡散項で解いてみた。結果はこちら。

下から時間が増えていっている。規則正しい構造からだんだん乱れていっているのがよくわかる。