Monte Carlo neutron slab model


by The_shadow
Tags: carlo, model, monte, neutron, slab
The_shadow
The_shadow is offline
#1
Dec16-12, 01:11 PM
P: 19
Hi All,

Trying to write a simple Monte Carlo code for mono energetic neutrons in a steady state slab geometry in one-dimension. Their is an incident source on the left hand side of the slab and its of length L.

I understand how to work out the path length and scattering angles based on the pseudo-random numbers, aswell as the type of scattering the particle undergoes, however I wish to calculate the flux and current and compare it to the deterministic model. I'm struggling to find much on this topic (that's useful).

For example the deterministic solutions for flux and current depend upon the source term (directly proportional), by contrast the Monte Carlo depends upon the initial numbers of particles, and a sufficiently large number is required to get good statistics - so presumably some scaling/normalization is required to compare directly to the deterministic equations? anyone know what that scaling factor is?
Phys.Org News Partner Science news on Phys.org
SensaBubble: It's a bubble, but not as we know it (w/ video)
The hemihelix: Scientists discover a new shape using rubber bands (w/ video)
Microbes provide insights into evolution of human language
mathman
mathman is offline
#2
Dec16-12, 03:45 PM
Sci Advisor
P: 5,941
Divide by the number of particles in the Monte Carlo to get results per particle. You should have a conversion factor for source strength to number of particles.
The_shadow
The_shadow is offline
#3
Dec17-12, 04:54 PM
P: 19
Hi Mathman! many thanks, I'm working on an Event-based Monte Carlo code - i.e the tracking of individual particles in a medium so I'm not sure where/how I derive the conversion factor for source strength to number of particles in order to compare it to the deterministic equations.

For example the current at a given surface is just the sum of net particles passing through that particular surface divided by the surface area and the flux is just the sum of track lengths divided by the volume of the element. The number of particles needs to be high enough (millions) to get good statistics - however that makes the flux and current completely arbitrary so how exactly would one compare this to the standard analytical solution for a neutron diffusion in a slab? i.e say:

{d^2\phi(x) \over d x^2} -\Sigma_a\Sigma_t\phi(x) + S = 0

Where S is the incident source. This source would produce N particles say, but given a source of say unity, what is N? Many thanks again,

mathman
mathman is offline
#4
Dec18-12, 05:55 PM
Sci Advisor
P: 5,941

Monte Carlo neutron slab model


What is your definition of source strength? When I was working in this area (a long time ago), the source strength was simply the number of neutrons, usually in the form of an energy spectrum.

Note:the equation looks like a Latex expression, without saying Latex.
Aaronvan
Aaronvan is offline
#5
Dec24-12, 04:57 PM
P: 36
Not sure if this will help but here you are.

http://people.sc.fsu.edu/~jburkardt/...imulation.html


Register to reply

Related Discussions
Ising Model Monte Carlo -- interpretations Classical Physics 1
Using a discrete Monte-Carlo technique in a multi-variable model Set Theory, Logic, Probability, Statistics 2
Monte-Carlo simulation for a matrix model General Physics 1
Monte Carlo Sim Science & Math Textbook Listings 1
Monte Carlo and Inverse Monte Carlo - Good Sources for Self Study? General Physics 1