Cette thèse propose l’étude d’une méthode numérique permettant de simuler un plasma. On considère un ensemble de particules, dont le mouvement est régi par l’équation de Vlasov, et qui est sensible aux forces électromagnétiques, qui proviennent des équations de Maxwell. La résolution numérique des équations de Vlasov-Maxwell est réalisée par une méthode Particle In Cell (PIC). La résolution des équations de Maxwell nécessite un maillage suffisamment fin afin de modéliser correctement les problémes multi-échelles que nous souhaitons traiter. Cependant, mailler finement tout le domaine de calcul a un coût. La nouveauté de cette thèse est de proposer un solveur PIC sur des maillages cartésiens localement raffinés, des maillages non conformes, afin de garantir la bonne modélisation du phénomène physique et d’éviter une trop forte pénalisation des temps de calcul.Nous utilisons une méthode Galerkin Discontinue en domaine temporelle (GDDT), qui offre l’avantage d’être d'une grande flexibilité dans le choix du maillage et qui est une méthode d’ordre élevé. Un point fondamental dans l’étude des solveurs PIC concerne le respect de la conservation de la charge. Nous proposons deux approches afin de traiter cet aspect. La première concerne les méthodes utilisant un système de Maxwell augmenté, dont la nouveauté a été de les étendre aux maillages non conformes. La seconde approche repose sur une méthode originale de pré-traitement du calcul du terme source de courant. / This thesis deals with the study of a numerical method to simulate a plasma. We consider a set of particles whose displacement is governed by the Vlasov equation and which creates an electromagnetic field thanks to Maxwell equations. The numerical resolution of the Vlasov-Maxwell system is performed by a Particle In Cell (PIC) method. The resolution of Maxwell equations needs a sufficiently fine mesh to correctly simulate the multi scaled problems that we have to face. Yet, a uniform fine mesh of the whole domain has a prohibitive cost. The novelty of this thesis is a PIC solver on locally refined Cartesian meshes : non conforming meshes, to guarantee the good modeling of the physical phenomena and to avoid too large CPU time. We use the Discontinuous Galerkin in Time Domain (DGTD) method which has the advantage of a great flexibility in the choice of the mesh and which is a high order method. A fundamental point in the study of PIC solvers is the respect of the charge conserving law. We propose two approaches to tackle this point. The first one deals with augmented Maxwell systems, that we have adapted to non conforming meshes. The second one deals with an original method of preprocessing of the calculation of the current source term.

Le travail présenté dans ce manuscrit est une contribution à l’étude des vibrations non-linéaires des poutres isotropes et en composite, en mouvement de rotation. Le modèle mathématique utilisé est basé sur la formulation intrinsèque et géométriquement exacte de Hodges, dédiée au traitement des poutres ayant des grands déplacements et de petites déformations. La résolution est faite dans le domaine fréquentiel suite à une discrétisation spatio-temporelle, en utilisant l’approximation de Galerkin et la méthode de l’équilibrage harmonique, avec des conditions aux limites correspondantes aux poutres encastrées-libres. Le systéme dynamique final est traité par des méthodes de continuation : la méthode asymptotique numérique et la méthode pseudo-longueur d’arc. Des algorithmes basés sur ces méthodes de continuation ont été développés et une étude comparative de convergence a été menée. Cette étude a cerné les aspects : statique, analyse modale linéaire, vibrations libres non-linéaires et les vibrations forcées non-linéaires des poutres rotatives. Ces algorithmes de continuations ont été testés pour le calculs des courbes de réponse sur des cas traités dans la littérature. La résonance interne et la stabilité des solutions obtenues sont étudiées / The work presented in this manuscript is a contribution to the non-linear vibrations of the isotropic beams and composite rotating beams study. The mathematical model used is based on the intrinsic formulation and geometrically exact of Hodges, developped for beams subjected to large displacements and small deformations. The resolution is done in the frequency domain after a spatial-temporal dicretisation, by using the Galerkin approximation and the the harmonic balance method, with boundary conditions corresponding to the clamped-free. The final dynamic system is treated by continuation methods : asymptotic numerical method and the pseudo-arc length method, whose algorithms based on these continuation methods were developed and a convergence study was carried out. This study surround the aspects : statics, linear modal analysis, non-linear free vibrations and the non-linear forced vibrations of the rotating beams. These continuation algorithms were tested for the response curves calculations on cases elaborated in the literature. Internal resonance and the stability of the solutions obtained are studied

