Detectors in transient simulations

Report all good and bad behavior here
Post Reply
Diego
Posts: 73
Joined: Wed Jun 01, 2011 8:49 pm

Detectors in transient simulations

Post by Diego » Tue Apr 17, 2018 12:53 pm

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 1723 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

User avatar
Jaakko Leppänen
Site Admin
Posts: 2387
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: Detectors in transient simulations

Post by Jaakko Leppänen » Tue Apr 17, 2018 2:35 pm

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?
- Jaakko

Ville Valtavirta
Posts: 468
Joined: Fri Sep 07, 2012 1:43 pm
Security question 1: No
Security question 2: 92

Re: Detectors in transient simulations

Post by Ville Valtavirta » Tue Apr 17, 2018 6:24 pm

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

Diego
Posts: 73
Joined: Wed Jun 01, 2011 8:49 pm

Re: Detectors in transient simulations

Post by Diego » Wed Apr 18, 2018 12:46 pm

Thanks! With the correction it works properly.
Best,
Diego

Post Reply