Fluid Dynamics
Finite Element Formulation for Fluid-11

■プログラムNSEQで使われている変数名■ 次に、プログラムNSEQで使われている変数名について説明します。NSEQでは、既にこれまでに紹介したプログラムで使われている変数名が多数、再使用されています。目的や意味は、いままでと同じですので、ここでの説明は省きます。

まず最初に、入力(subroutine INPUT)で使われる変数名を説明します。

変数名目的/意味データが格納されている
ファイル名
PROJECT計算のプロジェクト名(英数の8文字以内)NSDATA4.FLN
NF入力ファイル数NSDATA4.FLN
FNAME(I)各々入力ファイル名NSDATA4.FLN
VISCO流体の粘性係数プロジェクト名.DAT
FLMDA流体の第2粘性係数プロジェクト名.DAT
NE要素数プロジェクト名.DAT
NNODE節点数プロジェクト名.DAT
NBFx&y-方向のDirichlet境界条件の数プロジェクト名.DAT
DTタイムステッププロジェクト名.DAT
DTMAXタイムステップの最大値プロジェクト名.DAT
TMAX経過時間の最大値プロジェクト名.DAT
ITEFIX各タイムステップでの繰り返し計算の回数プロジェクト名.DAT
MAXITE繰り返し計算の最大回数プロジェクト名.DAT
ERMAX最大誤差(下を参照)プロジェクト名.DAT
NODEX(IEL,J)要素の節点情報プロジェクト名.DAT
XCOORD(1,NODE)x-座標値プロジェクト名.DAT
XCOORD(2,NODE)y-座標値プロジェクト名.DAT
IBNDF(I)x&y-方向のDirichlet境界条件の節点番号プロジェクト名.DAT
BV(I)x&y-方向のDirichlet境界条件の境界値プロジェクト名.DAT
TIME経過時間プロジェクト名.BIN
DTタイムステッププロジェクト名.BIN
U(1,I)x-方向の流速の初期値プロジェクト名.BIN
U(2,I)y-方向の流速の初期値プロジェクト名.BIN

入力ファイルには、基本的に3種類あります。ファイルは、ASCIIコードで書かれていますが、プロジェクト名.BINは、バイナリーファイルですから注意して下さい。
各タイムステップでは、Convective 項からの誤差を小さくするために、subroutine SOLVEFで繰り返し計算を行っています。その誤差判定にERMAXが使われています。ERMAXと比較される計算誤差は、次式で計算されています。

計算誤差=( U(t+Δt)-U(t)またはV(t+Δt)-V(t)の最大値 ) / (U(t)またはV(t)の最大値 )

BACK NEXT
Menu Govern FEM 2nd Visco DCF Cylinder /Heat Lub