JP2010049561A - 流動解析方法、流動解析装置、及び流動解析プログラム - Google Patents
流動解析方法、流動解析装置、及び流動解析プログラム Download PDFInfo
- Publication number
- JP2010049561A JP2010049561A JP2008214260A JP2008214260A JP2010049561A JP 2010049561 A JP2010049561 A JP 2010049561A JP 2008214260 A JP2008214260 A JP 2008214260A JP 2008214260 A JP2008214260 A JP 2008214260A JP 2010049561 A JP2010049561 A JP 2010049561A
- Authority
- JP
- Japan
- Prior art keywords
- flow analysis
- particle
- particles
- incompressible fluid
- behavior
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 90
- 238000005206 flow analysis Methods 0.000 title claims description 64
- 239000002245 particle Substances 0.000 claims abstract description 124
- 239000012530 fluid Substances 0.000 claims abstract description 31
- 238000004364 calculation method Methods 0.000 description 32
- 238000004458 analytical method Methods 0.000 description 19
- 238000004590 computer program Methods 0.000 description 6
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000010365 information processing Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
【解決手段】解像度が必要な領域内における粒子の粒子径を小さくする一方で解像度が必要でない領域内における粒子の粒子径を大きくすることにより、解析領域内で解像度を変化させる。これにより、解析領域内でスケールが大きく変化する非圧縮流体の挙動解析に要する時間を大幅に短縮することができる。
【選択図】図1
Description
本発明の実施形態となる流動解析装置1は、パーソナルコンピュータ,ワークステーション,汎用コンピュータ等の公知の情報処理装置により構成され、図1に示すように、RAM(Random Access Memory)2,ROM(Read Only Memory)3,初期設定データベース(DB)4,及びCPU(Central Processing Unit)5を主な構成要素として備える。RAM2は、CPU5が実行する流動解析装置1の起動プログラムや流動解析プログラム等のコンピュータプログラム及び各種データを一時的に格納するワーキングエリアを提供する。ROM3は、コンピュータプログラム及びこのコンピュータプログラムを実行するために必要な各種データを記憶する。なおROM3は、磁気的又は光学的な記憶媒体若しくは半導体メモリ等のCPU5により読み取り可能な記憶媒体を含む構成を有し、ROM3に記憶されているコンピュータプログラムや各種データの一部又は全部は電気通信回線を介してダウンロードされるように構成してもよい。
このような構成を有する流動解析装置1は、以下に示す流動解析処理を実行することによりスケールが大きく変化する解析領域内における非圧縮流体の挙動を現実的な時間で解析する。以下、図2に示すフローチャートを参照して、この流動解析処理を実行する際の流動解析装置1の動作について詳しく説明する。
本願発明の発明者は、従来のMPS法を利用した流動解析手法において粒子径が異なる複数の粒子により非圧縮流体を表現した場合、図3に示すように圧力値が空間的に振動することを確認した。そして鋭意研究を重ねてきた結果、数式7に示す重みを用いた場合には、図4に示すように滑らかな圧力分布が得られることを知見した。また本願発明の発明者は、粒子間で粒子径が1.5倍以上異なる場合には、図5(a)〜(d)に示すように計算が不安定になることを確認した。そして鋭意研究を重ねてきた結果、速度の発散が0であるとの非圧縮条件を与えて流動解析を行うことにより、図6(a)〜(c)に示すように計算の安定性が劇的に向上することを知見した。以上のことから、本願発明の発明者は、数式7に示す重みを用いて、速度の発散が0であるとの非圧縮条件を与えて流動解析を行うことにより、解析領域内で解像度を変化させることができ、結果として、スケールが大きく変化する非圧縮流体の挙動解析に要する時間を大幅に短縮することができた。
最後に、ダム崩壊問題の解析により本手法の優位性を示す。図7は、4種類の大きさの粒子が混在しているモデルである。粒子の種類を2種類から4種類まで変化させて解析を行ったところ、必要とした粒子数及び計算時間は図8に示すようになった。図8から明らかなように、粒子径を大きく変化させる程、計算効率が高くなり、計算時間は1/3以下となった。また容器右下の点における圧力値の変化を測定した所、図9に示すように粒子径が均一である場合(均一粒子)と粒子径が均一でない場合(非均一粒子)とではほぼ同じ結果になることが確認された。以上のことから、本実施形態の流動解析処理によれば、必要な部分の精度を保ったまま計算コストを下げることが可能であるといえる。
Claims (6)
- 非圧縮流体の挙動を解析するための流動解析方法であって、
MPS法により定義される粒子を利用して前記非圧縮流体を粒子径が異なる複数の粒子により表現して非圧縮流体の挙動を解析することを特徴とする流動解析方法。 - 請求項1に記載の流動解析方法において、速度の発散が0であることを非圧縮条件として非圧縮流体の挙動を解析することを特徴とする流動解析方法。
- 請求項2に記載の流動解析方法において、Vi及びVjをそれぞれ粒子i及び粒子jの体積、ri及びrjをそれぞれ粒子i及び粒子jの位置座標、Ri及びRjをそれぞれ粒子i及び粒子jの影響半径とした時、以下の数式1で表される粒子iと粒子j間の重みwijを用いて非圧縮流体の挙動を解析することを特徴とする流動解析方法。
- 請求項3に記載の流動解析方法において、前記複数の粒子を合体及び分裂させることにより非圧縮流体の解像度を空間的に変化させることを特徴とする流動解析方法。
- 非圧縮流体の挙動を解析するための流動解析装置であって、
MPS法により定義される粒子を利用して前記非圧縮流体を粒子径が異なる複数の粒子により表現する初期データ設定手段を備えることを特徴とする流動解析装置。 - 非圧縮流体の挙動を解析するための流動解析プログラムであって、
MPS法により定義される粒子を利用して前記非圧縮流体を粒子径が異なる複数の粒子により表現する処理をコンピュータに実行させることを特徴とする流動解析プログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008214260A JP5268496B2 (ja) | 2008-08-22 | 2008-08-22 | 流動解析方法、流動解析装置、及び流動解析プログラム |
US12/544,595 US8296112B2 (en) | 2008-08-22 | 2009-08-20 | Flow simulation method, flow simulation system, and computer program product |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008214260A JP5268496B2 (ja) | 2008-08-22 | 2008-08-22 | 流動解析方法、流動解析装置、及び流動解析プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010049561A true JP2010049561A (ja) | 2010-03-04 |
JP5268496B2 JP5268496B2 (ja) | 2013-08-21 |
Family
ID=41697164
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008214260A Expired - Fee Related JP5268496B2 (ja) | 2008-08-22 | 2008-08-22 | 流動解析方法、流動解析装置、及び流動解析プログラム |
Country Status (2)
Country | Link |
---|---|
US (1) | US8296112B2 (ja) |
JP (1) | JP5268496B2 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013038476A1 (ja) * | 2011-09-12 | 2013-03-21 | 富士通株式会社 | 運動解析装置、運動解析方法及び運動解析プログラム |
KR20150048960A (ko) * | 2013-10-28 | 2015-05-11 | 삼성전자주식회사 | 입자에 기반한 모델링 방법 및 장치 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101731968B1 (ko) * | 2010-11-01 | 2017-05-04 | 삼성전자주식회사 | 로봇의 리로케이션 장치 및 방법 |
JP5919784B2 (ja) * | 2011-12-09 | 2016-05-18 | 富士通株式会社 | 設計支援方法、設計支援プログラムおよび設計支援装置 |
US9070304B2 (en) * | 2012-03-28 | 2015-06-30 | Korea Institute Of Geoscience And Mineral Resources | Debris-flow simulation apparatus having variable flume |
KR101195409B1 (ko) * | 2012-04-10 | 2012-11-05 | 한국지질자원연구원 | 수막현상 재현형 토석류 모형시험장치 |
EP3525120A1 (en) * | 2013-07-24 | 2019-08-14 | Dassault Systemes Simulia Corp. | Lattice boltzmann collision operators enforcing isotropy and galilean invariance |
US10210287B2 (en) | 2013-12-31 | 2019-02-19 | Disney Enterprises, Inc. | Augmented material point method for simulating phase changes and varied materials |
US10019827B2 (en) * | 2013-12-31 | 2018-07-10 | Disney Enterprises, Inc. | Material point method for simulation of granular materials |
CN113807034B (zh) * | 2021-08-30 | 2023-05-16 | 西安交通大学 | 基于移动粒子半隐式法的轴对称流场二维模拟方法 |
CN115345088B (zh) * | 2022-08-22 | 2024-02-23 | 深圳十沣科技有限公司 | 多相流耦合计算方法、装置、设备及存储介质 |
CN118194744A (zh) * | 2024-03-07 | 2024-06-14 | 山东大学 | 考虑温度效应的富水破碎岩体cfd-dem注浆模拟方法及系统 |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07334484A (ja) | 1994-06-13 | 1995-12-22 | Seiichi Koshizuka | 流動解析装置及び流動解析方法 |
US6198489B1 (en) * | 1997-02-21 | 2001-03-06 | University Of Washington | Computer generated watercolor |
JP2000339292A (ja) | 1999-05-26 | 2000-12-08 | Canon Inc | シミュレーション方法及び装置 |
US20060169800A1 (en) * | 1999-06-11 | 2006-08-03 | Aradigm Corporation | Aerosol created by directed flow of fluids and devices and methods for producing same |
JP2002137272A (ja) | 2000-11-02 | 2002-05-14 | Mitsubishi Electric Corp | 樹脂成形シミュレーション方法およびその装置ならびに記録媒体 |
US7239990B2 (en) * | 2003-02-20 | 2007-07-03 | Robert Struijs | Method for the numerical simulation of a physical phenomenon with a preferential direction |
US7117138B2 (en) * | 2003-03-14 | 2006-10-03 | Seiko Epson Corporation | Coupled quadrilateral grid level set scheme for piezoelectric ink-jet simulation |
US7526415B2 (en) * | 2004-06-30 | 2009-04-28 | D. E. Shaw Research, Llc | Grid based computation for multiple body simulation. |
US7813903B2 (en) * | 2005-04-13 | 2010-10-12 | Autodesk, Inc. | Fixed time step dynamical solver for interacting particle systems |
US20060241928A1 (en) * | 2005-04-25 | 2006-10-26 | International Business Machines Corporation | Load balancing by spatial partitioning of interaction centers |
US20070233440A1 (en) * | 2006-03-29 | 2007-10-04 | International Business Machines Corporation | Reduced message count for interaction decomposition of N-body simulations |
US7565276B2 (en) * | 2006-04-05 | 2009-07-21 | Seoul National University Industry Foundation | Method of simulating detailed movements of fluids using derivative particles |
US7558714B2 (en) * | 2006-08-10 | 2009-07-07 | Exa Corporation | Computer simulation of physical processes |
US7921003B2 (en) * | 2007-01-23 | 2011-04-05 | Adobe Systems Incorporated | System and method for simulating shallow water effects on arbitrary surfaces |
CA2648441A1 (en) * | 2007-12-31 | 2009-06-30 | Exocortex Technologies, Inc. | Fast characterization of fluid dynamics |
US7930155B2 (en) * | 2008-04-22 | 2011-04-19 | Seiko Epson Corporation | Mass conserving algorithm for solving a solute advection diffusion equation inside an evaporating droplet |
US7949503B2 (en) * | 2008-07-07 | 2011-05-24 | King Fahd University Of Petroleum And Minerals | Facilities optimization method |
US8346711B2 (en) * | 2009-11-24 | 2013-01-01 | King Fahd University Of Petroleum And Minerals | Method for identifying multi-input multi-output Hammerstein models |
US20120065947A1 (en) * | 2010-09-09 | 2012-03-15 | Jiun-Der Yu | Collision Effect And Particle Information Update In Particulate Fluid Flow Simulations |
-
2008
- 2008-08-22 JP JP2008214260A patent/JP5268496B2/ja not_active Expired - Fee Related
-
2009
- 2009-08-20 US US12/544,595 patent/US8296112B2/en not_active Expired - Fee Related
Non-Patent Citations (5)
Title |
---|
CSNG200701210001; 越塚 誠一: '第4回 MPS粒子法による非圧縮性流れ/構造解析' シミュレーション Vol.24 No.3, 20050915, pp.35-42 * |
JPN6013000820; 豊田 英司: '空間解像度を可変とした粒子法による高速滑走体の流体数値シミュレーション' 日本船舶海洋工学会講演会論文集 No.4, 200705, pp.321-322 * |
JPN6013000823; 越塚 誠一: '第4回 MPS粒子法による非圧縮性流れ/構造解析' シミュレーション Vol.24 No.3, 20050915, pp.35-42 * |
JPN6013000825; 後藤 仁志: '砕波CGの写実性向上のための粒子法SPS表情モデルの構築' 海岸工学論文集 Vol.51, 20041008, pp.116-120 * |
JPN6013000827; 豊田 英司: '解像度可変型粒子法による水上滑走艇周りの流れのシミュレーション' 日本流体力学会年会2006講演アブストラクト集 , 20060905, p.57 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013038476A1 (ja) * | 2011-09-12 | 2013-03-21 | 富士通株式会社 | 運動解析装置、運動解析方法及び運動解析プログラム |
JPWO2013038476A1 (ja) * | 2011-09-12 | 2015-03-23 | 富士通株式会社 | 運動解析装置、運動解析方法及び運動解析プログラム |
KR20150048960A (ko) * | 2013-10-28 | 2015-05-11 | 삼성전자주식회사 | 입자에 기반한 모델링 방법 및 장치 |
KR102205845B1 (ko) * | 2013-10-28 | 2021-01-22 | 삼성전자주식회사 | 입자에 기반한 모델링 방법 및 장치 |
Also Published As
Publication number | Publication date |
---|---|
US8296112B2 (en) | 2012-10-23 |
US20100049489A1 (en) | 2010-02-25 |
JP5268496B2 (ja) | 2013-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5268496B2 (ja) | 流動解析方法、流動解析装置、及び流動解析プログラム | |
JP2008111675A (ja) | 流体解析装置、流体解析方法及び流体解析プログラム | |
JP2019125102A (ja) | 流体解析装置、流体解析方法、及び流体解析プログラム | |
JP2016024178A (ja) | 特定物質の解析用モデルの作成方法、特定物質の解析用モデルの作成用コンピュータプログラム、特定物質のシミュレーション方法及び特定物質のシミュレーション用コンピュータプログラム | |
Ghobadi et al. | Adapting SAFT-γ perturbation theory to site-based molecular dynamics simulation. II. Confined fluids and vapor-liquid interfaces | |
Guan et al. | Viscous fingering in a Hele-Shaw cell with finite viscosity ratio and interfacial tension | |
JP2016024787A (ja) | 特定物質の解析用モデルの作成方法、特定物質の解析用モデルの作成用コンピュータプログラム、特定物質のシミュレーション方法及び特定物質のシミュレーション用コンピュータプログラム | |
JP2017138123A (ja) | 気液二相流のシミュレーション方法 | |
JP5242224B2 (ja) | 流動解析方法、流動解析装置、及び流動解析プログラム | |
KR101652210B1 (ko) | 해석장치 | |
JP2006242893A (ja) | 地下水流動の到達範囲を求める数値解析方法及び数値解析プログラム | |
JP2017084081A (ja) | 磁性体シミュレーション装置、マイクロ磁化算出方法及びプログラム | |
Gera et al. | Three-dimensional multicomponent vesicles: dynamics and influence of material properties | |
JP5669589B2 (ja) | 解析装置 | |
JP2010122946A (ja) | 半導体デバイスモデルのモデルパラメータ抽出装置およびモデルパラメータ抽出プログラム | |
Zeifang et al. | An investigation of different splitting techniques for the isentropic Euler equations | |
Kuznetsov et al. | Estimating the domain of absolute stability of a numerical scheme based on the method of solution continuation with respect to a parameter for solving stiff initial value problems | |
JP5156516B2 (ja) | 材料特性の推定方法及び材料特性の推定プログラム | |
CN110244369A (zh) | 储层束缚和可动流体分布确定方法、装置及系统 | |
KR101807427B1 (ko) | 곡면을 구비한 노치의 최대 응력값 결정 방법 | |
KR102364390B1 (ko) | 매트랩을 활용한 모델링 자동화 장치 및 방법 | |
JP6183359B2 (ja) | 設計支援装置、設計支援方法及びプログラム | |
JP7285223B2 (ja) | シミュレーション装置、シミュレーション方法、及びプログラム | |
KR102478451B1 (ko) | 연산 장치 및 이를 이용한 시스템의 상태 추정 방법 | |
JP2013105259A (ja) | 流路断面形状の設計装置、その方法及びそのプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100924 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130115 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130318 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20130409 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130507 |
|
LAPS | Cancellation because of no payment of annual fees |