Currently, the NEMDriver class does the pin power reconstruction from its own simulation results. This is a very "quick and dirty" way of doing it, and the method doesn't even take into account things like the assembly gap, and just assumes that the assembly is only made of pins. Also, this means that one can't do pin power reconstruction with the finite-difference diffusion solver ! It also doesn't support pin power reconstruction with half/quarter/eighth symmetry, which is also quite limiting.
In the long term, a new class which takes information about the pin pitch, assembly pitch, assembly shapes, core shape, and symmetry of the problem should be created which can perform the pin power reconstruction. I think the best way to accomplish this would be for it to take a shared pointer to the solver class, and then it could be used to reconstruct pin powers on subsequent solves (i.e. depletion or transients). It should be able to use either the nodal diffusion solver or the finite-difference solver to accomplish this task, and both of those solvers should likely be modified so that they have a similar interface (i.e. can evaluate the flux / power at a given point) to accomplish this.
Currently, the
NEMDriverclass does the pin power reconstruction from its own simulation results. This is a very "quick and dirty" way of doing it, and the method doesn't even take into account things like the assembly gap, and just assumes that the assembly is only made of pins. Also, this means that one can't do pin power reconstruction with the finite-difference diffusion solver ! It also doesn't support pin power reconstruction with half/quarter/eighth symmetry, which is also quite limiting.In the long term, a new class which takes information about the pin pitch, assembly pitch, assembly shapes, core shape, and symmetry of the problem should be created which can perform the pin power reconstruction. I think the best way to accomplish this would be for it to take a shared pointer to the solver class, and then it could be used to reconstruct pin powers on subsequent solves (i.e. depletion or transients). It should be able to use either the nodal diffusion solver or the finite-difference solver to accomplish this task, and both of those solvers should likely be modified so that they have a similar interface (i.e. can evaluate the flux / power at a given point) to accomplish this.