Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Mecânica, 2018. / Em busca de melhor aproveitamento da energia eólica disponível, a demanda em relação a projetos de torres eólicas é por torres cada vez mais altas e pás maiores. Pás de torres eólicas são estruturas flexíveis que estão sujeitas a altas velocidades do vento e restritas pela movimentação da torre. As principais forças aplicadas no sistema são do tipo inercial, gravitacional e aerodinâmica, que são responsáveis pelas vibrações da pá. Além disso, o acoplamento das vibrações da pá e da torre leva à um sistema com elevado grau de complexidade. A torre, em função de suas dimensões, pode ser modelada como uma viga esbelta e pode estar sujeita a vibrações excessivas, que podem causar acidentes. A primeira frequência natural da torre está contida na faixa de atuação das pás, então um estudo do seu acoplamento com as pás é necessário. Este trabalho foca no estudo da resposta dinâmica da interação torre-pá através de modelos dinâmicos discretos. O modelo discreto considera efeitos de peso-próprio e de força centrífuga devido à rotação das pás. Evidencia-se a relação entre posição do rotor e frequência de resposta do sistema e, também, entre velocidade de rotação e frequência de resposta do sistema. / To better exploit available wind energy, the demand for wind tower design is for increasingly higher towers and larger blades. Wind tower blades are flexible structures that are subject to high forcing velocities and restricted by tower movement. The main forces applied in the system are inertial, gravitational and aerodynamic, which are responsible for blade vibrations. In addition, coupling vibrations of the blade and the tower leads to a system with a high degree of complexity. The tower, depending on its size, is considered a slender beam and may be subject to excessive vibrations, which may cause accidents. The first natural frequency of the tower is contained in the operating range of the blades, so a study of their coupling with the blades is necessary. This work presents a study on dynamic response of tower-blade interaction through discrete dynamical models. The discrete model considers self-weight effects and centrifugal force due to blade rotation. It shows the relationship between rotor position and system frequency response and also between rotation speed and frequency response of the system.

The time-dependent elastic wave equation is the foundation of seismology. It is very useful in studying the wave propagation in elastic solids. Simulation of Rayleigh waves, which is governed by the equation, requires high accuracy solutions. Finite difference method have been widely used for the simulation of Rayleigh waves. However, it is not obvious how to effectively impose the free surface boundary condition on a curved surface. On the other hand, discontinuous Galerkin methods are more flexible in handling complex geometries. / In this thesis, we develop a class of discontinuous Galerkin methods for time-dependentelastic wave equation in isotropic homogeneous medium. This method is explicit, locally and globally energy conserving. Also, the L² convergence of this method is optimal and the convergence in energy norm is independent of Poisson's ratio. / Besides, we apply our method to simulate Rayleigh waves on curved free surfaces. We also impose a perfectly matched layer to absorb the outward waves. Numerical examples show that our method can accurately capture features of Rayleigh waves even in a domain with high Poisson's ratio. / 時間依賴型彈性波動方程」是地震學的基礎。這組方程對於波在彈性固體中傳播的研究非常有用。雷利波是由這個方程所控制。模擬雷利波須要有非常準確的解。有限差分法廣泛地應用在雷利波的模擬上，可是如何有效地施加自由表面邊界條件於曲面上的方法並不明顯。另一方面，間斷伽遼金方法能更靈活地處理複雜的幾何形狀。 / 在這篇論文中，我們發展了一類間斷伽遼金方法去求「在均勻各向同性的介質上的時間依賴型彈性波動方程」的解。我們將表明，這種方法是顯式的，局部及全域能量守恆的，而它的收斂是最優的和獨立於泊松比的。 / 除此之外，我們運用這個方法來模擬雷利波在具有起伏的自由表面的半空間模型的傳播。我們會使用完美匹配層去吸收朝外的波動。數值算例反映，即使在高柏松比的介質中，這個方法也可以準確地捕捉雷利波的特徵。 / Lam, Chi Yeung. / Thesis M.Phil. Chinese University of Hong Kong 2014. / Includes bibliographical references (leaves 44-47). / Abstracts also in Chinese. / Title from PDF title page (viewed on 06, October, 2016). / Detailed summary in vernacular field only. / Detailed summary in vernacular field only. / Detailed summary in vernacular field only.

