昨日は休日出勤で、明日は学会出張のため、本日は休養日にあてる。 午前中は買い物と読書して、午後から並列化の勉強をしよう。

卒研のテーマとして、今年からOpenFOAMを使ってみることにした。 進学予定の学生なので、テーマそのもののアカデミックさを優先するよりも、 大学へ行っても使えるスキルを身につけて、「自分で調べてトライアンドエラー」を繰り返してもらうつもり。 これ以…

いがいと簡単そうで、結構大変そうだというのが判って来た。 これまで、あまり原理やモデルの特徴とか把握してこなかったので、良い勉強になると思ってがんばろう。

いまいち安定して計算できなかったMPSのプログラムだが、 粒子速度に合わせてCFL条件をきちんと守る可変タイムステップ導入（CFL Divergence-free形式を用いる圧力ソース項にする振動対策（＋粒子衝突対策など） 自由表面粒子の選択アルゴリズム（オリジナル…

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

「静止容器内に8割ほど水をいれておいて放置する」という数値実験をやってみた。 当然、理屈上は静止水面と、深さ方向に滑らかな圧力分布ができるはずなんだが、 どうも変な結果になってしまう。まるで沸騰しているみたいだ。むぅ。