流体解析ソフトウェア Particleworks
  • Home
  • Home
  • Case Examples
    • Particleworks解析事例
    • Particleworks Case Examples
  • Technical
    • 粒子法・MPS法
    • 技術コラム >
      • DX時代の製品開発プロセスとCAEの重要性 >
        • 第1回 序 略歴とコラム紹介
        • 第2回 DXとデジタルエンジニアリング
        • 第3回 製品開発プロセスの目指す姿
        • 第4回 DX時代のCAE
        • 第5回 評価CAEの概要と課題
        • 第6回 評価CAEの課題解決手法
        • 第7回 企画CAEの概要と課題
        • 第8回 企画CAEの運用と応用
        • 第9回 設計CAEの概要と課題
        • 第10回 設計CAEの課題解決の進め方
        • 第11回 開発プロセス運用の仕組み作り
        • 第12回 まとめと変革の時代に求められるエンジニア像
      • 粒子法のいま、そして未来へ >
        • 第1回 粒子法のいま
        • 第2回 SPH法におけるカーネル近似とカーネル関数の条件
        • 第3回 SPH法における空間離散化
      • 粒子法の非圧縮条件とは
      • 粒子法入門 >
        • 第1回 粒子法って何?
        • 第2回 粒子法は、他の方法とどう違うか
        • 第3回 粒子法の大きさと質量について
        • ​第4回 「粒子の動かし方」と「加速度の求め方」について
        • ​第5回 計算時間を短縮する方法について
    • Technical Column >
      • Growing the particle method, and its present state >
        • 1. Present State of the Particle Method
        • 2. Kernel Approximation and Kernel Function Conditions in the SPH Method (Preparation for Spatial Discretization)
      • Incompressibility of the particle method
      • Introduction to the particle method >
        • 1. What is a particle method?
        • 2. In what ways is the particle method different from other methods?
        • 3. Mass and volume of particles
        • 4. How to move particles and how to calculate accelerations of particles
        • 5. How to shorten the simulation time
    • 粒子法用語集
    • Particle Method Glossary
    • 参考文献・ウェブサイト
    • Reference Book/URL
    • 論文・講演
  • Contact
    • 導入の流れとライセンス形態
    • Particleworks / GranuleworksプリインストールGPU搭載ワークステーション
    • 開発元・パートナー
    • Developers, Partners
    • お問い合わせ
    • Contact Us

粒子法の非圧縮条件とは

粒子法の1つであるMPS(Moving Particle Semi-implicit)法では、非圧縮性流体を「密度一定」であるとして解いています。一方、有限体積法などの格子法では「速度の発散ゼロ」を使います。この違いがMPS法の特徴の1つなので、その違いについてここでは解説します。

 

粒子法の研究を開始した1993年に、既にSPH(Smoothed Particle Hydrodynamics)という粒子法がありましたが、陽的なアルゴリズムを用いて圧縮性流体を解く方法でした。工学的な流体力学の問題の多くは非圧縮性流体であり、格子法では半陰的なアルゴリズムを用いて解かれていました。そこで、半陰的アルゴリズムを用いて非圧縮性流体が解けるような粒子法を模索していました。そうしてたどり着いたのがMPS法です。ただ、これまでの格子法とは少し変えていまして、そこがMPS法の特徴であり、肝心かなめの中心部分になっています。

 

そもそもから話を始めると、物理の基本的な保存則は3個あり、質量保存則、運動量保存則、エネルギー保存則です。エネルギー保存則は質量保存則と運動量保存則から導くことができるので、実際には2個を解けば答えが出ます。流体力学では、質量保存則は「連続の式」、運動量保存則は「ナビエ・ストークス方程式」と呼ばれており、これら2個の式を解いています。粒子法でもそうです。非圧縮性が関係するのは質量保存則すなわち「連続の式」です。

 

圧縮性流体の連続の式は

$$ \frac{D \rho}{D t}+\rho \nabla \cdot \boldsymbol{u}=0 $$

です。非圧縮性流体では、式(1)の左辺の第1項がゼロであり、従って第2項もゼロになります。

$$ \frac{D \rho}{D t}=0 $$ $$ \nabla \cdot \boldsymbol{u}=0 $$

さらに、式(2)を時間で積分した式、すなわち密度\(\rho\)が時間的に一定であるという式

$$ \rho=\it{const} $$

も得られます(式(4)の右辺の\(\it{const}\)は一定値の意味)。

 

