Neste trabalho, estudamos propriedades de continuação única para as soluções da equação tipo Schrödinger com um ponto interação centrado em x=0, \\partial_tu=i(\\Delta_Z+V)u, onde V=V(x,t) é uma função de valor real e -\\Delta_Z é o operador escrito formalmente como \\[-\\Delta_Z=-\\frac\\frac{d^2}{dx^2}+Z\\delta_0,\\] sendo \\delta_0 a delta de Dirac centrada em zero e Z qualquer número real. Logo, usamos estes resultados para ver o possível fenômeno de concentração das soluções, que explodem, da equação de tipo Schrödinger não linear com um ponto de interação em x=0, \\[\\partial_tu=i(\\Delta_Zu+|u|^u),\\] com ho>5. Também, mostramos que para certas condições sobre o potencial dependente do tempo V, a equação linear em cima tem soluções não triviais. / In this work, we study unique continuation properties for solutions of the Schrödinger equations with an point interaction centered at $x=0$, \\begin\\label \\partial_tu=i(\\Delta_Z+V)u, \\end where $V=V(x,t)$ is real value function and $-\\Delta_Z$ is the operator formally written \\[-\\Delta_Z=-\\frac\\frac{d^2}{dx^2}+Z\\delta_0,\\] and $\\delta_0$ is Dirac\'s delta centered at zero and $Z$ is a real number. Next, we use these results in order to study the possible profile of the concentration of blow up solutions for the non linear Schrödinger equation with a point interaction at $x=0$, \\[\\partial_tu=i(\\Delta_Zu+|u|^u),\\] with $ho>5$. Besides, we show that the equation above has non trivial solutions for some conditions on the time dependent potencial $V$.

As equações diferencias parciais tem origem na modelagem do problemas nas ciências e engenharia, tais como a equação do calor, equação da onda, equação de Poisson, entre outras. Para muitas destas equações não é tão simples obter uma técnica analítica para achar sua solução e nestes casos é necessário uso de soluções aproximadas obtidas pelo computador. Existem técnicas tradicionais para solução numérica de uma grande classe de equações diferenciais, mas quando esta equação está na forma implícita, muitas destas técnicas já não podem ser aplicadas. Frequentemente as equações diferenciais parciais de segunda ordem tem maior estudo que as equações de primeira ordem sendo uma das razões que os modelos envolvem derivadas de segunda ordem. No caso das equações diferenciais parciais de primeira ordem implícitas a não linearidade em alguns casos não permite determinar uma solução de forma simples. O trabalho desenvolvido faz uma revisão do método das características para estabelecer as condições necessárias e suficientes, que permitam encontrar uma solução, ao mesmo tempo evidencia a complexidade de determinar uma solução clássica. Dentro das aplicações existentes relacionadas com as Equações Diferenciais Parciais Implícitas de Primeira Ordem, podemos mencionar a Equação cinemática e a Equação de Hamilton-Jacobi que podem-se associar com o movimento de partículas. Para a solução de uma Equação Diferencial Implícita de Primeira Ordem o método das características tem uma estrutura de solução que permite resolver a equação de forma analítica e numérica, desde que se verifique o Teorema de Cauchy. O objetivo deste trabalho de mestrado é obter um método numérico para a solução de equações diferenciais parciais de primeira ordem implícitas. Nós propomos um método numérico do tipo previsor-corretor que resolve uma EDP de primeira ordem implícita, utilizando o sistema característico em conjunto com as condições de banda, para reduzir o erro global nas iterações. / Partial differential equations arise in the modeling of problems in science and engineering, such as the heat equation, wave equation, Poisson equation, among others. For many of these equations it is not so simple to obtain an analytical technique to find a solution in these cases and it is necessary to use a computer to obtain approximate solutions. There are traditional techniques for numerical solution of a large class of differential equations, but when this equation is in implicit form, many of these techniques can no longer be applied. Often partial differential equations of second order are more studied than first order equations the reason being that one of the models involve secondorder derivatives. In the case of implicit partial differential equations of first order the non-linearity in some cases does not allow for a solution in simple from to be determined. The work reviews the method of characteristics to establish the necessary and sufficient conditions that will find a solution at the same time demonstrates the complexity of determining classical solution. Within existing applications related to Partial Differential Equations of First Order Implicit, we can mention the textit kinematic equation and textit equation Hamilton-Jacobi that can be associated with the movement of particles. For the solution of a differential equation First Implicit Order the method of characteristics has a solution framework that enables solve the equation analytically and numerically, provided there is the Cauchy theorem. The objective of this master thesis is to obtain a numerical method for the solution of partial differential equations first order implicit. We propose a numerical method of predictor-corrector type that resolves a EDP first implicate order, using the characteristic system in conjunction with the band conditions, to reduce the overall error in iterations.

