Xây Dựng Mô Hình Tính Toán Kết Hợp Vật Lý – Thủy Nhiệt Cho Bó Nhiên Liệu Lò Phản Ứng VVER

Tài liệu nghiên cứu Xây dựng mô hình tính toán kết hợp vật lý thủy nhiệt ho bó nhiên liệu lò phản ứng vver sử dụng, tổng hợp lý thuyết và thực hành, cung cấp kiến thức chuyên sâu

• Tác giả: Tạ Duy Long
• Người hướng dẫn: TS. Trần Hoài Nam
• Trường học: Trường Đại Học Bách Khoa Hà Nội
• Chuyên ngành: Kỹ Thuật Hạt Nhân
• Đề tài: Mô Hình Tính Toán Kết Hợp Vật Lý – Thủy Nhiệt Cho Bó Nhiên Liệu Lò Phản Ứng VVER
• Loại tài liệu: Luận Văn Thạc Sĩ
• Năm xuất bản: 2016
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CAM ĐOAN
• LỜI CẢM ƠN
• 1. CHƯƠNG 1: LÝ DO LỰA CHỌN ĐỀ TÀI
• 1.1. Mục tiêu nghiên cứu. Nội dung nghiên cứu
• 1.2. Xây dựng mô hình tính toán kết hợp Vật lý – Thủy nhiệt
• 1.3. Mô tả các chương trình tính toán
• 1.3.1. Chương trình tính toán MCNP5
• 1.3.2. Chương trình tính toán COBRA-EN
• 1.4. Tổng quan về lò phản ứng VVER-1000
• 1.5. Tổng quan về các hiệu ứng phản hồi nhiệt độ
• 1.6. Ảnh hưởng của mật độ chất tải nhiệt lên tính toán vùng hoạt lò phản ứng
• 2. CHƯƠNG 2: CÁC KẾT QUẢ TÍNH TOÁN
• 2.1. Tính toán benchmark với chương trình MCNP5
• 2.2. Các kết quả tính toán kết hợp Vật lý – Thủy nhiệt
• 2.3. Phân bố công suất trong bó nhiên liệu
• 2.4. Phân bố nhiệt độ chất tải nhiệt theo chiều cao bó nhiên liệu
• 2.5. Phân bố nhiệt độ trong thanh nhiên liệu theo chiều bán kính
• 3. CHƯƠNG 3: MÔ HÌNH HÓA BÓ NHIÊN LIỆU LÒ PHẢN ỨNG VVER-1000
• 3.1. Mô hình hóa bó nhiên liệu lò phản ứng VVER-1000 với chương trình MCNP5 và COBRA-EN
• 3.2. Mô hình COBRA-EN
• 3.3. Chương trình tính toán kết hợp
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Mô Hình Vật Lý Thủy Nhiệt VVER Bí Quyết

Hoạt động lò phản ứng hạt nhân là hệ rất phức tạp. Nó bao gồm tương tác vật lý như quá trình chuyên neutron, quá trình chuyên chất lỏng và truyền nhiệt. Việc mô phỏng các quá trình trong hoạt động lò phản ứng là một trong những nhiệm vụ quan trọng. Tính toán kết hợp Vật lý - Thủy nhiệt là hướng đi triển khai bởi nhiều nhóm nghiên cứu. Các mô hình tính toán kết hợp Vật lý và Thủy nhiệt đóng vai trò quan trọng trong việc đảm bảo độ chính xác của các tính toán. Các tính toán phụ thuộc vào trạng thái chuyên tiếp của lò phản ứng. Các tính toán ba chiều kết hợp Vật lý- Thủy nhiệt hiện được sử dụng để tính toán phân bố công suất trong một bó nhiên liệu lò phản ứng. Theo Tạ Duy Long (2016), nghiên cứu này xây dựng mô hình tính toán kết hợp cho bó nhiên liệu VVER sử dụng MCNP5 và COBRA-EN. Điều này nhằm cung cấp công cụ thiết yếu cho phân tích an toàn và hiệu suất lò.

1.1. Tầm quan trọng của Mô hình kết hợp VVER

Các tính toán kết hợp vật lý-thủy nhiệt đang ngày càng trở nên quan trọng trong việc đánh giá an toàn và hiệu suất của lò phản ứng. Việc sử dụng các mã mô phỏng tiên tiến như MCNP5 và COBRA-EN cho phép mô hình hóa chi tiết các quá trình neutron và thủy động lực học xảy ra trong bó nhiên liệu. Các mô hình này cung cấp thông tin chi tiết về phân bố nhiệt độ, mật độ và lưu lượng chất làm mát, cho phép các kỹ sư đánh giá chính xác hơn các điều kiện vận hành và dự đoán hành vi của lò phản ứng trong các kịch bản tai nạn khác nhau.

