Scribd
Upload
899 views

Monte Carlo Simulations of Radioactivity Using Python

In this presentation Monte Carlo simulation of radioactivity is done with the help if Python programming.

Uploaded by

Harsh Sharma
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PPTX, PDF, TXT or read online on Scribd
899 views

Monte Carlo Simulations of Radioactivity Using Python

In this presentation Monte Carlo simulation of radioactivity is done with the help if Python programming.

Uploaded by

Harsh Sharma
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PPTX, PDF, TXT or read online on Scribd
You are on page 1/ 12

Numerical Methods

Seminar

Monte Carlo Methods for


Simulation Of Radioactivity
What are Monte Carlo Methods?

 Monte Carlo methods are a broad class of computational algorithms that


rely on repeated random sampling to obtain numerical results.

 This technique is based on the random numbers and probability


statistics to investigate various problems in diverse areas.

 Used in physical and mathematical problems and are most useful when it
is difficult or impossible to use other approaches.

 Monte Carlo methods are mainly used in : optimization, numerical


integration, and generating draws from a probability distribution.
Usage in Physics

 Monte Carlo methods are useful for simulating systems with many coupled
degrees of freedom, such as fluids, disordered materials, strongly coupled
solids, and cellular structures, interacting particle systems , statistical
physics (Ising model). Evaluation of multidimensional definite integrals
with complicated boundary conditions.
Radioactivity Simulation using Monte
Carlo Methods
 What is radioactivity?
 The spontaneous decomposition of a nucleus into other nucleus with emission
of radiations is called radioactivity.
 In a radioactive source containing a very large number of radioactive nuclei, it is
not possible to predict when any one of the nuclei will decay. Although the decay
time for any one particular nucleus cannot be predicted, the average rate of
decay of a large sample of radioactive nuclei is highly predictable.
 Radioactive nuclei that emit α,β, and γ radiations. It is not possible to predict
when any one radioactive nucleus will decay and emit a particle.
 Since radioactive materials contains a large number of nuclei, it is possible to
predict the average rate of decay with high probability.
 There are a predictable number of decays per second even though it is not
possible to predict which nuclei among the sample will decay.
Simulation of Radioactivity using Monte
Carlo Methods
Algorithm:
Step 1:START
Step 2:Define No /*Total number of nuclei*/
time/*Total time of study*/
dT/*Time step*/
P/*Probability of decay*/
Initialize clock T=0
Step 3:Set Nundecayed=No/*Number of undecayed nuclei equal to No*/
Step 4: Loop: I to Nundecayed
4.1:Generate a random number x
Step 5:If |x-P|<0.01 and T<time:No =No -1
Update Clock:T=T+dt
Step 6:If No >0 and T<time:GOTO Step 3
Else:END
Python Program
Inputs:
Output:
Actual Graph of Radioactive decay
Thank You!

You might also like