iSALE - 3. 例. トレーサー粒子追跡 Diff
- Added parts are displayed like this.
- Deleted parts are displayed
like this.
== 概要
iSALEでは各計算グリッドについてトレーサー粒子を設定して,その粒子の移動,温度圧力履歴等の物理量を追跡することが出来る.
=== トレーサー粒子の追加
asteroid.inpにトレーサー粒子の計算を追加する.
* 例
---------------------Tracer Particle Parameter -------------------------
TR_SAVE Option for saving : 1 <- トレーサー粒子を計算するかどうか
TR_QUAL integration qual. : 1
TR_SPCH tracer spacing X : -1.0 : -1.0 <- x軸方向の粒子間隔
TR_SPCV tracer spacing Y : -1.0 : -1.0 <- y軸方向の粒子間隔
TR_VAR add. tracer fiels : #TrP-TrT-Trp-Trt-TrV-TrA# <- 計算する物理量
------------------------------------------------------------------------
* TR_SAVE:トレーサー粒子を計算するかどうかの選択
* 0 -> しない
* 1 -> する
* TR_SPCH, TRSPCV:トレーサー粒子の空間分布
* 正の値 -> 絶対値での指定.単位 [m]
* 負の値 -> 計算グリッド間隔を(1として)規格化した値での指定.
* TR_VAR:記録するトレーサー粒子の物理量
* TrT:最高衝撃温度
* TrP:最大衝撃圧力
* Trt:衝撃温度履歴
* Trp:衝撃圧力履歴
* TrV:体積ひずみ
* TrA:体積膨張
=== トレーサー粒子の描画
==== iSALEPlot上での描画
* コマンド
--------- TRACER PARAMETERS -------------------------------------------
TR_TYPE Tracer lines, grid or points (1,2,3) : -1 <- トレーサーの描画タイプ
TR_SPACE Spacing between tracers to be plotted : 10 <- トレーサーの描画間隔
TR_SIZE Size of tracer plot marker : 1 <- トレーサーのサイズ
-----------------------------------------------------------------------
* TR_TYPE:トレーサー粒子の描画タイプ
* 1 -> 「線」
* 2 -> 「グリッド」
* 3 -> 「点」
※現段階では,1-3のどれを入力しても「点」として出力される.原因は不明.
* TR_SPACE:トレーサーの描画間隔
* 0 -> 描画されない
*それ以外の整数 -> 現段階ではまだ反映されず,「1」として出力される.原因不明.
* TR_SIZE:トレーサーの大きさ
* 0 -> 描画されない
* それ以外 -> 値に従って大きくなる.
* 図
いくつかの描画例を下に示す.(左:圧力[Pa],右:温度[K])
{{attach_view(00000.png)}}
{{attach_view(00065.png)}}
{{attach_view(00112.png)}}
{{attach_view(00650.png)}}
{{attach_view(01950.png)}}
{{attach_view(02600.png)}}
{{attach_view(03250.png)}}
==== VIMoD上での描画
* パネルの説明
* 図
=== トレーサー粒子の解析
====トレーサー粒子ごとのデータの出力
iSALEPLotの描画コマンドにオプションとして「-T <描画する粒子の番号>」を追加することで.txtファイルを出力可能.
例:1番目のトレーサー粒子を出力.
../iSALEPlot -f ../iSALEPlot.inp -m jdata.dat -T 1
上記の操作により,「tracer-1000001.txt」というファイルが生成される.
中身は次の通り
Time xmark ymark TmpPre <- 左から時間,粒子のx座標,y座標,温度,圧力
0.0000000 50.000000 50.000000 753.21973 1.19651878E+09
1.5413910 206.66608 -2828.1284 753.21973 1.19651878E+09
3.0780520 258.27499 -3203.5825 753.21973 1.19651866E+09
4.6188130 261.67145 -3091.5713 753.21948 1.19648666E+09
6.1562600 270.00211 -2994.5752 753.22052 1.19664128E+09
7.6939030 276.39139 -2916.2378 753.23224 1.19834803E+09
9.2336054 279.68808 -2859.0605 754.13019 1.31637120E+09
10.770640 281.79733 -2828.4565 753.06561 1.18852262E+09
12.312140 283.42297 -2823.1621 753.07471 1.19253312E+09
13.848370 284.24585 -2840.7695 753.13367 1.20153984E+09
15.386140 271.85925 -2867.3828 753.03326 1.18787840E+09
16.924549 247.52042 -2886.7148 753.03821 1.18865293E+09
18.464939 246.17911 -2892.1021 752.95508 1.17706509E+09
20.003201 246.77846 -2901.5337 752.46234 1.10552653E+09
21.542770 245.76292 -2905.5068 753.19196 1.21384166E+09
23.077120 245.13913 -2909.4204 752.92352 1.17844365E+09
24.616541 245.11148 -2915.7446 752.61328 1.13912909E+09
26.154341 245.04880 -2924.6899 752.89874 1.18152627E+09
27.692949 245.04442 -2933.1025 752.97357 1.19366208E+09
29.233200 245.11763 -2937.0903 752.88684 1.18303552E+09
30.770411 245.26674 -2939.8745 752.48529 1.12413773E+09
32.310612 245.04799 -2940.2998 752.45294 1.12343501E+09
33.846161 244.43181 -2928.2373 752.54181 1.13237850E+09
35.388008 243.84512 -2930.9766 753.10828 1.22200755E+09
36.925468 243.14340 -2952.9873 753.79694 1.30881946E+09
38.464211 243.19815 -2970.5381 752.60376 1.15911859E+09
40.002048 243.25449 -2949.4482 752.75793 1.18333760E+09
41.543320 243.10425 -2934.8359 752.35272 1.12665702E+09
43.079689 243.10820 -2922.2207 752.18793 1.10040986E+09
44.615910 242.94418 -2914.3911 752.86719 1.20830656E+09
46.157761 242.23381 -2921.2798 752.64325 1.17738304E+09
47.693890 242.00459 -2923.3892 752.52234 1.15863322E+09
49.234959 241.93341 -2924.5425 752.70111 1.19225318E+09
50.769939 241.51730 -2938.4551 752.02716 1.08899290E+09
52.308472 241.61102 -2941.9395 751.87012 1.07313318E+09
53.849571 242.01355 -2936.9990 752.31018 1.14422451E+09
55.389370 240.86028 -2936.4121 752.29480 1.14128781E+09
56.926270 241.20326 -2933.1685 752.29950 1.14750758E+09
58.463509 241.56766 -2932.0737 752.55896 1.18605734E+09
60.003208 241.44037 -2935.7129 752.70404 1.20741389E+09
61.541512 242.22253 -2935.6523 752.05847 1.11615578E+09
* 左の3コラム:デフォルトで出力される値.
* 右の2コラム:「iSALEPlot.inp」で指定した描画パラメーター(詳細はこちらを参照)
==== 移動パス(流線)
深さ1 km付近に位置するトレーサー粒子の移動パス.
*下図は出力された各トレーサー粒子をコンパイルして描画したものである.
{{attach_view(trace_x.png)}}{{attach_view(trace.png)}}
==== 温度・圧力履歴
衝突地点下での複数粒子の経験する温度圧力(赤点)のコンパイル
* 線:各地点での温度圧力履歴
* 衝突地点直下(青線)
* 衝突地点下850 m 〜 2950 m(黒〜灰色線)
{{attach_view(pt_path_1.png)}}
iSALEでは各計算グリッドについてトレーサー粒子を設定して,その粒子の移動,温度圧力履歴等の物理量を追跡することが出来る.
=== トレーサー粒子の追加
asteroid.inpにトレーサー粒子の計算を追加する.
* 例
---------------------Tracer Particle Parameter -------------------------
TR_SAVE Option for saving : 1 <- トレーサー粒子を計算するかどうか
TR_QUAL integration qual. : 1
TR_SPCH tracer spacing X : -1.0 : -1.0 <- x軸方向の粒子間隔
TR_SPCV tracer spacing Y : -1.0 : -1.0 <- y軸方向の粒子間隔
TR_VAR add. tracer fiels : #TrP-TrT-Trp-Trt-TrV-TrA# <- 計算する物理量
------------------------------------------------------------------------
* TR_SAVE:トレーサー粒子を計算するかどうかの選択
* 0 -> しない
* 1 -> する
* TR_SPCH, TRSPCV:トレーサー粒子の空間分布
* 正の値 -> 絶対値での指定.単位 [m]
* 負の値 -> 計算グリッド間隔を(1として)規格化した値での指定.
* TR_VAR:記録するトレーサー粒子の物理量
* TrT:最高衝撃温度
* TrP:最大衝撃圧力
* Trt:衝撃温度履歴
* Trp:衝撃圧力履歴
* TrV:体積ひずみ
* TrA:体積膨張
=== トレーサー粒子の描画
==== iSALEPlot上での描画
* コマンド
--------- TRACER PARAMETERS -------------------------------------------
TR_TYPE Tracer lines, grid or points (1,2,3) : -1 <- トレーサーの描画タイプ
TR_SPACE Spacing between tracers to be plotted : 10 <- トレーサーの描画間隔
TR_SIZE Size of tracer plot marker : 1 <- トレーサーのサイズ
-----------------------------------------------------------------------
* TR_TYPE:トレーサー粒子の描画タイプ
* 1 -> 「線」
* 2 -> 「グリッド」
* 3 -> 「点」
※現段階では,1-3のどれを入力しても「点」として出力される.原因は不明.
* TR_SPACE:トレーサーの描画間隔
* 0 -> 描画されない
*それ以外の整数 -> 現段階ではまだ反映されず,「1」として出力される.原因不明.
* TR_SIZE:トレーサーの大きさ
* 0 -> 描画されない
* それ以外 -> 値に従って大きくなる.
* 図
いくつかの描画例を下に示す.(左:圧力[Pa],右:温度[K])
{{attach_view(00000.png)}}
{{attach_view(00065.png)}}
{{attach_view(00112.png)}}
{{attach_view(00650.png)}}
{{attach_view(01950.png)}}
{{attach_view(02600.png)}}
{{attach_view(03250.png)}}
==== VIMoD上での描画
* パネルの説明
* 図
=== トレーサー粒子の解析
====トレーサー粒子ごとのデータの出力
iSALEPLotの描画コマンドにオプションとして「-T <描画する粒子の番号>」を追加することで.txtファイルを出力可能.
例:1番目のトレーサー粒子を出力.
../iSALEPlot -f ../iSALEPlot.inp -m jdata.dat -T 1
上記の操作により,「tracer-1000001.txt」というファイルが生成される.
中身は次の通り
Time xmark ymark TmpPre <- 左から時間,粒子のx座標,y座標,温度,圧力
0.0000000 50.000000 50.000000 753.21973 1.19651878E+09
1.5413910 206.66608 -2828.1284 753.21973 1.19651878E+09
3.0780520 258.27499 -3203.5825 753.21973 1.19651866E+09
4.6188130 261.67145 -3091.5713 753.21948 1.19648666E+09
6.1562600 270.00211 -2994.5752 753.22052 1.19664128E+09
7.6939030 276.39139 -2916.2378 753.23224 1.19834803E+09
9.2336054 279.68808 -2859.0605 754.13019 1.31637120E+09
10.770640 281.79733 -2828.4565 753.06561 1.18852262E+09
12.312140 283.42297 -2823.1621 753.07471 1.19253312E+09
13.848370 284.24585 -2840.7695 753.13367 1.20153984E+09
15.386140 271.85925 -2867.3828 753.03326 1.18787840E+09
16.924549 247.52042 -2886.7148 753.03821 1.18865293E+09
18.464939 246.17911 -2892.1021 752.95508 1.17706509E+09
20.003201 246.77846 -2901.5337 752.46234 1.10552653E+09
21.542770 245.76292 -2905.5068 753.19196 1.21384166E+09
23.077120 245.13913 -2909.4204 752.92352 1.17844365E+09
24.616541 245.11148 -2915.7446 752.61328 1.13912909E+09
26.154341 245.04880 -2924.6899 752.89874 1.18152627E+09
27.692949 245.04442 -2933.1025 752.97357 1.19366208E+09
29.233200 245.11763 -2937.0903 752.88684 1.18303552E+09
30.770411 245.26674 -2939.8745 752.48529 1.12413773E+09
32.310612 245.04799 -2940.2998 752.45294 1.12343501E+09
33.846161 244.43181 -2928.2373 752.54181 1.13237850E+09
35.388008 243.84512 -2930.9766 753.10828 1.22200755E+09
36.925468 243.14340 -2952.9873 753.79694 1.30881946E+09
38.464211 243.19815 -2970.5381 752.60376 1.15911859E+09
40.002048 243.25449 -2949.4482 752.75793 1.18333760E+09
41.543320 243.10425 -2934.8359 752.35272 1.12665702E+09
43.079689 243.10820 -2922.2207 752.18793 1.10040986E+09
44.615910 242.94418 -2914.3911 752.86719 1.20830656E+09
46.157761 242.23381 -2921.2798 752.64325 1.17738304E+09
47.693890 242.00459 -2923.3892 752.52234 1.15863322E+09
49.234959 241.93341 -2924.5425 752.70111 1.19225318E+09
50.769939 241.51730 -2938.4551 752.02716 1.08899290E+09
52.308472 241.61102 -2941.9395 751.87012 1.07313318E+09
53.849571 242.01355 -2936.9990 752.31018 1.14422451E+09
55.389370 240.86028 -2936.4121 752.29480 1.14128781E+09
56.926270 241.20326 -2933.1685 752.29950 1.14750758E+09
58.463509 241.56766 -2932.0737 752.55896 1.18605734E+09
60.003208 241.44037 -2935.7129 752.70404 1.20741389E+09
61.541512 242.22253 -2935.6523 752.05847 1.11615578E+09
* 左の3コラム:デフォルトで出力される値.
* 右の2コラム:「iSALEPlot.inp」で指定した描画パラメーター(詳細はこちらを参照)
==== 移動パス(流線)
深さ1 km付近に位置するトレーサー粒子の移動パス.
*下図は出力された各トレーサー粒子をコンパイルして描画したものである.
==== 温度・圧力履歴
衝突地点下での複数粒子の経験する温度圧力(赤点)のコンパイル
* 線:各地点での温度圧力履歴
* 衝突地点直下(青線)
* 衝突地点下850 m 〜 2950 m(黒〜灰色線)
{{attach_view(pt_path_1.png)}}