1.2. Thách thức và Triển vọng Mô hình kết hợp

Mặc dù các mô hình kết hợp vật lý-thủy nhiệt mang lại nhiều lợi ích, nhưng việc phát triển và xác nhận chúng cũng đặt ra những thách thức đáng kể. Các mô hình này thường đòi hỏi tài nguyên tính toán lớn và cần được xác nhận bằng dữ liệu thực nghiệm để đảm bảo độ chính xác. Tuy nhiên, với sự tiến bộ không ngừng của công nghệ máy tính và sự sẵn có ngày càng tăng của dữ liệu thử nghiệm, các mô hình kết hợp vật lý-thủy nhiệt hứa hẹn sẽ đóng vai trò ngày càng quan trọng trong thiết kế, vận hành và đánh giá an toàn của các lò phản ứng hạt nhân trong tương lai.

II. Vấn Đề An Toàn Lò VVER Phân Tích và Giải Pháp Tính Toán

Các tính toán kết hợp Vật lý - Thủy nhiệt đặc biệt hữu ích cho các phân tích an toàn, cung cấp các giá trị hợp lý cho các thông số quan trọng. Mô hình này giúp xem xét ảnh hưởng của nhiệt độ trong các tính toán. Tính toán Vật lý lò phản ứng truyền thống xem xét các yếu tố nhiệt độ. Vì thế, việc xây dựng mô hình tính toán kết hợp Vật lý-Thủy nhiệt là cần thiết. Phân tích nhiệt độ trong các tính toán Vật lý của lò phản ứng. Tạ Duy Long (2016) đã xem xét ảnh hưởng nhiệt độ trong tính toán Vật lý lò phản ứng. Các kết quả tính toán kết hợp cho phép phân tích an toàn hơn, giảm độ bất định và các sai số.

2.1. Ảnh Hưởng Hiệu Ứng Nhiệt độ đến An Toàn Lò

Nhiệt độ là một trong những yếu tố quan trọng nhất ảnh hưởng đến hành vi của lò phản ứng hạt nhân. Sự thay đổi nhiệt độ có thể ảnh hưởng đến tốc độ phản ứng hạt nhân, phân bố công suất và tính chất của vật liệu lò phản ứng. Do đó, việc mô hình hóa chính xác các hiệu ứng nhiệt độ là rất quan trọng để đánh giá an toàn của lò phản ứng trong các điều kiện vận hành bình thường và các kịch bản tai nạn khác nhau.

2.2. Tính Toán Vật Lý và Thủy Nhiệt Liên Kết An Toàn

Các mô hình kết hợp vật lý-thủy nhiệt cho phép các kỹ sư mô phỏng đồng thời các quá trình neutron và thủy động lực học trong lò phản ứng. Điều này cho phép họ đánh giá tác động của các thay đổi nhiệt độ và lưu lượng chất làm mát đối với hiệu suất của lò phản ứng và sự ổn định của nó. Ví dụ, các mô hình này có thể được sử dụng để dự đoán sự hình thành điểm nóng trong bó nhiên liệu hoặc đánh giá hiệu quả của các hệ thống làm mát khẩn cấp.

2.3. Sai Số trong Tính Toán An Toàn Lò VVER

Phân tích độ nhạy cảm và đánh giá độ bất định là những khía cạnh quan trọng của phân tích an toàn lò phản ứng. Các mô hình kết hợp vật lý-thủy nhiệt có thể giúp xác định các nguồn độ bất định chính trong các tính toán an toàn và đánh giá tác động của chúng đối với kết quả. Điều này cho phép các kỹ sư đưa ra các quyết định sáng suốt hơn về thiết kế và vận hành của lò phản ứng.

III. MCNP5 COBRA EN Phương Pháp Tính Toán Kết Hợp VVER

Chương trình tính toán khởi chạy với tính toán vật lý bằng MCNP5, tạo ra phân bố công suất 3D. Giả sử phân bố nhiệt độ chất tải nhiệt đồng nhất. Phân bố công suất từ MCNP5 làm dữ liệu đầu vào cho COBRA-EN. Chương trình tính toán thủy nhiệt tạo ra phân bố 3D mật độ chất tải nhiệt. Kết quả này dùng làm dữ liệu cho MCNP5. Các tính toán được thực hiện đến khi điều kiện hội tụ thỏa mãn. Các chương trình hiện dùng tại Trung tâm Năng lượng Hạt nhân, Viện Khoa học và Kỹ thuật Hạt nhân. Để vận hành, tác giả phát triển các chương trình con để dùng dữ liệu đầu vào và một chương trình quản lý.

