10 - Propagation des rayonnements : méthodes et codes de calcul
327
L’inconvénient majeur de la simulation Monte Carlo est sa lenteur de convergence. Toutefois, on a vu que l’emploi des techniques de réduction de la variance permet de pallier
efficacement cette difficulté intrinsèque de la méthode.
Par ailleurs, la puissance accrue des ordinateurs et les techniques de programmation qui
« parallélisent » les logiciels de calcul, réduisent encore sensiblement les temps de calcul
et rendent de plus en plus attractif l’emploi de la méthode de Monte Carlo.
L’ensemble de ces dispositifs et les perspectives d’évolution qu’ils offrent conduisent à
faire évoluer le développement et l’usage des codes Monte Carlo dans deux directions :
celle de la R&D multi-physiques et celle du calcul industriel.
Plusieurs codes de transport tridimensionnel des neutrons et de γ fondés sur la méthode de
Monte Carlo existent : MCNP, TRIPOLI, MVP, MCU, MACBEND, une liste référencée en
est donnée dans le tableau 10.1. Selon les codes et leurs options de calcul, la représentation des sections efficaces en fonction de l’énergie est ponctuelle et/ou muligroupe, et/ou
multigroupe-tables de probabilité (Zheng et al., 1998).
Figure 10.28. Illustration de la convergence d’un débit d’équivalent de dose associé au flux des
neutrons d’énergie initiale de 2 MeV ayant traversé 70 cm d’une protection en béton. Les traits
verticaux représentent les barres d’erreur statistique qui sont de grande amplitude pour la simulation
analogue réduite pour les simulations biaisées (transformation exponentielle sur ces exemples). Ces
courbes montrent donc l’efficacité des techniques de réduction de la variance l’optimum étant
obtenu pour une valeur du paramètre de biaisage κ avoisinant 0,15 (Vergnaud, 2000).
327
L’inconvénient majeur de la simulation Monte Carlo est sa lenteur de convergence. Toutefois, on a vu que l’emploi des techniques de réduction de la variance permet de pallier
efficacement cette difficulté intrinsèque de la méthode.
Par ailleurs, la puissance accrue des ordinateurs et les techniques de programmation qui
« parallélisent » les logiciels de calcul, réduisent encore sensiblement les temps de calcul
et rendent de plus en plus attractif l’emploi de la méthode de Monte Carlo.
L’ensemble de ces dispositifs et les perspectives d’évolution qu’ils offrent conduisent à
faire évoluer le développement et l’usage des codes Monte Carlo dans deux directions :
celle de la R&D multi-physiques et celle du calcul industriel.
Plusieurs codes de transport tridimensionnel des neutrons et de γ fondés sur la méthode de
Monte Carlo existent : MCNP, TRIPOLI, MVP, MCU, MACBEND, une liste référencée en
est donnée dans le tableau 10.1. Selon les codes et leurs options de calcul, la représentation des sections efficaces en fonction de l’énergie est ponctuelle et/ou muligroupe, et/ou
multigroupe-tables de probabilité (Zheng et al., 1998).
Figure 10.28. Illustration de la convergence d’un débit d’équivalent de dose associé au flux des
neutrons d’énergie initiale de 2 MeV ayant traversé 70 cm d’une protection en béton. Les traits
verticaux représentent les barres d’erreur statistique qui sont de grande amplitude pour la simulation
analogue réduite pour les simulations biaisées (transformation exponentielle sur ces exemples). Ces
courbes montrent donc l’efficacité des techniques de réduction de la variance l’optimum étant
obtenu pour une valeur du paramètre de biaisage κ avoisinant 0,15 (Vergnaud, 2000).