Luận Văn Thạc Sĩ: Mô Phỏng Bài Toán Vận Chuyển Neutron

Luận văn thạc sĩ toán học nghiên cứu hus mô phỏng bài toán vận chuyển neutron 03, khảo sát thực trạng, phân tích nguyên nhân, đề xuất giải pháp cải thiện thực tiễn.

• Tác giả: Vũ Thị Linh
• Người hướng dẫn: Ts. Nguyễn Hoàng Oanh
• Trường học: Trường Đại Học Khoa Học Tự Nhiên
• Chuyên ngành: Vật lí lí thuyết và vật lí toán
• Đề tài: Mô Phỏng Bài Toán Vận Chuyển Neutron
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2015
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CẢM ƠN
• MỤC LỤC
• 1. CHƯƠNG I: GIỚI THIỆU VỀ CƠ BẢN VỀ BÀI TOÁN VẬN CHUYỂN NEUTRON
• 1.1. Tương tác của neutron với vật chất
• 2. CHƯƠNG II: GIỚI THIỆU VỀ PHƯƠNG PHÁP MONTE CARLO
• 3. CHƯƠNG III: MÔ PHỎNG BÀI TOÁN VẬN CHUYỂN NEUTRON
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Về Mô Phỏng Vận Chuyển Neutron Trong Lò Phản Ứng Hạt Nhân

Mô phỏng vận chuyển neutron trong lò phản ứng hạt nhân là một lĩnh vực nghiên cứu quan trọng trong vật lý hạt nhân. Nó giúp hiểu rõ hơn về cách thức neutron tương tác với vật chất và ảnh hưởng của chúng đến quá trình phân hạch. Việc mô phỏng này không chỉ giúp tối ưu hóa thiết kế lò phản ứng mà còn nâng cao hiệu quả an toàn trong vận hành. Các phương pháp mô phỏng hiện nay chủ yếu bao gồm phương pháp Monte Carlo và phương pháp tất định, mỗi phương pháp có những ưu điểm và nhược điểm riêng.

1.1. Khái Niệm Cơ Bản Về Vận Chuyển Neutron

Vận chuyển neutron là quá trình mà neutron di chuyển qua các vật liệu và tương tác với chúng. Các phản ứng này có thể bao gồm hấp thụ, phân hạch và tán xạ. Mỗi loại phản ứng có những đặc điểm riêng và ảnh hưởng đến sự phân bố neutron trong lò phản ứng.

1.2. Tầm Quan Trọng Của Mô Phỏng Neutron

Mô phỏng neutron giúp các nhà nghiên cứu dự đoán hành vi của neutron trong lò phản ứng, từ đó cải thiện thiết kế và tối ưu hóa hiệu suất. Điều này đặc biệt quan trọng trong việc đảm bảo an toàn và hiệu quả trong các ứng dụng năng lượng hạt nhân.

II. Vấn Đề Và Thách Thức Trong Mô Phỏng Vận Chuyển Neutron

Mô phỏng vận chuyển neutron đối mặt với nhiều thách thức, bao gồm độ phức tạp của các phương trình mô tả tương tác neutron và yêu cầu tính toán cao. Phương trình Boltzmann, mặc dù mạnh mẽ, nhưng rất khó giải quyết trong các bài toán thực tế. Hơn nữa, việc lựa chọn mô hình phù hợp để mô phỏng cũng là một thách thức lớn.

2.1. Khó Khăn Trong Việc Giải Phương Trình Boltzmann

Phương trình Boltzmann mô tả sự vận chuyển của neutron trong không gian và thời gian, nhưng việc giải nó trong các điều kiện thực tế là rất phức tạp. Cần phải sử dụng các phương pháp số để xấp xỉ nghiệm, điều này có thể dẫn đến sai số.

2.2. Tính Toán Tốn Thời Gian Và Tài Nguyên

Mô phỏng neutron yêu cầu tính toán rất lớn, đặc biệt là khi sử dụng phương pháp Monte Carlo. Điều này có thể dẫn đến thời gian xử lý dài và yêu cầu phần cứng mạnh mẽ, gây khó khăn cho nhiều nhà nghiên cứu.

III. Phương Pháp Monte Carlo Trong Mô Phỏng Vận Chuyển Neutron

Phương pháp Monte Carlo là một trong những kỹ thuật phổ biến nhất trong mô phỏng vận chuyển neutron. Nó sử dụng các mẫu ngẫu nhiên để mô phỏng quá trình tương tác của neutron với vật chất. Phương pháp này cho phép mô phỏng các tình huống phức tạp mà các phương pháp khác không thể thực hiện được.

3.1. Nguyên Tắc Cơ Bản Của Phương Pháp Monte Carlo

