流体問題における粘性係数依存性を克服する有限要素スキームとその高速求解法の確立

克服流体问题中粘度系数依赖性的有限元方案建立及其高速求解方法

• 批准号: 18K13461
• 负责人: 内海 晋弥
• 金额: $ 1.58万
• 依托单位: Gakushuin University
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Early-Career Scientists
• 财政年份: 2018
• 资助国家: 日本
• 起止时间: 2018-04-01 至 2024-03-31
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/en/grant/KAKENHI-PROJECT-18K13461/
• 关键词: ラグランジュ・ガレルキン法; ナヴィエ・ストークス方程式; 数値積分; 移流拡散方程式; L^2 収束次数; 射影法; 圧力安定化法; 曲線に囲まれた領域; 粘性係数依存性; 圧力の誤差評価; 数値解析; 変形エルミート要素; 混合有限要素近似; ストークス方程式; クリロフ部分空間法; 滑らかな領域; ラプラシアンの固有値の包含; オセーン方程式; 有限要素法; 誤差評価; 前処理

(i) 非圧縮粘性流の運動を記述するナヴィエ・ストークス方程式に対する数値解法の一つである，ラグランジュ・ガレルキンを考察している．本手法に特徴的な合成関数を含む積分は厳密に行うことが困難であり，そこに数値積分を用いる近似計算が使われた場合の理論解析は知られていなかった．田端正久氏との共同研究により，時間刻みが小さいという条件下で，三角形（四面体）の内部のみに積分点を持つ数値積分を用いるラグランジュ・ガレルキン法の解の収束性が証明された．また，辺上に積分点を持つ公式を使った時の収束性は証明されていないが，安定に計算できることが実験的に分かった．これらの結果は専門誌への掲載が決定した．(ii) 前年度までに結果を得ていた，オセーン問題のためのラグランジュ・ガレルキン，射影，圧力安定化の3種結合スキームに関して，圧力安定化項のパラメータの取り方を考察した．誤差評価の証明の過程でパラメータの指針はこれまでに分かっていたが，物理量の次元に着目して再考した．移流項の大きさは，パラメータの選択に影響を及ぼすことが分かり，ストークス問題に対してはパラメータ選択を再考する必要があることが分かった．(iii) 1次元移流拡散問題に対するラグランジュ・ガレルキン法を考え，L^2 ノルムでの収束次数を考察した．1次要素を用い，時間刻みが十分小さい場合，拡散係数に依存しない h^{3/2} の L^2 収束次数が得られた．ここで，h はメッシュ代表長を表す．これは，これまで得られていた h^1 の収束次数から改善されている．ある程度一般的な仮定の下で，h^{3/2} は他の手法も含めて現在知られている，拡散係数の小ささに依存しない最良次数の誤差評価である．

(I) not 圧 shrinkage of の viscous flow を account す る ナ ヴ ィ エ · ス ト ー ク ス equation に す seaborne る the numerical solution to a つ の で あ る, ラ グ ラ ン ジ ュ · ガ レ ル キ ン を investigation し て い る. This technique に な synthesis of 徴 masato number を は む integral line 厳 dense に う こ と が difficult で あ り, そ こ に を using the numerical integral い る approximate calculation が make わ れ た occasions の theory analytical は know ら れ て い な か っ た. Tian's regular long と の joint research に よ り, time carved み が small さ い と い う で, under the condition of the triangle (tetrahedron) の internal の み に integral point を hold つ を using the numerical integral い る ラ グ ラ ン ジ ュ · ガ レ ル キ ン method の solution の 収 beam sex が prove さ れ た. ま た, 辺 を に integral point on a つ formula を っ た when の 収 beam sex は prove さ れ て い な い が, settle に computing で き る こ と が be 験 に points か っ た. The れら れら results are announced in the へ journal of the へ that が decisions are made た た. (ii) before the annual ま で に results て を い た, オ セ ー ン problem の た め の ラ グ ラ ン ジ ュ · ガ レ ル キ ン, projection, pressure stabilization の three combination ス キ ー ム に masato し て, pressure stabilization item の パ ラ メ ー タ の take り party を investigation し た. Error evaluation 価 の prove の process で パ ラ メ ー タ の pointer は こ れ ま で に points か っ て い た が, quantities の dimensional に with mesh し て take exam し た. Advection term の big き さ は, パ ラ メ ー タ の sentaku に influence を and ぼ す こ と が points か り, ス ト ー ク ス problem に し seaborne て は パ ラ メ ー タ sentaku を take exam す る necessary が あ る こ と が points か っ た. (iii) 1 dimensional advection company, scattered issues に す seaborne る ラ グ ラ ン ジ ュ · ガ レ ル キ ン method を え, L ^ 2 ノ ル ム で の 収 bundle number を investigation し た. One element を い, time carved み が very small さ い occasions, company, dispersion coefficient に dependent し な い h ^ {L} の L ^ 2 収 bundle number が must ら れ た. Youdaoplaceholder0 で, h メッシュ メッシュ represents the long を table す. <s:1> れ, る れまで the られて た た h^1 <s:1> bunch count た ら improves されて る る る. あ る degree generally な 仮 under fixed の で, h ^ {l} は he の technique contains も め て now know ら れ て い る, small company, dispersion coefficient の さ さ に dependent し な い most good times の error evaluation 価 で あ る.

项目成果

曲有限要素空間の構成について

关于弯曲有限元空间的组成

• 发表时间: 2019
• 作者: Fukuizumi Reika;Hoshino Masato;Inui Takahisa;内海 晋弥
• 通讯作者: 内海 晋弥

内海 晋弥 の web ページ

内海慎也的网页

曲線で囲まれた領域におけるラプラシアンの固有値の包含について

关于在曲线所包围的区域中包含拉普拉斯算子的特征值

• 发表时间: 2020
• 作者: Koya Sakakibara;内海 晋弥
• 通讯作者: 内海 晋弥

Guaranteed bounds for the eigenvalues of Laplacian in curved domains

弯曲域中拉普拉斯算子特征值的保证界限

• 发表时间: 2019
• 作者: 榊原 航也;加藤孝盛;S. Uchiumi;S. Uchiumi
• 通讯作者: S. Uchiumi

数値積分を用いるLagrange-Galerkinスキームの収束性

使用数值积分的拉格朗日-伽辽金格式的收敛性

• 发表时间: 2018
• 作者: Gustafson Stephen;Inui Takahisa;加藤孝盛;榊原 航也;内海 晋弥，田端 正久
• 通讯作者: 内海 晋弥，田端 正久