偏微分方程式の無限精度数値シミュレーションに関する研究

2000 Fiscal Year Annual Research Report

• Project/Area Number: 10354001
• Research Institution: The University of Tokushima
• Principal Investigator: 今井 仁司 徳島大学, 工学部, 教授 (80203298)
• Co-Investigator(Kenkyū-buntansha): 友枝 謙二 大阪工業大学, 工学部, 教授 (60033916) ; 田端 正久 九州大学, 大学院・数理学研究院, 教授 (30093272) ; 池田 勉 龍谷大学, 理工学部, 教授 (50151296) ; 西田 孝明 京都大学, 大学院・理学研究科, 教授 (70026110) ; 中尾 充広 九州大学, 大学院・数理学研究院, 教授 (10136418)
• Keywords: 多倍長 / 並列計算 / PVM / シミュレーション / スペクトル法 / 精度保証 / 任意 / 無限

Research Abstract

本年度は研究期間の最終年度にあたる.新たにえられた主な成果を以下に示す.
1.無限精度数値シミュレーションの新たな応用と精度限界の更新.
一様流に対する無限精度数値シミュレーションの適用法を開発した.これはつい最近成功したもので,本格的なシミュレーションは今後の課題である.また,通常では実行不可能な,物理メモリ以上のメモリ容量を必要とする数値シミュレーションを,多倍長計算ライブラリのメモリ節約機能を用いて実行した.これによって,1次元の境界値問題を誤差が10^<-2300>で計算できた.研究開始初年度の2倍の計算精度を達成した.これは,有効数字の桁数でいうと,2300桁というとてつもない超高精度である.通常の数値計算精度である倍精度計算の,実に150倍以上という驚異的な計算精度である.
2.ライブラリの公開.
無限精度数値シミュレーションに必要なライブラリ(連立一次方程式を解くためのガウスの消去法の多倍長計算版とこのPVMによる並列化版)をネット上で公開した.これは,今井のホームページに,ライブラリを収めた科件費報告集を公開することで実現した.
3.可視化における無限拡大法の開発.
無限精度数値シミュレーションが,補間法であるスペクトル法と多倍長計算の組み合わせであることに着目して,データの可視化における無限拡大というまったく新しい手法を開発した.実際,10^<100>倍に拡大しても情報は劣化していなかった.この拡大率は,宇宙の果ての原子核を20cm先のモニター上に数cmの大きさに拡大することに匹敵する,想像を絶するものである.
4.複雑な応用問題ヘの適用のための基礎研究.
流体問題の誤差解析手法の確立や精度保証付き数値計算法の効率化の成功など,多くの関連する研究成果がえられた.

Research Products

• [Publications] H.Imai (共著者に T.Takeuchi): "Infinite Precision Numerical Simulation for PDE systems and Its Applications"数理解析研究所講究録. Vol.1147. 42-50 (2000)
• [Publications] Tarmizi (共著者に H.Imai,T.Takeuchi): "NUMERICAL SIMULATION OF SOME ONE-DIMENSIONAL FREE BOUNDARY PROBLEMS IN ARBITRARY PRECISION"GAKUTO International Series, Mathematical Sciences and Applications. Vol.14. 440-452 (2000)
• [Publications] Tarmizi (共著者に H.Imai,T.Takeuchi): "NUMERICAL SIMULATION OF ONE-DIMENSIONAL FREE BOUNDARY PROBLEMS IN INFINITE PRECISION"Advances in Mathematical Sciences and Applications. Vol.10,No.2(in press). (2000)
• [Publications] H.Ikeda, (共著者に T.Ikeda): "Bifurcation phenomena from standing pulse solutions of bistable reaction-diffusion systems"Journal of Dynamics and Differential Equations. Vol.12,No.1. 117-167 (2000)
• [Publications] T.Ikeda: "Hopf bifurcation of traveling pulses in some bistable reaction-diffusion systems"Methods and Applications of Analysis. Vol.7,No.1. 165-194 (2000)
• [Publications] M.Mimura (共著者に T.Ikeda): "Dynamics of travelling breathers arising in reaction-diffusion systems-ODE modelling approach-"Hiroshima Mathematical Journal. Vol.30,No.2. 221-256 (2000)
• [Publications] M.Tabata: "Error estimates for finite element approximations of drag and lift in nonstationary Navier-Stokes flows"Japan Journal of Industrial and Applied Mathematics. Vol.17. 371-389 (2000)
• [Publications] M.Tabata: "A stabilized finite element method for the Rayleigh-Benard equations with infinite Prandtl number in a spherical shell"Computational Methods in Applied Mechanics and Engineering. Vol.190. 387-402 (2000)
• [Publications] T.Nakaki (共著者に K.Tomoeda): "Numerical waiting time of interfaces in one-dimensional porous medium equations"Free Boundary Problems : Theory and Applications, Mathematical Sciences and Applications. Vol.14. 324-333 (2000)
• [Publications] Minamoto,T. (共著者に Nakao,M.T): "Numerical verification method for solutions of the perturbed Gelfand equation"Methoda and Applications of Analysis. Vol.7. 251-262 (2000)
• [Publications] Toyonaga,K. (共著者 Nakao,M.T): "Numerical enclosure for the optimal threshold probability in discounted Markov decision processes"Bulletin of Informatics and Cybernetics. Vol.32. 81-90 (2000)
• [Publications] T.Nishida (共著者に T.Ikeda): "Pattern Formation of Heat Convection Problems"Lecture Notes in Computational Science and Engineering, Springer, ed.by T.Miyoshi. 1-10 (2001)
• [Publications] Y.Watanabe (共著者に M.Nakao,T.Nishida): "A numerical verification of bifurcated solutions for the heat convection problems"Journal of Mathematical Fluid Mechanics. (to appear). (2001)
• [Publications] 藤間昌一: "A domain-decomposition finite-element scheme for flow problems-Choice of elements-"数理解析研究所講究録. Vol.1129. 9-22 (2000)