3.1. Mô tả chương trình MCNP5 trong Mô hình VVER

MCNP5 là một mã Monte Carlo linh hoạt có thể được sử dụng để mô phỏng vận chuyển neutron, photon và electron trong các hệ phức tạp. Trong mô hình VVER, MCNP5 được sử dụng để tính toán phân bố neutron trong vùng hoạt và xác định tỷ lệ phân hạch trong các thanh nhiên liệu khác nhau. Mã này có thể xử lý các hình học 3D phức tạp và có thể sử dụng nhiều thư viện dữ liệu hạt nhân khác nhau.

3.2. Ứng dụng COBRA EN cho Phân tích Thủy Nhiệt VVER

COBRA-EN là một mã thủy nhiệt có thể được sử dụng để mô phỏng dòng chảy và truyền nhiệt trong các bó nhiên liệu của lò phản ứng. Trong mô hình VVER, COBRA-EN được sử dụng để tính toán phân bố nhiệt độ và mật độ của chất làm mát trong các kênh khác nhau của bó nhiên liệu. Mã này có thể xử lý các hình học kênh phức tạp và có thể sử dụng nhiều mô hình tương quan truyền nhiệt và ma sát khác nhau.

3.3. Giao diện MCNP5 và COBRA EN Chi tiết kỹ thuật

Việc kết hợp MCNP5 và COBRA-EN đòi hỏi một giao diện cho phép trao đổi dữ liệu giữa hai mã. Giao diện này thường bao gồm một mã tập lệnh (script) chuyển đổi dữ liệu phân bố công suất từ MCNP5 sang định dạng mà COBRA-EN có thể sử dụng và sau đó chuyển đổi dữ liệu phân bố nhiệt độ và mật độ từ COBRA-EN sang định dạng mà MCNP5 có thể sử dụng. Quá trình này được lặp lại cho đến khi đạt được sự hội tụ giữa hai mã.

IV. Đánh Giá Hiệu Năng và An Toàn VVER qua Mô Hình

Nghiên cứu này giúp hiện thực hóa mô hình tính toán, kiểm định các thông số Vật lý và nhiệt độ trong các tính toán. Ảnh hưởng của mật độ chất tải nhiệt trong vùng hoạt lò là cần thiết. Các yếu tố có thể ảnh hưởng đến hệ số nhân gồm mật độ. Các tính toán giúp duy trì trạng thái hoạt động của lò. Trong nội dung luận văn, tác giả tập trung vào đánh giá ảnh hưởng của hiệu ứng mật độ chất làm mát.

4.1. Độ chính xác các Tính Toán Vật Lý VVER

Độ chính xác của các tính toán vật lý trong lò phản ứng VVER phụ thuộc vào độ chính xác của các thư viện dữ liệu hạt nhân được sử dụng và độ chi tiết của mô hình hình học. Các mô hình MCNP5 có thể sử dụng các thư viện dữ liệu hạt nhân khác nhau, chẳng hạn như ENDF/B-VII.1, để tính toán các tiết diện phản ứng và xác suất tương tác neutron.

4.2. Tính Thủy Nhiệt trong đánh giá An Toàn VVER

Các tính toán thủy nhiệt đóng vai trò quan trọng trong việc đánh giá an toàn của lò phản ứng VVER, vì chúng cung cấp thông tin chi tiết về phân bố nhiệt độ và lưu lượng chất làm mát trong bó nhiên liệu. Các mô hình COBRA-EN có thể sử dụng nhiều mô hình tương quan truyền nhiệt và ma sát khác nhau để dự đoán chính xác hành vi của chất làm mát trong các điều kiện vận hành khác nhau.

4.3. Kiểm Định Mô Hình Thủy Nhiệt và Kết quả Vật Lý

Việc kiểm định các mô hình thủy nhiệt và kết quả vật lý là rất quan trọng để đảm bảo độ chính xác của các tính toán an toàn. Các mô hình này có thể được kiểm định bằng cách so sánh kết quả tính toán với dữ liệu thực nghiệm thu được từ các lò phản ứng nghiên cứu hoặc các thí nghiệm tại chỗ.

V. Phân Bố Nhiệt và Mật Độ trong Bó Nhiên Liệu VVER Phân Tích

