luận văn thạc sĩ tính toán phổ năng lượng của kênh nơtron phin lọc từ lò phản ứng hạt nhân đà lạt

BỘ GIÁO DỤC VÀ ĐÀO TẠO VIỆN HÀN LÂM KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM HỌC VIỆN KHOA HỌC VÀ CÔNG NGHỆ - Huỳnh Thanh Sơn TÊN ĐỀ TÀI LUẬN VĂN TÍNH TỐN PHỔ NĂNG LƯỢNG CỦA KÊNH NƠTRON PHIN LỌC TỪ LÒ PHẢN ỨNG HẠT NHÂN ĐÀ LẠT LUẬN VĂN THẠC SĨ VẬT LÝ Khánh Hòa – 2020 BỘ GIÁO DỤC VÀ ĐÀO TẠO VIỆN HÀN LÂM KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM HỌC VIỆN KHOA HỌC VÀ CÔNG NGHỆ - Huỳnh Thanh Sơn TÍNH TỐN PHỔ NĂNG LƯỢNG CỦA KÊNH NƠTRON PHIN LỌC TỪ LÒ PHẢN ỨNG HẠT NHÂN ĐÀ LẠT Chuyên ngành: Vật lý kỹ thuật Mã số: 8520401 LUẬN VĂN THẠC SĨ VẬT LÝ NGƯỜI HƯỚNG DẪN KHOA HỌC: Hướng dẫn 1: PGS.TS Phù Chí Hịa- Trường Đại học Đà Lạt Hướng dẫn 2: TS Phạm Ngọc Sơn- Viện Nghiên cứu hạt nhân Đà Lạt Khánh Hịa – 2020 LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu tơi thực hướng dẫn PGS.TS Phù Chí Hịa TS Phạm Ngọc Sơn Các số liệu, kết nêu luận văn tiến hành Phòng Vật lý Điện tử hạt nhân, Viện Nghiên cứu hạt nhân Đà Lạt Các số liệu, kết trung thực chưa cơng bố cơng trình nghiên cứu khác Nếu không nêu trên, xin chịu hoàn toàn trách nhiệm đề tài Khánh Hịa, tháng 07 năm 2020 Tác giả Huỳnh Thanh Sơn i LỜI CẢM ƠN Với tình cảm chân thành lịng biết ơn sâu sắc, cho phép tơi gửi lời cảm ơn chân thành tới : Quý Thầy, Cô giáo Viện Hàm lâm Khoa học Công nghệ Việt Nam giảng dạy, truyền đạt cho tơi kiến thức q báu tận tình giúp đỡ tơi q trình học tập Q Thầy, Cơ giáo Viện Nghiên cứu Ứng dụng Công nghệ Nha Trang nhiệt tình tạo điều kiện thuận lợi để giúp tơi hồn thành khóa học Đặc biệt tơi xin bày tỏ lịng biết ơn sâu sắc tới Thầy PGS.TS Phù Chí Hịa, TS Phạm Ngọc Sơn tận tình hướng dẫn giúp đỡ tơi suốt trình thực luận văn Xin chân thành cảm ơn Ban giám hiệu Trường Trung học Phổ thơng Hồng Văn Thụ, động viên giúp đỡ suốt khóa học Cảm ơn bạn học viên lớp cao học ngành Vật lý kỹ thuật K6 – Viện Hàn lâm Khoa học Công nghệ Việt Nam tơi trao đổi kiến thức suốt khóa học Và cuối cùng, tơi xin bày tỏ lịng biết ơn tới gia đình tơi, người đã, luôn động viên – chia suốt đời Khánh Hòa, tháng 07 năm 2020 Huỳnh Thanh Sơn ii DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT STT Ký hiệu Giải thích ký hiệu PHITS Hệ thống mô hạt ion nặng phương pháp Monte – Carlo, (Particle and Heavy Ion Transport code System) MCNP Chương trình tính tốn mơ phương pháp Monte – Carlo (Monte Carlo N-Particle Transport Code) MCNP5 Chương trình tính tốn mơ phương pháp Monte – Carlo phiên (Monte Carlo N-Particle version 5) BNCT PGNAA IAEA EXFOR CFNB Chương trình tính tốn cho chùm nơtron lọc (Calculation for Filtered Neutron Beams) JAEA Cơ quan Năng lượng nguyên tử Nhật Bản (Japan Atomic Energy Agency) 10 GEANT4 Chương trình mơ hình học trình chuyển động (GEometry And Tracking) 11 JENDL-4.0 Phương pháp xạ trị kỹ thuật BNCT (Boron neutron capture theraphy) Phân tích kích hoạt nơtron gamma tức thời (Prompt Gamma-ray Neutron Activation Analysis) Cơ quan Năng lượng nguyên tử Quốc tế (International Atomic Energy Agency) Thư viện số liệu phản ứng hạt nhân thực nghiệm (Experimental Nuclear Reaction Data) Thư viện số liệu hạt nhân tính tốn (Japanes Evaluated Nuclear Data Library) iii DANH MỤC CÁC BẢNG Số bảng Nội dung bảng Trang 1.1 Bảng so sánh chương trình PHITS với chương trình máy tính khác 17 2.1 Các phần cấu trúc chương trình PHITS 22 2.2 Mơ tả chế độ tính tốn 24 2.3 Mô tả lịch sử tương tác hạt 24 2.4 Tên tập tin đầu vào-đầu 25 2.5 Mô tả thơng số cho nguồn xạ hình trụ 26 2.6 Bề mặt bao quanh vật thể 30 2.7 Định dạng khối vật thể 35 2.8 Mô tả chức Tally 39 3.1 Số liệu phân bố phổ lượng nơtron chùm tia nơtron ( với chuẩn trực hình cone, phin lọc Bismuth 6cm Sapphire 15cm) 49 3.2 Số liệu phân bố phổ lượng nơtron chùm tia nơtron ( với chuẩn trực hình trụ, phin lọc Bismuth 6cm Sapphire 15cm) 53 3.3 Các tham số phổ nơtron: Thông lượng nơtron nhiệt, nhiệt, hình cone tinh thể hình cone bình thường (Sử dụng phin lọc Bismuth-6cm, Sapphire-15cm) 58 3.4 Các tham số phổ nơtron: Thông lượng nơtron nhiệt, nhiệt, hình cone tinh thể hình trụ tinh thể (Sử dụng phin lọc Bismuth-6cm, Sapphire-15cm) 59 iv 3.5 Các tham số phổ nơtron: Thông lượng nơtron nhiệt, nhiệt, hình cone tinh thể với phin lọc sử dụng Bismuth-6cm chiều dài tinh thể Sapphire thay đổi từ đến 20cm 61 3.6 Các tham số phổ nơtron: Thông lượng nơtron nhiệt, nhiệt, hình cone tinh thể (Sử dụng phin lọc Bismuth-6cm, Sapphire-15cm khơng có phin lọc) 63 v DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Số hình Nội dung hình Trang 1.1 Sơ đồ thuật tốn Monte Carlo tính tốn q trình tương tác, hạt nơtron vật liệu 12 1.2 Thống kê số liệu phổ lượng hạt 16 1.3 Thống kê mật độ thông lượng hạt 16 2.1 Khai báo véctơ phương góc khối nguồn xạ 26 2.2 Nguồn hình trụ 27 2.3 Các trường hợp nguồn hình trụ 27 2.4 Kết khai báo nguồn xạ hình trụ 28 2.5 Kết khai báo hình cầu, hình hộp, mặt phẳng hình trụ 37 2.6 Kết khai báo chuẩn trực hình trụ 38 2.7 Thống kê mật độ thông lượng hạt sử dụng lưới xyz 43 2.8 Lỗi khơng có hạt thống kê lưới 43 2.9 Hiển thị thông lương hạt với phân bố lượng 45 3.1 Mặt cắt lò phản ứng hạt nhân Đà Lạt vị trí kênh ngang số 46 3.2 Mơ hình tính tốn hệ dẫn chùm tia nơtron phin lọc hình cone kênh ngang số chương trình PHITS 48 3.3 Mơ hình tính tốn hệ dẫn chùm tia nơtron phin lọc hình trụ kênh ngang số chương trình PHITS 48 3.4 Phổ nơtron đơn hình cone tinh thể hình cone bình thường theo thang logarit 58 3.5 Phổ nơtron đơn hình cone tinh thể hình cone bình thường theo thang tuyến tính 59 vi 3.6 Phổ nơtron đơn hình cone tinh thể hình trụ tinh thể theo thang tuyến tính 60 3.7 Phổ nơtron đơn hình cone tinh thể hình trụ tinh thể theo thang logarit 60 3.8 Phổ nơtron đơn hình cone tinh thể phin lọc Sapphire thay đổi chiều dài theo thang tuyến tính 62 3.9 Phổ nơtron đơn hình cone tinh thể (với phin lọc Bismuth-6cm, Sapphire-15cm khơng có phin lọc) theo thang logarit 63 3.10 Phổ nơtron đơn hình cone tinh thể (với phin lọc Bismuth-6cm, Sapphire-15cm khơng có phin lọc) theo thang tuyến tính 64 vii MỤC LỤC Trang Lời cam đoan i Lời cảm ơn ii Danh mục ký hiệu, chữ viết tắt iii Danh mục bảng iv Danh mục hình vẽ, đồ thị vi Mục lục viii Mở đầu ……………………………………………………………… CHƯƠNG TỔNG QUAN TÀI LIỆU……… …….……… 1.1 Tình hình nghiên cứu ngồi nước …………………… 1.2 Tình hình nghiên cứu nước ………………… … 1.3 Cơ sở lý thuyết kỹ thuật phin lọc nơtron …………… 1.4 Phương pháp tính tốn dịng nơtron phin lọc ………… 1.5 Kỹ thuật tính tốn mơ Monte Carlo ……… 1.5.1 Tìm hiểu chương trình mô giới nước …………………………….……… 1.5.2 Kỹ thuật tính tốn Monte Carlo ………………… 11 1.6 Giới thiệu chương trình tính tốn PHITS………….…… 15 CHƯƠNG 2PHƯƠNG PHÁP NGHIÊN CỨU VÀ THỰC HIỆN MÔ PHỎNG……… …… 20 2.1 Hướng dẫn cài đặt thực thi PHITS …………… …… 20 2.1.1 Cài đặt hệ điều hành Windows ……………… 20 2.1.2 Thực thi PHITS hệ điều hành Windows…… 21 2.1.3 Phần mềm đề xuất ………………………… 21 viii unit = material = all axis = eng # unit is [1/cm^2/source] # (D=all) number of specific material # axis of output file = tally_MeV.dat # file name of output for the above axis part = neutron epsout = y-txt = (n/cm2/MeV/s) factor = 1.6E+11 [End] (Dấu chấm dấu phẩy Tiếng Việt Do mặc định chương trình PHITS) 90 PHỤ LỤC 4: File input mơ chuẩn trực hình trụ [Material] mat[1] H2 O1 mat[2] Al mat[3] Bi mat[4] Al O mat[5] H C -0.023 # CONCRETE -0.0023 O -1.22 Na -0.0368 Mg -0.005 Al -0.078 Si -0.775 K -0.0299 Ca -0.1 Fe -0.032 mat[6] N 4.0586E-05 # Air O 1.0800E-05 Ar 2.4255E-07 mat[7] H -0.117 # SWX-201 B -0.05 C -0.833 mat[8] Pb mat[9] Bi mt9 bism.35t mat[10] Al mt10 $ Bismuth crystal at 20oc O3 alsap.20t osapp.20t $ Sapphire crystal at 20oC mat[11] Fe mat[12] H -0.049 # SWX-277 B -0.0163 91 C -0.9347 mat[13] H 2.0 C $ Polythylene (JIS-K6922-1 1997-3-1) 0.942 g/cm3 1.0 [Surface] 10 CZ 20 11 pz 12 pz 250 13 CZ 10.05 14 pz 50 15 pz 130 16 CZ 9.05 17 CZ 7.5 18 CZ 7.0 19 CZ 6.0 20 PZ 10 21 PZ 16 23 PZ 115 24 PZ 130 26 PZ 140 27 PZ 150 28 PZ 160 29 PZ 170 30 PZ 180 31 PZ 190 32 PZ 200 33 PZ 210 34 PZ 220 35 CZ 2.5 36 PZ 40 92 37 CZ 5.5 38 PZ 370 39 CZ 1.5 40 CZ 41 CZ 20 42 CZ 30 43 CZ 40 44 PZ 380 45 PZ 390 46 PZ 420 47 PZ 435 48 PZ 450 49 CZ 60 50 cz [Cell] 100 -1 49:-11:48 200 -49 43 11 -48 101 -2.302 -43 17 11 -36 201 -2.302 -43 13 14 -34 102 -2.7-13 16 14 -15 103 -0.95 -16 17 14 -15 104 -2.7-17 18 11 -15 105 -2.7-18 19 11 -15 106 -1.21e-3 -19 11 -20 107 -9.8-19 20 -21 109 -1.21e-3 -19 21 -23 110 10 -3.98 -35 23 -24 111 -1.21e-3 -39 24 -34 112 -0.95 -13 39 24 -26 93 113 -11.34 -13 39 26 -27 114 -0.95 -13 39 27 -28 115 -0.95 -13 39 28 -29 116 -11.34 -13 39 29 -30 117 -0.95 -13 39 30 -31 118 -0.95 -13 39 31 -32 119 -11.34 -13 39 32 -33 120 -0.95 -13 39 33 -34 121 -0.95 -19 35 23 -15 122 11 -7.874 -41 17 36 -14 123 -2.302 -43 41 36 -14 400 0-43 41 34 -38 401 -11.34 -41 50 34 -38 402 -0.95 -50 40 34 -38 403 -1.21e-3 -40 404 -11.34 -43 40 38 -44 405 12 -1.68 -42 40 44 -45 406 -11.34 -43 42 44 -47 407 12 -1.68 -42 41 45 -47 34 -38 408 12 -1.68 -41 46 -47 409 -11.34 -43 47 -48 410 -1.21e-3 -40 38 -45 411 -1.21e-3 -41 45 -46 (Dấu chấm dấu phẩy Tiếng Việt Do mặc định chương trình PHITS) 94 Phụ lục : Các nội dung hướng dẫn thiết kế file – input chương trình mơ PHITS Nội dung 1: Các phần PHITS STT 10 11 12 13 14 15 16 17 18 19 20 21 22 TÊN [title] [parameters] [source] [material] [surface] [cell] [transform] [temperature] [mat time change] [magnetic field] [electro magnetic field] [delta ray] MIÊU TẢ Tiêu đề Nhiều loại tham số Định nghĩa nguồn Định nghĩa vật liệu Định nghĩa bề mặt Định nghĩa khối vật thể Định nghĩa biến đổi tọa độ Định nghĩa nhiệt độ Định nghĩa vật liệu phụ thuộc thời gian Định nghĩa từ trường Định nghĩa điện từ trường Định nghĩa sản xuất tia  Định nghĩa cho việc sử dụng chế độ mô cấu [track structure] trúc theo dõi [super mirror] Định nghĩa siêu gương cho neutron lượng thấp Định nghĩa tùy chọn tán xạ đàn hồi cho neutron [elastic option] lượng thấp Định nghĩa lượng tối đa (dmax) hạt [data max] nhân để sử dụng thư viện liệu [frag data] Định nghĩa mặt cắt người dùng xác định [importance] Định nghĩa quan trọng khu vực [weight Định nghĩa cửa sổ trọng lượng khu vực window] Định nghĩa giá trị để ưu tiên [cửa sổ trọng [ww bias] lượng] [forced Định nghĩa va chạm cưỡng collisions] [volume] Định nghĩa thể tích khu vực 95 23 24 25 26 27 28 29 30 31 32 33 34 [multiplier] [mat name color] [reg name] [counter] [timer] [t-track] [t-cross] [t-point] [t-deposit] [t-deposit2] [t-heat] [t-yield] 35 36 [t-product] 37 38 [t-let] 39 40 41 [t-dpa] [t-sed] [t-time] [t-interact] (former name [t-star]) [t-dchain] 42 43 44 [t-wwg] [t-wwbg] 45 46 47 48 [t-userdefined] [t-gshow] [t-rshow] [t-3dshow] [t-volume] Định nghĩa số nhân Tên vật liệu định nghĩa màu cho đồ họa Định nghĩa tên vùng cho lô đồ họa Định nghĩa truy cập Định nghĩa thời gian đồng hồ Thông lượng hạt khu vực định Thông lượng hạt bề mặt định Thông lượng hạt điểm định Năng lượng hấp thụ thể tích định Năng lượng hấp thụ hai thể tích định Sản sinh nhiệt khu vực định Năng suất hạt nhân lại khu vực định Sản phẩm phản ứng khu vực định Dịch chuyển nguyên tử (DPA) khu vực định Phân phối LET khu vực định Phân phối số lượng vi mô khu vực định Thông tin thời gian hạt khu vực định Số lượng tương tác xảy khu vực định Năng suất hạt nhân lại (cho việc sử dụng DCHAIN) Thông số đầu cho [cửa sổ trọng lượng] Các tham số đầu cho [ww bias] Tự động tính toán khối lượng cho khu vực định Bất kỳ số lượng xác định người dùng Trực quan hình học 2D Trực quan hình học 2D với đại lượng vật lý Trực quan hình học 3D 96 49 [end] Kết thúc tập tin đầu vào Nội dung 2: Mơ tả chế độ tính tốn chương trình PHITS Tham số Giá trị icntl Chú thích (D=0) Tùy chọn điều khiển = 0Tính tốn PHITS bình thường = 1Tính tốn phản ứng hạt nhân, (đang phát triển) =3Chỉ xuất tiếng vang đầu vào để kiểm tra việc sử dụng nhớ, thư viện liên kết tệp =5Tính tốn giả sử tất vùng void (khơng có phản ứng, khơng ion hóa) để kiểm tra hình học, tính tốn khối lượng diện tích = 6Kiểm tra nguồn Vị trí nguồn tính [t-product] = 7Thực thi [t-gshow] cho đầu đồ họa =8Đầu hình học kiểm kê lưới xyz với tùy chọn gshow cho đầu đồ họa = 9Thực [t-rshow] cho đầu đồ họa = 10 Đầu hình học reg lưới tally với tùy chọn rshow cho đầu đồ họa = 11 Thực [t-3dshow] cho đầu đồ họa = 12 Tính tốn lại cách sử dụng tập tin dumpall Tệp Dumpall định tệp (15) = 13 Sử dụng hàm để tổng hợp kết kiểm đếm (Hàm Sumtally) = 14 Sử dụng tính tốn tự động thể tích khu vực = 15 Thực [t-wwbg] để thu tham số [ww bias] 97 Nội dung 3: Mô tả số lịch sử chương trình PHITS Tham số Giá trị Chú thích irskip Kiểm soát số ngẫu nhiên Bỏ qua số lượng lịch sử irskip hạt số ngẫu nhiên để gỡ lỗi (Tổng số lịch sử giảm irskip.) Bắt đầu tính tốn với số irskip + số hạt ngẫu nhiên để thực tính tốn song song tay (Tổng số lịch sử không thay đổi.) (D=0) Irskip > Irskip < rseed (D=0.0) rseed < rseed = rseed > Tùy chọn số ngẫu nhiên ban đầu Nhận số ngẫu nhiên ban đầu từ thời gian bắt đầu 6647299061401 (mặc định) Sử dụng rseed hạt giống ban đầu số ngẫu nhiên maxcas (D = 10) maxbch (D = 10) Số lượng lịch sử đợt Giới hạn 2147483647 Số lượng đợt Giới hạn 2147483647 maxbnk (D=10000) Số lượng hạt tạm thời lưu nhớ timeout (D = -1.0) Cắt thời gian CPU để tính tốn PHITS (tính theo đơn vị giây) istdev (D=0) Tham số điều khiển cho loại tính tốn độ khơng đảm bảo thống kê chế độ khởi động lại Chế độ tính tốn khởi động lại kích hoạt; khơng có kết kiểm đếm q khứ, phép = -2 = -1 =0 =1 tính bắt đầu với istdev = Chế độ tính tốn khởi động lại kích hoạt; khơng có kết kiểm đếm khứ, phép tính bắt đầu với istdev = istdev tự động đặt thành cho tính tốn song song nhớ chia sẻ cho trường hợp khác Một tính tốn bắt đầu Độ khơng đảm bảo 98 =2 ireschk (D=0) =0 =1 thống kê ước tính từ độ lệch chuẩn kết lơ Một tính tốn bắt đầu Độ khơng đảm bảo thống kê ước tính từ độ lệch chuẩn kết lịch sử Tùy chọn chọn chế độ song song chia sẻ Kiểm sốt tham số để kiểm tra tính quán kiểm đếm Kiểm tra tính quán cài đặt kiểm đếm cũ Không kiểm tra Tùy chọn hữu ích thiết lập kiểm đếm phức tạp có cài đặt khơng ghi đầy đủ tệp đầu kiểm đếm Nội dung 4: Tên tập tin đầu vào-đầu chương trình PHITS Tham số Giá trị Chú thích File(1) (D=c:/phits) Tên thư mục cài đặt PHITS Nếu tham số đặt đúng, không cần định tên tệp đầu vào khác, tức tệp (7, 20, 21, 24, 25 26), trừ cấu trúc thư mục PHITS không thay đổi Nếu tham số định trực tiếp tệp đầu vào, tham số định ghi đè cài đặt File(6) (D=phits.out) Tóm tắt tên tập tin đầu File(7) (D=file(1)/data/xsdir.jnd) Tên tập tin mặt cắt ngang File(11) (D=nuclcal.out) Tên tập tin đầu phản ứng hạt nhân File(12) (D=fort.12) Cắt tên tập tin đầu neutron File(13) (D=fort.13) Cắt tên tập tin đầu tia  99 File(10) (D=fort.10) Cắt tên tập tin đầu proton File(15) (D=dumpall.dat) tên tệp kết xuất cho tùy chọn dumpall = File(18) (D=voxel.bin) Tên tệp sử dụng ivoxel = 1, File(20) (D=file(1)/XS/egs/) Thư mục chứa liệu thư viện cho EGS5 File(21) (D=file(1)/dchain-sp/data/) Thư mục chứa liệu thư viện cho DCHAIN File(22) (D=batch.out) Tên tệp để thông tin đầu lô File(23) (D=pegs5) Tên tệp đầu cho PEGS5 File(24) (D=file(1)/data) Thư mục chứa RIsource.dat (cơ sở liệu DECDC) File(25) (D=file(1)/XS/tra) Thư mục chứa sở liệu mặt cắt sử dụng để mô cấu trúc theo dõi File(26) (D=file(1)/data/multiplier) Thư mục chứa phần [multiplier] xác định trước File(27) (D=file(1)/XS/yield/) Thư mục chứa liệu mặt cắt kích hoạt người dùng xác định Cần có thứ cuối đường dẫn File(28) (D=file(1)/data/aama.dat) Tên tệp đầu vào cho chương trình MUONIC ATOM CASCADE để sử dụng tham số người dùng xác định Nội dung 5: Các loại nguồn chương trình PHITS Loại nguồn Giải thích s-type = 1, (4) Hình trụ (hoặc hình trịn, bút chì) s-type = 2, (5) Hình chữ nhật đặc (hoặc hình chữ nhật) s-type = 3, (6) Gaussian (x, y, z độc lập) s-type = 7, (8) Parabola chung (x, y, z độc lập) 100 s-type = 9, (10) s-type = 11 Hình cầu vỏ hình cầu Phân bố đồng không gian pha thẳng đứng với hướng chùm s-type = 12 Đọc liệu từ đầu decay-turtle s-type =13, (14) Gaussian (mặt phẳng xyx) s-type = 15, (16) Parabola chung (mặt phẳng xy) s-type = 17 Đọc tập tin dump s-type = 18, (19) Hình nón s-type = 20, (21) Lăng kính tam giác s-type = 22, (23) phân phối lưới xyz s-type = 24 Nguồn Tetra-lưới s-type = 26 Nguồn định nghĩa bề mặt s-type = 100 Nguồn định nghĩa người dùng chỉnh sửa usrsors.f biên dịch PHITS Nội dung 6: Các thơng số nguồn chung chương trình PHITS Tham số Giải thích proj = Projectile: xem Bảng 2.8 để biết thông số kỹ thuật sx = sy = sz = (D=0) thành phần spin theo phương x (D=0) thành phần spin theo phương y (D=0) thành phần spin theo phương z reg = (D = all) Vùng nguồn giới hạn (các) chồng chéo vùng xác định s-type vùng định tham số Định dạng reg = {1 - 5} 10 34 Khung mạng vũ trụ sử dụng reg = (6