粒子法によるミルククラウンの解析

@tkoyama010

今回やったこと

PySPHを使用して粒子法でミルククラウンの解析をしました。
Paraviewを使用してミルククラウンをレンダリングしたアニメーションを作成しました。

タイトルが有限要素法ではなく粒子法になっていますが、これは間違えではありません。有限要素法でミルククラウン解析をすると届け出ていましたが間に合いませんでした。ごめんなさい。有限要素法で計算しないとだからといって何も出さないのももったいないですので今回は粒子法で問題を解いてきたという話をさせていただきます。使用したソフトはpysphというソフトでSPH法を使用したPythonのフレームワークの一つです。また、PyVistaで可視化をするつもりでしたが機能の調査が必要なため暫定的にParaviewで可視化を行いました。今回、SPHを可視化するにはPyVista使い方でまだ検討しなければいけない部分があることがわかりました。そこで、今回は可視化にはParaviewを使っています。

PySPHの解析事例

ミルククラウンの解析

解析の実行方法

$ sudo apt-get install python-pysph # pysphのインストール
$ git clone git@github.com:tkoyama010/pysph.git # 今回のスクリプトのチェックアウト
$ cd pysph
$ git checkout feat/milk-crown
$ cd ../
$ python3 ./pysph/pysph/examples/milk_crown.py # スクリプトの実行
$ pysph dump_vtk milk_crown_output # vtkへの出力