#/usr/env python from pyne.material import Material, MaterialLibrary, MultiMaterial import os matlib = MaterialLibrary() Lead = Material({'Pb': 1.00}) Lead = Lead.expand_elements() Lead.name = "Lead" Lead.density = 11.36 matlib["Lead"] = Lead print Lead Iron = Material({'Fe': 1.00}) Iron = Iron.expand_elements() Iron.name = "Iron" Iron.density = 7.87 matlib["Iron"] = Iron print Iron #espi metals ss316 = Material({'Fe': 65.34, 'Cr': 17, 'Ni': 12, 'Mo': 2.5, 'Mn': 2.0, 'Si':1.00, 'C12':0.08, 'P':0.05, 'S': 0.03 }) ss316 = ss316.expand_elements() #ss316.name = "ss316" ss316.name = "m1" ss316.density = 8.03 #matlib["ss316"] = ss316 matlib["m1"] = ss316 print ss316
230000:-2.0000000000E-05, 741820:-8.0100000000E-06, 741830:-4.3020000000E-06, 741840:-9.2160000000E-06, 741860:-8.5920000000E-06, 400910:-1.1019972084E-01, 400920:-1.6844253230E-01, 400940:-1.7070152836E-01, 400960:-2.7500821600E-02}) m1.normalize() m2.normalize() m3.normalize() m5.normalize() m1.density = -10.44 m2.density = -1.106 m3.density = -2.66 m5.density = -6.56 # from pixels, re-normalize m1_frac = 0.488 m2_frac = 0.238 m3_frac = 0.056 m5_frac = 0.100 #tot = m1_frac + m2_frac + m3_frac + m5_frac #m1_frac = m1_frac / tot #m2_frac = m2_frac / tot #m3_frac = m3_frac / tot #m5_frac = m5_frac / tot
#!/usr/bin/python # from pyne.material import Material, MaterialLibrary print "Welcome!" mat_lib = MaterialLibrary() # # define iron for a shield ironvec = {260000: 1} # pure iron iron = Material() iron.density = 7.86 iron.from_atom_frac(ironvec) iron = iron.expand_elements() # to write out an mcnp material card iron.metadata['mat_number'] = 101 iron.write_mcnp('twoslab_mcard.txt', 'atom') # # # define a simple water since O-18 not in mcnp xs libs watervec = {10010000: 2, 80160000: 1} # simple water water = Material() water.density = 1.0 water.from_atom_frac(watervec) # to write out an mcnp material card water.metadata['mat_number'] = 102 water.write_mcnp('twoslab_mcard.txt', 'atom') # # define a low density simple water since O-18 not in mcnp xs libs watervec = {10010000: 2, 80160000: 1} # simple water lowdensitywater = Material() lowdensitywater.density = 0.9 lowdensitywater.from_atom_frac(watervec)
'Si': 0.00026, 'P': 2e-05, 'S': 3e-05, 'Ti': 1e-05, 'V': 0.002, 'Cr': 0.09, 'Mn': 0.0055, 'Co': 5e-05, 'Ni': 0.0001, 'Cu': 3e-05, 'Nb': 5e-05, 'Mo': 3e-05, 'Ta': 0.0012, 'W': 0.011, }) eurofer_2016.density = 7.87 eurofer_2016 = eurofer_2016.expand_elements() #this needs changing to be the 2017 version eurofer = Material({ 'Fe': 0.88821, 'B': 1e-05, 'C': 0.00105, 'N': 0.0004, 'O': 1e-05, 'Al': 4e-05, 'Si': 0.00026, 'P': 2e-05, 'S': 3e-05, 'Ti': 1e-05, 'V': 0.002,