权益分类	功能权益	普通用户	{{item.name}}会员
{{category.name}}	{{benefitItem.name}}

衝撃波・岩石破砕を含む噴火シミュレータのための高精度数値計算スキームの開発

开发包括冲击波和岩石破裂在内的喷发模拟器高精度数值计算方案

基本信息

批准号：
17031008
负责人：
青木尊之
金额：
$ 2.62万
依托单位：
Tokyo Institute of Technology
依托单位国家：
日本
项目类别：
Grant-in-Aid for Scientific Research on Priority Areas
财政年份：
2005
资助国家：
日本
起止时间：
2005 至 2006
项目状态：
已结题

项目摘要

「噴火シミュレータ」をさらに高精度化するための以下の数値計算法の開発を行った。火山噴火は爆風の衝撃波伝播とかなり類似した部分があり圧縮性流体計算を必要とする。非常に狭い領域で大きなエネルギーを短時間に放出し、その影響が広範囲に及ぶ。数値計算手法に求められることとして、噴火口付近の領域に十分細かい計算格子を配置する、密度・圧力の急激な変化に対して安定に計算できるスキームである、長時間積分が可能である、ことが必要である。これらの全てを満足する数値計算手法として、保存形の局所補間微分オペレータ(IDO)法をさらに発展させた。(1)従来の物理量の値と微係数を従属変数として補間関数を構築するIDO法を値と区間積分値を従属変数として補間関数を構築する方式に変更し、区間積分値に対しては偏微分方程式を区間積分して得られる方程式で時間積分を行い、flux formの保存形IDO法を構築した。(2)保存形IDO法の精度検証を行い、風上補間の場合には補間関数は2次多項式になるにも関わらず、移流計算の精度は3次精度であることが明らかになった。また、中心補間に対しては4次多項式になり4次精度がでることが分かった。非保存形IDO法と比較して補間関数の多項式の次数は下がるが、フーリエ解析の結果はgainおよびphaseとも非保存形と区別できないほど一致した結果が得られた。(3)噴火の最中は火口付近の空間解像度を十分高くし、噴火終了後は噴煙を精度よく捕らえるために上空の空間解像度をAdaptiveに高くするAMR (Adaptive Mesh Refinement)法にIDO法を検討した。非保存形IDO法の場合と異なり、Cellのマージ過程に補間操作が無くなり、積分値は単純な足し算となるためこの部分の誤差がゼロであることが分かった。(4)岩石のモデル化として、複雑な形状を反映するために非球形モデルの個別要素法をさらに発展させ、球形粒子を複雑形状の岩石を多数(1000個)配置し、斜面を転がり落ちる落石過程をより高精度化した。また、岩石の慣性モーメントもより忠実に計算する方法に改善した。また、岩石破壊過程のモデルを検討し、計算量の問題から岩石全体の構成粒子を要素バネで連結することは困難であることが判明し、剛体連結したクラスター粒子の間を断裂するバネで連結し落石計算を行った。

The development of numerical value calculation method for high precision of "spray fire" is carried out. Volcanic eruption and blast wave propagation are necessary for calculation of compressible fluids. In a very narrow field, the large number of people in a short time, the impact of the large number of people and the small number of people. The calculation method of numerical value is very detailed in the field close to the nozzle. The calculation grid is arranged. The calculation method of density and pressure is very detailed. The calculation method of long time integration is possible and necessary. The method of calculating the total number of points and preserving the shape of the partial differential equation (IDO) is developed. (1)The IDO method for constructing interval integral values, interval integral values, interval integral values, (2)The accuracy of the preservation IDO method is verified in the middle of the line, and in the case of the wind supplement, the supplement is related to the second order polynomial, and the accuracy of the flow calculation is 3 degrees. The center of the space is the fourth order polynomial. Non-preserving IDO method and comparison of the degree of the polynomial of the number of complement, the result of the analysis, the result of the phase difference, the result of the consistency (3)The highest spatial resolution of the flame is very high, and the highest spatial resolution of the flame is very high. The AMR (Adaptive Mesh Refining) method is discussed after the end of the flame. The occasions of the non-saving IDO method are different, the interpolation operation in the Cell process is irrelevant, and the integral value is pure and sufficient. (4)The individual element method for the development of non-spherical particles, the arrangement of a large number of spherical particles, the formation of complex rocks, and the accurate process of falling rocks are described. The method of calculation is improved by the inertia of rock. The problem of rock fracture process is discussed and calculated. The problem of rock fracture process is difficult to determine. The problem of rock fracture process is difficult to determine. The problem of rock fracture process is difficult to determine.