Phương pháp Monte Carlo dựa trên việc tạo ra các mẫu ngẫu nhiên để mô phỏng các quá trình vật lý. Mỗi mẫu đại diện cho một kịch bản có thể xảy ra trong thực tế, từ đó cho phép tính toán xác suất và phân bố của các neutron.

3.2. Ưu Điểm Của Phương Pháp Monte Carlo

Phương pháp Monte Carlo có nhiều ưu điểm, bao gồm khả năng mô phỏng các tình huống phức tạp và tính chính xác cao. Nó cũng cho phép xử lý các biến đổi không gian và năng lượng một cách linh hoạt, điều này rất quan trọng trong nghiên cứu lò phản ứng hạt nhân.

IV. Ứng Dụng Thực Tiễn Của Mô Phỏng Vận Chuyển Neutron

Mô phỏng vận chuyển neutron có nhiều ứng dụng thực tiễn trong ngành công nghiệp hạt nhân. Nó được sử dụng để thiết kế lò phản ứng, tối ưu hóa quy trình vận hành và đảm bảo an toàn. Ngoài ra, mô phỏng cũng giúp nghiên cứu các vật liệu mới cho lò phản ứng.

4.1. Thiết Kế Lò Phản Ứng Hạt Nhân

Mô phỏng neutron giúp các kỹ sư thiết kế lò phản ứng hiệu quả hơn bằng cách dự đoán hành vi của neutron trong các điều kiện khác nhau. Điều này giúp tối ưu hóa hiệu suất và giảm thiểu rủi ro.

4.2. Nghiên Cứu Vật Liệu Mới

Mô phỏng cũng được sử dụng để nghiên cứu các vật liệu mới có khả năng hấp thụ neutron tốt hơn, từ đó cải thiện hiệu suất của lò phản ứng và giảm thiểu chất thải hạt nhân.

V. Kết Luận Về Mô Phỏng Vận Chuyển Neutron Trong Lò Phản Ứng Hạt Nhân

Mô phỏng vận chuyển neutron là một lĩnh vực nghiên cứu quan trọng với nhiều ứng dụng thực tiễn trong ngành công nghiệp hạt nhân. Mặc dù còn nhiều thách thức, nhưng với sự phát triển của công nghệ và phương pháp mô phỏng, tương lai của lĩnh vực này hứa hẹn sẽ mang lại nhiều tiến bộ đáng kể.

5.1. Tương Lai Của Nghiên Cứu Mô Phỏng Neutron

Với sự phát triển của công nghệ máy tính và các thuật toán mới, nghiên cứu mô phỏng neutron sẽ ngày càng trở nên chính xác và hiệu quả hơn. Điều này sẽ mở ra nhiều cơ hội mới trong việc tối ưu hóa lò phản ứng và nghiên cứu vật liệu.

5.2. Tầm Quan Trọng Của An Toàn Hạt Nhân

An toàn trong ngành công nghiệp hạt nhân là một yếu tố quan trọng. Mô phỏng vận chuyển neutron sẽ tiếp tục đóng vai trò quan trọng trong việc đảm bảo an toàn và hiệu quả trong vận hành lò phản ứng hạt nhân.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật hạt nhân và vật lý lò phản ứng, việc mô phỏng chính xác sự vận chuyển của dòng neutron đóng vai trò sống còn đối với hiệu suất vận hành và an toàn bức xạ. Năng lượng sinh ra trong một phản ứng phân hạch đạt mức trung bình khoảng 200 MeV, trong đó khoảng 80% được giải phóng trực tiếp dưới dạng động năng của các mảnh phân hạch. Để kiểm soát chuỗi phản ứng dây chuyền này, việc giải bài toán vận chuyển neutron là yêu cầu cốt lõi. Tuy nhiên, phương pháp tất định truyền thống gặp phải rào cản tính toán cực lớn khi phải giải phương trình vi tích phân Boltzmann trong không gian 7 chiều, bao gồm 3 tọa độ không gian, 2 biến góc định hướng chuyển động, 1 chiều năng lượng và 1 chiều thời gian.

Nghiên cứu này được thực hiện nhằm xây dựng chương trình mô phỏng sự vận chuyển neutron qua các lớp vật liệu dựa trên phương pháp thống kê Monte Carlo, ứng dụng nền tảng phần mềm tính toán số mã nguồn mở Scilab kết hợp ngôn ngữ lập trình C. Phạm vi nghiên cứu được triển khai tại Bộ môn Vật lý lý thuyết và Vật lý toán thuộc Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội trong giai đoạn 2011–2015. Đề tài tập trung khảo sát các tương tác va chạm đàn hồi và phi đàn hồi của chùm neutron tới với bia vật liệu phẳng có bề dày 1 đơn vị độ dài.