Model order reduction is necessary even in a time where the parallel processing is usual in almost any personal computer. The recent Model Reduction Methods are useful tools nowadays on reducing the problem processing. This work intends to describe a combination between POD (Proper Orthogonal Decomposition) and Ritz vectors that achieve an efficient Galerkin projection that changes during the processing, comparing the development of the error and the convergence rate between the full space and the projection space, in addition to check the stability of the projection space, leading to an adaptive model order reduction for nonlinear dynamical problems more efficient. This model reduction is supported by a secant formulation, which is updated by BFGS (Broyden - Fletcher - Goldfarb - Shanno) method to accelerate convergence of the model, and a tangent formulation to correct the projection space. Furthermore, this research shows that this method permits a correction of the reduced model at low cost, especially when the classical POD is no more efficient to represent accurately the solution. / A Redução de ordem de modelo é necessária, mesmo em uma época onde o processamento paralelo é usado em praticamente qualquer computador pessoal. Os recentes métodos de redução de modelo são ferramentas úteis nos dias de hoje para a redução de processamento de um problema. Este trabalho pretende descrever uma combinação entre POD (Proper Orthogonal Decomposition) e vetores de Ritz para uma projecção de Galerkin eficiente que sofre alterações durante o processamento, comparando o desenvolvimento do erro e a taxa de convergência entre o espaço total e o espaço de projeção, além da verificação de estabilidade do espaço de projeção, levando a uma redução de ordem do modelo adaptativo mais eficiente para problemas dinâmicos não-lineares. Esta redução de modelo é assistida por uma formulação secante, que é atualizado pela formula de BFGS (Broyden - Fletcher- Goldfarb - Shanno) com o intuito de acelerar a convergência do modelo, e uma formulação tangente para a correção do espaço de projeção. Além disso, esta pesquisa mostra que este método permite a correção do modelo reduzido com baixo custo, especialmente quando o clássico POD não é mais eficiente para representar com precisão a solução.

Dans cette thèse, nous nous intéressons à l’amélioration des performances numériques dela méthode Galerkin Discontinu en Domaine Temporel (GDDT), afin de valoriser son emploi industrielpour des problèmes de propagation d’ondes électromagnétiques. Pour ce faire, nous cherchons à réduire lenombre d’éléments des maillages utilisés en appliquant une stratégie de h-déraffinement/p-enrichissement.Dans un premier temps, nous montrons que si ce type de stratégie permet d’améliorer significativementl’efficacité numérique des résolutions dans un cadre conforme, son extension aux maillages non-conformespeut s’accompagner de contre-performances rédhibitoires limitant fortement leur intérêt pratique. Aprèsavoir identifié que ces dernières sont causées par le traitement des termes de flux non-conformes, nousproposons une méthode originale de condensation afin de retrouver des performances avantageuses. Cellecise base sur une redéfinition des flux non-conformes à partir d’un opérateur de reconstruction de traces,permettant de recréer une conformité d’espaces, et d’un produit scalaire condensé, assurant un calculapproché efficace. La stabilité et la consistance du schéma GDDT ainsi défini sont établies sous certainesconditions portant sur ces deux quantités. Dans un deuxième temps, nous détaillons la construction desopérateurs de trace et des produits scalaires associés. Nous proposons alors des flux condensés pourplusieurs configurations non-conformes, et validons numériquement la convergence du schéma GDDT résultant.Puis, nous cherchons à concevoir un algorithme de h-déraffinement/p-enrichissement automatisé,dans le but de générer des maillages hp minimisant les coûts de calcul du schéma. Ce processus est traduitsous la forme d’un problème d’optimisation combinatoire sous plusieurs contraintes de natures trèsdiverses. Nous présentons alors un algorithme de post-maillage basé sur un parcours efficace de l’arbrede recherche des configurations admissibles, associé à un processus de déraffinement hiérarchique. Enfin,nous mettons en œuvre la chaîne de calcul développée sur plusieurs cas-tests d’intérêt industriel, etévaluons son apport en termes de performances numériques. / This thesis is devoted to improving the numerical efficiency of the Discontinuous Galerkinin Time Domain (DGDT) method, in order to enhance its suitability for industrial use. One can noticethat, in an hp-conforming context, increasing correlatively the approximation order and the mesh sizeis a powerful strategy to reduce numerical costs. However, in complex geometries, the mesh can beconstrainted by the presence of small-scale inner elements, leading to hp-nonconforming configurationswith hanging nodes. The first issue we are dealing with is related to the nonconforming fluxes involvedin these configurations, whose high computational costs can deter the use of hp-coarsening strategies.In order to recover a satisfactory performance level, an original flux-lumping technique is set up. Thistechnique relies on recasting hybrid fluxes into conforming ones, and is performed by introducing twoingredients : a reconstruction operator designed to map traces from each side of a nonconforming interfaceinto the same functional space, and a lumped scalar product granting efficient integral computations.The resulting DGTD scheme is then proved to be stable and consistent, under some assumptions on thelatter two elements. Subsequently, we develop a lumped flux construction routine, and show numericalconvergence results on basic hybrid configurations. In a second part, we implement an automated strategyaiming at generating efficient hp-nonconforming meshes, well-suited to the previous DGDT scheme. To doso, a post-meshing process is formalized into a constrained optimization problem. We then put forward aheuristic hp-coarsening algorithm, based on a hierarchical coarsening approach coupled with an efficientsearch over the feasible configuration tree. Lastly, we present several numerical examples related toelectromagnetic wave propagation problems, and evaluate computational cost improvements.

