权益分类	功能权益	普通用户	{{item.name}}会员
{{category.name}}	{{benefitItem.name}}

hp-有限要素法による3次元計算コードの開発

使用 hp 有限元法开发 3D 计算代码

基本信息

批准号：
21540135
负责人：
鈴木厚
金额：
$ 1.66万
依托单位：
Kyushu University
依托单位国家：
日本
项目类别：
Grant-in-Aid for Scientific Research (C)
财政年份：
2009
资助国家：
日本
起止时间：
2009 至 2011
项目状态：
已结题

来源：
https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-21540135/
关键词：
数値流体力学有限要素法 hp-解適合近似並列計算

项目摘要

非圧縮流れの基礎方程式であるナビエ・ストークス方程式で記述される三次元流れ問題において,なるべく少ない計算量で高解像度の計算を実現でぎるアルゴリズムを開発することが本研究の目的である.ナビエ・ストークス方程式の物質微分項に特性曲線による近似を適用して時間微分の近似を行うと各ステップで線形のストークス方程式を解くことになる.本年度はこの定常のストークス方程式のhp-有限要素法アルゴリズムの開発を開始した.時間ステップを止めて,外力を既知とし,それに対応する流れ場を高精度に近似する適応的な要素細分(h)と高次要素(p)を導入する.計算効率を考え,六面体ではなく,四面体を細分する有限要素列を採用した.局所要素細分のアルゴリズムは,ある部分では要素細分が必要であり,またある部分では粗い要素分割に戻す操作が必要になる.通常の有限要素法では,ある要素を細分した際,その要素に隣接する要素も対応して細分する節点に関する適合性条件を満たす必要がある.このため,要素細分の局所性が弱くなる.これに対して,粗い要素の辺あるいは面の中点に細かい要素の節点を配置できる,不連続有限要素法(Discontinuous Galerkin (DG)法)がある.予備的なコードを作成したところ,DG法では離散化後の剛性行列が非対称なものが近似能力がすぐれているが,大規模非対称行列を解く困難さか残ることがわかった.このため,通常の要素細分アルゴリズムに基づき,また安定化手法を導入してストークス方程式を離散化する手法を採用した.高次要素においては,四面体の頂点だけでなく内部に節点を配置するが,これは細分を行う1次要素と同等の自由度の増加になるため,要素細分の手続きを共通化できることが分かった.

The fundamental equations of non-compression flow describe the problems of three-dimensional flow, and the calculation of high resolution is realized. This is the purpose of this study. The approximation of the characteristic curve of the material differential term of the equation is applicable to the approximation of the time differential. This year, the development of hp-finite element method for the steady state equations has begun. The time is short, the external force is known, and the flow field is approximated with high precision. The element is subdivided (h) and the higher order element (p) is introduced. The calculation rate of hexahedron is higher than tetrahedron. The segmentation of office elements is necessary, and the segmentation of elements is necessary, and the segmentation of elements is necessary. In general, when the finite element method is used to subdivide the elements, the elements are adjacent to each other, and the node is related to the suitability condition. The elements are subdivided and the nature of the bureau is weak. For this reason, the coarse elements are arranged at the midpoint of the plane, and the nodes of the fine elements are arranged at the midpoint of the plane, without continuous Galerkin (DG) method. The method of DG is difficult to solve the problem of large scale asymmetric arrays after discretization. This is a common element subdivision method, which is based on the introduction of a stabilization method and the discretization of an equation. Higher order elements, vertices of tetrahedra, internal nodes of tetrahedra, subdivision of first order elements, increase of degrees of freedom, subdivision of elements, commonalization of elements.