Ý nghĩa thực tiễn của công trình thể hiện ở việc giải quyết bài toán với cỡ mẫu 10.000 lịch sử hạt cho mỗi trường hợp khảo sát, giúp tối ưu hóa chi phí bản quyền phần mềm thương mại đắt đỏ về mức 0 đồng. Đồng thời, mô hình cung cấp một công cụ giảng dạy và nghiên cứu trực quan, làm sáng tỏ các quy luật tán xạ, hấp thụ và truyền qua của neutron với độ chính xác thống kê cao.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu xây dựng trên nền tảng lý thuyết vận chuyển neutron tổng quát và động học lò phản ứng hạt nhân. Ba trụ cột lý thuyết chính được vận dụng bao gồm:

Thứ nhất là lý thuyết vận chuyển hạt Boltzmann, mô tả sự cân bằng mật độ góc của neutron trong một yếu tố thể tích vô cùng nhỏ của không gian pha 6 chiều. Phương trình vi tích phân này cân bằng giữa tốc độ biến thiên mật độ hạt, dòng hạt di chuyển vào thể tích và các quá trình thất thoát do hấp thụ hoặc tán xạ ra ngoài chùm.

Thứ hai là lý thuyết tương tác hạt nhân vi mô và vĩ mô. Xác suất neutron tương tác với hạt nhân bia được xác định thông qua tiết diện vi mô $\sigma$ và tiết diện vĩ mô $\Sigma = N \cdot \sigma$, với $N$ là mật độ hạt nhân nguyên tử. Khoảng cách tương tác của lực hạt nhân mạnh được giới hạn trong phạm vi khoảng $10^{-15}$ m (cỡ bán kính hạt nhân). Các mode tương tác chính gồm có: phản ứng hấp thụ bắt bức xạ $(n,\gamma)$, phản ứng phân hạch sinh từ 2 đến 6 neutron mới, tán xạ đàn hồi bảo toàn động năng và tán xạ phi đàn hồi kích thích hạt nhân.

Thứ ba là mô hình phân bố năng lượng và thời gian sống của neutron. Phổ năng lượng neutron phân hạch tức thời tuân theo hàm phân bố Maxwell với đỉnh cực đại quanh 1 MeV và dải năng lượng kéo dài từ 100 eV đến 10 MeV. Thời gian sống trung bình của neutron trong lò phản ứng nhiệt cỡ 0,1 ms, trong đó quá trình làm chậm thông qua các tán xạ đàn hồi với hạt nhân nhẹ chỉ mất khoảng 0,01 ms.

Phương pháp nghiên cứu

Nghiên cứu kết hợp dữ liệu tiết diện hạt nhân thực nghiệm từ cơ sở dữ liệu quốc tế ENDF/B-VI ở nhiệt độ 300K của các đồng vị như Urani-235, Urani-238, Plutoni-239 và Americi-241. Để giải bài toán vận chuyển đa chiều phức tạp mà không cần giải tích phương trình vi phân, tác giả lựa chọn phương pháp thống kê Monte Carlo.

Cỡ mẫu mô phỏng được thiết lập chuẩn hóa ở mức 10.000 hạt neutron cho mỗi kịch bản tính toán. Phương pháp chọn mẫu dựa trên nguyên lý chuỗi Markov và kỹ thuật lấy mẫu giả ngẫu nhiên đồng dư tuyến tính. Để tạo phân bố góc cực $\cos\theta = 1 - 2r$ và phân bố quãng đường tự do $l = -\lambda \ln(r)$, thuật toán sử dụng phương pháp đổi biến xác suất và phương pháp loại trừ Von Neumann.

Lý do lựa chọn phương pháp Monte Carlo trên nền Scilab là vì sai số thống kê của Monte Carlo tỉ lệ thuận với nghịch đảo căn bậc hai của số phép thử ($1/\sqrt{N} = 1/\sqrt{10000} = 1%$), hoàn toàn độc lập với số chiều của bài toán tích phân. Điều này mang lại ưu thế vượt trội so với các phương pháp giải tích hình thang hay Simpson khi số chiều tăng lên. Tiến trình nghiên cứu được thực hiện tuần tự từ tháng 9 đến tháng 12 năm 2015, bao gồm các khâu thiết lập thuật toán, lập trình mã nguồn, chạy kiểm thử và phân tích đồ thị tương tác.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Mô phỏng bài toán vận chuyển 10.000 hạt neutron qua tấm chắn phẳng có bề dày chuẩn $t = 1$ đơn vị độ dài đã mang lại các phát hiện thực nghiệm số quan trọng:

