Page 1 of 1

Detectors in transient simulations

Posted: Tue Apr 17, 2018 12:53 pm
by Diego
Dear colleagues;
I'm running a transient simulation of a Minicore, starting from nominal power and then performing some rho insertions (namely moving some CR). To get my results I'm including two detectors to get the total power and some profiles (covering all the geometry), namely:

Code: Select all

det 1 dr -8 void di dettime	
det 2 dr -8 void dx  -31.72 31.72 45 dy  -31.72 31.72 45   dz -182.88 182.88 10 di dettime
I suppose that the sum over all bins in det 2 should be the same as det 1, but I found a difference such as:
transientdet.png
transientdet.png (21.88 KiB) Viewed 1744 times
It doesn't look like a lack of convergence (each bin of det 2 reports a stdev of ~4%). Maybe some CFE params?
I also checked the same for critical calculations and in such case both detectors agree (det1 1 vs sum over all bins in det2). I''m missing something in dx dy dz when using di? Anyone experienced similar problems?
Thanks,
Diego

Re: Detectors in transient simulations

Posted: Tue Apr 17, 2018 2:35 pm
by Jaakko Leppänen
It shouldn't be related to CFE. If the mesh covers the full geometry, you should get the exact same scores in both detectors.

Could you send me the input?

Re: Detectors in transient simulations

Posted: Tue Apr 17, 2018 6:24 pm
by Ville Valtavirta
The source of the error was tracked to lines 373-376 of detectoroutput.c. The statement

Code: Select all

                  /* Get index */

                  idx0 = DetIdx(det0, eb0, ub0, cb0, mb0, lb0, rb0, yb0, xb0,
                                tb0);
should be

Code: Select all

                  /* Get index */

                  idx0 = DetIdx(det0, eb0, ub0, cb0, mb0, lb0, zb0, yb0, xb0,
                                tb0);
i.e. "rb0" should be replaced with "zb0".

The bug affected detectors that had a time-binning specified with the "di" option and either axial binning "dz" or multiple response functions "dr".

-Ville

Re: Detectors in transient simulations

Posted: Wed Apr 18, 2018 12:46 pm
by Diego
Thanks! With the correction it works properly.
Best,
Diego