REACTOR_SIMULATION
Numerical Analysis Library

REACTOR_SIMULATION is a FORTRAN90 program which simulates the shielding effect of a slab intended to absorb or reflect most of the radiation from a neutron source. This program is intended as a simple demonstration of Monte Carlo simulation as applied to the analysis of complicated probabilistic systems.

REACTOR_SIMULATION is an example program that was developed in the book "Numerical Methods and Software".

Related Data and Programs:

HODGE, a C program which implements a 2D cellular automaton, that can be regarded as a model of the spread of an illness, and whose parameters can be tuned to exhibit stability, regular waves, or a variety of chaotic behavior. This is a simplified version of a program by Martin Gerhardt and Heike Schuster

ISING, a FORTRAN90 program which simulates a 3D Ising model of local variations in magnetism;

LIFE_OPEN_GL, a C program which uses OpenGL to display the evolution of John Conway's "Game of Life", by Simon Green.

MD1, a FORTRAN90 program which carries out a simple molecular dynamics simulation;

MD2, a FORTRAN90 program which carries out a simple molecular dynamics simulation;

MD3, a FORTRAN90 program which carries out a simple molecular dynamics simulation;

MD3GLUE, a FORTRAN90 program which carries out a simple molecular dynamics simulation;

MDDISK, a C program which simulates a fluid modeled as a collection of soft disks, by Dennis Rapaport;

NMS is a FORTRAN90 library which accompanies the book "Numerical Methods and Software", from which the REACTOR_SIMULATION program was extracted.

UNIFORM is a FORTRAN90 library which generates pseudorandom uniform numbers.

XFIRES, a C program which models the occurrence of fires in a forest, and displays the status of the forest using X Windows.

XISING, a C program which models the variations in ferromagnetism in a material, displaying the results using X Windows.

XWAVES, a C program which simulates the behavior of solutions of certain forms of the wave equation, displaying the results using X Windows.

Reference:

1. Paul Bratley, Bennett Fox, Linus Schrage,
   A Guide to Simulation,
   Second Edition,
   Springer, 1987,
   ISBN: 0387964673,
   LC: QA76.9.C65.B73.
2. David Kahaner, Cleve Moler, Steven Nash,
   Numerical Methods and Software,
   Prentice Hall, 1989,
   ISBN: 0-13-627258-4,
   LC: TA345.K34.
3. Pierre LEcuyer,
   Random Number Generation,
   in Handbook of Simulation,
   edited by Jerry Banks,
   Wiley, 1998,
   ISBN: 0471134031,
   LC: T57.62.H37.
4. Peter Lewis, Allen Goodman, James Miller,
   A Pseudo-Random Number Generator for the System/360,
   IBM Systems Journal,
   Volume 8, Number 2, 1969, pages 136-143.
5. George Marsaglia, Wai Wan Tsang,
   A fast, easily implemented method for sampling from decreasing or symmetric unimodal density functions,
   SIAM Journal of Scientific and Statistical Computing,
   Volume 5, Number 2, June 1984, pages 349-359.

Source Code:

• reactor_simulation.f90, the source code;
• reactor_simulation.csh, commands to compile, link and load the source code;
• reactor_simulation_output.txt, the output file from a sample run.

List of Routines:

• MAIN is the main program for the reactor shielding simulation.
• ABSORB returns TRUE if the particle is absorbed upon colliding.
• CROSS returns the "cross section" of a particle based on its energy.
• DIST2C returns the distance to collision.
• ENERGY assigns an energy to an emitted particle.
• OUTPUT prints the results of the reactor shielding simulation.
• SCATTER returns the new direction and energy of a particle that is scattered.
• SOURCE generates a new particle from the neutron source.
• TIMESTAMP prints the current YMDHMS date as a time stamp.
• UPDATE determines the position of the particle after it has traveled D units.

You can go up one level to the FORTRAN90 source codes.