28.08.2014 Views

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

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

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

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

viii<br />

ABSTRACT<br />

This project is <strong>de</strong>dicated to the implementation of the method of characteristics<br />

for <strong>la</strong>ttice calcu<strong>la</strong>tions in the neutronics mo<strong>de</strong>ling of nuclear reactors. A <strong>la</strong>ttice<br />

calcu<strong>la</strong>tion is a sequence of mo<strong>de</strong>ls that require the solution of the neutron transport<br />

equation at different stages. Within this framework, the collision probability<br />

method has been the <strong>de</strong>dicated tool for a long time. Nowadays, for the multigroup<br />

flux calcu<strong>la</strong>tion, the method of characteristics tends to supp<strong>la</strong>nt collision probability<br />

approaches thanks to its practical capability of treating configurations with a<br />

<strong>la</strong>rge number of regions and anisotropic scattering. However, for the other stages<br />

of a <strong>la</strong>ttice calcu<strong>la</strong>tion, the collision probability method remains the only solver<br />

used. In this work, beyond the computational <strong>de</strong>velopment re<strong>la</strong>ted to the method<br />

of characteristics in itself, the main i<strong>de</strong>a is to extend the usage of this method to<br />

all the calcu<strong>la</strong>tion stages that involve the solution of the transport equation.<br />

In this context, the method of characteristics was implemented for 2D geometries<br />

with exact and approximated boundary conditions. Special care was taken for<br />

the <strong>de</strong>velopment of acceleration techniques adapted to the different contexts. This<br />

acceleration strategy is based on two general c<strong>la</strong>sses of methods, synthetic preconditioning<br />

techniques and Krylov iterative methods, by analogy with the resolution<br />

of <strong>la</strong>rge linear systems. An already existing preconditioner was selected ; the work<br />

has consisted in <strong>de</strong>tailing the fundamental assumption of this method, performing a<br />

performances analysis prior and posterior to its implementation and, starting from<br />

there, improving it. A Krylov method was chosen and its implementation with the<br />

previous preconditioner in the special context of the characteristic method has been<br />

carried out. Besi<strong>de</strong>s, for comparison purpose, another preconditioning technique<br />

and another Krylov method were implemented.<br />

Apart from these computational <strong>de</strong>velopments, some re<strong>la</strong>ted investigations were<br />

carried out in or<strong>de</strong>r to c<strong>la</strong>rify some theoretical and practical aspects re<strong>la</strong>ted to the

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

Saved successfully!

Ooh no, something went wrong!