This doctoral thesis is concerned with the mathematical modeling of magnetoelastic materials and the analysis of PDE systems describing these materials and obtained from a variational approach. The purpose is to capture the behavior of elastic particles that are not only magnetic but exhibit a magnetic domain structure which is well described by the micromagnetic energy and the Landau-Lifshitz-Gilbert equation of the magnetization. The equation of motion for the material’s velocity is derived in a continuum mechanical setting from an energy ansatz. In the modeling process, the focus is on the interplay between Lagrangian and Eulerian coordinate systems to combine elasticity and magnetism in one model without the assumption of small deformations. The resulting general PDE system is simplified using special assumptions. Existence of weak solutions is proved for two variants of the PDE system, one including gradient flow dynamics on the magnetization, and the other featuring the Landau-Lifshitz-Gilbert equation. The proof is based on a Galerkin method and a fixed point argument. The analysis of the PDE system with the Landau-Lifshitz-Gilbert equation uses a more involved approach to obtain weak solutions based on G. Carbou and P. Fabrie 2001. / Die vorliegende Doktorarbeit beschäftigt sich mit der mathematischen Modellierung magnetoelastischer Materialien und der Analysis von Systemen partieller Differentialgleichungen für diese Materialien. Die Herleitung der partiellen Differentialgleichungen erfolgt mittels eines variationellen Zugangs. Ziel ist es, das Verhalten elastischer Teilchen zu beschreiben, welche nicht nur magnetisch sind, sondern sich durch eine magnetische Domänenstruktur auszeichnen. Diese Struktur wird beschrieben durch die mikromagnetische Energie und die Landau-Lifshitz-Gilbert Gleichung der Magnetisierung. Die Bewegungsgleichung für die Geschwindigkeit des Materials ist in einem kontinuumsmechanischen Setting von einer Energiegleichung abgeleitet. In der Modellierung liegt der Fokus auf dem Zusammenspiel von Lagrange’schen und Euler’schen Koordinaten, um Elastizität und Magnetismus in einem Modell zu kombinieren. Dies geschieht ohne die Annahme kleiner Deformationen. Das resultierende allgemeine System partieller Differentialgleichungen wird durch spezielle Annahmen vereinfacht und es wird die Existenz von schwachen Lösungen gezeigt. Der Beweis wird für zwei Varianten des Differentialgleichungssystems geführt. Das erste System enthält die Beschreibung der Dynamik der Magnetisierung mittels Gradientenfluss, im zweiten wird die Dynamik mittels Landau-Lifshitz-Gilbert Gleichung beschrieben. Schlüsselidee des Beweises ist ein Galerkin-Ansatz, kombiniert mit einem Fixpunkt-Argument. Zum Beweis der Existenz schwacher Lösungen des Systems mit Landau-Lifshitz-Gilbert Gleichung wird eine aufwändigere Methode herangezogen, welche auf einer Arbeit von G. Carbou und P. Fabrie aus 2001 beruht.

The present thesis considers the development and analysis of arbitrary Lagrangian-Eulerian discontinuous Galerkin (ALE-DG) methods with time-dependent approximation spaces for conservation laws and the Hamilton-Jacobi equations. Fundamentals about conservation laws, Hamilton-Jacobi equations and discontinuous Galerkin methods are presented. In particular, issues in the development of discontinuous Galerkin (DG) methods for the Hamilton-Jacobi equations are discussed. The development of the ALE-DG methods based on the assumption that the distribution of the grid points is explicitly given for an upcoming time level. This assumption allows to construct a time-dependent local affine linear mapping to a reference cell and a time-dependent finite element test function space. In addition, a version of Reynolds’ transport theorem can be proven. For the fully-discrete ALE-DG method for nonlinear scalar conservation laws the geometric conservation law and a local maximum principle are proven. Furthermore, conditions for slope limiters are stated. These conditions ensure the total variation stability of the method. In addition, entropy stability is discussed. For the corresponding semi-discrete ALE-DG method, error estimates are proven. If a piecewise $\mathcal{P}^{k}$ polynomial approximation space is used on the reference cell, the sub-optimal $\left(k+\frac{1}{2}\right)$ convergence for monotone fuxes and the optimal $(k+1)$ convergence for an upwind flux are proven in the $\mathrm{L}^{2}$-norm. The capability of the method is shown by numerical examples for nonlinear conservation laws. Likewise, for the semi-discrete ALE-DG method for nonlinear Hamilton-Jacobi equations, error estimates are proven. In the one dimensional case the optimal $\left(k+1\right)$ convergence and in the two dimensional case the sub-optimal $\left(k+\frac{1}{2}\right)$ convergence are proven in the $\mathrm{L}^{2}$-norm, if a piecewise $\mathcal{P}^{k}$ polynomial approximation space is used on the reference cell. For the fullydiscrete method, the geometric conservation is proven and for the piecewise constant forward Euler step the convergence of the method to the unique physical relevant solution is discussed. / Die vorliegende Arbeit beschäftigt sich mit der Entwicklung und Analyse von arbitrar Lagrangian-Eulerian discontinuous Galerkin (ALE-DG) Methoden mit zeitabhängigen Testfunktionen Räumen für Erhaltungs- und Hamilton-Jacobi Gleichungen. Grundlagen über Erhaltungsgleichungen, Hamilton-Jacobi Gleichungen und discontinuous Galerkin Methoden werden präsentiert. Insbesondere werden Probleme bei der Entwicklung von discontinuous Galerkin Methoden für die Hamilton-Jacobi Gleichungen untersucht. Die Entwicklung der ALE-DG Methode basiert auf der Annahme, dass die Verteilung der Gitterpunkte zu einem kommenden Zeitpunkt explizit gegeben ist. Diese Annahme ermöglicht die Konstruktion einer zeitabhängigen lokal affin-linearen Abbildung auf eine Referenzzelle und eines zeitabhängigen Testfunktionen Raums. Zusätzlich kann eine Version des Reynolds’schen Transportsatzes gezeigt werden. Für die vollständig diskretisierte ALE-DG Methode für nichtlineare Erhaltungsgleichungen werden der geometrischen Erhaltungssatz und ein lokales Maximumprinzip bewiesen. Des Weiteren werden Bedingungen für Limiter angegeben. Diese Bedingungen sichern die Stabilität der Methode im Sinne der totalen Variation. Zusätzlich wird die Entropie-Stabilität der Methode diskutiert. Für die zugehörige semi-diskretisierte ALE-DG Methode werden Fehlerabschätzungen gezeigt. Wenn auf der Referenzzelle ein Testfunktionen Raum, der stückweise Polynome vom Grad $k$ enthält verwendet wird, kann für einen monotonen Fluss die suboptimale Konvergenzordnung $\left(k+\frac{1}{2}\right)$ und für einen upwind Fluss die optimale Konvergenzordnung $\left(k+1\right)$ in der $\mathrm{L}^{2}$-Norm gezeigt werden. Die Leistungsfähigkeit der Methode wird anhand numerischer Beispiele für nichtlineare Erhaltungsgleichungen untersucht. Ebenso werden für die semi-diskretisierte ALE-DG Methode für nichtlineare Hamilton-Jacobi Gleichungen Fehlerabschätzungen gezeigt. Wenn auf der Referenzzelle ein Testfunktionen Raum, der stückweise Polynome vom Grad k enthält verwendet wird, kann im eindimensionalen Fall die optimale Konvergenzordnung $\left(k+1\right)$ und im zweidimensionalen Fall die suboptimale Konvergenzordnung $\left(k+\frac{1}{2}\right)$ in der $\mathrm{L}^{2}$-Norm gezeigt werden. Für die vollständig diskretisierte ALE-DG Methode werden der geometrischen Erhaltungssatz bewiesen und für die stückweise konstante explizite Euler Diskretisierung wird die Konvergenz gegen die eindeutige physikalisch relevante Lösung diskutiert.

Dans cette thèse, nous nous intéressons à un problème mathématique issu de la modélisation de taux d'érosion maximale dans la stratigraphie géologique. Une contrainte globale sur $\partial_t u$, la dérivée par rapport au temps de la solution, est la principale caractéristique de ce modèle. Ce qui nous amène à considérer une équation non linéaire pseudo-parabolique avec un coefficient de diffusion qui est une fonction non-linéaire de $\partial_t u$. En outre, le problème dégénère de telle sorte de tenir compte implicitement de la contrainte. Nous présentons un résultat de l'existence d'une solution au problème continu. Ensuite, une méthode DgFem (discontinuous Galerkin finite element method) pour son approximation numérique est développée. Notre objectif est d'utiliser les propriétéess d'approximation constante par morceaux pour tenir compte implicitement de la contrainte.

L'objet de cette thèse est d'étudier un problème de la magnétostatique tridimensionnel. On propose trois formulations mixtes couplant une méthode d'éléments finis pour tenir compte du milieu hétérogène et une méthode éléments de frontière pour le milieu extérieur homogène. Pour la méthode intégrale on a utilisé les équations de Calderon, l'opérateur de Neumann-Dirichlet ou d'autres opérateurs intégraux. L'utilisation des éléments d'arête de Nédélec pour le champ magnétique, et les éléments de face de Raviart pour l'induction magnétique permet d'utiliser des méthodes éléments finis conformes. Des résultats numériques ont permis de valider ces méthodes. La deuxième partie a porté sur la comparaison de diverses discrétisations pour l'opérateur de Poincaré-Steklov. Ces méthodes ont été comparées sur une formulation de la magnétostatique. Enfin, on propose des formulations discontinues du problème de la magnétostatique avec des conditions aux limites. On montre que ces formulations sont consistantes et des estimations d'erreur sont obtenues.

Ce travail de thèse a pour objet la mise en oeuvre d'une méthode de calcul des états liés dans les structures à multipuits quantiques. Il participe ainsi à l'amélioration des outils de simulation permettant d'optimiser les composants avant leur réalisation. Nous présentons le modèle physique utilisé ainsi que les différentes méthodes couramment employées pour le calcul de ces états. Une comparaison avec le calcul par éléments finis du premier ordre montre un avantage majeur de notre approche dans des cas limites usuels comme le couplage de deux puits identiques ou le calcul des bandes de valence d'un puits quantique large, ainsi qu'en terme de rapidité. La finalité de ce calcul est l'évaluation du gain matériau, élément de base de la simulation des composants. Ce nouveau module vient compléter le simulateur BCBV dont nous rappelons les principaux modèles tels que celui de dérive-diffusion et du couplage électrique-optique en semi-classique. Cependant, la présence de zones quantiques peut nécessiter une approche par la matrice de densité pour rendre compte, de manière plus précise, des phénomènes de transport. Pour finir, nous tentons de comparer les résultats de la simulation du gain avec des mesures effectuées à partir de lasers de type Fabry-Pérot.

The subject of this thesis is the development of an accurate, general and robust numerical method capable of predicting the flow behavior of two-phase immiscible fluids, separated by a well defined interface. In the quest of an accurate and robust numerical method for the modeling of two-phase flows, one has to keep in mind the intrinsic properties and difficulties associated with the problem: (i) those flows are mostly three-dimensional, (ii) some flows are steady, others unsteady, (iii) the interface might encounter a lot of topology changes (like merger or break-up), (iv) large jumps of density and viscosity might exist across the interface (e.g. ratio of density of 1/1000 for water and air), (v) surface tension forces may play a very important role in the interface dynamics. Hence, the influence of this force should be accurately evaluated and incorporated into the model, (vi) mass conservation is of primary importance. All these issues are addressed in this thesis, and special techniques are proposed for their treatment, which enables to construct the desired computational method. The chosen computational method combines a pressure stabilized finite element method for the Navier Stokes equations with a discontinuous Galerkin (DG) method for the level set equation. Such a combination of those two numerical methods results in a simple and effective algorithm that allows to simulate diverse flow regimes presenting large density and viscosity ratios (ratio up to 1/1000).

No description available.

In this thesis, we analyze theoretical properties of the discontinuous Galerkin method (DGM) and propose novel approaches to implementation with the aim to increase its efficiency. First, we derive explicit expressions for the eigenvalues (spectrum) of the discontinuous Galerkin spatial discretization applied to the linear advection equation. We show that the eigenvalues are related to the subdiagonal [p/p+1] Pade approximation of exp(-z) when the p-th degree basis functions are used. Then, we extend the analysis to nonuniform meshes where both the size of elements and the composition of the mesh influence the spectrum. We show that the spectrum depends on the ratio of the size of the largest to the smallest cell as well as the number of cells of different types. We find that the spectrum grows linearly as a function of the proportion of small cells present in the mesh when the size of small cells is greater than some critical value. When the smallest cells are smaller than this critical value, the corresponding eigenvalues lie outside of the main spectral curve. Numerical examples on nonuniform meshes are presented to show the improvement on the time step restriction. In particular, this result can be used to improve the time step restriction on Cartesian grids. Finally, we present a discontinuous Galerkin method for solutions of the Euler equations on Cartesian grids with embedded geometries. Cutting an embedded geometry out of the Cartesian grid creates cut cells, which are difficult to deal with for two reasons. One is the restrictive CFL number and the other is the integration on irregularly shaped cells. We use explicit time integration employing cell merging to avoid restrictively small time steps. We provide an algorithm for splitting complex cells into triangles and use standard quadrature rules on these for numerical integration. To avoid the loss of accuracy due to straight sided grids, we employ the curvature boundary conditions. We show that the proposed method is robust and high-order accurate.