ParaViewを用いたPHITS 計算結果の3次元表示

Presentation on theme: "ParaViewを用いたPHITS 計算結果の3次元表示"— Presentation transcript:

1 ParaViewを用いたPHITS 計算結果の3次元表示
Multi-Purpose Particle and Heavy Ion Transport code System ParaViewを用いたPHITS 計算結果の3次元表示 PHITS講習会　入門実習 2016年9月改訂 title 1

2 実習目的 フリーソフトParaViewを使ってPHITSの計算結果を3次元的に可視化する方法を実習します Purpose 2
「実習」　前の基本的な話から Purpose 2

3 ParaViewのインストール Install 3 ParaViewのサイト（http://www.paraview.org）にアクセスする
「Download Latest Release」（ の項目から選択して [Download] ボタンをクリックし、インストールファイルをダウンロードする。その際，バージョンは最新版として，ファイルの種類は，単純な使用のみであれば「ParaView Binary Installers」で問題ない ダウンロードしたファイルを実行してインストール ParaViewのウェブサイト ダウンロード用サイト Install 3

4 paraview.inpの確認 基本計算条件 deposit.eps（1ページ目） geometry.eps 入射粒子：
体系： タリー： 300MeV/uの4Heビーム（半径1.0cm） ５層の同心円（鉄→アルミ→炭素→水→ポリエチレン） [t-deposit]による吸収線量分布(MeV/source) [t-gshow]によるジオメトリの描画 deposit.eps（1ページ目） geometry.eps Check Input File 4

5 paraview用にタリーを変更 geometry.vtk deposit.vtk ファイルができたことを確認する 5
PHITSの2次元表示タリーは，１方向だけメッシュを粗くする場合が多い 　(例： axis = xzのときのY方向など） ParaViewでは，3次元空間をある程度均等に分割する必要がある epsout = 1のままだと描画に時間が掛かるのでepsout = 0とする その代わりParaView用ファイル(*.vtk)を出力するオプションvtkout = 1とする [ T - Gshow ] mesh = xyz x-type = 2 nx = 5 xmin = -30. xmax = 30. y-type = 2 ny = 30 ymin = -30. ymax = 30. z-type = 2 nz = 30 zmin = -30. zmax = 30. axis = yz vtkout = 1 [ T - Gshow ] mesh = xyz x-type = 2 nx = 1 xmin = -30. xmax = 30. y-type = 2 ny = 30 ymin = -30. ymax = 30. z-type = 2 nz = 30 zmin = -30. zmax = 30. axis = yz epsout = 1 [ T - Deposit ] mesh = xyz x-type = 2 nx = 50 xmin = -25. xmax = 25. y-type = 2 ny = 50 ymin = -25. ymax = 25. z-type = 2 nz = 50 zmin = -25. zmax = 25. axis = xz vtkout = 1 [ T - Deposit ] mesh = xyz x-type = 2 nx = 50 xmin = -25. xmax = 25. y-type = 2 ny = 1 ymin = -25. ymax = 25. z-type = 2 nz = 50 zmin = -25. zmax = 25. axis = xz epsout = 1 geometry.vtk deposit.vtk ファイルができたことを確認する 5

6 ParaViewでdeposit.vtkを開く
① ② ③ 6

7 タリー範囲をアウトラインで表示 ① 左クリックしながら回転 ホイールで拡大・縮小 7

8 タリー結果を色で表示 ① ② 8 表示する形式(outline, surface, wireframeなど）を選択
表示するデータ（material, p1(all), p2(He), p3(proton), region）を選択 （ただし，一番外側の値（０）のみが表示されてしまう） 8

9 ある断面で切り出す（Clip） ①Clip ② 切断面はマウス，もしくはOrigin, Normalパラメータで調整する 9

10 しきい値を設定 ①Threshold ④ ② ③ 10 しきい値は，表示データのみならず，表示していないデータでも設定可能
②で設定したパラメータの最小値・最大値を設定 Applyボタンを押すことにより適用 10

11 色を変更 ④ ③ ①Choose Preset ②Blue to Red Rainbow 11
Blue to Red Rainbowが最もPHITSの色使いに近い 11

12 Scaleを変更 ② Rescale to Custom Range ① 12
Log表示の場合，「０」のデータがあるとWarningが出てしまう（無視してよい） 12

13 体系のWireframeを表示 ①Open → geometry.vtkを選択 ③ ②Wireframe 13
Wireframeはepsファイルを重ね合わせたように表示される（この場合はyz平面） 13

14 物質を限定して表示 ②Threshold ①Threshold1を選択 ⑤ ③ Materialを選択 ④ Min, Maxともに1に設定
どの項目の下に新しい設定を追加するか選択する必要がある 複数の条件（例えばp1とmaterial）に対してThresholdをかけることができる 14

15 設定の保存 ①File 15 Fileを選択して「Save State」を選択し，*.pvsmファイルとして保存
読み込む場合は「Load State」。ただし，現在の設定に追加されるので注意 15

16 まとめ ParaViewを使えば，PHITSのタリー結果 （mesh = xyz）を３次元的に表示することがで きる
使いこなせばアニメーションの作成も可能 詳しくはParaViewマニュアルを参照ください 「ParaView　使い方」で検索すればたくさんあります （公開してよい）おもしろい動画や画像ができま したら，ぜひ事務局までお送りください Summary 16