Đồ án tập trung vào tìm hiểu phân bố nhiệt độ trong một lò thực tế trong quátrình hoạt động 100% công suất. Nội dung chính bao hàm trong ba chƣơng.Chƣơng I chú trọng về giới thiệu hai chƣơng trình cơ bản đƣợc sử dụng chonghiên cứu. Trong đó đề cập đến lý thuyết và các khai báo đã sử dụng ở mỗi chƣơngtrình để giải quyết bài toán.Chƣơng II cung cấp các thông số của lò nƣớc áp lực Tomari liên quan đến tínhtoán cũng nhƣ phƣơng pháp nghiên cứu để giải quyết bài toán.Chƣơng III cung cấp kết quả các bài toán và thảo luận.
ĐỒ ÁN TỐT NGHIỆP LỜI CẢM ƠN Tôi xin gửi lời cảm ơn chân thành đến ThS Trần Thùy Dƣơng Viện Kỹ thuật Hạt nhân Vật lý Môi trƣờng Cơ nhiệt tình bảo giúp đỡ tơi hồn thành đồ án Tơi xin gửi lời cảm ơn sâu sắc đến thầy cô viện Kỹ thuật Hạt nhân Vật lý Môi trƣờng tạo điều kiện giúp đỡ q trình tơi hoàn thành đồ án Đặc biệt TS Nguyễn Văn Thái, phó Viện trƣởng Viện Kỹ thuật Hạt nhân Vật lý Mơi trƣờng Thầy có góp ý tận tình giúp đỡ tơi q trình hồn thành đồ án Đồng thời, tơi xin cảm ơn ThS Trịnh Hữu Toản Thầy giúp đỡ để tơi hồn thành đồ án cách tốt Tơi xin gửi lời tri ân đến gia đình, ngƣời thân ln cổ vũ, khuyến khích tơi vƣơn lên suốt thời gian qua Cuối cùng, xin cảm ơn bạn bè lớp KTHN&VLMT K54 động viên, chia sẻ, giúp đỡ lúc gặp khó khăn q trình làm đồ án Sinh viên Nguyễn Hải Thìn Nguyễn Hải Thìn Lớp KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP TÓM TẮT NỘI DUNG ĐỒ ÁN Đồ án tập trung vào tìm hiểu phân bố nhiệt độ lò thực tế q trình hoạt động 100% cơng suất Nội dung bao hàm ba chƣơng Chƣơng I trọng giới thiệu hai chƣơng trình đƣợc sử dụng cho nghiên cứu Trong đề cập đến lý thuyết khai báo sử dụng chƣơng trình để giải tốn Chƣơng II cung cấp thơng số lò nƣớc áp lực Tomari liên quan đến tính tốn nhƣ phƣơng pháp nghiên cứu để giải toán Chƣơng III cung cấp kết toán thảo luận Nguyễn Hải Thìn Lớp KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP MỤC LỤC MỞ ĐẦU Chƣơng I TỔNG QUAN VỀ MCNP5 RELAP5 2 Tổng quan chung MCNP5 1.1 Giới thiệu chung 1.2 Dữ liệu hạt nhân phản ứng 1.3 Các khai báo cần thiết để giải toán Giới thiệu RELAP5 13 2.1 Cấu trúc RELAP5 13 2.2 Cơ sở lý thuyết 14 2.3 Sử dụng RELAP5 để giải toán 14 Chƣơng II THÔNG SỐ LỊ TOMARI, PHƢƠNG PHÁP NGHIÊN CỨU VÀ CÁC BÀI TỐN MÔ PHỎNG 19 Một số thông số lò Tomari 19 Phƣơng pháp nghiên cứu 30 Các tốn mơ 31 Chƣơng III KẾT QUẢ VÀ THẢO LUẬN 33 Bài toán 1: Bài toán sở 33 Bài toán 2: Phân bố nhiệt độ lò trƣờng hợp nguồn tâm lò trạng thái tới hạn, công suất 100% 37 Bài toán 3: Phân bố nhiệt độ chất làm mát lò trƣờng hợp nguồn vị trí đặt nguồn neutron ban đầu trạng thái tới hạn 40 KẾT LUẬN 46 TÀI LIỆU THAM KHẢO 47 PHỤ LỤC 48 Dữ liệu đầu vào cho MCNP5 48 Dữ liệu đầu vào cho RELAP5 79 Nguyễn Hải Thìn Lớp KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP Danh mục bảng biểu Bảng 1 Định nghĩa bề mặt mà MCNP sử dụng Bảng Các tally sử dụng MCNP5 11 Bảng Phân bố nhiên liệu đầu chu trình 24 Bảng 2 Các thơng số lò Tomari 26 Bảng Một số thành phần vật liệu 29 Bảng Khối lƣợng riêng chất sử dụng khai báo 30 Bảng Số lịch sử vòng 150 chu kỳ giá trị keff tƣơng ứng 34 Bảng Giá trị độ phản ứng so sánh với trƣờng hợp lò có chung nhiệt độ 40 Bảng 3 Giá trị độ phản ứng so sánh với trƣờng hợp lò có chung nhiệt độ 45 Nguyễn Hải Thìn Lớp KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP Danh mục hình vẽ đồ thị Hình 1 Lịch sử ngẫu nhiên neutron Hình Cấu trúc top-down chƣơng trình RELAP5 13 Hình Cấu trúc nhiên liệu 19 Hình 2 Ma trận nhiên liệu bó nhiên liệu 20 Hình Cấu tạo bó điều khiển 21 Hình Cấu tạo bó chất độc 22 Hình Cấu tạo nguồn ban đầu 23 Hình Ma trận cho nhiên liệu (U,Gd)O2 24 Hình Phân bố nhiên liệu ban đầu khởi động lò chu trình 25 Hình Mặt cắt ngang vùng hoạt lò phản ứng hạt nhân 26 Hình Sơ đồ thuật tốn mở rộng 31 Hình Lõi lò Tomari đƣợc mơ 33 Hình Nút hóa khu vực lõi lò Tomari 34 Hình 3 Kết khảo sát sơ nhiệt độ nhiên liệu 40 lần (lần lƣợt từ trái sang phải từ xuống dƣới) 35 Hình Phân bố nhiệt độ chất làm chậm 10 lần lặp 36 Hình Phân bố nhiệt độ chất làm chậm 37 Hình Phân bố nhiệt độ nhiên liệu 37 Hình Phân bố nhiệt độ nhiên liệu 38 Hình Phân bố nhiệt độ chất làm chậm 38 Hình Phân bố cơng suất mơ hình lò Tomari 39 Hình 10 Phân bố nhiệt độ chất làm chậm 41 Hình 11 Phân bố nhiệt độ nhiên liệu 41 Hình 12 Phân bố nhiệt độ nhiên liệu 42 Hình 13 Phân bố nhiệt độ chất làm chậm 42 Hình 14 Phân bố cơng suất mơ hình lò Tomari 43 Hình 15 Phân bố nhiệt độ theo chiều cao vùng hoạt 44 Hình 16 Phân bố cơng suất lò thực tế 44 Nguyễn Hải Thìn Lớp KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP MỞ ĐẦU Hiện nay, khủng hoảng lƣợng vấn đề trọng tâm với quốc gia, đặc biệt quốc gia phát triển nhƣ Việt Nam Tuy nhiên, việc phụ thuộc vào nguồn lƣợng hóa thạch hay thủy điện khơng gây ảnh hƣởng đến mơi sinh mà tạo nên nguy an ninh lƣợng cạn kiệt chúng Còn số giải pháp lƣợng thân thiện với môi trƣờng bền vững nhƣ lƣợng mặt trời, lƣợng gió,… ngồi tầm với chi phí đầu tƣ vận hành lớn Chính thế, điện hạt nhân trở thành giải pháp sáng giá khả dung hòa yêu cầu môi trƣờng, kinh tế lƣợng vận hành an tồn Do đó, u cầu thiết kỹ sƣ hạt nhân phải tính tốn đƣợc thơng số lò để kiểm tra đảm bảo hoạt động an toàn lò Trong đó, hệ số phản hồi âm có ảnh hƣởng nhiều tới kết xác tính tốn cho lò nhƣ hiệu ứng Doppler, hệ số chất làm chậm,… Trong phản hồi âm này, tính tốn nhiệt độ giải cách kết hợp chƣơng trình RELAP5 MCNP5 Do đó, tơi chọn đề tài nhằm tìm hiểu phân bố nhiệt độ và ảnh hƣởng đến kết tính tốn lò Mục đích nghiên cứu: Đồ án ƣu tiên việc tìm đƣợc phân bố nhiệt độ chất làm mát nhiên liệu lò Đối tƣợng nghiên cứu: Đối tƣợng đƣợc nghiên cứu vùng hoạt lò Tomari số Phạm vi nghiên cứu: Đồ án giới hạn nghiên cứu khn khổ tìm hiểu riêng phân bố nhiệt độ lò đầu chu kỳ thứ lò Tomari Phƣơng pháp nghiên cứu: đồ án kết hợp hai chƣơng trình RELAP5 MCNP5 Nguyễn Hải Thìn KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP Chƣơng I TỔNG QUAN VỀ MCNP5 RELAP5 Tổng quan chung MCNP5 1.1 Giới thiệu chung MCNP chƣơng trình tính tốn trình vận chuyển Monte Carlo đƣợc phát triển phòng thí nghiệm quốc gia Los Alamos, Mỹ Nó đƣợc sử dụng cho vài chế độ vận chuyển: neutron, photon, electron, kết hợp vận chuyển neutron/photon photon đƣợc sinh tƣơng tác neutron, neutron/photon/electron, photon/electron hay electron/photon Dải lƣợng neutron từ 10-11 MeV đến 20 MeV cho toàn đồng vị tới 150 MeV cho số đồng vị, dải lƣợng cho photon từ keV tới 100 GeV dải lƣợng electron từ 1keV tới GeV Khả tính tốn giá trị riêng keff cho hệ thống phân hạch tính chủ chốt Ngƣời dùng tạo tệp đầu vào đƣợc đọc sau MCNP Tệp bao gồm thông tin vấn đề mặt nhƣ: hình học, vật liệu, vị trí tính chất nguồn, tally mong muốn,… Monte Carlo nhân lý thuyết trình thống kê (nhƣ tƣơng tác hạt nhân với vật liệu) đặc biệt hữu ích trƣờng hợp hình học phức tạp Các kiện xác suất q trình đƣợc mơ Các phân bố xác suất kiểm soát kiện đƣợc mẫu hóa thống kê để mơ tả tồn tƣợng nói chung, mơ đƣợc biểu diễn máy tính số lần thử cần thiết để mô tả đầy đủ tƣợng thƣờng lớn Q trình mẫu hóa thống kê dựa lựa chọn số ngẫu nhiên – tƣơng tự nhƣ nhƣ ném xúc xắc – có tên “Monte Carlo” Trong vận chuyển hạt, kỹ thuật Monte Carlo có tính thực ƣu việt Nó theo dõi hạt từ lúc hạt sinh đến theo số phƣơng cách (hấp thụ, rò,…) Các phân bố thống kê đƣợc mẫu hóa ngẫu nhiên sử dụng liệu vận chuyển để xác định kết bƣớc thời gian Nguyễn Hải Thìn KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP Hình 1 Lịch sử ngẫu nhiên neutron Hình 1.1 đại điện cho lịch sử ngẫu nhiên hạt neutron tới phẳng vật liệu diễn q trình phân hạch Các số nằm đƣợc lựa chọn ngẫu nhiên để xác định tƣơng tác (nếu có) đâu xảy dựa quy luật (vật lý) thống kê (dữ liệu vận chuyển) kiểm soát q trình vật liệu liên quan Trong ví dụ này, hạt neutron va chạm xảy kiện Hạt neutron bị tán xạ theo nhƣ hƣớng hình vẽ, hƣớng đƣợc lựa chọn ngẫu nhiên từ phân bố tán xạ thực tế Một hạt photon đƣợc tạo tạm thời lƣu giữ “đóng băng” để phân tích sau Ở kiện 2, phân hạch xảy ra, dẫn đến biến neutron đến sinh hai neutron bay photon Một neutron photon đƣợc “đóng băng”để phân tích sau Hạt neutron sinh phân hạch bị bắt kiện biến Hạt neutron đƣợc “đóng băng” lúc phân hạch đƣợc phục hồi mẫu hóa ngẫu nhiên, rò khỏi mặt phẳng kiện Hạt photon sản phẩm phân hạch có va chạm kiện rò ngồi kiện Hạt photon lại đƣợc tạo kiện đƣợc theo dõi bị bắt kiện Chú ý MCNP khơi phục hạt bị “đóng băng” theo ngun tắc hạt lƣu trữ trạng thái “đóng băng” đƣợc đƣa Nguyễn Hải Thìn KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP Lịch sử neutron hoàn thành Càng nhiều lịch sử đƣợc theo dõi, phân bố neutron photon gần với dạng biết Khía cạnh mà ngƣời dùng quan tâm đƣợc đƣa tally kèm theo ƣớc tính sai số thống kê kết 1.2 Dữ liệu hạt nhân phản ứng MCNP sử dụng thƣ viện liệu nguyên tử hạt nhân lƣợng liên tục Dữ liệu đƣợc đánh giá đƣợc chuyển thành định dạng xấp xỉ cho MCNP mã nhƣ NJOY Các bảng số liệu hạt nhân có tƣơng tác neutron, tƣơng tác photon photon (đƣợc tạo neutron), phép đo liều hay kích hoạt neutron mơi trƣờng có tính đến chuyển động nhiệt phân tử Mỗi bảng số liệu sẵn có MCNP đƣợc lập danh sách tệp thƣ mục XSDIR Ngƣời sử dụng lựa chọn bảng số liệu hạt nhân qua ký hiệu nhận dạng ZAID bảng Nhìn chung, ký hiệu nhận dạng chứa số nguyên tử Z, số khối A thƣ viện riêng ID Trên 836 bảng tƣơng tác hạt nhân đƣợc sử dụng cho xấp xỉ 100 nguyên tố khác Có nhiều bảng dùng cho đồng vị đơn lẻ liệu đƣợc lấy từ giá trị khác nhau, khoảng nhiệt độ khác dung sai xử lý khác Dữ liệu phản ứng photon ( đƣợc tạo từ neutron có tƣơng tác) đƣợc đƣa vào nhƣ phần bảng tƣơng tác neutron Bảng tƣơng tác photon có cho tất nguyên tố có Z =1 tới Z=100 Dữ liệu bảng tƣơng tác photon cho phép MCNP đếm số tán xạ đàn hồi, không đàn hồi, hấp thụ quang điện với xác suất phát huỳnh quang tạo cặp.Tiết diện phản ứng cho gần 2000 phản ứng liên quan tới 400 hạt nhân bia trạng thái bền kích thích phần số liệu hạt nhân MCNP Các đơn vị đƣợc sử dụng chƣơng trình MCNP: Độ dài dùng centimeters (cm) Nguyễn Hải Thìn KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP Năng lƣợng theo MeV, Thời gian theo “shake” (1 shake = 10x10-9 s = 10-8 s ) Nhiệt độ theo MeV (kT), Mật độ nguyên tử có đơn vị: nguyên tử/barn-cm Khối lƣợng riêng theo g/cm3 Tiết diện phản ứng theo barns (10-24 cm2 ), Suất sinh nhiệt (heating numbers) theo MeV/va chạm Tỷ lệ khối lƣợng nguyên tử dựa khối lƣợng neutron 1.008664967 amu Trong đơn vị này, số Avogadro 0.59703109 x 10+24 1.3 Các khai báo cần thiết để giải tốn 1.3.1 Khai báo thẻ (cell cards) Các thẻ ô mục đứng sau tiêu đề Khơng có dòng trắng phân cách phần tiêu đề phần thẻ ô Trong phần đƣợc định nghĩa hình dạng vật liệu đƣợc điền Mẫu chung thẻ có dạng nhƣ sau: j j m d LIKE n g ds BUT ds Trong j giá trị nguyên có giá trị từ đến 99999 dùng để đặt tên cho ô m số vật liệu, xác định ô đƣợc điền đầy vật liệu thẻ khai báo vật liệu Mm Khi m=0 có nghĩa trống d khối lƣợng riêng vật liệu ô Giá trị dƣơng có nghĩa d đƣợc lấy theo mật độ nguyên tử theo đơn vị 1024 nguyên tử/cm3 Giá trị âm có nghĩa d đƣợc lấy theo khối lƣợng riêng với đơn vị g/cm3 Nguyễn Hải Thìn KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP * flag liq.mass flow vap.mass flow inter.veloc 3360200 3360201 -1.0 11791.38889 0.0 0.0 3360202 0.0 11791.38889 0.0 0.0 3360203 1.0e10 11791.38889 0.0 0.0 * ********************************************************************** ** * * upper plenum * 3400000 upvol1 tmdpvol 3400101 10.498 0.0 3400102 0.686 0.0 8.197 0.0 0.509016 90.0 0.0 3400200 3400201 0.0 3400202 1.0e+3 3400203 1.0e+10 15.41e+6 598.15 15.41e+6 15.41e+6 598.15 598.15 * * * * ********************************************************************** ** * Nguyễn Hải Thìn KTHN&VLMT K54 * 88 ĐỒ ÁN TỐT NGHIỆP * * * vessel heat slabs * * * * * ********************************************************************** ** * thermal shield * 13150000 12 13150100 13150101 1.9 13150201 13150301 0.0 2 1.85 13150401 549.8 13150501 315010000 10000 1 0.366 10 13150502 315020000 1 0.420 11 13150503 315030000 1 1.000 12 13150601 315010000 10000 1 0.366 10 13150602 315020000 1 0.420 11 13150603 315030000 1 1.000 12 13150701 0.0 0.0 0.0 12 13150801 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 12 13150901 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 12 * Nguyễn Hải Thìn KTHN&VLMT K54 89 ĐỒ ÁN TỐT NGHIỆP * downcomer vessel wall * * 13160000 12 13160100 13160101 2.1 13160201 13160301 0.0 13160401 549.8 2.0 13160501 315010000 10000 1 0.366 10 13160502 315110000 1 0.420 11 13160503 315120000 1 1.000 12 13160601 0 0.366 13160602 0 0.366 13160603 0 0.366 13160604 0 0.366 13160605 0 0.366 13160606 0 0.366 13160607 0 0.366 13160608 0 0.366 13160609 0 0.366 13160610 0 0.366 10 13160611 0 0.420 11 13160612 0 1.000 12 Nguyễn Hải Thìn KTHN&VLMT K54 90 ĐỒ ÁN TỐT NGHIỆP 13160701 0.0 0.0 0.0 12 13160801 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 12 * * barrel baffel region - heat slabs between barrel baffel and downcomer * 13200000 10 13200100 13200101 1.8 13200201 13200301 0.0 13200401 561.25 13200501 320010000 10000 1 0.366 10 13200601 315010000 10000 1 0.366 10 13200701 0.0 10 0.0 0.0 1.75 13200801 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 10 13200901 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 10 * * lower plenum - hemispherical shell,lower plenum wall, * core barrel wall,internals * * hemispherical shell * 13230000 13230100 2.0 Nguyễn Hải Thìn KTHN&VLMT K54 91 ĐỒ ÁN TỐT NGHIỆP 13230101 2.1 13230201 13230301 0.0 13230401 549.8 13230501 323010000 1 13230601 0 0.5 13230701 0.0 0.0 0.0 0.5 1 13230801 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 * * lower plenum wall * 13220000 2.0 13220100 13220101 2.1 13220201 13220301 0.0 13220401 549.8 13220501 322010000 1 13220601 0 1.81537 13220701 0.0 0.0 0.0 1.81537 1 13220801 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 * * * core barrel wall Nguyễn Hải Thìn KTHN&VLMT K54 92 ĐỒ ÁN TỐT NGHIỆP * 13250000 13250100 13250101 1.8 13250201 13250301 0.0 13250401 549.8 13250501 325010000 1 1.000 13250502 330010000 1 0.420 13250601 315120000 1 1.000 13250602 315110000 1 0.420 13250701 0.0 1.75 0.0 0.0 13250801 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 13250901 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 * * internals * 13270000 13270100 13270101 0.076029 13270201 13270301 0.0 13270401 549.8 13270501 0.0 20.4 0 Nguyễn Hải Thìn KTHN&VLMT K54 93 ĐỒ ÁN TỐT NGHIỆP 13270502 0 0 20.4 13270503 0 0 40.8 13270504 0 0 18.2 13270601 323010000 20.4 13270602 322010000 20.4 13270603 325010000 40.8 13270604 330010000 18.2 13270701 0.0 0.0 0.0 13270801 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 13270901 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 * * core baffle * 13280000 10 1.52 13280100 13280101 1.54 13280201 13280301 0.0 13280401 561.25 13280501 335010000 10000 13280601 320100000 1 0.366 13280602 320090000 1 0.366 13280603 320080000 1 0.366 13280604 320070000 1 0.366 Nguyễn Hải Thìn KTHN&VLMT K54 0.366 10 94 ĐỒ ÁN TỐT NGHIỆP 13280605 320060000 1 0.366 13280606 320050000 1 0.366 13280607 320040000 1 0.366 13280608 320030000 1 0.366 13280609 320020000 1 0.366 13280610 320010000 1 0.366 10 13280701 0.0 0.0 0.0 10 13280801 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 10 13280901 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 10 * * core fuel pin * 13360000 10 17 13360100 13360101 0.004095 13360102 0.00418 13360103 0.00475 13360201 13360202 -2 13360203 -3 16 13360301 1.0 13360302 0.0 16 13360401 866.5 13360402 755.4 0.0 Nguyễn Hải Thìn KTHN&VLMT K54 95 ĐỒ ÁN TỐT NGHIỆP 13360403 599.8 13360404 588.7 13360405 572 17 13360501 0 0 0.0 13360502 0 0 0.0 13360503 0 0 0.0 13360504 0 0 0.0 13360505 0 0 0.0 13360506 0 0 0.0 13360507 0 0 0.0 13360508 0 0 0.0 13360509 0 0 0.0 13360510 0 0 0.0 10 13360601 335010000 450.59 13360602 335020000 450.59 13360603 335030000 450.59 13360604 335040000 450.59 13360605 335050000 450.59 13360606 335060000 450.59 13360607 335070000 450.59 13360608 335080000 450.59 13360609 335090000 450.59 13360610 335100000 450.59 10 13360701 109 0.019445503 0.0 Nguyễn Hải Thìn KTHN&VLMT K54 0.003700 96 ĐỒ ÁN TỐT NGHIỆP 13360702 109 0.048320083 0.0 0.004823 13360703 109 0.088489085 0.0 0.005638 13360704 109 0.141790436 0.0 0.006144 13360705 109 0.195675335 0.0 0.006342 13360706 109 0.212708051 0.0 0.006231 13360707 109 0.174022014 0.0 0.005811 13360708 109 0.088800745 0.0 0.005083 13360709 109 0.025274471 0.0 0.004046 13360710 109 0.005474278 0.0 0.002700 10 13360801 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 10 13360901 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0 10 * * ********************************************************************** ** * * * heat structure thermal property data * * * ********************************************************************** * * * * composition type and data format * * * ********************************************************************** ** Nguyễn Hải Thìn KTHN&VLMT K54 97 ĐỒ ÁN TỐT NGHIỆP * * ************ Material property data **************** * - material property data * * composition (fuel uo2 pellet) 20100100 tbl/fctn 1 * thermal conductivity table * temperature thermal conduct * (k) (watts/m-k) 20100101 273.150 1.423079e+01 20100102 360.150 1.423079e+01 20100103 440.150 1.323077e+01 20100104 500.150 5.577050e+00 20100105 533.150 5.159348e+00 20100106 616.483 3.999700e+00 20100107 699.817 3.387041e+00 20100108 866.483 2.634192e+00 20100109 1366.483 2.090755e+00 20100110 1616.483 2.199791e+00 20100111 2255.372 2.319544e+00 20100112 2533.150 2.875816e+00 20100113 2810.928 3.409533e+00 20100114 4810.928 3.409533e+00 * volumetric heat capacity table Nguyễn Hải Thìn KTHN&VLMT K54 98 ĐỒ ÁN TỐT NGHIỆP * temperature heat capacity * (k) (j/m3-k) 20100151 273.150 2.042163e+06 20100152 323.150 2.303721e+06 20100153 473.150 2.987795e+06 20100154 673.150 3.064921e+06 20100155 1373.150 2.773183e+06 20100156 1973.150 4.463249e+06 20100157 2173.150 5.569840e+06 20100158 3113.150 6.465172e+06 20100159 4699.817 8.812486e+06 * * composition (fuel gap) 20100200 tbl/fctn 1 * thermal conductivity table * temperature thermal conduct * (k) (watts/m-k) 20100201 273.150 1.525747e+00 20100202 3255.372 1.525747e+00 * volumetric heat capacity table * temperature * (k) heat capacity (j/m3-k) 20100251 273.150 4.359297e+00 20100252 3255.372 4.359297e+00 Nguyễn Hải Thìn KTHN&VLMT K54 99 ĐỒ ÁN TỐT NGHIỆP * * composition (fuel cladding) 20100300 tbl/fctn 1 * thermal conductivity table * temperature thermal conduct * (k) (watts/m-k) 20100301 273.150 1.823523e+01 20100302 473.150 1.823523e+01 20100303 873.150 2.946908e+01 20100304 1073.150 3.146974e+01 20100305 1673.150 3.624304e+01 20100306 2473.150 5.664778e+01 * volumetric heat capacity table * * temperature (k) heat capacity (j/m3-k) 20100351 255.372 1.770009e+06 20100352 1077.761 2.379237e+06 20100353 1185.928 5.041761e+06 20100354 1248.428 2.975723e+06 20100355 2199.817 1.641510e+06 * * composition (inconel) 20100400 tbl/fctn 1 * thermal conductivity table Nguyễn Hải Thìn KTHN&VLMT K54 100 ĐỒ ÁN TỐT NGHIỆP * temperature * (k) thermal conduct (watts/m-k) 20100401 273.1500 1.318841e+01 20100402 838.7056 2.516807e+01 * volumetric heat capacity table * temperature * heat capacity (k) (j/m3-k) 20100451 273.1500 3.834840e+06 20100452 310.9278 3.834840e+06 20100453 477.5944 4.100422e+06 20100454 588.7056 4.276805e+06 20100455 699.8167 4.453860e+06 * * composition (stainless steel) 20100500 tbl/fctn 1 * thermal conductivity table * temperature * (k) thermal conduct (watts/m-k) 20100501 273.1500 1.298030e+01 20100502 1199.817 2.510576e+01 * volumetric heat capacity table * temperature * 20100551 (k) heat capacity (j/m3-k) 273.150 3.830413e+06 Nguyễn Hải Thìn KTHN&VLMT K54 101 ĐỒ ÁN TỐT NGHIỆP 20100552 366.483 3.830413e+06 20100553 422.039 3.964814e+06 20100554 477.594 4.099214e+06 20100555 533.150 4.233615e+06 20100556 588.706 4.334415e+06 20100557 644.261 4.435081e+06 20100558 699.817 4.502416e+06 20100559 810.928 4.636816e+06 20100560 1366.483 5.376019e+06 * ***************************************************** 20210900 power 20210901 0.0 20210902 1.0e+3 2600.e+6 20210903 1.0e+6 2600.e+6 2600.e+6 Nguyễn Hải Thìn KTHN&VLMT K54 102 ... tả lần ô bề mặt cấu trúc xuất lần hình học tốn Nguyễn Hải Thìn KTHN&VLMT K54 11 ĐỒ ÁN TỐT NGHIỆP Nhờ đó, khối lƣợng liệu đầu vào mà ngƣời dùng phải cung cấp lƣợng nhớ máy tính cần thiết để chạy... kiện biên bên phải 2.3.2.4 Các thẻ khai báo bảng chung Nhóm bao gồm bảng cung cấp tính chất nhiệt (độ dẫn nhiệt nhiệt dung) vật liệu bảng cơng suất lò theo thời gian Nguyễn Hải Thìn KTHN&VLMT... chung nhiệt độ 40 Bảng 3 Giá trị độ phản ứng so sánh với trƣờng hợp lò có chung nhiệt độ 45 Nguyễn Hải Thìn Lớp KTHN&VLMT K54 ĐỒ ÁN TỐT NGHIỆP Danh mục hình vẽ đồ thị Hình 1 Lịch sử ngẫu nhiên