26.03.2015 Views

De la fission aux nouvelles filières - Cenbg - IN2P3

De la fission aux nouvelles filières - Cenbg - IN2P3

De la fission aux nouvelles filières - Cenbg - IN2P3

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Sur <strong>la</strong> figure ci-dessus, le jeu consiste en <strong>la</strong> réalisation de trajectoires de particules et le score<br />

est de 1 à chaque fois qu’une particule traverse <strong>la</strong> surface. On peut alors calculer <strong>la</strong> moyenne<br />

du nombre de particules détectées (ayant traversé <strong>la</strong> surface).<br />

Dans le cadre du transport de particules, l’utilisation d’une méthode de Monte Carlo pour<br />

estimer des quantités physiques telles que flux, t<strong>aux</strong> de réaction ou facteur de multiplication<br />

repose sur <strong>la</strong> simu<strong>la</strong>tion d’histoires de particules dans <strong>la</strong> matière. Une histoire commence par<br />

<strong>la</strong> naissance de <strong>la</strong> particule (émission à partir d’une source) et se termine par sa mort<br />

(absorption, fuite).<br />

9.1.2 Les codes de calcul.<br />

Les codes de calcul, mettant en œuvre <strong>la</strong> méthode de Monte Carlo, pour <strong>la</strong> résolution des<br />

problèmes en physique des réacteurs sont principalement les logiciels MCNP et TRIPOLI.<br />

9.2 Méthodes déterministes.<br />

L’équation de Boltzmann est résolue numériquement en appliquant <strong>la</strong> discrétisation des<br />

variables :<br />

E, r, Ω → E i , r j , Ω k<br />

i = 1 … I<br />

j = 1 … J<br />

k = 1 … K<br />

⎬<br />

I x J x K Variables<br />

La fonction<br />

Φ<br />

(E,r,Ω) devient un vecteur Φ (Φ ijk ), où <strong>la</strong> composante Φ ijk est <strong>la</strong> valeur<br />

moyenne du flux dans l’élément △E i △r j △Ω k<br />

L’équation de Boltzmann peut se réécrire donc de façon discrétisée (en forme matricielle) :<br />

AΦ = FΦ A : opération de « disparition »<br />

⎨<br />

F : opération de « production »<br />

Formellement on cherche une solution du type (vecteur propre) :<br />

Φ = A -1 F Φ<br />

Différentes approximations sont utilisées dans l’opération de discrétisation en angle et<br />

espace :<br />

- Méthodes <strong>aux</strong> Ordonnées discrètes.<br />

- Méthodes Nodales<br />

- Etc.<br />

Ces approximations conduisent à des systèmes d’équations, résolus par des méthodes<br />

itératives, avec techniques d’accélération adaptées au type de matrice obtenu pour A et F.<br />

Un problème majeur est représenté par <strong>la</strong> variation en énergie des sections efficaces σ (E).<br />

Dans les matrices A et F, il faut introduire des valeurs moyennes ou « effectives » sur des<br />

« groupes » d’énergie △E i :<br />

280

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!