Optimizing launcher trajectories is essential for effective mission planning, and specialized software like ASTOS provide an initial, precise overview. However, as launcher development progresses, there is a growing need for the creation of an autonomous flight trajectory software that offers greaterflexibility in adjusting simulation parameters and better represents actual, real-life trajectories. Thisreport introduces an initial version of a comprehensive six-degree-of-freedom launcher trajectory calculation software developed using MATLAB and Simulink. The emphasis is on the development strategy, encompassing discussions on dynamics equations, essential features, and crucial models necessary for accurate simulations. Real-world scenarios often deviate from optimized trajectories, and the software addresses these deviations using sensitivity analysis through Monte Carlosimulations, enabling a thorough examination of uncertainties in input parameters and their impact on trajectories. The article delves into the establishment of the dispersion analysis tool and offers suggestions for further enhancements for both the Simulink model and this dispersion analysis tool. / Optimering av flygbanor är avgörande för effektiv uppdragsplanering, och specialiserad programvara som ASTOS ger en initial, exakt översikt. Men när flygbanans utveckling fortskrider finns det ett växande behov av att skapa en autonom flygbana som erbjuder större flexibilitet när det gäller attjustera simuleringsparametrar och bättre representerar faktiska, verkliga banor. Den här rapporten introducerar en initial version av en omfattande beräkningsprogramvara utvecklad med MATLAB och Simulink för sex frihetsgraders lanseringsbana. Tyngdpunkten ligger på utvecklingsstrategin, som omfattar diskussioner om dynamikekvationer, väsentliga funktioner och avgörande modeller som är nödvändiga för exakta simuleringar. Scenarier i verkligheten avviker ofta från optimerade banor, och programvaran adresserar dessa avvikelser med känslighetsanalys genom Monte Carlo-simuleringar,vilket möjliggör en grundlig undersökning av osäkerheter i inmatningsparametrar och deras påverkan på banor. Rapporten går in i skapandet av spridningsanalysverktyget och erbjuder förslag till ytterligare förbättringar för både Simulink-modellen och detta dispersionsanalysverktyg.

Electronic energy transfer has been investigated in pure donor systems by means of computer simulations. Calculated properties were the probability that the initially excited donor is excited at a time t after the excitation, Gs(t), the mean square displacement of the excitation and different fluorescence observables. For three dimensional systems the results obtained by Monte Carlo simulations were compared to the so-called GAF-theory {Gouchanour,C. R., Andersen, H. C. and Fayer, M. D., J. Chem. Phys. 81, 4380 (1984)}, and the agreement was found to be good. Anisotropic systems, i.e. mono-, bi- and multilayer systems, were compared to the two-particle model {Baumann,J. and Fayer, M. D., J. Chem. Phys. 85, 4087 (1986)}. The agreement between the Gs(t) calculated from the tp- model and the Monte Carlo simulations were good for all systems investigated. However, the agreement between the fluorescence observables obtained by MC and the tp-model were in general poor. A much better agreement was found when a phenomenological approach was used for calculating the fluorescence depolarization ratios. Three dimensional systems where the donors are rotating on the same time scale as the energy transfer takes place have also been studied and compared to analytical theories. The Molecular Dynamics simulations of decapeptide H142 shows that simulations in a continuum with a relative permeability do not provide a reliable alternative to simulations with explicit solvent molecules. / <p>Diss. (sammanfattning) Umeå : Umeå universitet, 1991, härtill 5 uppsatser</p> / digitalisering@umu

This thesis presents a methodology to incorporate reservoir uncertainties and estimate the loss in project value when facility planning decisions are based on erroneous estimates of input variables. We propose a tank model along with integrated asset development model to simulate the concept selection process. The model endogenizes drilling decisions and includes an option to expand. Key decision variables included in the model are number of pre-drill wells, initial facility capacity and number of well slots. Comparison is made between project value derived under erroneous estimates for reserve size and under an alternate hypothesis. The results suggest loss in project value of up to 40% when reservoir estimates are erroneous. Moreover, both optimistic and pessimistic reserve estimates results in a loss in project value. However, loss in project value is bigger when reserve size is underestimated than when it is overestimated. / text

In this thesis we introduce, implement and compare several multivariate goodness-of-fit tests. First of all, we will focus on universal mul- tivariate tests that do not place any assumptions on parametric families of null distributions. Thereafter, we will be concerned with testing of multi- variate normality and, by using Monte Carlo simulations, we will compare power of five different tests of bivariate normality against several alternati- ves. Then we describe multivariate skew-normal distribution and propose a new test of multivariate skew-normality based on empirical moment genera- ting functions. In the final analysis, we compare its power with other tests of multivariate skew-normality. 1

Les rayonnements ionisants sont connus pour induire des dommages critiques au sein de la matière biologique et spécialement au sein de l’ADN. Parmi ces dommages, les cassures doubles brins de l’ADN (DSB) sont considérées comme les principales responsables des effets létaux des rayonnements. Comprendre et prédire comment ces cassures sont créées et réparées dans les noyaux cellulaires demeure un défi dans la recherche en radiobiologie. Ce travail s’inscrit dans ce contexte, dans la modélisation des cassures double brin de l’ADN (DSB) à partir des dépôts d’énergie créés par l’irradiation au niveau intracellulaire. Le détail topologique au niveau nanométrique des dépôts d’énergie nécessaire à ce travail est obtenu par modélisation Monte Carlo à l’aide du code Geant4 et, en particulier son extension Geant4-DNA pour des processus à très faible énergie. Les dommages étudiés étant ceux localisés dans l’ADN, le premier objectif de ce travail a été de réaliser une géométrie détaillée de celui-ci afin de l’implémenter dans les calculs Monte Carlo. Deux types de noyaux cellulaires, représentant un fibroblaste et un endothélium, ont été décrits afin d’évaluer l’influence de la densité d’ADN dans les résultats sur la topologie des dépôts pouvant donner lieux à des cassures de la molécule. Cette géométrie nous permet d’effectuer une première sélection des dépôts d’énergie pouvant contribuer aux cassures car situées sur la chaîne sucre-phosphate. Ces dépôts sont ensuite analysés à l’aide d’un algorithme de clustérisation de manière à les regrouper sous forme d’agrégats afin d’étudier leur localisation et complexité. Néanmoins, dans cette étude, seule les interactions physiques entre les rayonnements ionisants et la cible sont modélisées, il n’est donc pas possible d’obtenir un nombre absolu de cassures de brins car cette modélisation n’inclue pas l’étape de création et de transport des radicaux libres pouvant donner lieu à des dommages indirects. Ainsi, le but de ce travail était d’évaluer la dépendance relative des dommages radio-induits directs avec la densité d’ADN, la qualité du rayonnement, la morphologie du noyau ou encore la condensation de la chromatine. Les différentes modélisations réalisées ont permis de quantifier l’influence de ces différents paramètres dans le nombre et la complexité des dommages directs induits dans l’ADN, pouvant ensuite contribuer aux effets tardifs sur le devenir cellulaire. / Ionizing radiations are known to induce critical damages on biological matter and especially on DNA. Among these damages, DNA double strand breaks (DSB) are considered as key precursor of lethal effects of ionizing radiations. Understand and predict how DNA double and simple strand breaks are created by ionising radiation and repaired in cell nucleus is nowadays a major challenge in radiobiology research. This work presents the results on the simulation of the DNA double strand breaks produced from the energy deposited by the irradiation at the intracellular level. At the nanometric scale, the only method to accurately simulate the topological details of energy deposited on the biological matter is the use of Monte Carlo codes. In this work, we used the Geant4 Monte Carlo code and, in particular, the low energy electromagnetic package extensions, referred as Geant4-DNA processes.In order to evaluate DNA radio-induced damages, the first objective of this work consisted in implementing a detailed geometry of the DNA on the Monte Carlo simulations. Two types of cell nuclei, representing a fibroblast and an endothelium, were described in order to evaluate the influence of the DNA density on the topology of the energy deposits contributing to strand breaks. Indeed, the implemented geometry allows the selection of energy transfer points that can lead to strand breaks because they are located on the backbone. Then, these energy transfer points were analysed with a clustering algorithm in order to reveal groups of aggregates and to study their location and complexity.In this work, only the physical interactions of ionizing radiations are simulated. Thus, it is not possible to achieve an absolute number of strand breaks as the creation and transportation of radical species which could lead to indirect DNA damages is not included. Nevertheless, the aim of this work was to evaluate the relative dependence of direct DNA damages with the DNA density, radiation quality, cell nuclei morphology or also chromatin condensation. The results presented in this work have allowed the quantification of the influence of these different parameters in the number and complexity of directs DNA damages which can then contribute to the late effects on cell fate.

Master of Science / Department of Electrical and Computer Engineering / Anil Pahwa / Unpredictable power outages due to environmental factors such as lighting, wind, trees, and animals, have always been a concern for utilities because they are often unavoidable. This research aims to study squirrel-related outages by modeling past real-life outage data and provide the optimal result which would assist utilities in increasing electric system reliability. This research is a novel approach to benchmark system performance in order to identify areas and durations with higher than expected outages. The model is illustrated with seven years (2005-2011) of animal-related outage data and 14 years of weather data (1998-2011) for four cities in Kansas, used as training data to predict future outages. The past data indicates that the number of outages on any day varies with the seasons and weather conditions on that day. The prediction is based on a Bayesian Model using conditional probability table, which is calculated based on training data. Since future weather conditions are unknown and random, Monte Carlo Simulation is used with the past 14 years of weather data to create different yearly scenarios. These scenarios are then used with the models to predict expected outages. Multiple runs of Monte Carlo analysis provide a probability distribution of expected outages. Further work discusses about cost-to-benefit analysis of implementation of outage mitigation methods. The analysis is performed by considering different combinations of outage reduction and mitigation levels. In this research, eight cases of outage reduction and nine cases of mitigation levels are defined. The probability of benefit is calculated by a statistical approach for every combination. Several optimal strategies are constructed using the probability values and outage history. The outcomes are compared with each other to propose the most beneficial outage mitigation strategy. This research will immensely assist utilities in reducing the outages due to squirrels more effectively with higher benefits and therefore improve reliability of the electricity supply to consumers.

Les méthodes de simulation Monte-Carlo s’étendent avec succès à différents domaines de la physique médicale mais aussi à différentes échelles, par exemple de la planification des traitements de radiothérapie jusqu’à une prévision des effets des rayonnements au niveau des cellules cancéreuses. La plateforme de simulation Monte-Carlo GATE, basée sur l’outil Geant4, propose des fonctionnalités dédiées aux simulations en physique médicale (médecine nucléaire et radiothérapie). Pour les applications en radiobiologie, les modèles physiques Geant4-DNA implémentés jusqu’à très basse énergie (eV) permettent d’estimer des quantités micro-dosimétriques d’intérêt. Dans le but d’implémenter une plateforme de simulation Monte-Carlo multi-échelles, nous nous sommes d’abord intéressés à la validation des modèles physiques de Geant4-DNA, puis à leur intégration dans la plateforme de simulation GATE et enfin à une validation de cette implémentation dans un contexte de radiothérapie et protonthérapie. De manière à valider les modèles physiques de Geant4-DNA, des points kernels de dose en électrons mono-énergétiques (de 10 keV à 100 keV) ont été simulés en utilisant les modèles physiques de Geant4 et de Geant4-DNA et ils ont comparés au code Monte-Carlo EGSnrc. Les parcours et pouvoirs d’arrêts des électrons (de 7,4 eV à 1 MeV) et des protons (de 1 keV à 100 MeV) calculés avec Geant4-DNA (processus et modèles préalablement intégrés dans GATE) ont ensuite été validés. Nous avons alors proposé de simuler avec la plateforme GATE l’impact de faisceaux cliniques et pré-cliniques sur l’ADN cellulaire. Nous avons ainsi modélisé un faisceau de protonthérapie de 193,1 MeV, un accélérateur linéaire en mode électrons de 6 MeV et un irradiateur RX de 250 kV. Ces simulations ont d’abord été validées en milieu aqueux par une comparaison de la dose macroscopique avec des mesures expérimentales. Les faisceaux ont ensuite été utilisés pour calculer, pour chacun d’entre eux, les fréquences de dépôts d’énergie à l’ADN. La molécule d’ADN a été simulée tout d’abord grâce à des cylindres équivalents en dimension à 10 paires de base (2 nm x 2 nm), équivalents à la taille d’un nucléosome (10 nm x 5 nm) et équivalents à la taille d’une fibre de chromatine (25 nm x 25 nm). Tous ces cylindres ont été placés aléatoirement dans un volume d’eau liquide (de rayon 500 nm). Nous avons ensuite reconstruit la molécule d’ADN dans Geant4 à partir de la lecture de fichiers PDB (Protein Data Bank) représentant douze paires de base de la molécule d’ADN et un dinucléosome (347 paires de base). Enfin, nous avons développé un outil permettant de corréler les positions de dépôts d’énergie directs dans l’eau liquide avec les coordonnées des paires de base de l’ADN, afin de calculer les nombres de cassures simple et double brin de l’ADN. Tous les calculs réalisés au cours de ce travail, ont été déployés sur l’Infrastructure de Grille Européenne ; des tests de performance sont proposés pour mesurer l’intérêt de ce type d’architecture pour les calculs Monte-Carlo. / The Monte Carlo simulation methods are successfully being used in various areas of medical physics but also at different scales, for example, from the radiation therapy treatment planning systems to the prediction of the effects of radiation in cancer cells. The Monte Carlo simulation platform GATE based on the Geant4 toolkit offers features dedicated to simulations in medical physics (nuclear medicine and radiotherapy). For radiobiology applications, the Geant4-DNA physical models are implemented to track particles till very low energy (eV) and are adapted for estimation of micro-dosimetric quantities. In order to implement a multi-scale Monte Carlo platform, we first validated the physical models of Geant4-DNA, and integrated them into GATE. Finally, we validated this implementation in the context of radiation therapy and proton therapy. In order to validate the Geant4-DNA physical models, dose point kernels for monoenergetic electrons (10 keV to 100 keV) were simulated using the physical models of Geant4-DNA and were compared to those simulated with Geant4 Standard physical models and another Monte Carlo code EGSnrc. The range and the stopping powers of electrons (7.4 eV to 1 MeV) and protons (1 keV to 100 MeV) calculated with GATE/Geant4-DNA were then compared with literature. We proposed to simulate with the GATE platform the impact of clinical and preclinical beams on cellular DNA. We modeled a clinical proton beam of 193.1 MeV, 6 MeV clinical electron beam and a X-ray irradiator beam. The beams models were validated by comparing absorbed dose computed and measured in liquid water. Then, the beams were used to calculate the frequency of energy deposits in DNA represented by different geometries. First, the DNA molecule was represented by small cylinders : 2 nm x 2 nm ( 10 bp), 5 nm x 10 nm ( nucleosome) and 25 nm x 25 nm ( chromatin fiber). All these cylinders were placed randomly in a sphere of liquid water (500 nm radius). Then we reconstructed the DNA molecule in Geant4 by reading PDB (Protein Data Bank) files representing twelve base pairs of the DNA molecule and a dinucleosome (347 base pairs). Finally, we developed a tool to correlate the positions of direct energy deposit in liquid water with the coordinates of the base pairs of DNA to calculate the number of single and double strand breaks in DNA. All calculations in this work were perfomed on the European Grid Infrastructure; performance tests are available to estimate the utility of this type of architecture for Monte Carlo calculations.

Bayer processing of bauxite employed for production of alumina yield a residue red-mud. The worldwide annual-rate of red-mud generation is approximately 120 million tons, and most of this is stockpiled. Red mud is rich in elements like aluminum, titanium and rare earths, in addition to the major iron-bearing constituents. The objective of this research is to explore such a strategy to extract Iron as fine particulate magnetite. Pyrometallurgical reduction experiments using carbon mixtures and a novel hydrometallurgical route are experimented. Reduction experiments performed with petroleum coke as a reductant resulted in incomplete reduction. The ‘optimal conditionsÂ’ for gaseous state reduction are determined to be: a processing temperature of 540oC ± 10oC, partial pressures CO (g) and CO2 (g) each of 0.070atm (bar) ± 0.001atm.(bar)/ inert diluent-gas: N2 (g), for a conversion-time of 30min. A mathematical-model was developed on the basis of unidirectional-diffusion of CO2 (g) within the CO2 (g)–CO (g)–N2 (g) gas-phase of the porous product-layer. Magnetic separation applied to the optimally reduced sample in: obtaining a magnetic portion with high iron and non-magnetic portion containing nonferrous (Al, Ti) is not successful. This finding was subsequently attributed to the nanometer length-scales crystallites of the predominant iron-containing phase, hydrated ferric-oxide(s) as determined by STEM micrographs. In addition, the presence of substitution for Fe3+ by Al3+ and Ti3+/4+ are determined with the help of MÖssbauer spectrograms. A hydrometallurgical route involving selective leaching and precipitation of iron in red-mud is tried. Red-mud is leached in oxalic acid at: 95˚C, 15 % Pulp density, 2.5 h leaching time, 1.5 pH. Kinetic studies yielded the leaching mechanism to be predominantly fluid film control. Ferric oxalate in the leach liquor is reduced to insoluble ferrous oxalate selectively using Iron powder. The ferrous oxalate formed is ~98 % pure. The precipitated ferrous oxalate is magnetically separated and reduced in Nitrogen atmosphere to form fine particulate magnetite. Additionally an economic feasibility study was conducted for the hydro and pyro alternatives to produce magnetite using Monte Carlo Simulations by imparting uncertainties in various input cost components. A traditional process was also compared to the proposed approaches for Total Capital Investment (TCI), Total Product Cost (TPC), Net Present Value (NPV) and sensitivity analysis. ~52 % reduction in Total product cost, 46 % reduction in Total Capital Investment was achieved for the hydrometallurgical process when compared to the traditional approach.

