université de montréal développement de la méthode des ...

28
Bien qu’un schéma de Jacobi converge moins vite qu’un schéma de Gauss-Seidel,
dans le cas de la méthode des caractéristiques, si le tracking est stocké dans un
fichier résidant sur le disque dur, il peut être très pénalisant de devoir lire le fichier
pour chaque groupe et une approche vectorielle est préférable.
Dans ce contexte, les itérations internes et thermiques peuvent être avantageusement
combinées dans un schéma de Jacobi à un niveau sous la forme
⃗Φ g(m+1) = L g (⃗ F g +
G∑
g ′ =1
∑ g←g ′
s
⃗Φ g′ (m)
)
. (2.38)
Cette approche, introduite dans le solveur MCI [Wu & Roy, 2003b], a été testée et
a été retenue pour ce projet. Même à peu de groupes (e.g. 20), séparer itérations
internes et itérations multigroupes n’est pas avantageux.
Dans les contextes de l’auto-protection des résonances ou du calcul SPH, les itérations
internes ne sont pas combinées avec un niveau supérieur d’itérations.
2.5 Tenants et aboutissants du développement de MOC
Dans ce paragraphe, nous allons insister sur les différences entre la méthode des
caractéristiques et la méthode des probabilités de collision car ces deux techniques
sont directement en concurrence dans le cadre du calcul de réseau.
D’abord, un point important qui explique le développement important que connait
MOC alors que les calculs se font avec des géométries de plus en plus détaillées,
concerne le stockage requis. En effet, comme le formalisme introduit précedemment
le montre, d’un point de vue du calcul de flux, le stockage est en O(N L + M) i.e.
des flux et des courants par régions. Ceci est à comparer avec CP qui, quant à
elle, nécessite un stockage en O(NL 2 ) i.e. les matrices denses qui contiennent les