Cahn-Hilliard方程式(2次元)をExcel VBAで計算してみる。

保存系のTDGL方程式であるCahn-Hilliard方程式を計算してみよう。式は

∂φ/∂=　∇²（φ³ - φ　-　∇²φ）

だ。128x128の格子で計算。これもあっという間に計算できる。

保存系というのは

∂φ/∂ｔ　+　∇・J　=0

J=　-δH/δφ

H=φ＾2/2　-　φ＾4/4　+（∇φ）＾2/2

から来てる。