式(4)は「密度一定」ということなので、各時刻で密度を計算して、それらを一定値にしなさいということを意味します。これをMPS法では課しています(Koshizuka and Oka, 1996)。一方、格子法では式(3)すなわち「速度の発散ゼロ」を使います。粒子法でも「速度の発散ゼロ」を使うこともできます。いったいどの式が非圧縮条件として本質的なのでしょうか。あるいは、どれを課しても同じと考えていいのでしょうか(図1)。

Picture
図1 非圧縮条件はどれが本質的なのか?

まず、式(4)と式(3)は何が違うのかを考えてみましょう。そもそもに話を戻すと、連続の式は質量保存則です。これは、質量が一定であるという意味であり、非圧縮性流体では本質的には式(4)です。式(4)を時間で微分すると、密度の時間変化はゼロであるという式(2)が得られ、式(1)を介して「速度の発散ゼロ」という式(3)が得られます。

 

貯金にたとえると、毎月、貯金の総額が一定であることを残高から確認するのが式(4)で、出金と入金の総計がゼロであることを確認するのが式(3)です。計算を間違えなければ、どちらも同じことですが、もし、計算に間違う可能性がある場合には、残高を確認する前者の方が直接的すなわち本質的ではないでしょうか。

Picture
図2 斜面上に拘束される玉の運動

また、図2(a)のような斜面上を転がる玉を考えます。玉は斜面上に存在しなければならないというのが直接的な拘束条件であって、流体力学における式(4)に対応します。玉が斜面上に存在しなければならないという条件を時間で微分すると、玉の速度\(\boldsymbol{v}\)が斜面と平行でなければならないという条件が得られます(図2(b))。こちらは式(3)に対応します。この時、次の時刻の玉の位置への移動を、\(\boldsymbol{v} \times \Delta t\)、とすると、図2(c)のように玉は地面に少しめり込むことも考えられます。つまり、瞬間的には、拘束条件=拘束条件の時間微分、なのですが、数値計算では誤差が入り込むので、拘束条件≠拘束条件の時間微分、となってしまいます。そう考えると、図2(a)の方が図2(b)よりも満たすべき拘束条件として本質的です。粒子法でも同様で、式(3)で非圧縮条件を課すと、数値計算の誤差によって流体の密度がわずかずつ変化してしまいます。そして時間とともに変化は蓄積されます。従って、非圧縮条件は式(4)が本質的で式(3)は派生的です。

 

MPS法より前に開発された粒子法であるSPHにおいても「密度一定」を非圧縮条件として課す方法があります(Shao and Lo, 2003)。この方法はISPH(Incompressible SPH)と呼ばれ、MPS法の非圧縮条件をSPHに取り入れたもので、元の論文(Shao and Lo, 2003)にはそのことが書かれています。

 

では、粒子法において「密度一定」と「速度の発散ゼロ」のどちらが優れているでしょうか。筆者は、「密度一定」の方が粒子法の特徴であるロバスト性(頑強さ)の点から優れていると考えています。詳細に見れば、それぞれ利点と欠点があります。「密度一定」の利点は、(1) 密度の誤差が蓄積しない、(2) 粒子配置が均等になる、(3) 自由表面境界条件が簡単、というものです。その一方で、圧力が振動しやすいという欠点があります。「速度の発散ゼロ」では利点と欠点は反対になります。そこで、非圧縮条件として式(2), (3), (4)を混ぜるという方法が考案されました(近藤, 越塚, 2008; Tanaka and Masunaga, 2010)。こうした方法ですと、「密度一定」の利点を残しつつ、欠点である圧力振動が改善されます。

参考文献:
S. Koshizuka and Y. Oka, "Moving-Particle Semi-implicit Method for Fragmentation of Incompressible Fluid," Nucl. Sci. Eng. 123, 421-434 (1996).
S. Shao and E. Y. M. Lo, "Incompressible SPH Method for Simulating Newtonian and non-Newtonian Flows with a Free Surface," Adv. Water Resources 26, 787-800 (2003).
近藤雅裕, 越塚誠一, "MPS法における不自然な数値振動の抑制," 日本計算工学会論文集, Paper No.20080015 (2008).
M. Tanaka and T. Masunaga, "Stabilization and Smoothing of Pressure in MPS Method by Quasi-Compressibility," J. Comput. Phys. 229, 4279-4290 (2010).
  • 粒子法の非圧縮条件とは
技術コラムTOPに戻る

 著者ご紹介

