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

KHUẾCH TÁN VÀ LÀM CHẬM NEUTRON

Tham khảo chính: “Introduction to Nuclear Engineering” của John R. Lamarsh

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

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.

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

(cid:0) (cid:0)

Trường hợp có vài chùm neutron

(

...)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

(cid:0)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) F n n n v ( ...)

F (cid:0)

(cid:0) (cid:0) (cid:0)

(cid:0) là thông lượng neutron (flux) [cm-2.s-1]

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

(cid:0) (cid:0)

(

dEEvEnE ()

()

)

(

(cid:0) ()

dEE )

(cid:0)

E (cid:0) )

(

)

(

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

trong đó, EvEn () là thông lượng phụ thuộc năng lượng (energy- dependent flux or flux per unit energy)

(cid:0) (cid:0)

dEE )

(cid:0) ()

(

dEE )

(

F a

F s

(cid:0) (cid:0) (cid:0) (cid:0) 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ự (cid:0) () (cid:0) (cid:0)

Đị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.

Định luật Fick (2)

(cid:0)

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

(cid:0) (x)

J x

d dx

(cid:0) (cid:0)

Jx là số neutron thực đi qua trên mỗi đơn vị thời gian 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 coefficient) [cm]

Định luật Fick (3)

Mật độ neutron tăng

Dòng

Gradient

J (r) = (cid:0) D(cid:0) (cid:0) (r)

(cid:0)

J (r)

(cid:0) (r)

Định luật Fick (4)

(cid:0)

Trong không gian 3 chiều và các neutron đơn năng

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

J

Dgrad

Định nghĩa toán tử gradient (tọa độ vuông góc)

J.n

nJ(cid:0)

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)

Đị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.

Đị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).

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 neutron production in volume V

Rate of neutron absorption in volume V

Rate of neutron leakage from V

Rate of change in the number of neutrons in a small volume V

- -

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à

ndV V

d dt

(cid:0)

The rate of change

(cid:0)

n t

(cid:0) hay (cid:0) (cid:0)

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

(cid:0)

(cid:0) (cid:0)

(cid:0)  Tốc độ mất neutron do hấp thụ trong mỗi cm3/s:

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

(cid:0) a

(cid:0) (cid:0) (cid:0)

Phương trình liên tục (4)

J.n

 Để ý 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à

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

(cid:0) (cid:0)

J.n A

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.n A

div dV J V

(cid:0) (cid:0) (cid:0)

Leakage rate

dVJdiv V

(cid:0) (cid:0)

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:

(cid:0)

sdV V

Jdiv V

n t

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

Dạng tổng quát

(cid:0)

Jdiv

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

n t

(cid:0)

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à:

0(cid:0)

Jdiv

(cid:0) (cid:0) (cid:0)

Phương trình khuếch tán (1)

(cid:0)

Nhớ lại định luật Fick:

(cid:0) (cid:0)

J

(cid:0) 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:

(cid:0)

(với neutron đơn năng và hệ số khuếch tán là hằng số)

(cid:0)2 (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) D s (cid:0) n t

= div grad is called Laplacian (cid:0)

Phương trình khuếch tán - laplacian (2)

Phương trình khuếch tán (3)

(cid:0)

Vì (cid:0) (cid:0) = nv

(cid:0)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

1 v

(cid:0) 2

(cid:0)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) 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): (cid:0) a

s D

1 2 L

(cid:0)2

(cid:0) (cid:0) (cid:0) (cid:0)

là diện tích khuếch tán [cm2] (cid:0)

L là độ dài khuếch tán [cm]

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.

Điều kiện biên (2)

 Thông lượng neutron phải thực, không âm và hữu hạn

Vacuum or air

Độ dài ngoại suy (?)

Diffusion theory

Surface

Exact flux

d (cid:0)

 Neutron không trở lại từ bề mặt tự do (chân không)

(cid:0)71,0(cid:0) 3(cid:0) 13,2(cid:0)

Diffusing medium

Thực nghiệm cho thấy các vật liệu không phải 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ệ

Điều kiện biên (3)

(cid:0)

 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ạ)

(cid:0)

(

) nA

) nB

(cid:0)

Các điều kiện biên này đôi khi được gọi là interface boundary conditions

 Tính đối xứng

Nghiệm của phương trình khuếch tán (1)

 Nguồn tấm phẳng vô hạn

(cid:0) 2

(cid:0)

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

s D

(cid:0)

(cid:0) (cid:0) (cid:0) (cid:0)

d dx (cid:0) 2 d dx

1 2 L 1 2 L

(cid:0) (cid:0) (cid:0)

x=0 plane

