L ỜI MỞ ĐẦ U
4.2.5 Các thông số mô phỏng cho MCNP5
Tally
Tally là phương thức tính toán được MCNP cung cấp để thu được các kết quả
ghi nhận của các hạt truyền trong môi trường vật chất. Tùy vào mục đích tính toán
32
MCNP cũ cung cấp tất cả là 7 loại tally gồm có F1, F2, F4, F5, F6, F7, F8. Bên cạnh
đó phiên bản MCNP5 còn cung cấp thêm một loại tally mới là tally FMESH. Trong phạm vi luận văn này, chỉ cần sử dụng hai loại tally là F8 và FMESH để xác định liều hấp thụ trong phantom như sau [6]:
Tally F8: là tally cung cấp phân bố của năng lượng của bức xạphát ra được thu
nhận trong cell và được gọi là tally độ cao xung. Ở luận văn này, tally F8 được
dùng để tính liều hấp thụ của các voxel trên trục trung tâm (phân bố liều hấp thụ
theo độ sâu) của phantom nước thông qua việc tính toán năng lượng để lại trong từng voxel.
Tally FMESH:là tally tính toán ước lượng chiều dài vết của thông lượng hạt tạo
ra trong từng cell, giá trị được tính trung bình trên một mạng lưới cell, đơn vị
của tally FMESH là hạt/cm2. Ở luận văn này, tally FMESH được dùng để chia nhỏ phantom nước và phantom CT thành một mạng lưới các voxel nhỏ và tính phân bố liều lên các voxel bằng cách chuyển đổi từthông lượng sang liều thông qua bộ hệ số chuyển đổi (flux-to-dose rate conversion factors).
Số hạt (NPS)
Quy định số hạt phát ra tại cửa thoát của máy gia tốc hay nói cách khác là số
biến cốđược sử dụng trong mô phỏng Monte Carlo và được người dùng thiết lập.
Thời gian (CTME)
Quy định thời gian tính toán mô phỏng, thời gian càng nhiều số lịch sử hạt phát ra càng lớn. Trong quá trình mô phỏng MCNP, chương trình sẽ dừng lại nếu một trong hai thông sốNPS và CTME đạt tới ngưỡng quy định trước.
Kĩ thuật giảm phương sai
Nhằm mục đích giảm thời gian tính toán của chương trình mà không ảnh hưởng
đến độ chính xác của bài toán. Trong luận văn này đã áp dụng hai kỹ thuật là năng lượng cắt (energy cut-off) và phân chia bremsstrahlung (bremsstrahlung splitting).
33
Do việc mô phỏng vận chuyển electron trong vật chất tốn rất nhiều thời gian chạy chương trình nên để rút ngắn thời gian, chọn mức năng lượng cắt của electron là vào khoảng 0,7 MeV và mức năng lượng cắt của photon là vào khoảng 0,01 MeV. Nếu trong quá trình mô phỏng, năng lượng của các electron và photon nào rơi xuống dưới
các ngưỡng cắt này thì chương trình MCNP sẽ ngưng khảo sát tiếp các hạt đó. Đồng thời, cũng áp dụng phương pháp phân chia bremsstrahlung có trong chương trình MCNP. Trong phương pháp này, ở mỗi tương tác phát bức xạ bremsstrahlung, electron sẽ phát ra N photon thay vì chỉ 1 photon như phương pháp mô phỏng thông thường, mỗi photon phát ra sẽ có trọng số 1/N, chúng ta chọn số photon phát ra ởđây là N = 80
[6].