注意書き

READ ME FIRST
1. 以下の項目を予め了解された方にソフトウェア利用を許可いたします。
・ダウンロード・インストール，ならびにお手持ちのコンピュータの利用環境下での
ソフトウェアの使用等の運用 (以下，運用等) については，読者の責任と判断によ
っておこなってください。
・読者の運用等の結果に際して，編者，著者及び丸善出版株式会社はいかなる責任も
負いません。
・また，運用等に関するいかなる質問に対しても，編者，著者及び丸善出版株式会社
はお答えできません。
・本プログラム提供サービスは予告なしに終了する可能性があります。
・本プログラム内容は予告なく修正を施す場合があります。
--------------------------------------------------------------------------------------------------------------------2. このフォルダ内のプログラムは過膨張及び不足膨張ジェットを計算します。また，フォ
ルダ内には下記のファイルが格納されています。
ReadMeFirst.pdf
本注意書き
laxf2v4.f
計算プログラム
input-pr03.dat
入力データ (過膨張，p-jet/pa=0.3，Mexit=3.0 の場合)
input-pr10.dat
入力データ (不足膨張，p-jet/pa=10.0，Mexit=1.0 の場合)
tpp4.f
簡易可視化プログラム
3. フォートランプログラム laxf2v4.f をコンパイルして実行モジュールを作成してくださ
い。例えば，インテルのフォートランを利用する場合は
> ifort laxf2v4.f
と入力して laxf2v4.exe という実行モジュールを作成いたします。インテルのフォートラン
の代わりにフリーのフォートランコンパイラ(g95 等)も利用可能です。
4. 例えば WINDOWS の場合に過膨張ジェットの計算を実行するには DOS 窓から
> laxf2v4 < input-pr03.dat
と入力してください。ダウンロードしたままの input-pr03.dat を用いると 200,000 ステップ計
算行った後に終了します。計算途中では 100 ステップ毎にステップ数，無次元時刻，時間刻
み幅，残差(L1 ノルム)，最大残差(L1 ノルム)を与える(i,j)座標，最大ノルム(L2)を出力しま
す。
5. 計算を実行すると
d15-pr03.dat
(フィールドデータ，約 6MB)
tec-pr03.plt (TECPLOT 形式の出力データ，約 6MB)
がディスクに出力されます。いずれもアスキー形式のデータで編集可能です。もしも
TECPLOT の利用が可能な場合は tec-pr03.plt を読み込んで流れ場を可視化して下さい。
TECPLOT 形式のデータが読める可視化ソフトも同様に利用可能です。フリーソフトでは
VisIt の利用実績があります(現時点で利用可能であることを保証するものではありませ
ん)。”VisIt, visualization” で検索してみてください。米国ローレンスリバモア研究所のサイ
トからダウンロードできることが分かります。
6. tec-pr03.plt は各格子点における密度(𝜌)，圧力(𝑝)，エントロピー(𝑠)，マッハ数(𝑀)を求めて
います。プロットする変数を選んでください。
7. 可視化ソフトを利用することが困難な方は，tpp4.f を利用してください。これは文字セ
ット(0,1,2,3,4,5,6,7,8,9,*)を用いた簡易可視化ソフトです。利用時は最初にフォートランコン
パイラーで tpp4.exe という実行モジュールを作成してください。次に
> tpp4.exe
と入力すると
> input: file name of field data
と出力されますので，ここで
> d15-pr03.dat
と入力します。すると
> select: ro=1,p=2,s=3,xm=4
と出力されますので，適当な変数の番号を入力してください。例えば 4 (マッハ数)を選んだ
場合，
>4
と入力してください。すると
emin= 7.479280979280894E-003
input: tmin and tmax
emax=
3.00000000000000
と出力されますので，例えば最小値 tmin=0.01, 最大値 tmax=2.99 を
> 0.01, 2.99
と入力します。すると次ページの結果を得ます。0.01 より小さな値を持つ計算セルには 0
が，2.99 より大きな値を持つ計算セルには*が記録され，その間を 1∼9 の数字で分割します。
かつて可視化ソフトが使えなかった時代には，このようにラインプリンターに文字を打ち出
すことによって流れ場の様子を確認することが行われました。
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111001111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111121
11111111111111111111111111111111111111111111111111111111111111111111121
11111111111111111111111111111111111111111111111111111111111111111111121
11111111111111111111111111111111111111111111111111111111111111111111122
42211111111111111111111111111111111111111111111111111111111111111111122
98643221111111111111111111111111222333333322222111111111111111111111222
**976654321111111111111111112233444555555544443322211111111111112223322
***98776654322111111111112234455666666666666655544332221111112223444543
****9977777655432211112234556667777777777777776666554432222223345566663
******98777776665432234456667777788888888888777777666554333344556667774
*******9877777776654455666777778888888888888888777776665444555666777774
*********98777777765556667777888888899999998888877777765555566667777884
**********9977776555555666777888889999999999988887776665555556667778884
************98655555555566777888899999999999999887766655555556667778884
**********9977776555555666777888889999999999988887776665555556667778884
*********98777777765556667777888888899999998888877777765555566667777884
*******9877777776654455666777778888888888888888777776665444555666777774
******98777776665432234456667777788888888888777777666554333344556667774
****9977777655432211112234556667777777777777776666554432222223345566663
***98776654322111111111112234455666666666666655544332221111112223444543
**976654321111111111111111112233444555555544443322211111111111112223322
98643221111111111111111111111111222333333322222111111111111111111111222
42211111111111111111111111111111111111111111111111111111111111111111122
11111111111111111111111111111111111111111111111111111111111111111111122
11111111111111111111111111111111111111111111111111111111111111111111121
11111111111111111111111111111111111111111111111111111111111111111111121
11111111111111111111111111111111111111111111111111111111111111111111121
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111001111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
過膨張ジェット流れ(マッハ数)の可視化例。第 9 章の 9.8.2 図と見比べてください。不足膨
張ジェットの場合は次ページのようになります。
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111100111
11111111111111111111111111111111111111111111111111111111111111000000111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111112222222222222222222211111111111111111111111
11111111111111111111111122233334444444444443333322222111111122222222222
11111111111111111111222334444455555555555555444444333222222233333444432
11111111111111111223344455555556666666666655555555544443233444455555553
11111111111111123344555566666666666666666666666665555544444555555555554
11111111111112344555666666666666666666666666666666665444444555566666664
01111111111234455666666666677777777776666666666666554445555555555666654
11111111123445566666677777777777777777777777666654444444444445555555553
11111111234556666777777788888888999998888777775433332222333333444444432
11111113455666677777888888889999999999999999852211111111111222222332222
11111134566667777888888888889999999999999999*41111111111111111111221111
11112345666777777888888888888999999999999999*51111111111111111111111111
11124556677777777788888888888899999999999999971111111111111111111111111
11245666677777777778888888888899999999999999981111111111111111111111111
12456666666777777777888888888889999999999999982111111111111111111111111
34555666666677777777888888888889999999999999992111111111111111111111111
34455555666667777777788888888889999999999999992111111111111111111111111
33444555666667777777788888888888999999999999992111111111111111111111111
23344555566666777777788888888888999999999999992111111111111111111111111
33444555666667777777788888888888999999999999992111111111111111111111111
34455555666667777777788888888889999999999999992111111111111111111111111
34555666666677777777888888888889999999999999992111111111111111111111111
12456666666777777777888888888889999999999999982111111111111111111111111
11245666677777777778888888888899999999999999981111111111111111111111111
11124556677777777788888888888899999999999999971111111111111111111111111
11112345666777777888888888888999999999999999*51111111111111111111111111
11111134566667777888888888889999999999999999*41111111111111111111221111
11111113455666677777888888889999999999999999852211111111111222222332222
11111111234556666777777788888888999998888777775433332222333333444444432
11111111123445566666677777777777777777777777666654444444444445555555553
01111111111234455666666666677777777776666666666666554445555555555666654
11111111111112344555666666666666666666666666666666665444444555566666664
11111111111111123344555566666666666666666666666665555544444555555555554
11111111111111111223344455555556666666666655555555544443233444455555553
11111111111111111111222334444455555555555555444444333222222233333444432
11111111111111111111111122233334444444444443333322222111111122222222222
11111111111111111111111111112222222222222222222211111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111111111111
11111111111111111111111111111111111111111111111111111111111111000000111
11111111111111111111111111111111111111111111111111111111111111111100111
11111111111111111111111111111111111111111111111111111111111111111111111
不足膨張ジェット流れ(マッハ数)の可視化例。第 9 章の 9.8.4 図と見比べてください。

Download Report