search

cyclus / cycamore

The CYClus Additional Module REpository
Other
15 stars 40 forks source link

Separations archetype: metastable nuclides not going to leftover stream #575

Open danielwojtaszek opened 9 months ago

danielwojtaszek commented 9 months ago

I've noticed that some metastable nuclides (such as 471080001 and 330820001) in the feed inventory of the separations archetype that are not placed in any product streams end up not going to the leftover stream, and thus disappear from the system altogether. Note that the decay control parameter of the simulation is set to never.

The disappearance of 471080001 (aka Ag-108m) is particularly odd since its half-life is > 400 years.

danielwojtaszek commented 9 months ago

I can work around this by creating a product stream for these metastable nuclides, but I found that some of these nuclides disappear when transferring them out of the separations facility to a sink facility.

gonuke commented 9 months ago

Thanks for reporting this @danielwojtaszek . Can you share the input file snippet that describes how you define you separations facility and its separation efficiencies?

danielwojtaszek commented 9 months ago

Hi Paul (@gonuke). Here's the prototype definition I was using when I discovered the issue.

<facility>
    <name>Reprocessing</name>
    <config>
      <Separations>
        <feedbuf_size>1e299</feedbuf_size>
        <feed_commods><val>Cooled_U_UNF_comm</val></feed_commods>
        <leftover_commod>HLW_comm</leftover_commod>
        <streams>
          <item><commod>TRU_comm</commod><info><buf_size>-1</buf_size>
            <efficiencies>
              <item><comp>942380000</comp><eff>0.990000</eff></item>
              <item><comp>942390000</comp><eff>0.990000</eff></item>
              <item><comp>942400000</comp><eff>0.990000</eff></item>
              <item><comp>942410000</comp><eff>0.990000</eff></item>
              <item><comp>942420000</comp><eff>0.990000</eff></item>
              <item><comp>942440000</comp><eff>0.990000</eff></item>
              <item><comp>932350000</comp><eff>0.990000</eff></item>
              <item><comp>932360000</comp><eff>0.990000</eff></item>
              <item><comp>932370000</comp><eff>0.990000</eff></item>
              <item><comp>932380000</comp><eff>0.990000</eff></item>
              <item><comp>932390000</comp><eff>0.990000</eff></item>
              <item><comp>952410000</comp><eff>0.990000</eff></item>
              <item><comp>952420000</comp><eff>0.990000</eff></item>
              <item><comp>952420001</comp><eff>0.990000</eff></item>
              <item><comp>952430000</comp><eff>0.990000</eff></item>
              <item><comp>952440000</comp><eff>0.990000</eff></item>
              <item><comp>952440001</comp><eff>0.990000</eff></item>
              <item><comp>962400000</comp><eff>0.990000</eff></item>
              <item><comp>962410000</comp><eff>0.990000</eff></item>
              <item><comp>962420000</comp><eff>0.990000</eff></item>
              <item><comp>962430000</comp><eff>0.990000</eff></item>
              <item><comp>962440000</comp><eff>0.990000</eff></item>
              <item><comp>962450000</comp><eff>0.990000</eff></item>
              <item><comp>962460000</comp><eff>0.990000</eff></item>
              <item><comp>962470000</comp><eff>0.990000</eff></item>
              <item><comp>962480000</comp><eff>0.990000</eff></item>
              <item><comp>962490000</comp><eff>0.990000</eff></item>
              <item><comp>962500000</comp><eff>0.990000</eff></item>
              <item><comp>952380000</comp><eff>0.990000</eff></item>
              <item><comp>952390000</comp><eff>0.990000</eff></item>
              <item><comp>952400000</comp><eff>0.990000</eff></item>
              <item><comp>952450000</comp><eff>0.990000</eff></item>
              <item><comp>952460000</comp><eff>0.990000</eff></item>
              <item><comp>962380000</comp><eff>0.990000</eff></item>
              <item><comp>962390000</comp><eff>0.990000</eff></item>
              <item><comp>962510000</comp><eff>0.990000</eff></item>
            </efficiencies>
          </info></item>
          <item><commod>RepU_comm</commod><info><buf_size>-1</buf_size>
            <efficiencies>
              <item><comp>922320000</comp><eff>1.000000</eff></item>
              <item><comp>922330000</comp><eff>1.000000</eff></item>
              <item><comp>922340000</comp><eff>1.000000</eff></item>
              <item><comp>922350000</comp><eff>1.000000</eff></item>
              <item><comp>922360000</comp><eff>1.000000</eff></item>
              <item><comp>922380000</comp><eff>1.000000</eff></item>
            </efficiencies>
          </info></item>
          <item><commod>REFP_comm</commod><info><buf_size>-1</buf_size>
            <efficiencies>
              <item><comp>571380000</comp><eff>0.990000</eff></item>
              <item><comp>571390000</comp><eff>0.990000</eff></item>
              <item><comp>581360000</comp><eff>0.990000</eff></item>
              <item><comp>581380000</comp><eff>0.990000</eff></item>
              <item><comp>581400000</comp><eff>0.990000</eff></item>
              <item><comp>581420000</comp><eff>0.990000</eff></item>
              <item><comp>591410000</comp><eff>0.990000</eff></item>
              <item><comp>601420000</comp><eff>0.990000</eff></item>
              <item><comp>601430000</comp><eff>0.990000</eff></item>
              <item><comp>601440000</comp><eff>0.990000</eff></item>
              <item><comp>601450000</comp><eff>0.990000</eff></item>
              <item><comp>601460000</comp><eff>0.990000</eff></item>
              <item><comp>601480000</comp><eff>0.990000</eff></item>
              <item><comp>601500000</comp><eff>0.990000</eff></item>
              <item><comp>641520000</comp><eff>0.990000</eff></item>
              <item><comp>641540000</comp><eff>0.990000</eff></item>
              <item><comp>641550000</comp><eff>0.990000</eff></item>
              <item><comp>641560000</comp><eff>0.990000</eff></item>
              <item><comp>641570000</comp><eff>0.990000</eff></item>
              <item><comp>641580000</comp><eff>0.990000</eff></item>
              <item><comp>641600000</comp><eff>0.990000</eff></item>
            </efficiencies>
          </info></item>
        </streams>
      </Separations>
    </config>
  </facility>