Điều kiện biên: thông lượng là hữu hạn khi x tiến đến vô cùng

Nghiệm của phương trình khuếch tán (2)

 Nguồn tấm phẳng vô hạn (tt)

(cid:0)

(cid:0) (cid:0) (cid:0) Nghiệm tổng quát:

(cid:0)

LxAe

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 (cid:0) (cid:0)

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

(Gọi là điều kiện nguồn)

xJ )(

lim x 0

S 2

Tính đối xứng

(cid:0) (cid:0)

Nghiệm của phương trình khuếch tán (3)

 Nguồn tấm phẳng vô hạn (tt)

(cid:0)

Từ định luật Fick,

Lxe

Dùng điều kiện nguồn

d dx

(cid:0) (cid:0) (cid:0) (cid:0)

(cid:0)

Lxe

DA L SL D 2 SL (cid:0) D 2

(cid:0)

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

(cid:0)

Lxe | |

SL D 2

(cid:0) (cid:0)

Nghiệm của phương trình khuếch tán (4)

 Nguồn điểm

(cid:0)

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,

(tại r khác zero)

d dr

1 2 r

1 2 L

(cid:0) (cid:0)

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à

rJr(cid:0) 4 )(

rJr )(

lim 2 r 0

S (cid:0)4

(cid:0) (điều kiện nguồn) (cid:0)

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

(cid:0) 2

(cid:0)

(cid:0) 1 L

(cid:0) d dr

(cid:0) (cid:0)

Nghiệm tổng quát

(cid:0)

(cid:0) (cid:0) (cid:0)

(cid:0)

(cid:0) (cid:0)

Nghiệm của phương trình khuếch tán (6)

 Nguồn điểm (tt)

(cid:0)

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

Lre

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

1 rL

d dr

1 2 r

(cid:0) (cid:0)

S (cid:0)4 D

(cid:0)

rJr )(

lim 2 r 0

S (cid:0)4

(cid:0) (cid:0)

(cid:0)

Se (cid:0) 4

(cid:0)

Nghiệm của phương trình khuếch tán (7)

(cid:0)

 Tấm trần

(

)

(

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) 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: (cid:0) da )

Đối với x dương

(cid:0) (cid:0) (cid:0)

Lda )

(

Lda )

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

) (2

Ce ( Ae Lda )

(cid:0) Ae da(cid:0) ( C Ae

Điều kiện biên

(cid:0) (cid:0) (cid:0)

Solutions of the diffusion equation (8)

 Tấm trần (tt)

Nghiệm phương trình trở thành

)

(cid:0) (cid:0) (cid:0)

(cid:0)Lda

(cid:0)

(cid:0) (cid:0) (cid:0)

Sử dụng điều kiện nguồn (tấm phẳng vô hạn)

Lda )

(cid:0) (cid:0) (cid:0)

(cid:0) 1

(cid:0) (cid:0) (cid:0)

Lda )

(cid:0)

Lda )

SL D 2 SL e D 2

e e

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

Lx |

Lda )

(cid:0)

Lda )

SL D 2

e e

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) Do đối xứng, ta có (cid:0) (cid:0) (cid:0)

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 ta được:

(

Lxda |

(cid:0)

Lda )

(

Lda )

( (cid:0) Ldae e e

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

e sinh

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

SL D 2 SL D 2

( cosh

da (

x | |) (cid:0)Lda )

(cid:0) (cid:0)

sinh và cosh là các hàm hyperbol được định nghĩa:

sinh

cosh

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

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

Nghiệm của phương trình khuếch tán (11)

(cid:0) 2

(cid:0)

 Bài toán nhiều vùng (tt) Ta có phương trình khuếch tán cho các vùng



for

1 2

1 2 L 1

d dx (cid:0) 2

(cid:0)

(cid:0) (cid:0) (cid:0)



for

2 2

d dx

1 2 L 2

(cid:0) (cid:0)

(cid:0) (cid:0) (cid:0) Điều kiện biên (cid:0) finite

xJ )(

(source condition)

S 2 (cid:0)

(cid:0) (cid:0)

lim x 0 (cid:0)

(

)2/

(

continuity

the flux

at the interface

)2/ (cid:0)

(cid:0) (cid:0) (cid:0)

continuity

current

at the interface

D 1

(cid:0)

a (cid:0) d 1 dx

d 2 dx

(cid:0) (cid:0) (cid:0) (cid:0)

Nghiệm của phương trình khuếch tán (12)

(cid:0)

 Bài toán nhiều vùng (tt) Nghiệm tổng quát

x )(

)

sinh(

)

Lx / 1

A 1

(cid:0) (cid:0)

cosh( Lx / 2

C 1 Lx / 2

(cid:0)

x )(

eA 2

eC 2

(cid:0) (cid:0) (cid:0)

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:

Lx / 1

)0(

D 1

(cid:0) d 1 dx

S 2

AD 1 1 L 1

CD 1 1 L 1

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

C 1

SL 1 D 2 1

(cid:0) (cid:0) (cid:0)

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:

L 2

cosh

sinh

A 1

eA 2

a L 2 1

(cid:0) (cid:0) (cid:0) (cid:0)

L 2

sinh

cosh

AD 2 2 L

SL 1 D 2 1 S 2

AD 1 1 L 1

a L 2 1

a L 2 1 a L 2 1

2 ph/tr chứa 2 ẩn số. Giải cho A1 và A2

(cid:0) (cid:0) (cid:0) (cid:0)

A 1

(cid:0) (cid:0)

a a

cosh( cosh(

sinh( sinh(

LD 2 1 LD 1 2

L )2/ 1 L )2/ 1

LD 2 1 LD 1

L )2/ 1 L )2/ 1

L 1

(cid:0)

A 2

(cid:0)

SL 1 D 2 1 LSL 1 2

cosh(

sinh(

LD 1 2

LD 1

L )2/ 1

e L )2/ 1

(cid:0)

Do đó ta biết được thông lượng ở 2 vùng.

Độ dài khuếch tán (1)

Neutron emitted here

Neutron absorbed here

Actual path of neutron

Để 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

Độ dài khuếch tán (2)

(cid:0)

Se (cid:0) 4

(cid:0)

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 Gọi dn là số neutron bị hấp thu trong khoảng r và r+dr

dVr

)((cid:0)

(cid:0) (cid:0)

4(cid:0)

(cid:0)

S D

S 2 L

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

Độ dài khuếch tán (3)

Lr(cid:0)

Xác suất để một neutron từ nguồn bị hấp thụ trong khoảng dr

rp )(

(cid:0)

1 2 L 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

(cid:0)

2 rpr )(

(cid:0) (cid:0)

3 er

2 L

0 1 2 L

L (cid:0)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

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).

Phương pháp khuếch tán nhóm (2)

Group 1

Group 2

Group g

Group N-1 Group N

In c re a s in g E n e rg y

Các nhóm năng lượng cho tính toán khuếch tán nhóm

Phương pháp khuếch tán nhóm (3)

(cid:0)

Thông lượng neutron nhóm g

((cid:0)

dEE)

(cid:0) (cid:0)

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ạ

(

(cid:0) ()

dEE )

(cid:0) (cid:0) Absorption rate (cid:0)

(

(cid:0) ()

dEE )

(cid:0) (cid:0)

(cid:0)

(

dEE )

(cid:0) (cid:0)

(cid:0)

(cid:0) (cid:0) Absorption rate

Phương pháp khuếch tán nhóm (4)

g (cid:0)

(cid:0)

Gọi tiết diện chuyển nhóm từ nhóm g sang nhóm h

(cid:0) gh

(cid:0) (cid:0) (cid:0) (cid:0) Total transfer rate out of g

(cid:0) (cid:0)

(cid:0)

(cid:0) hg

(cid:0) (cid:0) (cid:0) (cid:0) Total transfer rate into g (cid:0)

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

(cid:0) 2

(cid:0)

(cid:0) gh

(cid:0) hg

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

gD là hệ số khuếch tán nhóm được cho bởi

(cid:0)

(cid:0) (cid:0) (cid:0)

ED (

(cid:0) ()

dEE )

(

dEE )

(cid:0)

1 (cid:0) 3

1 (

)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

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.

(cid:0)

Khuếch tán neutron nhiệt (1)

kTEeE

En (

)

(cid:0) (cid:0)

n 23)

(cid:0) (

Phân bố năng lượng của neutron nhiệt được cho bởi hàm Maxwellian (cid:0) 2 kT

(cid:0)

Thông lượng phụ thuộc năng lượng của neutron nhiệt

(

)

EvEn ()

(

)

(cid:0)

kTEEe

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

2 m

(cid:0) 2 kT

n 23 )

(cid:0) (

(cid:0) (cid:0)

Ev (

)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) Bằng cách sử dụng công thức

E 2 m

(cid:0) (cid:0)

Khuếch tán neutron nhiệt (2)

(cid:0)

dEE)

(cid:0) Kí hiệu thông lượng neutron nhiệt một nhóm là (cid:0) (cid:0) ((cid:0) (cid:0)

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 ((cid:0) 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

kTE

(cid:0)

(cid:0) 2 kT

2 m

n 23 )

(cid:0) (

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

n 2 (cid:0)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

kT 2 m

(cid:0) (cid:0)

Khuếch tán neutron nhiệt (3)

mv (cid:0)2

1 2

(cid:0) Sẽ rất thuận tiện nếu thay và

Giá trị của ET và vT có thể được tính từ các công thức sau

617

eV 510

(cid:0) (cid:0) (cid:0)

cm/s

284

(cid:0)

(cid:0) (cid:0) T[K]

2 (cid:0)

(cid:0)

Khuếch tán neutron nhiệt (4)

T và thông lượng tại 2200 m/s (cid:0)

0 = nv0

0 cho(cid:0) (cid:0)

(cid:0)

Quan hệ giữa (cid:0) có thể tìm được bằng cách chia(cid:0) (cid:0)

(cid:0) (cid:0)

v0 v T

(cid:0)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

(cid:0) (cid:0)

T là tất cả thông lượng neutron

0 là giả-thông lượng được tính bằng cách

T 0 T Về ý nghĩa vật lý thì (cid:0) nhiệt gộp lại bằng cách tính tích phân phân bố Maxwell trong khi (cid:0) (cid:0) giả thiết tất cả các neutron nhiệt chỉ có một năng lượng.

T0 = 293.61 K (cid:0) (cid:0)

Khuếch tán neutron nhiệt (5)

(cid:0)

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

ED (

(cid:0) ()

dEE )

(

dEE )

(cid:0)

1 (cid:0) 3

1 (

)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

(

(cid:0) ()

dEE )

(cid:0)

T Tích phân là tốc độ hấp thụ toàn phần và bằng

Tiết diện hấp thụ một nhóm neutron nhiệt 1 (cid:0) (cid:0) (cid:0) (cid:0)

(

Tg )( a

(cid:0)E 0 )

với ga(T) h/s non-1/v hay h/s Westcott (cid:0)

(cid:0)

(

)

Tg )( a (cid:0)

(cid:0) (cid:0) (cid:0)

(

)

Tg )( a

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

T 0 T

(cid:0) (cid:0)

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

Moderator

TL ,

, 0197

(cid:0)

4,9 (cid:0)

(cid:0) (cid:0)

H2O D2O Be Graphit

00,1 10,1 85,1 60,1

Density, g/cm

D , cm 16,0 87,0 50,0 84,0

a ,0 3,9 04,1 4,2

cm(cid:0) 510 310 410

TL , cm 85,2 97 21 59

2cm 1,8 310 480 3500

[Argonne National Laboratory ANL-5800, 1963]

(cid:0) (cid:0)

Khuếch tán neutron nhiệt (7)

(cid:0)2

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

(cid:0) Ta

(cid:0)

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

s T D

1 2 L T

(cid:0)2 L T

(cid:0) (cid:0) (cid:0) (cid:0)

(cid:0) là diện tích khuếch tán nhiệt (thermal diffusion area);

LT là độ dài khuếch tán nhiệt (thermal diffusion length)

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.

Tính toán làm chậm neutron hai nhóm (2)

(cid:0) 2

(cid:0)

(cid:0) Nhớ lại phương trình khuếch tán nhóm tổng quát N

(cid:0) gh

(cid:0) hg

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

x khác 0, S=0

Giả sử (cid:0)

(cid:0)D

a = 0 (good moderator) 0

(cid:0) 11

(cid:0) (cid:0) (cid:0) (cid:0)

(cid:0) 2

(cid:0)

(cid:0) 11 D

1 (cid:0) L T 2 T

(cid:0) (cid:0) (cid:0) (cid:0)

21(cid:0)

(cid:0) (cid:0) (cid:0)

Tính toán làm chậm neutron hai nhóm (3)

T từ phương trình khuếch tán nhiệt, đầu tiên

(cid:0)

Để tìm (cid:0) cần giải phương trình khuếch tán nhanh để có (cid:0)

(cid:0)

(cid:0) (cid:0)

d 1 dr

1 2 r

(cid:0) D

(cid:0)

Dùng laplacian của tọa độ cầu d dr

(cid:0)

là tuổi neutron (age) [cm2] (cid:0)

(cid:0)

Se 4(cid:0)

rD 1

(cid:0)

Lr T

Tương tự như trong trường hợp giải khuếch tán neutron từ một nguồn điểm

(cid:0)

(

)

(cid:0)

(cid:0) (cid:0) (cid:0) (cid:0)

(cid:0) 4

)

2 SL T 2 LDr ( T

(cid:0)

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ó

(cid:0)

J x

sinh

aJ )(

(cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0) (cid:0)

d dx

SL D 2

( cosh

da (

Lx ) Lda )

sinh