MCNP5 (Monte Carlo N-Particle Transport Code) là chương trình mô phỏng vật lý hạt tổng quát sử dụng phương pháp Monte Carlo để tính toán sự truyền dẫn bức xạ neutron, photon và electron trong các môi trường vật chất có cấu trúc hình học không gian ba chiều phức tạp.
Lịch sử phát triển và kiến trúc phần mềm MCNP5
Hệ mã mô phỏng MCNP có cội nguồn từ Dự án Manhattan tại Phòng thí nghiệm Quốc gia Los Alamos (LANL), nơi các nhà vật lý Stanislaw Ulam, John von Neumann và Enrico Fermi lần đầu tiên áp dụng phương pháp thống kê ngẫu nhiên Monte Carlo để giải quyết các bài toán khuếch tán neutron trong vũ khí hạt nhân. Trải qua nhiều thập kỷ hoàn thiện, từ các phiên bản MCNP3, MCNP4 đến MCNP5 (phát hành đầu những năm 2000), hệ mã đã được hiện đại hóa hoàn toàn sang ngôn ngữ Fortran 90 kết hợp C, hỗ trợ tính toán song song đa lõi thông qua giao thức OpenMP trên bộ nhớ dùng chung và MPI trên hệ thống cụm máy tính phân tán (clusters).
MCNP5 cung cấp khả năng mô phỏng liên tục theo năng lượng cho ba loại hạt cơ bản: neutron (từ 10^-11 MeV đến 20 MeV, và mở rộng tới 150 MeV khi liên kết thư viện vật lý năng lượng cao), photon (từ 1 keV đến 100 GeV) và electron (từ 1 keV đến 1 GeV), đồng thời giải quyết bài toán vận chuyển ghép cặp neutron - photon - electron một cách nhất quán.
Nguyên lý toán học của thuật toán Monte Carlo trong truyền dẫn bức xạ
Khác biệt căn bản với các phương pháp tất định giải phương trình vi tích phân vận chuyển Boltzmann (như phương pháp tọa độ rời rạc Discrete Ordinates Sn hoặc phương pháp phần tử hữu hạn FEM), phương pháp Monte Carlo không đòi hỏi phải rời rạc hóa không gian, góc hay phổ năng lượng. Thay vào đó, nó mô phỏng trực tiếp lịch sử sống (particle history) của từng hạt riêng lẻ từ khi được sinh ra ở nguồn bức xạ cho đến khi bị hấp thụ hoàn toàn hoặc thoát khỏi hệ thống hình học.
Tại mỗi bước đi tự do của hạt, quãng đường đi tự do trung bình, loại phản ứng tương tác hạt nhân (tán xạ đàn hồi, tán xạ không đàn hồi, phân hạch, bắt bức xạ hay hiệu ứng quang điện, tán xạ Compton, tạo cặp) và năng lượng cũng như góc tán xạ sau va chạm đều được chọn ngẫu nhiên dựa trên các hàm mật độ xác suất (PDF) rút ra từ các tập tiết diện vi phân thực nghiệm. Theo Định lý Giới hạn Trung tâm, giá trị kỳ vọng của đại lượng vật lý cần đo (ví dụ thông lượng neutron hay liều hấp thụ) sẽ hội tụ về nghiệm thực khi số lượng lịch sử hạt mô phỏng N tiến tới vô cùng, với sai số thống kê tỷ lệ nghịch với căn bậc hai của N (1/sqrt(N)).
Thư viện dữ liệu tiết diện hạt nhân (Nuclear Cross-Section Libraries)
Độ chính xác tuyệt đối của kết quả tính toán MCNP5 phụ thuộc chặt chẽ vào độ tin cậy của các thư viện dữ liệu hạt nhân được đánh giá quốc tế (Evaluated Nuclear Data File - ENDF). MCNP5 tích hợp và sử dụng định dạng bảng dữ liệu ACE (A Compact ENDF) được tiền xử lý bởi mã NJOY, bao gồm:
- Thư viện ENDF/B-VII và ENDF/B-VIII: Cung cấp tiết diện tán xạ, phân hạch và bắt neutron cho hàng trăm đồng vị với độ phân giải năng lượng liên tục cực cao, tái hiện chính xác các đỉnh cộng hưởng nhiệt và cộng hưởng giải được.
- Dữ liệu tán xạ nhiệt S(alpha, beta): Cần thiết để mô phỏng tương tác của neutron nhiệt chậm với các phân tử liên kết hóa học trong chất làm chậm (như nước nhẹ H2O, nước nặng D2O, graphite, polyethylene), nơi chuyển động dao động và quay của mạng tinh thể ảnh hưởng lớn đến quá trình nhiệt hóa.
- Thư viện photon và electron: Tích hợp dữ liệu tiết diện tương tác nguyên tử EPDL97 và EEDL mô tả chính xác hiệu ứng quang điện, phát xạ tia hãm (bremsstrahlung) và ion hóa vỏ nguyên tử.
| Loại hạt | Dải năng lượng mô phỏng | Các quá trình tương tác chính | Thư viện tiết diện chuẩn |
|---|---|---|---|
| Neutron | 10^-11 MeV đến 20 MeV (mở rộng 150 MeV) | Tán xạ đàn hồi/không đàn hồi, phân hạch sinh neutron (n,f), bắt bức xạ (n,gamma), phản ứng sinh hạt tích điện (n,p), (n,alpha) | ENDF/B-VII.1, JEFF-3.3, JENDL-4.0 |
| Photon | 1 keV đến 100 GeV | Tán xạ kết hợp (Thomson/Rayleigh), tán xạ không kết hợp (Compton), hấp thụ quang điện, tạo cặp electron-positron | EPDL97, MCPLIB04 |
| Electron | 1 keV đến 1 GeV | Tán xạ đàn hồi hạt nhân (Molière), ion hóa va chạm, kích thích nguyên tử, phát xạ bức xạ hãm (Bremsstrahlung) | EEDL, EL03 (mô hình ngẫu nhiên hóa bước dừng condensed history) |
Cấu trúc tệp đầu vào (Input File Structure)
Một tệp đầu vào chuẩn của MCNP5 (thường có phần mở rộng .inp) được thiết kế theo cấu trúc ba khối thẻ phân cách nhau bằng một dòng trắng:
- Khối thẻ ô không gian (Cell Cards): Định nghĩa các miền vật lý trong không gian ba chiều bằng các phép toán logic đại số Boole (giao, hợp, phần bù) của các mặt biên hình học, gán số hiệu vật liệu, mật độ nguyên tử và tham số tầm quan trọng (importance).
- Khối thẻ mặt biên (Surface Cards): Khai báo các mặt phẳng, mặt cầu, mặt trụ, mặt nón hoặc các đa thức bậc hai tổng quát với các thông số tọa độ cụ thể để giới hạn các ô không gian.
- Khối thẻ dữ liệu (Data Cards): Khai báo định nghĩa thành phần đồng vị vật liệu (thẻ M), phân bố nguồn bức xạ tổng quát (thẻ SDEF), chế độ vận chuyển hạt (thẻ MODE N P E), cấu hình bộ ghi nhận đại lượng đo (thẻ Tally F1 đến F8), và các lệnh điều khiển kiểm tra hội tụ số.
Các kỹ thuật giảm phương sai (Variance Reduction Techniques)
Trong các bài toán che chắn bức xạ dày (ví dụ tường bê tông lò phản ứng dày 2 mét) hoặc phát hiện hạt ở khoảng cách xa, xác suất để một hạt tự nhiên xuyên qua che chắn chỉ vào khoảng 10^-8 đến 10^-12. Nếu mô phỏng tự nhiên, máy tính cần hàng triệu tỷ lịch sử hạt để thu được một kết quả có ý nghĩa thống kê. MCNP5 cung cấp các kỹ thuật giảm phương sai tinh vi giúp tăng tốc tính toán hàng nghìn lần:
- Tách ô hình học và Roulette Nga (Geometry Splitting & Russian Roulette): Gán tham số tầm quan trọng (importance) tăng dần khi hạt tiến gần về phía vùng đo quan tâm; khi hạt vượt qua ranh giới vào ô có importance gấp đôi, nó được tách thành 2 hạt với trọng số giảm một nửa. Ngược lại, khi hạt di chuyển ra xa, nó chịu bài toán roulette Nga để triệt tiêu các hạt không có giá trị thống kê.
- Trọng số cửa sổ (Weight Window): Thiết lập ngưỡng chặn trên và chặn dưới cho trọng số hạt theo từng vùng không gian và dải năng lượng, đảm bảo trọng số của các hạt luôn dao động trong khoảng tối ưu.
- Va chạm cưỡng bức (Forced Collisions): Bắt buộc hạt phải xảy ra va chạm bên trong một ô vật chất cụ thể, hữu ích trong các bài toán đầu dò khí hoặc chất nhấp nháy có mật độ mỏng.
- Thiên vị nguồn (Source Biasing): Phát hạt thiên lệch theo các hướng, năng lượng hoặc vị trí có khả năng đến được vùng thu nhận cao hơn.
Ứng dụng thực tiễn trong kỹ thuật hạt nhân và y học bức xạ
MCNP5 là công cụ tiêu chuẩn vàng được kiểm chứng trong nhiều ngành công nghiệp và nghiên cứu mũi nhọn:
- Phân tích an toàn và độ tới hạn lò phản ứng: Tính toán hệ số nhân neutron hiệu dụng (k-effective) theo chu trình lặp Shannon entropy, phân bố công suất vùng lõi và thiết kế thanh điều khiển an toàn.
- Thiết kế che chắn bảo vệ bức xạ: Tối ưu hóa chiều dày và vật liệu che chắn bức xạ hỗn hợp neutron - gamma cho container lưu giữ nhiên liệu hạt nhân đã qua sử dụng và cơ sở bức xạ công nghiệp.
- Y học bức xạ và xạ trị ung thư: Tính toán bản đồ phân bố liều hấp thụ 3D cho các kỹ thuật xạ trị phức tạp, thiết kế đầu máy xạ trị gia tốc tuyến tính (LINAC), và tối ưu hóa liều trong kỹ thuật xạ trị bắt giữ neutron bằng hợp chất boron (BNCT).