As equações diferencias parciais tem origem na modelagem do problemas nas ciências e engenharia, tais como a equação do calor, equação da onda, equação de Poisson, entre outras. Para muitas destas equações não é tão simples obter uma técnica analítica para achar sua solução e nestes casos é necessário uso de soluções aproximadas obtidas pelo computador. Existem técnicas tradicionais para solução numérica de uma grande classe de equações diferenciais, mas quando esta equação está na forma implícita, muitas destas técnicas já não podem ser aplicadas. Frequentemente as equações diferenciais parciais de segunda ordem tem maior estudo que as equações de primeira ordem sendo uma das razões que os modelos envolvem derivadas de segunda ordem. No caso das equações diferenciais parciais de primeira ordem implícitas a não linearidade em alguns casos não permite determinar uma solução de forma simples. O trabalho desenvolvido faz uma revisão do método das características para estabelecer as condições necessárias e suficientes, que permitam encontrar uma solução, ao mesmo tempo evidencia a complexidade de determinar uma solução clássica. Dentro das aplicações existentes relacionadas com as Equações Diferenciais Parciais Implícitas de Primeira Ordem, podemos mencionar a Equação cinemática e a Equação de Hamilton-Jacobi que podem-se associar com o movimento de partículas. Para a solução de uma Equação Diferencial Implícita de Primeira Ordem o método das características tem uma estrutura de solução que permite resolver a equação de forma analítica e numérica, desde que se verifique o Teorema de Cauchy. O objetivo deste trabalho de mestrado é obter um método numérico para a solução de equações diferenciais parciais de primeira ordem implícitas. Nós propomos um método numérico do tipo previsor-corretor que resolve uma EDP de primeira ordem implícita, utilizando o sistema característico em conjunto com as condições de banda, para reduzir o erro global nas iterações. / Partial differential equations arise in the modeling of problems in science and engineering, such as the heat equation, wave equation, Poisson equation, among others. For many of these equations it is not so simple to obtain an analytical technique to find a solution in these cases and it is necessary to use a computer to obtain approximate solutions. There are traditional techniques for numerical solution of a large class of differential equations, but when this equation is in implicit form, many of these techniques can no longer be applied. Often partial differential equations of second order are more studied than first order equations the reason being that one of the models involve secondorder derivatives. In the case of implicit partial differential equations of first order the non-linearity in some cases does not allow for a solution in simple from to be determined. The work reviews the method of characteristics to establish the necessary and sufficient conditions that will find a solution at the same time demonstrates the complexity of determining classical solution. Within existing applications related to Partial Differential Equations of First Order Implicit, we can mention the textit kinematic equation and textit equation Hamilton-Jacobi that can be associated with the movement of particles. For the solution of a differential equation First Implicit Order the method of characteristics has a solution framework that enables solve the equation analytically and numerically, provided there is the Cauchy theorem. The objective of this master thesis is to obtain a numerical method for the solution of partial differential equations first order implicit. We propose a numerical method of predictor-corrector type that resolves a EDP first implicate order, using the characteristic system in conjunction with the band conditions, to reduce the overall error in iterations.

Orientador: Dílson Amâncio Alves / Banca: Luiz Fernando Bovolato / Banca: Luis Carlos Origa de Oliveira / Banca: Vivaldo Fernando da Costa / Banca: Eduardo Nobuhiro Asada / Resumo: Este trabalho apresenta uma metodologia alternativa para a melhoria da margem de carregamento e redução da perda total de potência ativa com base no método da continuação. Para atingir esta meta, uma equação de parametrização baseada na perda de potência ativa total e as equações da potência reativa nas barras de geração são acrescentadas às equações de fluxo de carga convencional. As tensões nas barras PV são consideradas como variáveis de controle e um novo parâmetro é escolhido para reduzir as perdas de potência ativa nas linhas de transmissão. Os resultados mostram que este procedimento, em geral, conduz a um aumento no ponto de máximo carregamento e por conseguinte, melhoria na margem estática da estabilidade de tensão. Este procedimento também leva a uma redução nos custos operacionais e, simultaneamente, uma melhoria no perfil da tensão. / Abstract: This work presents an alternative methodology for loading margin improvement and total real power losses reduction by using a continuation method. In order to attain this goal, a parameterizing equation based on the total real power losses and the equations of the reactive power at the slack and generation buses are added to the conventional Power Flow equations. The voltages at these buses are considered as control variables and a new parameter is chosen with to reduce the real power losses in the transmission lines. The results show that this procedure leads to maximum loading point increase and consequently, in static voltage stability margin improvement. Besides, this procedure also takes to a reduction in the operational costs and, simultaneously, to voltage profile improvement. / Doutor

Made available in DSpace on 2014-06-11T19:22:33Z (GMT). No. of bitstreams: 0 Previous issue date: 2011-01-14Bitstream added on 2014-06-13T18:08:29Z : No. of bitstreams: 1 amorim_es_me_ilha.pdf: 2160628 bytes, checksum: 3be5e9c0540f3af6109d0248dc52ede1 (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / O objetivo do trabalho proposto é o de apresentar aplicações práticas da técnica de análise modal na avaliação da estabilidade de tensão e à redução da perda total de potência reativa na transmissão. A análise modal tem demonstrado ser uma ferramenta de análise estática muito útil que pode ser utilizada para identificar áreas propensas à instabilidade de tensão. Seus resultados são usados para identificar as medidas de reforço mais adequadas para aumentar a margem de estabilidade de tensão de sistemas elétricos de potência. A melhor localização para a instalação de compensação reativa shunt e série é determinada baseada nos fatores de participação de barra e de ramos, respectivamente. Outros objetivos deste trabalho são os de apresentar uma metodologia alternativa para a alocação de bancos de capacitores shunt, e o efeito da compensação de bancos de capacitores série na margem de carregamento. Os objetivos são à redução da perda total de potência reativa na transmissão e o aumento na margem de carregamento do sistema. Os valores dos montantes de compensação reativa shunt a serem alocados são determinados com base nas curvas de perda total de potência reativa na transmissão versus magnitude de tensão da barra. No caso da compensação série, o valor da reatância capacitiva escolhido é limitado pelo valor prático de compensação que é da ordem de 80%. Também são investigados os efeitos individuais que cada um destes procedimentos traz para o sistema. Entre os efeitos investigados estão: a redução das perdas reativas, a variação na margem de estabilidade estática de tensão, e a melhoria do perfil de tensão. Os resultados obtidos para os sistemas do IEEE (14, 57 e 118 barras) mostram que o procedimento conduz a uma sensível redução da perda total de potência... / The objective of the proposed work is to present practical applications of the modal analysis technique applied to voltage stability assessment and reactive power losses reduction. Modal analysis has proven to be a useful steady-state analysis approach which can be applied to identify areas prone to voltage instability. Its results are used to identify the most effective remedial actions to increase system voltage stability margin. The best location for placing shunt and series capacitor is determined based on the bus and branch participations factors respectively. Other goals of the proposed work are to present an alternative procedure to the allocation of shunt capacitor, and the effects of series capacitor compensation on the system loading margin. The objectives are the total reactive power losses reduction in the transmission and the loading margin increase. The values of shunt reactive compensation to be allocated are determined based on the curve of total transmission reactive power losses versus bus voltage magnitude. In case of series compensation, the chosen reactive capacitance values are limited by a practical upper limit of series compensation degree of 80%. The individual effects that each one of these procedures brings for the system are also investigated. Among the investigated effects are: the reduction in the total reactive power losses, the changes in the static voltage stability margin, and the improvement of voltage profile. The results obtained for the IEEE system (14, 57 and 118 bus) show that the procedure leads to a sensible reduction of total reactive power losses and simultaneously, an improvement in the voltage profile and load margin increase.

Made available in DSpace on 2015-05-15T11:46:02Z (GMT). No. of bitstreams: 1 arquivototal.pdf: 425597 bytes, checksum: 2cb5c42e376ba77764139fbc3fcdaf15 (MD5) Previous issue date: 2011-07-02 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / We consider the semilinear heat equation involving gradient terms in a bounded domain of Rn. It is assumed the non-linearity is globally Lipschitz. We prove that the system is approximately controllable when the control acts on a bounded subset of the domain. The proof uses a variant of a classical fixed point method and is a simpler alternative to the earlier proof existing in the literature by means of the penalization of an optimal control problem. We also prove that the control may be built so that, in addition to the approximate controllability requirement, it ensures that the state reaches exactly a finite number of constraints. / Consideremos a equação do calor semilinear envolvendo termos do gradiente em um domínio limitado do Rn. Assumimos que a não-linearidade é globalmente Lipschtz. Usando o método do ponto fixo, provamos que o sistema é finito-aproximadamente controlável e nulamente controlável, quando o controle age em um subconjunto não vazio do domínio.

Nas últimas décadas, considerável evidência tem sido direcionada no desenvolvimento de métodos computacionais que analisam os comportamentos não lineares das estruturas. O método da corda é considerado o tipo de método de seguimento de trajetórias que possui a técnica mais comum e versátil para analisar tais comportamentos não lineares. No entanto, testes realizados pelos autores concluem que tais métodos ainda apresentam algumas dificuldades quando as estruturas possuem formas complexas e comportamentos fortemente não lineares. Devido a isso, o objetivo do trabalho em questão é acrescentar duas modificações ao método da corda, a fim de desenvolver um modelo mais seguro e estável. A esses métodos, inclui-se (i) uma equação alternativa para o passo previsor, que tornará o método mais estável ao ultrapassar pontos críticos; e (ii) um parâmetro de controle, que através do conceito de soma dos comprimentos de corda, é chamado de método da corda acumulado. Ao final, por consequência destas alterações, um procedimento mais robusto para o método da corda é obtido. Sua validação é apresentada, pelos autores, através de exemplos utilizando estruturas treliçadas (bidimensionais e tridimensionais) em um programa de simulação numérica desenvolvida. / In the last decades, a considerable effort has been directed to develop of computational methods to analyze the nonlinear behavior of structures. The arc-length method is considered the type of path following method which have the most common and versatile techniques for analyzing such nonlinear behaviors. Nevertheless, tests performed by the authors conclude that such methods present some difficulties when the structures possess complex shapes and strongly nonlinear behaviors. Due to this, the purpose of the present work is to add two amendments to the arc-length method, in order to develop a model safer and more stable. In such methods, it is included (i) an alternative equation for the predictor step, that will make the method most stable to overcoming critical points, and (ii) a control parameter, through the concept of sum of the arc length up to now, referred to as accumulated arc-length method. By the end, as a result of these modifications, a more robust procedure for the method of the string is obtained. And its validation is presented by the authors through examples using truss structures in a two-dimensional and three-dimensional numerical simulation program developed.

Orientador: Alberto Luiz Serpa / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-05T19:36:12Z (GMT). No. of bitstreams: 1 Paulino_LeonardoPereira_M.pdf: 2501868 bytes, checksum: 87eef8b309aa247d3e5c6c4f572b1323 (MD5) Previous issue date: 2004 / Resumo: Em muitos problemas não lineares na mecânica estrutural é necessário o emprego de métodos de solução que controlam os deslocamentos e o nível de carregamento simultaneamente. Um dos métodos mais conhecidos para estas situações é o método do arc-length. Este método consiste na introdução de uma restrição (" esfera generalizada", por exemplo) à equação de equilíbrio não linear do problema. A solução do problema de equilíbrio juntamente com a restrição" esférica" mantém o deslocamento limitado, evitando a divergência durante o processo de resolução. O método clássico do arc-length apresenta certa complexidade por envolver duas fases de cálculo, predição e correção da solução, além do fato de que mais de uma solução satisfaz as equações envolvidas, sendo necessário o emprego de alguns critérios ainda não suficientemente robustos para a escolha da solução adequada. Propõe-se neste trabalho o estudo de uma formulação alternativa para o método do arc-length através do emprego dos conceitos de programação matemática, onde um problema de minimização com restrições é formulado, minimizando-se uma função objetivo que caracterize o problema de equilíbrio e esteja sujeita a restrições que garantam o controle dos deslocamentos. Além disso, uma restrição adicional é imposta ao problema para buscar a unicidade da solução. A verificação do desempenho da metodologia proposta é feita através da solução de problemas de estruturas reticuladas envolvendo efeitos de grandes deslocamentos. A implementação foi feita em Matlab pela facilidade e disponibilidade de ferramentas computacionais tais como as funções do "toolbox" de otimização / Abstract: Non-linear problems in structural mechanics in general require the use of solution methods that control the displacements and the load leveI simultaneously. One of the most popular method used in these cases is the arc-length method. The arc-length method introduces one additional constraint equation to the non-linear equilibrium equations of the problem. The solution of the equilibrium equations with the additional constraints should keep a limited displacement to avoid numerical divergence during the solution procedure. The classical arc length method involves two solution phases, prediction and correction. A difficulty of this method is that it can re-compute a solution already determined, requiring the use of some cri teria that are not robust enough to choose the adequated solution. This work proposes an alternative formulation for the arc-length method using the concepts of mathematical programming, where a constrained minimization problem is formulated. The objective function is established in terms of the equilibrium residue and the arc-length constraint should ensure displacements and load leveI control. Besides this, an additional constraint equation is employed to guarantee the uniqueness of the solution, improving the solution path compared to the classical arc-length method. To verify the proposed methodology computational performance some problems involving the effects of large displacements in truss structures are presented. The implementation was done in the software Matlab due to its facility and the computational tools that are available, such as the optimization toolbox functions / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Neste trabalho; estudamos a dinâmica das órbitas de uma partícula infinitesimal P contida no espaço euclidiano tridimensional; atraída unicamente pela força gravitacional induzida por um corpo maciço com densidade de massa constante na forma de anel ou disco circular. Este problema denominamos: Problema do anel ou disco circular homogêneo. Este problema apresenta alguns subproblemas: a linha vertical perpendicular ao plano que contém anel ou disco circular e passa pelo centro de massa; a qual denominamos de eixo-z ; o plano que contém o anel ou disco circular; que denominamos de plano horizontal e o plano perpendicular ao plano horizontal e que contém o eixo-z. Obtivemos resultados importantes da dinâmica em cada um destes subproblemas. Mostramos que o problema do anel ou disco circular é invariante por rotações em torno do eixo-z e desta forma podemos formular o problema em coordenadas giratórias. Escrevemos o problema como uma perturbação do problema de Kepler onde o parâmetro perturbador está associado a espessura do anel ou raio do disco. Utilizando a formulação do problema em coordenadas convenientes conseguimos obter uma grande quantidade de órbitas periódicas simétricas; como continuação de órbitas circulares e também elípticas no caso espacial e no caso dos subproblemas planares. A técnica empregada para conseguir tais órbitas foi o método de continuação analítica de Poincaré primeiramente aplicado a um problema geral que consiste em uma perturbação do problema de Kepler analítica e depois aplicamos ao nosso problema particular. Ainda estudamos as órbitas de escapes associadas a nosso problema; concentrando o escape nas direções dos eixos coordenados. _________________________________________________________________________________________ ABSTRACT: Our main concern, in this work, is describe the dynamics of the orbits of an in¯nitesimal particle moving in the space R3 under the in°uence of the gravitational force induced by a homogeneous annulus disk ¯xed on a plane. The aspects of the dynamics that we have interest are namely: Describe the di®erent ways the gravitational potential associated to the problem in each case; Characterize homogeneity properties of the potential; Describe the con¯guration space of these problems; Determine the symmetries of the vectorial ¯eld associate; Identify sub-problems associates according the dimension of the ambient space; In each sub-case,describe the dynamics and compares them to each other; Relate the potential singularities with the singularities of the vectorial ¯eld of each problem at issue; Introduce a convenient parameter; Determine a great diversity of families of periodic orbits in the di®erent sub-problems. Study the escape orbits in the di®erent cases; Compare the results obtained with the n-body problem in Celestial Mechanics.

Neste trabalho, estudamos propriedades de continuação única para as soluções da equação tipo Schrödinger com um ponto interação centrado em x=0, \\partial_tu=i(\\Delta_Z+V)u, onde V=V(x,t) é uma função de valor real e -\\Delta_Z é o operador escrito formalmente como \\[-\\Delta_Z=-\\frac\\frac{d^2}{dx^2}+Z\\delta_0,\\] sendo \\delta_0 a delta de Dirac centrada em zero e Z qualquer número real. Logo, usamos estes resultados para ver o possível fenômeno de concentração das soluções, que explodem, da equação de tipo Schrödinger não linear com um ponto de interação em x=0, \\[\\partial_tu=i(\\Delta_Zu+|u|^u),\\] com ho>5. Também, mostramos que para certas condições sobre o potencial dependente do tempo V, a equação linear em cima tem soluções não triviais. / In this work, we study unique continuation properties for solutions of the Schrödinger equations with an point interaction centered at $x=0$, \\begin\\label \\partial_tu=i(\\Delta_Z+V)u, \\end where $V=V(x,t)$ is real value function and $-\\Delta_Z$ is the operator formally written \\[-\\Delta_Z=-\\frac\\frac{d^2}{dx^2}+Z\\delta_0,\\] and $\\delta_0$ is Dirac\'s delta centered at zero and $Z$ is a real number. Next, we use these results in order to study the possible profile of the concentration of blow up solutions for the non linear Schrödinger equation with a point interaction at $x=0$, \\[\\partial_tu=i(\\Delta_Zu+|u|^u),\\] with $ho>5$. Besides, we show that the equation above has non trivial solutions for some conditions on the time dependent potencial $V$.