In this study, the mapped Galerkin method was applied to solve the vector wave equation based on H&amp / #8722 / field and to obtain the propagation constant in x &amp / #8722 / y space. The vector wave equation was solved by the transformation of the infinite x &amp / #8722 / y plane onto a unit square. Two-dimensional Fourier series expansions were used in the solutions. Modal fields and propagation constants of dielectric waveguides were calculated. In the first part of the study, all of the calculations were made in step index fibers. Transverse magnetic fields were obtained in the u &amp / #8722 / v and x &amp / #8722 / y space through the solution of the matrix eigenvalue equation. Some graphics were plotted in the light of the results obtained. The results are found to be in accord with the results of other numerical techniques and exact solutions. After that, the propagation constant in x&amp / #8722 / y space was calculated with ease using the solution of the modal field components. In the second part of the study, the similar calculations were made in graded index fibers.

Incompressible flows -- flows in which variations in the density of a fluid are negligible -- arise in a wide variety of applications, from hydraulics to aerodynamics. The incompressible Navier-Stokes equations which govern such flows are also of fundamental physical and mathematical interest. They are believed to hold the key to understanding turbulent phenomena; precise conditions for the existence and uniqueness of solutions remain unknown -- and establishing such conditions is the subject of one of the Clay Mathematics Institute's Millennium Prize Problems. Typical solutions of incompressible flow problems involve both fine- and large-scale phenomena, so that a uniform finite element mesh of sufficient granularity will at best be wasteful of computational resources, and at worst be infeasible because of resource limitations. Thus adaptive mesh refinements are required. In industry, the adaptivity schemes used are ad hoc, requiring a domain expert to predict features of the solution. A badly chosen mesh may cause the code to take considerably longer to converge, or fail to converge altogether. Typically, the Navier-Stokes solve will be just one component in an optimization loop, which means that any failure requiring human intervention is costly. Therefore, I pursue technological foundations for a solver of the incompressible Navier-Stokes equations that provides robust adaptivity starting with a coarse mesh. By robust, I mean both that the solver always converges to a solution in predictable time, and that the adaptive scheme is independent of the problem -- no special expertise is required for adaptivity. The cornerstone of my approach is the discontinuous Petrov-Galerkin (DPG) finite element methodology developed by Leszek Demkowicz and Jay Gopalakrishnan. For a large class of problems, DPG can be shown to converge at optimal rates. DPG also provides an accurate mechanism for measuring the error, and this can be used to drive adaptive mesh refinements. Several approximations to Navier-Stokes are of interest, and I study each of these in turn, culminating in the study of the steady 2D incompressible Navier-Stokes equations. The Stokes equations can be obtained by neglecting the convective term; these are accurate for "creeping" viscous flows. The Oseen equations replace the convective term, which is nonlinear, with a linear approximation. The steady-state incompressible Navier-Stokes equations approximate the transient equations by neglecting time variations. Crucial to this work is Camellia, a toolbox I developed for solving DPG problems which uses the Trilinos numerical libraries. Camellia supports 2D meshes of triangles and quads of variable polynomial order, allows simple specification of variational forms, supports h- and p-refinements, and distributes the computation of the stiffness matrix, among other features. The central contribution of this dissertation is design and development of mathematical techniques and software, based on the DPG method, for solving the 2D incompressible Navier-Stokes equations in the laminar regime (Reynolds numbers up to about 1000). Along the way, I investigate approximations to these equations -- the Stokes equations and the Oseen equations -- followed by the steady-state Navier-Stokes equations. / text

