# Prompt neutron lifetime

**URL:** <https://openmc.discourse.group/t/prompt-neutron-lifetime/4853>\
**Category:** User Support\
**Created:** [October 1, 2024, 10:27am UTC](https://openmc.discourse.group/t/prompt-neutron-lifetime/4853 "2024-10-01T10:27:14Z")\
**Posts on this page:** 4\
**Page:** 1

<div class="post-metadata">

**Author:** ![fyh\_sinap](https://avatars.discourse-cdn.com/v4/letter/f/ed655f/32.png) [@fyh\_sinap](https://openmc.discourse.group/u/fyh_sinap)\
**Post date:** [October 1, 2024, 10:27am UTC](https://openmc.discourse.group/t/prompt-neutron-lifetime/4853/1 "2024-10-01T10:27:14Z")

</div>

Hello, everyone!  
I want to use point reactor kinetics equations for calculating transient problems. The delayed neutron lifetime, fraction can be obtained by mgxs class. However, I can’t find a way to process the neutron generation time or prompt neutron lifetime. Does anyone know it,very gratefully!

---

<div class="post-metadata">

**Author:** ![mkreher13](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/mkreher13/32/4093_2.png) [@mkreher13](https://openmc.discourse.group/u/mkreher13)\
**Post date:** [October 4, 2024, 11:44pm UTC](https://openmc.discourse.group/t/prompt-neutron-lifetime/4853/2 "2024-10-04T23:44:43Z")

</div>

Hello, although I’m not super fresh on the topic, I did my dissertation on transient problems in OpenMC. If you see [this paper](https://www.tandfonline.com/doi/full/10.1080/00295639.2021.1980363), equation 10 shows one way to calculate prompt neutron lifetime. Please see the PDF version of the paper, as the online version distorts equation 10.

You might also find it useful to poke around [this branch](https://github.com/mkreher13/openmc/tree/transient-coupled-two-way-extrap/openmc/kinetics) of OpenMC on my GitHub for more.

---

<div class="post-metadata">

**Author:** ![fyh\_sinap](https://avatars.discourse-cdn.com/v4/letter/f/ed655f/32.png) [@fyh\_sinap](https://openmc.discourse.group/u/fyh_sinap)\
**Post date:** [October 17, 2024, 7:18am UTC](https://openmc.discourse.group/t/prompt-neutron-lifetime/4853/3 "2024-10-17T07:18:31Z")

</div>

Thank you for answering this question. That means that the adjoint flux is required for calculation neutron prompt life. It seems that OpenMC doesn’t have that function at the moment. I will try to generate the dynamic parameters using deterministic code. Thank you very much!

---

<div class="post-metadata">

**Author:** ![mkreher13](https://yyz2.discourse-cdn.com/free1/user_avatar/openmc.discourse.group/mkreher13/32/4093_2.png) [@mkreher13](https://openmc.discourse.group/u/mkreher13)\
**Post date:** [October 22, 2024, 7:23am UTC](https://openmc.discourse.group/t/prompt-neutron-lifetime/4853/4 "2024-10-22T07:23:58Z")

</div>

Hi again,  
In the [Github link](https://github.com/mkreher13/openmc/tree/transient-coupled-two-way-extrap/openmc/kinetics) I mentioned, an example of prompt neutron lifetime being calculated is shown in state.py in a function called pnl().

That refers to an adjoint flux which is calculated using a source iteration with matrices of Monte Carlo tallied entries. The code for that eigenvalue calculation can be found in aux\_functions() in a function called compute\_eigenvalue(). To obtain an adjoint, call the function with the transpose of the matrices.
