Re: depletion problem

51 views
Skip to first unread message
Message has been deleted

Ariful Islam Pranto

unread,
Jul 18, 2020, 11:06:12 AM7/18/20
to OpenMC Users Group
Hi Shurabil,

1. For a 2D problem, the unit of power can be given as [w/cm] and volume assigned in depletion material as [cm2] because of infinite axial length. For a 3D problem, the unit of power will be [w] and volume will be assigned in depletion material as [cm3].
  Hints: fuel.volume = math.pi * radii[0] ** 2*height

2. You have to provide either [Power] or [Power Density] in the depletion algorithm class. Power or Power density
 
3. Prev_res. If you are using openmc-0.11, you will see the simulation start from scratch which is a bug and solved in develop branch.

On Saturday, July 18, 2020 at 8:03:48 PM UTC+6 Shurabil wrote:
Hello guys,
1)In the pincell depletion problem,a pincell with infinite axial length has been modeled but before exporting the file to materials.xml the volume is assigned as if the pincell were only 1 cm in height.Why is that?Why didnt the person who wrote this just modeled the pincell as 1 cm from the beginning.Please explain this
2)In the lwr next gen benchmark problem,an infinite(axial length) pincell/assembly and power density(w/cm3,w/cm,w/gHM) is given.I dont have power.There is no volume.How do i normalize the power density to power for a pincell/assembly of infinite axial length
3)How do I use the prev_result parametre of openmc depletion operator to start calculation from where i left.For example,if i calculate depletion for 10 days and then i want to start from 10 days and go to 20 days how do i do that?I tried to use prev_result and loaded the depletion results.h5 for 10 days and insert it in openmc deplete operator but for some reasons the calculation starts from zero second again.What am I doing wrong?
Thanks in advance
Message has been deleted

Ariful Islam Pranto

unread,
Jul 18, 2020, 11:50:34 AM7/18/20
to OpenMC Users Group
Details of your  2nd question,

power  = x [w/cm3] * volume
power = x[w/cm] * Area >> Your first question Why 1 cm?
Power = x[w/gHM] * operator.heavy_metal[g]

On Saturday, July 18, 2020 at 9:34:58 PM UTC+6 Shurabil wrote:
1)that wasnt my question
2)i i already knew that
3)thanks
Reply all
Reply to author
Forward
Message has been deleted
0 new messages