画像
東京大学大学院 工学系研究科 システム創成学専攻 教授
プロメテック・ソフトウェア創業者・社外取締役
粒子法(MPS法)考案者
越塚 誠一 先生
​<著者プロフィール>
2004年 東京大学大学院 工学系研究科 システム量子工学専攻 教授
2005年 プロメテック・ソフトウェア株式会社 取締役
2008年 東京大学大学院 工学系研究科 システム創成学専攻 教授
(Sitemap)
​Home
Case Examples
 - Particleworks Case Examples
Technical
 - Technical Column
 - References
 - Papers, Lectures
 - Particle Method Glossary
Contact
 - Developer, Partners
 - Contact Us
Privacy Policy
Terms of Use
GDPR PRIVACY POLICY
(Related Sites)
Prometech Software Site
Granuleworks Site
Prometech Simulation
​Conference Site
GDEP Solutions Site
Contact Us
[Developer, Main Domestic / Global Dealer​]
  Prometech Software, Inc.
Prometech Software, Inc.
​URL: www.prometech.co.jp
​
E-mail: web@prometech.co.jp
(Sitemap)
Home
事例
 - 解析事例
Learning
 - 粒子法・MPS法
​ - 技術コラム
 - 粒子法用語集
 - 参考文献・ウェブサイト
 - 論文・講演
お問い合わせ・ご相談
 - 導入の流れとライセンス形態
 - Particleworks / Granuleworks
   プリインストールGPU搭載
   ワークステーション
 - 開発元・パートナー
 - お問い合わせ
プライバシーポリシー
利用規約
GDPR プライバシーポリシー
(Related Sites)
Prometech Software サイト
​Granuleworksサイト
​Prometech Simulation
Conference サイト
Particleworks Europe サイト
プロメテックCGリサーチ サイト
GDEP Solutions サイト
-  動作確認済み GPU搭載ワークステーション
HPC WORLD サイト
お問い合わせフォーム
[開発元・国内、海外総販売店]
  プロメテック・ソフトウェア株式会社
Prometech Software, Inc.
URL: www.prometech.co.jp
E-mail: web@prometech.co.jp

ⓒPrometech Software, Inc.
  • Home
  • Home
  • Case Examples
    • Particleworks解析事例
    • Particleworks Case Examples
  • Technical
    • 粒子法・MPS法
    • 技術コラム >
      • DX時代の製品開発プロセスとCAEの重要性 >
        • 第1回 序 略歴とコラム紹介
        • 第2回 DXとデジタルエンジニアリング
        • 第3回 製品開発プロセスの目指す姿
        • 第4回 DX時代のCAE
        • 第5回 評価CAEの概要と課題
        • 第6回 評価CAEの課題解決手法
        • 第7回 企画CAEの概要と課題
        • 第8回 企画CAEの運用と応用
        • 第9回 設計CAEの概要と課題
        • 第10回 設計CAEの課題解決の進め方
        • 第11回 開発プロセス運用の仕組み作り
        • 第12回 まとめと変革の時代に求められるエンジニア像
      • 粒子法のいま、そして未来へ >
        • 第1回 粒子法のいま
        • 第2回 SPH法におけるカーネル近似とカーネル関数の条件
        • 第3回 SPH法における空間離散化
      • 粒子法の非圧縮条件とは
      • 粒子法入門 >
        • 第1回 粒子法って何?
        • 第2回 粒子法は、他の方法とどう違うか
        • 第3回 粒子法の大きさと質量について
        • ​第4回 「粒子の動かし方」と「加速度の求め方」について
        • ​第5回 計算時間を短縮する方法について
    • Technical Column >
      • Growing the particle method, and its present state >
        • 1. Present State of the Particle Method
        • 2. Kernel Approximation and Kernel Function Conditions in the SPH Method (Preparation for Spatial Discretization)
      • Incompressibility of the particle method
      • Introduction to the particle method >
        • 1. What is a particle method?
        • 2. In what ways is the particle method different from other methods?
        • 3. Mass and volume of particles
        • 4. How to move particles and how to calculate accelerations of particles
        • 5. How to shorten the simulation time
    • 粒子法用語集
    • Particle Method Glossary
    • 参考文献・ウェブサイト
    • Reference Book/URL
    • 論文・講演
  • Contact
    • 導入の流れとライセンス形態
    • Particleworks / GranuleworksプリインストールGPU搭載ワークステーション
    • 開発元・パートナー
    • Developers, Partners
    • お問い合わせ
    • Contact Us