Le travail de thèse porte sur de nouveaux compléments et améliorations pour la théorie de la péridynamique concernant la modélisation numérique de structures minces telles que les poutres et les plaques, les composites isotropes et multicouches soumis à un chargement dynamique. Nos développements ont principalement porté sur l'exploration des possibilités offertes par la méthode péridynamique, largement appliquée dans divers domaines de l'ingénierie où des discontinuités fortes ou faibles peuvent se produire, telles que des fissures. La procédure de généralisation de la méthode Peridynamics pour la modélisation des structures de poutres de Timoshenko et des structures de plaques de Reissner-Mindlin avec une large plage de rapport épaisseur sur longueur allant de structures épaisses à très minces est indiquée. Et un impact avec une faible vitesse simplifié basé sur le modèle péridynamique développé pour la poutre de Timoshenko et la plaque de Reissner-Mindlin a été proposé en utilisant une procédure de contact spécifique pour l'estimation « naturelle » de la charge d'impact. L’originalité de la méthode actuelle réside dans l’introduction avec deux techniques permettant de réduire le problème de blocage par cisaillement qui se pose dans les structures à poutres et à plaques minces, à savoir la méthode d’intégration réduite (ou sélective) et la formulation mixte. Le modèle péridynamique résultant pour les structures de poutre de Timoshenko et les structures de plaque de Reissner-Mindlin est efficace et ne souffre d'aucun phénomène de verrouillage par cisaillement. En outre, la procédure de généralisation de la méthode péridynamique pour la modélisation de structures composites minces renforcées par des fibres est introduite. L’approche péridynamique pour la modélisation d’une couche est d’abord validée en quasi-statique, ce qui inclut des problèmes de prévision de la propagation de fissures soumis à des conditions de chargement mécaniques. La méthode péridynamique a ensuite été étendue à l’analyse de structures composites minces renforcées par des fibres utilisant la théorie fondamentale d’une couche. Enfin, plusieurs applications impliquant des structures composites minces renforcées par des fibres et des résultats numériques ont été validées par comparaison à la solution FEM obtenue à l'aide d'un logiciel commercial ou à des solutions de référence de la littérature. Dans toutes les applications, Péridynamics montre que les résultats correspondent parfaitement aux solutions de référence, ce qui prouve son potentiel d’efficacité, en particulier pour la simulation de chemins de fissures dans les structures isotropes et composites. / This thesis introduces some new complements and improvments for the Bond-Based Peridynamics theory concerning the numerical modeling of thin structures such as beams and plates, isotropic and multilayer composites subjected to dynamic loading. Our developments have been focused mainly on exploring the possibilities offered by the Peridynamic method, which has been widely applied in various engineering domains where strong or weak discontinuities may occur such as cracks or heterogeneous media. The generalization procedure of the Peridynamics method for the modeling of Timoshenko beam structures and Reissner-Mindlin plate structures respectively with a wide range of thickness to length ratio starting from thick structures to very thin structures is given. And A simplified low velocity impact based on the developed Peridynamic model for Timoshenko beam and ReissnerMindlin plate has been proposed by using a specific contact procedure for the estimation of the impact load. The originality of the present method was the introduction for the first time of two techniques for the alleviation of the shear locking problem which arises in thin beam and plate structures, namely the reduced (or selective) integration method and mixed formulation. The resulting Peridynamic model for Timoshenko beam structures and Reissner-Mindlin plate structures is efficient and does not suffer from any shear locking phenomenon. Besides, the generalization procedure of Peridynamic method for the modeling of fiber-reinforced thin composite structures is introduced. The Peridynamic approach for the modeling of a lamina is firstly validated in the quasi-statics including a crack propagation prediction problems subjected to mechanical loading conditions and then the Peridynamic method was further extended to analyze fiber-reinforced thin composite structures using the fundamental lamina theory. Finally, several applications involving fiber-reinforced thin composite structures and numerical results were validated by comparison to the FEM solution obtained using commercial software or to reference solutions from the literature. In all applications, the Peridynamics shows that results are matching perfectly the reference solutions, which proves its efficiency potentiality especially for crack paths simulation in isotropic and composite structures.

