• Search Research Projects
  • Search Researchers
  • How to Use
  1. Back to previous page

部分構造反復法による3次元ストークス問題ソルバーの開発と実装

Research Project

Project/Area Number 15740065
Research Category

Grant-in-Aid for Young Scientists (B)

Allocation TypeSingle-year Grants
Research Field General mathematics (including Probability theory/Statistical mathematics)
Research InstitutionKyushu University

Principal Investigator

鈴木 厚  九州大学, 大学院・数理学研究院, 助手 (60284155)

Project Period (FY) 2003 – 2004
Project Status Completed (Fiscal Year 2004)
Budget Amount *help
¥2,000,000 (Direct Cost: ¥2,000,000)
Fiscal Year 2004: ¥800,000 (Direct Cost: ¥800,000)
Fiscal Year 2003: ¥1,200,000 (Direct Cost: ¥1,200,000)
Keywords部分構造反復法 / バランシング前処理 / 正射影付き有限要素方程式 / 安定化有限要素法 / 共役勾配法 / 合同な部分領域
Research Abstract

ストークス問題では,滑り境界条件を課す場合があり,解は剛体回転の自由度と圧力の不定性自由度を含む.この自由度を取り除くために線形拘束を用いるが,これは関数空間の制約として取り扱う.全節点自由度の空間からの制限を実現する正射影によりアルゴリズムを構成することができる.有限要素剛性行列は全節点に自由度を持つ有限要素基底を用いて生成することができ,手続きが簡単である.この行列自身の正則性は考える必要がなく,線形拘束をみたす制約付き空間の上での正則性,すなわち正射影を行列の両辺に作用させたものの正則性を議論することで十分であり,これは有限要素解の一意性から自然に導かれる.得られた連立一次方程式は正射影付きクリロフ部分空間法(特に共役勾配法)を用いて解くことができる.行列・ベクトル積とベクトル加減算に正射影を付加するプログラム実装を行うことで,数値誤差の蓄積を防ぎ,制約付き空間のなかでクリロフ部分空間法を実行することができる.
部分構造反復法では部分領域内部の自由度を消去し,内部境界上の自由度の方程式に変換する.このためには部分領域での係数行列の成分を知る必要があるが,正射影を伴う有限要素方程式では,成分の値を陽に得ることは難しい.正射影作用素を領域分割に対応して分割することで,有限要素方程式に対する汎用的な部分構造反復法を構成した.各部分領域で,領域毎に線形拘束を課すことで,行列成分を陽に必要とする問題を回避できるが,過剰な制約となってしまうため,全体領域で拘束条件の自由度に対応する問題を付加する.これにより内部境界上での問題の正則性が保たれることがわかった.
この手法によるストークス方程式の部分構造反復法に,弾性体問題に対して優れた収束特性をもつバランシング前処理(BDD法)を導入した.BDD法はPavarin-Widlund[2002]によって圧力に不連続要素を用いるQ2/QO要素の有限要素近似に拡張されているが,圧力に連続要素を用いるP1/P1安定化近似に対してアルゴリズムを構築した.

Report

(2 results)
  • 2004 Annual Research Report
  • 2003 Annual Research Report
  • Research Products

    (5 results)

All 2005 2004 Other

All Journal Article (2 results) Publications (3 results)

  • [Journal Article] Finite element matrices in congruent subdomains and their effective use for large scale computations2005

    • Author(s)
      A.Suzuki, M.Tabata
    • Journal Title

      International Journal for Numerical Methods in Enginnering 62(印刷中)

    • NAID

      120000981469

    • Related Report
      2004 Annual Research Report
  • [Journal Article] Finite element computation of Earth's mantle convection problem on a vector-parallel computer2004

    • Author(s)
      A.Suzuki
    • Journal Title

      Computational Mechanics : Proceedings of Sixth World Congress on Computational Mechanics in conjunction with the Second Asian-Pacific Congress on Computational Mechanics

      Pages: 645-652

    • Related Report
      2004 Annual Research Report
  • [Publications] 鈴木 厚: "地球マントル対流問題のための有限要素コードのベクトル化と並列化"計算工学講演会論文集. 8・2. 729-730 (2003)

    • Related Report
      2003 Annual Research Report
  • [Publications] 鈴木 厚: "有限要素剛性行列に対する汎用的部分構造反復法"第53回理論応用力学講演会 講演論文集. 403-404 (2004)

    • Related Report
      2003 Annual Research Report
  • [Publications] A.Suzuki, M.Tabata: "Finite element matrices in congruent subdomains and their effective use for large-scale computations."MHF preprint series, Kyushu University. 2003-6. 16 (2003)

    • Related Report
      2003 Annual Research Report

URL: 

Published: 2003-04-01   Modified: 2016-04-21  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi