No description available.

This is Part B. / <p> Another approach is adopted for deriving the moments equations in spherical geometry using a spherical harmonics expansion of the neutron transport equation over a variable range of the direction cosine. Because of complications and uncertainties in establishing boundary conditions for the equations, only the zero'th order equations are solved, in an idealized situation, in order that a feel for equations and boundary conditions may be obtained.</p> <p> The equations are compared to equations given in a paper 'Directionally Discontinuous Harmonic Solutions of the Neutron Transport Equation in Spherical Geometry', by A. A. Harms and E. A. Attia. Analytical solutions for the zero'th order equations are given for equations developed there and to the equations developed in this paper. Numerical values are presented to give an idea of what accuracies might be expected. It is hoped that similar techniques can be used to solve the higher order equations analytically, and that appropriate boundary conditions can be found.</p> / Thesis / Master of Engineering (MEngr)

No description available.

No description available.

In this thesis, the parallel performance of the finite element-spherical harmonics (FE-PN) method implemented in the general-purpose radiation transport code EVENT is studied both analytically and empirically. EVENT solves the coupled set of space-angle discretized FE-PN equations using a parallel block-Jacobi domain decomposition method. As part of the analytical study, the thesis presents complexity results for EVENT when solving for a 3D criticality benchmark radiation transport problem in parallel. The empirical analysis is concerned with the impact of the main algorithmic factors affecting performance. Firstly, EVENT supports two solution strategies, namely MOD (Moments Over Domains) and DOM (Domains Over Moments), to solve the transport equation in parallel. The two strategies differ in the way they solve the multi-level space-angle coupled systems of equations. The thesis presents empirical evidence of which of the two solution strategies is more efficient. Secondly, different preconditioners are used in the Preconditioned Conjugate Gradient (PCG) inside EVENT. Performance of EVENT is compared when using three preconditioners, namely diagonal, SSOR(Symmetric Successive Over-Relaxation) and ILU. The other two factors, angular and spatial resolutions of the problem affect both the performance and precision of EVENT. The thesis presents comparative results on EVENTs performance as these two resolutions are increased. From the empirical performance study of EVENT, a bottleneck is identified that limits the improvement in performance as number of processors used by EVENT is increased. In some experiments, it is observed that uneven assignment of computational load among processors causes a significant portion of the total time being spent in synchronization among processors. The thesis presents two indicators that identify when such inefficiency occur; and in such a case, a load rebalancing strategy is applied that computes a new partition of the problem so that each partition corresponds to equal amount of computational load.

Diffusion MRI (dMRI) measurements allow us to infer the microstructural properties of white matter and to reconstruct fiber pathways in-vivo. High angular diffusion imaging (HARDI) allows for the creation of more and more complex local models connecting the microstructure to the measured signal. One of the challenges is the derivation of meaningful metrics describing the underlying structure from the local models. The aim hereby is to increase the specificity of the widely used metric fractional anisotropy (FA) by using the additional information contained within the HARDI data. A local model which is connected directly to the underlying microstructure through the model of a single fiber population is spherical deconvolution. It produces a fiber orientation density function (fODF), which can often be interpreted as superposition of multiple peaks, each associated to one relatively coherent fiber population (bundle). Parameterizing these peaks one is able to disentangle and characterize these bundles. In this work, the fODF peaks are approximated by Bingham distributions, capturing first and second order statistics of the fiber orientations, from which metrics for the parametric quantification of fiber bundles are derived. Meaningful relationships between these measures and the underlying microstructural properties are proposed. The focus lies on metrics derived directly from properties of the Bingham distribution, such as peak length, peak direction, peak spread, integral over the peak, as well as a metric derived from the comparison of the largest peaks, which probes the complexity of the underlying microstructure. These metrics are compared to the conventionally used fractional anisotropy (FA) and it is shown how they may help to increase the specificity of the characterization of microstructural properties. Visualization of the micro-structural arrangement is another application of dMRI. This is done by using tractography to propagate the fiber layout, extracted from the local model, in each voxel. In practice most tractography algorithms use little of the additional information gained from HARDI based local models aside from the reconstructed fiber bundle directions. In this work an approach to tractography based on the Bingham parameterization of the fODF is introduced. For each of the fiber populations present in a voxel the diffusion signal and tensor are computed. Then tensor deflection tractography is performed. This allows incorporating the complete bundle information, performing local interpolation as well as using multiple directions per voxel for generating tracts. Another aspect of this work is the investigation of the spherical harmonic representation which is used most commonly for the fODF by means of the parameters derived from the Bingham distribution fit. Here a strong connection between the approximation errors in the spherical representation of the Dirac delta function and the distribution of crossing angles recovered from the fODF was discovered. The final aspect of this work is the application of the metrics derived from the Bingham fit to a number of fetal datasets for quantifying the brain’s development. This is done by introducing the Gini-coefficient as a metric describing the brain’s age.

Call number: LD2668 .T4 1964 C15 / Master of Science

An important aspect for precision guidance, navigation, and control for lunar operations is environmental modeling. In particular, consider gravity field modeling. Available gravity field models for the Moon reach degree and order 165 requiring the use and storage of approximately 26,000 spherical harmonic coefficients. Although the high degree and order provide a means by which to accurately predict trajectories within the influence of the Moon's gravitational field, the size of these models makes using them computationally expensive and restricts their use in design environments with limited computer memory and storage. It is desirable to determine reduced complexity realizations of the gravitational models to lower the computational burden while retaining the structure of the original gravitational field for use in rapid design environments. The extended Kalman filter and the unscented Kalman filter are used to create reduced order models and are compared against a simple truncation based reduction method. Both variations of the Kalman filter out perform the truncation based method as a means by which to reduce the complexity of the gravitational field. The extended Kalman filter and unscented Kalman filter were able to achieve good estimates of position while reducing the number of spherical harmonic coefficients used in gravitational acceleration calculations by approximately 5,400, greatly increasing the speed of the calculations while reducing the required computer allocation. / text

A significant component of recent space exploration has been unmanned mission to comets and asteroids. The increase in interest for these bodies necessitates an increase in demand for higher fidelity trajectory simulations in order to assure mission success. Most available methods for simulating trajectories about asymmetric bodies assume they are of uniform density. This thesis examines a hybrid method that merges a mass concentration ("mascon") model and a spherical harmonic model using the "Brillouin sphere" as the interface. This joint model will be used for simulating trajectories about variable density bodies and, in particular, contact binary asteroids and comets. The scope of this thesis is confined to the analysis and characterization of the spherical harmonic modeling method in which three bodies of increasing asymmetrical severity are used as test cases: Earth, asteroid 101955 Bennu, and asteroid 25143 Itokawa. Since the zonal harmonics are well defined for Earth, it is used as the initial baseline for the method. Trajectories in the equatorial plane and inclined to this plane are simulated to analyze the dynamical behavior of the environment around each of the three bodies. There are multiple degrees of freedom in the spherical harmonic modeling method which are characterized as follows: (1) The radius of the Brillouin sphere is varied as a function of the altitude of the simulated orbit, (2) The truncation degree of the series is chosen based upon the error incurred in the acceleration field on the chosen Brillouin sphere, and (3) The gravitational potential and acceleration field are calculated using the determined radial location of the Brillouin sphere and the truncation degree. An ideal Brillouin sphere radius and truncation degree are able to be determined as a function of a given orbit where the error in the acceleration field is locally minimized. The dual-density model for a contact binary is found to more accurately describe the dynamical environment around Asteroid 25143 Itokawa compared to the single density model.

Avec les progrès de la technologie, la densité de puissance des appareils électriques augmente, ce qui conduit à une augmentation des interférences électromagnétiques entre deux systèmes d’électroniques de puissance. Ces interférences peuvent engendrer un dysfonctionnement des appareillages. Afin d’éviter de tels problèmes, une grande quantité de vérifications expérimentaux sont nécessaires après le développement d’un prototype, ce qui conduira à un surcoût important dans son cycle industriel. Cette thèse a pour objectif de prendre en compte des problèmes de compatibilité électromagnétique en champ proche dès la phase de conception d’un produit. Ceci est basé sur le développement en harmoniques sphériques du rayonnement parasite engendré par les dispositifs ou les composants d’un dispositif. Les coefficients harmoniques d’un dispositif permettent non seulement de modéliser ses champs proches par une source ponctuelle, mais également de déterminer le couplage inductif avec les autres. Dans ce contexte, la méthode du développement multipolaire est étudiée théoriquement, numériquement (par le logiciel Flux3D) et expérimentalement dans ce travail. Dans ce document, le nouveau système de mesure automatisé est présenté et différentes approches liées à la solution de ce problème sont étudiées, tels que le choix optimal d’origine du développement, la compensation d’erreurs des mesures dues aux sources secondaires, l’étude de propagation d’incertitude dans le problème inverse et prise en compte d’informations a priori, etc. / With the advances in technology, the power density of electrical equipment is increasing, which leads to an increase in electromagnetic interference between two power electronic systems. This interference may cause the electronic equipment to malfunction. In order to avoid such problems, a large number of experimental verifications are necessary after the development of a prototype, which leads to a significant additional cost in its industrial cycle. This thesis aims at taking into account the problems of near-field electromagnetic compatibility in the design phase of a product. The method is based on the development in spherical harmonics of the parasitic radiation generated by the devices or the components of the device. The harmonic coefficients of a device make it possible not only to model its near field by a punctual source, but also to determine the inductive coupling with the others. In this context, the multipole expansion method is studied theoretically, numerically (by the software Flux3D) and experimentally in this work. In this document, the new automated measuring system is presented and different approaches related to the solution of this problem are then studied, such as the optimal choice of the development origin, the error compensation of the measurements due to the secondary sources, the study of the propagation of uncertainty in the inverse problem and how to take into account the a priori information, etc.