どうにも安定しなかったMPS水槽実験だけど、圧力ポアソン方程式のソース項の取り扱いが原因だったようで（さんざん論文で指摘されていること）、 田中・益永(JSCES,2008)の通り、オリジナルのPND条件形式から、Divergence-free条件形式に変更したところかな…