solve double mach reflection problem x方向格子点数 y方向格子点数 時間刻み幅 比熱比と関連定数 定数(計算速度性能に寄与するかはよくわからない) パラメータ等設定 includeを用いるのは本来好ましくない 最大反復数 終了時間
格子点座標 0とmax+1は便宜上 計算領域は0-max
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | x(0:jmax+1,0:kmax+1) | ||||
real | :: | y(0:jmax+1,0:kmax+1) |
逆数を定義しておく intgがちょっと速くなるはず
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | dltx | ||||
real | :: | dlty |
保存量
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | cons(4,0:jmax+1,0:kmax+1) |
基本量
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | pris(4,0:jmax+1,0:kmax+1) |
t:現在時間, tt:サブステップの時間
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | tn | ||||
real | :: | tt |
反復数
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
integer | :: | im |
before shock state
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | r0 | ||||
real | :: | p0 | ||||
real | :: | u0 | ||||
real | :: | v0 | ||||
real | :: | ru0 | ||||
real | :: | rv0 | ||||
real | :: | e0 |
after shock state
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | r1 | ||||
real | :: | p1 | ||||
real | :: | u1 | ||||
real | :: | v1 | ||||
real | :: | ru1 | ||||
real | :: | rv1 | ||||
real | :: | e1 |
非粘性流束
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | xflx(4,jmax,kmax) | ||||
real | :: | yflx(4,jmax,kmax) |
流束評価面に補間された値 l:left, r:right
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | pxl(4,0:jmax+1,0:kmax+1) | ||||
real | :: | pxr(4,0:jmax+1,0:kmax+1) | ||||
real | :: | pyl(4,0:jmax+1,0:kmax+1) | ||||
real | :: | pyr(4,0:jmax+1,0:kmax+1) |
時間増分
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | df(4,jmax,kmax) |
非粘性流束評価で用いる値
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | conL(4) | ||||
real | :: | conR(4) | ||||
real | :: | dnflx(4) |
流束評価面の法線ベクトル
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real | :: | dnx | ||||
real | :: | dny |
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
double precision | :: | conn(4,0:jmax+1,0:kmax+1) | n stepの保存量の一時的置き場 |