Can you provide more of your input? Also, this can occur if you’re running a older version of OpenMC that doesn’t have IFP implemented yet.
thw code is given below…My docker image for my mac is markmit1908/openmc-arm:latest
How do i know whether this image is ifp implemented or not…if not, then what can i do
import openmc
####pure MATERIALS
m5 = openmc.Material(name=“M5 Tubing”)
m5.add_element(‘Zr’, 0.9825, percent_type=‘wo’)
m5.add_element(‘Nb’, 0.01, percent_type=‘wo’)
m5.add_element(‘O’, 0.007, percent_type=‘wo’)
m5.add_element(‘Fe’, 0.0005, percent_type=‘wo’)
m5.set_density(‘g/cm3’, 6.55)
#m5.temperature = 622.
helium = openmc.Material(name=“Helium Gas”)
helium.set_density(‘g/cm3’, 0.0001786)
helium.add_nuclide(‘He4’, 1.0)
#helium.temperature = 622.
water = openmc.Material(name=‘Water’)
water.add_element(‘H’, 2)
water.add_element(‘O’, 1)
water.add_s_alpha_beta(‘c_H_in_H2O’)
water.set_density(‘g/cm3’, 1.0)
#water.temperature = 565.
A14B14C14D13 = openmc.Material(name=“1.87% u235”)##
A14B14C14D13.add_element(“U”, 1.0, enrichment=1.87)
A14B14C14D13.add_element(“O”, 2.0)
A14B14C14D13.set_density(‘g/cm3’, 10.53)
#A14B14C14D13.temperature = 565.
A233B233C233 = openmc.Material(name=“4.55% u235”)##
A233B233C233.add_element(“U”, 1.0, enrichment=4.55)
A233B233C233.add_element(“O”, 2.0)
A233B233C233.set_density(‘g/cm3’, 10.53)
#A233B233C233.temperature = 565.
D2 = openmc.Material(name=“4.33% u235”)##
D2.add_element(“U”, 1.0, enrichment=4.33)
D2.add_element(“O”, 2.0)
D2.set_density(‘g/cm3’, 10.53)
#D2.temperature = 565.
#mix1
Gd2O3 = openmc.Material(name=“A3_Gd2O3”)
Gd2O3.add_element(“Gd”, 2.0)
Gd2O3.add_element(“O”, 3.0)
Gd2O3.set_density(‘g/cm3’, 7.41)
A3_UO2 = openmc.Material(name=“4.55% u235”)##
A3_UO2.add_element(“U”, 1, enrichment=4.55)
A3_UO2.add_element(“O”, 2.0)
A3_UO2.set_density(‘g/cm3’, 10.53)
A3_UO2_Gd2O3 = openmc.Material.mix_materials([A3_UO2, Gd2O3], [0.98, 0.02], ‘wo’)##
A3_UO2_Gd2O3.set_density(‘g/cm3’, 10.53)
#A3_UO2_Gd2O3.temperature = 565.
#mix2
B3_UO2 = openmc.Material(name=“4.3% u235”)##
B3_UO2.add_element(“U”, 1, enrichment=4.3)
B3_UO2.add_element(“O”, 2.0)
B3_UO2.set_density(‘g/cm3’, 10.53)
B3_UO2_Gd2O3 = openmc.Material.mix_materials([B3_UO2, Gd2O3], [0.94, 0.06], ‘wo’)##
B3_UO2_Gd2O3.set_density(‘g/cm3’, 10.53)
#B3_UO2_Gd2O3.temperature = 565.
#mix3
C3_UO2 = openmc.Material(name=“4.29% u235”)##
C3_UO2.add_element(“U”, 1, enrichment=4.29)
C3_UO2.add_element(“O”, 2.0)
C3_UO2.set_density(‘g/cm3’, 10.53)
C3_UO2_Gd2O3 = openmc.Material.mix_materials([C3_UO2, Gd2O3], [0.92, 0.08], ‘wo’)##
C3_UO2_Gd2O3.set_density(‘g/cm3’, 10.53)
#C3_UO2_Gd2O3.temperature = 565.
materials = openmc.Materials([m5, helium, water, A14B14C14D13, A233B233C233, D2, Gd2O3, A3_UO2, A3_UO2_Gd2O3, B3_UO2, B3_UO2_Gd2O3, C3_UO2, C3_UO2_Gd2O3, ])
#materials.cross_sections = ‘root/endfb-viii.0-hdf5/cross_sections.xml’
materials.export_to_xml()
#boundary
Aplane1 = openmc.ZPlane(z0=-100, boundary_type=‘reflective’)
Aplane2 = openmc.ZPlane(z0=-92)
Aplane3 = openmc.ZPlane(z0=84)
Aplane4 = openmc.ZPlane(z0=92)
Aplane5 = openmc.ZPlane(z0=100, boundary_type=‘reflective’)
Apallet_out_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.405765)
Acladout_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.47498)
Acladin_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.41402)
Ar1 = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.5715)
Ar2 = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.6120)
Bplane1 = openmc.ZPlane(z0=-100, boundary_type=‘reflective’)
Bplane2 = openmc.ZPlane(z0=-92)
Bplane3 = openmc.ZPlane(z0=84)
Bplane4 = openmc.ZPlane(z0=92)
Bplane5 = openmc.ZPlane(z0=100, boundary_type=‘reflective’)
Bpallet_out_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.405765)
Bcladout_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.47498)
Bcladin_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.41402)
Br1 = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.5715)
Br2 = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.6120)
Cplane1 = openmc.ZPlane(z0=-100, boundary_type=‘reflective’)
Cplane2 = openmc.ZPlane(z0=-92)
Cplane3 = openmc.ZPlane(z0=84)
Cplane4 = openmc.ZPlane(z0=92)
Cplane5 = openmc.ZPlane(z0=100, boundary_type=‘reflective’)
Cpallet_out_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.405765)
Ccladout_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.47498)
Ccladin_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.41402)
Cr1 = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.5715)
Cr2 = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.6120)
Dplane1 = openmc.ZPlane(z0=-100, boundary_type=‘reflective’)
Dplane2 = openmc.ZPlane(z0=-92)
Dplane3 = openmc.ZPlane(z0=84)
Dplane4 = openmc.ZPlane(z0=92)
Dplane5 = openmc.ZPlane(z0=100, boundary_type=‘reflective’)
Dpallet_out_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.405765)
Dcladout_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.47498)
Dcladin_rad = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.41402)
Dr1 = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.5715)
Dr2 = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.6120)
AcellA1B1C1D1 = openmc.Cell(fill=A14B14C14D13, region=-Apallet_out_rad & +Aplane4 & -Aplane5)
AcellA1B1C1D1m = openmc.Cell(fill=A14B14C14D13, region=-Apallet_out_rad & +Aplane4 & -Aplane5)
BcellA1B1C1D1 = openmc.Cell(fill=A14B14C14D13, region=-Bpallet_out_rad & +Bplane4 & -Bplane5)
BcellA1B1C1D1m = openmc.Cell(fill=A14B14C14D13, region=-Bpallet_out_rad & +Bplane4 & -Bplane5)
CcellA1B1C1D1 = openmc.Cell(fill=A14B14C14D13, region=-Cpallet_out_rad & +Cplane4 & -Cplane5)
CcellA1B1C1D1m = openmc.Cell(fill=A14B14C14D13, region=-Cpallet_out_rad & +Cplane4 & -Cplane5)
DcellA1B1C1D1 = openmc.Cell(fill=A14B14C14D13, region=-Dpallet_out_rad & +Dplane4 & -Dplane5)
AcellA2B2C2m = openmc.Cell(fill=A233B233C233, region=-Apallet_out_rad & +Aplane3 & -Aplane4)
BcellA2B2C2m = openmc.Cell(fill=A233B233C233, region=-Bpallet_out_rad & +Bplane3 & -Bplane4)
CcellA2B2C2m = openmc.Cell(fill=A233B233C233, region=-Cpallet_out_rad & +Cplane3 & -Cplane4)
AcellA3m = openmc.Cell(fill=A3_UO2_Gd2O3, region=-Apallet_out_rad & +Aplane2 & -Aplane3)
BcellB3m = openmc.Cell(fill=B3_UO2_Gd2O3, region=-Bpallet_out_rad & +Bplane2 & -Bplane3)
CcellC3m = openmc.Cell(fill=C3_UO2_Gd2O3, region=-Cpallet_out_rad & +Cplane2 & -Cplane3)
AcellA4B4C4D4 = openmc.Cell(fill=A14B14C14D13, region=-Apallet_out_rad & +Aplane1 & -Aplane2)
BcellA4B4C4D4 = openmc.Cell(fill=A14B14C14D13, region=-Bpallet_out_rad & +Bplane1 & -Bplane2)
CcellA4B4C4D4 = openmc.Cell(fill=A14B14C14D13, region=-Cpallet_out_rad & +Cplane1 & -Cplane2)
AcellA4B4C4D4m = openmc.Cell(fill=A14B14C14D13, region=-Apallet_out_rad & +Aplane1 & -Aplane2)
BcellA4B4C4D4m = openmc.Cell(fill=A14B14C14D13, region=-Bpallet_out_rad & +Bplane1 & -Bplane2)
CcellA4B4C4D4m = openmc.Cell(fill=A14B14C14D13, region=-Cpallet_out_rad & +Cplane1 & -Cplane2)
DcellA4B4C4D4 = openmc.Cell(fill=A14B14C14D13, region=-Dpallet_out_rad & +Dplane1 & -Dplane2)
AcellA23B23C23 = openmc.Cell(fill=A14B14C14D13, region=-Apallet_out_rad & +Aplane2 & -Aplane4)
BcellA23B23C23 = openmc.Cell(fill=A14B14C14D13, region=-Bpallet_out_rad & +Bplane2 & -Bplane4)
CcellA23B23C23 = openmc.Cell(fill=A14B14C14D13, region=-Cpallet_out_rad & +Cplane2 & -Cplane4)
DcellD23 = openmc.Cell(fill=D2, region=-Dpallet_out_rad & +Dplane2 & -Dplane4)
Agap_cell = openmc.Cell(fill=helium, region=+Apallet_out_rad & -Acladin_rad & +Aplane1 & -Aplane5)
Agap_cellm = openmc.Cell(fill=helium, region=+Apallet_out_rad & -Acladin_rad & +Aplane1 & -Aplane5)
Bgap_cell = openmc.Cell(fill=helium, region=+Bpallet_out_rad & -Bcladin_rad & +Bplane1 & -Bplane5)
Bgap_cellm = openmc.Cell(fill=helium, region=+Bpallet_out_rad & -Bcladin_rad & +Bplane1 & -Bplane5)
Cgap_cell = openmc.Cell(fill=helium, region=+Cpallet_out_rad & -Ccladin_rad & +Cplane1 & -Cplane5)
Cgap_cellm = openmc.Cell(fill=helium, region=+Cpallet_out_rad & -Ccladin_rad & +Cplane1 & -Cplane5)
Dgap_cell = openmc.Cell(fill=helium, region=+Dpallet_out_rad & -Dcladin_rad & +Dplane1 & -Dplane5)
Aclad_cell = openmc.Cell(fill=m5, region=+Acladin_rad & -Acladout_rad & +Aplane1 & -Aplane5)
Aclad_cellm = openmc.Cell(fill=m5, region=+Acladin_rad & -Acladout_rad & +Aplane1 & -Aplane5)
Bclad_cell = openmc.Cell(fill=m5, region=+Bcladin_rad & -Bcladout_rad & +Bplane1 & -Bplane5)
Bclad_cellm = openmc.Cell(fill=m5, region=+Bcladin_rad & -Bcladout_rad & +Bplane1 & -Bplane5)
Cclad_cell = openmc.Cell(fill=m5, region=+Ccladin_rad & -Ccladout_rad & +Cplane1 & -Cplane5)
Cclad_cellm = openmc.Cell(fill=m5, region=+Ccladin_rad & -Ccladout_rad & +Cplane1 & -Cplane5)
Dclad_cell = openmc.Cell(fill=m5, region=+Dcladin_rad & -Dcladout_rad & +Dplane1 & -Dplane5)
Agi1 = openmc.Cell(fill=water, region=-Ar1)
Bgi1 = openmc.Cell(fill=water, region=-Br1)
Cgi1 = openmc.Cell(fill=water, region=-Cr1)
Dgi1 = openmc.Cell(fill=water, region=-Dr1)
Agi2 = openmc.Cell(fill=m5, region=-Ar2 & +Ar1)
Bgi2 = openmc.Cell(fill=m5, region=-Br2 & +Br1)
Cgi2 = openmc.Cell(fill=m5, region=-Cr2 & +Cr1)
Dgi2 = openmc.Cell(fill=m5, region=-Dr2 & +Dr1)
Awatercellp = openmc.Cell(fill=water, region=+Acladout_rad)
Bwatercellp = openmc.Cell(fill=water, region=+Bcladout_rad)
Cwatercellp = openmc.Cell(fill=water, region=+Ccladout_rad)
Dwatercellp = openmc.Cell(fill=water, region=+Dcladout_rad)
Awatercellpm = openmc.Cell(fill=water, region=+Acladout_rad)
Bwatercellpm = openmc.Cell(fill=water, region=+Bcladout_rad)
Cwatercellpm = openmc.Cell(fill=water, region=+Ccladout_rad)
Dwatercellpm = openmc.Cell(fill=water, region=+Dcladout_rad)
Awatercellg = openmc.Cell(fill=water, region=+Ar2)
Bwatercellg = openmc.Cell(fill=water, region=+Br2)
Cwatercellg = openmc.Cell(fill=water, region=+Cr2)
Dwatercellg = openmc.Cell(fill=water, region=+Dr2)
rodAuni1 = openmc.Universe(cells=[AcellA1B1C1D1, AcellA23B23C23, AcellA4B4C4D4, Agap_cell, Aclad_cell, Awatercellp])
rodAuni2mix = openmc.Universe(cells=[AcellA1B1C1D1m, AcellA2B2C2m, AcellA3m, AcellA4B4C4D4m, Agap_cellm, Aclad_cellm, Awatercellpm])
rodBuni1 = openmc.Universe(cells=[BcellA1B1C1D1, BcellA23B23C23, BcellA4B4C4D4, Bgap_cell, Bclad_cell, Bwatercellp])
rodBuni2mix = openmc.Universe(cells=[BcellA1B1C1D1m, BcellA2B2C2m, BcellB3m, BcellA4B4C4D4m, Bgap_cellm, Bclad_cellm, Bwatercellpm])
rodCuni1 = openmc.Universe(cells=[CcellA1B1C1D1, CcellA23B23C23, CcellA4B4C4D4, Cgap_cell, Cclad_cell, Cwatercellp])
rodCuni2mix = openmc.Universe(cells=[CcellA1B1C1D1m, CcellA2B2C2m, CcellC3m, CcellA4B4C4D4m, Cgap_cellm, Cclad_cellm, Cwatercellpm])
rodDuni = openmc.Universe(cells=[DcellA1B1C1D1, DcellD23, DcellA4B4C4D4, Dgap_cell, Dclad_cell, Dwatercellp])
Agiuni = openmc.Universe(cells=[Agi1, Agi2, Awatercellg])
Bgiuni = openmc.Universe(cells=[Bgi1, Bgi2, Bwatercellg])
Cgiuni = openmc.Universe(cells=[Cgi1, Cgi2, Cwatercellg])
Dgiuni = openmc.Universe(cells=[Dgi1, Dgi2, Dwatercellg])
#rodAuni2mix.plot(color_by=‘material’, basis=‘xz’, width=(2, 200))
Aouter = openmc.Cell(fill=water)
Aouter_uni = openmc.Universe(cells=[Aouter])
Bouter = openmc.Cell(fill=water)
Bouter_uni = openmc.Universe(cells=[Bouter])
Couter = openmc.Cell(fill=water)
Couter_uni = openmc.Universe(cells=[Couter])
Douter = openmc.Cell(fill=water)
Douter_uni = openmc.Universe(cells=[Douter])
#A_assembly.pitch = (1.259,)
A_assembly = openmc.RectLattice()
A_assembly.center = (0, 0)
A_assembly.pitch = (1.259, 1.259)
A_assembly.lower_left = [-1.259 * 17. / 2.0] * 2
A_assembly.outer = Aouter_uni
A_assembly.universes = [
[rodAuni1]*17,
[rodAuni1]*17,
[rodAuni1]*4 + [rodAuni2mix] + [Agiuni] + [rodAuni1]*2 + [Agiuni] + [rodAuni1]*2 + [Agiuni] + [rodAuni2mix] + [rodAuni1]*4,
[rodAuni1]*3 + [Agiuni] + [rodAuni1]*2 + [rodAuni2mix] + [rodAuni1] + [rodAuni2mix] + [rodAuni1] + [rodAuni2mix] + [rodAuni1]*2 + [Agiuni] + [rodAuni1]*3,
[rodAuni1]*2 + [rodAuni2mix] + [rodAuni1] + [rodAuni2mix] + [rodAuni1]*7 + [rodAuni2mix] + [rodAuni1] + [rodAuni2mix] + [rodAuni1]*2,
[rodAuni1] + [rodAuni1] + [Agiuni] + [rodAuni1] + [rodAuni1] + [Agiuni] + [rodAuni1] + [rodAuni2mix] + [Agiuni] + [rodAuni2mix] + [rodAuni1] + [Agiuni] + [rodAuni1] + [rodAuni1] + [Agiuni] + [rodAuni1] + [rodAuni1],
[rodAuni1]*3 + [rodAuni2mix] + [rodAuni1]*9 + [rodAuni2mix] + [rodAuni1]*3,
[rodAuni1]*5 + [rodAuni2mix] + [rodAuni1]*5 + [rodAuni2mix] + [rodAuni1]*5,
[rodAuni1]*2 + [Agiuni] + [rodAuni2mix] + [rodAuni1] + [Agiuni] + [rodAuni1]*2 + [Agiuni] + [rodAuni1]*2 + [Agiuni] + [rodAuni1] + [rodAuni2mix] + [Agiuni] + [rodAuni1]*2,
[rodAuni1]*5 + [rodAuni2mix] + [rodAuni1]*5 + [rodAuni2mix] + [rodAuni1]*5,
[rodAuni1]*3 + [rodAuni2mix] + [rodAuni1]*9 + [rodAuni2mix] + [rodAuni1]*3,
[rodAuni1] + [rodAuni1] + [Agiuni] + [rodAuni1] + [rodAuni1] + [Agiuni] + [rodAuni1] + [rodAuni2mix] + [Agiuni] + [rodAuni2mix] + [rodAuni1] + [Agiuni] + [rodAuni1] + [rodAuni1] + [Agiuni] + [rodAuni1] + [rodAuni1],
[rodAuni1]*2 + [rodAuni2mix] + [rodAuni1] + [rodAuni2mix] + [rodAuni1]*7 + [rodAuni2mix] + [rodAuni1] + [rodAuni2mix] + [rodAuni1]*2,
[rodAuni1]*3 + [Agiuni] + [rodAuni1]*2 + [rodAuni2mix] + [rodAuni1] + [rodAuni2mix] + [rodAuni1] + [rodAuni2mix] + [rodAuni1]*2 + [Agiuni] + [rodAuni1]*3,
[rodAuni1]*4 + [rodAuni2mix] + [Agiuni] + [rodAuni1]*2 + [Agiuni] + [rodAuni1]*2 + [Agiuni] + [rodAuni2mix] + [rodAuni1]*4,
[rodAuni1]*17,
[rodAuni1]*17,
]
B_assembly = openmc.RectLattice()
B_assembly.center = (0, 0)
B_assembly.pitch = (1.259, 1.259)
B_assembly.lower_left = [-1.259 * 17. / 2.0] * 2
B_assembly.outer = Bouter_uni
B_assembly.universes = [
[rodBuni1]*17,
[rodBuni1]*17,
[rodBuni1]*4 + [rodBuni2mix] + [Bgiuni] + [rodBuni1]*2 + [Bgiuni] + [rodBuni1]*2 + [Bgiuni] + [rodBuni2mix] + [rodBuni1]*4,
[rodBuni1]*3 + [Bgiuni] + [rodBuni1]*2 + [rodBuni2mix] + [rodBuni1] + [rodBuni2mix] + [rodBuni1] + [rodBuni2mix] + [rodBuni1]*2 + [Bgiuni] + [rodBuni1]*3,
[rodBuni1]*2 + [rodBuni2mix] + [rodBuni1] + [rodBuni2mix] + [rodBuni1]*7 + [rodBuni2mix] + [rodBuni1] + [rodBuni2mix] + [rodBuni1]*2,
[rodBuni1] + [rodBuni1] + [Bgiuni] + [rodBuni1] + [rodBuni1] + [Bgiuni] + [rodBuni1] + [rodBuni2mix] + [Bgiuni] + [rodBuni2mix] + [rodBuni1] + [Bgiuni] + [rodBuni1] + [rodBuni1] + [Bgiuni] + [rodBuni1] + [rodBuni1],
[rodBuni1]*3 + [rodBuni2mix] + [rodBuni1]*9 + [rodBuni2mix] + [rodBuni1]*3,
[rodBuni1]*5 + [rodBuni2mix] + [rodBuni1]*5 + [rodBuni2mix] + [rodBuni1]*5,
[rodBuni1]*2 + [Bgiuni] + [rodBuni2mix] + [rodBuni1] + [Bgiuni] + [rodBuni1]*2 + [Bgiuni] + [rodBuni1]*2 + [Bgiuni] + [rodBuni1] + [rodBuni2mix] + [Bgiuni] + [rodBuni1]*2,
[rodBuni1]*5 + [rodBuni2mix] + [rodBuni1]*5 + [rodBuni2mix] + [rodBuni1]*5,
[rodBuni1]*3 + [rodBuni2mix] + [rodBuni1]*9 + [rodBuni2mix] + [rodBuni1]*3,
[rodBuni1] + [rodBuni1] + [Bgiuni] + [rodBuni1] + [rodBuni1] + [Bgiuni] + [rodBuni1] + [rodBuni2mix] + [Bgiuni] + [rodBuni2mix] + [rodBuni1] + [Bgiuni] + [rodBuni1] + [rodBuni1] + [Bgiuni] + [rodBuni1] + [rodBuni1],
[rodBuni1]*2 + [rodBuni2mix] + [rodBuni1] + [rodBuni2mix] + [rodBuni1]*7 + [rodBuni2mix] + [rodBuni1] + [rodBuni2mix] + [rodBuni1]*2,
[rodBuni1]*3 + [Bgiuni] + [rodBuni1]*2 + [rodBuni2mix] + [rodBuni1] + [rodBuni2mix] + [rodBuni1] + [rodBuni2mix] + [rodBuni1]*2 + [Bgiuni] + [rodBuni1]*3,
[rodBuni1]*4 + [rodBuni2mix] + [Bgiuni] + [rodBuni1]*2 + [Bgiuni] + [rodBuni1]*2 + [Bgiuni] + [rodBuni2mix] + [rodBuni1]*4,
[rodBuni1]*17,
[rodBuni1]*17,
]
C_assembly = openmc.RectLattice()
C_assembly.center = (0, 0)
C_assembly.pitch = (1.259, 1.259)
C_assembly.lower_left = [-1.259 * 17. / 2.0] * 2
C_assembly.outer = Couter_uni
C_assembly.universes = [
[rodCuni1]*17,
[rodCuni1]*17,
[rodCuni1]*4 + [rodCuni2mix] + [Cgiuni] + [rodCuni1]*2 + [Cgiuni] + [rodCuni1]*2 + [Cgiuni] + [rodCuni2mix] + [rodCuni1]*4,
[rodCuni1]*3 + [Cgiuni] + [rodCuni1]*2 + [rodCuni2mix] + [rodCuni1] + [rodCuni2mix] + [rodCuni1] + [rodCuni2mix] + [rodCuni1]*2 + [Cgiuni] + [rodCuni1]*3,
[rodCuni1]*2 + [rodCuni2mix] + [rodCuni1] + [rodCuni2mix] + [rodCuni1]*7 + [rodCuni2mix] + [rodCuni1] + [rodCuni2mix] + [rodCuni1]*2,
[rodCuni1] + [rodCuni1] + [Cgiuni] + [rodCuni1] + [rodCuni1] + [Cgiuni] + [rodCuni1] + [rodCuni2mix] + [Cgiuni] + [rodCuni2mix] + [rodCuni1] + [Cgiuni] + [rodCuni1] + [rodCuni1] + [Cgiuni] + [rodCuni1] + [rodCuni1],
[rodCuni1]*3 + [rodCuni2mix] + [rodCuni1]*9 + [rodCuni2mix] + [rodCuni1]*3,
[rodCuni1]*5 + [rodCuni2mix] + [rodCuni1]*5 + [rodCuni2mix] + [rodCuni1]*5,
[rodCuni1]*2 + [Cgiuni] + [rodCuni2mix] + [rodCuni1] + [Cgiuni] + [rodCuni1]*2 + [Cgiuni] + [rodCuni1]*2 + [Cgiuni] + [rodCuni1] + [rodCuni2mix] + [Cgiuni] + [rodCuni1]*2,
[rodCuni1]*5 + [rodCuni2mix] + [rodCuni1]*5 + [rodCuni2mix] + [rodCuni1]*5,
[rodCuni1]*3 + [rodCuni2mix] + [rodCuni1]*9 + [rodCuni2mix] + [rodCuni1]*3,
[rodCuni1] + [rodCuni1] + [Cgiuni] + [rodCuni1] + [rodCuni1] + [Cgiuni] + [rodCuni1] + [rodCuni2mix] + [Cgiuni] + [rodCuni2mix] + [rodCuni1] + [Cgiuni] + [rodCuni1] + [rodCuni1] + [Cgiuni] + [rodCuni1] + [rodCuni1],
[rodCuni1]*2 + [rodCuni2mix] + [rodCuni1] + [rodCuni2mix] + [rodCuni1]*7 + [rodCuni2mix] + [rodCuni1] + [rodCuni2mix] + [rodCuni1]*2,
[rodCuni1]*3 + [Cgiuni] + [rodCuni1]*2 + [rodCuni2mix] + [rodCuni1] + [rodCuni2mix] + [rodCuni1] + [rodCuni2mix] + [rodCuni1]*2 + [Cgiuni] + [rodCuni1]*3,
[rodCuni1]*4 + [rodCuni2mix] + [Cgiuni] + [rodCuni1]*2 + [Cgiuni] + [rodCuni1]*2 + [Cgiuni] + [rodCuni2mix] + [rodCuni1]*4,
[rodCuni1]*17,
[rodCuni1]*17,
]
D_assembly = openmc.RectLattice()
D_assembly.center = (0, 0)
D_assembly.pitch = (1.259, 1.259)
D_assembly.lower_left = [-1.259 * 17. / 2.0] * 2
D_assembly.outer = Douter_uni
D_assembly.universes = [
[rodDuni]*17,
[rodDuni]*17,
[rodDuni]*4 + [rodDuni] + [Dgiuni] + [rodDuni]*2 + [Dgiuni] + [rodDuni]*2 + [Dgiuni] + [rodDuni] + [rodDuni]*4,
[rodDuni]*3 + [Dgiuni] + [rodDuni]*2 + [rodDuni] + [rodDuni] + [rodDuni] + [rodDuni] + [rodDuni] + [rodDuni]*2 + [Dgiuni] + [rodDuni]*3,
[rodDuni]*2 + [rodDuni] + [rodDuni] + [rodDuni] + [rodDuni]*7 + [rodDuni] + [rodDuni] + [rodDuni] + [rodDuni]*2,
[rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni] + [rodDuni],
[rodDuni]*3 + [rodDuni] + [rodDuni]*9 + [rodDuni] + [rodDuni]*3,
[rodDuni]*5 + [rodDuni] + [rodDuni]*5 + [rodDuni] + [rodDuni]*5,
[rodDuni]*2 + [Dgiuni] + [rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni]*2 + [Dgiuni] + [rodDuni]*2 + [Dgiuni] + [rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni]*2,
[rodDuni]*5 + [rodDuni] + [rodDuni]*5 + [rodDuni] + [rodDuni]*5,
[rodDuni]*3 + [rodDuni] + [rodDuni]*9 + [rodDuni] + [rodDuni]*3,
[rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni] + [rodDuni] + [Dgiuni] + [rodDuni] + [rodDuni],
[rodDuni]*2 + [rodDuni] + [rodDuni] + [rodDuni] + [rodDuni]*7 + [rodDuni] + [rodDuni] + [rodDuni] + [rodDuni]*2,
[rodDuni]*3 + [Dgiuni] + [rodDuni]*2 + [rodDuni] + [rodDuni] + [rodDuni] + [rodDuni] + [rodDuni] + [rodDuni]*2 + [Dgiuni] + [rodDuni]*3,
[rodDuni]*4 + [rodDuni] + [Dgiuni] + [rodDuni]*2 + [Dgiuni] + [rodDuni]*2 + [Dgiuni] + [rodDuni] + [rodDuni]*4,
[rodDuni]*17,
[rodDuni]*17,
]
box = -openmc.model.RectangularPrism(height=21.403, width=21.403, boundary_type=‘reflective’) & -openmc.ZPlane(z0=100, boundary_type=‘reflective’) & +openmc.ZPlane(z0=-100, boundary_type=‘reflective’)
#cell1 = openmc.Cell(fill=A_assembly, region=box)
#cell2 = openmc.Cell(fill=B_assembly, region=box)
#cell3 = openmc.Cell(fill=C_assembly, region=box)
cell4 = openmc.Cell(fill=D_assembly, region=box)
#gemA = openmc.Geometry([cell1])
#gemB = openmc.Geometry([cell2])
#gemC = openmc.Geometry([cell3])
gemD = openmc.Geometry([cell4])
#gemA.export_to_xml()
#gemB.export_to_xml()
#gemC.export_to_xml()
gemD.export_to_xml()
src = openmc.Source()
src.space = openmc.stats.Point((0.0, 0.0, 0.0))
src.angle = openmc.stats.Isotropic()
src.energy = openmc.stats.Discrete([2.0e6], [1.0])
settings = openmc.Settings()
#settings.run_mode = ‘fixed source’
settings.batches = 300
settings.inactive = 50
settings.particles = 100000
settings.source = src
settings.temperature = {‘method’: ‘interpolation’}
#settings.create_delayed_neutrons=False
settings.ifp_n_generation = 40
settings.export_to_xml()
tally = openmc.Tally(name=“ifp-scores”)
tally.scores = [
“ifp-time-numerator”,
“ifp-beta-numerator”,
“ifp-denominator”
]
tallies = openmc.Tallies([tally])
tallies.export_to_xml()
openmc.run()
So from the look of it, the IFP method within OpenMC might only be included within development branches. The documentation for it is version 0.15.3.dev29+ga7d1ceba3 which is not the newest main stable version of 0.15.2, and seems to be a newer dev version. I can’t find anything on IFP method within the 0.15.2 documentation. Previous development branches developed by individuals using the code were discussed in these previous user support threads here,
I don’t know if they have formally been consolidated into the main branch. I tried to play around with it on my pc running 0.15.2, and I was able to get the tallies to be exported to XML, but I don’t know how to run the IFP simulation, and the documentation doesn’t state how to actually run the simulation. I tried defining ifp_n_generation in the settings to try to have it run the sim but that didn’t seem to work. Therefore, it might be in the newest version, the tally definition is at least, but I don’t know if the full process is integrated. @paulromano can probably give a confirmation if it has been implemented.
@paulromano wanna say something?
openmc/openmc:develop can be used to use “ifp method”, i’m using this docker image and the issue is being solved.
Hi, that is correct. In theory, the IFP method will be implemented in the version 0.15.3 which is the next one. To utilize it you have to clone the develop branch of OpenMC.