Many real-life applications such as the shape optimization of technological devices, the identification of parameters in environmental processes and flow control problems lead to optimization problems governed by systems of convection diusion partial dierential equations (PDEs). When convection dominates diusion, the solutions of these PDEs typically exhibit layers on small regions where the solution has large gradients. Hence, it requires special numerical techniques, which take into account the structure of the convection. The integration of discretization and optimization is important for the overall eciency of the solution process. Discontinuous Galerkin (DG) methods became recently as an alternative to the finite dierence, finite volume and continuous finite element methods for solving wave dominated problems like convection diusion equations since they possess higher accuracy. This thesis will focus on analysis and application of DG methods for linear-quadratic convection dominated optimal control problems. Because of the inconsistencies of the standard stabilized methods such as streamline upwind Petrov Galerkin (SUPG) on convection diusion optimal control problems, the discretize-then-optimize and the optimize-then-discretize do not commute. However, the upwind symmetric interior penalty Galerkin (SIPG) method leads to the same discrete optimality systems. The other DG methods such as nonsymmetric interior penalty Galerkin (NIPG) and incomplete interior penalty Galerkin (IIPG) method also yield the same discrete optimality systems when penalization constant is taken large enough. We will study a posteriori error estimates of the upwind SIPG method for the distributed unconstrained and control constrained optimal control problems. In convection dominated optimal control problems with boundary and/or interior layers, the oscillations are propagated downwind and upwind direction in the interior domain, due the opposite sign of convection terms in state and adjoint equations. Hence, we will use residual based a posteriori error estimators to reduce these oscillations around the boundary and/or interior layers. Finally, theoretical analysis will be confirmed by several numerical examples with and without control constraints

The study of wave production and propagation is a common phenomenon seen within a variety of math and physics problems. This thesis in particular will investigate the production and propagation of sound waves through musical instruments. Although this field of work has been examined since the late 1800s, approaching these types of problems can be very difficult. With the exception of the last fifty years, we have only been able to approach such problems by linearizing the necessary equations of gas dynamics. Without the use of a computer, one can only get so far in studying nonlinear acoustic problems. In addition, the numerical theory for nonlinear problems is incomplete. Proving stability is challenging and there are a variety of open problems within this field. This thesis will be examining the propagation of sound waves specifically through brass instruments. However, we will not be able to fully examine this problem in a master’s thesis because of the complexity. Instead, the objective is to provide a foundation and global picture of this problem by merge the fields of nonlinear acoustics as well as computational and analytical gas dynamics. To study the general behaviour of nonlinear wave propagation (and to verify previous findings), experiments have been carried on a trumpet. The purpose of these experiments is take measurements of the sound pressure waves at various locations along the instrument in order to understand the evolution of the wave propagation. In particular, we want to establish if the nonlinear distortion is strong enough to have musical consequences; and if there are such outcomes, what prerequisites are required for the observable behaviour. Additionally, by using the discontinuous Galerkin numerical method, a model of the system will be presented in this thesis. It will then be compared with the experimental data to verify how well we were able to describe the nonlinear wave motion within a trumpet.

The physical processes in lakes remain only partially understood despite successful data collection from a variety of sources spanning several decades. Although numerical models are already frequently employed to simulate the physics of lakes, especially in the context of water quality management, improved methods are necessary to better capture the wide array of dynamically important physical processes, spanning length scales from ~ 10 km (basin-scale oscillations) - 1 m (short internal waves). In this thesis, high-order numerical methods are explored for specialized model equations of lakes, so that their use can be taken into consideration in the next generation of more sophisticated models that will better capture important small scale features than their present day counterparts. The full three-dimensional incompressible density-stratified Navier-Stokes equations remain too computationally expensive to be solved for situations that involve both complicated geometries and require resolution of features at length-scales spanning four orders of magnitude. The main source of computational expense lay with the requirement of having to solve a three-dimensional Poisson equation for pressure at every time-step. Simplified model equations are thus the only way that numerical lake modelling can be carried out at present time, and progress can be made by seeking intelligent parameterizations as a means of capturing more physics within the framework of such simplified equation sets. In this thesis, we employ the long-accepted practice of sub-dividing the lake into vertical layers of different constant densities as an approximation to continuous vertical stratification. We build on this approach by including weakly non-hydrostatic dispersive correction terms in the model equations in order to parameterize the effects of small vertical accelerations that are often disregarded by operational models. Favouring the inclusion of weakly non-hydrostatic effects over the more popular hydrostatic approximation allows these models to capture the emergence of small-scale internal wave phenomena, such as internal solitary waves and undular bores, that are missed by purely hydrostatic models. The Fourier and Chebyshev pseudospectral methods are employed for these weakly non-hydrostatic layered models in simple idealized lake geometries, e.g., doubly periodic domains, periodic channels, and annular domains, for a set of test problems relevant to lake dynamics since they offer excellent resolution characteristics at minimal memory costs. This feature makes them an excellent benchmark to compare other methods against. The Discontinuous Galerkin Finite Element Method (DG-FEM) is then explored as a mid- to high-order method that can be used in arbitrary lake geometries. The DG-FEM can be interpreted as a domain-decomposition extension of a polynomial pseudospectral method and shares many of the same attractive features, such as fast convergence rates and the ability to resolve small-scale features with a relatively low number of grid points when compared to a low-order method. The DG-FEM is further complemented by certain desirable attributes it shares with the finite volume method, such as the freedom to specify upwind-biased numerical flux functions for advection-dominated flows, the flexibility to deal with complicated geometries, and the notion that each element (or cell) can be regarded as a control volume for conserved fluid quantities. Practical implementation details of the numerical methods used in this thesis are discussed, and the various modelling and methodology choices that have been made in the course of this work are justified as the difficulties that these choices address are revealed to the reader. Theoretical calculations are intermittently carried out throughout the thesis to help improve intuition in situations where numerical methods alone fall short of giving complete explanations of the physical processes under consideration. The utility of the DG-FEM method beyond purely hyperbolic systems is also a recurring theme in this thesis. The DG-FEM method is applied to dispersive shallow water type systems as well as incompressible flow situations. Furthermore, it is employed for eigenvalue problems where orthogonal bases must be constructed from the eigenspaces of elliptic operators. The technique is applied to the problem calculating the free modes of oscillation in rotating basins with irregular geometries where the corresponding linear operator is not self-adjoint.

Thesis (Doctora).

Ce travail de thèse s’inscrit dans l’effort de réduction des nuisances sonores dues à la soufflante d’unréacteur double-flux à l’aide de matériaux absorbants acoustiques, appelés communément «liners». Afind’optimiser ces traitements acoustiques, il convient d’étudier en détail la physique de la propagationacoustique en présence de liner. De plus, il s’agit d’améliorer la compréhension des instabilités hydrodynamiquespouvant se développer sur un liner sous des conditions particulières et possiblement génératricesde bruit. Ce travail de thèse a consisté à développer un code de calcul en formulation Galerkin discontinuepour l’analyse modale et la stabilité dans un conduit traité acoustiquement, code qui a été appliqué à desconfigurations réalistes, en considérant une section transverse ou longitudinale d’un conduit. Les étudesmodales réalisées dans la section transverse ont apporté des informations sur la propagation acoustiquedans une nacelle de turbofan avec des discontinuités du traitement acoustique («splices»), ainsi que dansle banc B2A de l’ONERA. Les calculs dans la section longitudinale ont nécessité l’implantation de conditionsaux limites PML pour tronquer le domaine de calcul, ainsi que d’une condition aux limites sur leliner, modélisée en domaine temporel à partir d’une extension de travaux existants dans la littérature.Avec ces outils, le code a permis de mettre en évidence une dynamique de type amplificateur de bruit dueau développement d’une instabilité hydrodynamique sur le liner en présence d’écoulement cisaillé ainsiqu’un rayonnement acoustique en amont et en aval du conduit dû à cette instabilité. / The current work deals with the reduction of aircraft engine fan noise using acoustic lining. In orderto optimise these liners, it is necessary to deeply understand the physics of acoustic wave propagation in lined ducts and to have a better knowledge of the hydrodynamic instabilities existing under particular conditions and likely to radiate noise. This work is about the development of a discontinuous Galerkin solver for modal and stability analysis in lined flow duct and the application of this solver to realistic configurations by considering the transverse or longitudinal section of a duct. The modal studies in the transverse section brought informations on acoustic propagation in a turbofan nacelle with lining discontinuities (“splices”) and in the B2A bench of ONERA. The computation in the longitudinal section of a duct required the implementation of PML boundary conditions in order to truncate the computational domain and of a boundary condition at the lined wall, modeled in temporal domain by the enhancement of a method published in the literature. With these features, the application of the solver highlighted a noise amplifier dynamics caused by the development of a hydrodynamic instability on the liner with sheared flow and a noise radiation mechanism upstream and downstream the lined section.

O estudo da dinâmica estrutural de risers oceânicos apresenta instigantes desafios aos pesquisadores da área da engenharia de estruturas, uma vez que os meios tradicionais de análises dinâmicas lineares nem sempre se ajustam às suas complexas particularidades. No atual estágio do desenvolvimento científico da área de engenharia de estruturas, a aplicação de técnicas de análise dinâmica não linear, dentro de determinadas hipóteses, mostra-se como uma das alternativas possíveis e viáveis à tradicional análise dinâmica linear. Com vistas a uma nova abordagem do problema, o presente trabalho adota uma metodologia de análise não linear dinâmica de risers oceânicos em configuração de lançamento de catenária, conjugada a uma técnica de processamento de Modelos de Ordem Reduzida para o estudo dos fenômenos dinâmicos manifestados por risers. Trata-se de um método de modelagem local, restrito à região de contato unilateral do riser com o solo, considerado este último um meio viscoelástico. Os resultados da aplicação desta metodologia são demonstrados nos estudos de caso apresentados com comparações com modelos numéricos (Método dos Elementos Finitos) e modelos físicos. / The dynamic study of offshore risers still demands large efforts from structural engineering researchers, since these systems may behave in a way that is not well modeled and understood using simply linear dynamic theories. Nevertheless, the current development stage of non linear dynamic theories gives hope that their use for the analyses of such systems can be of great value, even though, this must be carefully done specially by the analyst. The present work refers to a non linear dynamic methodology application to offshore risers, particularly steel catenary risers, by a technique known as reduced-order modeling, in the study of dynamic phenomena that these structures may present. The model is local, which means that it represents the touch-down zone of the riser-soil system. The soil modeling was presumed to be viscoelastic. The results obtained in case studies are compared with those from numerical (Finite Element Method) and small scale physical models.

Submitted by Irene Nascimento (irene.kessia@ufpe.br) on 2017-04-11T18:59:11Z No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) DissertaçãoOmar.pdf: 629619 bytes, checksum: f018416fe978f2e27de6abfe2542c60c (MD5) / Made available in DSpace on 2017-04-11T18:59:11Z (GMT). No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) DissertaçãoOmar.pdf: 629619 bytes, checksum: f018416fe978f2e27de6abfe2542c60c (MD5) Previous issue date: 2016-02-19 / CNPQ / Estudamos aspectos teóricos de um sistema que modela o comportamento dos unidos micro polares incompressíveis num domínio limitado _ Rn (n = 2 ou 3). Especificamente, utilizamos o método espectral de Galerkin para mostrar a existência de soluções fortes e com determinadas condições mostramos a unicidade das soluções / We study theoretical aspects of a system that models the behavior of incompressible micropolar uids in a bounded domain _ Rn (n = 2 or 3). Speci cally, we use the spectral Galerkin method to show the existence of strong solutions and under certain conditions show the uniqueness of solutions.

Orientador: Philippe Remy Bernard Devloo / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica e Instituto de Geociências / Made available in DSpace on 2018-08-23T03:59:08Z (GMT). No. of bitstreams: 1 DuranTriana_OmarYesid_M.pdf: 4774391 bytes, checksum: f29f87826d760e9d63ab1efb621493f4 (MD5) Previous issue date: 2013 / Resumo: A reativação de falhas, resultante de variações na pressão de poros, pode ocasionar atividades sísmicas, subsidência, dano nos poços e criação de caminhos de escape dos fluidos contidos nos reservatórios. Para se garantir uma produção de hidrocarbonetos eficiente, mostra-se um fator crítico a avaliação das tendências de reativação das falhas existentes no meio poroso. Neste trabalho, apresenta-se uma aproximação numérica para uma análise de deformação quase-estática com escoamento monofásico, considerando a compressibilidade da rocha e dos fluidos. Um modelo bidimensional foi empregado considerando a teoria de poroelasticidade linear e um novo tratamento da poroelasticidade através de estruturas de dados multifísicos. Formas adimensionais das equações de poroelasticidade são apresentadas, juntamente com a reprodução de diversas soluções analíticas e semi-analíticas das mesmas em semiespaços, com o propósito de se validar o algoritmo desenvolvido. O modelo computacional foi utilizado para avaliar as mudanças de tensão, no reservatório e em suas fronteiras, com o objetivo de se estudar as tendências reativação de falhas em diferentes cenários. As tendências de reativação de falhas, resultantes da indução de variações de tensão na rocha, foram calculadas através do método de variação de tensão de ruptura de Coulomb para a definição das seções com potencial de deslocamento por tensões cisalhante das falhas pre-existentes em tempos específicos, associados com as alterações na pressão de poros. Mostrou-se que a reativação de falhas depende da geometria de reservatório, das propriedades poroelásticas da rocha, coeficiente de atrito e a distribuição da pressão de poros. Um estudo sobre precisão dos cálculos baseado na dimensão do material circundante é apresentada e vários cenários com diferentes programas depleção foram avaliados para determinar a influência das taxas de produção sobre a tendência de reativação das falhas / Abstract: Fault reactivation resulting from pore pressure changes may be accompanied by seismic activity, subsidence, well damage and the creation of fluid leakage paths. To ensure acceptable reservoir performance in hydrocarbon production, it is critical to assess the reactivation tendencies of existing faults. In this work, a numerical approximation is presented that allows quasi-static deformation coupled with monophasic flow considering compressible constituents. Two dimensional modeling is carried out using the theory of linear poroelasticity and a new treatment of poroelastic equations defined into a multiphysics data structure. Dimensionless forms of poroelasticity equations are presented and several analytic and semi analytic solutions, as well as poroelastic inclusion theory were reproduced with the proposed implementation in order to validate it. The computational model is used to evaluate the stress changes around and into the reservoir in order to assess the fault reactivation tendency at different scenarios. Fault reactivation tendency resulting from induced stress changes was calculated using the Coulomb failure stress change method for definition of the shear slip potential along pre-existing faults at one specific time associated to pore pressure change. It was found that fault reactivation tendency depends on the reservoir geometry, poroelastic properties of the reservoir and surrounding rocks, reservoir geometry, static friction coefficient, and pore pressure distribution. A numerical study about the accuracy of surrounding material dimensions is presented and several scenarios with different depletion programs were evaluated to determine the influence of the production rates over fault reactivation tendency / Mestrado / Explotação / Mestre em Ciências e Engenharia de Petróleo

The thesis deals with the numerical solution of partial differential equati- ons describing the flow of the so-called shallow water neglecting the flow in the vertical direction. These equations are of hyperbolical type of the first or- der with a reactive term representing the bottom topology. We discretize the resulting system of equations by the implicit space-time discontinuous Ga- lerkin method (STDGM). In the literature, the explicit techniques are used most of the time. The implicit approach is suitable especially for adaptive methods, because it allows the usage of different meshes for different time niveaus. In the thesis we derive the corresponding method and an adaptive algorithm. Finally, we present usage of the method in several examples. 1

Wave propagation problems can be modeled by partial differential equations. In this thesis, we study wave propagation in fluids and in solids, modeled by the acoustic wave equation and the elastic wave equation, respectively. In real-world applications, waves often propagate in heterogeneous media with complex geometries, which makes it impossible to derive exact solutions to the governing equations. Alternatively, we seek approximated solutions by constructing numerical methods and implementing on modern computers. An efficient numerical method produces accurate approximations at low computational cost. There are many choices of numerical methods for solving partial differential equations. Which method is more efficient than the others depends on the particular problem we consider. In this thesis, we study two numerical methods: the finite difference method and the discontinuous Galerkin method. The finite difference method is conceptually simple and easy to implement, but has difficulties in handling complex geometries of the computational domain. We construct high order finite difference methods for wave propagation in heterogeneous media with complex geometries. In addition, we derive error estimates to a class of finite difference operators applied to the acoustic wave equation. The discontinuous Galerkin method is flexible with complex geometries. Moreover, the discontinuous nature between elements makes the method suitable for multiphysics problems. We use an energy based discontinuous Galerkin method to solve a coupled acoustic-elastic problem.