• JoomlaWorks Simple Image Rotator
  • JoomlaWorks Simple Image Rotator
  • JoomlaWorks Simple Image Rotator
  • JoomlaWorks Simple Image Rotator
  • JoomlaWorks Simple Image Rotator
  • JoomlaWorks Simple Image Rotator
  • JoomlaWorks Simple Image Rotator
  • JoomlaWorks Simple Image Rotator
  • JoomlaWorks Simple Image Rotator
  • JoomlaWorks Simple Image Rotator
 
  Bookmark and Share
 
 
Mémoire de Maîtrise
DOI
https://doi.org/10.11606/D.45.2016.tde-25082016-221140
Document
Auteur
Nom complet
Guilherme Galina Loch
Adresse Mail
Unité de l'USP
Domain de Connaissance
Date de Soutenance
Editeur
São Paulo, 2016
Directeur
Jury
Bevilacqua, Joyce da Silva (Président)
Kuhl, Nelson Mugayar
Rodrigues Junior, Orlando
Titre en portugais
Sistemas rígidos associados a cadeias de decaimento radioativo
Mots-clés en portugais
Cadeias de decaimento radioativo
Equações diferenciais ordinárias rígidas
Método de Radau II
Método de Rosenbrock
Resumé en portugais
Os progressos computacionais nas últimas décadas e a teoria matemática cada vez mais sólida têm possibilitado a resolução de problemas de alta complexidade, permitindo uma modelagem cada vez mais detalhada da realidade. Tal verdade aplica-se inclusive para os sistemas rígidos de Equações Diferencias Ordinárias (EDOs): existem métodos numéricos altamente performáticos para este tipo de problema, que permitem uma grande variação no tamanho do passo de integração sem impactar na sua convergência. Este trabalho apresenta um estudo sobre o conceito de rigidez e técnicas numéricas para resolução de problemas rígidos de EDOs. O que nos motivou a estudar tais técnicas foram problemas oriundos da Física Nuclear que envolvem cadeias de decaimento radioativo. Estes problemas podem ser modelados por uma cadeia fechada de compartimentos que se traduz em um sistema de EDOs. Os elementos destas cadeias podem possuir constantes de decaimento com ordens de grandeza muito distintas, caracterizando a sua rigidez e exigindo cautela na resolução das equações que as modelam. Embora seja possível determinar a solução analítica para estes problemas, o uso de métodos numéricos facilita a obtenção da solução quando consideramos sistemas com um número elevado de equações. Além disso, soluções numéricas permitem adaptações na modelagem ou em ajustes de dados com mais facilidade. Métodos implícitos são indicados para a resolução deste tipo de problema, pois possuem uma região de estabilidade ilimitada. Neste trabalho, implementamos dois métodos numéricos que possuem esta característica: o método de Radau II e o método de Rosenbrock. Estes métodos foram utilizados para obtenção de soluções numéricas robustas para problemas rígidos de decaimento radioativo envolvendo cadeias naturais e artificiais, considerando retiradas de elementos das cadeias durante o processo de decaimento e quando queremos determinar qual era o estado inicial de uma cadeia que está em decaimento. Ambos os métodos foram implementados com estratégias de controle do tamanho do passo de integração e produziram resultados consistentes dentro de uma precisão pré-fixada.
Titre en anglais
Stiff systems associated with radioactive decay chains
Mots-clés en anglais
Radau II method
Radioactive decay chains
Rosenbrock method
Stiff ordinary differential equations
Resumé en anglais
The computational progress in the last decades and the increasingly solid mathematical theory have made possible the resolution of highly complex problems allowing an ever more detailed modelling of reality. This is true even for the systems of stiff Ordinary Differential Equations (ODEs): there are highly performative numerical methods for this kind of problem which allow a wide variation in the size of integration step without impacting on their convergence. This thesis presents a study about the concept of stiffness and numerical techniques to solve stiff problems of ODEs. What motivated us to study these techniques were problems from the Nuclear Physics involving radioactive decay chains. These problems could be modelled by a closed chain of compartments which is translated into a system of ODEs. The elements of these chains could have decay constants with very different orders of magnitude which characterizes the stiffness of the problem and requires caution in solving the model equations. Although it is possible to determine the analytical solution to these problems when we consider systems with a high number of equations, calculate the solution by numerical methods becomes easier. Furthermore, numerical solutions allow adaptations in modelling or data adjustments more easily. Implicit methods are indicated to solve this kind of problem because they have an unlimited region of stability. In this study, we implemented two numerical methods which have this feature: Radau II method and Rosenbrock method. These methods were used to obtain robust numerical solutions for stiff problems of radioactive decay involving natural and artificial chains, considering the removal of elements during the decay process and when we want to determine what was the initial state of a chain which is decaying. Both methods were implemented with control strategies for integration step size providing consistent results within a pre-established accuracy.
 
AVERTISSEMENT - Regarde ce document est soumise à votre acceptation des conditions d'utilisation suivantes:
Ce document est uniquement à des fins privées pour la recherche et l'enseignement. Reproduction à des fins commerciales est interdite. Cette droits couvrent l'ensemble des données sur ce document ainsi que son contenu. Toute utilisation ou de copie de ce document, en totalité ou en partie, doit inclure le nom de l'auteur.
Date de Publication
2016-09-12
 
AVERTISSEMENT: Apprenez ce que sont des œvres dérivées cliquant ici.
Tous droits de la thèse/dissertation appartiennent aux auteurs
CeTI-SC/STI
Bibliothèque Numérique de Thèses et Mémoires de l'USP. Copyright © 2001-2024. Tous droits réservés.