Thứ nhất, trong bài toán tán xạ đàn hồi với xác suất hấp thụ $p_c = 0.3333$ (xác suất tán xạ $p_s = 0.6667$), sự gia tăng của quãng đường tự do trung bình $\lambda$ từ 0.1 đến 10 đơn vị độ dài làm thay đổi căn bản phân bố trạng thái hạt. Khi $\lambda = 10$, tỉ lệ neutron truyền qua tấm vật liệu đạt 39,72% (tương ứng 3.972 hạt), tỉ lệ hạt bị phản xạ ngược về phía nguồn phát chiếm 55,10% (5.510 hạt), và tỉ lệ hạt bị hấp thụ giảm mạnh xuống chỉ còn 5,18% (518 hạt). Khi $\lambda \le 0.2$, tỉ lệ hạt truyền qua ghi nhận bằng 0%.

Thứ hai, khi nâng xác suất hấp thụ lên $p_c = 0.5$ (với $p_s = 0.5$), tại giá trị $\lambda = 10$, tỉ lệ neutron bị hấp thụ tăng lên mức 8,18% (tăng tương đối 57,9% so với trường hợp $p_c = 0.3333$). Tỉ lệ phản xạ ghi nhận giảm nhẹ còn 52,88%, trong khi tỉ lệ truyền qua đạt mức 38,94%.

Thứ ba, trong bài toán tán xạ phi đàn hồi với $\lambda = 0.3$ và $p_c = 0.3333$, sự mất mát năng lượng $f$ sau mỗi va chạm đóng vai trò quyết định đến khả năng đâm xuyên. Khi năng lượng bị triệt tiêu hoàn toàn sau va chạm ($f = 1,0$ hay mất mát 100%), tỉ lệ hạt bị hấp thụ tuyệt đối đạt 100% (10.000 hạt), triệt tiêu hoàn toàn tỉ lệ phản xạ (0%) và tỉ lệ truyền qua (0%). Ngược lại, khi không có sự mất mát năng lượng ($f = 0$), tỉ lệ hấp thụ chỉ là 34,78%, tỉ lệ phản xạ đạt 63,50% và tỉ lệ truyền qua đạt 1,72%.

Thảo luận kết quả

Các kết quả mô phỏng phản ánh chính xác bản chất vật lý vi mô của tương tác neutron. Khi năng lượng neutron tới tăng cao, vận tốc $v$ của hạt tăng lên, dẫn đến quãng đường tự do trung bình $\lambda$ (tỉ lệ thuận với $\sqrt{E}$) tăng theo. Theo quy luật cơ học lượng tử, thời gian neutron lưu lại trong vùng tương tác của lực hạt nhân tỉ lệ với đại lượng nghịch đảo vận tốc $1/v$. Do đó, tiết diện hấp thụ giảm nhanh, làm cho tỉ lệ hạt bị bắt giữ giảm mạnh và tỉ lệ hạt truyền qua tấm vật liệu tăng tuyến tính theo năng lượng.

Tỉ lệ phản xạ ngược về phía nguồn phát luôn duy trì ở mức cao (thường lớn hơn 50% ở dải năng lượng thấp và trung bình) do nguồn phát ban đầu có tính chất đẳng hướng trong không gian góc $4\pi$. Một số lượng lớn neutron sau khi đi vào lớp mỏng ở bề mặt bia đã trải qua các va chạm tán xạ đàn hồi đổi hướng và bay ngược trở lại nửa không gian phát xạ trước khi có cơ hội đi sâu vào cấu trúc bên trong.

Đối với trường hợp tán xạ phi đàn hồi, sự mất mát năng lượng $f$ làm cho neutron sau mỗi lần va chạm bị chậm lại đáng kể, đưa hạt rơi vào miền năng lượng thấp nơi có tiết diện bắt bức xạ cộng hưởng rất lớn. Các dữ liệu này được thể hiện trực quan qua 9 biểu đồ biến thiên liên tục và 3 bảng số liệu tổng hợp trong luận văn, mô tả chi tiết đường cong phân rã của dòng hạt theo bề dày và năng lượng. Kết quả này hoàn toàn tương thích với các mô hình lý thuyết của Viện Năng lượng Nguyên tử Việt Nam và dữ liệu đo đạc chuẩn quốc tế ENDF/B-VI.

Đề xuất và khuyến nghị

Dựa trên các kết quả đạt được từ chương trình mô phỏng, nghiên cứu đưa ra 4 khuyến nghị và giải pháp phát triển kỹ thuật cụ thể:

Thứ nhất, nâng cấp mô hình mô phỏng từ cấu trúc phẳng 2 chiều lên không gian 3 chiều hoàn chỉnh. Nhóm nghiên cứu vật lý tính toán cần bổ sung việc lấy mẫu góc phương vị $\phi$ trong khoảng từ 0 đến $2\pi$ và cập nhật tọa độ không gian $(x, y, z)$. Mục tiêu đặt ra là mở rộng khả năng mô phỏng với cỡ mẫu 100.000 hạt trong vòng 6 tháng tới.

