Project

General

Profile

Forums » Support »

Too many steps in routine RungeKutta

Added by YANG Richard Yi almost 9 years ago

Je veux simuler un sénario simple : 4 REP utilisent toujours 3,2% UOX et un REP utilise les URE sortent de 4 REP précédents. Comme les compositions de tous les réacteurs ne varient pas, seulement 2 .dat fichiers sont offerts : un fichier pour le database de 4 REP de UOX et un fichier pour XSM_CLOSEST.
Mais quand CLASS_exec est exécuté, il y a des messages dans le terminal "Too many steps in routine RungeKutta".

De plus, il y a des isotopes qui doient exister pendant l'évolution dans le REP de URE disparaissent dans le résultat (par exemple, Pu243).

Je joins notre fichier d'entrée et le .log fichier.

Download all files
SimpleExample.cxx (6.21 KB) SimpleExample.cxx
CLASS_OUTPUT.log (4.33 MB) CLASS_OUTPUT.log

Replies (9)

RE: Too many steps in routine RungeKutta - Added by MOUGINOT Baptiste almost 9 years ago

Je ne sais pas exactement comment fonctionne tes URE mais :

"Too many steps in routine RungeKutta"

vient d'un problème de pas en temps dans tes calculs : généralement la puissance du réacteur est mal défini. (doit être donné en Watt)

J'ai effectivement l'impression que tu t'es tromper dans ton calcul de puissance : tu prend la puissance spécifique en MW/tHM mais dans ton calcul tu fais : 

power[i]=power[i]*HMM[i]*1.0e6;

je pense que ca serait correct avec :

power[i]=power[i]*HMM[i]*1.0e3;

!!WARNING!! [IrradiationModel.cxx:Getn2nXsMatrix] Can't have n,2n reaction on this nuclei, ZAI : 64 156 0
!!WARNING!! [IrradiationModel.cxx:GetNuclearProcessMatrix] Can't have nuclear process on this nucleus, ZAI 94 243 0 : its halflife seems to be below the threshold!

tu n'as pas à t'inquiéter à propos de ses messages : CLASS t'informe juste que (1) il n'as pas de section efficaces pour cette reaction : elle est donc ignorée ou que (2) le noyau créé a un temps de vie inférieur à la limite fixée ( il est normalement traiter correctement!)

Bon courage
BaM

RE: Too many steps in routine RungeKutta - Added by YANG Richard Yi almost 9 years ago

1MW=1,0e6W, donc je pense que la puissance est correcte, n'est-ce pas ?

C'est étrange parce que tous ces sections efficaces existent dans .dat fichier, et quand le .dat fichier de RU est offert à un réacteur de "fixed fuel", il n'y a pas de problème.

RE: Too many steps in routine RungeKutta - Added by LENIAU Baptiste almost 9 years ago

Bonjour Richard,

Tu pourrais envoyer tes fichiers .dat de réacteurs que j'essaye de reproduire le bug.

merci

RE: Too many steps in routine RungeKutta - Added by YANG Richard Yi almost 9 years ago

Bonjour Monsieur,
Voici des databases de réacteurs. Le fichier de UOX est pour les réacteurs de "fixed fuel" et le fichier de RU est pour le calcul de XSM_CLOSEST.
Merci.

Download all files
UOX_3.2_40.17_33.5.info (199 Bytes) UOX_3.2_40.17_33.5.info
UOX_3.2_40.17_33.5.dat (355 KB) UOX_3.2_40.17_33.5.dat
RU_3.2.info (199 Bytes) RU_3.2.info
RU_3.2.dat (355 KB) RU_3.2.dat

RE: Too many steps in routine RungeKutta - Added by MOUGINOT Baptiste almost 9 years ago

effectivement 1MW=1e6W... donc c'est correct... désolé!

J'ai vérifié la valeur de ses sections efficaces : l'ordre de grandeur est bon, tes .info semble correct aussi... ton input semble correct aussi à première vue....

je regarde ca dans la journée si Baptiste L. n'a pas trouvé ...

RE: Too many steps in routine RungeKutta - Added by MOUGINOT Baptiste almost 9 years ago

peux-tu poster également ton EQM_Simple.hxx/cxx ?

RE: Too many steps in routine RungeKutta - Added by YANG Richard Yi almost 9 years ago

Ce EQM prend simplement une propotion des isotopes fissiles fixée que j'ai calculé.

Download all files
EQM_Simple.cxx (673 Bytes) EQM_Simple.cxx
EQM_Simple.hxx (309 Bytes) EQM_Simple.hxx

RE: Too many steps in routine RungeKutta - Added by LENIAU Baptiste over 8 years ago

je réécris le mail que je t'ai envoyé pour avoir une trace de la solution :

j’ai trouvé le problème.
Il vient du fait que le fichier d'EvolutionData "RU_3.2.dat" possède un espace à la fin de chaque ligne alors qu’il ne faut pas.
La lecture de ces fichiers est un peu sensible.

SOLUTION : enlever l’espace à la fin de la ligne du temps dans RU_3.2.dat

à bientôt ,

    (1-9/9)