Nous présentons dans cette thèse les évolutions apportées au solveur Galerkin Discontinu Teta-CLAC, issu de la collaboration IRMA-AxesSim, au cours du projet HOROCH (2015-2018). Ce solveur permet de résoudre les équations de Maxwell en 3D, en parallèle sur un grand nombre d'accélérateurs OpenCL. L'objectif du projet HOROCH était d'effectuer des simulations de grande envergure sur un modèle numérique complet de corps humain. Ce modèle comporte 24 millions de mailles hexaédriques pour des calculs dans la bande de fréquences des objets connectés allant de 1 à 3 GHz (Bluetooth). Les applications sont nombreuses : téléphonie et accessoires, sport (maillots connectés), médecine (sondes : gélules, patchs), etc. Les évolutions ainsi apportées comprennent, entre autres : l'optimisation des kernels OpenCL à destination des CPU dans le but d'utiliser au mieux les architectures hybrides ; l'expérimentation du runtime StarPU ; le design d'un schéma d'intégration à pas de temps local ; et bon nombre d'optimisations permettant au solveur de traiter des simulations de plusieurs millions de mailles. / In this thesis, we present the evolutions made to the Discontinuous Galerkin solver Teta-CLAC – resulting from the IRMA-AxesSim collaboration – during the HOROCH project (2015-2018). This solver allows to solve the Maxwell equations in 3D and in parallel on a large amount of OpenCL accelerators. The goal of the HOROCH project was to perform large-scale simulations on a complete digital human body model. This model is composed of 24 million hexahedral cells in order to perform calculations in the frequency band of connected objects going from 1 to 3 GHz (Bluetooth). The applications are numerous: telephony and accessories, sport (connected shirts), medicine (probes: capsules, patches), etc. The changes thus made include, among others: optimization of OpenCL kernels for CPUs in order to make the best use of hybrid architectures; StarPU runtime experimentation; the design of an integration scheme using local time steps; and many optimizations allowing the solver to process simulations of several millions of cells.

Elektrische und magnetische Emitter können zur Erforschung unterschiedlicher plasmonischer Nanostrukturen genutzt werden. Indem wir die Änderung der Abstrahldynamik und in der Lebensdauer bestimmen, detektieren wir die photonische lokale Zustandsdichte. Diese Zustandsdichte, die eine Eigenschaft der Umgebung ist, ermöglicht uns nicht nur Rückschlüsse auf die elektronischen und andere physikalische Eigenschaften dieser zu treffen sondern auch die allgemeinen Eigenschaften der plasmonischen Nanostruktur im Bezug auf Licht-Materie Kopplung zu bestimmen. Eine starke Licht-Materie-Kopplung ist für die zukünftige Anwendung im Bereich der Quantentechnologien wichtig. Wenn Emitter hierbei mit plasmonischen Nanostrukturen koppeln, fokussieren letztere nicht nur das emittierte Lichts an der Oberfläche im Subwellenlängenbereich sondern ermöglichen durch die Feldüberhöhung an der Oberfläche auch eine starke Licht-Materie-Kopplung. In der Arbeit konzentrieren wir uns auf zwei grundlegend unterschiedliche plasmonische Systeme: zunächst untersuchen wir analytisch den Einfluss von Graphen auf elektrische und magnetische Emitter und diskutieren dann die Lebensdaueränderungen und Strahlungsdynamiken in der Nähe von Silber- und Goldnanostrukturen. Im ersten Teil der Arbeit analysieren wir den Einfluss von Graphen mit einer Bandlücke auf den Emitter und zeigen Möglichkeiten zur experimentellen Bestimmung der Bandlücke auf. Im zweiten Teil modellieren wir die Propagation elektromagnetischer Felder im dreidimensionalen Raum mit Hilfe der Diskontinuierlichen Galerkin Zeitraum Methode mit erweiterten Funktionalitäten. Diese verwenden wir sowohl zur theoretischen Modellierung des ersten dreidimensionalen Fluoreszenlebensdauerabbildungsmikroskopie mit einem einzelnen Quantenemitter als auch zur selbstkonsistent Beschreibung von Emittern in der Nähe eines Goldpentamers. Die Kombination der Studien betont die Stärke von Emittern elektrische, optische und magnetische Eigenschaften zu detektieren. / Electric and magnetic emitters can be used to probe different plasmonic nanostructures. By determining the modification of the radiation dynamics and the lifetimes, we can measure the photonic local density of states. This, being a property of the enviroment, does not only allow us to draw conclusions regarding the electronic and other physical properties of the latter but also regarding the general light-matter coupling properties of the plasmonic nanostructure. A strong light-matter coupling is important for future applications in quantum technology. If emitters couple specifically to plasmonic nanostructure, the latter do not only focus the emitted light at the sub-wavelength scale at the surface of the structure but also allow for such a strong light-matter coupling due to the field enhancement at the surface. In this work, we focus on two different basic plasmonic systems: first, we study analytically the influence of graphene on electric and magnetic emitters, and second we discuss lifetime modifications and radiation dynamics close to silver and gold nanostructures. In the first part of this work, we specifically focus on the influence of graphene exhibiting a finite band gap on the emitter. In the second part, we model the propagation of electromagnetic fields in three-dimensional space making use of the discontinuous Galerkin time-domain method with extended functionalities. This framework we apply to model the first three-dimensional scanning-probe fluorescence-lifetime imaging microscopy by use of a single quantum-emitter as well as for a self-consistent description of emitters in the proximity of a gold pentamer. The combination of these studies stress that the strength of emitters lies in the detection of electronic, optical and magnetic properties.

In dieser Arbeit werden die nichtlokalen sowie nichtlinearen Eigenschaften plasmonischer Nanopartikel behandelt, wie sie im hydrodynamischen Modell enthalten sind. Das hydrodynamische Materialmodell stellt eine Erweiterung des Drude Modells dar, in der Korrekturen in der Beschreibung des Elektronenplasmas berücksichtigt werden. Einer ausführlichen Einführung des Materialmodells folgt eine analytische Diskussion der Auswirkungen der Nichtlokalität am Beispiel eines einzelnen Zylinders. Hierbei werden die durch die Nichtlokalität herbeigeführten Frequenzverschiebungen in den Streu- und Absorptionsspektren quantifiziert und asymptotisch behandelt. Des Weiteren wird mit Hilfe einer konformen Abbildung das Problem eines zylindrischen Dimers in der Elektrostatischen Näherung gelöst und die Moden der Struktur bestimmt. Diese Untersuchungen dienen als maßgebliche Grundlage für weiterführende numerische Studien die mit der diskontinuierlichen Galerkin Zeitraummethode durchgeführt werden. Die durch die analytischen Betrachtungen gewonnene Kenntnis der Moden ermöglicht es, im Zusammenhang mit gruppentheoretischen Betrachtungen und numerischen Untersuchungen, rigorose Auswahlregeln für die Anregung der Moden durch lineare und nichtlineare Prozesse aufzustellen. In weiterführenden numerischen Simulationen werden außerdem Strukturen niedrigerer Symmetrie, auf die sich die Auswahlregeln übertragen lassen, untersucht. Zudem werden numerische Studien präsentiert in denen der Einfluss der Nichtlokalität auf Feldüberhöhungen in Dimeren und doppel-resonantes Verhalten (es liegt sowohl bei der Frequenz des eingestrahlten Lichtes als auch bei der zweiten harmonischen eine Resonanz vor) untersucht werden. / This thesis deals with the nonlocal and nonlinear properties of plasmonic nanoparticles, as described by the hydrodynamic model. The hydrodynamic material model represents an extension of the Drude model that contains corrections to the descriptions of the electron plasma. After a thorough derivation of the material model, analytical discussions of nonlocality are presented for the example of a single cylinder. The frequency shifts in the scattering and absorption spectra are quantified and treated asymptotically. Furthermore, by applying a conformal map, the problem of a cylindrical dimer is solved in the electrostatic limit and the modes of the structure are determined. These investigations lay the foundations for numerical investigations which are performed employing the discontinuous Galerkin time domain method. The analytical knowledge of the modes, in conjunction with group theoretical considerations and numerical analysis, enables the formulation of rigorous selection rules for the excitation of modes by linear and nonlinear processes. In further numerical studies, the influence of nonlocality on the field enhancement in dimer structures and double-resonant behavior (a resonance is found at the frequency of the incoming light and at the second harmonic) are investigated.

L’imagerie du sous-sol par équations d’onde est une application de l’ingénierie pétrolière qui mobilise des ressources de calcul très importantes. On dispose aujourd’hui de calculateurs puissants qui rendent accessible l’imagerie de régions complexes mais des progrès sont encore nécessaires pour réduire les coûts de calcul et améliorer la qualité des simulations. Les méthodes utilisées aujourd’hui ne permettent toujours pas d’imager correctement des régions très hétérogènes 3D parce qu’elles sont trop coûteuses et /ou pas assez précises. Les méthodes d’éléments finis sont reconnues pour leur efficacité à produire des simulations de qualité dans des milieux hétérogènes. Dans cette thèse, on a fait le choix d’utiliser une méthode de Galerkine discontinue (DG) d’ordre élevé à flux centrés pour résoudre l’équation des ondes acoustiques et on développe un schéma d’ordre élevé pour l’intégration en temps qui peut se coupler avec la technique de discrétisation en espace, sans générer des coûts de calcul plus élevés qu’avec le schéma d’ordre deux Leap-Frog qui est le plus couramment employé. Le nouveau schéma est comparé au schéma d’ordre élevé ADER qui s’avère plus coûteux car il requiert un plus grand nombre d’opérations pour un niveau de précision fixé. De plus, le schéma ADER utilise plus de mémoire, ce qui joue aussi en faveur du nouveau schéma car la production d’images du sous-sol consomme beaucoup de mémoire et justifie de développer des méthodes numériques qui utilisent la mémoire au minimum. On analyse également la précision des deux schémas intégrés dans un code industriel et appliqués à des cas test réalistes. On met en évidence des phénomènes de pollution numériques liés à la mise en oeuvre d'une source ponctuelle dans le schéma DG et on montre qu'on peut éliminer ces ondes parasites en introduisant un terme de pénalisation non dissipatif dans la formulation DG. On finit cette thèse en discutant les difficultés engendrées par l'utilisation de schémas numériques dans un contexte industriel, et en particulier l'effet des calculs en simple précision. / Oil engineering uses a wide variety of technologies including imaging wave equation which involves very large computing resources. Very powerful computers are now available that make imaging of complex areas possible, but further progress is needed both to reduce the computational cost and improve the simulation accuracy. The current methods still do not allow to image properly heterogeneous 3D regions because they are too expensive and / or not accurate enough. Finite element methods turn out to be efficient for producing good simulations in heterogeneous media. In this thesis, we thus chose to use a high order Discontinuous Galerkin (DG) method based upon centered fluxes to solve the acoustic wave equation and developed a high-order scheme for time integration which can be coupled with the space discretization technique, without generating higher computational cost than the second-order Leap Frog scheme which is the most widely used . The new scheme is compared to the high order ADER scheme which is more expensive because it requires a larger number of computations for a fixed level of accuracy. In addition, the ADER scheme uses more memory, which also works in favor of the new scheme since producing subsurface images consumes lots of memory and justifies the development of low-memory numerical methods. The accuracy of both schemes is then analyzed when they are included in an industrial code and applied to realistic problems. The comparison highlights the phenomena of numerical pollution that occur when injecting a point source in the DG scheme and shows that spurious waves can be eliminated by introducing a non-dissipative penalty term in the DG formulation. This work ends by discussing the difficulties induced by using numerical methods in an industrial framework, and in particular the effect of single precision calculations.

La caractérisation d’objets enfouis à partir de mesures d’ondes diffractées est un problème présent dans de nombreuses applications comme l’exploration géophysique, le contrôle non-destructif, l’imagerie médicale, etc. Elle peut être obtenue numériquement par la résolution d’un problème inverse. Néanmoins, c’est un problème non linéaire et mal posé, ce qui rend la tâche difficile. Une reconstruction précise nécessite un choix judicieux de plusieurs paramètres très différents, dépendant des données de la méthode numérique d’optimisation choisie.La contribution principale de cette thèse est une étude de la reconstruction complète d’obstacles élastiques immergés à partir de mesures du champ lointain diffracté. Les paramètres à reconstruire sont la frontière, les coefficients de Lamé, la densité et la position de l’obstacle. On établit tout d’abord des résultats d’existence et d’unicité pour un problème aux limites généralisé englobant le problème direct d’élasto-acoustique. On analyse la sensibilité du champ diffracté par rapport aux différents paramètres du solide, ce qui nous conduit à caractériser les dérivées partielles de Fréchet comme des solutions du problème direct avec des seconds membres modifiés. Les dérivées sont calculées numériquement grâce à la méthode de Galerkine discontinue avec pénalité intérieure et le code est validé par des comparaisons avec des solutions analytiques. Ensuite, deux méthodologies sont introduites pour résoudre le problème inverse. Toutes deux reposent sur une méthode itérative de type Newton généralisée et la première consiste à retrouver les paramètres de nature différente indépendamment, alors que la seconde reconstruit tous les paramètre en même temps. À cause du comportement différent des paramètres, on réalise des tests de sensibilité pour évaluer l’influence de ces paramètres sur les mesures. On conclut que les paramètres matériels ont une influence plus faible sur les mesures que les paramètres de forme et, ainsi, qu’une stratégie efficace pour retrouver des paramètres de nature distincte doit prendre en compte ces différents niveaux de sensibilité. On a effectué de nombreuses expériences à différents niveaux de bruit, avec des données partielles ou complètes pour retrouver certains paramètres, par exemple les coefficients de Lamé et les paramètres de forme, la densité, les paramètres de forme et la localisation. Cet ensemble de tests contribue à la mise en place d’une stratégie pour la reconstruction complète des conditions plus proches de la réalité. Dans la dernière partie de la thèse, on étend ces résultats à des matériaux plus complexes, en particulier élastiques anisotropes. / The characterization of hidden objects from scattered wave measurements arises in many applications such as geophysical exploration, non destructive testing, medical imaging, etc. It can be achieved numerically by solving an Inverse Problem. However, this is a nonlinear and ill-posed problem, thus a difficult task. A successful reconstruction requires careful selection of very different parameters depending on the data and the chosen optimization numerical method.The main contribution of this thesis is an investigation of the full reconstruction of immersed elastic scatterers from far-field pattern measurements. The sought-after parameters are the boundary, the Lamé coefficients, the density and the location of the obstacle. First, existence and uniqueness results of a generalized Boundary Value Problem including the direct elasto-acoustic problem are established. The sensitivity of the scattered field with respect to the different parametersdescribing the solid is analyzed and we end up with the characterization of the corresponding partial Fréchet derivatives as solutions to the direct problem with modified right-hand sides. These Fréchet derivatives are computed numerically thanks to the Interior Penalty Discontinuous Galerkin method and the code is validated thanks to comparison with analytical solutions. Then, two solution methodologies are introduced for solving the inverse problem. Both are based on an iterative regularized Newton-type methodology and the first one consists in retrieving the parameters of different nature independently, while the second one reconstructs all parameters together. Due to the different behavior of the parameters, sensitivity tests are performed to assess the impact of the parameters on the measurements. We conclude that material parameters have a weaker influence on the measurements than shape parameters, and therefore, a successful strategy to retrieve parameters of distinct nature should take into account these different levels of sensitivity. Various experiments at different noise levels and with full or limited aperture data are carried out to retrieve some of the physical properties, e.g. Lamé coefficients with shape parameters, density with shape parameters a, density, shape and location. This set of tests contributes to a final strategy for the full reconstruction and in more realistic conditions. In the final part of the thesis, we extend the results to more complex material parameters, in particular anisotropic elastic.

A new combination of techniques for the numerical computation of incompressible flow is presented. The temporal discretization bases on the discontinuous Galerkin-formulation. Both constant (DG(0)) and linear approximation (DG(1)) in time is discussed. In case of DG(1) an iterative method reduces the problem to a sequence of problems each with the dimension of the DG(0) approach. For the semi-discrete problems a Galerkin/least-squares method is applied. Furthermore a non-overlapping domain decomposition method can be used for a parallelized computation. The main advantage of this approach is the low amount of information which must be exchanged between the subdomains. Due to the slight bandwidth a workstation-cluster is a suitable platform. Otherwise this method is efficient only for a small number of subdomains. The interface condition is of the Robin/Robin-type and for the Navier-Stokes equation a formulation introducing a further pressure interface condition is used. Additionally a suggestion for the implementation of the standard k-epsilon turbulence model with special wall function is done in this context. All the features mentioned above are implemented in a code called ParallelNS. Using this code the verification of this approach was done on a large number of examples ranging from simple advection-diffusion problems to turbulent convection in a closed cavity.