Thứ hai, tích hợp cơ sở dữ liệu tiết diện hạt nhân liên tục ENDF/B-VIII thay vì sử dụng các giá trị xác suất hấp thụ trung bình cố định. Chủ thể thực hiện là các chuyên gia vật lý hạt nhân, hướng tới mục tiêu giảm sai số tính toán thông lượng neutron xuống dưới 0,5% trong thời hạn 12 tháng.

Thứ ba, tối ưu hóa thuật toán bằng kỹ thuật tính toán song song đa luồng OpenMP hoặc GPU CUDA trên nền tảng C. Giải pháp này do các kỹ sư phần mềm khoa học phối hợp thực hiện nhằm tăng tốc độ xử lý lên gấp 500% (tăng 5 lần hiệu suất thời gian thực) trong lộ trình 9 tháng.

Thứ tư, đóng gói mã nguồn Scilab thành bộ công cụ mở có giao diện đồ họa trực quan (GUI) phục vụ công tác đào tạo thực hành vật lý hạt nhân tại các trường đại học. Đơn vị chủ trì là Bộ môn Vật lý lý thuyết với mục tiêu hoàn thiện và chuyển giao bộ tài liệu hướng dẫn trong vòng 18 tháng.

Đối tượng nên tham khảo luận văn

Luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng trọng tâm sau:

Nhóm học viên cao học và sinh viên chuyên ngành Vật lý hạt nhân, Vật lý lý thuyết: Tài liệu cung cấp ví dụ thực hành mẫu mực về cách thức chuyển đổi phương trình vi tích phân phức tạp sang thuật toán lấy mẫu thống kê, hỗ trợ nắm vững cơ chế phân hạch, bắt nơtron và tán xạ hạt nhân.

Nhóm kỹ sư thiết kế lò phản ứng và che chắn an toàn bức xạ: Cung cấp các công thức ước lượng độ suy giảm dòng hạt, giúp tính toán độ dày vật liệu che chắn tối ưu nhằm giảm thiểu thông lượng neutron rò rỉ ra ngoài môi trường.

Nhóm giảng viên và nhà nghiên cứu Vật lý tính toán: Luận văn là tài liệu tham khảo giá trị về phương pháp Monte Carlo, thuật toán sinh số ngẫu nhiên đồng dư tuyến tính và phương pháp loại trừ Von Neumann ứng dụng trong ngôn ngữ Scilab và C.

Nhóm chuyên gia phát triển detector và thiết bị đo phóng xạ: Giúp hiểu rõ cơ chế tương tác đa lớp và sự phân bố của dòng neutron phản xạ, hỗ trợ tối ưu hóa cấu trúc các đầu dò neutron nhiệt và neutron nhanh.

Câu hỏi thường gặp

Tại sao phương pháp Monte Carlo lại vượt trội hơn phương pháp giải tích tất định trong bài toán vận chuyển neutron?
Phương pháp tất định đòi hỏi giải phương trình vi tích phân Boltzmann 7 chiều vô cùng phức tạp và buộc phải dùng nhiều giả thiết xấp xỉ thô. Trong khi đó, Monte Carlo theo dõi trực tiếp từng lịch sử va chạm ngẫu nhiên của hạt nhân, có sai số giảm dần theo $1/\sqrt{N}$ và hoàn toàn không bị ảnh hưởng bởi sự gia tăng số chiều hình học.

Quãng đường tự do trung bình của neutron phụ thuộc vào yếu tố nào?
Quãng đường tự do trung bình $\lambda$ tỉ lệ nghịch với tiết diện tương tác vĩ mô toàn phần ($\lambda = 1/\Sigma_t$) và tỉ lệ thuận với căn bậc hai của động năng neutron ($\sqrt{E}$). Năng lượng hạt càng cao, vận tốc càng lớn thì quãng đường tự do trung bình càng dài, giúp neutron đâm xuyên sâu hơn qua vật liệu.

Thuật toán loại trừ Von Neumann được ứng dụng như thế nào trong mô phỏng?
Thuật toán này được dùng để sinh biến ngẫu nhiên theo một hàm phân bố mật độ xác suất bất kỳ mà không cần tìm hàm tích phân ngược. Bằng cách gieo cặp số ngẫu nhiên đều và chỉ chấp nhận các điểm nằm dưới đường cong phân bố, thuật toán đảm bảo mô phỏng chính xác góc tán xạ vi mô.

Mất mát năng lượng trong tán xạ phi đàn hồi ảnh hưởng thế nào đến khả năng che chắn bức xạ?
Tán xạ phi đàn hồi làm suy giảm $f%$ động năng của neutron sau va chạm, nhanh chóng đưa neutron về vùng năng lượng thấp hoặc miền cộng hưởng. Tại đây, tiết diện bắt nơtron tăng vọt từ 100 đến 1000 lần, khiến tỉ lệ hấp thụ đạt tới 100% khi $f=1,0$, giúp tăng hiệu quả che chắn bảo vệ.

