https://github.com/openmsr/pincell_msr
openmc msr depletion capabilities example
https://github.com/openmsr/pincell_msr
depletion msr openmc removal-rate
Last synced: 6 months ago
JSON representation
openmc msr depletion capabilities example
- Host: GitHub
- URL: https://github.com/openmsr/pincell_msr
- Owner: openmsr
- Created: 2022-11-01T11:02:50.000Z (over 3 years ago)
- Default Branch: main
- Last Pushed: 2024-04-19T09:14:51.000Z (about 2 years ago)
- Last Synced: 2024-04-20T09:38:23.152Z (about 2 years ago)
- Topics: depletion, msr, openmc, removal-rate
- Language: Jupyter Notebook
- Homepage:
- Size: 1.79 MB
- Stars: 3
- Watchers: 1
- Forks: 3
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
Awesome Lists containing this project
README
# Pincell MSR depletion
This is a slightly modified version of the [OpenMC pincell depletion](https://github.com/openmc-dev/openmc/tree/main/examples/pincell_depletion) example, aimed to show some recent depletion capabilities developed primarily for Molten Salt Reactor (MSR) simulations.
## Functionalities:
### Transfer Rates
`TransferRates` methodology allows continuous material reprocessing (removal/feed) during a depletion or coupled transport-depletion simulation.
It has been integrated into OpenMC main code from version 0.14.0 and more info can be found on the user manual.
### Reactivity control
Reactivity control is currently an openmc-dev [open PR](https://github.com/openmc-dev/openmc/pull/2693). It allows to add reactivity controls during a transport-depletion simulation with the aim of keep `keff`equals to 1 as a function of some user defined parameters, such as the geometrical position of a control rod or some batchwise material addition or removal.
### Redox control
Redox control is another [open PR](https://github.com/openmc-dev/openmc/pull/2783) that adds the capability to maintain the redox of the fuel salt constant during a depletion simulation, by continuously adding or removing a user-defined buffer nuclides vector. It makes use of the `TransferRates` implementation to add a further term to the Bateman equation.
For example, UF4 frees 4 atoms of fluorine every time a fission event occurs, which will then bind with fission products with lower oxidation states than uranium, giving rise to an excess of free fluorine atoms and thus increasing the redox potential of the salt. On the other hand, neutron capture reactions, will likely increase the oxidation state of the transmuting nuclide, thus creating a deficiency of fluorine atoms and reducing the redox potential of the salt.
Assuming the oxidation states of the nuclides present in the salt are known, the redox can be computed multiplying the number of atoms $N_i$ of each nuclide present in the salt, by their oxidation states $Ox_i$:
```math
\begin{equation}
\sum_i N_i Ox_i
\end{equation}
```
The added term for a generic buffer nuclide can be expressed mathematically as:
```math
\frac{dN_b(t)}{dt} = \cdots + \frac{1}{Ox_b}\sum_i N_i\left( L_{ii}Ox_i - \sum_j G_{i\rightarrow j } Ox_j\right)
```
where $`Ox_b`$ is the oxidation state of the buffer nuclide, and $`Ox_i`$, $`Ox_j`$ are the oxidation states of the i-th nuclide undergoing a reaction or a decay and of the j-th product nuclide, respectively.
Therefore, the first term in the right hand side represents the losses of the i-th nuclide (i.e. the diagonal terms of the Bateman matrix), and the second term the gains of the j-th nuclides (i.e. the off-diagonal terms of the Bateman matrix).
## Examples:
### pincell_msr_transfer_rates
In this example `TransferRates` methodology is demonstrated for a slightly modified version of the existing [openmc pincell depletion example](https://github.com/openmc-dev/openmc-notebooks/blob/main/depletion.ipynb).
### pincell_msr_reactivity_control
The same pincell model will be used here to show some `ReactiviyControl` functionalities, where we set a geometrical parameter of the pincell to keep it critical during the transport-depletion run.
### pincell_msr_redox
In this example, we will run the pincell with molten salt (UF4) as fuel and add a redox control to keep the redox potential constant during the transport-depletion run by adding a nuclide buffer.
### pincell_msr_lumped
Here we show how to set up `TransferRates` to simulate the circulation of fuel salt from in-core to out-of-core regions and vice versa, as a function of the regions' volumes and flow rate.
### pincell_msr_dnf
The same settings of the **pincell_msr_lumped** cam be used to estimate the Delayed Neutron Fraction (DNF) out-of-core and in-core.