This thesis is concerned with motion planning for the classical two-phase Stefan problem in level set formulation. The interface separating the fluid phases from the solid phases is represented as the zero level set of a continuous function whose evolution is described by the level set equation. Heat conduction in the two phases is modeled by the heat equation. A quadratic tracking-type cost functional that incorporates temperature tracking terms and a control cost term that expresses the desire to have the interface follow a prescribed trajectory by adjusting the heat flux through part of the boundary of the computational domain. The formal Lagrange approach is used to establish a first-order optimality system by applying shape calculus tools. For the numerical solution, the level set equation and its adjoint are discretized in space by discontinuous Galerkin methods that are combined with suitable explicit Runge-Kutta time stepping schemes, while the temperature and its adjoint are approximated in space by the extended finite element method (which accounts for the weak discontinuity of the temperature by a dynamic local modification of the underlying finite element spaces) combined with the implicit Euler method for the temporal discretization. The curvature of the interface which arises in the adjoint system is discretized by a finite element method as well. The projected gradient method, and, in the absence of control constraints, the limited memory BFGS method are used to solve the arising optimization problems. Several numerical examples highlight the potential of the proposed optimal control approach. In particular, they show that it inherits the geometric flexibility of the level set method. Thus, in addition to unidirectional solidification, closed interfaces and changes of topology can be tracked. Finally, the Moreau-Yosida regularization is applied to transform a state constraint on the position of the interface into a penalty term that is added to the cost functional. The optimality conditions for this penalized optimal control problem and its numerical solution are discussed. An example confirms the efficacy of the state constraint. / Die vorliegende Arbeit beschäftigt sich mit einem Optimalsteuerungsproblem für das klassische Stefan-Problem in zwei Phasen. Die Phasengrenze wird als Niveaulinie einer stetigen Funktion modelliert, was die Lösung der so genannten Level-Set-Gleichung erfordert. Durch Anpassen des Wärmeflusses am Rand des betrachteten Gebiets soll ein gewünschter Verlauf der Phasengrenze angesteuert werden. Zusammen mit dem Wunsch, ein vorgegebenes Temperaturprofil zu approximieren, wird dieses Ziel in einem quadratischen Zielfunktional formuliert. Die notwendigen Optimalitätsbedingungen erster Ordnung werden formal mit Hilfe der entsprechenden Lagrange-Funktion und unter Benutzung von Techniken aus der Formoptimierung hergeleitet. Für die numerische Lösung müssen die auftretenden partiellen Differentialgleichungen diskretisiert werden. Dies geschieht im Falle der Level-Set-Gleichung und ihrer Adjungierten auf Basis von unstetigen Galerkin-Verfahren und expliziten Runge-Kutta-Methoden. Die Wärmeleitungsgleichung und die entsprechende Gleichung im adjungierten System werden mit einer erweiterten Finite-Elemente-Methode im Ort sowie dem impliziten Euler-Verfahren in der Zeit diskretisiert. Dieser Zugang umgeht die aufwändige Adaption des Gitters, die normalerweise bei der FE-Diskretisierung von Phasenübergangsproblemen unvermeidbar ist. Auch die Krümmung der Phasengrenze wird numerisch mit Hilfe der Methode der finiten Elemente angenähert. Zur Lösung der auftretenden Optimierungsprobleme werden ein Gradienten-Projektionsverfahren und, im Fall dass keine Kontrollschranken vorliegen, die BFGS-Methode mit beschränktem Speicherbedarf eingesetzt. Numerische Beispiele beleuchten die Stärken des vorgeschlagenen Zugangs. Es stellt sich insbesondere heraus, dass sich die geometrische Flexibilität der Level-Set-Methode auf den vorgeschlagenen Zugang zur optimalen Steuerung vererbt. Zusätzlich zur gerichteten Bewegung einer flachen Phasengrenze können somit auch geschlossene Phasengrenzen sowie topologische Veränderungen angesteuert werden. Exemplarisch, und zwar an Hand einer Beschränkung an die Lage der Phasengrenze, wird auch noch die Behandlung von Zustandsbeschränkungen mittels der Moreau-Yosida-Regularisierung diskutiert. Ein numerisches Beispiel demonstriert die Wirkung der Zustandsbeschränkung.

Title: Fluid-structure interaction of compressible flow Author: RNDr. Jaroslava Hasnedlová Department: Department of Numerical Mathematics, Institute of Applied Mathematics Supervisors: Prof. RNDr. Miloslav Feistauer, DrSc., Dr. h. c., Prof. Dr. Dr. h. c. Rolf Rannacher Supervisors' e-mail addresses: feist@karlin.mff.cuni.cz, rannacher@iwr.uni-heidelberg.de Abstract: The presented work is split into two parts. The first part is devoted to the theory of the discontinuous Galerkin finite element (DGFE) method for the space-time discretization of a nonstationary convection-diffusion initial-boundary value problem with nonlinear convection and linear diffusion. The DGFE method is applied sep- arately in space and time using, in general, different space grids on different time levels and different polynomial degrees p and q in space and time discretization. The main result is the proof of error estimates in L2 (L2 )-norm and in DG-norm formed by the L2 (H1 )-seminorm and penalty terms. The second part of the thesis deals with the realization of fluid-structure interaction problem of the compressible viscous flow with the elastic structure. The time-dependence of the domain occupied by the fluid is treated by the ALE (Arbitrary Lagrangian-Eulerian) method, when the compress- ible Navier-Stokes equations are formulated in...