[ES] Esta tesis trata sobre la integración numérica de sistemas hamiltonianos con potenciales explícitamente dependientes del tiempo. Los problemas de este tipo son comunes en la física matemática, porque provienen de la mecánica cuántica, clásica y celestial. La meta de la tesis es construir integradores para unos problemas relevantes no autónomos: la ecuación de Schrödinger, que es el fundamento de la mecánica cuántica; las ecuaciones de Hill y de onda, que describen sistemas oscilatorios; el problema de Kepler con la masa variante en el tiempo. El Capítulo 1 describe la motivación y los objetivos de la obra en el contexto histórico de la integración numérica. En el Capítulo 2 se introducen los conceptos esenciales y unas herramientas fundamentales utilizadas a lo largo de la tesis. El diseño de los integradores propuestos se basa en los métodos de composición y escisión y en el desarrollo de Magnus. En el Capítulo 3 se describe el primero. Su idea principal consta de una recombinación de unos integradores sencillos para obtener la solución del problema. El concepto importante de las condiciones de orden se describe en ese capítulo. En el Capítulo 4 se hace un resumen de las álgebras de Lie y del desarrollo de Magnus que son las herramientas algebraicas que permiten expresar la solución de ecuaciones diferenciales dependientes del tiempo. La ecuación lineal de Schrödinger con potencial dependiente del tiempo está examinada en el Capítulo 5. Dado su estructura particular, nuevos métodos casi sin conmutadores, basados en el desarrollo de Magnus, son construidos. Su eficiencia es demostrada en unos experimentos numéricos con el modelo de Walker-Preston de una molécula dentro de un campo electromagnético. En el Capítulo 6, se diseñan los métodos de Magnus-escisión para las ecuaciones de onda y de Hill. Su eficiencia está demostrada en los experimentos numéricos con varios sistemas oscilatorios: con la ecuación de Mathieu, la ec. de Hill matricial, las ecuaciones de onda y de Klein-Gordon-Fock. El Capítulo 7 explica cómo el enfoque algebraico y el desarrollo de Magnus pueden generalizarse a los problemas no lineales. El ejemplo utilizado es el problema de Kepler con masa decreciente. El Capítulo 8 concluye la tesis, reseña los resultados y traza las posibles direcciones de la investigación futura. / [CAT] Aquesta tesi tracta de la integració numèrica de sistemes hamiltonians amb potencials explícitament dependents del temps. Els problemes d'aquest tipus són comuns en la física matemàtica, perquè provenen de la mecànica quàntica, clàssica i celest. L'objectiu de la tesi és construir integradors per a uns problemes rellevants no autònoms: l'equació de Schrödinger, que és el fonament de la mecànica quàntica; les equacions de Hill i d'ona, que descriuen sistemes oscil·latoris; el problema de Kepler amb la massa variant en el temps. El Capítol 1 descriu la motivació i els objectius de l'obra en el context històric de la integració numèrica. En Capítol 2 s'introdueixen els conceptes essencials i unes ferramentes fonamentals utilitzades al llarg de la tesi. El disseny dels integradors proposats es basa en els mètodes de composició i escissió i en el desenvolupament de Magnus. En el Capítol 3, es descriu el primer. La seua idea principal consta d'una recombinació d'uns integradors senzills per a obtenir la solució del problema. El concepte important de les condicions d'orde es descriu en eixe capítol. El Capítol 4 fa un resum de les àlgebres de Lie i del desenvolupament de Magnus que són les ferramentes algebraiques que permeten expressar la solució d'equacions diferencials dependents del temps. L'equació lineal de Schrödinger amb potencial dependent del temps està examinada en el Capítol 5. Donat la seua estructura particular, nous mètodes quasi sense commutadors, basats en el desenvolupament de Magnus, són construïts. La seua eficiència és demostrada en uns experiments numèrics amb el model de Walker-Preston d'una molècula dins d'un camp electromagnètic. En el Capítol 6 es dissenyen els mètodes de Magnus-escissió per a les equacions d'onda i de Hill. El seu rendiment està demostrat en els experiments numèrics amb diversos sistemes oscil·latoris: amb l'equació de Mathieu, l'ec. de Hill matricial, les equacions d'onda i de Klein-Gordon-Fock. El Capítol 7 explica com l'enfocament algebraic i el desenvolupament de Magnus poden generalitzar-se als problemes no lineals. L'exemple utilitzat és el problema de Kepler amb massa decreixent. El Capítol 8 conclou la tesi, ressenya els resultats i traça les possibles direccions de la investigació futura. / [EN] The present thesis addresses the numerical integration of Hamiltonian systems with explicitly time-dependent potentials. These problems are common in mathematical physics because they come from quantum, classical and celestial mechanics. The goal of the thesis is to construct integrators for several import ant non-autonomous problems: the Schrödinger equation, which is the cornerstone of quantum mechanics; the Hill and the wave equations, that describe oscillating systems; the Kepler problem with time-variant mass. Chapter 1 describes the motivation and the aims of the work in the historical context of numerical integration. In Chapter 2 essential concepts and some fundamental tools used throughout the thesis are introduced. The design of the proposed integrators is based on the composition and splitting methods and the Magnus expansion. In Chapter 3, the former is described. Their main idea is to recombine some simpler integrators to obtain the solution. The salient concept of order conditions is described in that chapter. Chapter 4 summarises Lie algebras and the Magnus expansion ¿ algebraic tools that help to express the solution of time-dependent differential equations. The linear Schrödinger equation with time-dependent potential is considered in Chapter 5. Given its particular structure, new, Magnus-based quasi-commutator-free integrators are build. Their efficiency is shown in numerical experiments with the Walker-Preston model of a molecule in an electromagnetic field. In Chapter 6, Magnus-splitting methods for the wave and the Hill equations are designed. Their performance is demonstrated in numerical experiments with various oscillatory systems: the Mathieu equation, the matrix Hill eq., the wave and the Klein-Gordon-Fock eq. Chapter 7 shows how the algebraic approach and the Magnus expansion can be generalised to non-linear problems. The example used is the Kepler problem with decreasing mass. The thesis is concluded by Chapter 8, in which the results are reviewed and possible directions of future work are outlined. / Kopylov, N. (2019). Magnus-based geometric integrators for dynamical systems with time-dependent potentials [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/118798 / TESIS

