This thesis is a blend of neutron transport theory and numerical analysis. We start with the study of the problem of the Mika/Case eigenexpansion used in the solution process of the homogeneous one-speed Boltzmann neutron transport equation with anisotropic scattering for plane symmetry. The anisotropic scattering is expressed as a finite Legendre series in which the coefficients are the ``scattering coefficients'. This eigenexpansion consists of a discrete spectrum of eigenvalues with its corresponding eigenfunctions and the continuous spectrum [-1,+1] with its corresponding eigendistributions. In the general case where the anisotropic scattering can be of any (finite) order, multiple discrete eigenvalues exist and these have to be located to have the complete spectrum. We have devised a stable and robust method that locates all these discrete eigenvalues. The method is a two-step process: first the number of discrete eigenvalues is calculated and this is followed by the calculation of the discrete eigenvalues themselves, now being able to count them down and make sure none are forgotten.
During our numerical experiments, we came across what we called near-singular eigenvalues: discrete eigenvalues that are located extremely close to the continuum and hence lead to near-singular behaviour in the eigenfunction. Our solution method has been adapted and allows for the automatic detection of such a near-singular eigenvalue.
For the elements of the continuous spectrum [-1,+1], there is no non-zero function satisfying the associated eigenequation but there is a non-zero distribution that does satisfy it. It is not feasible to compute a distribution as such but one can evaluate integrals in which this distribution appears. The continuum part of the eigenexpansion can hence only be characterised by its (angular) moments. Accurate and fast numerical quadrature is needed to evaluate these integrals. Several quadrature methods have been evaluated on a representative test function.
The eigenexpansion was proved to be orthogonal and complete and hence can be used to represent the infinite medium Green's function. The latter is the building block of the Boundary Sources Method, an integral solution method for the neutron transport equation. Using angular and angular/spatial moments of the Green's function, it is possible to solve with high accuracy slab problems. We have written a one-dimensional slab code implementing this Boundary Sources Method allowing for media with arbitrary order anisotropic scattering. Our results are very good and the code can be considered as a benchmark code for others.
As a final application, we have used our code to study the discrete spectrum of a well-known scattering kernel in radiative transfer, the Henyey-Greenstein kernel. This kernel has one free parameter which is used to fit the kernel to experimental data. Since the kernel is a continuous function, a finite Legendre approximation needs to be adopted. Depending on the free parameter, the approximation order and the number of secondaries per collision, the number of discrete eigenvalues ranges from two to thirty and even more. Bounds for the minimum approximation order are derived for different requirements on the approximation: non-negativity, an absolute and relative error tolerance.
Identifer | oai:union.ndltd.org:BICfB/oai:ulb.ac.be:ETDULB:ULBetd-04252005-114848 |
Date | 12 May 2005 |
Creators | Van den Eynde, Gert |
Contributors | Beauwens, Robert, Mund, Ernest, Mika, Janusz, Ganapol, Barry, Notay, Yvan, Labeau, Pierre-Etienne, Louis, Yves, Aït Abderrahim, Hamid |
Publisher | Universite Libre de Bruxelles |
Source Sets | Bibliothèque interuniversitaire de la Communauté française de Belgique |
Language | English |
Detected Language | English |
Type | text |
Format | application/pdf |
Source | http://theses.ulb.ac.be/ETD-db/collection/available/ULBetd-04252005-114848/ |
Rights | unrestricted, J'accepte que le texte de la thèse (ci-après l'oeuvre), sous réserve des parties couvertes par la confidentialité, soit publié dans le recueil électronique des thèses ULB. A cette fin, je donne licence à ULB : - le droit de fixer et de reproduire l'oeuvre sur support électronique : logiciel ETD/db - le droit de communiquer l'oeuvre au public Cette licence, gratuite et non exclusive, est valable pour toute la durée de la propriété littéraire et artistique, y compris ses éventuelles prolongations, et pour le monde entier. Je conserve tous les autres droits pour la reproduction et la communication de la thèse, ainsi que le droit de l'utiliser dans de futurs travaux. Je certifie avoir obtenu, conformément à la législation sur le droit d'auteur et aux exigences du droit à l'image, toutes les autorisations nécessaires à la reproduction dans ma thèse d'images, de textes, et/ou de toute oeuvre protégés par le droit d'auteur, et avoir obtenu les autorisations nécessaires à leur communication à des tiers. Au cas où un tiers est titulaire d'un droit de propriété intellectuelle sur tout ou partie de ma thèse, je certifie avoir obtenu son autorisation écrite pour l'exercice des droits mentionnés ci-dessus. |
Page generated in 0.0017 seconds