Monte Carlo Methods for Isotopic Inventory Analysis

1
Monte Carlo Methods for Isotopic Inventory
Analysis

• Professor Paul Wilson
• Phiphat Phruksarojanakun (Nor)

2
Outline

• Problems and Motivations
• Methodology
• Improvements on Principle
• Results
• Future Work

3
Problems and Motivations

• Restrictions to Traditional Methods
• Static System

4
Problems and Motivations

• Future Nuclear System
• Dynamic
• Streams of Materials

Heat Exchanger
Reactor
5
Methodology

• 0-D Steady State Calculation
• A control volume

6
Methodology

• 0-D Steady State Calculation
• A control volume

Neutron Flux
7
Methodology

• 0-D Steady State Calculation
• A control volume

Neutron Flux
tr
8
Methodology

• 0-D Steady State Calculation
• A control volume
• An initial isotope

9
Methodology

• 0-D Steady State Calculation
• A control volume
• An initial isotope
• Atomic, Mass and Isomeric State ? KZA
• KZA (A1000M)10I

10
Methodology
tr
11
Methodology
12
Methodology

• Pseudo Code
• Determine if there is a reaction
• before the atom exits.

• Calculation
• ? 1 e-?t 1 e-N
• N -ln(1 ?)
• ? ?decay ?sf
• Nr ? tr

13
Methodology

• Pseudo Code
• Determine if there is a reaction
• before the atom exits.

• Calculation
• ? 1 e-?t 1 e-N
• N -ln(1 ?)
• ? ?decay ?sf
• Nr ? tr

14
Methodology

• Pseudo Code
• Determine if there is a reaction
• before the atom exits.

• Calculation
• ? 1 e-?t 1 e-N
• N -ln(1 ?)
• ? ?decay ?sf
• Nr ? tr

. . .
1
2
3
Nr
1/ ?
15
Methodology

• Pseudo Code
• Determine if there is a reaction
• before the atom exits.
• If N gt Nr
• -Record the presence of the atom
• -An atom exits
• -N N - Nr

• Calculation
• ? 1 e-?t 1 e-N
• N -ln(1 ?)
• ? ?decay ?sf
• Nr ? tr

16
Methodology

• Pseudo Code
• Determine if there is a reaction
• before the atom exits.
• If N lt Nr
• - Record the presence of the atom
• - Sample a new atom
• - tr tr N/ ?

• Calculation
• ? 1 e-?t 1 e-N
• N -ln(1 ?)
• ? ?decay ?sf
• Nr ? tr

17
Methodology

• Pseudo Code
• Determine if there is a reaction
• before the atom exits.
• If N lt Nr
• - Record the presence of the atom
• - Sample a new atom
• - tr tr N/ ?

• Calculation
• ? 1 e-?t 1 e-N
• N -ln(1 ?)
• ? ?decay ?sf
• Nr ? tr

18
Improvements on Principle

• 0-D Calculation
• Simple Flow
• Complex Flow
• Arbitrary Source Sampling
• Loop

19
0-D Benchmark

• MATLAB
• 32000 NPS
• 10 isotopes
• t1/2 1 s
• Tally at discrete locations in CV

20
0-D Benchmark
21
0-D Benchmark

• C
• O-14 pure decay, 1010 NPS, Error 0.001

22
0-D Benchmark
12.04
7.83
0.89
3.29
23
Flow Benchmark
24
Flow Benchmark
25
MPI Performance
26
Future Work

• Material Sampling
• Variance Reduction Techniques
• Force Reaction
• Uniform Sampling
• User Input Interface

27
Thats it !!!

• QUESTIONS???