This thesis deals with the concepts of numerical integrator using floating point arithmetic for solving partial differential equations. The integrator uses Euler method and Taylor series. Thesis shows parallel and serial approach to computing with exponents and significands. There is also a comparison between modern parallel systems and the proposed concepts.

The application of the finite integral of multiple variable functions is penetrating into more and more industries and science disciplines. The demands placed on solutions to these problems (such as high accuracy or high speed) are often quite contradictory. Therefore, it is not always possible to apply analytical approaches to these problems; numerical methods provide a suitable alternative. However, the ever-growing complexity of these problems places too high a demand on many of these numerical methods, and so neither of these methods are useful for solving such problems. The goal of this thesis is to design and implement a new numerical method that provides highly accurate and very fast computation of finite integrals of multiple variable functions. This new method combines pre-existing approaches in the field of numerical mathematics.

Gitter QCD ist ein erfolgreiches Instrument zur nicht-perturbativen Berechnung von QCD Observablen. Die hierfür notwendige Auswertung des QCD Pfadintegrals besteht aus zwei Teilen: Zuerst werden Stützstellen generiert, an denen danach das Pfadintegral ausgewertet wird. In der Regel werden für den ersten Teil Markov-chain Monte Carlo (MCMC) Methoden verwendet, die für die meisten Anwendungen sehr gute Ergebnisse liefern, aber auch Probleme wie eine langsame Fehlerskalierung und das numerische Vorzeichenproblem bergen. Der zweite Teil beinhaltet die Berechnung von Quark zusammenhängenden und unzusammenhängenden Diagrammen. Letztere tragen maßgeblich zu physikalischen Observablen bei, jedoch leidet deren Berechnung an großen Fehlerabschätzungen. In dieser Arbeit werden Methoden präsentiert, um die beschriebenen Schwierigkeiten in beiden Auswertungsteilen des QCD Pfadintegrals anzugehen und somit Observablen effizienter beziehungsweise genauer abschätzen zu können. Für die Berechnung der unzusammenhängenden Diagramme haben wir die Methode der exakten Eigenmodenrekonstruktion mit Deflation getestet und konnten eine 5.5 fache Verbesserung der Laufzeit erreichen. Um die Probleme von MCMC Methoden zu adressieren haben wir die rekursive numerische Integration zur Vereinfachung von Integralauswertungen getestet. Wir haben diese Methode, kominiert mit einer Gauß-Quadraturregel, auf den eindimensionalen quantenmechanischen Rotor angewandt und konnten exponentiell skalierende Fehlerabschätzungen erreichen. Der nächste Schritt ist eine Verallgemeinerung zu höheren Raumzeit Dimensionen. Außerdem haben wir symmetrisierte Quadraturregeln entwickelt, um das Vorzeichenproblem zu umgehen. Wir haben diese Regeln auf die eindimensionale QCD mit chemischem Potential angewandt und konnten zeigen, dass sie das Vorzeichenproblem beseitigen und sehr effizient auf Modelle mit einer Variablen angewendet werden können. Zukünftig kann die Effizienz für mehr Variablen verbessert werden. / Lattice QCD is a very successful tool to compute QCD observables non-perturbatively from first principles. The therefore needed evaluation of the QCD path integral consists of two parts: first, sampling points are generated at which second, the path integral is evaluated. The first part is typically achieved by Markov-chain Monte Carlo (MCMC) methods which work very well for most applications but also have some issues as their slow error scaling and the numerical sign-problem. The second part includes the computation of quark connected and disconnected diagrams. Improvements of the signal-to-noise ratio have to be found since the disconnected diagrams, though their estimation being very noisy, contribute significantly to physical observables. Methods are proposed to overcome the aforementioned difficulties in both parts of the evaluation of the lattice QCD path integral and therefore to estimate observables more efficiently and more accurately. For the computation of quark disconnected diagrams we tested the exact eigenmode reconstruction with deflation method and found that this method resulted in a 5.5-fold reduction of runtime. To address the difficulties of MCMC methods, we tested the recursive numerical integration method, which simplifies the evaluation of the integral. We applied the method in combination with a Gauss quadrature rule to the one-dimensional quantum-mechanical rotor and found that we can compute error estimates that scale exponentially to the correct result. A generalization to higher space-time dimensions can be done in the future. Additionally, we developed the symmetrized quadrature rules to address the sign-problem. We applied them to the one-dimensional QCD with a chemical potential and found that this method is capable of overcoming the sign-problem completely and is very efficient for models with one variable. Improvements of the efficiency for multi-variable scenarios can be made in the future.