Việc sử dụng Scilab mang lại lợi ích gì so với các phần mềm thương mại như MATLAB hay MCNP?
Scilab là phần mềm mã nguồn mở hoàn toàn miễn phí, dung lượng nhẹ và cú pháp tính toán ma trận tương tự MATLAB. Việc triển khai trên Scilab giúp người học dễ dàng đọc hiểu thuật toán, chỉnh sửa mã nguồn và mô phỏng trực tiếp mà không tốn chi phí bản quyền đắt đỏ.

Kết luận

• Luận văn đã hệ thống hóa toàn diện cơ sở vật lý tương tác neutron và xây dựng thành công chương trình mô phỏng Monte Carlo 2D trên Scilab và C.
• Kết quả mô phỏng với 10.000 hạt khẳng định rõ quy luật: khi quãng đường tự do $\lambda$ tăng từ 0.1 lên 10, tỉ lệ hạt truyền qua tăng từ 0% lên 39,72% với $p_c = 0.3333$.
• Đã định lượng chính xác vai trò của tán xạ phi đàn hồi, chứng minh sự mất mát năng lượng 100% sẽ làm tỉ lệ hạt bị hấp thụ đạt mức tuyệt đối 100%.
• Giải pháp mã nguồn mở Scilab giúp tiết kiệm 100% chi phí bản quyền, tạo công cụ học tập trực quan và hiệu quả cho sinh viên ngành vật lý.
• Đề ra lộ trình nâng cấp thuật toán sang không gian 3 chiều và tích hợp thư viện hạt nhân ENDF/B-VIII trong giai đoạn tiếp theo.

Công trình là bước đi nền tảng vững chắc trong việc làm chủ kỹ thuật mô phỏng bức xạ hạt nhân. Quý độc giả, học viên và các nhóm nghiên cứu quan tâm có thể khai thác mã nguồn thuật toán để tiếp tục mở rộng mô hình cho các cấu hình lò phản ứng phức tạp hơn.

Trích đoạn nội dung tài liệu

mở đầu này chúng tôi tìm hiểu cơ sở vật lý của các tương tác neutron, những kiến thức cần thiết để nghiên cứu mô phỏng Monte Carlo bài toán vận chuyển neutron. Tương tác vật lý của neutron được mô tả cả bằng cơ học cổ điển và cơ học lượng tử. Các phản ứng tán xạ đàn hồi ở mức năng lượng thấp thuộc phạm vi ứng dụng của cơ học cổ điển. Cơ học lượng tử đóng vai trò quan trọng trong các hiện tượng như: quá trình neutron đâm xuyên sâu vào trong hạt nhân bia và tương tác với các cấu phần của hạt nhân; một vài hiện tượng thể hiện tính sóng của neutron như nhiễu xạ neutron xuất hiện khi các neutron có năng lượng thấp tương tác với cấu trúc tinh thể của vật liệu,v.

Ngoài tán xạ đàn hồi, các tương tác neutron khác thường gặp trong các quá trình phân hạch không quá khác nhau. Chúng đều tạo ra các hạt nhân hợp phần ở trang thái kích thích với thời gian tồn tại ngắn và nhanh chóng phân rã để trở lại trạng thái cơ bản (ground state). Mode phân rã sẽ quy định loại phản ứng, trong quá trình tính toán bài toán vận chuyển neutron, chúng có thể được chia ra thành 3 loại phản ứng: hấp thụ, phân hạch và tán xạ. Mỗi loại phản ứng xảy ra với một giá trị cận dưới của động năng của neutron tới, nếu động năng nhỏ hơn giá trị này thì phản ứng không xảy ra.

Khi hấp thụ neutron hạt nhân trở thành hạt nhân hợp phần với năng lượng kích thích tối thiểu bằng năng lượng liên kết của neutron trong hạt nhân đó. Nếu năng lượng liên kết này lớn hơn một giá trị ngưỡng nào đó (phụ thuộc vào mỗi đồng vị) thì hạt nhân ban đầu có thể bị phân hạch khi hấp thụ neutron với năng lượng bất kỳ. Nếu năng lượng liên kết nhỏ hơn giá trị tới hạn này thì quá trình phân hạch chỉ xảy 3 LUAN VAN CHAT LUONG download : add [e-mail] ra khi động năng của neutron đủ lớn để cho năng lượng kích thích vượt quá năng lượng ngưỡng. Xác suất tương tác giữa neutron tới và hạt nhân bia được đặc trưng bởi tiết diện tương tác vi mô, phụ thuộc vào hạt nhân bia, loại tương tác và năng lượng của neutron tới.

