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.

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.