双曲型偏微分方程式の有限要素解析における衝撃捕獲法の開発

1995 Fiscal Year Annual Research Report

• Project/Area Number: 07855010
• Research Institution: Tokyo Institute of Technology
• Principal Investigator: 山田 貴博 東京工業大学, 工業材料研究所, 助手 (40240022)
• Keywords: 双曲型偏微分方程式 / 数値解析 / 有限要素法 / 気泡関数要素 / 安定化法 / 衝撃捕獲 / 人工粘性 / 圧縮性オイラー方程式

Research Abstract

本研究は,圧縮性オイラー方程式などの双曲型偏微分方程式の有限要素解析において、衝撃波を、鋭敏かつ安定に捕えることのできる衝撃捕獲手法の開発を行うものである.本年度は,以下のように研究を行った.
1.衝撃波を含まない双曲型偏微分方程式に対する気泡関数要素の有効性について検討を行い、移流方向の情報すなわち特性曲線の方向を使うことなく、移流の大きさのみの情報で適切な上流方向の数値粘性を導入することができることが確認された。また、このような性質は特に連立の双曲型方程式である圧縮性オイラー方程式に対して有効であることが分かった。
2.基本的な1変数の双曲型偏微分方程式の気泡関数要素を用いた有限要素解析において、衝撃波の不連続性を要素毎に評価する測度とその測度を用いた人工粘性について数値実験により検討を行った。その結果、定常衝撃波の問題には、物理量の2階微分と要素代表長さの3乗に比例する人工粘性項が有効であることが分かった。一方、非定常衝撃波に対しては物理量の1階微分と要素代表長さの2乗に比例する人工粘性により、移動する衝撃波を鋭敏に捕獲することができることが示された。
3.以上のような成果をもとに、大域的な方程式の双曲性に対しては気泡関数による上流化を行い、衝撃波に対しては局所的な人工粘性を用いる圧縮性オイラー方程式に対する有限要素解析手法を開発し、数値実験により手法の妥当性について検討を行った。その結果、本手法においては、衝撃波が3〜5要素程度の幅で捕獲され、また、要素分割による数値解の依存性も少ないことが明かとなり、解析手法の有効性が確認された。

Research Products

• [Publications] 山田貴博: "非圧縮体に対するMINI要素の数値不安定性" 構造工学における数値解析法シンポジウム論文集. 19. 101-106 (1995)
• [Publications] 山田貴博: "圧縮性オイラー方程式に体する気泡関数要素" 構造工学における数値解析法シンポジウム論文集. 19. 445-450 (1995)
• [Publications] 山田貴博: "Stokes方程式に対する安定化有限要素法とMINI要素" 日本応用数理学会平成7年度年会研究発表予稿集. 96-97 (1995)
• [Publications] Takahiro Yamada: "A bubble element for inviscid flows" Proc. of the 9th International Conference on Finite Elements in Fluids. part II. 1567-1576 (1995)
• [Publications] 山田貴博: "圧縮性オイラー方程式に対する気泡関数要素を用いた有限要素スキーム" 第9回数値流体力学シンポジウム講演論文集. 139-140 (1995)