Beaucoup de méthodes numériques ont été développées pour résoudre l'équation de Vlasov, car obtenir des simulations numériques précises en un temps raisonnable pour cette équation est un véritable défi. Cette équation décrit en effet l'évolution de la fonction de distribution de particules (électrons/ions) qui dépend de 3 variables d'espace, 3 variables de vitesse et du temps. L'idée principale de cette thèse est de réécrire l'équation de Vlasov sous forme d'un système hyperbolique par semi-discrétisation en vitesse. Cette semi-discrétisation est effectuée par méthode d'éléments finis. Le modèle ainsi obtenu est appelé équation de Vlasov réduite. Nous proposons différentes méthodes numériques pour résoudre efficacement ce modèle: méthodes des volumes finis, méthodes semi-Lagrangiennes et méthodes Galerkin discontinus. / Many numerical methods have been developed in order to selve the Vlasov equation, because computing precise simulations in a reasonable time is a real challenge. This equation describes the time evolution of the distribution function of charged particles (electrons/ions), which depends on 3 variables in space, 3 in velocity and time. The main idea of this thesis is to rewrite the Vlasov equation in the form of a hyperbolic system using a semi-discretization of the velocity. This semi-discretization is achieved using the finite element method. The resulting model is called the reduced Vlasov equation. We propose different numerical methods to salve this new model efficiently: finite volume methods, semi-Lagrangian methods and discontinuous Galerkin methods.

Cette thèse étudie le raffinement local de maillage à la fois en espace et en temps pour l’équation de l’elastodynamique du second ordre pour le calcul haute performance. L’objectif est de mettre en place des méthodes numériques pour traiter des hétérogénéités de petite taille ayant un impact important sur la propagation des ondes. Nous utilisons une approche par éléments finis de Galerkin discontinus avec pénalisation pour leur flexibilité et facilité de parallélisation. La formulation éléments finis que nous proposons a pour particularité d’être élasto-acoustique, pour pouvoir prendre en compte des hétérogénéités acoustiques de petite taille. Par ailleurs, nous proposons un terme de pénalisation optimisé qui est mieux adapté à l’équation de l’élastodynamique, conduisant en particulier à une meilleure condition CFL. Nous avons aussi amélioré une formulation PML du second ordre pour laquelle nous avons proposé une nouvelle discrétisation temporelle qui rend la formulation plus stable. En tirant parti de la p-adaptivité et des maillages non-conformes des méthodes de Galerkin discontinues combiné à une méthode de pas de temps local, nous avons grandement réduit le coût du raffinement local. Ces méthodes ont été implémentées en C++, en utilisant des techniques de template metaprogramming, au sein d’un code parallèle à mémoire distribuée (MPI) et partagée (OpenMP). Enfin, nous montrons le potentiel de notre approche sur des cas tests de validation et sur des cas plus réalistes avec des milieux présentant des hydrofractures. / This thesis studies local mesh refinement both in time and space for the second order elastodynamic equation in a high performance computing context. The objective is to develop numerical methods to treat small heterogeneities that have global impact on wave propagation. We use an internal penalty discontinuous Galerkin finite element approach for its flexibity and parallelization capabilities. The elasto-acoustic finite element formulation we discuss is elasto-acoustic in order to handle local acoustic heterogeneities. We also propose an optimized penalty term more suited to the elastodynamic equation that results in better CFL condition. We improve a second order PML formulation with an original time discretization that results in a more stable formulation. Using the p-adaptivity and nonconforming mesh capabilities of discontinuous Galerkin methods combined with a local time stepping method, we greatly reduce the high computational cost of local refinements. These methods have been implemented in C++, using template metaprogramming, in a distributed memory (MPI) and shared memory (OpenMP) parallel code. Finally, we show the potential of our methods on validation test cases and on more realistic test cases with medium including hydrofractures.