Luận văn thạc sĩ nghiên cứu chiếu xạ thanh long trên thiết bị gia tốc điện tử UERL 10 15s2

THƯ VIỆN BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH Nguyễn Nguyệt Lệ NGHIÊN CỨU CHIẾU XẠ THANH LONG TRÊN THIẾT BỊ GIA TỐC ĐIỆN TỬ UERL-10-15S2 Chuyên ngành: Vật lý nguyên tử, hạt nhân lượng cao Mã số: 604405 LUẬN VĂN THẠC SĨ VẬT LÝ NGƯỜI HƯỚNG DẪN KHOA HỌC: TS TRẦN VĂN HÙNG Thành phố Hồ Chí Minh - 2010 LỜI CÁM ƠN Trong trình học tập trường Đại học Sư phạm Thành phố Hồ Chí Minh, Quý Thầy Cô cung cấp cho kiến thức chuyên sâu, giúp trưởng thành học tập nghiên cứu khoa học Tôi xin gửi lời biết ơn đến tất Quý Thầy Cô tận tình giảng dạy tơi suốt thời gian học trường Tôi xin gửi lời biết ơn sâu sắc đến TS Trần Văn Hùng, Nghiên cứu viên Trung tâm nghiên cứu triển khai công nghệ xạ - Viện Năng lượng Nguyên tử Việt Nam định hình cho tơi lựa chọn đề tài tận tình hướng dẫn tơi suốt thời gian thực luận văn Đặc biệt, học Thầy phương pháp làm việc khoa học học có từ thực tiễn, từ vốn sống, từ am hiểu thấu đáo riêng Thầy mà khó có sách diễn đạt hết điều Tơi chân thành cám ơn Ths Trần Khắc Ân, Giám đốc Trung tâm nghiên cứu triển khai công nghệ xạ - Viện Năng lượng Nguyên tử Việt Nam tạo điều kiện tốt cho tơi suốt q trình thực đề tài Xin phép gửi lời cám ơn đến Quý Thầy Hội đồng Bảo vệ Luận văn Thạc sĩ đọc, đóng góp ý kiến, nhận xét đánh giá luận văn Tôi gửi lời cám ơn Cử nhân Nguyễn Anh Tuấn cử nhân Cao Văn Chung - Phòng vật lý Trung tâm nghiên cứu triển khai cơng nghệ xạ - có ý kiến đóng góp q báu dành cho tơi nhiều hỗ trợ nhiệt tình trình thực luận văn Tôi xin phép gửi lời cám ơn đến Sở Giáo dục – Đào tạo TP Hồ Chí Minh, Ban Giám Hiệu trường THPT Bùi Thị Xuân đồng nghiệp tạo nhiều điều kiện thuận lợi giúp đỡ tơi q trình học tập, thực luận văn Cuối cùng, xin khắc sâu công ơn Cha Mẹ, người thân, bạn bè ủng hộ, động viên giúp đỡ tơi suốt khóa học DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT Các ký hiệu D1 : Liều chiếu từ đầu quét D2 : Liều chiếu từ đầu quét DT : Liều chiếu tổng Các chữ viết tắt ACTL : Activation Library ADN : Acid Deoxyribonucleic APHIS : Animal Plant Health Inspection Service ARN : Axít ribonucleic EB : Electron beam EDB : Ethylene dibromide EtO : Ethylene oxi ENDF : Evaluated Nuclear Data File ENDL : Evaluated Nuclear Data Library IAEA : International Atomic Energy Agency LET : Linear energy transfer MCNP : Monte Carlo N-Particle NFI : National Fisheries Institute NFPA : National Food Processors Association EUREP : Euro-Retailer Produce Working Group- EUREP) GAP : Good Agricultural Practice R Dose : Relative Dose WHO : World Health Organization MỞ ĐẦU Chiếu xạ thực phẩm rau sử dụng xạ ion hóa, chẳng hạn chùm điện tử, tia gamma tia X để giảm ngăn cản sinh trưởng tiêu diệt vi sinh vật có hại vật phẩm Qua nhiều thập kỷ nghiên cứu khẳng định chiếu xạ có nhiều ứng dụng hữu ích, ví dụ tiêu diệt trùng hoa hạt ngũ cốc, chống nảy mầm khoai tây, hành tây, làm chậm chín loại tươi rau củ, gia tăng tính an toàn khử trùng sản phẩm thịt tươi đông lạnh, hải sản trứng sữa.v.v Lịch sử chiếu xạ thực phẩm khám phá tia X Roentgen 1895 chất phóng xạ Becquerel 1896 Theo sau khám phá có nhiều nghiên cứu ảnh hưởng xạ lên quan sinh học Đầu kỷ 20, nghiên cứu cho thấy xạ ion hóa hữu ích ứng dụng chiếu xạ thực phẩm Nguồn xạ sử dụng máy gia tốc hạt, tạo chùm điện tử tới lượng 24 MeV, vào cuối năm 40 kỷ 20, đồng vị phóng xạ nhân tạo Co60 Cs137 (phát xạ gamma) ứng dụng chiếu xạ công nghiệp cách phổ biến Tuy nhiên, thiết bị gia tốc điện tử ngày có tính chất ưu việt mà thiết bị sử dụng nguồn Co60 Cs137 Chẳng hạn khơng để lại chất thải phóng xạ (vì nguồn Co60 Cs137 hoạt độ thấp không sử dụng ứng dụng chiếu xạ phải chơn cất chất thải phóng xạ), hoạt động phát xạ ion hóa, cịn khơng sử dụng, tắt nguồn điện, chùm xạ ion hóa tắt.v.v Chính vậy, ngày giới có khoảng 200 thiết bị gia tốc điện tử hoạt động, phục vụ cho chiếu xạ khử trùng thực phẩm, dụng cụ y tế, nghiên cứu chế tạo vật liệu mới.v.v Ở Việt Nam, cơng ty Sơn Sơn (Bình Chánh) sử dụng chùm tia X từ máy gia tốc điện tử MeV để xử lý thực phẩm Chùm tia X tạo từ chùm điện tử MeV đập lên bia W Vào cuối năm (2010), Trung tâm Nghiên cứu Triển khai Công nghệ Bức xạ đưa vào vận hành thiết bị chiếu xạ dùng trực tiếp chùm điện tử 10 MeV từ máy gia tốc để xử lý hoa quả, thực phẩm đông lạnh phục vụ xuất khử trùng dụng cụ y tế Cũng cần nói thêm, tình hình có tính chất thời việc xuất trái Việt Nam là: Mới Hoa Kỳ chấp nhận cho phép nhập trái Thanh long Việt Nam Tuy nhiên, yêu cầu bắt buộc từ phía Hoa Kỳ Thanh long phải qua chiếu xạ để đảm bảo kiểm dịch côn trùng Chính vậy, tác giả lựa chọn đề tài “Nghiên cứu chiếu xạ Thanh long thiết bị gia tốc điện tử UERL- 10-15S2” nhằm góp phần giải vấn đề chiếu xạ trái Thanh Long xuất Việt Nam thiết bị gia tốc điện tử UERL- 10-15S2 vào hoạt động Mục đích đề tài xuất phát từ đặc điểm khả xuyên sâu chùm điện tử thấp tia X tia gamma nên khuôn khổ luận văn chủ yếu tập trung nghiên cứu phân bố liều trái Thanh long nhằm tiêu diệt côn trùng, ấu trùng trứng côn trùng, sâu bệnh trái Thanh long, đảm bảo u cầu kiểm dịch Chính vậy, luận văn mang ý nghĩa thực tế cao, nhằm đáp ứng xuất trái Thanh long vào thị trường Mỹ, Châu Âu nước khác, góp phần nâng cao giá trị xuất trái Thanh long Việt Nam Để đạt mục tiêu trên, luận văn có nhiệm vụ sau: - Tìm hiểu kích cỡ trái Thanh long, cách đóng gói xuất - Tính tốn phân bố liều bề mặt trái Thanh long - Tính tốn phân bố liều bên trái Thanh long - Độ bất đồng liều - Các giải pháp giảm hệ số bất đồng liều, đặc biệt biên độ bất đồng liều bề mặt trái Thanh long nhằm đáp ứng tiêu diệt côn trùng bề mặt - Đánh giá suất xử lý thiết bị Để thực luận văn, chúng tơi dùng chương trình MCNP, phần mềm vận chuyển xạ đa dựa phương pháp Monte-Carlo xây dựng phòng thí nghiệm quốc gia LosAlamos, Mỹ Đây cơng cụ tính tốn mạnh, mơ vận chuyển nơtron, photon electron, giải pháp toán vận chuyển xạ chiều dùng nhiều lĩnh vực tính tốn Vật lý hạt nhân Trong đề tài MCNP sử dụng để tính tốn phân bố liều xử lý trái Thanh long Luận văn xếp thành ba chương theo cấu trúc sau: CHƯƠNG 1: HÓA BỨC XẠ, HIỆU ỨNG CỦA BỨC XẠ LÊN CƠ QUAN SINH HỌC VÀ THÀNH PHẦN THỰC PHẨM 1.1 Tương tác hạt tích điện với vật chất 1.2 Cơ sở hóa xạ 1.3 Định nghĩa liều, đơn vị liều 1.4 Hiệu ứng xạ lên quan sinh học 1.4.1 Các quan vi sinh 1.4.2 Hiệu ứng xạ ion hóa 1.5 Hiệu ứng xạ lên thành phần thực phẩm 1.6 Ứng dụng chiếu xạ thực phẩm CHƯƠNG 2: THIẾT BỊ GIA TỐC ĐIỆN TỬ UERL 1015S VÀ CODE MCNP 2.1 Sơ sử dụng máy gia tốc điện tử chiếu xạ 2.2 Phân bố chùm tia liều bên vật chất chiếu xạ 2.3 Thiết bị gia tốc điện tử UERL-10-15S2 2.4 Code MCNP CHƯƠNG 3: TÍNH TỐN PHÂN BỐ LIỀU CHIẾU XẠ THANH LONG 3.1 Một số nhìn nhận chung 3.2 Phân bố liều bề mặt 3.3 Phân bố liều theo độ sâu 3.4 Hệ số bất đồng 3.5 Các kỹ thuật làm giảm độ bất đồng 3.6 Đánh giá suất thiết bị CHƯƠNG HÓA BỨC XẠ, HIỆU ỨNG CỦA BỨC XẠ LÊN CƠ QUAN SINH HỌC VÀ THÀNH PHẦN THỰC PHẨM Bức xạ dạng lượng phát trình vận động biến đổi vật chất Về mặt vật lý thể dạng sóng, hạt sóng hạt Mỗi dạng xạ đặc trưng dải lượng hay tương ứng với nó, dải bước sóng xác định Mối tương quan lượng E bước sóng λ xa mơ tả biểu thức: E  h  c 2 (1.1) Trong h - số Planck; c- vận tốc ánh sáng chân không Hiện nay, dạng xạ phổ biến áp dụng xạ electron, tia gamma, xạ hãm, xạ tử ngoại, chùm ion, xạ nơtron Cịn nguồn xạ thơng dụng bao gồm nguồn xạ thụ động (nguồn đồng vị phóng xạ Co60, Cs137…) nguồn xạ chủ động ( máy gia tốc, thiết bị phát chùm tia) [4] Để sâu nghiên cứu trình truyền lượng xạ cho vật chất xem xét đặc trưng sau đây: 1.1 Tương tác hạt tích điện với vật chất [5],[13] Các hạt mang điện bao gồm hạt nhẹ mang điện (electron, positron) hạt nặng mang điện (proton, hạt  , mảnh vỡ hạt nhân ) Các hạt mang điện vào môi trường chúng tham gia loại tương tác khác với nguyên tử toàn bộ, với electron riêng lẻ vỏ nguyên tử hay với hạt nhân nguyên tử truyền lượng cho nguyên tử, cho electron hay cho hạt nhân Các tương tác là: - Tán xạ đàn hồi với nguyên tử tồn bộ, electron khơng lượng, đổi hướng chuyển động - Va chạm không đàn hồi với electron vỏ nguyên tử: + Kích thích: xạ lượng + Ion hóa: xạ lượng - Tán xạ nhân + Tán xạ đàn hồi: xạ không lượng, đổi hướng bay + Tán xạ không đàn hồi: xạ lượng phát xạ hãm - Gây phản ứng hạt nhân Tùy theo loại xạ (nặng hay nhẹ) tùy thuộc vào lượng xạ vào môi trường mà xác xuất xảy trình lớn hay bé Các tương tác hạt mang điện với nguyên tử môi trường minh họa hình 1.1 Hình 1.1 Các tương tác hạt mang điện với nguyên tử môi trường Đối với ứng dụng chiếu xạ thực phẩm rau quả, xạ quan tâm chủ yếu điện tử Phạm vi nghiên cứu luận văn dùng trực tiếp chùm điện tử 10 MeV từ thiết bị gia tốc để xử lý hoa quả, cụ thể trái long, chúng tơi xin sâu vào khảo sát tương tác điện tử với vật chất Khi điện tử nhanh (điện tử lượng cao) qua vật chất, chúng bị va chạm Coulomb với điện tử hạt nhân nguyên tử, va chạm gây nhiều kiểu mát lượng góc lệch điện tử Đối với điện tử tới có lượng khoảng -10 MeV, khối lượng tương đối tính lớn khối lượng điện tử nguyên tử, nhỏ so với khối lượng hạt nhân nguyên tử Vì vậy, xét túy động học, va chạm tán xạ với hạt nhân nặng không dẫn đến truyền lượng, làm thay đổi đáng kể đường điện tử (sự hãm điện tử cách đột ngột) làm phát sinh tia X Ngược lại, va chạm với điện tử nguyên tử dẫn đến truyền lượng đáng kể, bứt điện tử khỏi quỹ đạo chúng ( ion hóa) Đến lượt mình, điện tử thứ cấp nhanh bị va chạm Coulomb, sinh điện tử tam cấp, tứ cấp… động điện tử tới bị hấp thụ gần hết Xét phương diện biến thiên lượng theo chiều sâu tương tác điện tử với vật chất xảy q trình vật lý sau đây: - Tán xạ khơng đàn hồi điện tử nguyên tử - Phát xạ điện tử thứ cấp - Sự va chạm (tán xạ) đàn hồi với hạt nhân làm mở rộng chùm tia theo chiều ngang - Phát xạ hãm( tia X) chủ yếu va chạm hạt nhân 1.1.1 Tán xạ không đàn hồi điện tử nguyên tử Sơ đồ tương tác hình 1.2 cho thấy điện tử nhanh tới với vận tốc v làm chuyển động điện tử nguyên tử Khoảng cách gần hạt gọi thơng số va chạm, kí hiệu b Hình 1.2 Sơ đồ va chạm điên tử tới điện tử nguyên tử Độ biến đổi vi phân momen dp chuyển đến điện tử nguyên tử thời gian dt lực điện trường F hạt xác định theo phương trình: F dp e  (đơn vị cgs) dt r (1.2) r : khoảng cách tức thời điện tử Nếu giả thuyết lượng truyền E , phần tương đối nhỏ lượng điện tử tới E hướng điện tử tới thay đổi không đáng kể va chạm, tổng biến thiên momen p nguyên tử điện tử thu cách lấy tích phân theo thời gian thành phần vng góc Vì dt  dx nên p xác định theo công thức: v p  ( e2 dx 2e2 )  cos  ( )  v r bv (1.3) Điện tử tới bị lượng E Mà E xác định công thức: E  (p)2 2m (1.4) với m khối lượng nghỉ điện tử Thay (1.3) vào (1.4) ta được: E  2e mv 2b2 (1.5) Nếu gọi N mật độ nguyên tử có nguyên tử số Z , có (2 bdb) NZ s điện tử nguyên tử gia số độ dài s có thơng số va chạm nằm b b  db biễu diễn hình 1.3 Do đó, mát lượng va chạm trung bình đơn vị chiều dài điện tử tới ( E )coll , s tính cách lấy tích phân tất giá trị thông số va chạm, công thức: (E / s )coll  4 NZ (e4 / mv )  (db / b)  4 NZ (e4 / mv )ln(bmax / bmin ) (1.6) Vì hạt tới hạt bia điện tử, nên lượng truyền cực đại E / ; thơng số va chạm tối thiểu tính từ phương trình (1.5) Giá trị cực đại thơng số va chạm cần phải kể đến thực tế điện tử ngun tử khơng hồn tồn tự mà bị liên kết trạng thái nguyên tử khác Hình 1.3 Một điện tử tới truyền lượng E cho điện tử thể tích hình trụ ( (2 bdb / s ) Dựa mơ hình thống kê ngun tử xác định lượng kích thích ion hóa trung bình Iav dạng vi phân ( dE e4 E )coll  4 NZ ( )ln( ) (1.7) ds mv I av I av xấp xỉ với nguyên tử số vật liệu, với số tỷ lệ bậc với lượng Rydberg 13.5 eV Chú ý lượng điện tử phi tương đối tính, E  mv cơng suất hãm va chạm giảm xấp xỉ theo E 1 Tuy nhiên, lượng tương đối tính vận tốc v gần tốc độ ánh sáng công suất hãm tăng logaric với lượng 1.1.2 Phát xạ điện tử thứ cấp Phổ điện tử thứ cấp, q trình va chạm sinh ra, thu ý xác suất tương tác vi phân d cần phải tỷ lệ với 2 bdb Việc lấy vi phân phương trình (1.5) cho thấy phổ lượng biến thiên biểu thức Hình 3.22 Kết tính tốn phân bố liều che chắn hai lớp nhơm tồn trái long 3.6 Đánh giá suất thiết bị.[4] Máy gia tốc UERL 10-15S2 thuộc nhóm máy gia tốc lượng cao Máy có ưu điểm sau: - Suất liều lớn: ưu điểm giúp thời gian xử lý nhanh, cho sản lượng cao, có giá thành giảm tiết kiệm lượng - Tác động theo hướng định: Nếu nguồn gamma phát xạ theo hướng, kể hướng khơng có sản phẩm cần chiếu xạ, chùm hạt gia tốc ln hướng theo phía có sản phẩm Do hiệu suất sử dụng lượng tăng đáng kể so với nguồn gamma - Hiệu suất sử dụng lượng cao Theo định nghĩa Năng lượng hấp thụ sản phẩm Hiệu suất sử dụng lượng  Năng lượng nguồn phát Hiệu suất sử dụng lượng để xử lý thực phẩm giới thiệu bảng 3.12 Bảng 3.12 So sánh hiệu suất sử dụng lượng nguồn xạ [4] Nguồn xạ Hiệu suất % Máy gia tốc electron nhanh 10 MeV 60 Nguồn xạ hãm Emax=5 MeV 50 Nguồn 60Co 30 Nguồn 137Cs 20 Từ bảng ta thấy hiệu suất sử dụng lượng nguồn xạ electron nhanh lớn Bên cạnh ưu điểm nêu trên, máy gia tốc UERL 10-15S2 nhược điểm sau: - Độ xuyên thấp: Nhược điểm chủ yếu xạ electron quan điểm công nghệ xạ độ xuyên thấp so với xạ gamma - Tính khơng đồng liều: Khả xun sâu thấp electron nhanh cịn gây tính khơng đồng liều vật bị chiếu KẾT LUẬN Việc ứng dụng máy gia tốc chùm tia điện tử để chiếu xạ trái nói chung trái Thanh long nói riêng khó khăn khả xuyên sâu chùm điện tử thấp Trong năm gần việc nghiên cứu sử dụng code MCNP q trình mơ tính tốn phân bố liều chiếu xạ để xử lý hoa quả, thực phẩm đông lạnh phục vụ xuất khử trùng dụng cụ y tế trở nên phổ biến ưu điểm chương trình mang lại Trong luận văn này, nghiên cứu chiếu xạ long thiết bị gia tốc điện tử UERL10-15S2, thu kết sau: - Tìm hiểu kích cỡ trái Thanh long, cách đóng gói xuất - Tính tốn phân bố liều bề mặt trái Thanh long - Tính tốn phân bố liều bên trái Thanh long - Độ bất đồng liều - Các giải pháp giảm hệ số bất đồng liều, đặc biệt biên độ bất đồng liều bề mặt trái Thanh long nhằm đáp ứng tiêu diệt côn trùng bề mặt - Đánh giá suất xử lý thiết bị Qua tính tốn phân bố liều ta thấy chiếu trực tiếp chùm tia điện tử lượng 10 MeV xuống trái long độ bất đồng 2.44 Để khai thác máy gia tốc chùm tia điện tử lượng 10 MeV vào việc chiếu xạ trái long ta cần phải tính tốn phương pháp giảm độ bất đồng liều Các tính tốn luận văn cho thấy che chắn hai lớp nhôm (chiều dày lớp 0.2 cm) cho độ bất đồng liều tốt (1.67) Đây xem phương pháp hữu hiệu để khai thác máy gia tốc chùm tia điện tử chiếu xạ long phát đề tài Với kết nhận ban đầu đề tài xem tài liệu tham khảo hữu ích cho nhân viên vận hành trình tìm hiểu máy gia tốc xác định liều cho trái Thanh long Đồng thời sở khuyến cáo khách hàng nên đóng gói kích thước thùng hàng để xử lý hiệu Qua đó, đề tài giải nhiệm vụ trung tâm nghiên cứu triển khai cơng nghệ xạ dựa sở tính tốn đề tài, trung tâm đưa vào kiểm chứng thực tiễn máy gia tốc lắp đặt HƯỚNG PHÁT TRIỂN Đề tài tính tốn mơ chương trình MCNP, chưa đo phân bố liều thực tế trung tâm “Nghiên cứu triển khai công nghệ xạ” trình xây dựng lắp đặt máy gia tốc Nếu có thêm số liệu đo đạc số liệu tính tốn đề tài khẳng định tính đắn phương pháp tính tốn có giá trị nhiều Với hạn chế nêu, tác giả xin đề nghị luận văn sau tiếp tục nghiên cứu tìm hiểu sâu thêm đề tài tiến hành lấy số liệu đo đạc từ thực nghiệm sau trung tâm hoàn thành việc lắp đặt máy gia tốc Từ kiểm chứng đối chiếu với số liệu tính tốn để ngày hồn chĩnh luận văn hơn, góp phần đưa chiếu xạ Thanh long cách sử dụng máy gia tốc ứng dụng rộng rãi Bình Thuận đáp ứng nhu cầu xuất Đây công việc nên tiếp tục tương lai TÀI LIỆU THAM KHẢO Tiếng Việt Trần Văn Hùng (2008), Báo cáo tổng kết tính tốn liều chiếu xạ chế độ chiếu xạ phục vụ xử lý hàng đông lạnh dụng cụ y tế máy gia tốc UERL-10-15T, TP Hồ Chí Minh Ngơ Quang Huy (2006), Cơ sở vật lý hạt nhân, Nxb khoa học kỹ thuật, Hà Nội Nguyễn Quang Miên (2008), Giáo trình ghi nhận đo lường xạ hạt nhân, trường Đại học Sư phạm TP Hồ Chí Minh Trần Đại Nghiệp (2002), Giáo trình cơng nghệ xạ, Nxb khoa học kỹ thuật, Hà Nội Nguyễn Đơng Sơn (2009), Giáo trình ứng dụng xạ ion hóa kỹ thuật hạt nhân y tế,trường Đại học Sư phạm TP Hồ Chí Minh Châu Văn Tạo (2004), An tồn xạ ion hóa, Nxb Đại học quốc gia TP Hồ Chí Minh, TP Hồ Chí Minh http://www.xaluan.com/modules.php?name=News&file=article&sid=113150 Tiếng Anh Corad Service Ltd, office 206, 27, Partizanskaya Street, St, Petersburg, 195248, Russia – Design and Manufacture in X-Ray and accelerators technology (2009), Technical requirements of the electron beam system delivered under contract No 01/12-08-1 and rooms for this electron beam system Council for Agricultural Sience and Technology,(1989), Tạp chí Ionizing enery in food processing and pest control, Task Force Report, No 115, ii.applications 10 J.F.Briesmeister, (1997), MCNP- A General Monte Carlo N-Particle Transport Code Version 4C2, Transport Methods Group, Los Alamos National Laboratory 11 Loiseau (2007), IAEA / RCA Regional training course e-beam & X-ray applications fruits & frozen food, High – performance X-ray sustems for radiation processing 12 NCRP report No.51 Radiation protection design guidelines for 0.1 – 100 MeV paricle accelerator facilities 13 R.B.Miller (2005) Electronic Irradiation of Food, An Introduction to the Technology, Springer, New Mexico PHỤ LỤC Phụ lục A: Input tính tốn phân bố liều toàn trái long -0.97 -46 #21 #22 #23 #24 #25 #26 #27 #28 #29 & #30 #31 #32 #33 #34 #35 #36 #37 #38 #39 #40 #41 #42 & #43 #44 #45 #46 #47 #48 #49 #50 #51 #52 #53 #54 #55 & #56 #57 #58 #59 #60 #61 #62 #63 #64 #65 #66 #67 #68 & #69 #70 #71 #72 #73 trcl=(0 -125) imp:e=2 $trai long 10 -0.001293 -6 23 24 -25 26 -27 #1 #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 #45 #46 #47 #48 #49 #50 #51 #52 #53 #54 #55 & #56 #57 #58 #59 #60 #61 #62 #63 #64 #65 #66 #67 #68 & #69 #70 #71 #72 #73 imp:e=1 $kk 20 10 -4.5 -6 13 -8 11 -12 imp:e=1 $cua so tital 21 -0.97 -39 trcl=(0 -120.5) imp:e=2 22 like 21 but trcl=(0 -121.5) 23 like 21 but trcl=(0 -122.5) 24 like 21 but trcl=(0 -123.75) 25 like 21 but trcl=(0 -125) 26 like 21 but trcl=(0 -126.25) 27 like 21 but trcl=(0 -127.5) 28 like 21 but trcl=(0 -128.5) 29 like 21 but trcl=(0 -129.5) 30 like 21 but trcl=(1 -120.6) 31 like 21 but trcl=(1 -121.6) 32 like 21 but trcl=(1 -122.6) 33 like 21 but trcl=(1 -123.76) 34 like 21 but trcl=(1 -125) 35 like 21 but trcl=(1 -126.24) 36 like 21 but trcl=(1 -127.4) 37 like 21 but trcl=(1 -128.4) 38 like 21 but trcl=(1 -129.4) 39 like 21 but trcl=(2 -120.8) 40 like 21 but trcl=(2 -121.8) 41 like 21 but trcl=(2 -122.8) 42 like 21 but trcl=(2 -123.8) 43 like 21 but trcl=(2 -125) 44 like 21 but trcl=(2 -126.1) 45 like 21 but trcl=(2 -127.1) 46 like 21 but trcl=(2 -128.1) 47 like 21 but trcl=(2 -129.18) 48 like 21 but trcl=(3 -121.2) 49 like 21 but trcl=(3 -122) 50 like 21 but trcl=(3 -123) 51 like 21 but trcl=(3 -124) 52 like 21 but trcl=(3 -125) 53 like 21 but trcl=(3 -126) 54 like 21 but trcl=(3 -127) 55 like 21 but trcl=(3 -128) 56 like 21 but trcl=(3 -128.8) 57 like 21 but trcl=(4 -121.8) 58 like 21 but trcl=(4 -122.59) 59 like 21 but trcl=(4 -123.38) 60 like 21 but trcl=(4 -124.2) 61 like 21 but trcl=(4 -125) 62 like 21 but trcl=(4 -125.8) 63 like 21 but trcl=(4 -126.6) 64 like 21 but trcl=(4 -127.36) 65 like 21 but trcl=(4 -128.1) 66 like 21 but trcl=(5 -122.9) 67 like 21 but trcl=(5 -123.65) 68 like 21 but trcl=(5 -124.35) 69 like 21 but trcl=(5 -125) 70 like 21 but trcl=(5 -125.67) 71 like 21 but trcl=(5 -126.35) 72 like 21 but trcl=(5 -127.1) 73 like 21 but trcl=(5.75 -125) c 31 -2.7 40 -41 42 -43 44 -45 trcl=(0 -120.9) imp:e=1 $lop nhom c 32 -2.7 47 -48 49 -50 51 -52 trcl=(5.75 -120.3) imp:e=1 $lop nhom 100 (6:-23:-24:25:-26:27) imp:e=0 p -0.61875 -1.5 p -0.61875 -0.5 p 0.61875 0.5 p 0.61875 1.5 pz pz -80 px -10 px 10 px -1.05 10 px 1.05 11 py -40 12 py 40 13 pz -80.0135 14 px -20 15 px 20 16 py -23 17 py 23 18 pz -110.0135 19 px -60 20 px 60 21 pz -150.0235 22 pz -130.0135 23 pz -140 24 px -25 25 px 25 26 py -50 27 py 50 28 px -220 29 px 220 30 py -220 31 py 220 32 pz -250 33 px -0.25 34 px 0.25 35 py -0.25 36 py 0.25 37 pz -0.25 38 pz 0.25 39 so 0.25 40 px -10 41 px 10 42 py -40 43 py 40 44 pz -0.3 45 pz 0.3 46 ell -2.34 0 2.34 0 47 px -2 48 px 49 py -2 50 py 51 pz -0.3 52 pz 0.3 mode e m2 7000 -0.755 8000 -0.232 & 18000 -0.013$Khong m6 1000 -0.0069 6000 -0.04542 7000 -0.00383 8000 -0.004385$Dummy m5 13000 -1.0$ Nhom m10 22000 $tital sdef erg=10 pos=0 -80.0001 X=d2 Y=d1 Z=-80.0001 par=3 vec=0 -1 dir=1 SI1 -30 30 SP1 si2 -10 10 sp2 f4:e 21 22 23 24 25 26 27 28 29 f14:e 30 31 32 33 34 35 36 37 38 f24:e 39 40 41 42 43 44 45 46 47 f34:e 48 49 50 51 52 53 54 55 56 f44:e 57 58 59 60 61 62 63 64 65 f54:e 66 67 68 69 70 71 72 73 DE 0.01 0.03 0.05 0.07 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 & 0.55 0.6 0.65 0.7 0.8 1.4 1.8 2.3 2.5 2.8 3.3 3.8 4.3 4.8 & 5.3 5.8 6.3 6.8 7.3 7.8 8.3 8.8 9.3 9.8 10.3 DF 6.22E+09 1.99E+09 1.26E+09 9.64E+08 7.41E+08 5.66E+08 4.76E+08 & 4.21E+08 3.83E+08 3.56E+08 3.35E+08 3.19E+08 3.05E+08 2.94E+08 & 2.85E+08 2.77E+08 2.70E+08 2.59E+08 2.43E+08 2.24E+08 2.15E+08 & 2.11E+08 2.08E+08 2.06E+08 2.03E+08 2.01E+08 1.99E+08 1.98E+08 & 1.97E+08 1.96E+08 1.95E+08 1.95E+08 1.94E+08 1.94E+08 1.94E+08 & 1.93E+08 1.93E+08 1.93E+08 1.93E+08 1.92E+08 ctme 1440 Phụ lục B: Input tính tốn phân bố liều toàn trái long che hai lớp nhôm -0.97 -46 #21 #22 #23 #24 #25 #26 #27 #28 #29 & #30 #31 #32 #33 #34 #35 #36 #37 #38 #39 #40 #41 #42 & #43 #44 #45 #46 #47 #48 #49 #50 #51 #52 #53 #54 #55 & #56 #57 #58 #59 #60 #61 #62 #63 #64 #65 #66 #67 #68 & #69 #70 #71 #72 #73 trcl=(0 -125) imp:e=2 $hang 10 -0.001293 -6 23 24 -25 26 -27 #1 #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 #45 #46 #47 #48 #49 #50 #51 #52 #53 #54 #55 & #56 #57 #58 #59 #60 #61 #62 #63 #64 #65 #66 #67 #68 & #69 #70 #71 #72 #73 #74 #75 imp:e=1 $kk 20 10 -4.5 -6 13 -8 11 -12 imp:e=1 $cua so tital 21 -0.97 -39 trcl=(0 -120.5) imp:e=2 22 like 21 but trcl=(0 -121.5) 23 like 21 but trcl=(0 -122.5) 24 like 21 but trcl=(0 -123.75) 25 like 21 but trcl=(0 -125) 26 like 21 but trcl=(0 -126.25) 27 like 21 but trcl=(0 -127.5) 28 like 21 but trcl=(0 -128.5) 29 like 21 but trcl=(0 -129.5) 30 like 21 but trcl=(1 -120.6) 31 like 21 but trcl=(1 -121.6) 32 like 21 but trcl=(1 -122.6) 33 like 21 but trcl=(1 -123.76) 34 like 21 but trcl=(1 -125) 35 like 21 but trcl=(1 -126.24) 36 like 21 but trcl=(1 -127.4) 37 like 21 but trcl=(1 -128.4) 38 like 21 but trcl=(1 -129.4) 39 like 21 but trcl=(2 -120.8) 40 like 21 but trcl=(2 -121.8) 41 like 21 but trcl=(2 -122.8) 42 like 21 but trcl=(2 -123.8) 43 like 21 but trcl=(2 -125) 44 like 21 but trcl=(2 -126.1) 45 like 21 but trcl=(2 -127.1) 46 like 21 but trcl=(2 -128.1) 47 like 21 but trcl=(2 -129.18) 48 like 21 but trcl=(3 -121.2) 49 like 21 but trcl=(3 -122) 50 like 21 but trcl=(3 -123) 51 like 21 but trcl=(3 -124) 52 like 21 but trcl=(3 -125) 53 like 21 but trcl=(3 -126) 54 like 21 but trcl=(3 -127) 55 like 21 but trcl=(3 -128) 56 like 21 but trcl=(3 -128.8) 57 like 21 but trcl=(4 -121.8) 58 like 21 but trcl=(4 -122.59) 59 like 21 but trcl=(4 -123.38) 60 like 21 but trcl=(4 -124.2) 61 like 21 but trcl=(4 -125) 62 like 21 but trcl=(4 -125.8) 63 like 21 but trcl=(4 -126.6) 64 like 21 but trcl=(4 -127.36) 65 like 21 but trcl=(4 -128.1) 66 like 21 but trcl=(5 -122.9) 67 like 21 but trcl=(5 -123.65) 68 like 21 but trcl=(5 -124.35) 69 like 21 but trcl=(5 -125) 70 like 21 but trcl=(5 -125.67) 71 like 21 but trcl=(5 -126.35) 72 like 21 but trcl=(5 -127.1) 73 like 21 but trcl=(5.75 -125) 74 -2.7 40 -41 42 -43 44 -45 trcl=(0 -120.1) imp:e=1 $lop nhom 75 -2.7 47 -48 49 -50 51 -52 trcl=(5.75 -119.85) imp:e=1 $lop nhom 100 (6:-23:-24:25:-26:27) imp:e=0 p -0.61875 -1.5 p -0.61875 -0.5 p 0.61875 0.5 p 0.61875 1.5 pz pz -80 px -10 px 10 px -1.05 10 px 1.05 11 py -40 12 py 40 13 pz -80.0135 14 px -20 15 px 20 16 py -23 17 py 23 18 pz -110.0135 19 px -60 20 px 60 21 pz -150.0235 22 pz -130.0135 23 pz -140 24 px -25 25 px 25 26 py -50 27 py 50 28 px -220 29 px 220 30 py -220 31 py 220 32 pz -250 33 px -0.25 34 px 0.25 35 py -0.25 36 py 0.25 37 pz -0.25 38 pz 0.25 39 so 0.25 40 px -10 41 px 10 42 py -40 43 py 40 44 pz -0.01 45 pz 0.01 46 ell -2.34 0 2.34 0 47 px -2.5 48 px 2.5 49 py -2 50 py 51 pz -0.2 52 pz 0.2 mode e m2 7000 -0.755 8000 -0.232 & 18000 -0.013$Khong m6 1000 -0.0069 6000 -0.04542 7000 -0.00383 8000 -0.004385$Dummy m5 13000 -1.0$ Nhom m10 22000 $tital sdef erg=10 pos=0 -80.0001 X=d2 Y=d1 Z=-80.0001 par=3 vec=0 -1 dir=1 SI1 -30 30 SP1 si2 -10 10 sp2 f4:e 21 22 23 24 25 26 27 28 29 f14:e 30 31 32 33 34 35 36 37 38 f24:e 39 40 41 42 43 44 45 46 47 f34:e 48 49 50 51 52 53 54 55 56 f44:e 57 58 59 60 61 62 63 64 65 f54:e 66 67 68 69 70 71 72 73 DE 0.01 0.03 0.05 0.07 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 & 0.55 0.6 0.65 0.7 0.8 1.4 1.8 2.3 2.5 2.8 3.3 3.8 4.3 4.8 & 5.3 5.8 6.3 6.8 7.3 7.8 8.3 8.8 9.3 9.8 10.3 DF 6.22E+09 1.99E+09 1.26E+09 9.64E+08 7.41E+08 5.66E+08 4.76E+08 & 4.21E+08 3.83E+08 3.56E+08 3.35E+08 3.19E+08 3.05E+08 2.94E+08 & 2.85E+08 2.77E+08 2.70E+08 2.59E+08 2.43E+08 2.24E+08 2.15E+08 & 2.11E+08 2.08E+08 2.06E+08 2.03E+08 2.01E+08 1.99E+08 1.98E+08 & 1.97E+08 1.96E+08 1.95E+08 1.95E+08 1.94E+08 1.94E+08 1.94E+08 & 1.93E+08 1.93E+08 1.93E+08 1.93E+08 1.92E+08 ctme 500 ... CHƯƠNG 2: THIẾT BỊ GIA TỐC ĐIỆN TỬ UERL 101 5S VÀ CODE MCNP 2.1 Sơ sử dụng máy gia tốc điện tử chiếu xạ 2.2 Phân bố chùm tia liều bên vật chất chiếu xạ 2.3 Thiết bị gia tốc điện tử UERL- 10- 15S2 2.4... cứu chiếu xạ Thanh long thiết bị gia tốc điện tử UERL- 10- 15S2? ?? nhằm góp phần giải vấn đề chiếu xạ trái Thanh Long xuất Việt Nam thiết bị gia tốc điện tử UERL- 10- 15S2 vào hoạt động Mục đích đề... dopt-vật liệu = dopt-nước / 2.3 Thiết bị gia tốc điện tử UERL 10- 15S2 [8] 2.3.1 Giới thiệu máy gia tốc UERL- 10- 15S2 Máy gia tốc electron tuyến tính UERL- 10- 15S2 lắp đặt Vinagamma thuộc loại