Chương trình khởi chạy với tính toán vật lý. Giả sử phân bố nhiệt độ chất tải nhiệt đồng nhất. Phân bố công suất từ chương trình dùng làm dữ liệu đầu vào. chương trình tính toán thủy nhiệt tạo ra phân bố mật độ chất tải nhiệt. Kết quả này làm dữ liệu cho tính toán vật lý. Các tính toán được thực hiện khi điều kiện hội tụ thỏa mãn.

5.1. Phân bố công suất Ảnh hưởng của Vật lý Nhiệt

Phân bố công suất là một yếu tố quan trọng cần được xem xét trong thiết kế và vận hành của lò phản ứng VVER. Các mô hình kết hợp vật lý-thủy nhiệt có thể giúp dự đoán chính xác hơn phân bố công suất trong bó nhiên liệu và đánh giá tác động của các thay đổi nhiệt độ và lưu lượng chất làm mát đối với phân bố này.

5.2. Mật độ Chất làm mát yếu tố An toàn quan trọng

Mật độ chất làm mát là một yếu tố quan trọng cần được xem xét trong phân tích an toàn của lò phản ứng VVER. Các mô hình kết hợp vật lý-thủy nhiệt có thể giúp dự đoán chính xác hơn phân bố mật độ chất làm mát trong bó nhiên liệu và đánh giá tác động của các thay đổi mật độ đối với tính ổn định của lò phản ứng.

5.3. Tương tác Vật Lý và Thủy Nhiệt trong Bó Nhiên Liệu

Các mô hình kết hợp vật lý-thủy nhiệt cho phép các kỹ sư mô phỏng đồng thời các quá trình neutron và thủy động lực học trong bó nhiên liệu. Điều này cho phép họ đánh giá tác động của các thay đổi nhiệt độ và lưu lượng chất làm mát đối với phân bố công suất, phân bố mật độ chất làm mát và các thông số quan trọng khác.

VI. Kết Luận và Tương Lai Mô Hình Hóa VVER Toàn Diện

Mô hình Vật lý-Thủy Nhiệt kết hợp cung cấp cách tiếp cận mạnh mẽ để mô hình hóa các thông số của lò phản ứng VVER. Các mô hình có thể được kiểm định chéo và đánh giá độ tin cậy. Mô hình hóa Vật lý và Thủy Nhiệt rất quan trọng để phân tích. Các mô hình có thể được sử dụng cho các hoạt động quản lý lò phản ứng trong tương lai.

6.1. Phát triển mô hình Vật lý và Thủy nhiệt VVER

Phát triển mô hình cho phép tính toán chuyên sâu hơn.

6.2. Ứng dụng và Đánh Giá

Các kết quả đánh giá giúp hiểu chuyên sâu hơn.

Nội dung chính

Tổng quan nghiên cứu

Lò phản ứng VVER-1000 là một trong những loại lò phản ứng hạt nhân thế hệ hiện đại, được sử dụng rộng rãi tại nhiều quốc gia với công suất khoảng 1000 MW điện. Việc tính toán chính xác các đặc tính vật lý và thủy nhiệt của bó nhiên liệu trong lò phản ứng là yếu tố then chốt đảm bảo an toàn và hiệu quả vận hành. Nghiên cứu tập trung xây dựng mô hình tính toán kết hợp vật lý – thủy nhiệt cho bó nhiên liệu lõi phản ứng VVER sử dụng chương trình MCNP5 và COBRA-EN, nhằm mô phỏng phân bố công suất, nhiệt độ và các hiệu ứng phản hồi nhiệt độ trong vùng hoạt của lò. Mục tiêu cụ thể là phát triển mô hình tính toán kết hợp, đánh giá phân bố công suất 3 chiều, phân bố nhiệt độ chất tải nhiệt và nhiệt độ trong thanh nhiên liệu, đồng thời phân tích ảnh hưởng của mật độ chất tải nhiệt lên tính toán vùng hoạt. Nghiên cứu được thực hiện trên mô hình bó nhiên liệu VVER-1000 tại điều kiện hoạt động Hot Full Power (575 K), với dữ liệu benchmark từ các nghiên cứu quốc tế. Kết quả nghiên cứu có ý nghĩa quan trọng trong việc nâng cao độ chính xác của các mô hình tính toán an toàn, hỗ trợ thiết kế và vận hành lò phản ứng hạt nhân hiệu quả hơn.

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 dựa trên hai khung lý thuyết chính: vật lý hạt nhân và thủy nhiệt. Vật lý hạt nhân được mô hình hóa bằng chương trình MCNP5, sử dụng phương pháp Monte Carlo để tính toán phân bố neutron và công suất trong bó nhiên liệu. Thủy nhiệt được mô phỏng bằng chương trình COBRA-EN, mô hình thủy động lực học đa kênh, tính toán phân bố nhiệt độ và các hiệu ứng phản hồi nhiệt độ trong kênh chất tải nhiệt. Các khái niệm chính bao gồm: hệ số nhân vô hạn (k-inf), phân bố công suất 3 chiều, phân bố nhiệt độ chất tải nhiệt, hiệu ứng phản hồi nhiệt độ, và mật độ chất tải nhiệt. Mô hình kết hợp vật lý – thủy nhiệt được xây dựng dựa trên thuật toán vòng lặp, trong đó kết quả phân bố công suất từ MCNP5 được dùng làm dữ liệu đầu vào cho COBRA-EN, và ngược lại, nhiệt độ từ COBRA-EN được cập nhật cho MCNP5 để tính toán lại phân bố neutron, đảm bảo điều kiện hội tụ.

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

