Die Dynamik der Atmosphäre der Erde umfasst einen Bereich von mikrophysikalischer Turbulenz über konvektive Prozesse und Wolkenbildung bis zu planetaren Wellenmustern. Für Wettervorhersage und zur Betrachtung des Klimas über Jahrzehnte und Jahrhunderte ist diese Gegenstand der Modellierung mit numerischen Verfahren. Mit voranschreitender Entwicklung der Rechentechnik sind Neuentwicklungen der dynamischen Kerne von Klimamodellen, die mit der feiner werdenden Auflösung auch entsprechende Prozesse auflösen können, notwendig. Der dynamische Kern eines Modells besteht in der Umsetzung (Diskretisierung) der grundlegenden dynamischen Gleichungen für die Entwicklung von Masse, Energie und Impuls, so dass sie mit Computern numerisch gelöst werden können. Die vorliegende Arbeit untersucht die Eignung eines unstetigen Galerkin-Verfahrens niedriger Ordnung für atmosphärische Anwendungen. Diese Eignung für Gleichungen mit Wirkungen von externen Kräften wie Erdanziehungskraft und Corioliskraft ist aus der Theorie nicht selbstverständlich. Es werden nötige Anpassungen beschrieben, die das Verfahren stabilisieren, ohne sogenannte „slope limiter” einzusetzen. Für das unmodifizierte Verfahren wird belegt, dass es nicht geeignet ist, atmosphärische Gleichgewichte stabil darzustellen. Das entwickelte stabilisierte Modell reproduziert eine Reihe von Standard-Testfällen der atmosphärischen Dynamik mit Euler- und Flachwassergleichungen in einem weiten Bereich von räumlichen und zeitlichen Skalen. Die Lösung der thermischen Windgleichung entlang der mit den Isobaren identischen charakteristischen Kurven liefert atmosphärische Gleichgewichtszustände mit durch vorgegebenem Grundstrom einstellbarer Neigung zu(barotropen und baroklinen)Instabilitäten, die für die Entwicklung von Zyklonen wesentlich sind. Im Gegensatz zu früheren Arbeiten sind diese Zustände direkt im z-System(Höhe in Metern)definiert und müssen nicht aus Druckkoordinaten übertragen werden.Mit diesen Zuständen, sowohl als Referenzzustand, von dem lediglich die Abweichungen numerisch betrachtet werden, und insbesondere auch als Startzustand, der einer kleinen Störung unterliegt, werden verschiedene Studien der Simulation von barotroper und barokliner Instabilität durchgeführt. Hervorzuheben ist dabei die durch die Formulierung von Grundströmen mit einstellbarer Baroklinität ermöglichte simulationsgestützte Studie des Grades der baroklinen Instabilität verschiedener Wellenlängen in Abhängigkeit von statischer Stabilität und vertikalem Windgradient als Entsprechung zu Stabilitätskarten aus theoretischen Betrachtungen in der Literatur. / The dynamics of the Earth’s atmosphere encompass a range from microphysical turbulence over convective processes and cloud formation up to planetary wave patterns. For weather forecasting and the investigation of climate over decades and centuries, these are subject to modelling with numerical methods. With progressing development of computer technology, re-development of the dynamical cores of climate models is in order to properly handle processes covered by the increasing resolution. The dynamical core of a model consists of the adaptation(discretization)of the basic equations for the dynamics of mass, energy and momentum for solving them numerically employing computers. The presented work investigates the applicability of a low-order Discontinuous Galerkin (DG) method for atmospheric applications. With equations that include external forces like gravitation and the Coriolis force, that is not given by theory. Necessary changes for stabilizing the method without resorting to slope limiters are presented. For the unmodified method, the basic inability to properly keep atmospheric balances is demonstrated. The developed stabilized model reproduces a set of standard test cases in a wide range of spatial and temporal scales. The solution of the termal wind equation along its characteristics curves, those being identical to the isobars, produces balanced atmospheric states with tunable (barotropic and baroclinic) instability via a prescribed zonal wind field. The constructed instability directly relates to the generation of cyclones. In contrast to earlier works, these balanced states are directly given in the z system (height in meters), without need for elaborate conversion from pressure coordinates. With these constructed states, both as reference state, the deviations from which being considered numerically, and as especially as initial condition subject to a small perturbation, several studies of barotropic and baroclinic instability are conducted via simulations. Particularily, the construction of steady states with configurable zonal flows of certain baroclinity facilitates a simulation-based study of baroclinic instability of differing wavelengths, depending on static stability and vertical wind gradient, in correspondence with stability maps from theoretical considerations in the literature.

This study is concerned with the numerical solution of certain stochastic models of chemical reactions. Such descriptions have been shown to be useful tools when studying biochemical processes inside living cells where classical deterministic rate equations fail to reproduce actual behavior. The main contribution of this thesis lies in its theoretical and practical investigation of different methods for obtaining numerical solutions to such descriptions. In a preliminary study, a simple but often quite effective approach to the moment closure problem is examined. A more advanced program is then developed for obtaining a consistent representation of the high dimensional probability density of the solution. The proposed method gains efficiency by utilizing a rapidly converging representation of certain functions defined over the semi-infinite integer lattice. Another contribution of this study, where the focus instead is on the spatially distributed case, is a suggestion for how to obtain a consistent stochastic reaction-diffusion model over an unstructured grid. Here it is also shown how to efficiently collect samples from the resulting model by making use of a hybrid method. In a final study, a time-parallel stochastic simulation algorithm is suggested and analyzed. Efficiency is here achieved by moving parts of the solution phase into the deterministic regime given that a parallel architecture is available. Necessary background material is developed in three chapters in this summary. An introductory chapter on an accessible level motivates the purpose of considering stochastic models in applied physics. In a second chapter the actual stochastic models considered are developed in a multi-faceted way. Finally, the current state-of-the-art in numerical solution methods is summarized and commented upon.

Cette thèse se concentre sur le développement d'une méthode précise et efficace pour la simulation des grandes échelles (LES) des écoulements turbulents. Une approche de la LES basée sur la méthode variationnelle multi-échelles (VMS) est considérée. La VMS applique aux équations de la dynamique des fluides une séparation d'échelles a priori sans recours à des hypothèses sur les conditions aux limites ou sur l'uniformité du maillage. Afin d'assurer effectivement une séparation d'échelles dans l'espace des nombres d'onde associé, nous choisissons d'utiliser les ondelettes de deuxième génération (SGW), une base polynomiale qui présente des propriétés de localisation spatiale-fréquence optimales. A partir de la séparation d'échelles ainsi réalisée, l'action du modèle sous-maille est limitée à un intervalle de nombres d'onde proche de la coupure spectrale. Cette approche VMS-LES basée sur les ondelettes est désignée par WAVVMS-LES. Elle est incorporée dans un solveur d'ordre élevé pour la simulation des écoulements incompressibles sur la base d'une méthode de Galerkin discontinue (DG-FEM) stabilisée pour la pression. La méthode est évaluée par réalisation de LES sur des maillages fortement sous-résolus pour le cas test du tourbillon de Taylor-Green 3D à deux nombres de Reynolds différents. / This thesis focuses on the development of an accurate and efficient method for performing Large-Eddy Simulation (LES) of turbulent flows. An LES approach based upon the Variational Multiscale (VMS) method is considered. VMS produces an a priori scale-separation of the governing equations, in a manner which makes no assumptions on the boundary conditions and mesh uniformity. In order to ensure that scale-separation in wavenumber is achieved, we have chosen to make use of the Second Generation Wavelets (SGW), a polynomial basis which exhibits optimal space-frequency localisation properties. Once scale-separation has been achieved, the action of the subgrid model is restricted to the wavenumber band closest to the cutoff. We call this approach wavelet-based VMS-LES (WAV-VMS-LES). This approach has been incorporated within the framework of a high-order incompressible flow solver based upon pressure-stabilised discontinuous Galerkin FEM (DG-FEM). The method has been assessed by performing highly under-resolved LES upon the 3D Taylor-Green Vortex test case at two different Reynolds numbers.

L'étude de la susceptibilité électromagnétique de circuits électroniques nécessite l'utilisation d'un outil de simulation rapide, précis et suffisamment flexible pour intégrer les dernières innovations technologiques. La méthode itérative basée sur le concept d'onde (notée WCIP pour Wave Concept Iterative Procedure) initialement proposée par H. Baudrand est particulièrement adaptée pour la modélisation numérique de circuits multicouches à plusieurs niveaux de métallisation. Pour ce type de circuits, elle se révèle être l'une des méthodes qui utilise le plus petit nombre d'inconnues pour atteindre une précision donnée. Néanmoins, la WCIP n'est pas adaptée à la prise en compte des diélectriques inhomogènes et des trous d'interconnexion. L'objectif de la thèse est de s'affranchir de ces limitations par un couplage avec des méthodes numériques volumiques. En premier lieu, l'hybridation a été mise en œuvre avec une méthode basée sur la théorie des lignes de transmission pour des raisons de correspondance de maillages. Par la suite, le couplage avec une technique d'éléments finis de type Galerkin Discontinu (notée GD) Hybridée permet d'atteindre des objectifs de précision et de rapidité car GD apporte une flexibilité dans la discrétisation. En effet, c'est une méthode d'éléments finis non conforme qui permet notamment de faire varier d'un élément à l'autre l'ordre polynomial d'approximation. On a ainsi développé une nouvelle méthode numérique hybride couplant la WCIP avec des méthodes volumiques qui offrent plus de souplesse pour la prise en compte des milieux complexes. Enfin, une stratégie de résolution par décomposition de domaines est également abordée à la fin du manuscrit. / Electromagnetic susceptibility study of electronic circuits requires the use of a simulation tool which is fast, accurate and flexible enough to incorporate the latest technological innovations. The Wave Concept Iterative Procedure (WCIP) initially proposed by H. Baudrand is particularly adapted for numerical modeling of multilayered circuits with multilevel metallization. For this kind of circuits, it turns out to be one of the methods that uses the smallest number of unknowns to reach a given accuracy. However, the WCIP is not appropriate for inhomogeneous dielectric substrates and metallized via holes. The aim of this PhD thesis is to overcome these limitations coupling the WCIP with volume numerical methods. First, hybridization is carried out with the Frequency Domain Transmission Line Matrix (denoted FDTLM) assuming matching meshes at the interface between computational domains of both methods. Subsequently, the coupling with a finite element technique like a Hybridized Discontinuous Galerkin (denoted DG) method is considered to achieve the objectives of accuracy and speed because DG brings flexibility in the discretization. Indeed, it is a nonconforming finite element method which allows in particular changing the polynomial approximation order from one element to another. Therefore, a new hybrid method is developed coupling the WCIP with volume numerical methods which offer more flexibility for dealing with complex environments. Finally, a domain decomposition solution strategy is also discussed at the end of the manuscript.

Cette thèse se penche sur la propagation d’ondes au sein du système coupléTerre-océan-atmosphère. La compréhension de ces phénomènes a une importance majeure pour l’étude de perturbations sismiques et d’explosions atmosphériques notamment dans le cadre de missions spatiales planétaires. Les formes d’ondes issues du couplage fluide-solide permettent d’obtenir de précieuses informations sur la source du signal ou les propriétés des milieux de propagation. On développe donc deux outils numériques d’ordre élevé pour la propagation d’ondes acoustiques et de gravité. L'u en différences finies et se concentre sur le milieu atmosphérique et la propagation d’ondes linéaires dans un milieu stratifié visqueux et avec du vent. Cette méthode linéaire est validée par des solutions quasi-analytiques reposant sur les équations de dispersion dans une atmosphère stratifiée. Elle est aussi appliquée à deux cas d’études : la propagation d’ondes liée à l’impact d’une météorite à la surface de Mars (mission NASA INSIGHT), et la propagation d’ondes atmosphériques liées au tsunami de Sumatra en 2004. La seconde méthode résout la propagation non-linéaire d’ondes gravito-acoustiques dans une atmosphère couplée, avec topographie, à la propagation d’ondes élastiques dans un solide visco-élastique. Cette méthode repose sur sur le couplage d’une formulation en éléments finis discontinus, pour résoudre les équations de Navier-Stokes la partie fluide, par éléments finis continus pour résoudre les équations de l’élastodynamique dans la partie solide. Elle a été validée grâce à des solutions analytiques ainsi que par des comparaisons avec les résultats de la méthode par différences finies. / This thesis deals with the wave propagation problem within the Earth-ocean-atmosphere coupled system. A good understanding of the these phenomena has a major importance for seismic and atmospheric explosion studies, especially for planetary missions. Atmospheric wave-forms generated by explosions or surface oscillations can bring valuable information about the source mechanism or the properties of the various propagation media. We develop two new numerical full-wave high-order modeling tools to model the propagation of acoustic and gravity waves in realistic atmospheres. The first one relies on a high-order staggered finite difference method and focus only on the atmosphere. It enables the simultaneous propagation of linear acoustic and gravity waves in stratified viscous and windy atmosphere. This method is validated against quasi-analytical solutions based on the dispersion equations for a stratified atmosphere. It has also been employed to investigate two cases : the atmospheric propagation generated by a meteor impact on Mars for the INSIGHT NASA mission and for the study of tsunami-induced acoutic and gravity waves following the 2004 Sumatra tsunami. The second numerical method resolves the non-linear acoustic and gravity wave propagation in a realistic atmosphere coupled, with topography, to the elastic wave propagation in a visco-elastic solid. This numerical tool relies on a discontinuous Galerkin method to solve the full Navier-Stokes equations in the fluid domain and a continuous Galerkin method to solve the elastodynamics equations in the solid domain. It is validated against analytical solutions and numerical results provided by the finite-difference method.

Neste trabalho, é apresentada uma formulação apropriada à análise de guias de ondas eletromagnéticos, cobrindo do espectro de microondas até o da óptica. Nas regiões a partir do ultravioleta, os comprimentos de onda são equivalentes às dimensões atômicas e a formulação necessita de uma abordagem quântica, que não é considerada neste estudo. A formulação é fundamentada nos métodos vetorial magnético e dos elementos finitos (MEF), em meios não homogêneos, anisotrópicos e não dissipativos, embora a dissipação possa ser facilmente introduzida na análise. Deu-se preferência à formulação com o campo magnético em vez do elétrico, pelo fato do campo magnético ignorar descontinuidades elétricas. Ele é contínuo em regiões de permeabilidade homogênea, propriedade dos meios dielétricos em geral ( = 0), independente da permissividade dos respectivos meios, conquanto os campos elétricos sejam descontínuos entre regiões de permissividades diferentes. / This work presents a suitable formulation to the analysis of electromagnetic waveguide, covering the spectrum of the microwave to optics. In regions from the ultraviolet, the wavelengths are equivalent to atomic dimensions and the design needs a quantum approach, which is not considered in this study. The formulation is based on the magnetic vector and the finite element methods (FEM), in non-homogeneous, anisotropic and non-dissipative dielectric materials, while the dissipation can be easily introduced in the analysis. Preference was given to the formulation with the magnetic field, because the magnetic field ignores electrical discontinuities. It is continuous in regions of homogeneous permeability, property of all dielectric materials (=0), independent of the permissiveness of respective regions, while electric fields are discontinuous between regions of different permittivities.

Made available in DSpace on 2015-05-15T11:46:21Z (GMT). No. of bitstreams: 1 arquivototal.pdf: 1219498 bytes, checksum: d931fbd819e737ae31a593f5c64c55e0 (MD5) Previous issue date: 2014-09-26 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / In this work we study the Navier-Stokes equations in bounded domains of Rn. Initially we the case n = 2 and we show that its variational formulation is well put (in case the Hadamard). We show the existence of solution for the case n  4 . In both cases we use the Faedo-Galerkin method. / Neste trabalho estudamos as equações de Navier-Stokes em domínios limitados do Rn. Inicialmente consideramos o caso n = 2emostramos que sua formulação variacional est´a bem posta (no sentido de Hadamard). Em seguida, mostramos a existência de solução para o caso n  4. Em ambos os casos utilizamos o método de Faedo-Galerkin.

Este trabalho apresenta uma metodologia de análise da heterogeneidade do concreto a partir de modelos computacionais desenvolvidos com o método de Galerkin livre de elementos. Esse método se caracteriza pela discretização de um domínio de interesse por um conjunto de partículas sem que exista explicitamente uma malha de elementos no sentido convencional. O objetivo é a previsão das propriedades mecânicas macroscópicas do material resultante a partir das fases individuais e do arranjo geométrico. O concreto foi admitido, na escala mesoscópica, como um composto formado por inclusões (agregado graúdo) imersas em uma matriz (argamassa). Para a simulação foi desenvolvida uma formulação multiregiões onde se admitiu que cada agregado e a argamassa são domínios distintos interligados nas suas interfaces. Para isto foram utilizadas técnicas de subdivisões do domínio (elemento representativo) ao ponto que os seus comportamentos mecânicos não foram comprometidos. Para simular o processo das perdas de rigidez com a formação da fissuração no concreto foi admitido o efeito da mecânica do dano contínuo através do modelo de Mazars. Para as análises foram desenvolvidos modelos computacionais bidimensionais e tridimensionais da heterogeneidade do concreto. A geometria dos agregados foi aproximada por circunferências e elipses no caso 2D e por esferas e elipsoides no caso 3D. Como conclusão a metodologia de multiregiões com o método de Galerkin livre de elementos foi satisfatória e os modelos apresentaram caminhos preferenciais de ruptura adequados durante a evolução da danificação. / This thesis presents a methodology for analyzing the heterogeneity of concrete from computational models developed with the element free Galerkin method. This method is characterized by discretization of a domain of interest by a set of particles with no explicit mesh in the conventional sense. The goal is to predict the macroscopic mechanical properties of the material resulting from the individual phases and the geometric arrangement. The concrete was assumed, in the mesoscopic scale, as a compound formed by inclusions (coarse aggregate) embedded in a matrix (mortar). For the simulation, a formulation was developed where multi regions were admitted, assuming that each aggregate and mortar are distinct domains connected by their interfaces. For this we used techniques of subdivisions of the domain (representative elements) to the point that their mechanical behaviors were not compromised. To simulate the process of loss of stiffness with the formation of cracks in the concrete, continuum damage mechanics was admitted through Mazars model. For the analysis, two-dimensional and three-dimensional computer models of the heterogeneity of the concrete were developed. The shape of the aggregates was approximated by circles and ellipses in the two-dimensional case, and by spheres and ellipsoids for the 3D problems. In conclusion the multi region methodology with the element free Galerkin methods was satisfactory and the models presented suitable preferred paths for the rupture during the evolution of damage.

In this work we study the existence and uniqueness of solutions of a boundary value problem for equations of non-homogeneous wave. We will use the Faedo - Galerkin method to ensure the existence of solutions and also prove the exponential decay of the solution. / Neste trabalho estudaremos a existência e unicidade de soluções de um problema de valor de fronteira para equações da onda não homogênea. Usaremos o método de Faedo-Galerkin para garantir a existência de soluções e também provaremos o decaimento exponencial da solução do problema.

Submitted by Erika Demachki (erikademachki@gmail.com) on 2017-03-09T21:07:51Z No. of bitstreams: 2 Dissertação - Steffânio Moreno de Sousa - 2017.pdf: 1402439 bytes, checksum: 6a4985baeae0454d7088588c3c87e295 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Approved for entry into archive by Luciana Ferreira (lucgeral@gmail.com) on 2017-03-10T10:48:27Z (GMT) No. of bitstreams: 2 Dissertação - Steffânio Moreno de Sousa - 2017.pdf: 1402439 bytes, checksum: 6a4985baeae0454d7088588c3c87e295 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Made available in DSpace on 2017-03-10T10:48:28Z (GMT). No. of bitstreams: 2 Dissertação - Steffânio Moreno de Sousa - 2017.pdf: 1402439 bytes, checksum: 6a4985baeae0454d7088588c3c87e295 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Previous issue date: 2017-03-03 / Conselho Nacional de Pesquisa e Desenvolvimento Científico e Tecnológico - CNPq / In this work we study existence and regularity of non-negative solution of elliptical systems of the type 8>< >: 􀀀Dpu = f (x;u;v); 􀀀Dqv = g(x;u;v) in W; u;v > 0 in W; where 1 < p;q <N, W IRN is a bounded domain with smooth boundary ¶W, and f ;g are of the type singular-convex or W=IRN and f ;g are concave-convex. In case W we will find solutions that cancel out ¶W while in the case W = IRN solutions in C1(IRN)\L¥(IRN).We will use the Galerkin method and the comparison principle. In case W = IRN we will use the method of sub and super solutions, variational methods and principles of maximum. / Neste trabalho estudaremos existência e regularidade de soluções positivas de sistemas elípticos do tipo 8>< >: 􀀀Dpu = f (x;u;v); 􀀀Dqv = g(x;u;v) in W u;v > 0 in W; onde 1 < p;q < N, W IRN é um domínio limitado com fronteira ¶W regular, e f , g são do tipo convexo-singular ou W = IRN e f , g são do tipo côncavo-convexo. No caso W limitado encontraremos soluções que se anulam em ¶W, enquanto que, no caso W = IRN as soluções em C1(IRN) \ L¥(IRN). No caso f , g singulares utilizaremos o método de Galerkin e princípio de comparação. No caso W = IRN utilizaremos o método de sub e super-soluções, métodos variacionais e princípios de máximo.