Over the last three decades, CFD simulations have become commonplace as a tool in the engineering and design of high-speed aircraft. Experiments are often complemented by computational simulations, and CFD technologies have proved very useful in both the reduction of aircraft development cycles, and in the simulation of conditions difficult to reproduce experimentally. Great advances have been made in the field since its introduction, especially in areas of meshing, computer architecture, and solution strategies. Despite this, there still exist many computational limitations in existing CFD methods; in particular, reliable higher order and hp-adaptive methods for the Navier-Stokes equations that govern viscous compressible flow. Solutions to the equations of viscous flow can display shocks and boundary layers, which are characterized by localized regions of rapid change and high gradients. The use of adaptive meshes is crucial in such settings -- good resolution for such problems under uniform meshes is computationally prohibitive and impractical for most physical regimes of interest. However, the construction of "good" meshes is a difficult task, usually requiring a-priori knowledge of the form of the solution. An alternative to such is the construction of automatically adaptive schemes; such methods begin with a coarse mesh and refine based on the minimization of error. However, this task is difficult, as the convergence of numerical methods for problems in CFD is notoriously sensitive to mesh quality. Additionally, the use of adaptivity becomes more difficult in the context of higher order and hp methods. Many of the above issues are tied to the notion of robustness, which we define loosely for CFD applications as the degradation of the quality of numerical solutions on a coarse mesh with respect to the Reynolds number, or nondimensional viscosity. For typical physical conditions of interest for the compressible Navier-Stokes equations, the Reynolds number dictates the scale of shock and boundary layer phenomena, and can be extremely high -- on the order of 10⁷ in a unit domain. For an under-resolved mesh, the Galerkin finite element method develops large oscillations which prevent convergence and pollute the solution. The issue of robustness for finite element methods was addressed early on by Brooks and Hughes in the SUPG method, which introduced the idea of residual-based stabilization to combat such oscillations. Residual-based stabilizations can alternatively be viewed as modifying the standard finite element test space, and consequently the norm in which the finite element method converges. Demkowicz and Gopalakrishnan generalized this idea in 2009 by introducing the Discontinous Petrov-Galerkin (DPG) method with optimal test functions, where test functions are determined such that they minimize the discrete linear residual in a dual space. Under the ultra-weak variational formulation, these test functions can be computed locally to yield a symmetric, positive-definite system. The main theoretical thrust of this research is to develop a DPG method that is provably robust for singular perturbation problems in CFD, but does not suffer from discretization error in the approximation of test functions. Such a method is developed for the prototypical singular perturbation problem of convection-diffusion, where it is demonstrated that the method does not suffer from error in the approximation of test functions, and that the L² error is robustly bounded by the energy error in which DPG is optimal -- in other words, as the energy error decreases, the L² error of the solution is guaranteed to decrease as well. The method is then extended to the linearized Navier-Stokes equations, and applied to the solution of the nonlinear compressible Navier-Stokes equations. The numerical work in this dissertation has focused on the development of a 2D compressible flow code under the Camellia library, developed and maintained by Nathan Roberts at ICES. In particular, we have developed a framework allowing for rapid implementation of problems and the easy application of higher order and hp-adaptive schemes based on a natural error representation function that stems from the DPG residual. Finally, the DPG method is applied to several convection diffusion problems which mimic difficult problems in compressible flow simulations, including problems exhibiting both boundary layers and singularities in stresses. A viscous Burgers' equation is solved as an extension of DPG to nonlinear problems, and the effectiveness of DPG as a numerical method for compressible flow is assessed with the application of DPG to two benchmark problems in supersonic flow. In particular, DPG is used to solve the Carter flat plate problem and the Holden compression corner problem over a range of Mach numbers and laminar Reynolds numbers using automatically adaptive schemes, beginning with very under-resolved/coarse initial meshes. / text

This paper presents three new coupling methods for interior penalty discontinuous Galerkin finite element methods and boundary element methods. The new methods allow one to use discontinuous basis functions on the interface between the subdomains represented by the finite element and boundary element methods. This feature is particularly important when discontinuous Galerkin finite element methods are used. Error and stability analysis is presented for some of the methods. Numerical examples suggest that all three methods exhibit very similar convergence properties, consistent with available theoretical results.

Neste trabalho, utilizando o método espectral de Galerkin, provamos a existência de soluções fracas (quando a dimensão n é maior que 2) e existência e unicidade de soluções fracas (quando a dimensão é 2) para um sistema de equações diferenciais parciais que descrevem o movimento de um fluido quimicamente ativo em um domínio limitado em Rn, n 2≥2. / In this work, by using the spectral Galerkin method, we prove the existence of weak solutions (when the dimension n is great than 2) and existence and uniqueness of weak solutions (when the dimension is 2) for a system of partial differential equations that describes the motion of a chemical active fluid in a bounded domain in Rn, n≥2.

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.