Ở năng lượng thấp, xác xuất xảy ra các phản ứng hợp phần tỉ lệ trực tiếp với thời gian mà neutron đi qua vùng tương tác của lực tương tác mạnh. Điều này dẫn đến sự phụ thuộc của các tiết diện tán xạ vào đại lượng nghịch đảo vận tốc của neutron, 1/v. Sự phụ thuộc của các tiết diện phóng xạ vào các đại lượng vật lý trở nên phức tạp hơn khi năng lượng tăng và các tương tác ở cấp độ vi mô cần đến sự vận dụng của cơ học lượng tử xuất hiện, dẫn đến sự thay đổi phức tạp trong tiết diện phản ứng. Nếu năng lượng toàn phần của neutron xấp xỉ bằng năng lượng của trạng thái hợp phần kích thích, xác suất tương tác có thể tăng lên vài bậc.

Điều này có thể được quan sát với các đỉnh cộng hưởng trên đồ thị sự phụ thuộc của tiết diện tương tác vào năng lượng. Khoảng cách giữa các đỉnh cộng hưởng thu hẹp lại khi năng lượng tăng. Đến một lúc nào đó, không thể quan sát được các đỉnh riêng rẽ. Vùng này là miền cộng hưởng không thể phân giải.

Ở mức năng lượng cao hơn nữa, sự cộng hưởng trở nên trùng lặp và chúng ta quan sát được dạng đường cong trơn. Sự hấp thụ Hình 1. Tiết diện phản ứng (n, γ) của 238U và tiết diện phân hạch của 235U, 238 U và 239Pu (Dữ liệu của thí nghiệm ENDF/B-VI, 300K) 4 LUAN VAN CHAT LUONG download : add [e-mail] Sự hấp thụ đề cập đến tất cả các phản ứng mà không có neutron thứ cấp được phát ra. Hiện tượng hấp thụ đơn giản nhất là hiện tượng bắt neutron bức xạ hay phản ứng (n,γ).

Neutron tới được giữ lại trong các hạt nhân bia và năng lượng dự thừa được giải phóng bằng sự phát xạ photon năng lượng cao. Một ví dụ cho phản ứng bắt bức xạ neutron của 238U : 238 239 U+n  U* + γ (1.1) Trong thực tế phản ứng (n,γ) là phương thức hấp thụ duy nhất với các đồng vị actini giàu neutron. Các phản ứng phân hạch và tán xạ chiếm ưu thế ở vùng năng lượng cao. Tiết diện tương tác phản ứng (n, γ) của 238U được biểu diễn trong hình 1.

Các phản ứng bắt quan trọng khác gồm (n,p), (n,α) và (n,3He). Chúng là các phản ứng điển hình chỉ xảy ra khi với một ngưỡng năng lượng của neutron đến. Một ví dụ là phản ứng (n,α) của 10B: 10 11 B+n  B*  7 Li + 4He (1.2) Urani và các actini là những đồng vị nặng với hơn 90 proton và 140–160 neutron trong hạt nhân. Các phản ứng hợp phần rất phức tạp và các tiết diện tương tác của actini được đặc trưng bởi số lượng lớn các vùng cộng hưởng.

Đỉnh cộng hưởng đầu tiên được phân tách khá rõ ràng trong miền năng lượng tương đối thấp. Dải cộng hưởng không phân giải được bắt đầu từ khoảng 10 KeV. Miền liên tục này được thể hiện trong hình là một đường trơn liên tục. Sự phân hạch Sự phân hạch là tương tác trong đó các hạt nhân hợp phần ở trạng thái kích thích sau khi hấp thụ neutron chậm bị phân chia thành 2 mảnh: 235 236 U+n U*  135 Xe + 98Sr + 3n (1.3) Trong phản ứng này, sản phẩm trung gian là 236U sẽ bị phân rã thành các hạt nhân nhỏ hơn.

Thông thường thì sản phẩm tạo ra sẽ là 2 đồng vị có số khối trung bình với một đồng vị nặng hơn đồng vị kia một chút. Sản phẩm của phản ứng phân hạch thường có nhiều neutron và bị phân rã β- để đi đến trạng thái bền vững. 5 LUAN VAN CHAT LUONG download : add [e-mail] Phản ứng phân hạch có thể xảy ra đối với tất cả các nguyên tốt thuộc họ actini và chiếm ưu thế trong vùng năng lượng cao cỡ MeV. Trong thực tế, phản ứng cần năng lượng khoảng 6 MeV để có thể xảy ra trong khi năng lượng sinh ra trong quá trình tạo thành hạt nhân hợp phần là khoảng 5 MeV.

