Project/Area Number |
13304007
|
Research Category |
Grant-in-Aid for Scientific Research (A)
|
Allocation Type | Single-year Grants |
Section | 一般 |
Research Field |
General mathematics (including Probability theory/Statistical mathematics)
|
Research Institution | Kyushu University |
Principal Investigator |
TABATA Masahisa Kyushu University, Faculty of Mathematics, Professor, 大学院・数理学研究院, 教授 (30093272)
|
Co-Investigator(Kenkyū-buntansha) |
KANAYAMA Hiroshi Kyushu University, Faculty of Engineering, Professor, 大学院・工学研究院, 教授 (90294884)
USHIJIMA Teruo Univ.of Electro-Communications, Professor, 電気通信学部, 教授 (10012410)
IMAI Hitoshi Tokushima Univ., Faculty of Engineering, Professor, 工学部, 教授 (80203298)
NAKAO Mitsuhiro Kyushu University, Faculty of Mathematics, Professor, 大学院・数理学研究院, 教授 (10136418)
KIKUCHI Fumio Univ.of Tokyo, Graduate School of Mathematical Sciences, Professor, 大学院・数理科学研究科, 教授 (40013734)
友枝 謙二 大阪工業大学, 工学部, 教授 (60033916)
大塚 厚二 広島国際学院大学, 工学部, 教授 (30141683)
|
Project Period (FY) |
2001 – 2003
|
Project Status |
Completed (Fiscal Year 2003)
|
Budget Amount *help |
¥40,820,000 (Direct Cost: ¥31,400,000、Indirect Cost: ¥9,420,000)
Fiscal Year 2003: ¥12,480,000 (Direct Cost: ¥9,600,000、Indirect Cost: ¥2,880,000)
Fiscal Year 2002: ¥12,610,000 (Direct Cost: ¥9,700,000、Indirect Cost: ¥2,910,000)
Fiscal Year 2001: ¥15,730,000 (Direct Cost: ¥12,100,000、Indirect Cost: ¥3,630,000)
|
Keywords | finite element method / error estimate / numerical simulation / parallel computation / heat convection problems / moving boundary problems / validated computation / Nedelec's edge elements / 温度依存粘性 / 無限精度数値計算 / ヘルムホルツ方程式 / 板曲げ要素 / 非線形拡散 / 特性有限要素法 / 特異極限法 / 精度保証付き数値計算法 / 多孔質媒体流に対する差分法 / ライスナー・ミンドリン型要素 / 自由境界 / 円外帰着波動問題の基本解近似解法 / 一様可解性 / 地球マントル対流 |
Research Abstract |
1.In devising numerical schemes for flow problems, how to approximate the convection torn is a crucial point. Characteristic finite element approximation is based on the approximation of the material derivative, which is the sum of the time derivative term and the convection term. So far finite element schemes of characteristic method of the first-order accuracy in time increment have been used. We have developed a finite element scheme of the second-order accuracy in time increment and obtained the best possible error estimate. This scheme is more robust than the first-order scheme with respect to numerical integration error and can solve flow problems more stably and accurately. 2.We have developed a finite element scheme and established an error estimate for heat convection problems with temperature-dependent viscosity. The viscosity of heat conduction problems such as mantle convection in the Earth and melting glass convection in the furnace is strongly dependent on the temperature.
… More
The dependence plays an important role in die formation of convection patterns. Our scheme is applicable for the general Rayleigh-Benard problems with temperature-dependent viscosity, thermal conductivity, and thermal expansion coefficient. Using this scheme we have carried out large-scale numerical simulation of Earth's mantle convection in three-dimensional spherical shell and succeeded in obtaining complex heat convection patterns. 3.In the infinite precision computation we have succeeded a large-scale parallel computation using a cluster of high-performance computers with 10CPU and 20GB memory. For one-dimensional boundary-value problems very precise results with precision 4995 digits have been obtained. We have used this system to perform direct numerical simulation of inverse problems and made possible a numerical analysis of inverse problems. 4.Formulating eddy current problems in magnetic vector potential and electric scalar potential, we have solved them using a hierarchical domain decomposition method. This solution has been shown to be effective under the environment of parallel computation. By this method we have carried out large-scale numerical simulation of nonlinear static magnetic problems in magnetic vector potential. Less
|