Neste trabalho estudamos um modelo estocástico com simetria Ising e dois estados absorventes em três dimensões com uma rede cúbica e em duas dimensões através de uma rede triangular. O estudo levou em conta cálculos de aproximação de campo médio e simulações de Monte Carlo. Os resultados mostraram que o modelo tem transição de segunda ordem de uma fase paramagnética para uma fase ferromagnética, uma transição da fase ferromagnética para uma fase absorvente, também de segunda ordem, e ainda uma transição de primeira ordem da fase paramagnética para a fase absorvente. No espaço de parâmetros as três linhas de transição se encontram no diagrama de fases em um ponto onde o modelo se comporta como o modelo do votante. / In this work we studied a stochastic model with ising symmetry and two simmetric absorbing configurations in a three-dimensional cubic lattice and in two dimensions using a triangular lattice. The study took into account simple mean-field approximations and Monte Carlo simulations. The results showed that the model has a second-order transition from a paramagnetic phase to a ferromagnetic phase and second-order transition from ferromagnetic phase to the absorbing one. A first-order phase transition from the paramagnetic phase to the absorbing phase is observed too. In the phase diagram the two second-order transition lines aproaches to the point where the model behaves like the voter model.

Au cours du traitement d'une tumeur avec des ions légers, la position du pic de Bragg doit être connue avec précision. Une fraction importante des ions incidents subissent des collisions nucléaires avec les noyaux cibles conduisant à l'émission de particules promptes qui peuvent être porteuses d'informations sur le parcours des ions. Ce travail, qui se concentre sur les gamma prompts, montre que le rendement en profondeur de ces émissions est fortement corrélé au parcours des ions et que les taux de comptage mesurés permettent d'envisager un système d'imagerie réaliste, fonctionnant en temps réel. Dans ce but, nous avons réalisé des expériences au GANIL et au GSI avec un détecteur collimaté placé perpendiculairement à l'axe du faisceau et la technique du temps de vol a été utilisée pour réduire le bruit de fond induit par les neutrons et les particules chargées. Des simulations Geant4 ont été réalisées pour concevoir le dispositif expérimental et interpréter les données. Un accord qualitatif entre les simulations et l'expérience est observé pour la quantité d'énergie déposée dans le détecteur et pour la forme du spectre de temps de vol. Cependant, des divergences apparaissent pour le rendement de gamma prompts et la distribution en profondeur des gamma détectés. Ces divergences sont discutées, principalement en termes de modèles de physique nucléaire qui doivent être améliorés. Après avoir sélectionné les modèles physiques offrant les simulations les plus en accord avec les mesures, des études concernant les lieux d'émissions des gamma prompts et l'influence de la diffusion dans la cible ont été réalisés afin de déterminer l'impact sur la corrélation avec le parcours des ions / During the treatment of a tumor with light ions, the Bragg peak location must be accurately known. A significant fraction of the incident ions undergo nuclear collisions with the target nuclei leading to the prompt emission of particles which may carry information on the ion path. This work, which focuses on prompt gamma, shows that the depth profile of these emissions is highly correlated to the ions path and the measured counting rates allow to consider a realistic imaging system, operating in real time. For that purpose, we performed experiments at GANIL and at GSI with a collimated detector placed perpendicular to the beam axis and the time of flight technique was used in order to reduce the noise induced by neutrons and charged particles. Geant4 simulations were performed for the experimental design and data interpretation. A qualitative agreement between simulations and experiment is observed for the amount of energy deposited in the detector and the shape of the time of flight spectrum. However, discrepancies appear for the prompt gamma yield and the depth distribution of gamma detected. These discrepancies are discussed, mainly in terms of nuclear physics models that must be improved. After selecting the physical models which lead to the best agreement between simulations and measurements, studies on the location of prompt gamma emission and on the influence of diffusion in the target were performed to determine the impact on the correlation with the ion path