流体解析ソフトウェア 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

​第2回 粒子法は、他の方法とどう違うか

前回のコラムでは、粒子法を使えば様々な難しい問題を容易にシミュレーションできることを紹介しました。粒子法はなぜそのようなことをできるのでしょう。今回のコラムでは、粒子法を他の計算手法と比較することでその理由を説明したいと思います。

1. 差分法

流体シミュレーションの伝統的な計算手法として、差分法があります。差分法を使ってパイプの中の水の流れ(図1)を計算する際の様子を図2に示します。計算領域を網目のような格子(計算格子)に分割して考え、格子のある位置(四辺形の中心や、1辺の中央など)に計算点(計算上の観測点のようなもの)を置きます。それらの計算点に、変数(圧力、速度、流体密度、温度などの値)の情報を持たせ、時間が経過するとともに変数の値を更新していきます。差分法では、格子を移動させたり変形したりすることは基本的には行いません注1ので、計算点も移動しません。つまり空間に固定された同じ位置で流れを観測することになります。温度の時間変化を求める場合で考えると、差分法ではある場所に温度計を固定して温度を測り続けることになります。場所を固定していますので、毎時刻においてその位置を通過する別の流体の温度を測定(計算)していることになります。
画像
図1 パイプ内の流体の流れ
画像
図2 差分法での計算格子(メッシュ)による流体の表現
(矢印は流速を表しています。本当は全ての格子の各辺に垂直な矢印を、また各格子の中央に圧力の計算点を描くべきですが、簡略化して描画しています。)

注1

例外として、Arbitrary Lagrangian Eulerian (ALE)法のように格子を変形させる方法もあります。

2. 粒子法

パイプの中の水の流れ(図1)を粒子法で計算する際の模式図を図3に示します。粒子(計算点)を使って流体を分割します。格子は使いません。ある小さな流体領域を1つの流体粒子で表し、流体の全体を粒子の集合で表すということです。差分法とは異なり、粒子法の計算点は移動します。計算点の移動速度は、その流体粒子の速度となります。図4のような複雑に表面形状が変わっても、この計算点の移動により動画1に示すように簡単に対応できるのです​。
​粒子法が簡単なもう1つの理由は、差分法で必要な移流項注2の計算をする必要が無いからです。移流項とは、流体のように位置が時々刻々移動する物体のある物理量(速度、温度など)の時間変化を、空間に固定された計測点での情報のみを用いて求める際に必要となる項のことです。この移流項の計算は複雑になることが多いのですが、粒子法の場合は(計算点を流れとともに移動させるおかげで)移流項を計算する必要がなく、簡単なのです。
具体的な例として、流体の温度の時間変化を求める場合で差分法と粒子法を比較してみましょう。差分法では場所を固定して温度を計測し続けるのに対して、粒子法ではある流体粒子に着目しその粒子の温度を計測し続けるという違いがあります。この観測方法の違いにより、粒子法では移流項を計算する必要がなく簡単なのです。ここでは温度の時間変化の場合で説明しましたが、流速の時間変化でも同様です。差分法では空間に固定した位置で流速の変化を考え、粒子法ではある流体粒子の流速の変化を見ます。
画像
図3 粒子法での粒子による流体表現
(矢印は流速を表しています。本当は全ての粒子に流速を表す矢印を、また粒子の中心に圧力を表す点を描く​​べきですが、簡略化して描画しています。)
画像
図4 粒子による複雑な表面形状の表現の模式図
動画1 複雑な表面形状をもつ流体の粒子法による計算例

注2

移流項は、ある位置の上を今まさに通過しようとする流体の物理量の時間変化を、空間に固定された位置での情報のみを基に求めようとする際に現れる項です。次の時刻に流体の位置が変わる(下流に流れる)ことに起因する物理量の変化を移流項は表します。もう少し具体的に言うと、流体の微小な塊が現在の場所から下流の場所に移動しようとするときに、同時刻において両者の位置で物理量にあらかじめ差がある場合の影響を表します。
移流項を数式でしっかり理解したい人は、連続体力学や流体力学の教科書[1,2]、粒子法入門[3]を読むことをお勧めします。多くの教科書では、実質微分やラグランジュ微分などの節で移流項を説明しています。ただし、粒子法では移流項は計算する必要がないので、しっかり理解できなくても心配する必要はありません。「粒子法では移流項という面倒な項を計算する必要がなくて簡単だ。」ということを理解いただければ、まずは十分だと思います。
参考文献:
1. 佐野理,"連続体の力学 (基礎物理学選書)", 裳華房, 2000年, ISBN-13: 978-4785321376
2. 日野幹雄, "流体力学",朝倉書店, 1992年, ISBN-13: 978-4254200669
3. 越塚誠一、柴田和也、室谷浩平、"粒子法入門"、丸善出版、2014年6月25日
  • 第1回 粒子法って何?
  • 第2回 粒子法は、ほかの方法とどう違うか
  • 第3回 粒子法の大きさと質量について
  • 第4回 「粒子法の動かし方」と「加速度の求め方」について​​
  • 第5回 計算時間を短縮する方法について
技術コラムTOPに戻る

 著者ご紹介

画像
東京大学大学院 工学系研究科 システム創成学専攻 准教授
プロメテック・ソフトウェア技術顧問
柴田 和也 先生
​<著者プロフィール>
2007年 東京大学大学院 工学系研究科
           システム量子工学専攻 博士課程修了 博士(工学)
2007年 (独)海上技術安全研究所 入所
           ​海の10モードプロジェクトチーム研究員
2009年 東京大学大学院 工学系研究科 システム創成学専攻 助教
2013年 東京大学大学院 工学系研究科 システム創成学専攻 講師
2017年 東京大学大学院 工学系研究科 システム創成学専攻 准教授
Picture
Picture
(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