Nguồn dữ liệu chính là mô hình bó nhiên liệu VVER-1000 với 66 thanh nhiên liệu, trong đó nhiên liệu UO2 làm giàu 3.7% 235U và U-Gd2O3 làm giàu 3.6% 235U, cùng các thanh điều khiển và thanh trung tâm bằng hợp kim Zr-Nb. Cỡ mẫu mô hình là toàn bộ bó nhiên liệu, mô hình hóa chi tiết 1/6 bó nhiên liệu trong COBRA-EN để giảm thiểu tính toán nhưng vẫn đảm bảo độ chính xác. Phương pháp phân tích sử dụng vòng lặp tính toán kết hợp MCNP5 và COBRA-EN, với thuật toán tự động cập nhật dữ liệu đầu vào và đầu ra giữa hai chương trình cho đến khi hội tụ các tham số như hệ số nhân vô hạn và phân bố nhiệt độ. Timeline nghiên cứu kéo dài trong khoảng thời gian thực hiện mô hình hóa, tính toán benchmark, và phân tích kết quả, dựa trên dữ liệu benchmark quốc tế và các điều kiện vận hành thực tế của lò phản ứng.

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

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

1. Phân bố công suất 3 chiều trong bó nhiên liệu: Kết quả tính toán cho thấy phân bố công suất theo chiều cao và bán kính bó nhiên liệu có sự biến đổi rõ rệt, với công suất tương đối dao động từ 0.5 đến 3.0 m chiều cao. Phân bố công suất theo bán kính thể hiện sự tập trung công suất ở vùng trung tâm, giảm dần về phía biên, với sai số dưới 5% so với benchmark quốc tế.

2. Phân bố nhiệt độ chất tải nhiệt: Nhiệt độ chất tải nhiệt phân bố theo chiều cao bó nhiên liệu dao động từ khoảng 560 K đến 610 K, với nhiệt độ trung bình tại đáy là 560 K và tăng dần lên đỉnh bó nhiên liệu đạt 610 K. Sự phân bố này phù hợp với điều kiện hoạt động Hot Full Power và được cập nhật liên tục trong vòng lặp tính toán.

3. Hệ số nhân vô hạn (k-inf): Hệ số nhân vô hạn được tính toán qua các vòng lặp hội tụ với giá trị ổn định khoảng 1.0001, sai số giữa các vòng lặp giảm dần và đạt mức dưới 0.0001 sau vòng lặp thứ 6, chứng tỏ mô hình tính toán kết hợp có độ chính xác cao và ổn định.

4. Ảnh hưởng của mật độ chất tải nhiệt: Mật độ chất tải nhiệt ảnh hưởng trực tiếp đến phân bố nhiệt độ và hiệu ứng phản hồi nhiệt độ trong vùng hoạt. Mô hình cho thấy khi mật độ chất tải nhiệt giảm, nhiệt độ trong thanh nhiên liệu và chất tải nhiệt tăng lên, làm giảm hệ số nhân hiệu dụng, ảnh hưởng đến an toàn vận hành.

Thảo luận kết quả

Nguyên nhân của sự phân bố công suất và nhiệt độ như trên là do đặc tính vật lý của nhiên liệu và cấu trúc thanh nhiên liệu, cũng như sự tương tác giữa neutron và chất tải nhiệt. So sánh với các nghiên cứu quốc tế, kết quả phân bố công suất và nhiệt độ trong nghiên cứu này có độ lệch nhỏ, chứng tỏ mô hình tính toán kết hợp MCNP5 và COBRA-EN phù hợp và tin cậy. Việc sử dụng vòng lặp cập nhật dữ liệu giữa hai chương trình giúp mô hình phản ánh chính xác các hiệu ứng phản hồi nhiệt độ, điều mà các mô hình tính toán riêng lẻ khó đạt được. Dữ liệu có thể được trình bày qua biểu đồ phân bố công suất theo chiều cao và bán kính (Hình 11, 12), biểu đồ phân bố nhiệt độ chất tải nhiệt theo chiều cao (Hình 13), và bảng hệ số nhân vô hạn theo vòng lặp (Bảng 7, 8) để minh họa sự hội tụ và độ chính xác của mô hình.

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

