保存/散逸型微分方程式に対する離散変分法の拡張に関する研究

Research Project

Project/Area Number	15760043
Research Category	Grant-in-Aid for Young Scientists (B)
Allocation Type	Single-year Grants
Research Field	Engineering fundamentals
Research Institution	The University of Tokyo (2004) Nagoya University (2003)
Principal Investigator	松尾宇泰東京大学, 大学院・情報学環, 講師 (90293670)
Project Period (FY)	2003 – 2004
Project Status	Completed (Fiscal Year 2004)
Budget Amount *help	¥3,500,000 (Direct Cost: ¥3,500,000) Fiscal Year 2004: ¥1,200,000 (Direct Cost: ¥1,200,000) Fiscal Year 2003: ¥2,300,000 (Direct Cost: ¥2,300,000)
Keywords	微分方程式 / 数値解法 / 保存・散逸系 / 保存 / 散逸型微分方程式 / 離散変分法 / composition method / 散逸数値解法 / 一般化BDF解法
Research Abstract	研究実績の概要は以下のとおりである. 本研究は,国産の数値手法である「離散変分法」を拡張することを目的としていたが,その中で特に「時間方向の精度の向上」について重点的な研究を行い,種々の知見を得た. まず,従来,6次までであった次数を,原理的には任意次数に上げる手法を開発した.これは微分方程式の数値解法におけるGeneralized Backward Difference Formula(GBDF)を援用することで達成された.簡単な例題において,実際に高い精度と安定性を持つことが示された. この手法では,離散スキーム全体を一気に,Bounary Value Methodとして解くことが必要となるが,一般にスキームは非線形であるため,それは高次元の非線形代数方程式の数値解法に帰着する.これを解く方法として単純な関数反復や,数値ライブラリによる数値Newton法などを試し,数値Newton法が様々な意味でコストが低く,実用的であることを確認した. さらに,上記の手法が保存/散逸型偏微分方程式の数値解法にも適用可能であることを示した.これは保存/散逸型偏微分方程式を,いったん空間方向のみ適切に離散化して保存/散逸型常微分方程式にし,それに対して上の手法を適用することで達成される.しかしGBDFは非常に計算コストの高い方法であるため,偏微分方程式の場合には本研究で得られた手法が必ずしも理想的ではないことも,各種数値実験により検証された.