Power Distribution for each burnup step using SIE method

Questions and discussion about applications, input, output and general user topics
Post Reply
Yousef1986
Posts: 73
Joined: Fri Sep 04, 2015 3:33 pm
Security question 1: No
Security question 2: 92
Location: Karlsruhe, Germany
Contact:

Power Distribution for each burnup step using SIE method

Post by Yousef1986 » Wed Jun 21, 2017 2:48 pm

Dear Jaakko,

I am trying to print power distribution using "set cpd" card at each burnup step. I am using the SIE method for the depletion and got the following problem:
The power distribution is printed only at fresh calculation (i.e. zero burnups and at the predictor step only) and as you know there is only a single predictor step in the whole burnup calculation when SIE method is used.

So, I made a change in printcoredistr.c after line 41 to allow for printing the power for each corrector step/iteration:

/* Check corrector step */

if (((long)RDB[DATA_BURN_STEP_PC] == CORRECTOR_STEP) &&
(RDB[DATA_BURN_SIE] == (double)NO))
return;

/* Set file name */

/* Check if running SIE burnup calculation */
if (RDB[DATA_BURN_SIE] == (double)YES)
{
if ((long)RDB[DATA_BURN_STEP_PC] == CORRECTOR_STEP)
sprintf(outfile, "%s_core%ld_CI_%ld.m", GetText(DATA_PTR_INPUT_FNAME),
(long)RDB[DATA_BURN_STEP], (long)RDB[DATA_BURN_CI_I]+1);
else
sprintf(outfile, "%s_core%ld.m", GetText(DATA_PTR_INPUT_FNAME),
(long)RDB[DATA_BURN_STEP]);
}
else
sprintf(outfile, "%s_core%ld.m", GetText(DATA_PTR_INPUT_FNAME),
(long)RDB[DATA_BURN_STEP]);



After that modification, now I could see a file created after each corrector iteration. However, all the power values in those files are zeros. It seems that the power is scored only at the predictor step.

Could you please help me out to get the power values after each corrector iteration? Which subroutine should be modified?

Thanks
Yousef Alzaben
Ph.D. Student
Karlsruhe Institute of Technology (KIT)

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

Re: Power Distribution for each burnup step using SIE method

Post by Ville Valtavirta » Thu Jun 22, 2017 8:53 am

Hi Yousef,

try modifying the if-statement on line 551 of score.c to your liking. Currently it skips scoring things like core power distribution, meshplots etc. on corrector steps.

-Ville

Yousef1986
Posts: 73
Joined: Fri Sep 04, 2015 3:33 pm
Security question 1: No
Security question 2: 92
Location: Karlsruhe, Germany
Contact:

Re: Power Distribution for each burnup step using SIE method

Post by Yousef1986 » Thu Jun 22, 2017 12:57 pm

Hi Ville,

I have modified the if-statement on line 531 and 551 of score.c to skip the return statement in the case of SIE method and now it works fine.
Thanks for the advice.
Yousef Alzaben
Ph.D. Student
Karlsruhe Institute of Technology (KIT)

Post Reply