From 5a5f8fbf0ddd2c3425bc15a64e1428c1ca29d0db Mon Sep 17 00:00:00 2001 From: Gavin Ridley Date: Wed, 9 Aug 2017 11:14:15 -0500 Subject: [PATCH] added notes to controlRodWorthMatcher script --- python/controlRodWorthMatcher.py | 14 +++++++++++++- 1 file changed, 13 insertions(+), 1 deletion(-) diff --git a/python/controlRodWorthMatcher.py b/python/controlRodWorthMatcher.py index 2067a35076..d86caaa8bc 100644 --- a/python/controlRodWorthMatcher.py +++ b/python/controlRodWorthMatcher.py @@ -10,7 +10,19 @@ # as compared to the state with no control rods inserted. # All preceding columns are assumed to belong to control rods 1, 2, 3, all # the way up to n. This script will find the moltres absorbtion factor for -# each rod that most closely matches reactivities +# each rod that most closely matches reactivities. +# +# NOTE: +# This script is unfinished. I'm putting it into the python directory so +# that others have an idea of the methodology I had in mind to determine +# the rod multiplicative factor for absorbtion for systems with many rods. +# +# The general idea would be to form a residual term that is squared reactivity +# differences between detailed transport (e.g. Serpent) and Moltres' results. +# A Jacobian is formed from perturbing each rod multiplicative factor. Then +# a Newton iteration is done to drive the residual to zero. Repeat until results +# are within uncertainty of the Monte Carlo neutronics. Note that forming the Jacobian +# for this residual should probably happen in parallel on a computing cluster. # import numpy import subprocess