1. Triển khai mô hình tính toán kết hợp trong thiết kế lò phản ứng: Áp dụng mô hình MCNP5 và COBRA-EN để mô phỏng chính xác phân bố công suất và nhiệt độ, nâng cao độ an toàn và hiệu quả vận hành trong vòng 1-2 năm, do các đơn vị thiết kế và vận hành lò phản ứng thực hiện.

2. Cập nhật dữ liệu vận hành thực tế: Thu thập và cập nhật dữ liệu nhiệt độ, áp suất và mật độ chất tải nhiệt từ các lò phản ứng đang vận hành để hiệu chỉnh mô hình, đảm bảo tính thực tiễn và độ tin cậy, thực hiện định kỳ hàng năm bởi các trung tâm nghiên cứu hạt nhân.

3. Phát triển phần mềm tích hợp tự động: Xây dựng phần mềm tích hợp tự động kết nối MCNP5 và COBRA-EN, giảm thiểu sai sót và tăng tốc độ tính toán, hoàn thành trong 3 năm, do các nhóm phát triển phần mềm chuyên ngành thực hiện.

4. Đào tạo và nâng cao năng lực chuyên môn: Tổ chức các khóa đào tạo về mô hình tính toán kết hợp cho kỹ sư và nhà khoa học hạt nhân, nâng cao năng lực phân tích và vận hành, triển khai liên tục, do các trường đại học và viện nghiên cứu đảm nhiệm.

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

1. Kỹ sư và nhà thiết kế lò phản ứng hạt nhân: Nghiên cứu cung cấp công cụ mô hình hóa chính xác, hỗ trợ thiết kế vùng hoạt và bó nhiên liệu, giúp tối ưu hóa hiệu suất và an toàn.

2. Nhà nghiên cứu vật lý hạt nhân và thủy nhiệt: Luận văn cung cấp phương pháp kết hợp mô hình vật lý và thủy nhiệt, mở rộng kiến thức và ứng dụng trong nghiên cứu chuyên sâu.

3. Cơ quan quản lý an toàn hạt nhân: Tham khảo để đánh giá các mô hình tính toán an toàn, kiểm soát và giám sát vận hành lò phản ứng.

4. Sinh viên và học viên cao học chuyên ngành kỹ thuật hạt nhân: Tài liệu tham khảo quý giá cho việc học tập, nghiên cứu và phát triển các mô hình tính toán trong lĩnh vực hạt nhân.

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

1. Mô hình tính toán kết hợp vật lý – thủy nhiệt là gì?
   Là phương pháp sử dụng đồng thời hai chương trình MCNP5 (vật lý hạt nhân) và COBRA-EN (thủy nhiệt) để mô phỏng phân bố neutron, công suất và nhiệt độ trong bó nhiên liệu, cập nhật dữ liệu qua vòng lặp cho đến khi hội tụ.

2. Tại sao cần mô hình kết hợp thay vì mô hình riêng lẻ?
   Mô hình kết hợp phản ánh chính xác các hiệu ứng phản hồi nhiệt độ và tương tác giữa neutron và chất tải nhiệt, giúp nâng cao độ chính xác và tin cậy của kết quả tính toán.

3. Cỡ mẫu và phạm vi mô hình hóa như thế nào?
   Mô hình hóa toàn bộ bó nhiên liệu VVER-1000 với 66 thanh nhiên liệu, trong COBRA-EN mô hình hóa chi tiết 1/6 bó nhiên liệu để giảm thiểu tính toán nhưng vẫn đảm bảo độ chính xác.

4. Kết quả mô hình có thể ứng dụng thực tế ra sao?
   Kết quả giúp tối ưu hóa thiết kế và vận hành lò phản ứng, nâng cao an toàn, đồng thời hỗ trợ đánh giá và kiểm soát các thông số vận hành quan trọng như phân bố công suất và nhiệt độ.

5. Mô hình có thể áp dụng cho các loại lò phản ứng khác không?
   Phương pháp mô hình kết hợp có thể điều chỉnh và áp dụng cho các loại lò phản ứng khác có cấu trúc tương tự, tuy nhiên cần hiệu chỉnh các tham số đặc thù của từng loại lò.