Sự chênh lệch năng lượng này làm cho phản ứng phụ thuộc lớn vào hiệu ứng chẵn lẽ: sự thay đổi số lượng neutron từ lẻ sang chẵn giải phóng thêm vài MeV, đủ để nâng năng lượng lên trên ngưỡng phân hạch. Các đồng vị có số neutron lẻ (235U, 239Pu, 241Pu, 242Am,…) là những đồng vị phân hạch, chúng có thể bị phân hạch bởi sự hấp thụ các neutron có năng lượng thấp. Các đồng vị phân hạch có số neutron chẵn (238U, 237Np, 241Am…) chỉ phân hạch khi hấp thụ các neutron có động năng lớn hơn giá trị ngưỡng, khoảng 1 MeV. Các tiết diện phân hạch của 235U, 238U, và 239Pu cũng được vẽ trong hình 1.

Phản ứng phân hạch tạo ra 1 đến 6 neutron mới. Số neutron trung bình sinh ra là hàm tuyến tính của năng lượng neutron tới và nó cũng phụ thuộc vào số khối của đồng vị. Ngoài các neutron tức thời phát ra trong sự kiện phân hạch, còn có nhiều neutron được tạo ra từ các chuỗi phân hạch của các đồng vị sản phẩm có số lượng neutron dư. Sự phóng xạ neutron dư hầu như là tức thời nhưng nó cũng có thể xảy ra sau sự phân rã của đồng vị phóng xạ khác có thời gian sống dài hơn.

Một ví dụ của chuỗi phản ứng là sự phân rã của sản phẩm phân hạch 87Br: 87 Br  87Kr*  86Kr + n (1.4) Hạt nhân 87Br ở trạng thái kích thích với chu kì bán rã 56 s và phân rã tạo nên hạt nhân 87Kr. Hạt nhân 87Kr được tạo nên ở trạng thái kích thích cao với năng lượng 6 MeV đủ để phóng ra tức thời một neutron và chuyển thành hạt nhân 86Kr bền. Thời gian sống của chuỗi phân rã được chi phối bởi phản ứng đầu tiên với chu kì bán rã là 56s. Phản ứng này quyết định thời gian mà neutron có thể được phát ra sau sự kiện phân hạch.

Các neutron trễ này có vai trò quan trọng đối với hằng số thời gian trong các lò phản ứng hạt nhân. 6 LUAN VAN CHAT LUONG download : add luanvanchat@agmail. Số lượng trung bình các neutron phân hạch tức thời của 235U, 238U, 239 Pu và 241Am (Dữ liệu của thí nghiệm ENDF/B-VI).2 cho thấy số lượng neutron phát ra tăng lên theo số khối của đồng vị và năng lượng neutron tới. Một đại lượng quan trọng trong vật lý lò là tỉ lệ của số neutron trễ được vẽ trên hình 1.3 với 4 đồng vị.

Tỉ lệ này với các đồng vị urani cao khiến cho những những nhiên liệu làm từ urani tốt hơn nhiều so với các dạng nhiên liệu khác. Tỉ lệ neutron trễ với các đồng vị 235U, 238U, 239Pu và 241Am 7 LUAN VAN CHAT LUONG download : add [e-mail] Neutron phân hạch được phát ra đẳng hướng do có rất nhiều hạt được sinh ra trong phản ứng mà các định luật bảo toàn năng lượng, động lượng không chi phối hướng tới và phát. Hơn nữa thời gian sống của trạng thái trung gian đủ dài để các hạt neutron tới “quên” đi hướng cũ của nó. Phân bố năng lượng của neutron phân hạch tức thời có dạng hàm phân bố Maxwell, có đỉnh ở quanh giá trị 1 MeV.

Đoạn cuối của đường phân bố kéo dài tới khoảng năng lượng cao nhất là 10 MeV và kết thúc ở khoảng năng lượng thấp nhất là 100 eV. Phổ của phân bố phân hạch tức thời của 235U, 238U, 239Pu, 241Am được vẽ trong hình 1. Các phân bố hơi khác nhau đối với các đồng vị khác nhau, nhưng hầu như không phụ thuộc vào năng lượng neutron tới. Các neutron phân rã được sinh ra ở các mức năng lượng thấp, có vai trò quan trọng trong tính toán lò phản ứng.

Năng lượng toàn phần sinh ra trong phản ứng phân hạch (Q) khoảng 200 MeV và nó có sự phụ thuộc nhỏ vào đồng vị sử dụng làm bia. Phần lớn năng lượng, khoảng 80% được giải phóng dưới dạng động năng của các mảnh phân hạch. Năng lượng này được sinh ra ngay tại vị trí phân hạch. Tuy nhiên ta không thể thu được toàn bộ 200 MeV năng lượng phân hạch do các neutron và photon bị rò ra bên ngoài lõi lò.

Khoảng 4% năng lượng bị mất do các hạt neutrino không tương tác với vật chất. Không phải tất cả năng lượng được giải phóng ngay. Các phân rã phóng xạ của các sản phẩm phân hạch chiếm khoảng 6% tổng năng lượng phân hạch với hằng số phóng xạ biến thiên từ vài phần nghìn giây cho đến hàng nghìn năm.