You might want to check the depletion notebook from openmc for transport-coupled depletion.
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Pincell Depletion\n",
"This notebook is intended to introduce the reader to the depletion interface contained in OpenMC. It is recommended that you are moderately familiar with building models using the OpenMC Python API. The earlier examples are excellent starting points, as this notebook will not focus heavily on model building.\n",
"\n",
"If you have a real power reactor, the fuel composition is constantly changing as fission events produce energy, remove some fissile isotopes, and produce fission products. Other reactions, like $(n, \\alpha)$ and $(n, \\gamma)$ will alter the composition as well. Furthermore, some nuclides undergo spontaneous decay with widely ranging frequencies. Depletion is the process of modeling this behavior.\n",
"\n",
"In this notebook, we will model a simple fuel pin in an infinite lattice using the Python API. We will then build and examine some of the necessary components for performing depletion analysis. Then, we will use the depletion interface in OpenMC to simulate the fuel pin producing power over several months. Lastly, we will wrap up with some helpful tips to improve the fidelity of depletion simulations."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
show original
The details on python api for depletion calculation could be seen in Depletion - Python API .Transport Independent Depletion section has described this feature. So you can try it before setting up the integrator and do the integrator.integrate().
Also, you might want to ask this kind of question on the User support section of the forum, which will catch more attention and show the related topics, such as this one Negative atoms when using openmc.deplete.IndependentOperator from Shimwell and others.
1 Like