Kết luận

• Xây dựng thành công mô hình tính toán kết hợp vật lý – thủy nhiệt cho bó nhiên liệu lò phản ứng VVER-1000 sử dụng MCNP5 và COBRA-EN.
• Mô hình cho kết quả phân bố công suất và nhiệt độ chính xác, hội tụ ổn định với sai số nhỏ.
• Phân tích ảnh hưởng mật độ chất tải nhiệt giúp nâng cao hiểu biết về hiệu ứng phản hồi nhiệt độ trong vùng hoạt.
• Mô hình có thể ứng dụng trong thiết kế, vận hành và đánh giá an toàn lò phản ứng hạt nhân.
• Đề xuất phát triển phần mềm tích hợp và đào tạo chuyên môn để ứng dụng rộng rãi trong ngành hạt nhân.

Triển khai mô hình trong các dự án thực tế, cập nhật dữ liệu vận hành, phát triển phần mềm tích hợp tự động.

Các nhà nghiên cứu và kỹ sư trong lĩnh vực kỹ thuật hạt nhân nên áp dụng và phát triển mô hình này để nâng cao hiệu quả và an toàn trong vận hành lò phản ứng.

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

Chương 1. Lý do lựa chọn đề tài  lý                   H    -  (Neutronics-Thermal hydraulics)         -   àn (safety margins)  t vài         [2]                          -           Zwermann und W. Verdú [4]  -TF (CTF) cho các tính toán vùng                     8             Tuy nhiên, các tính to                                  lý    1. Mục tiêu nghiên cứu.

Nội dung nghiên cứu. t hp khi chy vi mt tính toán vt lý b t qu phân b công sut 3 chiu trong bó nhiên liu vi u kin gi s u là phân b nhi cht ti nhing nht trong toàn bó nhiên liu. Phân b công suc t c s dng làm d liu vào cho tính toán thy nhit v-EN. y nhit qu phân b 3 chiu m cht ti nhit trong bó nhiên liu, kt qu này tip tc s dng làm d liu u vào cho tính toán vt lý b vòng lp tip theo.

c thc hiu kin hi t c thõa mãn.   10       trình PERL.           -  cho các tính toán cháy và tính toán phân tích an toàn. - : Xây dng mô hình tính toán kt hp Vt lý  Thy nhit.

- : Các kt qu tính toán.     mô hình tính toán     . Tổng quan về lò phản ứng VVER-1000    -             -     -   Hình 1       - 1000 Bushehr. Hình 1: Cấu hình vùng hoạt lò phản ứng VVER-1000 (Bushehr) 12      -          bó nhiên -1000/V392 là 353 cm.

Hình 2 mô   -      Hình 2: Bó nhiên liệu lò phản ứng VVER-1000.    - Các thông s thit k c trình bày trong ni dung bài toán benchmark. - c nhiu nhóm nghiên cu trên th gii tính toán li và kim tra kt qut qu  tin cy cao và là   khnh tính chính xác cho các tính toán ca tác gi. Tổng quan về các hiệu ứng phản hồi nhiệt độ                     -238 do  2                   -             14                               Hình 3 .60 550 560 570 580 590 600 610 Nhiệt độ (K) Hình 3: Mối tương quan giữa nhiệt độ và mật độ nước tại áp suất 15.7 Mpa                   VVER-                -  15    -EN.

Ảnh hưởng của mật độ chất tải nhiệt lên tính toán vùng hoạt lò phản ứng.                                 i i+1        : -   -    -    -    i vi vùng hot ca lò phn c hu hn, ta cn quan tâm n h s nhân hiu dng vì s phát trin ca phn ng dây chuyn ph thuc vào h s nhân hiu dng k eff. P (2) vi P < 1 là xác sung hu hn. 16 i vi lò nhic gii hn thì có hi trình làm cht do quá trình khu s nhân tr thành:    (3)  L     - Xác sut tránh hp th cng: 1s 2 p (4)  a1   1s2 :  s   s.

Xây dựng mô hình tính toán kết hợp Vật lý – Thủy nhiệt 3. Mô tả các chương trình tính toán 3. Chương trình tính toán MCNP5 -                             -            -                    nh tính toán MCNP5   -  chính xác ca kt qu cao. - c kinh bi nhiu nhóm nghiên cu ln trên th gii, u công b trên các tp chí khoa hc uy tín s d tính toán này.

