Skip to content

Add support for OpenMC MGXS delayed neutron data.#1027

Open
wdhawkins wants to merge 1 commit intoOpen-Sn:mainfrom
wdhawkins:openmc_delayed_neutron_xs
Open

Add support for OpenMC MGXS delayed neutron data.#1027
wdhawkins wants to merge 1 commit intoOpen-Sn:mainfrom
wdhawkins:openmc_delayed_neutron_xs

Conversation

@wdhawkins
Copy link
Copy Markdown
Collaborator

PR Checklist

  • I have updated the user guide and/or Python API documentation if necessary.

@wdhawkins wdhawkins self-assigned this Apr 13, 2026
@wdhawkins wdhawkins added the enhancement New feature or request label Apr 13, 2026
@wdhawkins wdhawkins force-pushed the openmc_delayed_neutron_xs branch from ae85618 to 9c057c5 Compare April 13, 2026 15:56
@wdhawkins wdhawkins changed the title Add support for OpenMC MGXS dealyed neutron data. Add support for OpenMC MGXS delayed neutron data. Apr 13, 2026
@wdhawkins wdhawkins force-pushed the openmc_delayed_neutron_xs branch from 9c057c5 to c2141d7 Compare April 13, 2026 16:03
@wdhawkins wdhawkins force-pushed the openmc_delayed_neutron_xs branch from c2141d7 to 5751176 Compare April 14, 2026 18:43
Comment on lines +313 to +317
When delayed-neutron data is present, OpenSn also imports precursor-group
information, prompt and delayed fission production, precursor decay
constants, and delayed emission spectra. If ``chi-prompt`` is not present,
OpenSn falls back to the steady fission spectrum ``chi`` as the prompt
spectrum.
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

how do we ensure that we have a null transient when going from steady state to time-dependent? i think my comment has a lot of ramification that should be best discussed in person

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yes, we should probably discuss this in person. But, if I understand your question correctly, using initial_state=existing with the transient solver should prevent any artificial transient jump from the mode transition itself. So, the workflow would be:

  1. Solve the steady-state problem to convergence.
  2. Switch to time-dependent mode.
  3. Initialize TransientSolver with initial_state=existing.

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

for a null transient test to pass, nu prompt + nu delayed must be = to nu tot. furthermore, chi is not easily deduced from chi prompt and chi delayed (there is an approximation that happens to yield chi used in steady-state, so consistency is also an issue). finally, when coming from MC, there is always statistical noise in the values. So, all in all, it is not obvious to how a null transient will happen. but we definitely need to add some consistency checks between steady-state data and transient data and report back to the user what those are, so that they are informed.

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I see. We should talk more about this.

@wdhawkins wdhawkins force-pushed the openmc_delayed_neutron_xs branch from 5751176 to 232c16b Compare April 15, 2026 04:08
@wdhawkins wdhawkins force-pushed the openmc_delayed_neutron_xs branch from 232c16b to 8ca4bea Compare April 15, 2026 14:53
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

enhancement New feature or request

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants