#
# 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])
# 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/>.
u4_c = 0.1771
u5_f = 1.219
u5_c = 0.9519e-1
u8_f = 0.2998
u8_c = 0.7019e-1
print 'u5: ', u5_f / (u5_f + u5_c)
print 'u8: ', u8_f / (u8_f + u8_c)
from pirs.mcnp import Material
u = Material('U')
a4 = u.how_much(1, 92234).v
a5 = u.how_much(1, 92235).v
a8 = u.how_much(1, 92238).v
print a5, a8
print u.report()
kinf = (a4 * u4_f + a5 * u5_f + a8 * u8_f) / (a4 * (u4_f + u4_c) + a5 *
(u5_f + u5_c) + a8 *
(u8_f + u8_c))
print kinf
# Check "emply" mixture as void.
from pirs.mcnp import Material
v = Material() # This should be void. It can be added only with vol. fractions
u = Material('U')
u.dens = 18.0
uu = Material(u, (1, 3), v, (1, 3))
print(uu.report())