#
# PIRS-2 is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
u = Material((92235, 4, 2), (92238, 96, 2))
p = Material((94239, 90, 2), (94240, 10, 2))
o = Material(8016)
uox = u + 2 * o
pox = p + 2 * o
mox = Material((uox, 1), (pox, 1))
print mox.report()
def of(m):
a1 = m.how_much(1, ZAID=[92235, 94239])
a2 = m.how_much(1, Z=[92, 94])
return a1 / a2 - 0.10
mox.tune(of, [uox, pox])
print mox.report()
print mox.how_much(1, ZAID=[92235, 94239])
print mox.how_much(1, Z=[92, 94])
# PIRS-2 is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
# depleted u, mass fractions
ud = Material((92235, 0.2, 2),
(92238, 99.8, 2))
# simplified o for oxide
o = Material(8016)
# pu
pu = Material((94239, 93.6, 2),
(94240, 5.9, 2),
(94241, 0.4, 2),
(94242, 0.1, 2))
# oxides and mox
ux = Material((ud, 1), (o, 2))
px = Material((pu, 1), (o, 2))
mox = Material((ux, 0.5), (px, 0.5))
print mox.report()
# objective function
def obj(mix):
a1 = mix.how_much(2, ZAID=[92235, 94239, 94241, 94243])
a2 = mix.how_much(2, Z=[92, 94])
return a1/a2 - 0.025
mox.tune(obj, [ux, px])
print mox.report()