3. 例. トレーサー粒子追跡
概要
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])
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(黒〜灰色線)
Keyword(s):
References:[SideMenu]