Für die Simulation technischer Bauteile mit Hilfe der Finite-Elemente-Methode (FEM) werden tensorielle Stoffgesetze benötigt, die zu einem beliebigen dreidimensionalen Verzerrungszustand und gegebenenfalls der Belastungsvorgeschichte und -geschwindigkeit des Materials die zugehörige Spannungsantwort liefern. Die Entwicklung derart komplexer Materialmodelle verläuft oftmals über Zwischenstufen, die zunächst nur Vorhersagen für den einachsigen Zug-/Druckversuch erlauben. Zur automatischen Verallgemeinerung solcher eindimensionaler Materialbeschreibungen zu vollständig dreidimensionalen Stoffgesetzen für die Finite-Elemente-Methode wird im Rahmen dieser Arbeit das Konzept repräsentativer Raumrichtungen vorgeschlagen, welches auf der Integration einachsiger Spannungszustände über eine diskrete Anzahl gleichmäßig verteilter (repräsentativer) Raumrichtungen basiert. Zur Untersuchung der grundlegenden Eigenschaften des Algorithmus wurden verschiedene inelastische tensorielle Beispielstoffgesetze herangezogen, deren eindimensionale Formulierung als Eingangsmodell für die repräsentativen Raumrichtungen dient. Hierbei zeigt sich, dass die wesentlichen Materialeigenschaften des jeweiligen uniaxialen Eingangsmodells bei der Verallgemeinerung vollständig erhalten bleiben. Weiterhin werden einige wichtige Effekte vom Konzept automatisch generiert, wie z. B. die anisotrope Entfestigung technischer Gummiwerkstoffe oder die formative Verfestigung metallischer Werkstoffe, was eine realitätsnahe Simulation dieser Materialklassen ohne zusätzlichen Arbeitsaufwand erlaubt. Das Konzept wurde zusätzlich auf Stoffgesetze angewendet, die ausschließlich in Form einer eindimensionalen Materialbeschreibung vorliegen und somit konkrete Anwendungsfällle darstellen. Darüber hinaus wurden für einige ausgewählte Stoffgesetze in repräsentativen Raumrichtungen Vergleiche mit Ergebnissen aus experimentellen Versuchen vorgenommen, wobei sich stets eine gute Übereinstimmung zwischen Experiment und Simulation ergibt. Das Konzept repräsentativer Raumrichtungen wurde in die zwei kommerziellen Finite-Elemente-Programme MSC.Marc und ABAQUS implementiert. Hiermit können Simulationen inhomogener Verzerrungs- und Spannungsverteilungen durchgeführt werden, obwohl das zugrunde liegende Stoffgesetz lediglich einachsige Spannungszustände beschreibt. In diesem Zusammenhang werden verschiedene Methoden vorgestellt, mit deren Hilfe die Effizienz einer FEM-Simulation erheblich gesteigert werden kann. Dies betrifft zum einen die Generierung einer gleichmäßigen Verteilung von repräsentativen Raumrichtungen mit Hilfe eines numerischen Algorithmus zur Simulation sich abstoßender elektrischer Punktladungen auf der Kugeloberfläche. Zum anderen besteht die Möglichkeit, die einzelnen Sätze von repräsentativen Raumrichungen in den Gaußpunkten eines finiten Elementes unterschiedlich zueinander auszurichten, was bei gleichbleibendem Rechenaufwand eine beträchtliche Erhöhung der Rechengenauigkeit erlaubt. / The simulation of technical components using the finite element method (FEM) requires tensorial constitutive models which describe the complete relation between a given three-dimensional state of strain (in some cases also the loading history and strain rate) and the corresponding state of stress. The development of such complex material models often leads to an intermediate stage that enables the prediction of uniaxial tension and compression only. The automatic generalization of those one-dimensional material descriptions to complete three-dimensional constitutive models for the finite element method can be accomplished by using the concept of representative directions which is based on the integration of uniaxial stresses over a discrete number of uniformly distributed (representative) directions in space. In order to investigate the fundamental characteristics of the algorithm several inelastic tensorial constitutive models were used, whose one-dimensional formulation serves as the input model for the use within the representative directions. In this context it becomes evident that the essential material properties of the respective uniaxial input model are completely preserved during the process of generalization. Furthermore, some important effects are produced automatically by the concept such as the anisotropic stress softening of technical rubber materials or the distortional hardening of metallic materials, which enables a realistic simulation of those material classes without spending additional effort. The concept was also applied to material models that are available in form of a one-dimensional material description only, so that these can be regarded as concrete applications. In addition, some of the material models in representative directions were compared to experimental data, whereas a good agreement between measurement and simulation can be noticed. The concept of representative directions has been implemented into the commercial finite element programs MSC.Marc and ABAQUS. This enables simulations of inhomogeneous strain and stress distributions even though the underlying material model describes uniaxial loading processes only. In this context, several methods are introduced which can be applied to increase the efficiency of a finite element simulation to a great extent. On the one hand this affects the generation of a uniform distribution of representative directions using a numerical algorithm simulating the repulsion of electric charges on the surface of a sphere. On the other hand, it is possible to adjust the sets of representative directions at the integration points of a finite element differently, which leads to an increasing computational accuracy at constant computational effort.

