JP6098330B2 - 数値計算プログラム、数値計算方法及び情報処理装置 - Google Patents
数値計算プログラム、数値計算方法及び情報処理装置 Download PDFInfo
- Publication number
- JP6098330B2 JP6098330B2 JP2013088355A JP2013088355A JP6098330B2 JP 6098330 B2 JP6098330 B2 JP 6098330B2 JP 2013088355 A JP2013088355 A JP 2013088355A JP 2013088355 A JP2013088355 A JP 2013088355A JP 6098330 B2 JP6098330 B2 JP 6098330B2
- Authority
- JP
- Japan
- Prior art keywords
- particle
- particles
- internal energy
- continuum
- numerical calculation
- 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.)
- Active
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Power Engineering (AREA)
Description
vab=va−vb
rab(太字)=ra−rb
rab(イタリック)=|rab(太字)|
Tab=Ta−Tb
(A)粒子が連続体内部に存在する場合には放射冷却が起こらず、
(B)粒子が平面上に一様分布している場合には放射冷却が、当該粒子と等価な体積を有する、各辺が等長である平行六面体が有する六面のうち一面からの放射冷却が起こり、且つ
(C)粒子が孤立している場合には当該粒子と等価な体積の球と同じ放射冷却が起こる
ように補正するための連続関数である。なお、シュテファン・ボルツマンの法則は、(8)式において、係数ga以外の部分で表される。また、σは、物質に依存しないシュテファン・ボルツマン係数であり、αは物体の、黒体からのずれを表す係数(以下、放射率と呼ぶ)を表し、dは空間次元を表す。dは、1、2又は3である。さらに、(ma/ρa)1-(1/d)は、粒子aの表面積を表している。ここで、「当該粒子と等価な体積を有する、各辺が等長である平行六面体が有する六面のうち一面」は、当該粒子と等価な体積の立方体が有する六面のうち一面と等しい面積を有することになる。また、係数gaの満たすべき(B)の条件の代わりとして、表面形状に応じて粒子を、当該粒子と等価な体積の多面体や半球と仮定した場合において、当該粒子と等価な体積の立方体が有する六面のうち一面と等しい面積を有する立体であっても良い。
(A)粒子が物体内部に存在するときには放射冷却は起こらない。
これはga(0)=0と等価である。なぜならば、内部粒子の時Qa=0であり、この時ga(0)=0であれば、放射冷却によって失われる内部エネルギーが0となり、放射冷却が起きていないことになるからである。
(B)平面上に一様に分布しているとき、当該粒子と等価な体積の立方体のうち一面からの放射冷却が起こる(従来手法)。
これはga(1/2)=1と等価である。なぜならば、表面粒子の場合には、Qa=1/2であり、この時ga(1/2)=1であれば、放射冷却によって失われる内部エネルギーは、シュテファン・ボルツマンの法則に従って計算した量((8)式においてgaを除外した式の値)に一致するためである。
(C)孤立した粒子は等価な体積の球と同様に放射冷却が生ずる。
この条件は、以下のような式で表される。
連続体を粒子の集合として表現する場合における各粒子について、当該粒子が前記連続体表面にさらされている度合いに対応する係数に応じ且つ放射冷却に基づく内部エネルギーの時間微分を算出し、
各粒子について算出された前記内部エネルギーの時間微分に基づき、単位時間後の前記内部エネルギーを算出する
処理を、コンピュータに実行させるための数値計算プログラム。
前記係数が、
粒子が前記連続体内部に存在する場合には放射冷却が起こらず、
粒子が平面上に一様分布している場合には当該粒子と等価な体積を有する、各辺が等長である平行六面体が有する六面のうち一面からの放射冷却が起こり、且つ
粒子が孤立している場合には当該粒子と等価な体積の球と同じ放射冷却が起こる
ように補正するための連続関数である
付記1記載の数値計算プログラム。
前記連続関数が、規格化された数密度saの関数となっている
付記2記載の数値計算プログラム。
連続体を粒子の集合として表現する場合における各粒子について、当該粒子が前記連続体表面にさらされている度合いに対応する係数に応じ且つ放射冷却に基づく内部エネルギーの時間微分を算出し、
各粒子について算出された前記内部エネルギーの時間微分に基づき、単位時間後の前記内部エネルギーを算出する
処理を含み、コンピュータにより実行される数値計算方法。
連続体を粒子の集合として表現する場合における各粒子について、当該粒子が前記連続体表面にさらされている度合いに対応する係数に応じ且つ放射冷却に基づく内部エネルギーの時間微分を算出する手段と、
各粒子について算出された前記内部エネルギーの時間微分に基づき、単位時間後の前記内部エネルギーを算出する手段と、
を有する情報処理装置。
120 第1データ格納部
130 物理量算出部
140 第2データ格納部
150 出力部
131 近傍リスト生成部
132 放射冷却算出部
133 積分処理部
Claims (6)
- 連続体を粒子の集合として表現する場合における各粒子について、当該粒子が前記連続体表面にさらされている度合いに対応する係数に応じ且つ放射冷却に基づく内部エネルギーの時間微分を算出し、
各粒子について算出された前記内部エネルギーの時間微分に基づき、単位時間後の前記内部エネルギーを算出する
処理を、コンピュータに実行させるための数値計算プログラムであって、
前記係数が、
粒子が前記連続体内部に存在する場合には放射冷却が起こらず、
粒子が平面上に一様分布している場合には当該粒子と等価な体積を有する、各辺が等長である平行六面体が有する六面のうち一面からの放射冷却が起こり、且つ
粒子が孤立している場合には当該粒子と等価な体積の球と同じ放射冷却が起こる
ように補正するための連続関数である
数値計算プログラム。 - 連続体を粒子の集合として表現する場合における各粒子について、当該粒子が前記連続体表面にさらされている度合いに対応する係数に応じ且つ放射冷却に基づく内部エネルギーの時間微分を算出し、
各粒子について算出された前記内部エネルギーの時間微分に基づき、単位時間後の前記内部エネルギーを算出する
処理を含み、
前記係数が、
粒子が前記連続体内部に存在する場合には放射冷却が起こらず、
粒子が平面上に一様分布している場合には当該粒子と等価な体積を有する、各辺が等長である平行六面体が有する六面のうち一面からの放射冷却が起こり、且つ
粒子が孤立している場合には当該粒子と等価な体積の球と同じ放射冷却が起こる
ように補正するための連続関数である、
コンピュータにより実行される数値計算方法。 - 連続体を粒子の集合として表現する場合における各粒子について、当該粒子が前記連続体表面にさらされている度合いに対応する係数に応じ且つ放射冷却に基づく内部エネルギーの時間微分を算出する手段と、
各粒子について算出された前記内部エネルギーの時間微分に基づき、単位時間後の前記内部エネルギーを算出する手段と、
を有し、
前記係数が、
粒子が前記連続体内部に存在する場合には放射冷却が起こらず、
粒子が平面上に一様分布している場合には当該粒子と等価な体積を有する、各辺が等長である平行六面体が有する六面のうち一面からの放射冷却が起こり、且つ
粒子が孤立している場合には当該粒子と等価な体積の球と同じ放射冷却が起こる
ように補正するための連続関数である、
情報処理装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013088355A JP6098330B2 (ja) | 2013-04-19 | 2013-04-19 | 数値計算プログラム、数値計算方法及び情報処理装置 |
US14/191,919 US20140316596A1 (en) | 2013-04-19 | 2014-02-27 | Information processing method and information processing system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013088355A JP6098330B2 (ja) | 2013-04-19 | 2013-04-19 | 数値計算プログラム、数値計算方法及び情報処理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014211798A JP2014211798A (ja) | 2014-11-13 |
JP6098330B2 true JP6098330B2 (ja) | 2017-03-22 |
Family
ID=51729627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013088355A Active JP6098330B2 (ja) | 2013-04-19 | 2013-04-19 | 数値計算プログラム、数値計算方法及び情報処理装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20140316596A1 (ja) |
JP (1) | JP6098330B2 (ja) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7197404B2 (en) * | 2004-03-01 | 2007-03-27 | Richard Andrew Holland | Computation of radiating particle and wave distributions using a generalized discrete field constructed from representative ray sets |
JP4278610B2 (ja) * | 2004-12-28 | 2009-06-17 | 富士通株式会社 | 数値解析支援装置,数値解析支援方法,数値解析支援プログラムおよび同プログラムを記録したコンピュータ読取可能な記録媒体 |
US7473028B1 (en) * | 2005-04-22 | 2009-01-06 | The Ohio State University | Method and device for investigation of phase transformations in metals and alloys |
WO2007072863A1 (ja) * | 2005-12-20 | 2007-06-28 | Sintokogio, Ltd. | 投射機による投射条件情報の推定方法及びその装置 |
US8783337B2 (en) * | 2006-12-01 | 2014-07-22 | The Invention Science Fund I Llc | System for changing the convective heat transfer coefficient for a surface |
US8155939B2 (en) * | 2008-06-12 | 2012-04-10 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Particle-surface interaction model and method of determining particle-surface interactions |
IL197176A0 (en) * | 2009-02-23 | 2009-12-24 | Yuli Lozinski Dr | New heat flow measuring system |
-
2013
- 2013-04-19 JP JP2013088355A patent/JP6098330B2/ja active Active
-
2014
- 2014-02-27 US US14/191,919 patent/US20140316596A1/en not_active Abandoned
Also Published As
Publication number | Publication date |
---|---|
JP2014211798A (ja) | 2014-11-13 |
US20140316596A1 (en) | 2014-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6657359B2 (ja) | ハイブリッド熱格子ボルツマン法のための温度結合アルゴリズム | |
Adami et al. | A new surface-tension formulation for multi-phase SPH using a reproducing divergence approximation | |
Ding et al. | Numerical computation of three-dimensional incompressible viscous flows in the primitive variable form by local multiquadric differential quadrature method | |
US20190368344A1 (en) | Mass exchange model for relative permeability simulation | |
JP5724814B2 (ja) | 熱流体シミュレーションプログラム,熱流体シミュレーション装置および熱流体シミュレーション方法 | |
JP6098190B2 (ja) | シミュレーションプログラム、シミュレーション方法及びシミュレーション装置 | |
US20200394277A1 (en) | Computer simulation of physical fluids on irregular spatial grids stabilized for explicit numerical diffusion problems | |
JP5704246B2 (ja) | 物体運動解析装置、物体運動解析方法、及び物体運動解析プログラム | |
JP2008152423A (ja) | 粒子モデルを用いた変形挙動シミュレーション方法およびプログラム | |
JP5892257B2 (ja) | シミュレーションプログラム、シミュレーション方法及びシミュレーション装置 | |
Ortega et al. | A meshless finite point method for three‐dimensional analysis of compressible flow problems involving moving boundaries and adaptivity | |
TWI614687B (zh) | 構造解析方法及構造解析程式 | |
JP6098330B2 (ja) | 数値計算プログラム、数値計算方法及び情報処理装置 | |
Cruchaga et al. | A surface remeshing technique for a Lagrangian description of 3D two‐fluid flow problems | |
JP6163897B2 (ja) | 数値計算プログラム、数値計算方法及び情報処理装置 | |
JP6808195B2 (ja) | 流体シミュレーションプログラム、流体シミュレーション装置および流体シミュレーション方法 | |
JP5761355B2 (ja) | 運動解析装置、運動解析方法及び運動解析プログラム | |
US20240054268A1 (en) | Methods and systems for physics-based reduced-order modeling of local dynamics in additive manufacturing | |
Shahane et al. | Virtually-guided certification with uncertainty quantification applied to die casting | |
JP7395456B2 (ja) | シミュレーション装置、及びプログラム | |
Pino Muñoz et al. | Direct 3D simulation of powder sintering by surface and volume diffusion | |
Vakhrushev et al. | Modeling of turbulent melt flow and solidification processes in steel continuous caster with the open source software package openFOAM | |
JP2019101939A (ja) | 流体中の繊維状物質の運動状態の解析方法及びその解析装置 | |
JP4145078B2 (ja) | 鋳造シミュレーション方法 | |
Bharadwaj et al. | An Eulerian meshless method for two-phase flows with embedded geometries |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160113 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20161206 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20161227 |
|
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: 20170124 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170206 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6098330 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |