[pt] O presente trabalho tem como finalidade avaliar os diferentes modelos de mistura para o cálculo da combustão de reagentes pré- misturados utilizando a abordagem de Reator Parcialmente Misturado (PaSR). Os modelos de mistura considerados neste trabalho foram os modelos IEM estendido, Langevin e Langevin estendido. Investiga-se aqui o grau de mistura previsto por tais modelos e sua influência sobre as propriedades termoquímicas em um processo de combustão. A primeira parte deste trabalho consiste na apresentação e avaliação destes modelos de mistura, considerando-se um campo escalar inerte em presença de um campo turbulento homogêneo e isotrópico. Uma vez que estes modelos de mistura envolvem formulações do tipo estocástico, sua implementação foi realizada utilizando o método de Monte Carlo, mediante a utilização de esquemas numéricos adequados à resolução de equações diferenciais estocásticas. Assim, estuda-se a evolução da Função Densidade de Probabilidade (PDF) e das principais propriedades do campo escalar para cada modelo implementado. Os resultados obtidos também são comparados com simulação numérica direta e com resultados analáticos disponsáveis. Um ótimo acordo em termos qualitativos e quantitativos é obtido. A segunda parte deste trabalho utiliza estes modelos para o estudo numérico de um PaSR no qual são modelados os processos difusivos e reativos presentes durante a combustão. O PaSR é usado para avaliar a influência dos modelos de mistura nas propriedades termoquímicas da mistura em uma situação de combustão de tipo pré-misturada, que é modelada utilizando-se uma variável de progresso de uma reação. Os resultados obtidos com os diferentes modelos de mistura são comparados para diferentes regimes de funcionamento do PaSR, mostrando que, em situações de mistura rápida e reação intensa, os diferentes modelos apresentam resultados similares. Porém, nos casos de mistura lenta e reação moderada, discrepancias importantes são observadas entre os resultados dos modelos; as quais atingem até 65% para o valor médio da variável de progresso da reação. / [en] The present work evaluates several mixing models for the prediction of premixed combustion in a Partially Stirred Reactor (PaSR). The models considered in this work were the extended IEM, Langevin and extended Langevin models. The degree of mixing and its influence on the termochemical properties in a combustion process are investigated here. The first part of this work consists on the presentation and the assesment of these mixing models in which a single scalar field was considered in presence of a homogeneous and isotropic turbulent field. Since these mixing models involve stochastic terms, their implementation is performed by the Monte Carlo method using numerical schemes which solve the corresponding Stochastic Differential Equations (SDE). The evolution of the Probability Density Function (PDF) and the main properties for a single scalar field are studied for each mixing model. The numerical results are compared with Direct Numerical Simulation and available analytical results. Excellent qualitative and quantitative agreements are obtained. In the second part of this work, mixing models are used for numerical simulation of a PaSR where the diffusive and reactive processes occur. The PaSR is used to assess the mixing model influence on the termochemical properties of the mixture in a premixed combustion process, which is modeled using a reaction progress variable. The results obtained with the different mixing models are compared in several operating regimes of the PaSR, showing that when mixing is fast and reaction is intense, the different models lead to similar results. However, when mixing is slow and reaction is weak, important discrepancies are observed between the model results, which reach 65%, as far as the averaged reaction progress variable is concerned

Este trabalho apresenta um estudo qualitativo das equações diferenciais nãolineares que descrevem o sincronismo de fase nos PLLs de 3ª ordem que compõem redes OWMS de topologia mista, Estrela Simples e Cadeia Simples. O objetivo é determinar, através da Teoria de Bifurcações, os valores ou relações entre os parâmetros constitutivos da rede que permitam a existência e a estabilidade do estado síncrono, quando são aplicadas, no oscilador mestre, duas funções de excitação muito comuns na prática: o degrau e a rampa de fase. Na determinação da estabilidade dos pontos de equilíbrio, sob o ponto de vista de Lyapunov, a existência de pontos de equilíbrio não-hiperbólicos não permite uma aproximação linear e, nesses casos, é aplicado o Teorema da Variedade Central. Essa técnica de simplificação de sistemas dinâmicos permite fazer uma aproximação homeomórfica em torno desses pontos, preservando a orientação no espaço de fases e possibilitando determinar localmente suas estabilidades. / This work presents a qualitative study of the non-linear differential equations that describe the synchronous state in 3rd order PLLs that compose One-way masterslave time distribution networks with Single Star and Single Chain topologies. Using bifurcation theory, the dynamical behavior of third-order phase-locked loops employed to extract the syncronous state in each node is analyzed depending on constitutive node parameters when two usual inputs, the step and the ramp phase pertubations, are supposed to appear in the master node. When parameter combinations result in non hyperbolic synchronous states, from Lyapunov point of view, the linear approximation does not provide any information about the local behavior of the system. In this case, the center manifold theorem permits the construction of an equivalent vector field representing the asymptotic behavior of the original system in the neighborhood of these points. Thus, the local stability can be determined.

A tese apresenta um método para a solução dos problemas do valor inicial (PVIs) com margens de erro comparáveis às de métodos numéricos consagrados (MN), tanto para a função quanto para suas derivadas. O método é aplicável a equações diferenciais (EDs) lineares ou não, sendo o ferramental desenvolvido até a quarta ordem, que pode ser expandido para ordens superiores. A solução é uma expressão polinomial de alto grau com coeficientes expressos pela razão entre dois inteiros. O método se mostra eficaz mesmo em alguns casos em que os MN não conseguiram dar a partida. As resoluções são obtidas considerando que o espaço de soluções é um espaço de Hilbert, equipado com a base completa dos polinômios de Legendre. Em decorrência do método aqui desenvolvido, os majorantes de erros para a função e derivadas são determinados analiticamente por um cálculo matricial também deduzido nesta tese. Paralelamente a toda fundamentação analítica, foi desenvolvido o software SAM, que automatiza todas as tarefas na busca de soluções dos PVIs. A tese propõe e verifica a validade de um novo critério de erro no qual pesam tanto os erros locais quanto os erros globais, simultaneamente. Como subprodutos dos resultados já descritos, igualmente integrados ao SAM, obtiveram-se também: (1) Um critério objetivo para analisar a qualidade de um MN, sem necessidade do conhecimento de seu algoritmo; (2) Uma ferramenta para aproximações polinomiais de alta precisão para funções de quadrado integrável em determinado intervalo limitado, com um majorante de erro; (3) Um ferramental analítico para transposição genérica (linear ou não) dos PVIs até 4ª ordem, nas mudanças de domínio; (4) As matrizes de integração e diferenciação genéricas para todas as bases polinomiais do espaço de Hilbert. / This thesis shows a new method to get polynomial solutions to the initial value problems (IVP), with an error margin comparable to the consecrate numerical methods (NM), for both the function and its derivatives. The method works with differential equations (DEs) linear or not, beeing the developed tolls available until 4th order, whose can be expanded to higher orders. The solution is a polynomial high degree expression with coefficients expressed by the ratio between two integers. The method behaves efficiently even in some cases that NM cannot get started. The resolutions are gotten considering that, the solution space is a Hilbert space, equipped with a complete set basis of Legendre Polynomials. Due the method here developed, the errors majoratives for the function and its derivatives are found analytically by a matrix calculus, also derived in this thesis. Beside all analytical foundation, a software (SAM) was developed to automate the whole process, joining all the tasks involved in the search for solutions to the IVP. This thesis proposes, verifies and validates a new error criterion, which takes in account simultaneously the local and global errors. As sub-products of the results described before, also integrated to the SAM, the following achievements should be highlighted: (1) An objective criterion to analyze the quality of any NM, despite of the knowledge of its algorithm; (2) A tool for a polynomial approximation, of high precision, for functions whose square is integrable in a given limited domain, with an errors majorative; (3) A tool-kit for a generically transpose (linear or not) of the IVPs domain and form, taking into account its derivatives, until the 4th order; (4) The generic matrices for integration and differentiation for all the polynomial basis of the Hilbert space.

Reatores tubulares de polimerização podem apresentar um perfil de velocidade bastante distorcido. Partindo desta observação, um modelo estocástico baseado no modelo de dispersão axial foi proposto para a representação matemática da fluidodinâmica de um reator tubular para produção de poliestireno. A equação diferencial foi obtida inserindo a aleatoriedade no parâmetro de dispersão, resultando na adição de um termo estocástico ao modelo capaz de simular as oscilações observadas experimentalmente. A equação diferencial estocástica foi discretizada e resolvida pelo método Euler-Maruyama de forma satisfatória. Uma função estimadora foi desenvolvida para a obtenção do parâmetro do termo estocástico e o parâmetro do termo determinístico foi calculado pelo método dos mínimos quadrados. Uma análise de convergência foi conduzida para determinar o número de elementos da discretização e o modelo foi validado através da comparação de trajetórias e de intervalos de confiança computacionais com dados experimentais. O resultado obtido foi satisfatório, o que auxilia na compreensão do comportamento fluidodinâmico complexo do reator estudado. / The velocity profile of polymerization tubular reactors may be very distorted. Based on this observation, a stochastic model based on the axial dispersion model was proposed for the mathematical representation of the fluid dynamics of a tubular reactor for polystyrene production. The differential equation was built by inserting randomness in the dipersion coefficient, which added a stochastic term to the model. This term was capable of simulating the experimentally observed fluctuations. The stochastic differential equation was discretized and solved by the Euler-Maruyama method adequately. An estimator function has been developed to calculate the parameter of the stochastic term, while the parameter of the deterministic term was estimated by a least squares method. A convergence analysis was carried out in order to determine the number of elements needed for the time discretization. The model was validated through comparisons of sample paths and computational confidence intervals with experimental data. The result was considered satisfactory, allowing a better understanding of the complex fluid dynamic behaviour of the analised reactor. Key-words: modelling, simulation, stochastic differential equation, polymerization tubular reactor, time residence distribution.

Les systèmes décrits par les équations aux dérivées partielles, appartiennent à la classe des systèmes dynamiques impliquant des fonctions dépendantes de plusieurs variables, comme le temps et l'espace. Déjà fortement répandus pour la modélisation mathématique de phénomènes physiques et environnementaux, ces systèmes ont un rôle croissant dans les domaines de l'automatique. Cette expansion, provoquée par les avancées technologiques au niveau des capteurs facilitant l'acquisition de données et par les nouveaux enjeux environnementaux, incite au développement de nouvelles thématiques de recherche. L'une de ces thématiques, est l'étude des problèmes inverses et plus particulièrement l'identification paramétrique des équations aux dérivées partielles. Tout abord, une description détaillée des différentes classes de systèmes décrits par ces équations est présentée puis les problèmes d'identification qui leur sont associés sont soulevés. L'accent est mis sur l'estimation paramétrique des équations linéaires, homogènes ou non, et sur les équations linéaires à paramètres variant. Un point commun à ces problèmes d'identification réside dans le caractère bruité et échantillonné des mesures de la sortie. Pour ce faire, deux types d'outils principaux ont été élaborés. Certaines techniques de discrétisation spatio-temporelle ont été utilisées pour faire face au caractère échantillonné des données; les méthodes de variable instrumentale, pour traiter le problème lié à la présence de bruit de mesure. Les performances de ces méthodes ont été évaluées selon des protocoles de simulation numérique reproduisant des situations réalistes de phénomènes physique et environnementaux, comme la diffusion de polluant dans une rivière / A large variety of natural, industrial, and environmental systems involves phenomena that are continuous functions not only of time, but also of other independent variables, such as space coordinates. Typical examples are transportation phenomena of mass or energy, such as heat transmission and/or exchange, humidity diffusion or concentration distributions. These systems are intrinsically distributed parameter systems whose description usually requires the introduction of partial differential equations. There is a significant number of phenomena that can be simulated and explained by partial differential equations. Unfortunately all phenomena are not likely to be represented by a single equation. Also, it is necessary to model the largest possible number of behaviors to consider several classes of partial differential equations. The most common are linear equations, but the most representative are non-linear equations. The nonlinear equations can be formulated in many different ways, the interest in nonlinear equations with linear parameters varying is studied. The aim of the thesis is to develop new estimators to identify the systems described by these partial differential equations. These estimators must be adapted with the actual data obtained in experiments. It is therefore necessary to develop estimators that provide convergent estimates when one is in the presence of missing data and are robust to measurement noise. In this thesis, identification methods are proposed for partial differential equation parameter estimation. These methods involve the introduction of estimators based on the instrumental variable technique

Der erste Teil der Arbeit beschreibt eine Methode, Auszahlungen zu bewerten, die einem auf dem Finanzmarkt nicht absicherbaren Risiken ausgesetzt sind. Im zweiten Teil berechnen wir den maximalen Nutzen und optimale Handelsstrategien auf unvollständigen Märkten mit Hilfe von stochastischen Rückwärtsgleichungen. Wir betrachten Händler, deren Einkommen einer externen Risikoquelle ausgesetzt sind. Diese vervollständigen den Markt, indem sie entweder einen Bond schaffen oder gegenseitig Verträge schliessen. Eine andere Moeglichkeit ist eine Anleihe, die von einer Versicherung herausgegeben wird. Die Risikoquellen, die wir in Betracht ziehen, können Versicherungs-, Wetter-oder Klimarisiko sein. Aktienpreise sind exogen gegeben. Wir berechnen Preise für die zusätzlichen Anlagen so dass Angebot und Nachfrage dafür gleich sind. Wir haben partielle Markträumung. Die Präferenzen der Händler sind durch erwarteten Nutzen gegeben. In Kapitel 2 bis Kapitel 4 haben die Händler exponentielle Nutzenfunktionen. Um den Gleichgewichtspreis zu finden, wenden wir stochastische Rückwärtsgleichungen an. In Kapitel 5 beschreiben wir ein Einperiodenmodell mit Nutzenfunktionen, die die Inada-Bedingungen erfüllen. Der zweite Teil dieser Arbeit beschäftigt sich mit dem robusten Nutzenmaximierungsproblem auf einem unvollständigen Finanzmarkt. Entweder das Wahrscheinlichkeitsmass oder die Koeffizienten des Aktienmarktes sind ungewiss. Die Lösung der Rückwärtsgleichung beschreibt die nutzenmaximierende Handelsstrategie und das Wahrscheinlichkeitsmass, das in der Auswertung des robusten Nutzens benutzt wird. Für die exponentielle Nutzenfunktion berechnen wir Nutzenindifferenzpreise. Ausserdem wenden wir diese Techniken auf die Maximierung des erwarteten Nutzens bezüglich eines festen Wahrscheinlichkeitsmasses an. Dafür betrachten wir abgeschlossene, im allgemeinen nicht konvexe zulässige Mengen für die Handelsstrategien. / The first part of the thesis proposes a method to find prices and hedging strategies for risky claims exposed to a risk factor that is not hedgeable on a financial market. In the second part we calculate the maximal utility and optimal trading strategies on incomplete markets using Backward Stochastic Differential Equations. We consider agents with incomes exposed to a non-hedgeable external source of risk by creating either a bond or by signing contracts. The sources of risk we think of may be insurance, weather or climate risk. Stock prices are seen as exogenuosly given. We calculate prices for the additional securities such that supply is equal to demand, the market clears partially. The preferences of the agents are described by expected utility. In Chapter 2 through Chapter 4 the agents use exponential utility functions, the model is placed in a Brownian filtration. In order to find the equilibrium price, we use Backward Stochastic Differential Equations. Chapter 5 provides a one--period model where the agents use utility functions satisfying the Inada condition. The second part of this thesis considers the robust utility maximization problem on an incomplete financial market. Either the probability measure or drift and volatility of the stock price process are uncertain. We apply a martingale argument and solve a saddle point problem. The solution of a Backward Stochastic Differential Equation describes the maximizing trading strategy as well as the probability measure that is used in the robust utility. We consider the exponential, the power and the logarithmic utility functions. For the exponential utility function we calculate utility indifference prices of not perfectly hedgeable claims. Finally, we maximize the expected utility with respect to a single probability measure. We apply a martingale argument and solve maximization problems. This allows us to consider closed, in general non--convex constraints on the values of trading strategies.

Für zufällige dynamische Systeme mit stetiger Zeit existieren zwei wichtige Klassen von Generatoren: Zum einen stationäre zufällige ifferentialgleichungen, i.e. gewöhnliche Differentialgleichungen, die von einem stationärer zufälligen Vektorfeld getrieben werden, und zum anderen stochastische Stratonovichdifferentialgleichungen mit weißem Rauschen. Während die erste Klasse sich gut in den ergodentheoretischen Rahmen der Theorie der zufälligen dynamischen Systeme einfügt, widersetzte sich die zweite Klasse lange Zeit der dynamischen Untersuchung aufgrund des "Konflikts zwischen Ergodentheorie und stochastischer Analysis". In dieser Arbeit wird gezeigt, daß beide Klassen von zufälligen dynamischen Systemen nicht wesentlich verschieden sind, genauer: Zu jeder stochastischen Stratonovichdifferentialgleichung mit weißem Rauschen (unter den üblichen Regularitätsforderungen an die Vektorfelder, die die Existenz von Flüssen garantieren) existiert eine stationäre zufällige Differentialgleichung derart, daß die erzeugten zufälligen dynamischen Systeme konjugiert sind. Als Anwendung wird eine Version des lokalen Linearisierungssatzes von Hartman/Grobman für stochastische Stratonovichdifferentialgleichungen bewiesen. / For continuous time random dynamical systems there exist two important classes of generators: on the one hand stationary random differential quations, i.e. ordinary differential equations driven by a stationary random vector field, and on the other hand stochastic Stratonovich differential equations with white noise. While the first class fits well into the framework of the theory of random dynamical systems, the second class resisted for a long time the dynamical investigation due to the "conflict between ergodic theory and stochastic analysis". The main result of this thesis is that both classes of random dynamical systems are not essentially distinct, more precisely: For each stochastic Stratonovich differential equation with white noise (under usual regularity assumptions) there exists a stationary random differential equation such that the corresponding random dynamical systems are conjugate. As an application a version of the local Hartman/Grobman theorem for stochastic differential equations is proved.

Using a Lie symmetry group generator and a generalized form of Manale's formula for solving second order ordinary di erential equations, we determine new symmetries for the one and two dimensional heat equations, leading to new solutions. As an application, we test a formula resulting from this approach on thin plate heat conduction. / Applied Mathematics / M. Sc. (Applied Mathematics)

The purpose of this thesis in to introduce the systematic study of symmetries in binary differential equations (BDEs). We formalize the concept of a symmetric BDE, under the linear action of a compact Lie group. One of the main results establishes a formula that relates the algebraic and geometric effects of the occurrence of the symmetry in the problem. Using tools from invariant theory and representation theory for compact Lie groups we deduce the general forms of equivariant binary differential equations under compact subgroups of O(2). A study about the behavior of the invariant straight lines on the configuration of homogeneous BDEs of degree n is done with emphasis on cases in which n = 0 and n = 1. Also for the linear case (n = 1) the equivariant normal forms are presented. Symmetries of linear 1-forms are also studied and related with symmetries of tangent orthogonal vectors fields associated with it. / O objetivo desta tese é introduzir o estudo sistemático de simetrias em equações diferenciais binárias (EDBs). Neste trabalho formalizamos o conceito de EDB simétrica sobre a ação de um grupo de Lie compacto. Um dos principais resultados é uma fórmula que relaciona o efeito geométrico e algébrico das simetrias presentes no problema. Utilizando ferramentas da teoria invariante e de representação para grupos compactos deduzimos as formas gerais para EDBs equivariantes. Um estudo sobre o comportamento das retas invariantes na configuração de EDBs com coeficientes homogêneos de grau n é feito com ênfase nos casos de grau 0 e 1, ainda no caso de grau 1 são apresentadas suas formas normais. Simetrias de 1-formas lineares são também estudadas e relacionadas com as simetrias dos seus campos tangente e ortogonal.

The purpose of this thesis in to introduce the systematic study of symmetries in binary differential equations (BDEs). We formalize the concept of a symmetric BDE, under the linear action of a compact Lie group. One of the main results establishes a formula that relates the algebraic and geometric effects of the occurrence of the symmetry in the problem. Using tools from invariant theory and representation theory for compact Lie groups we deduce the general forms of equivariant binary differential equations under compact subgroups of O(2). A study about the behavior of the invariant straight lines on the configuration of homogeneous BDEs of degree n is done with emphasis on cases in which n = 0 and n = 1. Also for the linear case (n = 1) the equivariant normal forms are presented. Symmetries of linear 1-forms are also studied and related with symmetries of tangent orthogonal vectors fields associated with it. / O objetivo desta tese é introduzir o estudo sistemático de simetrias em equações diferenciais binárias (EDBs). Neste trabalho formalizamos o conceito de EDB simétrica sobre a ação de um grupo de Lie compacto. Um dos principais resultados é uma fórmula que relaciona o efeito geométrico e algébrico das simetrias presentes no problema. Utilizando ferramentas da teoria invariante e de representação para grupos compactos deduzimos as formas gerais para EDBs equivariantes. Um estudo sobre o comportamento das retas invariantes na configuração de EDBs com coeficientes homogêneos de grau n é feito com ênfase nos casos de grau 0 e 1, ainda no caso de grau 1 são apresentadas suas formas normais. Simetrias de 1-formas lineares são também estudadas e relacionadas com as simetrias dos seus campos tangente e ortogonal.