This PhD thesis consists of a summary and seven papers, where various applications of auto-validated computations are studied. In the first paper we describe a rigorous method to determine unknown parameters in a system of ordinary differential equations from measured data with known bounds on the noise of the measurements. Papers II, III, IV, and V are concerned with Abelian integrals. In Paper II, we construct an auto-validated algorithm to compute Abelian integrals. In Paper III we investigate, via an example, how one can use this algorithm to determine the possible configurations of limit cycles that can bifurcate from a given Hamiltonian vector field. In Paper IV we construct an example of a perturbation of degree five of a Hamiltonian vector field of degree five, with 27 limit cycles, and in Paper V we construct an example of a perturbation of degree seven of a Hamiltonian vector field of degree seven, with 53 limit cycles. These are new lower bounds for the maximum number of limit cycles that can bifurcate from a Hamiltonian vector field for those degrees. In Papers VI, and VII, we study a certain kind of normal form for real hyperbolic saddles, which is numerically robust. In Paper VI we describe an algorithm how to automatically compute these normal forms in the planar case. In Paper VII we use the properties of the normal form to compute local invariant manifolds in a neighbourhood of the saddle.

Der Strahlungstransport stellt eine von drei Arten des Wärmetransports zwischen Gebieten unterschiedlicher Temperatur dar. Eine der einfachsten Formen bildet der Strahlungstransport im Vakuum, ein Vorgang, der im kosmischen Umfeld, beispielsweise bei der Energieübertragung von einem Stern auf seine Planeten, beobachtbar ist. Hierbei ist es hinreichend, sich auf die Betrachtung von Oberflächen zu beschränken. Strahlungstransport kann jedoch auch in semitransparenten Medien, wie biologischem Gewebe oder Glas, beobachtet werden. Das Medium, in dem der Strahlungstransport erfolgt, wirkt sich durch Vorgänge wie Absorption, Emission, Reflexion oder Streuung auf den Strahlungstransport aus. Für die Modellierung des Strahlungstransports in einem solchen Umfeld können verschiedene Modelle, darunter das Strahlenmodell, genutzt werden. Dieses Modell beschreibt den Wärmetransport anhand einer skalaren Größe, die Strahlungsintensität genannt wird. Betrachtet wird die Strahlungsintensität in diesem Modell entlang eines Strahls in eine vorgegebene Richtung. Die mathematische Darstellung des Strahlenmodells des Strahlungstransports in partizipierenden Medien führt auf eine richtungsabhängige Integro-Differentialgleichung. Ist die Richtungsabhängigkeit nicht von Interesse, so kann der Übergang zu einer winkelintegrierten Form erfolgen. Dieser Übergang führt schließlich auf ein System schwach singulärer fredholmscher Integralgleichungen zweiter Art. Dieses charakterisiert nun nicht mehr die erwähnte Strahlungsintensität, sondern beschreibt die sogenannte Einstrahlung sowie den Strahlungsfluss. Das System singulärer Integralgleichungen kann mittels eines Galerkin-Ansatzes numerisch gelöst werden. Geht man von einer hinreichenden Glattheit des Randes aus, kann die Kompaktheit des Operators der Integralgleichungen gezeigt werden. Dies wiederum erlaubt Rückschlüsse auf die Existenz und Eindeutigkeit einer Lösung. Ein Augenmerk bei der Ermittlung der Galerkin-Näherung ist auf die Bestimmung der singulären Integrale der Galerkin-Diskretisierung zu richten. Für die Bestimmung multidimensionaler, singulärer Integrale stellt die Arbeit das Verfahren der hierarchischen Integration vor. Basierend auf einer Zerlegung des Integrationsgebietes, erfolgt die Beschreibung singulärer Integrale durch ein Gleichungssystem, dessen rechte Seite nur von regulären Integralen abhängig ist. Können diese regulären Integrale sowie die Lösung des Gleichungssystems exakt bestimmt werden, so sind auch die singulären Integrale exakt bestimmt. Bei einer numerischen Bestimmung der regulären Integrale ist die Fehlerordnung ausschlaggebend für den Fehler der singulären Integrale. Als Integrationsgebiete werden Hyperwürfel beliebiger Dimension sowie Simplizes bis einschließlich Dimension 3 als Integrationsgebiete betrachtet. Als Voraussetzungen an den Kern des Doppelintegrals sind nur die Eigenschaften der Translationsinvarianz sowie der Homogenität zu richten. Kann ein nicht translationsinvarianter oder nicht homogener Kern eines Integrals in Summanden zerlegt werden, die selbst translationsinvariant und homogen sind, ist auch die Bestimmung solcher Integrale möglich. Darüber hinaus stellt die Arbeit Verbindungen zu dem Begriff des Hadamard partie finie her. Auf diese Weise lässt sich das Verfahren der hierarchischen Integration für beliebige Dimensionen und beliebige Singularitätsordnungen anwenden. Die Strahlungstransportgleichung ist im Allgemeinen mittels eines Galerkin-Ansatzes lösbar, führt jedoch auf eine voll besetzte Systemmatrix. Numerische Beispiele beleuchten daher Methoden der Matrixkompression mittels hierarchischer Matrizen sowie der direkten Erzeugung schwach besetzter Matrizen über regulären Gittern und Gittern mit hängenden Knoten und skizziert Ansätze zur Parallelisierung auf entsprechenden Computersystemen.

The kinetics of drugs photodegradation have traditionally been treated using thermal kinetic analysis methods consisting most commonly in zero and first order kinetics. These treatment strategies were shown to lack specificity and present a number of limitations when applied to photoreactions kinetics. Nevertheless, these methods have widely been used due to a lack of integrated rate-laws for the majority of photoreactions types, in turn, due to the presence of a variable time-dependent factor in most photoreactions rate-laws that prevents their mathematical integration. To address these limitations, a new methodology for the development and validation of semi-empirical integrated rate-laws that faithfully describe photoreactions kinetics and photoreactions simulated cases generated by numerical integration methods (NIMs), is hereby presented. Using this methodology, a new kinetic order was ascribed to photoreactions namely the Φ-order kinetics. Semi-empirical integrated rate-laws were, thus, developed for three photoreaction types namely, unimolecular, AB(1Φ), photoreversible ,AB(2Φ), and consecutive, AB4(4Φ), photoreactions. The proposed models were further tested experimentally on drugs following these photodegradation mechanisms using; nifedipine and nisoldipine for unimolecular photoreactions; montelukast and fluvoxamine for photoreversible reactions; and riboflavin for consecutive photoreactions. The developed models not only accurately described the photoreaction kinetics of these drugs but also allowed the determination of all the kinetic parameters that characterise them. Furthermore, the above studied drugs were shown to act as precise and simple actinometers when analytically treated with the Φ-order kinetic methods, hereby presented. A universal standard method for the precise and worldwide reproducible study of drugs stability and compounds photoreactions, based on monochromatic irradiation and Φ-kinetics data analysis, is also detailed and adopted throughout the thesis. Finally, two new kinetic parameters namely, the pseudo-rate-constant and pseudo-initial velocity have been identified and shown to be more reliable and accurate in the description and universal comparison of photoreactions kinetics.

