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 TỐN CHE CHẮN AN TỒ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 TỐN CHE CHẮN 1.1 Các đại lượng tính tốn an tồ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 tố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 TỐ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 tố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 tốn phương pháp tính 54 3.2.2.1 Bài tố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 tố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 tố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 tốn suất liều vị trí xung quanh phòng dọc theo Maze 63 3.3.3 Kết tính tố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 tố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ơ tố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ơ tố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 tố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 tố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 tốn: 55 Bảng Kết tính tốn trường hợp khơng tính đến ảnh hưởng tán xạ 57 Bảng 4: Kết tính tố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 tốn cho trường hợp chùm tia chiếu phía tường trước 65 Bảng Kết tính tốn cho trường hợp chùm tia chiếu phía tường sau 66 Bảng 10 Kết tính tốn cho trường hợp chùm tia chiếu lên trần 66 Bảng 11 Kết tính tố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 tốn cho trường hợp 67 Bảng 13 Kết tính tốn nguồn I-192, vật liệu bê tông 73 Bảng 14 Kết tính tốn nguồn I-192, vật liệu sắt 73 Bảng 15 Kết tính tốn nguồn I-192, vật liệu chì 74 Bảng 16 Kết tính tốn nguồn Cs-137, vật liệu bê tơng 75 Bảng 17 Kết tính tốn nguồn Cs-137, vật liệu sắt 76 Bảng 18 Kết tính tốn nguồn Cs-137, vật liệu chì 77 Bảng 19 Kết tính tốn nguồn Co-60, vật liệu bê tơng 78 Bảng 20 Kết tính tốn nguồn Co-60, vật liệu sắt 79 Bảng 21 Kết tính tốn nguồn Co-60, vật liệu chì 80 Bảng 22 Kết tính tố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 tố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 tốn nguồn tia X có lượng 150 keV, vật liệu chì 83 Bảng 25 Kết tính tố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 tố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 tốn nguồn tia X có lượng 200 keV, vật liệu sắt 85 Bảng 28 Kết tính tốn nguồn tia X có lượng 200 keV, vật liệu chì 86 Bảng 29 Kết tính tố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 tố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 tố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 tố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 tố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 tố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 tố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 tố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 tố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 tốn trường hợp khơng tính đến tán xạ 56 Hình 12 Hình học tính tố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 tốn che chắn 62 Hình 14 Hình học tốn tính tố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 tố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 tố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 TỐ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 TỐ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 TỐ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 tố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 ... điểm phương pháp mô Monte Carlo so với phương pháp tính tốn giải tích khác Từ vận dụng cách tối ưu hai phương pháp cho toán thiết kế che chắn nhằm đảm bảo an toàn xạ thực tế Áp dụng thử phương pháp. .. dụng phương pháp mơ MCNP tính tốn che chắn an tồn xạ Do hạn chế thời gian nên luận văn tập trung nghiên cứu tính tố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. .. trình cho việc tính tố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 tố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ơ