Nghiên cứu phương pháp mô phỏng monte carlo trong tính toán che chắn an toàn bức xạ đối với bức xạ gamma

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI LẠI TIẾN THỊNH LẠI TIẾN THỊNH KỸ THUẬT HẠT NHÂN NGHIÊN CỨU PHƯƠNG PHÁP MÔ PHỎNG MONTE - CARLO TRONG TÍNH TOÁN CHE CHẮN AN TOÀN BỨC XẠ ĐỐI VỚI BỨC XẠ GAMMA LUẬN VĂN THẠC SĨ KHOA HỌC KỸ THUẬT HẠT NHÂN 2008 – 2010 Hà Nội – Năm 2011 MỤC LỤC Lời cam đoan i Danh mục bảng ii Danh mục hình vẽ, đồ thị iv Mở đầu vii Chương 1: TỔNG QUAN PHƯƠNG PHÁP TÍNH TOÁN CHE CHẮN 1.1 Các đại lượng tính toán an toàn xạ 1 1.1.1 Năng lượng 1.1.2 Hoạt độ phóng xạ 1.1.3 Liều hấp thụ xạ 1.1.4 Liều tương đương 1.1.5 Liều hiệu dụng 1.1.6 Chiều dày làm yếu nửa (HVL) 1.1.7 Chiều dày làm yếu mười lần (TVL) 1.2 Lịch sử nghiên cứu phương pháp tính toán che chắn xạ gamma 1.2.1 Những năm gần 1.2.1.1 Che chắn máy phát tia X 1.2.1.2 Che chắn nguồn phóng xạ radium 1.2.2 Dự án Manhattan thời gian đầu sau chiến tranh 1.2.2.1 Thời kỳ chiến tranh 1.2.2.2 Lò phản ứng hạt nhân cho động đẩy 1.2.3 Thập niên năm 1950 1.2.3.1 Những tiến phương pháp che chắn tia gamma 10 1.2.3.2 Những tiến phương pháp tính toán Monte Carlo 11 1.2.4 Thập niên năm 1960 12 1.2.4.1 Che chắn không gian 13 1.1.4.2 Cấu trúc che chắn 13 1.2.4.3 Các ứng dụng máy tính kỹ thuật số 14 1.2.5 Những năm thập niên 1970 15 1.2.6 Những năm 1980 1990 1.2.6.1 Những tiến hế số tích lũy 17 1.2.6.2 Tiết diện hệ số chuyển đổi liều 18 1.2.6.3 Ứng dụng máy tính 20 1.3 Che chắn xạ thực tiễn 21 1.3.1 Khái niệm 21 1.3.1.1 Đặc điểm nguồn xạ 21 1.3.1.2 Đặc tính suy giảm 22 1.3.1.3 Đặc trưng cường độ 23 1.3.1.4 Các hệ số chuyển đổi từ thông lượng thành liều 23 1.3.1.5 Liều nội 24 1.3.1.6 Liều phantom 25 1.3.2 Các phương pháp phân tích 27 1.3.2.1 Liều xạ không tương tác 27 1.3.2.2 Phương pháp điểm hạt nhân cho liều xạ không tương tác 1.3.2.3 Ứng dụng lựa chọn hình học 27 28 1.3.3 Các phương pháp trung gian dùng tính toán che chắn gamma 29 1.3.3.1 Khái niệm hệ số tích lũy 30 1.3.3.2 Hệ số tích lũy hình học 30 1.3.3.3 Tích lũy môi trường không phân lớp 31 1.3.3.4 Tích lũy môi trường phân lớp 31 1.3.3.5 Suy giảm chùm tia rộng 32 1.3.3.6 Tia tới xiên 32 1.3.3.7 Sự suy giảm chùm tia X 33 1.3.4 Tán xạ tia gamma 33 1.3.4.1 Phương pháp suất phản chiếu (Albedo method) 34 1.3.4.2 Áp dụng phương pháp suất phản chiếu 34 1.3.5 Tính xấp xỉ suất phản chiếu liều tia gamma 35 1.3.5.1 Bức xạ truyền qua ống dẫn 36 1.3.5.2 Tia gamma skyshine 37 1.3.6 Thuyết vận chuyển 39 1.3.6.1 Thuyết vận chuyển tất định 39 1.3.6.2 Lý thuyết vận chuyển Monte Carlo 41 Chương 2: PHƯƠNG PHÁP MONTE-CARLO VÀ CHƯƠNG TRÌNH MCNP5 2.1 Phương pháp Monte-Carlo 2.1.1 Giới thiệu chung 44 44 44 2.1.2 Mô đại lượng ngẫu nhiên phương pháp Mont – Carlo 44 2.2 Chương trình MCNP5 46 2.2.1 Số liệu phản ứng hạt nhân 46 2.2.2 Vật lý 47 2.2.3 Đánh giá sai số chương trình MCNP5 47 2.2.3.1 Sai số thống kê 47 2.2.3.2 Các nguồn gây sai số khác 48 Chương 3: KẾT QUẢ TÍNH TOÁN VÀ THẢO LUẬN 51 3.1 Ứng dụng mcnp đánh giá suy giảm cường độ chùm xạ qua vật liệu che chắn 51 3.1.1 Bài toán 51 3.1.2 Kết tính toán 52 3.1.3 Nhận xét 53 3.2 Ứng dụng mcnp đánh giá liều gây vật liệu xây dựng có chứa phóng xạ 53 3.2.1 Bài toán phương pháp tính 54 3.2.2.1 Bài toán 54 3.2.2.2 Phương pháp giải tích 55 3.2.2.3 Phương pháp mô sử dụng chương trình MCNP5 56 3.2.3 Kết 57 3.2.3.1 Kết tính MCNP so với kết tính phương pháp giải tích trường hợp không tính đến ảnh hưởng tán xạ 57 3.2.3.2 Kết tính chiều dày tường điển hình (20cm) trường hợp có tính đến ảnh hưởng tán xạ 58 3.2.3.3 Kết tính toán suất liều tương đương đề xuất mức giới hạn hoạt độ cho nhân phóng xạ có vật liệu xây dựng 3.2.3.4 Nhận xét 59 60 3.3 Tính toán che chắn cho máy gia tốc tuyến tính lượng 10 MeV 61 3.3.1 Bài toán 61 3.3.2 Phương pháp đánh giá 62 3.3.2.1 Chuẩn liều vị trí điểm đồng tâm 62 3.3.2.2 Tính toán suất liều vị trí xung quanh phòng dọc theo Maze 63 3.3.3 Kết tính toán 65 3.3.4 Nhận xét 67 Kết luận 68 Tài liệu tham khảo 70 Phụ lục I Kết tính toán suy giảm cường độ chùm xạ qua vật liệu che chắn 73 Phụ lục II Mã chương trình MCNP mô toán đánh giá liều từ vật liệu xây dựng có chứa phóng xạ 95 Phụ lục III Mã chương trình mô toán đánh giá suy giảm cường độ chùm xạ qua vật liệu che chắn 99 Phụ lục IV Mã chương trình MCNP mô toán đánh giá che chắn cho phòng đặt máy giá tốc phát lượng 10 MeV 104 i LỜI CAM ĐOAN Tác giả xin cam đoan công trình khoa học chưa cá nhân tổ chức công bố Tất số liệu luận văn trung thực, khách quan tác giả trực tiếp tính toán, mô chương trình MCNP5 Cục An toàn xạ hạt nhân Hà Nội, ngày 15 tháng năm 2011 Người cam đoan ký tên ii DANH MỤC CÁC BẢNG Bảng Kết tính toán chiều dày làm yếu 10 lần chương trình mô MCNP5 cho số vật liệu che chắn điển hình 52 Bảng 2: Năng lượng tia gamma trung bình, tỷ số phân nhánh () số thông số sử dụng tính toán: 55 Bảng Kết tính toán trường hợp không tính đến ảnh hưởng tán xạ 57 Bảng 4: Kết tính toán với trường hợp tính đến ảnh hưởng tán xạ tường dày 20 cm 58 Bảng Suất liều tương đương tâm phòng 59 Bảng Suất liều tương đương tổng cộng tâm phòng 59 Bảng Giới hạn hoạt độ đồng vị phóng xạ 59 Bảng Kết tính toán cho trường hợp chùm tia chiếu phía tường trước 65 Bảng Kết tính toán cho trường hợp chùm tia chiếu phía tường sau 66 Bảng 10 Kết tính toán cho trường hợp chùm tia chiếu lên trần 66 Bảng 11 Kết tính toán cho trường hợp chùm tia chiếu xuống sàn 66 Bảng 12 Tổng hợp kết tính toán cho trường hợp 67 Bảng 13 Kết tính toán nguồn I-192, vật liệu bê tông 73 Bảng 14 Kết tính toán nguồn I-192, vật liệu sắt 73 Bảng 15 Kết tính toán nguồn I-192, vật liệu chì 74 Bảng 16 Kết tính toán nguồn Cs-137, vật liệu bê tông 75 Bảng 17 Kết tính toán nguồn Cs-137, vật liệu sắt 76 Bảng 18 Kết tính toán nguồn Cs-137, vật liệu chì 77 Bảng 19 Kết tính toán nguồn Co-60, vật liệu bê tông 78 Bảng 20 Kết tính toán nguồn Co-60, vật liệu sắt 79 Bảng 21 Kết tính toán nguồn Co-60, vật liệu chì 80 Bảng 22 Kết tính toán nguồn tia X có lượng 150 keV, vật liệu bê tông 81 Bảng 23 Kết tính toán nguồn tia X có lượng 150 keV, vật liệu sắt 82 iii Bảng 24 Kết tính toán nguồn tia X có lượng 150 keV, vật liệu chì 83 Bảng 25 Kết tính toán nguồn tia X có lượng 200 keV, vật liệu bê tông Barit 84 Bảng 26 Kết tính toán nguồn tia X có lượng 200 keV, vật liệu bê tông 85 Bảng 27 Kết tính toán nguồn tia X có lượng 200 keV, vật liệu sắt 85 Bảng 28 Kết tính toán nguồn tia X có lượng 200 keV, vật liệu chì 86 Bảng 29 Kết tính toán nguồn tia X có lượng 200 keV, vật liệu bê tông Barit 87 Bảng 30 Kết tính toán nguồn photon phát từ máy gia tốc có lượng MeV, vật liệu bê tông 88 Bảng 31 Kết tính toán nguồn photon phát từ máy gia tốc có lượng MeV, vật liệu sắt 89 Bảng 32 Kết tính toán nguồn photon phát từ máy gia tốc có lượng MeV, vật liệu chì 90 Bảng 33 Kết tính toán nguồn photon phát từ máy gia tốc có lượng MeV, vật liệu bê tông 91 Bảng 34 Kết tính toán nguồn photon phát từ máy gia tốc có lượng MeV, vật liệu sắt 92 Bảng 35 Kết tính toán nguồn photon phát từ máy gia tốc có lượng MeV, vật liệu chì 93 iv DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Hình Máy tính kỹ số 12 Hình So sánh hệ số tích lũy yếu tố tích tụ tính mã Pallas phù hợp theo phương pháp cấp số nhân 19 Hình So sánh hàm đáp ứng photon Dữ liệu trích dẫn từ ICRP (1987) 26 Hình Nguồn thẳng đẳng hướng với che chắn 29 Hình Hệ số điều chỉnh cho hệ số tích lũy Bx ranh giới môi trường hữu hạn phương tiện hữu hạn nằm môi trường vô hạn có hệ số tích lũy B độ sâu thâm nhập 30 Hình Mối quan hệ góc lượng công thức suất phản chiếu 35 Hình Suất phản chiếu tương đương liều môi trường photon tán xạ có lượng 1,25 MeV từ bêtông, tính toán cách sử dụng công thức xấp xỉ bảy thông số Chilton-Huddleston 36 Hình Hình học sử dụng để tính toán khả che chắn xạ vật liệu theo chiều dày 52 Hình Hình học sử dụng để tính toán suất liều hấp thụ gamma không khí tâm phòng gây nhân phóng xạ có vật liệu xây dựng 54 Hình 10 Hình học sử dụng để tính suất liều hấp thụ không khí điểm P (xp,yp,z p) 56 Hình 11 Hình học tính toán trường hợp không tính đến tán xạ 56 Hình 12 Hình học tính toán trường hợp có tính đến ảnh hưởng tán xạ 57 Hình 13 Mặt phòng đặt máy gia tốc tuyến tính cần tính toán che chắn 62 Hình 14 Hình học toán tính toán che chắn cho máy gia tốc lượng 10 MeV – chùm tia hướng phía tường trước 63 Hình 39 Đồ thị thể suy giảm chùm tia theo chiều dày nguồn photon phát từ máy gia tốc có lượng 10 MeV, vật liệu che chắn sắt - Kết tính toán nguồn photon phát từ máy gia tốc có lượng 10 MeV, vật liệu chì đưa Bảng 30: Bảng 35 Kết tính toán nguồn photon phát từ máy gia tốc có lượng 10 MeV, vật liệu chì: STT 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 Chiều dày, cm 0.5 1.5 2.5 3.5 4.5 5.5 6.5 7.5 10 11 12 13 14 15 Kerma, MeV/g 4.18E-08 3.48E-08 2.86E-08 2.34E-08 1.94E-08 1.59E-08 1.27E-08 1.03E-08 8.37E-09 6.62E-09 5.42E-09 4.38E-09 3.51E-09 2.83E-09 2.25E-09 1.87E-09 1.55E-09 9.76E-10 6.64E-10 4.41E-10 2.87E-10 1.68E-10 1.09E-10 7.17E-11 93 Sai số 1.60% 1.74% 1.90% 2.09% 2.28% 2.52% 2.80% 3.09% 3.42% 3.84% 4.16% 4.61% 5.13% 4.01% 4.44% 4.85% 5.31% 6.55% 6.59% 7.87% 9.88% 13.14% 16.92% 19.68% Hình 40 Đồ thị thể suy giảm chùm tia theo chiều dày nguồn photon phát từ máy gia tốc có lượng 10 MeV, vật liệu che chắn chì 94 PHỤ LỤC II MÃ CHƯƠNG TRÌNH MCNP MÔ PHỎNG BÀI TOÁN ĐÁNH GIÁ LIỀU TỪ VẬT LIỆU XÂY DỰNG CÓ CHỨA PHÓNG XẠ - Không tính đến tán xạ: continue mode m1 m2 -2.35 -1 -0.001293 -2 -0.001293 -3 box 0 700 0 1200 0 12.766 s 350 600 152.766 15 s 350 600 10 3000 p 6000.04p 7000.04p 1000.04p 8000.04p 13000.04p 19000.04p -0.000124 -0.755268 -0.005558 -0.4980765 -0.045746 -0.019239 2r $Air 8000.04p $Concrete 11000.04p 14000.04p 20000.04p -0.231781 18000.04p -0.012827 -0.017101 12000.04p -0.3150923 16000.04p -0.0829411 26000.04p $ 1, -0.002565 -0.001283 -0.012398 imp:p phys:p nps 5e7 c ctme sdef x=d1 y=d2 z=d3 erg=d4 cel=1 par=2 si1 700 sp1 si2 1200 sp2 si3 12.766 sp3 c si4 L 0.0532 0.2420 0.2952 0.3519 0.7860 0.6093 0.7684 0.8062 0.9341 c 1.1552 1.2381 1.2810 1.3777 1.4015 1.4080 1.5092 1.7296 1.7645 c 2.1185 2.2042 2.4478 $ U-238 c sp4 0.0110 0.0712 0.1820 0.3510 0.0104 0.4460 0.0476 0.0120 0.0307 c 0.0190 0.0578 0.0160 0.0480 0.0140 0.0250 0.0208 0.0290 0.1510 c 0.0117 0.0498 0.0155 c si4 L 0.0995 0.1291 0.2093 0.2702 0.3276 0.3383 0.4098 0.4630 0.7721 c 0.8356 0.9111 0.9644 0.9688 1.5883 1.6307 0.2410 0.7273 0.7855 c 0.2386 0.3001 0.2774 0.5834 0.8604 2.6145 $Th-232 c sp4 0.0130 0.0223 0.0381 0.0344 0.0310 0.1126 0.0195 0.0450 0.0145 c 0.0153 0.2660 0.0505 0.1623 0.0326 0.0153 0.0404 0.0665 0.0107 c 0.4350 0.0327 0.0230 0.3060 0.0455 0.3560 c si4 L 1.461 $K-40 c sp4 c si4 L 0.031817 0.032194 0.036304 0.036378 0.037255 0.2835 0.661657 $ c sp4 0.019600 0.035800 0.003420 0.006600 0.002040 0.0000058 0.8510 c si4 L 0.662 $ Cs-137 c sp4 0.852 si4 L 0.810 $ U-238 sp4 c si4 L 0.587 2.615 $ Th-232 c sp4 0.859179 0.140821 f6:p 1.1203 1.8474 0.1470 0.0220 0.7948 1.6206 0.0434 0.0138 Cs-137 Tính đến tán xạ: continue 95 c Created on: Wednesday, April 18, 2007 at 15:26 -2.35 -1 :-2 :-3 :-4 :-5 :-6 -0.001293 -7 -0.001293 -9 -0.001293 -8 #1 8 mode m1 m2 box box box box box box s s box p 6000.04p 7000.04p 1000.04p 8000.04p 13000.04p 19000.04p 0 700 0 -21.2766 0 280 $W2 1200 700 0 21.2766 0 280 $W2 0 700 0 1200 0 -21.2766 $Floor-ceiling 0 280 700 0 1200 0 21.2766 $Floor-ceiling 0 -21.2766 0 1200 0 280 $W1 700 0 21.2766 0 1200 0 280 $W1 350 600 140 15 350 600 140 3000 0 700 0 1200 0 280 -0.000124 -0.755268 -0.005558 -0.4980765 -0.045746 -0.019239 3r $Air 8000.04p $Concrete 11000.04p 14000.04p 20000.04p -0.231781 18000.04p -0.012827 -0.017101 12000.04p -0.3150923 16000.04p -0.0829411 26000.04p $ 1, -0.002565 -0.001283 -0.012398 imp:p phys:p nps 5e7 c ctme 180 sdef x=d1 y=d2 z=d3 erg=d4 cel=1 par=2 si1 700 sp1 si2 1200 sp2 si3 -21.2766 sp3 c si4 L 0.0532 0.2420 0.2952 0.3519 0.7860 0.6093 0.7684 0.8062 0.9341 1.1203 c 1.1552 1.2381 1.2810 1.3777 1.4015 1.4080 1.5092 1.7296 1.7645 1.8474 c 2.1185 2.2042 2.4478 $ U-238 c sp4 0.0110 0.0712 0.1820 0.3510 0.0104 0.4460 0.0476 0.0120 0.0307 0.1470 c 0.0190 0.0578 0.0160 0.0480 0.0140 0.0250 0.0208 0.0290 0.1510 0.0220 c 0.0117 0.0498 0.0155 si4 L 0.0995 0.1291 0.2093 0.2702 0.3276 0.3383 0.4098 0.4630 0.7721 0.7948 0.8356 0.9111 0.9644 0.9688 1.5883 1.6307 0.2410 0.7273 0.7855 1.6206 0.2386 0.3001 0.2774 0.5834 0.8604 2.6145 $Th-232 sp4 0.0130 0.0223 0.0381 0.0344 0.0310 0.1126 0.0195 0.0450 0.0145 0.0434 0.0153 0.2660 0.0505 0.1623 0.0326 0.0153 0.0404 0.0665 0.0107 0.0138 0.4350 0.0327 0.0230 0.3060 0.0455 0.3560 c si4 L 1.461 $K-40 c sp4 c si4 L 0.662 $ Cs-137 c sp4 0.852 c si4 L 0.810 $ U-238 c sp4 c si4 L 0.587 2.615 $ Th-232 c sp4 0.859179 0.140821 f6:p Sử dụng cầu chuẩn ICRU để đánh giá liều hiệu dụng: continue 1 1 -1 -1 -1 -1 -1 -3 -4 -5 -6 -7 96 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 21 22 23 30 31 32 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 -1 -8 -1 -9 -1 -10 -1 -11 -1 -12 -1 -13 -1 -14 -1 -15 -1 -16 -1 -17 -1 -18 -1 -19 -1 -20 -1 -21 -1 -22 -1 -23 -1 -24 -1 -25 -1 -26 -1 -27 -1 -28 -1 -29 -0.001293 -1 10 11 12 13 14 15 16 17 18 19 20 24 25 26 27 28 29 31 -0.001293 -30 30 -2.35 -31 box box s s s s s s s s s s s s s s s s s s s s s s s s s s s s box 0 700 0 1200 0 280 -21.27 -21.27 -21.27 742.54 0 1242.54 0 322.54 110 200 40 15 110 600 40 15 110 1000 40 15 350 200 40 15 350 600 40 15 350 1000 40 15 590 200 40 15 590 600 40 15 590 1000 40 15 110 200 140 15 110 600 140 15 110 1000 140 15 350 200 140 15 350 600 140 15 350 1000 140 15 590 200 140 15 590 600 140 15 590 1000 140 15 110 200 240 15 110 600 240 15 110 1000 240 15 350 200 240 15 350 600 240 15 350 1000 240 15 590 200 240 15 590 600 240 15 590 1000 240 15 350 600 10 3000 0 700 0 1200 0 97 mode m1 p 1000.04p 6000.04p -0.101 $ Phantom ICRU -0.111 7000.04p -0.026 8000.04p 0.762 m2 1000.04p -0.005558 $Concrete 8000.04p -0.4980765 11000.04p -0.017101 12000.04p 0.002565 13000.04p -0.045746 14000.04p -0.3150923 16000.04p 0.001283 19000.04p -0.019239 20000.04p -0.0829411 26000.04p 0.012398 m3 6000.04p -0.000124 $Air 7000.04p -0.755268 8000.04p -0.231781 18000.04p 0.012827 imp:p 29r $ 1, 32 phys:p nps 2e7 c ctme 180 sdef x=d1 y=d2 z=d3 erg=d4 cel=32 par=2 si1 700 sp1 si2 1200 sp2 si3 sp3 c si4 L 0.0633 0.0924 0.0928 0.1862 0.0533 0.2420 0.2952 0.3519 0.6093 c 0.6655 0.7684 0.8062 0.9341 1.1203 1.1552 1.2381 1.2810 1.3777 c 1.4015 1.4080 1.5092 1.6613 1.7296 1.7645 1.8474 2.1185 2.2040 c 2.4477 0.0465 $U-238,Emission probability total = 1.964 c sp4 0.0370 0.0262 0.0259 0.0353 0.0107 0.0719 0.1828 0.3534 0.4516 c 0.0152 0.0485 0.0126 0.0307 0.1478 0.0162 0.0579 0.0143 0.0395 c 0.0132 0.0237 0.0211 0.0104 0.0282 0.1517 0.0200 0.0115 0.0489 c 0.0154 0.0425 c si4 L 0.0844 0.2410 0.2386 0.3001 0.7273 0.7854 1.6207 0.2774 0.5832 c 0.8606 2.6145 0.0995 0.1291 0.2093 0.2702 0.3274 0.3383 0.4095 c 0.4630 0.7553 0.7723 0.7949 0.8357 0.9112 0.9648 0.9690 1.5882 c 1.6306 $Th-232, Emission probability total = 2.261 c sp4 0.0117 0.0412 0.4360 0.0318 0.0674 0.0111 0.0151 0.0237 0.3055 c 0.0449 0.3586 0.0126 0.0242 0.0389 0.0346 0.0307 0.1127 0.0192 c 0.0440 0.0100 0.0149 0.0425 0.0161 0.2580 0.0499 0.1580 0.0322 c 0.0151 si4 L 1.461 $K-40, Emission probability total = 0.107 sp4 c si4 L 0.810 $U-238, Emission probability total = 2.12 c sp4 c si4 L 0.587 2.615 $Th-232, Emission probability total = 2.406 c sp4 2.05 0.356 f6:p 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 98 PHỤ LỤC III MÃ CHƯƠNG TRÌNH MCNP MÔ PHỎNG BÀI TOÁN ĐÁNH GIÁ SỰ SUY GIẢM CƯỜNG ĐỘ CHÙM BỨC XẠ KHI ĐI QUA VẬT LIỆU CHE CHẮN - Nguồn xạ Ir-192: continue c Created on: Wednesday, April 18, 2007 at 15:26 -2.35 -1 -0.001293 -2 -0.001293 -3 -0.001293 -4 -0.001293 -5 -0.001293 -6 -0.001293 -7 -0.001293 -8 -0.001293 -9 10 -0.001293 -10 11 -0.001293 -11 12 -0.001293 -12 13 -0.001293 -13 14 -0.001293 -14 15 -0.001293 -15 16 -0.001293 -16 17 -0.001293 -17 18 -0.001293 -18 19 -0.001293 -19 20 -0.001293 -20 21 -0.001293 -21 22 -0.001293 -22 10 11 12 13 14 15 16 17 18 19 20 21 23 22 10 11 12 13 14 15 16 17 18 19 20 21 22 mode m1 m2 m3 box s s s s s s s s s s s s s s s s s s s s s p 6000.04p 7000.04p 1000.04p 8000.04p 13000.04p 19000.04p 6000.04p 0 65 0 2000 0 2000 -15 950 1000 15 -15 1000 1000 15 -15 1050 1000 15 -15 1100 1000 15 -15 1150 1000 15 -15 1200 1000 15 -15 1250 1000 15 -15 1300 1000 15 -15 1350 1000 15 -15 1400 1000 15 150 1000 1000 15 150 1100 1000 15 150 1200 1000 15 150 1300 1000 15 150 1400 1000 15 150 1500 1000 15 150 1600 1000 15 150 1700 1000 15 150 1800 1000 15 150 1900 1000 15 10 1000 1000 1500 -0.000124 -0.755268 -0.005558 -0.4980765 -0.045746 -0.019239 -0.0006 $Air 8000.04p -0.231781 $Concrete 11000.04p -0.017101 14000.04p -0.3150923 20000.04p -0.0829411 $carbon steel 99 18000.04p -0.012827 12000.04p 16000.04p 26000.04p -0.002565 -0.001283 -0.012398 25000.04p -0.0035 26000.04p -0.9959 m4 82000.04p -1 $Pb imp:p 21r $ 1, 23 phys:p nps 1e7 c ctme c sdef pos=300 1000 1000 erg=d1 par=2 $ nguon Ir-191 sdef pos=300 1000 1000 erg=d1 vec=-1 0 dir=d2 par=2 si1 L 0.06149 0.06300 0.07130 0.07340 0.11009 $ nguon Ir-192, 81.91Ci 0.13634 0.17698 0.20131 0.20580 0.28004 $ he so nhan 2.24343427 0.28327 0.29596 0.30847 0.31651 0.32931 0.37449 0.41647 0.42053 0.46807 0.48458 0.48530 0.48904 0.58859 0.59337 0.59940 0.60442 0.61247 0.70398 0.76600 0.88454 1.06148 1.08970 1.37830 sp1 7.13192e-3 9.04863e-3 2.95484e-3 7.72031e-4 5.65651e-5 8.18388e-4 1.91136e-5 2.10347e-3 1.47230e-2 1.03680e-4 1.17097e-3 1.27795e-1 1.45714e-1 3.69344e-1 8.27303e-5 3.21293e-3 2.69153e-3 3.27088e-4 2.13200e-1 1.42059e-2 9.80639e-6 1.97599e-3 2.01254e-2 1.89531e-4 1.72860e-5 3.66937e-2 2.36646e-2 2.38206e-5 6.65052e-6 1.30113e-3 2.35353e-4 4.78285e-6 5.53616e-6 si2 -1 0.94 sp2 0.97 0.03 sb2 0 f6:p 10 11 12 13 14 15 16 17 18 19 20 21 Mô tả nguồn Cs-137: sdef pos=300 1000 1000 erg=d1 vec=-1 0 dir=d2 par=2 $ nguon Cs-137 si1 L 0.66162 0.03219 0.03182 0.03640 $ nguon Cs-137 he so nhan 0.9128 sp1 8.4620e-1 3.7000e-2 1.92004e-2 1.0400e-2 si2 -1 0.94 sp2 0.97 0.03 sb2 0 Mô tả nguồn Co-60: sdef pos=300 1000 1000 erg=d1 vec=-1 0 dir=d2 par=2 $ nguon Co-60 si1 L 1.33251 1.17323$ nguon Co-60 he so nhan 1.9984 sp1 9.998e-1 9.986e-1 si2 -1 0.94 sp2 0.97 0.03 sb2 0 Mô tả nguồn tia X lượng 150 keV: c sdef pos=300 1000 1000 erg=d1 par=2 $ nguon tia X dai sdef pos=300 1000 1000 erg=d1 vec=-1 0 dir=d2 par=2 $ 150keV si1 L 4.950E-02 5.000E-02 5.050E-02 5.100E-02 5.150E-02 5.300e-02 5.350E-02 5.400E-02 5.450E-02 5.500E-02 5.650e-02 5.700E-02 5.750E-02 5.800E-02 5.850E-02 6.000e-02 6.050E-02 6.100E-02 6.150E-02 6.200E-02 6.350e-02 6.400E-02 6.450E-02 6.500E-02 6.550E-02 6.700e-02 6.750E-02 6.800E-02 6.850E-02 6.900E-02 100 hep 150keV nguon tia X dai hep 5.200E-02 5.550E-02 5.900E-02 6.250E-02 6.600E-02 6.950E-02 5.250E-02 5.600E-02 5.950E-02 6.300E-02 6.650E-02 7.000E-02 sp1 si2 sp2 sb2 7.050e-02 7.400e-02 7.750e-02 8.100e-02 8.450e-02 8.800e-02 9.150e-02 9.500e-02 9.850e-02 1.020e-01 1.055e-01 1.090e-01 1.125e-01 1.160e-01 1.195e-01 1.230e-01 1.265e-01 1.300e-01 1.335e-01 1.370e-01 1.405e-01 1.440e-01 1.475e-01 1.510e-01 D 0.000E+00 4.403e-01 3.528e-01 4.828e+00 2.315e+00 7.779e+00 4.881e+00 6.817e+00 8.671e+00 1.089e+01 1.327e+01 1.387e+01 1.450e+01 1.510e+01 1.599e+01 1.692e+01 1.685e+01 1.658e+01 1.581e+01 1.613e+01 1.545e+01 1.393e+01 1.325e+01 1.105e+01 9.339e+00 7.957e+00 6.779e+00 4.493e+00 2.598e+00 0.000e+00 -1 0.94 0.97 0.03 0 7.100E-02 7.450E-02 7.800E-02 8.150E-02 8.500E-02 8.850E-02 9.200E-02 9.550E-02 9.900E-02 1.025E-01 1.060E-01 1.095E-01 1.130E-01 1.165E-01 1.200E-01 1.235E-01 1.270E-01 1.305E-01 1.340E-01 1.375E-01 1.410E-01 1.445E-01 1.480E-01 7.150E-02 7.500E-02 7.850E-02 8.200E-02 8.550E-02 8.900E-02 9.250E-02 9.600E-02 9.950E-02 1.030E-01 1.065E-01 1.100E-01 1.135E-01 1.170E-01 1.205E-01 1.240E-01 1.275E-01 1.310E-01 1.345E-01 1.380E-01 1.415E-01 1.450E-01 1.485E-01 7.200E-02 7.550E-02 7.900E-02 8.250E-02 8.600E-02 8.950E-02 9.300E-02 9.650E-02 1.000E-01 1.035E-01 1.070E-01 1.105E-01 1.140E-01 1.175E-01 1.210E-01 1.245E-01 1.280E-01 1.315E-01 1.350E-01 1.385E-01 1.420E-01 1.455E-01 1.490E-01 7.250E-02 7.600E-02 7.950E-02 8.300E-02 8.650E-02 9.000E-02 9.350E-02 9.700E-02 1.005E-01 1.040E-01 1.075E-01 1.110E-01 1.145E-01 1.180E-01 1.215E-01 1.250E-01 1.285E-01 1.320E-01 1.355E-01 1.390E-01 1.425E-01 1.460E-01 1.495E-01 7.300E-02 7.650E-02 8.000E-02 8.350E-02 8.700E-02 9.050E-02 9.400E-02 9.750E-02 1.010E-01 1.045E-01 1.080E-01 1.115E-01 1.150E-01 1.185E-01 1.220E-01 1.255E-01 1.290E-01 1.325E-01 1.360E-01 1.395E-01 1.430E-01 1.465E-01 1.500E-01 7.350E-02 7.700E-02 8.050E-02 8.400E-02 8.750E-02 9.100E-02 9.450E-02 9.800E-02 1.015E-01 1.050E-01 1.085E-01 1.120E-01 1.155E-01 1.190E-01 1.225E-01 1.260E-01 1.295E-01 1.330E-01 1.365E-01 1.400E-01 1.435E-01 1.470E-01 1.505E-01 3.499E-02 3.645E-01 8.281E-01 1.662E+00 2.417E+00 1.499E+01 5.079E+00 6.750E+00 8.470E+00 1.086E+01 1.320E+01 1.358E+01 1.483E+01 1.634E+01 1.596E+01 1.687E+01 1.646E+01 1.674E+01 1.626E+01 1.561E+01 1.526E+01 1.372E+01 1.261E+01 1.112E+01 9.823E+00 8.208E+00 6.199E+00 4.502E+00 2.003E+00 1.400E-01 5.569E-01 9.359E-01 1.458E+00 2.618E+00 6.173E+00 5.155E+00 7.592E+00 9.120E+00 1.106E+01 1.316E+01 1.385E+01 1.479E+01 1.599E+01 1.653E+01 1.672E+01 1.714E+01 1.650E+01 1.588E+01 1.556E+01 1.479E+01 1.329E+01 1.232E+01 1.106E+01 8.709E+00 8.193E+00 6.202E+00 3.869E+00 1.691E+00 2.682E-01 5.802E-01 2.429E+00 1.621E+00 2.846E+00 4.388E+00 5.502E+00 8.027E+00 9.584E+00 1.144E+01 1.358E+01 1.413E+01 1.527E+01 1.564E+01 1.674E+01 1.735E+01 1.670E+01 1.721E+01 1.622E+01 1.579E+01 1.430E+01 1.293E+01 1.240E+01 1.048E+01 9.339E+00 7.228E+00 5.645E+00 3.630E+00 1.420E+00 2.566E-01 6.181E-01 3.770E+00 1.735E+00 3.032E+00 6.537E+00 5.992E+00 7.546E+00 9.750E+00 1.201E+01 1.355E+01 1.417E+01 1.525E+01 1.575E+01 1.607E+01 1.700E+01 1.673E+01 1.641E+01 1.616E+01 1.509E+01 1.460E+01 1.329E+01 1.209E+01 1.020E+01 9.321E+00 7.178E+00 5.584E+00 3.376E+00 1.079E+00 3.324E-01 7.464E-01 1.639E+00 2.114E+00 3.373E+00 8.852E+00 6.528E+00 7.887E+00 9.695E+00 1.213E+01 1.397E+01 1.426E+01 1.538E+01 1.601E+01 1.646E+01 1.637E+01 1.667E+01 1.636E+01 1.603E+01 1.560E+01 1.455E+01 1.309E+01 1.205E+01 1.019E+01 9.097E+00 6.770E+00 4.895E+00 3.015E+00 6.677E-01 2.828E-01 2.391E-01 8.730E+00 2.271E+00 3.598E+00 5.467E+00 6.485E+00 8.257E+00 1.067E+01 1.242E+01 1.415E+01 1.458E+01 1.582E+01 1.621E+01 1.717E+01 1.649E+01 1.706E+01 1.686E+01 1.578E+01 1.548E+01 1.396E+01 1.311E+01 1.173E+01 9.910E+00 8.744E+00 6.613E+00 4.858E+00 2.607E+00 2.712E-01 Mô tả nguồn tia X lượng 200 keV: sdef pos=300 1000 1000 erg=d1 vec=-1 0 dir=d2 par=2 $ 200keV si1 L 1.030E-01 1.035E-01 1.040E-01 1.045E-01 1.050E-01 1.065e-01 1.070E-01 1.075E-01 1.080E-01 1.085E-01 1.100e-01 1.105E-01 1.110E-01 1.115E-01 1.120E-01 101 nguon tia X dai hep 1.055E-01 1.060E-01 1.090E-01 1.095E-01 1.125E-01 1.130E-01 sp1 si2 sp2 sb2 1.135e-01 1.170e-01 1.205e-01 1.240e-01 1.275e-01 1.310e-01 1.345e-01 1.380e-01 1.415e-01 1.450e-01 1.485e-01 1.520e-01 1.555e-01 1.590e-01 1.625e-01 1.660e-01 1.695e-01 1.730e-01 1.765e-01 1.800e-01 1.835e-01 1.870e-01 1.905e-01 1.940e-01 1.975e-01 2.010e-01 D 0.000E+00 5.033e-01 8.266e-01 1.121e+00 1.205e+00 1.952e+00 2.658e+00 3.556e+00 4.784e+00 6.183e+00 7.453e+00 9.036e+00 1.100e+01 1.259e+01 1.414e+01 1.576e+01 1.700e+01 1.829e+01 1.937e+01 1.933e+01 1.953e+01 1.915e+01 1.886e+01 1.722e+01 1.499e+01 1.342e+01 8.565e+00 4.455e+00 5.172e-01 -1 0.94 0.97 0.03 0 1.140E-01 1.175E-01 1.210E-01 1.245E-01 1.280E-01 1.315E-01 1.350E-01 1.385E-01 1.420E-01 1.455E-01 1.490E-01 1.525E-01 1.560E-01 1.595E-01 1.630E-01 1.665E-01 1.700E-01 1.735E-01 1.770E-01 1.805E-01 1.840E-01 1.875E-01 1.910E-01 1.945E-01 1.980E-01 2.015E-01 1.847E-01 4.664E-01 7.820E-01 1.191E+00 1.375E+00 2.192E+00 2.846E+00 3.597E+00 5.115E+00 6.300E+00 7.353E+00 9.103E+00 1.068E+01 1.289E+01 1.493E+01 1.644E+01 1.716E+01 1.847E+01 1.897E+01 2.023E+01 2.070E+01 1.884E+01 1.845E+01 1.731E+01 1.531E+01 1.183E+01 8.321E+00 3.885E+00 2.124E-01 1.145E-01 1.180E-01 1.215E-01 1.250E-01 1.285E-01 1.320E-01 1.355E-01 1.390E-01 1.425E-01 1.460E-01 1.495E-01 1.530E-01 1.565E-01 1.600E-01 1.635E-01 1.670E-01 1.705E-01 1.740E-01 1.775E-01 1.810E-01 1.845E-01 1.880E-01 1.915E-01 1.950E-01 1.985E-01 2.020E-01 2.263E-01 5.218E-01 8.743E-01 1.288E+00 1.662E+00 2.164E+00 2.946E+00 3.853E+00 5.221E+00 6.428E+00 7.884E+00 9.964E+00 1.146E+01 1.291E+01 1.497E+01 1.668E+01 1.748E+01 1.809E+01 1.991E+01 1.975E+01 2.052E+01 2.016E+01 1.836E+01 1.660E+01 1.487E+01 1.174E+01 7.846E+00 3.240E+00 1.185E-01 1.150E-01 1.185E-01 1.220E-01 1.255E-01 1.290E-01 1.325E-01 1.360E-01 1.395E-01 1.430E-01 1.465E-01 1.500E-01 1.535E-01 1.570E-01 1.605E-01 1.640E-01 1.675E-01 1.710E-01 1.745E-01 1.780E-01 1.815E-01 1.850E-01 1.885E-01 1.920E-01 1.955E-01 1.990E-01 2.025E-01 2.940E-01 5.464E-01 9.236E-01 1.415E+00 1.838E+00 2.189E+00 3.142E+00 3.908E+00 5.461E+00 6.853E+00 8.140E+00 9.905E+00 1.127E+01 1.346E+01 1.559E+01 1.661E+01 1.730E+01 1.904E+01 1.945E+01 1.921E+01 1.985E+01 1.905E+01 1.801E+01 1.658E+01 1.529E+01 1.120E+01 6.893E+00 2.666E+00 0.000E+00 1.155E-01 1.190E-01 1.225E-01 1.260E-01 1.295E-01 1.330E-01 1.365E-01 1.400E-01 1.435E-01 1.470E-01 1.505E-01 1.540E-01 1.575E-01 1.610E-01 1.645E-01 1.680E-01 1.715E-01 1.750E-01 1.785E-01 1.820E-01 1.855E-01 1.890E-01 1.925E-01 1.960E-01 1.995E-01 1.160E-01 1.195E-01 1.230E-01 1.265E-01 1.300E-01 1.335E-01 1.370E-01 1.405E-01 1.440E-01 1.475E-01 1.510E-01 1.545E-01 1.580E-01 1.615E-01 1.650E-01 1.685E-01 1.720E-01 1.755E-01 1.790E-01 1.825E-01 1.860E-01 1.895E-01 1.930E-01 1.965E-01 2.000E-01 1.165E-01 1.200E-01 1.235E-01 1.270E-01 1.305E-01 1.340E-01 1.375E-01 1.410E-01 1.445E-01 1.480E-01 1.515E-01 1.550E-01 1.585E-01 1.620E-01 1.655E-01 1.690E-01 1.725E-01 1.760E-01 1.795E-01 1.830E-01 1.865E-01 1.900E-01 1.935E-01 1.970E-01 2.005E-01 3.694E-01 6.096E-01 9.836E-01 1.447E+00 1.826E+00 2.389E+00 3.316E+00 4.308E+00 5.458E+00 6.487E+00 8.691E+00 1.046E+01 1.178E+01 1.370E+01 1.510E+01 1.686E+01 1.720E+01 1.893E+01 2.012E+01 1.980E+01 1.911E+01 1.865E+01 1.787E+01 1.590E+01 1.335E+01 1.072E+01 6.390E+00 1.889E+00 3.217E-01 7.081E-01 9.913E-01 1.387E+00 1.778E+00 2.387E+00 3.416E+00 4.567E+00 5.879E+00 6.885E+00 8.879E+00 1.035E+01 1.190E+01 1.369E+01 1.616E+01 1.701E+01 1.823E+01 1.852E+01 1.945E+01 2.016E+01 2.005E+01 1.922E+01 1.786E+01 1.545E+01 1.340E+01 1.011E+01 5.846E+00 1.379E+00 4.110E-01 6.881E-01 1.056E+00 1.599E+00 1.967E+00 2.463E+00 3.499E+00 4.749E+00 5.723E+00 7.319E+00 9.191E+00 1.027E+01 1.263E+01 1.432E+01 1.566E+01 1.684E+01 1.854E+01 1.900E+01 1.964E+01 2.003E+01 1.978E+01 2.003E+01 1.763E+01 1.567E+01 1.370E+01 9.283E+00 5.406E+00 7.466E-01 Mô tả nguồn photon lượng MeV: c sdef pos=300 1000 1000 erg=d1 par=2 $ nguon MGT 6MeV sdef pos=300 1000 1000 erg=d1 vec=-1 0 dir=d2 par=2 si1 L 0.25 0.35 0.45 0.55 0.65 0.75 0.85 0.95 102 sp1 si2 sp2 sb2 1.05 1.15 1.25 1.35 1.45 1.55 1.65 1.75 1.85 1.95 2.05 2.15 2.25 2.35 2.45 2.55 2.65 2.75 2.85 2.95 3.05 3.15 3.25 3.35 3.45 3.55 3.65 3.75 3.85 3.95 4.05 4.15 4.25 4.35 4.45 4.55 4.65 4.75 4.85 4.95 5.05 5.15 5.25 5.35 5.45 5.55 5.65 5.75 5.85 5.95 6.05 2.047E-2 2.373E-2 2.608E-2 2.778E-2 2.897E-2 3.044e-2 3.027E-2 2.996E-2 2.953E-2 2.901E-2 2.631e-2 2.554E-2 2.475E-2 2.396E-2 2.315E-2 1.997e-2 1.920E-2 1.845E-2 1.771E-2 1.698E-2 1.429e-2 1.366E-2 1.306E-2 1.248E-2 1.192E-2 9.872e-3 9.409E-3 8.966E-3 8.540E-3 8.133E-3 6.669e-3 6.342E-3 6.030E-3 5.732E-3 5.447E-3 4.434e-3 4.210E-3 3.996E-3 -1 0.94 0.97 0.03 0 2.976E-2 2.841E-2 2.235E-2 1.628E-2 1.137E-2 7.742E-3 5.176E-3 3.023E-2 2.775E-2 2.155E-2 1.560E-2 1.085E-2 7.368E-3 4.917E-3 3.044E-2 2.705E-2 2.076E-2 1.493E-2 1.035E-3 7.011E-3 4.670E-3 Mô tả nguồn photon lượng 10 MeV: c sdef pos=300 1000 1000 erg=d1 par=2 $ nguon MGT 10MeV sdef pos=300 1000 1000 erg=d1 vec=-1 0 dir=d2 par=2 si1 L 0.25 0.35 0.45 0.55 0.65 0.75 0.85 0.95 1.05 1.15 1.25 1.35 1.45 1.55 1.65 1.75 1.85 1.95 2.05 2.15 2.25 2.35 2.45 2.55 2.65 2.75 2.85 2.95 3.05 3.15 3.25 3.35 3.45 3.55 3.65 3.75 3.85 3.95 4.05 4.15 4.25 4.35 4.45 4.55 4.65 4.75 4.85 4.95 5.05 5.15 5.25 5.35 5.45 5.55 5.65 5.75 5.85 5.95 6.05 6.15 6.25 6.35 6.45 6.55 6.65 6.75 6.85 6.95 7.05 7.15 7.25 7.35 7.45 7.55 7.65 7.75 7.85 7.95 8.05 8.15 8.25 8.35 8.45 8.55 8.65 8.75 8.85 8.95 9.05 9.15 9.25 9.35 9.45 9.55 9.65 9.75 9.85 9.95 10.05 sp1 1.924E-2 2.230E-2 2.451E-2 2.610E-2 2.722E-2 2.797E-2 2.861e-2 2.845E-2 2.815E-2 2.775E-2 2.726E-2 2.670E-2 2.472e-2 2.400E-2 2.326E-2 2.251E-2 2.176E-2 2.101E-2 1.877e-2 1.805E-2 1.734E-2 1.664E-2 1.593E-2 1.530E-2 1.343e-2 1.284E-2 1.227E-2 1.173E-2 1.120E-2 1.069E-2 9.277e-3 8.843E-3 8.426E-3 8.026E-3 7.643E-3 7.276E-3 6.267e-3 5.960E-3 5.667E-3 5.387E-3 5.119E-3 4.864E-3 4.167e-3 3.956E-3 3.756E-3 3.564E-3 3.382E-3 3.209E-3 2.738e-3 2.597E-3 2.462E-3 2.334E-3 2.212E-3 2.096E-3 1.783e-3 1.689E-3 1.600E-3 1.515E-3 1.435E-3 1.359E-3 1.153e-3 1.091E-3 1.033E-3 9.770E-4 9.250E-4 8.750E-4 7.410e-4 7.010E-4 6.630E-4 6.270E-4 5.930E-4 5.610E-4 4.740e-4 4.480E-4 4.230E-4 si2 -1 0.94 sp2 0.97 0.03 sb2 0 2.841E-2 2.608E-2 2.025E-2 1.466E-2 1.020E-2 6.925E-3 4.621E-3 3.044E-3 1.986E-3 1.286E-3 8.280E-4 5.300E-4 2.861E-2 2.542E-2 1.951E-2 1.403E-2 9.730E-3 6.589E-3 4.388E-3 2.888E-3 1.882E-3 1.218E-3 7.830E-4 5.010E-4 103 PHỤ LỤC IV MÃ CHƯƠNG TRÌNH MCNP MÔ PHỎNG BÀI TOÁN ĐÁNH GIÁ CHE CHẮN CHO PHÒNG ĐẶT MÁY GIA TỐC PHÁT NĂNG LƯỢNG 10 MeV - Mã chương trình chuẩn liều vị trí điểm đồng tâm: c Created on: Sunday, October 14, 2007 at 16:08 -0.001293 -1 -0.001293 -2 s 0 100 15 s 0 1000 mode m1 p 6000.04p -0.000124 $Air 7000.04p -0.755268 8000.04p -0.231781 18000.04p -0.012827 imp:p 1r $ 1, phys:p nps 2e6 sdef pos=0 0 erg=d1 vec=0 dir=d2 par=2$ nguon may gia toc 10MeV si1 L 0.25 0.35 0.45 0.55 0.65 0.75 0.85 0.95 1.05 1.15 1.25 1.35 1.45 1.55 1.65 1.75 1.85 1.95 2.05 2.15 2.25 2.35 2.45 2.55 2.65 2.75 2.85 2.95 3.05 3.15 3.25 3.35 3.45 3.55 3.65 3.75 3.85 3.95 4.05 4.15 4.25 4.35 4.45 4.55 4.65 4.75 4.85 4.95 5.05 5.15 5.25 5.35 5.45 5.55 5.65 5.75 5.85 5.95 6.05 6.15 6.25 6.35 6.45 6.55 6.65 6.75 6.85 6.95 7.05 7.15 7.25 7.35 7.45 7.55 7.65 7.75 7.85 7.95 8.05 8.15 8.25 8.35 8.45 8.55 8.65 8.75 8.85 8.95 9.05 9.15 9.25 9.35 9.45 9.55 9.65 9.75 9.85 9.95 10.05 sp1 1.924E-2 2.230E-2 2.451E-2 2.610E-2 2.722E-2 2.797E-2 2.841E-2 2.861E-2 2.861e-2 2.845E-2 2.815E-2 2.775E-2 2.726E-2 2.670E-2 2.608E-2 2.542E-2 2.472e-2 2.400E-2 2.326E-2 2.251E-2 2.176E-2 2.101E-2 2.025E-2 1.951E-2 1.877e-2 1.805E-2 1.734E-2 1.664E-2 1.593E-2 1.530E-2 1.466E-2 1.403E-2 1.343e-2 1.284E-2 1.227E-2 1.173E-2 1.120E-2 1.069E-2 1.020E-2 9.730E-3 9.277e-3 8.843E-3 8.426E-3 8.026E-3 7.643E-3 7.276E-3 6.925E-3 6.589E-3 6.267e-3 5.960E-3 5.667E-3 5.387E-3 5.119E-3 4.864E-3 4.621E-3 4.388E-3 4.167e-3 3.956E-3 3.756E-3 3.564E-3 3.382E-3 3.209E-3 3.044E-3 2.888E-3 104 2.738e-3 2.597E-3 1.986E-3 1.882E-3 1.783e-3 1.689E-3 1.286E-3 1.218E-3 1.153e-3 1.091E-3 8.280E-4 7.830E-4 7.410e-4 7.010E-4 5.300E-4 5.010E-4 4.740e-4 4.480E-4 si2 -1 0.97 sp2 0.985 0.015 sb2 0 f6:p 2.462E-3 2.334E-3 2.212E-3 2.096E-3 1.600E-3 1.515E-3 1.435E-3 1.359E-3 1.033E-3 9.770E-4 9.250E-4 8.750E-4 6.630E-4 6.270E-4 5.930E-4 5.610E-4 4.230E-4 Mã chương trình tính toán che chắn cho phòng đặt máy gia tốc: continue 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 -2.35 -3 -2.35 -4 -2.35 -5 -2.35 -6 -2.35 -7 -2.35 -8 -2.35 -9 -2.35 -10 -2.35 -11 -2.35 -12 -2.35 -13 -2.35 -14 -2.35 -15 -2.35 -16 -2.35 -17 -2.35 -18 -2.35 -19 -2.35 -20 -2.35 -21 -2.35 -22 -2.35 -23 -2.35 -24 -2.35 -25 -2.35 -26 -2.35 -27 -2.35 -28 -1 -2 -0.001293 -29 :-30 :-31 :-32 :-33 35 36 37 39 -0.001293 -35 -0.001293 -36 -0.001293 -37 -0.001293 -38 -0.001293 -39 -0.001293 -40 -0.001293 -41 -0.001293 -42 -0.001293 -43 -0.001293 -44 -0.001293 -34 38 40 41 42 43 44 34 box box box box box 0 1005 0 1191.5 -500 -500 -200 1505 0 0 1005 0 75 0 75 80 0 165 0 80 75 810 0 165 0 0 0 0 553 2191.5 0 661 553 553 30 105 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 mode m1 m2 m3 imp:p box box box box box box box box box box box box box box box box box box box box box box box box box box box box s s s s s s s s s s s p 6000.04p 7000.04p 1000.04p 8000.04p 13000.04p 19000.04p 8000.04p 11000.04p 14000.04p 22000.04p 80 75 461 810 0 165 0 92 890 75 115 0 165 0 553 240 160 0 379 0 553 160 240 630 0 379 0 30 160 240 335 630 0 379 0 36.4 $Slip tran 160 240 371.4 630 0 379 0 36.4 $Slip tran 160 240 407.8 630 0 379 0 36.4 $Slip tran 160 240 444.2 630 0 379 0 36.4 $Slip tran 160 240 480.6 630 0 379 0 36.4 $Slip tran 160 240 517 630 0 379 0 36 $Slip tran 790 240 215 0 379 0 553 619 80 0 213 0 553 80 619 810 0 213 0 30 80 619 461 810 0 213 0 92 890 619 115 0 213 0 553 832 80 0 61 0 553 80 832 352 0 61 0 30 80 832 461 352 0 61 0 92 432 832 573 0 61 0 553 893 80 0 233.5 0 553 80 893 925 0 233.5 0 30 80 893 461 925 0 233.5 0 92 1126.5 1005 0 65 0 553 80 75 30 810 0 165 0 431 160 240 30 630 0 379 0 305 80 619 30 810 0 213 0 431 80 832 30 352 0 61 0 431 80 893 30 925 0 233.5 0 431 550 429.5 182.5 3000 95 1009.75 182.5 15 352 1009.75 182.5 15 700 1009.75 182.5 15 1020 1009.75 182.5 15 500 1009.75 182.5 15 475 429.5 568 15 475 529.5 568 15 475 579.5 568 15 1020 700 182.5 15 1020 440 182.5 15 -0.000124 -0.755268 -0.005558 -0.4980765 -0.045746 -0.019239 -0.5134 -0.006 -0.2764 -0.0045 7r 16 $Air 8000.04p $Concrete 11000.04p 14000.04p 20000.04p $Soil 12000.04p 19000.04p 25000.04p 32 phys:p nps 1e8 c ctme 420 sdef pos=475 429.5 82.5 erg=d1 si1 L 0.25 0.35 0.45 0.55 0.65 1.05 1.15 1.25 1.35 1.45 1.85 1.95 2.05 2.15 2.25 2.65 2.75 2.85 2.95 3.05 3.45 3.55 3.65 3.75 3.85 4.25 4.35 4.45 4.55 4.65 5.05 5.15 5.25 5.35 5.45 5.85 5.95 6.05 6.15 6.25 6.65 6.75 6.85 6.95 7.05 7.45 7.55 7.65 7.75 7.85 8.25 8.35 8.45 8.55 8.65 vec=0 0.75 0.85 1.55 1.65 2.35 2.45 3.15 3.25 3.95 4.05 4.75 4.85 5.55 5.65 6.35 6.45 7.15 7.25 7.95 8.05 8.75 8.85 -0.231781 18000.04p -0.012827 -0.017101 12000.04p -0.3150923 16000.04p -0.0829411 26000.04p -0.002565 -0.001283 -0.012398 -0.013 13000.04p -0.014 20000.04p -0.0007 25r -0.067 -0.05 dir=d2 par=2 $nguon 0.95 1.75 2.55 3.35 4.15 4.95 5.75 6.55 7.35 8.15 8.95 106 $ 1, 11 $ 12, 40 9.05 9.15 9.25 9.35 9.45 9.55 9.65 9.75 9.85 9.95 10.05 sp1 1.924E-2 2.230E-2 2.451E-2 2.610E-2 2.722E-2 2.861e-2 2.845E-2 2.815E-2 2.775E-2 2.726E-2 2.472e-2 2.400E-2 2.326E-2 2.251E-2 2.176E-2 1.877e-2 1.805E-2 1.734E-2 1.664E-2 1.593E-2 1.343e-2 1.284E-2 1.227E-2 1.173E-2 1.120E-2 9.277e-3 8.843E-3 8.426E-3 8.026E-3 7.643E-3 6.267e-3 5.960E-3 5.667E-3 5.387E-3 5.119E-3 4.167e-3 3.956E-3 3.756E-3 3.564E-3 3.382E-3 2.738e-3 2.597E-3 2.462E-3 2.334E-3 2.212E-3 1.783e-3 1.689E-3 1.600E-3 1.515E-3 1.435E-3 1.153e-3 1.091E-3 1.033E-3 9.770E-4 9.250E-4 7.410e-4 7.010E-4 6.630E-4 6.270E-4 5.930E-4 4.740e-4 4.480E-4 4.230E-4 si2 -1 0.97 sp2 0.985 0.015 sb2 0 f6:p 29 30 31 32 33 34 35 36 37 38 2.797E-2 2.670E-2 2.101E-2 1.530E-2 1.069E-2 7.276E-3 4.864E-3 3.209E-3 2.096E-3 1.359E-3 8.750E-4 5.610E-4 2.841E-2 2.608E-2 2.025E-2 1.466E-2 1.020E-2 6.925E-3 4.621E-3 3.044E-3 1.986E-3 1.286E-3 8.280E-4 5.300E-4 2.861E-2 2.542E-2 1.951E-2 1.403E-2 9.730E-3 6.589E-3 4.388E-3 2.888E-3 1.882E-3 1.218E-3 7.830E-4 5.010E-4 107 ... Chương 3- Kết tính toán thảo luận CHƯƠNG I TỔNG QUAN PHƯƠNG PHÁP TÍNH TOÁN CHE CHẮN 1.1 CÁC ĐẠI LƯỢNG CƠ BẢN TRONG TÍNH TOÁN AN TOÀN BỨC XẠ Trong tính toán thiết kế bảo vệ an toàn xạ gamma (bức xạ. .. cho việc tính toán che chắn an toàn xạ gamma hay tia-X, luận văn nghiên cứu sở lý thuyết yếu tố ảnh hưởng đến tính toán che chắn loại nguồn xạ, vật liệu che chắn v.v nghiên cứu phần mềm mô MCNP5,... dụng phương pháp mô MCNP tính toán che chắn an toàn xạ Do hạn chế thời gian nên luận văn tập trung nghiên cứu tính toán che chắn cho xạ gamma viii B Mục tiêu luận văn Mục tiêu luận văn nghiên cứu