Orientador: Maria Cristina Dias Tavares / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e Computação / Made available in DSpace on 2018-08-07T13:37:43Z (GMT). No. of bitstreams: 1 Paz_MarcosdeAraujo_D.pdf: 6127820 bytes, checksum: 21f6dc485bcf097846769f086c2362ad (MD5) Previous issue date: 2005 / Resumo: O objetivo deste trabalho é contribuir para o estudo de linhas de transmissão e o desenvolvimento de modelos de linhas voltados a programas do tipo EMTP. O trabalho desenvolvido consiste em aproveitar as propriedades e semelhanças dos modos não homopolares de propagação das linhas de transmissão, sejam elas idealmente transpostas ou que apresentem plano de simetria vertical. Por meio de manipulações da matriz de transformação fase-modo-fase através do uso de vetores ortogonais e elementos complexos no domínio modal, pode-se reduzir a dimensão dos vetores de estado. O procedimento proposto é direcionado às modelagens que utilizam transformação fase-modo-fase na determinação das variáveis de estado (corrente e tensão) nas linhas de transmissão. Os recursos usados levam a uma redução tanto em número de operações por iteração quanto em alocação de estados, otimizando o procedimento de cá1culo e aumentando a velocidade de simulação no processamento. A economia computacional que o procedimento proposto proporciona pode ser usada em ferramentas de simulação em tempo real / Abstract: The research objective is to contribute with the study of the transmission line and the development of line models in EMTP-type programs. The proposed model takes advantage of some similarities between non homopolar modes, applying orthogonal vectors and complex number theory to manipulate the state vectors in modal domain reducing the state vectors dimension. The proposed procedure is aimed to the modelling that use phase-mode-phase transformation in the determination of state variables (current and voltage) in transmission lines. The resources used lead to a reduction both in the number of operations per iteration and state allocation, optimizing the procedure and increasing the processing simulation speed. The computational economy, which the proposed procedure provide, can be applied to real-time simulation tools / Doutorado / Energia Eletrica / Doutor em Engenharia Elétrica