JP6163897B2 - 数値計算プログラム、数値計算方法及び情報処理装置 - Google Patents
数値計算プログラム、数値計算方法及び情報処理装置 Download PDFInfo
- Publication number
- JP6163897B2 JP6163897B2 JP2013122330A JP2013122330A JP6163897B2 JP 6163897 B2 JP6163897 B2 JP 6163897B2 JP 2013122330 A JP2013122330 A JP 2013122330A JP 2013122330 A JP2013122330 A JP 2013122330A JP 6163897 B2 JP6163897 B2 JP 6163897B2
- Authority
- JP
- Japan
- Prior art keywords
- particles
- physical quantity
- particle
- density
- calculated
- 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.)
- Expired - Fee Related
Links
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Description
xij n+1/2=xi n+1/2−xj n+1/2
xij=xi−xj
vij=vi−vj
A)より具体的には、変数q+及びq-を以下のようにして求める。
B)また、変数q+及びq-の勾配を以下のように計算する。
C)変数q+及びq-について、粒子ij間の中間値を以下のように計算する。
第1の基準密度の流体に係る第1の粒子の密度の代わりに、当該第1の粒子の密度の、前記第1の基準密度に対する比率をリーマン不変量において用いることで得られる第1の物理量を算出し、
第2の基準密度の流体に係る第2の粒子の密度の代わりに、当該第2の粒子の密度の、前記第2の基準密度に対する比率をリーマン不変量において用いることで得られる第2の物理量を算出し、
前記第1の物理量の勾配及び前記第2の物理量の勾配を算出し、
前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値を、前記第1の物理量、前記第2の物理量、前記第1の物理量の勾配及び前記第2の物理量の勾配を用いて算出し、
前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する中間値から、前記第1の粒子と前記第2の粒子との間における圧力の時空間中間値を算出し、
前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する中間値から、前記第1の粒子と前記第2の粒子との間における速度の時空間中間値を算出する
第1の処理を、前記第1の粒子との位置関係から特定される第2の粒子の各々について実行し、
算出された前記圧力の時空間中間値を用いて前記第1の粒子の速度を更新し、
算出された前記速度の時空間中間値に基づき前記第1の粒子の密度を更新する
処理を、コンピュータに実行させるための数値計算プログラム。
前記勾配を算出する処理が、
前記第1の基準密度の空間勾配を用いて、前記第1の物理量の勾配を算出し、
前記第2の基準密度の空間勾配を用いて、前記第2の物理量の勾配を算出する
処理を含む付記1記載の数値計算プログラム。
前記第1の粒子の速度を更新する処理が、
前記圧力の時空間中間値から、前記第1の粒子の加速度を更新し、
前記第1の粒子の加速度を用いて、前記第1の粒子の速度を更新する
処理を含む付記1又は2記載の数値計算プログラム。
前記第1の粒子の密度を算出する処理が、
前記第1の粒子の速度と前記速度の時空間中間値とを用いて、前記第1の粒子の密度時間変化を算出し、
前記第1の粒子の密度時間変化を用いて、前記第1の粒子の密度を更新する
処理を含む付記1乃至3のいずれか1つ記載の数値計算プログラム。
第1の基準密度の流体に係る第1の粒子の密度の代わりに、当該第1の粒子の密度の、前記第1の基準密度に対する比率をリーマン不変量において用いることで得られる第1の物理量を算出し、
第2の基準密度の流体に係る第2の粒子の密度の代わりに、当該第2の粒子の密度の、前記第2の基準密度に対する比率をリーマン不変量において用いることで得られる第2の物理量を算出し、
前記第1の物理量の勾配及び前記第2の物理量の勾配を算出し、
前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値を、前記第1の物理量、前記第2の物理量、前記第1の物理量の勾配及び前記第2の物理量の勾配を用いて算出し、
前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する中間値から、前記第1の粒子と前記第2の粒子との間における圧力の時空間中間値を算出し、
前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する中間値から、前記第1の粒子と前記第2の粒子との間における速度の時空間中間値を算出する
第1の処理を、前記第1の粒子との位置関係から特定される第2の粒子の各々について実行し、
算出された前記圧力の時空間中間値を用いて前記第1の粒子の速度を更新し、
算出された前記速度の時空間中間値に基づき前記第1の粒子の密度を更新する
処理を含み、コンピュータにより実行される数値計算方法。
データ格納部と、
前記データ格納部に格納されているデータを用いて処理を行う処理部と、
を有し、
前記処理部は、
第1の基準密度の流体に係る第1の粒子の密度の代わりに、当該第1の粒子の密度の、前記第1の基準密度に対する比率をリーマン不変量において用いることで得られる第1の物理量を算出して、前記データ格納部に格納し、
第2の基準密度の流体に係る第2の粒子の密度の代わりに、当該第2の粒子の密度の、前記第2の基準密度に対する比率をリーマン不変量において用いることで得られる第2の物理量を算出して、前記データ格納部に格納し、
前記第1の物理量の勾配及び前記第2の物理量の勾配を算出して、前記データ格納部に格納し、
前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値を、前記第1の物理量、前記第2の物理量、前記第1の物理量の勾配及び前記第2の物理量の勾配を用いて算出して、前記データ格納部に格納し、
前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する中間値から、前記第1の粒子と前記第2の粒子との間における圧力の時空間中間値を算出して、前記データ格納部に格納し、
前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する中間値から、前記第1の粒子と前記第2の粒子との間における速度の時空間中間値を算出して、前記データ格納部に格納する
第1の処理を、前記第1の粒子との位置関係から特定される第2の粒子の各々について実行し、
算出された前記圧力の時空間中間値を用いて前記第1の粒子の速度を更新して、前記データ格納部に格納し、
算出された前記速度の時空間中間値に基づき前記第1の粒子の密度を更新して、前記データ格納部に格納する
処理を実行する情報処理装置。
110 初期条件データ格納部
120 処理部
130 データ格納部
140 処理結果格納部
150 出力部
Claims (6)
- 第1の基準密度の流体に係る第1の粒子の密度の代わりに、当該第1の粒子の密度の、前記第1の基準密度に対する比率をリーマン不変量において用いることで得られる第1の物理量を算出し、第2の基準密度の流体に係る第2の粒子の密度の代わりに、当該第2の粒子の密度の、前記第2の基準密度に対する比率をリーマン不変量において用いることで得られる第2の物理量を算出し、前記第1の物理量の勾配及び前記第2の物理量の勾配を算出し、前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値を、前記第1の物理量、前記第2の物理量、前記第1の物理量の勾配及び前記第2の物理量の勾配を用いて算出し、前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値から、前記第1の粒子と前記第2の粒子との間における圧力の時空間中間値を算出し、前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値から、前記第1の粒子と前記第2の粒子との間における速度の時空間中間値を算出する第1の処理を、前記第1の粒子との位置関係から特定される第2の粒子の各々について実行し、
算出された前記圧力の時空間中間値を用いて前記第1の粒子の速度を更新し、
算出された前記速度の時空間中間値に基づき前記第1の粒子の密度を更新する
処理を、コンピュータに実行させるための数値計算プログラム。 - 前記第1の物理量の勾配及び前記第2の物理量の勾配を算出する処理が、
前記第1の基準密度の空間勾配を用いて、前記第1の物理量の勾配を算出し、
前記第2の基準密度の空間勾配を用いて、前記第2の物理量の勾配を算出する
処理を含む請求項1記載の数値計算プログラム。 - 前記第1の粒子の速度を更新する処理が、
前記圧力の時空間中間値から、前記第1の粒子の加速度を更新し、
前記第1の粒子の加速度を用いて、前記第1の粒子の速度を更新する
処理を含む請求項1又は2記載の数値計算プログラム。 - 前記第1の粒子の密度を更新する処理が、
前記第1の粒子の速度と前記速度の時空間中間値とを用いて、前記第1の粒子の密度時間変化を算出し、
前記第1の粒子の密度時間変化を用いて、前記第1の粒子の密度を更新する
処理を含む請求項1乃至3のいずれか1つ記載の数値計算プログラム。 - 第1の基準密度の流体に係る第1の粒子の密度の代わりに、当該第1の粒子の密度の、前記第1の基準密度に対する比率をリーマン不変量において用いることで得られる第1の物理量を算出し、第2の基準密度の流体に係る第2の粒子の密度の代わりに、当該第2の粒子の密度の、前記第2の基準密度に対する比率をリーマン不変量において用いることで得られる第2の物理量を算出し、前記第1の物理量の勾配及び前記第2の物理量の勾配を算出し、前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値を、前記第1の物理量、前記第2の物理量、前記第1の物理量の勾配及び前記第2の物理量の勾配を用いて算出し、前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値から、前記第1の粒子と前記第2の粒子との間における圧力の時空間中間値を算出し、前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値から、前記第1の粒子と前記第2の粒子との間における速度の時空間中間値を算出する第1の処理を、前記第1の粒子との位置関係から特定される第2の粒子の各々について実行し、
算出された前記圧力の時空間中間値を用いて前記第1の粒子の速度を更新し、
算出された前記速度の時空間中間値に基づき前記第1の粒子の密度を更新する
処理を含み、コンピュータにより実行される数値計算方法。 - データ格納部と、
前記データ格納部に格納されているデータを用いて処理を行う処理部と、
を有し、
前記処理部は、
第1の基準密度の流体に係る第1の粒子の密度の代わりに、当該第1の粒子の密度の、前記第1の基準密度に対する比率をリーマン不変量において用いることで得られる第1の物理量を算出して、前記データ格納部に格納し、第2の基準密度の流体に係る第2の粒子の密度の代わりに、当該第2の粒子の密度の、前記第2の基準密度に対する比率をリーマン不変量において用いることで得られる第2の物理量を算出して、前記データ格納部に格納し、前記第1の物理量の勾配及び前記第2の物理量の勾配を算出して、前記データ格納部に格納し、前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値を、前記第1の物理量、前記第2の物理量、前記第1の物理量の勾配及び前記第2の物理量の勾配を用いて算出して、前記データ格納部に格納し、前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値から、前記第1の粒子と前記第2の粒子との間における圧力の時空間中間値を算出して、前記データ格納部に格納し、前記第1の粒子と前記第2の粒子との間における前記第1及び第2の物理量に関する時空間中間値から、前記第1の粒子と前記第2の粒子との間における速度の時空間中間値を算出して、前記データ格納部に格納する第1の処理を、前記第1の粒子との位置関係から特定される第2の粒子の各々について実行し、
算出された前記圧力の時空間中間値を用いて前記第1の粒子の速度を更新して、前記データ格納部に格納し、
算出された前記速度の時空間中間値に基づき前記第1の粒子の密度を更新して、前記データ格納部に格納する
処理を実行する情報処理装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013122330A JP6163897B2 (ja) | 2013-06-11 | 2013-06-11 | 数値計算プログラム、数値計算方法及び情報処理装置 |
US14/256,164 US20140365185A1 (en) | 2013-06-11 | 2014-04-18 | Numerical calculation method and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013122330A JP6163897B2 (ja) | 2013-06-11 | 2013-06-11 | 数値計算プログラム、数値計算方法及び情報処理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014241002A JP2014241002A (ja) | 2014-12-25 |
JP6163897B2 true JP6163897B2 (ja) | 2017-07-19 |
Family
ID=52006194
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013122330A Expired - Fee Related JP6163897B2 (ja) | 2013-06-11 | 2013-06-11 | 数値計算プログラム、数値計算方法及び情報処理装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20140365185A1 (ja) |
JP (1) | JP6163897B2 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6547547B2 (ja) * | 2015-09-25 | 2019-07-24 | 富士通株式会社 | 粒子シミュレーションプログラム、計算機資源配分方法、および粒子シミュレーション装置 |
CN105787162B (zh) * | 2015-11-23 | 2019-04-02 | 南京航空航天大学 | 基于非结构rkdg实现多介质界面追踪的数值模拟方法 |
CN116992796B (zh) * | 2023-09-27 | 2023-12-22 | 中国科学技术大学 | 一种自适应低耗散的sph-hllc黎曼求解器耦合方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012031398A1 (en) * | 2010-09-09 | 2012-03-15 | Tianjin Aerocode Engineering Application Software Development Inc. | Numerical method for simulating subsonic flows based on euler equations in lagrangian formulation |
WO2012111082A1 (ja) * | 2011-02-15 | 2012-08-23 | 富士通株式会社 | シミュレーション装置、シミュレーション方法、及びプログラム |
WO2013042234A1 (ja) * | 2011-09-21 | 2013-03-28 | 富士通株式会社 | 物体運動解析装置、物体運動解析方法、及び物体運動解析プログラム |
-
2013
- 2013-06-11 JP JP2013122330A patent/JP6163897B2/ja not_active Expired - Fee Related
-
2014
- 2014-04-18 US US14/256,164 patent/US20140365185A1/en not_active Abandoned
Also Published As
Publication number | Publication date |
---|---|
JP2014241002A (ja) | 2014-12-25 |
US20140365185A1 (en) | 2014-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lian et al. | Implementation of regularized isogeometric boundary element methods for gradient‐based shape optimization in two‐dimensional linear elasticity | |
Wang et al. | Algorithms for interface treatment and load computation in embedded boundary methods for fluid and fluid–structure interaction problems | |
Li | An overview of the immersed interface method and its applications | |
Alauzet et al. | 3D transient fixed point mesh adaptation for time-dependent problems: Application to CFD simulations | |
Oñate et al. | Lagrangian formulation for finite element analysis of quasi‐incompressible fluids with reduced mass losses | |
Marchandise et al. | A stabilized finite element method using a discontinuous level set approach for the computation of bubble dynamics | |
Liang et al. | An efficient staggered grid material point method | |
Pino Muñoz et al. | A finite element‐based level set method for fluid–elastic solid interaction with surface tension | |
Horgue et al. | A penalization technique applied to the “Volume-Of-Fluid” method: Wettability condition on immersed boundaries | |
Shamanskiy et al. | Mesh moving techniques in fluid-structure interaction: robustness, accumulated distortion and computational efficiency | |
Weiler et al. | Projective fluids | |
Baiges et al. | The Fixed‐Mesh ALE approach for the numerical simulation of floating solids | |
Cruchaga et al. | Finite element computation and experimental validation of sloshing in rectangular tanks | |
Evje et al. | Numerical treatment of two-phase flow in capillary heterogeneous porous media by finite-volume approximations | |
JP6163897B2 (ja) | 数値計算プログラム、数値計算方法及び情報処理装置 | |
Karakus et al. | An adaptive fully discontinuous Galerkin level set method for incompressible multiphase flows | |
Diehl et al. | Bond-based peridynamics: a quantitative study of mode i crack opening | |
Subber et al. | Asynchronous space–time algorithm based on a domain decomposition method for structural dynamics problems on non-matching meshes | |
Franci et al. | On the effect of the bulk tangent matrix in partitioned solution schemes for nearly incompressible fluids | |
De Corato et al. | Finite element formulation of fluctuating hydrodynamics for fluids filled with rigid particles using boundary fitted meshes | |
US20230008706A1 (en) | Particles-based fluid analysis simulation method using dummy particles, and fluid analysis simulation device | |
Khezri et al. | Application of RKP‐FSM in the buckling and free vibration analysis of thin plates with abrupt thickness changes and internal supports | |
Pu et al. | An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows | |
JP2021101329A (ja) | K−ω乱流モデルについての普遍的な壁面境界条件処理 | |
Zhang et al. | A moving mesh finite difference method for non-monotone solutions of non-equilibrium equations in porous media |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160310 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170314 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170511 |
|
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: 20170523 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170605 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6163897 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |