Khóa luận: Xác Định Bề Dày Vật Liệu Bằng Mô Phỏng Monte Carlo

Khám phá phương pháp đo bề dày vật liệu chính xác qua mô phỏng Monte Carlo vật lý, ứng dụng hiệu quả trong phân tích và kiểm định.

• Tác giả: Dương Thái Đương
• Người hướng dẫn: ThS. Hoàng Đức Tâm
• Trường học: Đại học Sư Phạm Tp.HCM
• Chuyên ngành: Vật lý hạt nhân
• Đề tài: Xác định bề dày vật liệu bằng phương pháp mô phỏng Monte-Carlo
• Loại tài liệu: Luận văn tốt nghiệp
• Năm xuất bản: 2014
• Địa điểm: Thành phố Hồ Chí Minh

Mục lục chi tiết

• LỜI CẢM ƠN
• MỤC LỤC
• LỜI NÓI ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. TAN XẠ COMPTON
• 1.1.1. Tan xạ Compton
• 1.2. SỰ HẤP THỤ TIA GAMMA QUA MÔI TRƯỜNG VẬT CHẤT
• 1.2.1. Định lý Beer-Boucher
• 1.3. PHƯƠNG PHÁP MONTE CARLO
• 1.3.1. Phương pháp Monte Carlo
• 1.3.2. Chương trình mô phỏng MCNP
• 2. CHƯƠNG 2: KỸ THUẬT GAMMA TÁN XẠ VÀ PHƯƠNG PHÁP XÁC ĐỊNH BỀ DÀY
• 2.1. KỸ THUẬT GAMMA TÁN XẠ
• 2.1.1. Giới thiệu kỹ thuật gamma tán xạ
• 2.1.2. Công thức giải tích xác định bề dày vật liệu
• 2.1.3. Bề dày bão hoà
• 2.2. MÔ PHỎNG KỸ THUẬT GAMMA TÁN XẠ
• 2.2.1. Mô phỏng detector, nguồn và vật liệu
• 3. CHƯƠNG 3: KẾT QUẢ
• 3.1. SỰ THAY ĐỔI CƯỜNG ĐỘ TÁN XẠ VÀO THỂ TÍCH TÁN XẠ
• 3.2. KẾT QUẢ XÁC ĐỊNH BỀ DÀY
• 3.2.1. Kết quả ứng với từng hệ đo
• 4. CHƯƠNG 4: KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN ĐỀ TÀI
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan phương pháp Monte Carlo đo bề dày vật liệu

Bài viết này giới thiệu một cách tiếp cận hiện đại để xác định độ dày vật liệu: sử dụng phương pháp Monte Carlo trong vật lý. Đây là một công cụ kiểm tra không phá hủy (NDT) mạnh mẽ, cho phép phân tích cấu trúc bên trong mà không cần can thiệp vật lý. Phương pháp này dựa trên việc mô phỏng hành vi của hàng triệu hạt bức xạ khi chúng tương tác với vật chất, từ đó suy ra các đặc tính quan trọng như bề dày. Luận văn của Dương Thái Dương (2014) là cơ sở chính cho các phân tích kỹ thuật trong bài viết, tập trung vào kỹ thuật tán xạ gamma ngược và ứng dụng phần mềm MCNP.

1.1. Phương pháp Monte Carlo Từ lý thuyết đến ứng dụng

Về bản chất, phương pháp Monte Carlo là một kỹ thuật tính toán sử dụng các số ngẫu nhiên để giải quyết các bài toán phức tạp không thể giải bằng phương trình giải tích. Trong vật lý hạt nhân ứng dụng, phương pháp này mô phỏng quỹ đạo và tương tác bức xạ vật chất của từng hạt riêng lẻ (photon, neutron, electron). Bằng cách lặp lại quá trình này hàng triệu lần, chúng ta thu được một bức tranh thống kê chính xác về hành vi tổng thể của chùm bức xạ. Đặc điểm chính của phương pháp này là thuật toán đơn giản cho từng sự kiện và sai số tỷ lệ nghịch với căn bậc hai của số lượng sự kiện mô phỏng (1/√N). Điều này làm cho nó trở thành một công cụ cực kỳ mạnh mẽ để giải quyết các bài toán mô phỏng vận chuyển hạt.

1.2. Thách thức trong kiểm tra không phá hủy NDT hiện nay

Các phương pháp kiểm tra không phá hủy (NDT) truyền thống như siêu âm hay chụp ảnh phóng xạ truyền qua gặp nhiều hạn chế khi chỉ có thể tiếp cận vật liệu từ một phía. Ví dụ, việc kiểm tra độ ăn mòn của đường ống dẫn dầu dưới lòng đất hoặc lớp vỏ lò phản ứng là bất khả thi với kỹ thuật truyền qua. Kỹ thuật gamma tán xạ ngược, được mô phỏng bằng Monte Carlo, giải quyết triệt để vấn đề này. Nó cho phép đặt cả nguồn phóng xạ và đầu dò bức xạ ở cùng một phía so với vật thể. Thách thức chính của kỹ thuật này là xử lý tín hiệu nhiễu từ tán xạ đa cấp và phông môi trường, đòi hỏi các mô hình mô phỏng số trong vật lý phải có độ chính xác cao.

II. Nguyên lý vật lý của kỹ thuật đo bề dày bằng gamma

Cơ sở của kỹ thuật này nằm ở cách các hạt gamma tương tác với electron trong vật liệu. Khi một chùm tia gamma đi vào vật chất, một phần năng lượng của nó bị mất đi do các hiệu ứng vật lý, chủ yếu là tán xạ Compton và hiệu ứng quang điện. Bằng cách đo năng lượng và cường độ của chùm tia tán xạ ngược, chúng ta có thể suy ra quãng đường mà chúng đã đi qua bên trong vật liệu, từ đó xác định được bề dày. Định luật Beer-Boucher cung cấp nền tảng toán học cho sự suy giảm cường độ này, là chìa khóa để xây dựng các công thức giải tích và mô phỏng.

2.1. Tán xạ Compton Chìa khóa của phương pháp tán xạ ngược

Hiệu ứng tán xạ Compton xảy ra khi một photon gamma va chạm không đàn hồi với một electron và bị lệch hướng, đồng thời mất đi một phần năng lượng. Năng lượng của photon tán xạ phụ thuộc trực tiếp vào góc tán xạ. Trong kỹ thuật đo bề dày, hệ thống được thiết kế với một góc cố định giữa nguồn và đầu dò (ví dụ 120° trong nghiên cứu của Dương Thái Dương). Điều này đảm bảo rằng các photon đến được đầu dò phải có một mức năng lượng xác định. Cường độ của chùm tia tại mức năng lượng này tỷ lệ thuận với số lượng nguyên tử trong vùng tương tác, và do đó, liên quan trực tiếp đến bề dày của vật liệu.

2.2. Sự suy giảm tia gamma và hệ số suy giảm tuyến tính

Khi di chuyển trong vật chất, cường độ của chùm tia gamma giảm theo hàm mũ. Hiện tượng này được mô tả bởi định luật Beer-Boucher: I = I₀ * e^(-µd), trong đó µ là hệ số suy giảm tuyến tính và d là bề dày. Hệ số µ phụ thuộc vào năng lượng của tia gamma và bản chất của vật liệu. Quá trình suy giảm tia gamma bao gồm cả tán xạ Compton, hiệu ứng quang điện và hiệu ứng tạo cặp. Trong mô phỏng, các thư viện dữ liệu hạt nhân cung cấp tiết diện tương tác chính xác cho các quá trình này, cho phép MCNP hay Geant4 tính toán chính xác sự suy giảm và tán xạ của chùm tia.

III. Hướng dẫn thiết lập mô phỏng MCNP đo bề dày vật liệu

Việc xây dựng một mô hình mô phỏng Monte Carlo chính xác đòi hỏi sự tỉ mỉ trong việc định nghĩa từng thành phần của hệ đo. Quá trình này bao gồm việc mô hình hóa hình học của nguồn, vật liệu bia và đầu dò, cũng như khai báo các đặc tính vật lý của chúng. Phần mềm MCNP (Monte Carlo N-Particle) là một công cụ tiêu chuẩn trong ngành, cho phép người dùng xây dựng các mô hình phức tạp và theo dõi hành trình của từng hạt. Việc tối ưu hóa mô phỏng để giảm thời gian tính toán và tăng độ chính xác thống kê là một bước quan trọng trong quy trình.

3.1. Mô hình hóa nguồn ¹³⁷Cs và detector nhấp nháy NaI Tl

Trong luận văn tham khảo, nguồn phóng xạ được sử dụng là ¹³⁷Cs, phát ra tia gamma đơn năng lượng 662 keV. Nguồn được mô hình hóa dưới dạng hình trụ với các thông số kích thước chính xác. Thành phần quan trọng không kém là detector nhấp nháy NaI(Tl). Đây là loại đầu dò phổ biến nhờ hiệu suất cao và khả năng hoạt động ở điều kiện môi trường bình thường. Mô hình detector trong MCNP bao gồm tinh thể NaI(Tl), lớp phản xạ Al₂O₃, vỏ nhôm và các thành phần khác. Để kết quả mô phỏng khớp với thực nghiệm, hàm đáp ứng năng lượng (FWHM) của detector được hiệu chuẩn thực nghiệm và đưa vào mô phỏng, đảm bảo phổ năng lượng mô phỏng có dạng phân bố Gauss tương tự thực tế.

3.2. Xác định vị trí đặt bia tối ưu qua thể tích tán xạ

Một trong những yếu tố quyết định độ chính xác của phép đo là vị trí đặt bia vật liệu. Nghiên cứu của Hoàng Đức Tâm và cộng sự (2013), được trích dẫn trong tài liệu, chỉ ra rằng cường độ tín hiệu tán xạ phụ thuộc mạnh vào 'thể tích tán xạ' (vùng không gian giao nhau giữa trường nhìn của nguồn và đầu dò). Vị trí tối ưu là nơi thể tích này đạt cực đại. Bằng cách sử dụng mô phỏng Monte Carlo, các nhà nghiên cứu đã tính toán và xác định được vị trí đặt bia lý tưởng, nơi mà số đếm photon tán xạ thu được tại phổ kế năng lượng là lớn nhất, từ đó giảm thiểu sai số thống kê và nâng cao độ tin cậy của kết quả đo độ dày lớp phủ.

IV. Phân tích kết quả mô phỏng và so sánh với thực tế

Sau khi chạy mô phỏng, bước tiếp theo là phân tích dữ liệu đầu ra, cụ thể là phổ năng lượng ghi nhận tại đầu dò. Quá trình này bao gồm việc xử lý phổ để loại bỏ nhiễu nền, xác định diện tích đỉnh tán xạ, và cuối cùng áp dụng công thức giải tích để tính toán bề dày. Kết quả từ mô phỏng sau đó được so sánh trực tiếp với giá trị bề dày thực tế của vật liệu mẫu để đánh giá độ chính xác của phương pháp. Nghiên cứu trên thép C45 cho thấy sự phù hợp cao giữa lý thuyết và mô phỏng.

4.1. Kỹ thuật xử lý phổ và tính toán cường độ tán xạ

Phổ năng lượng thô từ MCNP chứa cả tín hiệu từ tán xạ một lần (thông tin hữu ích) và phông nền từ tán xạ nhiều lần và môi trường. Để trích xuất thông tin chính xác, các chương trình như Genie2K được sử dụng. Phông nền được mô hình hóa bằng một hàm đa thức bậc ba và trừ đi khỏi phổ tổng. Sau đó, một hàm Gauss được khớp vào đỉnh năng lượng tán xạ (khoảng 225 keV cho góc 120° từ nguồn 662 keV) để tính diện tích đỉnh. Diện tích này chính là cường độ chùm tia tán xạ, đại lượng đầu vào quan trọng cho công thức tính bề dày. Độ chính xác của bước này ảnh hưởng trực tiếp đến kết quả cuối cùng.

4.2. Đánh giá độ chính xác Sai lệch dưới 5 so với thực tế

Trong luận văn của Dương Thái Dương (2014), phương pháp này được áp dụng để xác định bề dày các tấm thép C45 có độ dày đã biết. Các mô phỏng được thực hiện với nhiều góc đặt bia khác nhau (0°, 15°, 30°, 45°, 60°). Kết quả cho thấy sự phù hợp rất tốt: sai số tương đối giữa giá trị tính từ mô phỏng và giá trị thực tế hầu hết đều dưới 5%, một số trường hợp dưới 10%. Kết quả này khẳng định tính khả thi và độ tin cậy cao của việc kết hợp phương pháp Monte Carlo và công thức giải tích để đo bề dày vật liệu, mở đường cho các nghiên cứu thực nghiệm trong tương lai. Điều này cũng đảm bảo các tiêu chuẩn về an toàn bức xạ được mô phỏng trước khi triển khai.

V. Kết luận và hướng phát triển của kỹ thuật mô phỏng

Nghiên cứu mô phỏng xác định bề dày vật liệu bằng phương pháp Monte Carlo đã chứng minh được hiệu quả và độ chính xác cao, là một tiền đề lý thuyết vững chắc cho các ứng dụng thực tiễn. Kỹ thuật này không chỉ giải quyết được những hạn chế của phương pháp đo truyền thống mà còn mở ra nhiều hướng phát triển mới trong lĩnh vực vật lý hạt nhân ứng dụng và khoa học vật liệu. Tuy nhiên, vẫn còn những thách thức cần được khắc phục để hoàn thiện phương pháp.

5.1. Tiềm năng ứng dụng trong công nghiệp và nghiên cứu

Kỹ thuật gamma tán xạ ngược mô phỏng bằng Monte Carlo có tiềm năng ứng dụng rộng rãi. Trong công nghiệp, nó có thể được dùng để đo độ dày lớp phủ, kiểm tra độ ăn mòn của đường ống, bồn chứa mà không cần dừng hoạt động. Trong nghiên cứu vật liệu, phương pháp này giúp xác định mật độ và thành phần của các vật liệu đa lớp. Trong y học, các nguyên tắc của mô phỏng vận chuyển hạt được áp dụng trong lập kế hoạch xạ trị và thiết kế các thiết bị che chắn bức xạ, góp phần nâng cao hiệu quả điều trị và an toàn bức xạ.

5.2. Hướng phát triển Từ mô phỏng đến kiểm chứng thực nghiệm

Mặc dù kết quả mô phỏng rất hứa hẹn, hướng phát triển tất yếu là tiến hành kiểm chứng bằng thực nghiệm. Các nghiên cứu trong tương lai cần tập trung vào việc xây dựng một hệ đo thực tế dựa trên các thông số đã được tối ưu hóa mô phỏng. Ngoài ra, cần mở rộng nghiên cứu trên nhiều loại vật liệu khác nhau (nhôm, đồng, hợp kim) để kiểm tra tính tổng quát của phương pháp. Việc nghiên cứu sâu hơn về ảnh hưởng của tán xạ nhiều lần và phát triển các thuật toán xử lý phổ tiên tiến hơn cũng là một hướng đi quan trọng để tăng độ chính xác của kỹ thuật kiểm tra không phá hủy này.

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

CHƯƠNG I. Chương nảy sẽ trình bay tán xạ Compton vả phương pháp Monte Carlo. CHUONG IL KỸ THUAT GAMMA TAN XA VA PHƯƠNG PHÁP XÁC ĐỊNH BE DAY. Trong chương này chúng tôi sẽ trình bay công thức giải tích xác định bẻ day vật liệu và phương pháp mô phỏng dé xác định bẻ dày vật liệu.

CHUONG IIL KET QUA. Chương nay sẽ trình bảy toàn bộ kết quả nghiên cửu. kết luận, hưởng phát trién đẻ tài. TONG QUAN Nam 1905, Albert Einstein dựa trên y tưởng của Max Plank vẻ lượng tử nâng lượng đã đưa ra thuyết photon ánh sáng.

nhằm giải thích ba định luật quang điện mà Hertz đã phát hiện ra từ năm 1887. Cũng như nhiều ly thuyết khác, thuyết photon mới ra đời cũng có nhiều tranh cãi. Đến năm 1923, thí nghiệm tan xạ của chùm tia X lên tam graphite của Compton đã củng cổ cho tính đúng đắn của thuyết photon. Năm 1923, Athur Holly Compton, một chuyên gia tia X ở đại học Washington da làm một thi nghiệm, cho tia X đơn sắc phát ra từ đối âm cực của một catốt của dng tia X.

Chim tia X được chuẩn trực và gọi vào tắm graphite, chùm tia tán xạ được thu bởi detector đặt ở các góc khác nhau. Kết quả thu được tử thí nghiệm được cho trên hình HF Mad ~Â= 2 (I~eosØ)= 34, sin! Ø (1.1) Điều nay hoàn toan trai ngược với thuyết sóng ánh sang. Vi vậy, dé giải thích hiện tượng trên ta cần thừa nhận thuyết photon ánh sáng của Einstein. TAN XA COMPTON 1.

Tan xa Compton Các electron liên kết trong các nguyên tử graphite rất nhỏ so với năng lượng của lượng tử gamma nên có thé xem như các electron tự do. Áp dụng các phương trinh bảo toan năng lượng và động lượng ta thu được công thức (1.—h__— (12) I+ ¬ (I~cos) Trong đó £, vả £ lần lượt là năng lượng của lượng tử gamma trước va sau tán xa, m là khối lượng của electron, Hinh 1. Tan xạ Compton, Từ công thức (1.1) ta thay, AA không phụ thuộc vào bước sóng tới 4, ma chi phụ thuộc vào góc tán xạ Ø. Với góc tán xạ cảng lớn thi.

độ thay đối bước sóng càng lớn, tức năng lượng ị gamma tán xạ càng nhỏ. năng lượng của gamma tán xạ theo góc tần xạ. i Với góc tấn xạ từ 120° trở lên thi với năng H lượng ban đầu £, bat ky đều cho năng lượng tản xạ gan như là bằng nhau. Dé thị hình 1.3 vẽ cho bốn trường hợp với năng lượng ban dau lan % ? FO 42 00 SG 129 129 140 142 +60 lượt là 4MeV, 2MeV, IMeV và 0,5MeV.

Tiết điện tán xạ Compton Hình 1. Nang lượng tan xạ Compton thay đổi theo góc tán xạ Tiết diện vi phân tán xạ của hiệu ứng Compton được tinh theo công thức thực nghiệm Klein-Nishina: do, — r` l+co8 a’(1-cos0y da 2[I+z(I ~ c esoØs) ] '* Ti se0s " [ 1 + a( 1- c0 s8)] (1.3) > ˆ T đồ r= + là ban kính electron cô điển. ta 7 te Lay tích phân biểu thức (1.3) trên toan góc khôi ta được công thức tỉnh tiết diện tan xạ toàn phan [3]: »|l+a|2(l+a) 1 l l Xét hai trường hợp đặc biệt sau đây: Nếu a =5; «1, áp dụng khai triển gin đúng, công thức (1.5) Ta thấy đối với năng lượng thắp (a nhỏ), tiết điện tan xa Compton ting tuyến tinh theo sự giảm năng lượng va giảm tới giá tr giới hạn ơ,„. Nếu a = -^*_ >1, công thức (1.4) trở thành: me ø,.6) Như vậy đối với năng lượng lớn thi tiết điện tan xạ Compton đối với một electron gin như tý lệ nghịch với năng lượng của lượng tử gamma.

Đối với nguyên tử có Z electron thì tiết điện tán xạ tăng lên Z lần. đối với năng lượng của lượng tử gamma lớn thì ta có: a.7) Sự thay đổi tiết điện tán xa Compton theo góc tán xạ được biểu diễn trên hình 1.4, Tiết điện tán xạ Compton thay đối theo góc tán xạ. Với năng lượng nhỏ (khoảng vải trăm keV) thi tiết diện tán xạ sé nhỏ nhất tại góc tan xạ 90. Đôi năng lượng lớn (>1MeV) thi tiết điện tán xạ giảm đơn điệu theo sự tăng của góc tán xạ.

Mặt khác từ hình 1.4 ta cũng thấy ứng với góc tán xạ bat kỷ, với năng lượng của lượng tứ gamma tới cảng lớn thi tiết diện tan xạ cảng nhỏ. Điều nay phủ hợp với kết quá tỉnh toán ly thuyết trong các công thức (1. SỰ HAP THY TIA GAMMA QUA MOI TRUONG VAT CHAT 1. Dinh ly Beer-Boucher Lượng tử gamma tương tác với môi trường vật chất ngoài tán xạ Compton thi còn phải kẻ đến hai loại tương tác chính nữa là hiệu ứng quang điện và sinh-huỷ cặp.

Chính những tương tác đó là nguyên nhân làm suy giảm cưởng độ chùm tia gamma tới. Cường độ chùm tia gamma còn lại sau khi xuyên qua bẻ day d vật chất được tính theo công thức; l=len (1.8) Trong đó 7, /, lần lượt la cường độ chùm tia gamma sau khi xuyên qua bẻ day d và cường độ chim tia ban đầu, „ được gọi là hệ số suy giảm tuyến tính.8) được gọi là định ly Beer-Boucher. Hệ số suy giảm Hệ số suy giảm tuyến tính đặc trưng cho sự suy giảm của chim tia gamma qua môi trường vật chất. Y nghĩa cúa hệ số suy giảm tuyến tính là thông lượng chim tia gamma giảm đi ¢ lần khi đi được quãng đường rà Đại lượng P chính lả quãng đường tự do trung bình của lượng tử gamma trong môi trường vật chất.

Điều này dé thấy nếu chú ý rằng e *“ chính là xác suất dé chùm gamma xuyên qua quảng đường d ma không tương tác. Quãng đường tự đo trung bình sẽ được tính theo công thức: (1.9) Hệ số suy giảm tuyến tính „ phụ thuộc vào bản chất từng môi trường vả ca ning lượng của chim gamma tới. Nếu gọi o là tiết điện tương tác toan phân của lượng tử gamma với một nguyên tứ.10) Trong đó N là số nguyên tử trong Icm’ vật chất. Tiết diện tương tác toản phan được tính theo công thức: FF „+0 „+Ø,„ (I.

lần lượt là tiết điện tương tác của hiệu ứng quang điện, tán xa Compton, hiệu ứng sinh-huy cặp. Ngoài hệ số suy giám tuyến tính người ta còn định nghĩa hệ số suy giảm khối được tính như sau: 2 Nam, 20n= 2 TM 2 ‘g) (1.12) Ũ Trong trường hợp vật chat la một hợp chat thì hệ số suy giám khối được tinh bằng: “(E)=G + Cup, ++ Cos, = Cu (1.13) - Trong đó C, la thành phan phan trim của chất thd i. PHƯƠNG PHÁP MONTE CARLO 1. Phương pháp Monte Carlo Phương pháp Monte Carlo hay còn gọi là phương pháp thử thống ké được định nghĩa như là phương pháp tinh bảng cách biêu diễn nghiệm của bài toán đưới dạng các tham số của một đám đông lý thuyết va sử dung day số ngẫu nhiên dé xây dung mẫu đám đông mà từ đó ta thu được ước lượng thống kê của các tham số.

Nói cách khác, phương pháp Monte Carlo cung cấp những lời giải gan đúng bảng cách thực hiện các thi nghiệm lay mẫu thông kê sử dụng số ngẫu nhiên [2]. Đẻ hiểu rd hơn vé phương pháp nay ta hãy xét một vi du đơn giản sau. Tính tích phan: f={ f(x)de.14) Thông thưởng. dé tinh tích phân nay người ta chia đoạn [a,] thành n đoạn nhỏ bằng nhau vả khi đó tích phân trên được tính gan đúng bing công thức: Js 22 s(x).

trong đó x, =a+(¡-š ||°=*) " n Trong khi đó, nếu dang phương pháp Monte Carlo thi các giá trị x được chon hoàn hoàn ngẫu nhiên trong khoảng [a. Gia trị của x, được tỉnh như sau: x =a+q,(b-a) (1.15) với g, là một số ngẫu nhiên có phân bố đều trong khoảng [0,1]. Như vậy, khi tinh bằng phương pháp Monte Carlo chúng ta đã sử dụng biển ngẫu nhiền. biến ngẫu nhiên này cho ta giá trị của đại lượng ngẫu nhiên ¢, có phân bế đều trong khoảng [0,1].

Hai phương pháp tính cho kết quả càng gan nhau khi N cảng lớn. Với những phân tích trên thi phương pháp Monte Carlo có hai đặc điểm chính sau đây: Phương pháp Monte Carlo có thuật toán đơn gián, khi xây dựng thuật toán ta chi xây dựng cho một sự kiện và sau đó lập cho tất ca sự kiện còn lại. Chính vi vậy ma người ta còn gọi phương pháp Monte Carlo là phương pháp thử thống kê. Sai số của phương pháp Monte Carlo tỷ lệ với l/VX với N là số sự kiện ngẫu nhiên dùng dé mô phỏng.

Như vay, nêu muốn giảm sai số của kết quả 10 lan thi phải 12 tăng số sự kiện mô phỏng lên 100 lần, Rõ ring phương pháp Monte Carlo không phù hợp cho bài toán đôi hỏi độ chính xác cao. Tuy nhiên trong vật lý thực nghiệm nói chung và vật lý hạt nhân thực nghiệm nói riêng, đa số kết quả không đòi hỏi độ chính xác quá cao. Nguyên nhân lả vi kết quả của thí nghiệm luôn chửa đựng những sai số. Sai số này chủ yếu xuất phát tir thiết bị đo và đôi khi do cả bản chất vật lý của hiện tượng đang nghiên cứu.

Những sai số này thường lớn hơn nhiều so với sai số bắt nguồn từ phương pháp tính. Chính vi vậy ma phương pháp Monte Carlo trở thành công cụ toán học mạnh đối với các nhả vật lý, nó đã và đang được sử dụng rộng rải cả trong vật lý thực nghiệm và vật lý lý thuyết. Chương trình mô phỏng MCNP MCNP (Monte Carlo N-Partical) là chương trình ửng dụng phương pháp Monte Carlo dé mỏ phỏng các quá trình vật lý hạt nhân đối với neutron, photon va electron (các qua trình phân ra hạt nhân, tương tác cúa các hạt với môi trưởng vật chất, thông lượng neutron. Chương trình này được thiết lập rat tốt cho phép người sứ dụng xây dựng các dạng hình học phức tạp và mô phỏng dựa trên các thư viện hạt nhắn.

Sự phức tap của tương tác photon với môi trưởng vat chat cũng được xử lý trong MCNP. Chương trình điều khiển các quá trình này bằng cách gieo số thống ké theo quy luật cho trước và mỏ phỏng được thực hiện trên máy tinh vi số lần thử can thiết thưởng rat lớn. MCNP ban đầu được phát triển bởi nhóm Monte Carlo và sau này là nhỏm Radiation Transport (nhóm X-6) của phòng thi nghiệm Vật ly Lý thuyết Ứng dụng thuộc phòng thí nghiệm quốc gia Los Almos (M9). MCNP được cung cấp tới người dùng thông qua Trung tâm Thông tin Che chắn Bức xạ (RSICC) ở Oak Ridge, Tennessee (Mỹ) và ngân hàng dữ liệu OECD/NEA ở Paris (Pháp).

KỸ THUẬT GAMMA TÁN XẠ VÀ PHƯƠNG PHÁP XÁC ĐỊNH BE DAY 2. KY THUAT GAMMA TAN XA 2.