Chương trình tính toán COBRA-EN - -  COBRA-         -3C/MIT cho các phân              -            .   - momen cho h           (Electric Power Research Institute            (Simplified Boiling Water Reactor   và AP600  trong các phân      -           22 - S d thích hp cho các tính toán kênh cht ti nhit trong mô hình bó nhiên liu. - c s dng bi mt s nhóm nghiên cu trên th ging công b s dp chí khoa hc uy tín.      -                    -    thành các nhóm: - Các th n hình hc ca mô hình.

- Các th n thông s công sut ca các thanh hay bó nhiên liu. - Các th   n các thông s thy nhit ca mô hình (áp sut, nhi   - Các th n các hiu chnh vi mô hình tính toán. -    Hình 4 d  -                     Bảng 1 23                   COBRA-EN. Hình 4 : Mô hình hóa trong chương trình tính toán COBRA-EN.

Bảng 1: Một số đặc trưng của thanh nhiên liệu và kênh chất tải nhiệt         Chu vi   Gap     24 -  - Phân b nhi cht ti nhit. - Phân b void fraction (pht ti nhit). - Nhi nhiên liu - H s DNB (DNBR). Mô hình bó nhiên liệu lò phản ứng VVER-1000   benchmark  -1000 [1].

Các thanh nhiên liu UO2   làm giàu 3.7 wt% 235U, các thanh nhiên liu U- làm giàu 3.6 wt% 235U và cha 4. V bc các thanh nhiên liu và các ng du khin, thanh trung tâm s dng vt liu là hp kim Zr-Nb. M hng v thành phn trong nhiên liu và các vt liu c Bảng 2. 25 Hình 5: Tiết diện mặt cắt ngang bó nhiên liệu lò VVER-1000.

Hình 6: Cấu hình mặt cắt ngang ô mạng các thanh nhiên liệu và các ống dẫn 26 Bảng 2: Mật độ hạt nhân các đồng vị thành phần trong nhiên liệu và vật liệu cấu trúc.   UO2 UO2-Gd2O3   U235 8.94E-05 27 Bảng 3: Kích thước vỏ thanh nhiên liệu và ống dẫn thanh điều khiển, thanh trung tâm Bán kính trong (R1) (cm) Bán kính ngoài (R2) (cm)  0. Mô hình hóa bó nhiên liệu lò phản ứng VVER-1000 với chương trình MCNP5 và COBRA-EN 3. Mô hình hóa bó nhiên liệu với chương trình MCNP5        benchmark  -1000 [1]  2O3),                   trong Hình 5.

Các thanh nhiên liu UO 2  làm giàu 3.7 wt% 235U, các thanh nhiên liu U-Gd  làm giàu 3.6 wt% 235 U và cha 4. V bc các thanh nhiên liu và các ng du khin, thanh trung tâm s dng vt liu là hp kim Zr-Nb. Mt  hng v thành phn trong nhiên liu và các vt liu c ra trong Bảng 3. 28      -  th  29         tính toán  MeV/g.

Hình 7: Mô hình bó nhiên liệu lò phản ứng VVER-1000 với chương trình MCNP5 (a – Tiết diện theo chiều cao; b – Tiết diện mắt cắt ngang bó nhiên liệu) 31 3. Mô hình COBRA-EN  - -ình tính  Hình 8. Hình 8: Mô hình hóa 1/6 bó nhiên liệu lò phản ứng VVER-1000 với chương trình COBRA-EN -EN s dng mô hình 1/6 bó nhiên liu lò phn ng VVER-1000 vi 66 thanh nhiên liu    a mô hình tính toán trong COBRA-EN nên 1/6 bó nhiên lic mô hình hóa các thanh nhiên liu ti c hai biên ci xng) và các kênh cht ti nhit bao quanh thanh nhiên li trình tính toán này không s dng các d liu hình h   , hình dng ca kênh cht ti nhit trong vic xây dng mô hình. Vic xây d  c thc 32 hi-EN s dng các d liu v a kênh cht ti nhit (din tích kênh cht ti nhit, chu vi kênh cht ti nhi và ma các kênh cht ti nhit lin k.

-c mô hình  - M  công sut các thanh (nhiên liu và các thanh cu trúc): Trong th này  chiu cao ca mi lp (bó nhiên lic chia thành 10 lu nhau theo chiu cao). M công sut ti mi thanh nhiên lic tính toán t mt  công sut trung bình ca bó nhiên liu và phân b công sut ca bó nhiên liu. - Mi liên kt gia các kênh cht ti nhit lin k: Trong th này cn khai báo mi liên h gia các kênh cht ti nhit lin k vt này là: S kênh, chu vi phn tip xúc gia các kênh và khong cách gia tâm ca 2 kênh cht ti nhit lin k.