›
›

Bài giảng Khuếch tán và làm chậm Neutron

ppt 58 trang phuongnguyen 3941
Download
Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Khuếch tán và làm chậm Neutron", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pptbai_giang_khuech_tan_va_lam_cham_neutron.ppt

Nội dung text: Bài giảng Khuếch tán và làm chậm Neutron

  1. KHUẾCH TÁN VÀ LÀM CHẬM NEUTRON Tham khảo chính: “Introduction to Nuclear Engineering” của John R. Lamarsh
  2. Nội dung ➢ Mở đầu ➢ Thông lượng Neutron ➢ Định luật Fick ➢ Phương trình liên tục ➢ Phương trình khuếch tán ➢ Điều kiện biên ➢ Các nghiệm của phương trình khuếch tán ➢ Độ dài khuếch tán ➢ Phương pháp khuếch tán nhóm ➢ Khuếch tán neutron nhiệt ➢ Tính toán làm chậm neutron hai nhóm
  3. Mở đầu • Để thiết kế một lò phản ứng hạt nhân một cách đúng đắn cần phải tiên đoán neutron phân bố như thế nào trong hệ thống. • Đây là một bài toán khó do lộ trình của các neutron là rất phức tạp do kết quả của các va chạm hạt nhân liên tiếp. • Một trong các phép gần đúng là coi neutron “khuếch tán” (diffusion) như là một chất khí trong một chất khí khác. • Phân bố neutron có thể tìm được bằng cách giải phương trình khuếch tán. • Cách làm này gọi là gần đúng khuếch tán và hiện nay vẫn còn được dùng rộng rãi.
  4. Thông lượng neutron (1) Số va chạm trên mỗi cm3/s khi một chùm neutron đi vào một bia mỏng F = t I Trường hợp có vài chùm neutron F = t (I A + I B + IC + ) F = t (nA + nB + nC + )v F = t nv F = t  = nv là thông lượng neutron (flux) [cm-2.s-1]
  5. Thông lượng neutron (2) Dễ dàng mở rộng kết quả này cho các neutron có một phân bố năng lượng F = t (E)n(E)v(E)dE = t (E)(E)dE 0 0 trong đó, (E) = n(E)v(E) là thông lượng phụ thuộc năng lượng (energy- dependent flux or flux per unit energy) Phương trình trên là tốc độ tương tác tổng. Tốc độ tương tác từng phần (tán xạ, hấp thụ) có thể xác định bằng các biểu thức tương tự Fs = s (E)(E)dE Fa = a (E)(E)dE 0 0
  6. Định luật Fick (1) Lý thuyết khuếch tán dựa trên định luật Fick mà ban đầu được dùng cho khuếch tán hóa học: Một chất tan khuếch tán từ vùng có nồng độ cao hơn tới vùng có nồng độ thấp hơn; Lưu lượng chất tan tỷ lệ với gradient âm của nồng độ chất tan. Cách di chuyển của các neutron cũng giống như cách của chất tan trong dung dịch.
  7. Định luật Fick (2) Giả sử thông lượng neutron biến đổi dọc theo trục x như trong hình vẽ. Định luật Fick được viết như sau d J = −D (x) x dx J là số neutron thực đi qua x J trên mỗi đơn vị thời gian x qua mỗi đơn vị diện tích vuông góc với hướng x (neutrons/cm2.s); D là hệ số khuếch tán (diffusion x coefficient) [cm]
  8. Định luật Fick (3) Mật độ neutron tăng J (r) = −D(r) Dòng Gradient J (r) (r)
  9. Định luật Fick (4) Trong không gian 3 chiều và các neutron đơn năng J = −Dgrad = −D Định nghĩa toán tử gradient (tọa độ vuông góc) Số hạt đi qua bề mặt theo hướng n trong mỗi đơn vị thời gian trên mỗi đơn vị diện tích (Dòng pháp tuyến-normal current) J.n = J n
  10. Định luật Fick (5) Là quãng đường tự do trung bình vận chuyển -transport mean free paths Cosine góc tán xạ trung bình trong đó A là nguyên tử khối.
  11. Định luật Fick (6) Lưu ý: Định luật Fick không phải là mối quan hệ chính xác. Nó không có hiệu lực trong một số trường hợp: 1) Trong một môi trường hấp thụ neutron mạnh; 2) Trong khoảng 3 lần quãng đường tự do trung bình từ nguồn neutron hoặc từ bề mặt của một môi trường; 3) Khi tán xạ neutron là không đẳng hướng (mạnh).
  12. Phương trình liên tục (1) Để ý một thể tích V bất kì trong một môi trường chứa neutron. Theo thời gian số neutron trong V có thể thay đổi. Phương trình liên tục (equation of continuity) biểu diễn sự bảo toàn neutron: Rate of Rate of Rate of Rate of change in neutron neutron neutron the number = production - absorption - leakage of neutrons in volume in volume from V in a small V V volume V
  13. Phương trình liên tục (2) ➢ Gọi n là mật độ neutron tại tại một điểm và thời điểm nào đó trong V. Tốc độ thay đổi (the rate of change) số neutron là d ndV dt V n hay The rate of change = dV V t
  14. Phương trình liên tục (3) ➢ Gọi s là tốc độ phát ra neutron từ nguồn trong mỗi cm3 trong V. Tốc độ neutron sinh ra trong toàn bộ thể tích V được cho bởi Production rate = sdV V 3 ➢ Tốc độ mất neutron do hấp thụ trong mỗi cm /s:  a Trong toàn bộ thể tích V tổng số neutron mất đi trong mỗi giây do hấp thụ là Absorption rate = adV V
  15. Phương trình liên tục (4) ➢ Để ý dòng neutron đi vào và đi ra khỏi V. Nếu J là vector mật độ dòng neutron trên bề mặt của V và n là vector pháp tuyến hướng ra ngoài từ bề mặt. Số neutron thực (net) đi qua bề mặt ra ngoài trên mỗi cm2/s là J.n Tốc độ rò neutron toàn phần (có thể âm hoặc dương) qua bề mặt A của thể tích là Leakage rate = J.ndA A
  16. Phương trình liên tục (5) Tích phân theo bề mặt này có thể chuyển thành một tích phân theo thể tích bằng cách dùng lý thuyết phân kì (divergence) J.ndA = divJdV A V Leakage rate = divJdV V
  17. Phương trình liên tục (6) Sử dụng những kết quả trên đây chúng ta thu được phương trình liên tục như sau: n dV = sdV −  adV − divJdV V t V V V n Dạng tổng quát = s −   − divJ t a Nếu mật độ neutron không phụ thuộc thời gian thì phương trình liên tục ở trạng thái dừng (steady- state equation of continuity) sẽ là: divJ + a − s = 0
  18. Phương trình khuếch tán (1) Nhớ lại định luật Fick: J = −D Thay định luật Fick và phương trình liên tục ta có phương trình khuếch tán neutron (diffusion equation) phụ thuộc thời gian: (với neutron đơn năng 2 n D  −  a + s = và hệ số khuếch tán là t hằng số) 2 = div grad is called Laplacian
  19. Phương trình khuếch tán - laplacian (2)
  20. Phương trình khuếch tán (3) Vì  = nv 1  D 2 −   + s = a v t Trường hợp không phụ thuộc thời gian, phương trình khuếch tán dừng (steady-state diffusion): 2 D  − a + s = 0 1 s  2 −  = − L2 D D L2 = là diện tích khuếch tán [cm2]  a L là độ dài khuếch tán [cm]
  21. Các điều kiện biên (1) Thông lượng neutron có thể thu được bằng cách giải phương trình khuếch tán. Trước khi giải phương trình khuếch tán cần phải chỉ rõ các điều kiện biên (boundary conditions) mà nghiệm của phương trình phải thỏa mãn.
  22. Điều kiện biên (2) ➢ Thông lượng neutron phải thực, không âm và hữu hạn ➢ Neutron không trở lại từ bề mặt tự do (chân không) Diffusion Vacuum or air Độ dài ngoại suy (?) theory d = 0,71tr Exact Surface flux tr = 3D d = 2,13D Diffusing medium Thực nghiệm cho thấy các vật liệu không phải d là khí thì D<1cm nên thực tế có thể giả thiết thông lượng triệt tiêu ở bề mặt thật của hệ
  23. Điều kiện biên (3) ➢ Thông lượng và dòng neutron là liên tục tại bề mặt giữa hai môi trường khác nhau (ví dụ như vùng hoạt và vành phản xạ) A = B (J A )n = (J B )n Các điều kiện biên này đôi khi được gọi là interface boundary conditions ➢ Tính đối xứng
  24. Nghiệm của phương trình khuếch tán (1) ❑ Nguồn tấm phẳng vô hạn Xét một nguồn tấm phẳng vô hạn phát ra S neutrons trong mỗi cm2/s trong môi trường khuếch tán vô hạn d 2 1 s −  = − , x = 0 dx2 L2 D 2x d 2 1 −  = 0, x 0 J dx 2 L2 J Điều kiện biên: thông lượng là hữu hạn khi x tiến đến vô cùng x=0 plane
  25. Nghiệm của phương trình khuếch tán (2) ❑ Nguồn tấm phẳng vô hạn (tt) Nghiệm tổng quát:  = Ae−x L + Cex L A và C là hằng số được xác định nhờ các điều kiện biên: C=0 để thông lượng hữu hạn khi x tiến đến vô cùng, nên  = Ae−x L Trong giới hạn khi x tiến đến zero, dòng thực đi ra từ box (h/vẽ) phải tiến dần đến mật độ nguồn của tấm S lim J (x) = (Gọi là điều kiện nguồn) x→0 2 Tính đối xứng
  26. Nghiệm của phương trình khuếch tán (3) ❑ Nguồn tấm phẳng vô hạn (tt) Từ định luật Fick, d DA J = −D = e −x L Dùng điều kiện dx L nguồn SL A = 2D SL  = e− x L 2D Do tính đối xứng của bài toán, thông lượng phải giống nhau tại –x cũng như +x SL  = e −|x| L 2D
  27. Nghiệm của phương trình khuếch tán (4) ❑ Nguồn điểm Xét nguồn điểm phát S n/s đẳng hướng trong môi trường vô hạn. Với Laplacian trong trục tọa độ cầu, phương trình khuếch tán trở thành, 1 d d 1 r 2 −  = 0 (tại r khác zero) r 2 dr dr L2 Số neutron thực đi qua bề mặt một hình cầu nhỏ bán kính r bao quanh nguồn trong một đơn vị thời gian là 4 r 2 J(r) S lim r 2 J (r) = (điều kiện nguồn) r→0 4
  28. Nghiệm của phương trình khuếch tán (5) ❑ Nguồn điểm (tt) Sử dụng một biến mới  = r d 2 1 −  = 0 dr 2 L Nghiệm tổng quát  = Ae−r L + Cer L e−r L er L  = A + C r r
  29. Nghiệm của phương trình khuếch tán (6) ❑ Nguồn điểm (tt) Vì thông lượng phải là hữu hạn khi r tiến đến vô cùng, C phải bằng không. Hằng số A thu được từ điều kiện nguồn d 1 1 J = −D = DA + e −r L dr rL r 2 S A = S 4 D lim r 2 J (r) = r→0 4 Se −r L  = 4 Dr
  30. Nghiệm của phương trình khuếch tán (7) ❑ Tấm trần Khảo sát một tấm vô hạn có độ dày 2a có một nguồn tấm phẳng vô hạn ở trung tâm phát ra S n/cm2.s. Điều kiện biên: (a + d) = (−a − d) = 0 Đối với x dương  = Ae−x L + Cex L (a + d) = Ae−(a+d ) L + Ce(a+d ) L = 0 −2(a+d ) L Điều kiện C = Ae biên
  31. Solutions of the diffusion equation (8) ❑ Tấm trần (tt) Nghiệm phương trình trở thành  = Ae−x L − e x L−2(a+d ) L  Sử dụng điều kiện nguồn (tấm phẳng vô hạn) SL −1 A = (1+ e −2(a+d ) L ) 2D SL e−x L − e x L−2(a+d ) L  = 2D 1+ e−2(a+d ) L SL e−|x| L − e|x| L−2(a+d ) L Do đối xứng, ta có  = 2D 1+ e−2(a+d ) L
  32. Nghiệm của phương trình khuếch tán (9) ❑ Tấm trần (tt) Nhân tử và mẫu với e(a+d ) L ta được: SL e(a+d −|x|) L − e −(a+d −|x|) L  = 2D e(a+d ) L + e −(a+d ) L SL sinh (a + d − | x |) L = 2D cosh(a + d) L sinh và cosh là các hàm hyperbol được định nghĩa: e x − e−x e x + e−x sinh x = cosh x = 2 2
  33. Nghiệm của phương trình khuếch tán (10) ❑ Bài toán nhiều vùng (multi-region) Xét một tấm vô hạn như vừa rồi nhưng được bao bọc bởi một môi trường khác dày vô hạn
  34. Nghiệm của phương trình khuếch tán (11) ❑ Bài toán nhiều vùng (tt) Ta có phương trình khuếch tán cho các vùng d 2 1 1 −  = 0 for x a , x 0 2 2 1 2 dx L1 d 2 1 2 −  = 0 for x  a 2 2 2 2 dx L2 Điều kiện biên 1. 2 = finite as x → S 2. lim J (x) = (source condition) x→0 2 3. 1 ( a / 2) = 2 ( a / 2) continuity of the flux at the interface d1 d2 4. D1 = D2 continuity of current at the interface dx x= a / 2 dx x= a / 2
  35. Nghiệm của phương trình khuếch tán (12) ❑ Bài toán nhiều vùng (tt) Nghiệm tổng quát 1 (x) = A1 cosh( x / L1 ) + C1 sinh( x / L1 ) −x / L2 x / L2 2 (x) = A2e + C2e Như vậy ta có 4 ẩn số và 4 điều kiện biên. Từ B.C. 1 ta có ngay C2=0 Từ B.C. 2: x / L1 −x / L1 x / L1 −x / L1 S d1 D1 A1 e − e D1C1 e + e = J(0) = −D1 = − − 2 dx L 2 L 2 0 1 x=0 1 x=0 SL1 C1 = − 2D1
  36. Nghiệm của phương trình khuếch tán (13) ❑ Bài toán nhiều vùng (tt) Từ B.C. 3 và 4: SL a 1 a −a / 2L2 − A1 cosh + sinh = A2e 2 ph/tr 2L1 2D1 2L1 chứa 2 D A D A ẩn số. 1 1 a S a 2 2 −a / 2L2 − sinh + cosh = e Giải cho L1 2L1 2 2L1 L2 A1 và A2 SL1 D1L2 cosh(a / 2L1 ) + D2 L1 sinh( a / 2L1 ) A1 = 2D1 D2 L1 cosh(a / 2L1 ) + D1L2 sinh( a / 2L1 ) a / 2L1 SL1L2 e A2 = 2 D2 L1 cosh(a / 2L1 ) + D1L2 sinh( a / 2L1 ) Do đó ta biết được thông lượng ở 2 vùng.
  37. Độ dài khuếch tán (1) Để xem xét ý nghĩa vật lý của độ dài khuếch tán chúng ta khảo sát một nguồn điểm đơn năng phát ra S n/s trong một môi trường làm chậm đồng nhất vô hạn. Hành trình của một neutron trong một môi trường làm chậm Neutron Neutron emitted here r absorbed here Actual path of neutron
  38. Độ dài khuếch tán (2) Se −r L Biểu thức thông lượng từ một nguồn điểm,  = L đặc trưng cho sự suy giảm thông lượng 4 Dr Gọi dn là số neutron bị hấp thu trong khoảng r và r+dr dn = a(r)dV dV = 4 r 2 dr S S dn = a re −r L dr = re −r L dr D L2
  39. Độ dài khuếch tán (3) Xác suất để một neutron từ nguồn bị hấp thụ trong khoảng dr 1 p(r)dr = re −r L dr L2 Trung bình bình phương khoảng cách từ nguồn đến điểm mà tại đó một neutron bị hấp thụ có thể thu được bằng cách lấy trung bình r2 theo phân bố xác suất p(r)dr r 2 = r 2 p(r)dr 0 1 = r 3e −r L dr = 6L2 L2 0 2 1 2 L = 6 r
  40. Phương pháp khuếch tán nhóm (1) ▪ Chúng ta vừa giải phương trình khuếch tán cho một số trường hợp đặc biệt với giả thiết neutron khuếch tán là đơn năng. ▪ Các neutron phát ra từ phản ứng phân hạch với một phổ năng lượng liên tục sau đó mở rộng ra khi neutron khuếch tán trong hệ thống, mất năng lượng trong các va chạm đàn hồi và không đàn hồi. ▪ Một trong các cách hiệu quả nhất để tính sự làm chậm và khuếch tán neutron là phương pháp khuếch tán nhóm (group diffusion method).
  41. Phương pháp khuếch tán nhóm (2) Các nhóm năng lượng cho tính toán khuếch tán nhóm Group 1 Group 2 Group g Group N-1 Increasing Energy Increasing Group N
  42. Phương pháp khuếch tán nhóm (3) Thông lượng neutron nhóm g g = (E)dE g Các neutron biến khỏi nhóm gth: (1) phản ứng hấp thụ; và (2) va chạm tán xạ Absorption rate = a (E)(E)dE g a (E)(E)dE g ag = (E)dE g Absorption rate = agg
  43. Phương pháp khuếch tán nhóm (4) Gọi tiết diện chuyển nhóm từ  g→h nhóm g sang nhóm h N Total transfer rate out of g = g→hg h=g+1 g−1 Total transfer rate into g = h→gh h=1
  44. Phương pháp khuếch tán nhóm (5) Kết hợp các tốc độ hấp thụ, tốc độ chuyển nhóm chúng ta thu được phương trình khuếch tán dừng cho nhóm neutron g như sau N g−1 2 Dg g − agg − g→hg + h→gh = −sg h=g+1 h=1 Dg là hệ số khuếch tán nhóm được cho bởi 1 1 1 Dg = D(E)(E)dE = (E)dE g g g 3g tr (E) tr là tiết diện vận chuyển, và sg là tổng số neutron phát ra trong mỗi cm3/s vào trong nhóm.
  45. Khuếch tán neutron nhiệt (1) Phân bố năng lượng của neutron nhiệt được cho bởi hàm Maxwellian 2 n n(E) = E1 2e −E kT ( kT)3 2 Thông lượng phụ thuộc năng lượng của neutron nhiệt (E) = n(E)v(E) 1 2 2 n 2 = Ee−E kT ( kT)3 2 m 1 2 2E Bằng cách sử dụng công thức v(E) = m
  46. Khuếch tán neutron nhiệt (2) Kí hiệu thông lượng neutron nhiệt một nhóm là T T = (E)dE T Kí hiệu T có nghĩa là tích phân được lấy trên vùng năng lượng nhiệt thường cho đến khoảng 5kT ( 0.1 eV) nhưng do trên năng lượng này thông lượng giảm rất nhanh nên có thể lấy tích phân đến vô cùng 1 2 2 n 2 −E kT T = Ee dE ( kT)3 2 m 0 1 2 2n 2kT = m
  47. Khuếch tán neutron nhiệt (3) 1 2 Sẽ rất thuận tiện nếu thay E = kT và mvT = ET T 2 Giá trị của ET và vT có thể được tính từ các công thức sau −5 ET = 8,617T 10 eV 1 2 4 vT =1,284T 10 cm/s T[K] 2 T = nvT
  48. Khuếch tán neutron nhiệt (4) Quan hệ giữa T và thông lượng tại 2200 m/s 0 = nv0 có thể tìm được bằng cách chia 0 cho T  v 0 = 0 T 2 vT 1 2  T 0 = 0 T0 = 293.61 K T 2 T Về ý nghĩa vật lý thì T là tất cả thông lượng neutron nhiệt gộp lại bằng cách tính tích phân phân bố Maxwell trong khi 0 là giả-thông lượng được tính bằng cách giả thiết tất cả các neutron nhiệt chỉ có một năng lượng.
  49. Khuếch tán neutron nhiệt (5) Hệ số khuếch tán trong tính toán một nhóm khuếch tán neutron nhiệt được kí hiệu là D 1 1 1 Dg = D(E)(E)dE = (E)dE g g g 3g tr (E) Tiết diện hấp thụ một nhóm neutron nhiệt 1 a =  a (E)(E)dE T T Tích phân là tốc độ hấp thụ toàn phần và bằng ga (T)a (E0 )0 với ga(T) h/s non-1/v hay h/s Westcott a = g a (T)a (E0 )0 T 1 2 T0 = g a (T)a (E0 ) 2 T
  50. Khuếch tán neutron nhiệt (6) Các thông số khuếch tán neutron nhiệt của các chất làm chậm thông thường ở 20oC Density, 3 Moderator g/cm −1 2 2 D , cm a , cm LT , cm LT , H2O 1,00 0,16 0,0197 8,1 2,85 −5 3 D2O 1,10 0,87 9,3 10 9,4 10 97 −3 Be 1,85 0,50 1,04 10 480 21 Graphit 1,60 0,84 2,4 10−4 3500 59 [Argonne National Laboratory ANL-5800, 1963]
  51. Khuếch tán neutron nhiệt (7) Với thông lượng nhiệt đã được định nghĩa và các tham số (hệ số khuếch tán, tiết diện hấp thụ), phương trình khuếch tán một nhóm đối với neutron nhiệt có thể được viết như sau 2 D T − aT = −sT 2 1 sT  T − 2 T = − LT D 2 D LT = là diện tích khuếch tán nhiệt a (thermal diffusion area); LT là độ dài khuếch tán nhiệt (thermal diffusion length)
  52. Tính toán làm chậm neutron hai nhóm (1) Trong một số tính toán lò, ít nhất hai nhóm neutron phải được sử dụng để thu được kết quả đủ chính xác. Một nhóm mô tả neutron nhiệt còn nhóm nhanh gồm tất cả các neutron có năng lượng trên năng lượng nhiệt. Khảo sát một bài toán tính thông lượng hai nhóm như là một hàm của khoảng cách từ một nguồn điểm phát ra S neutron phân hạch (neutron nhanh) trong mỗi giây trong một môi trường làm chậm đồng nhất vô hạn.
  53. Tính toán làm chậm neutron hai nhóm (2) Nhớ lại phương trình khuếch tán nhóm tổng quát N g−1 2 Dg g − agg − g→hg + h→gh = −sg h=g+1 h=1 Giả sử  = 0 (good moderator) a x khác 0, 2 S=0 D1 1 − 11 = 0 2 1 11  T − 2 T = − LT D 1 = 1→2
  54. Tính toán làm chậm neutron hai nhóm (3) Để tìm T từ phương trình khuếch tán nhiệt, đầu tiên cần giải phương trình khuếch tán nhanh để có 1 Dùng laplacian của tọa độ cầu 1 d 2 d1 1 2 r − 1 = 0 r dr dr  T D 1 2  T = là tuổi neutron (age) [cm ] 1 Tương tự như Se −r T trong trường hợp giải 1 = 4 rD1r khuếch tán 2 neutron từ một SL −r  T −r LT T nguồn điểm T = 2 (e − e ) 4 rD(LT −T )
  55. Ví dụ 5.1 (1) Bài toán: Một tấm vô hạn dày 2a có một nguồn tấm phẳng vô hạn tại trung tâm phát ra S neutrons trên mỗi cm2/s. (a) Tìm một biểu thức cho số neutron thoát ra trong một giây từ 1cm2 trên hai mặt của tấm; (b) Hỏi xác suất để một neutron từ nguồn thoát ra khỏi tấm? Lời giải: a) Để ý trên bề mặt tại x=a. Sử dụng định luật Fick ta có d d SL sinh (a + d − x) L J x = −D J(a) = − D dx dx 2D cosh(a + d) L x=a SL sinh (a + d− | x |) L S sinh( d L)  = = 2D cosh(a + d) L 2 cosh(a + d) L Do đối xứng neutron thoát ra hai mặt là cosh(d L) 2J(a) = S [ans] cosh(a + d) L
  56. Ví dụ 5.1 (2) b) Xác suất để một neutron nguồn thoát khỏi tấm bằng số neutron thoát trên mỗi cm2/s chia cho số neutron phát ra trên mỗi cm2/s từ nguồn 2J (a) cosh(d L) = [ans] S cosh(a + d) L
  57. Ví dụ 5.2 (1) Bài toán: Một nguồn điểm phát ra 107 neutron nhiệt trong mỗi giây được đặt trong một vật thể vô hạn chứa nước mật độ đơn vị tại nhiệt độ phòng. Hỏi thông lượng neutron nhiệt cách nguồn 15cm? LT = 2,85 cm D = 0,16 cm Lời giải: Thông lượng được cho bởi Se −r LT  = T 4 Dr Sử dụng các giá trị của S, r, LT và D đã cho 7 −15 2,85 10 e 3 2  = = 1,72 10 neutron/cm s [ans.] T 4 0,16 15
  58. Chicago Pile 1