粒子法

粒子法(りゅうしほう)とは、連続体に関する方程式を数値的に解くための離散化手法の一つで、計算対象物を粒子の集まりとして表すことからこのように呼ばれる。

主に流体解析や構造解析に用いられる手法で、代表的なものとして, SPH(Smoothed Particle Hydrodynamics)法や、MPS(Moving Particle Semi-implicit)法などがある。また、連続体の離散化手法とは異なるものであるがDEM(Distinct Element Method)法も粒子法と呼ばれる場合がある。

特徴と利点

流体解析においては、ラグランジュ法に属し、粒子の位置変化が対流に相当するため対流項の計算が無いという特徴がある。(有限体積法有限要素法に代表されるオイラー法では対流項の扱いが最も煩雑かつ難解な部分である。)

その他の主なメリットとして、質量が保存すること、水面等の扱いが容易であること、計算格子の作成が不要な点が挙げられる。粒子集団で構成される計算空間は接続情報が無いため、格子に比べて作成や扱いが容易である。また、格子解法のように微小な水滴が消えてしまうことはない。よって、さざ波や水しぶきなどの流体表面の細部の挙動を解析しやすい。構造解析においては、変形する母材の扱いが容易。

CG領域ではSPH法をアレンジした粒子法や粒子/格子ハイブリッド手法が使われることがある。 → FLIP、Hybrido等

現状の課題と実用状況

有限要素法や有限体積法に比べると歴史が浅く、専用の解析ソフトウェアの数は極めて少ない。

現実世界に極めて近い計算モデル[?]で解析を行うため万能にも思われる粒子法であるが、従来の手法と比較した場合の明らかな短所が幾つか存在する。

*空間分解能を調整する技術が一般的でない。格子手法のように要所に計算点を集中させることができない。

*乱流モデルはじめとする諸々の物理モデルの整備が不十分。工学的に適用できる対象は限定的。

*流入/流出や大気開放といった境界条件の扱いが厄介。

*粒子系の特性に依存する非物理挙動が顕著。

粒子法は格子法に比べると研究分野として未成熟である。粒子群の安定性といった根本的な課題が未解決のままである(張力不安定性等)。

著名なソフトウェア等

産業界に向けては、例えば東大発ベンチャーのプロメテック・ソフトウェア株式会社がParticleworks/Granuleworks(MPS/DEMに基づく流体/粉体解析ソフトウェア)を開発・販売しており、粒子法のエンジニアリング実務への応用も進行している。

  • SPlisHSplasH(SPH)
  • DualSPHysics (SPH)
  • GPUSPH (SPH)
  • CFDEM/LIGGGHTS (DEM/SPH)
  • プロメテック/Particleworks (MPS)
  • 株式会社富士テクニカルリサーチ/MPS-RYUJIN (MPS)

パッケージの一部に粒子法を含む製品

  • Ansys/LS-DYNA (SPH)
  • RealFlow(Cinema4D)(SPH/Hybrido)

主な手法

参考文献

和書

  • 越塚誠一、日本計算工学会(編):「粒子法」、丸善株式会社、 ISBN 978-4-621-07522-7 (2005年2月)。※ MPS法
  • 越塚誠一(編,著): 「粒子法シミュレーション-物理ベースCG入門」、培風館、ISBN 978-4-563-01574-9(2008年2月28日)。
  • W.G.Hoover, 志田 晃一郎 (訳):「粒子法による力学 - 連続体シミュレーションへの展開」、森北出版、ISBN 978-4627919310(2008年12月9日)。※ SPH法
  • 越塚誠一、柴田和也、室谷浩平:「粒子法入門 - 流体シミュレーションの基礎から並列計算と可視化まで」、丸善、ISBN 978-4621088340(2014年6月28日)。※ MPS法
  • 後藤仁志:「粒子法」、森北出版、ISBN 978-4627922310(2018年1月17日)。※ 代表的な手法(SPH法・MPS法)を統一的に解説。

その他

  • 越塚誠一:「粒子法シミュレーションの大規模化と高速化」、スーパーコンピューティングニュース(東京大学情報基盤センタ発行)、Vol.11(No.1), pp.125-138 (2009)。 (url=http://www.cc.u-tokyo.ac.jp/support/press/news/VOL11/special/200902SP-koshizuka.pdf )
有限差分法
放物型偏微分方程式
  • FTCSスキーム(英語版)
  • クランク・ニコルソン法(英語版)
双曲型偏微分方程式
  • ラックス・フリードリヒ法(英語版)
  • ラックス・ウェンドロフ法(英語版)
  • マコマック法(英語版)
  • 風上スキーム(英語版)
  • 特性曲線法
その他
有限体積法
  • ゴドノフスキーム(英語版)
  • 高分解能スキーム(英語版)
  • 保存法則用単調性上流中心差分スキーム(英語版)
  • 移流上流分離法(英語版)
  • リーマン解法(英語版)
有限要素法
  • hp-FEM(英語版)
  • 拡張型有限要素法(英語版) (XFEM)
  • 不連続ガラーキン法(英語版) (DG)
  • スペクトル要素法(英語版) (SEM)
  • モルタル有限要素法(英語版)
メッシュフリー法・粒子法
  • SPH法
  • MPS法(英語版)
  • MPM法(英語版)
領域分割法
  • シューア補元法(英語版)
  • 仮想領域法(英語版)
  • シュヴァルツ交代法加法シュヴァルツ法(英語版)抽象加法シュヴァルツ法(英語版)
  • ノイマン・ディレクレ法(英語版)
  • ノイマン・ノイマン法(英語版)
  • ポアンカレ・ステクロフ法(英語版)
  • バランシング領域分割法(英語版) (BDD)
  • BDDC法(英語版)
  • FETI法(英語版)
  • FETI-DP法(英語版)
その他