Numerical dispersion resulting from using the second-order central-difference operation to approximate the differential operation is the main error source of the FDTD method. The effect of numerical dispersion can be minimized if the spatial grid size is small than£f/10. It is difficultly to analyze the modeling of electrically large structures since a huge amount of computer memory will be needed if using a very fine grid to discretize the structure. Using higher-order FDTD is the effective alternative to reduce the effect of numerical dispersion. In this paper will discuss the handling of the discontinuous PEC boundary condition in four-order FDTD and its applications to antenna pattern analysis. Using the fourth-order FDTD can enlarge the spatial grid size and reduce the requirement of computer¡¦s memory. The far field range of small size antenna operating at higher frequency is shorter enough to directly derive the far field pattern by enlarging the spatial size of fourth-order FDTD. It will compare the far field pattern derived by four-order FDTD with near-to-far field transformation and analyze their characteristic individually.

Résumé : La surveillance de l’intégrité des structures (Structural Health Monitoring - SHM) est une nouvelle technologie, et comme toute nouvelle avancée technologique, elle n’a pas encore réalisé son plein potentiel. Le SHM s’appuie sur des avancées dans plusieurs disciplines, dont l’évaluation non-desctructive, les matériaux intelligents, et les capteurs et actionneurs intégrés. Une des disciplines qui permet son déploiement est la simulation numérique. Le SHM englobe une variété de techniques basées sur la génération d’ondes vibratoires et d’ondes ultrasonores guidées. L’utilisation d’ondes guidées offre en particulier une vaste gamme d’avantages. Le défi majeur associé à la pleine utilisation de la simulation numérique dans la conception d’un système SHM basé sur l’utilisation d’ondes guidées réside dans les ressources de calcul requises pour une simulation précise. La principale raison pour ces exigences est la dispersion induite par la discrétisation numérique, tel qu’indiqué dans la littérature. La méthodes des éléments spectraux (SEM) est une variante de la p-version de la méthode des éléments finis (FEM) qui offre certains outils pour solutionner le problème des erreurs de dispersion, mais la littérature souffre toujours d’une lacune dans l’étude systématique des erreurs de dispersion numérique et de sa dépendance sur les paramètres de simulation. Le présent ouvrage tente de combler cette lacune pour les théories d’ingénierie en vibrations. Il présente d’abord le développement de la formulation des éléments spectraux pour différentes théories d’ingénierie pertinentes pour la propagation des ondes vibratoires dans différents types de structures, comme des tiges et des plaques. Puis, une nouvelle technique pour le calcul des erreurs de dispersion numériques est présentée et appliquée systématiquement dans le but d’évaluer la dispersion numérique induite en termes d’erreurs dans les vitesses de propagation. Cette technique est utilisable pour les différentes formes de propagation des ondes vibratoires dans les éléments structuraux visés dans la présente thèse afin d’évaluer quantitativement les exigences de précision en termes de paramètres de maillage. Les ondes de Lamb constituent un cas particulier de la déformation plane des ondes élastiques, en raison de la présence des doubles frontières à traction libre qui couplent les ondes longitudinales et de cisaillement et qui conduisent à une infinité de modes propagatifs qui sont dispersifs par nature. La simulation des ondes de Lamb n’a pas fait l’objet d’analyse systématique de la dispersion numérique dans la littérature autant pour la SEM que la FEM. Nous rapportons ici pour la première fois les résultats de l’analyse de dispersion numérique pour la propagation des ondes Lamb. Pour toutes les analyses de dispersion numérique présentées ici, l’analyse a été effectuée à˘ala fois dans le domaine fréquentiel et dans le domaine temporel. En se basant sur la nouvelle compréhension des effets de discrétisation numérique de la propagation des ondes guidées, nous étudions l’application de la SEM à la simulation numérique pour des applications de conception en SHM. Pour ce faire, l’excitation piézoélectrique est développée, et une nouvelle technique de condensation statique est développée et mise en œuvre pour les équations de la matrice semi-discrète, qui élimine le besoin de solution itérative, ainsi surnommée fortement couplée ou entièrement couplée. Cet élément piézoélectrique précis est ensuite utilisé pour étudier en détails les subtilités de la conception d’un système SHM en mettant l’accent sur la propagation des ondes de Lamb. Afin d’éviter la contamination des résultats par les réflexions sur les bords une nouvelle forme particulière d’élément absorbant a été développée et mise en œuvre. Les résultats de simulation dans le domaine fréquentiel jettent un éclairage nouveau sur les limites des modèles théoriques actuels pour l’excitation des ondes de Lamb par piézoélectriques. L’excitation par un élément piézoélectrique couplé est ensuite entièrement simulée dans le domaine temporel, et les résultats de simulation sont validés par deux cas de mesures expérimentales ainsi que par la simulation classique avec des éléments finis en utilisant le logiciel commercial ANSYS. // Abstract : Structural health monitoring (SHM) is a novel technology, and like any new technological advancement it has yet not realized its full potential. It builds on advancements in several disciplines including nondestructive evaluation, smart materials, and embedded sensors and actuators. One of the enabling disciplines is the numerical simulation. SHM encompasses a variety of techniques, vibration based, impedance and guided ultrasonic waves. Guided waves offers a wide repertoire of advantages. The major challenge facing the full utilization of the numerical simulation in designing a viable guided waves based SHM System is the formidable computational requirements for accurate simulation. The main reason for these requirements is the dispersion induced by numerical discretization as explained in the literature review. The spectral element (SEM) is a variant of the p-version finite element (FEM) that offers certain remedies to the numerical dispersion errors problem, yet it lacks a systematic study of the numerical dispersion errors and its dependence on the meshing parameters. The present work attempts to fill that gap for engineering theories. It starts by developing the formulation of the spectral element for different relevant engineering theories for guided waves propagation in various structural elements, like rods and plates. Then, extending the utility of a novel technique for computing the numerical dispersion errors, we systematically apply it in order to evaluate the numerically induced dispersion in terms of errors in the propagation speeds. This technique is employed for the various forms of guided waves propagation in structural elements covered in the present thesis in order to quantitatively assess the accuracy requirements in terms of the meshing parameters. The Lamb guided waves constitute a special case of the plane strain elastic waves, that is due to the presence of the double traction free boundaries, couple in the section plane and this coupling leads to an infinitude of propagating modes that are dispersive in nature. Lamb waves simulation have not been a subject of numerical dispersion analysis in the open literature neither for SEM nor FEM for that matter. We report here for the first time the numerical dispersion analysis results for Lamb waves propagation. For all the numerical dispersion analysis presented here, the analysis was done for both the frequency domain and time domain analysis. Based on the established understanding of the numerical discretization effects on the guided waves propagation, we utilize this knowledge to study the application of SEM to SHM simulations. In order to do so the piezoelectric excitation is developed, and a new static condensation technique is developed for the semidiscrete matrix equations, that eliminate the need for iterative solution, thus dubbed strongly coupled or fully coupled implementation. This accurate piezoelectric element are then used to study in details the intricacies of the design of an SHM system with specific emphasis on the Lamb waves propagation. In order to avoid the contamination of the results by the reflections from the edges a new special form of absorbing boundary was developed and implemented. The Simulation results in the frequency domain illuminated the limitations of the current theoretical models for piezoelectric excitation of Lamb waves. The piezoelectric excitation of a fully coupled element is then simulated in the time domain, and the results of simulation was verified against two cases of experimental measurements as well as conventional finite element simulation using the commercial software ANSYS.

In this thesis, we produce a rigorous and quantitative analysis of the errors introduced by finite difference schemes into strong constraint 4D-Variational (4D-Var) data assimilation. Strong constraint 4D-Var data assimilation is a method that solves a particular kind of inverse problem; given a set of observations and a numerical model for a physical system together with a priori information on the initial condition, estimate an improved initial condition for the numerical model, known as the analysis vector. This method has many forms of error affecting the accuracy of the analysis vector, and is derived under the assumption that the numerical model is perfect, when in reality this is not true. Therefore it is important to assess whether this assumption is realistic and if not, how the method should be modified to account for model error. Here we analyse how the errors introduced by finite difference schemes used as the numerical model, affect the accuracy of the analysis vector. Initially the 1D linear advection equation is considered as our physical system. All forms of error, other than those introduced by finite difference schemes, are initially removed. The error introduced by `representative schemes' is considered in terms of numerical dissipation and numerical dispersion. A spectral approach is successfully implemented to analyse the impact on the analysis vector, examining the effects on unresolvable wavenumber components and the l2-norm of the error. Subsequently, a similar also successful analysis is conducted when observation errors are re-introduced to the problem. We then explore how the results can be extended to weak constraint 4D-Var. The 2D linear advection equation is then considered as our physical system, demonstrating how the results from the 1D problem extend to 2D. The linearised shallow water equations extend the problem further, highlighting the difficulties associated with analysing a coupled system of PDEs.

In the thesis, we propose a newly-developed method called the method of Connected Local Fields (CLF) to analyze opto-electromagnetic passive devices. The method of CLF somewhat resembles a hybrid between the finite difference and pseudo-spectral methods. For opto-electromagnetic passive devices, our primary concern is their steady state behavior, or narrow-band characteristics, so we use a frequency-domain method, in which the system is governed by the Helmholtz equation. The essence of CLF is to use the intrinsic general solution of the Helmholtz equation to expand the local fields on the compact stencil. The original equation can then be transformed into the discretized form called LFE-9 (in 2-D case), and the intrinsic reconstruction formulae describing each overlapping local region can be obtained. Further, we present rigorous analysis of the numerical dispersion equation of LFE-9, by means of first-order approximation, and acquire the closed-form formula of the relative numerical dispersion error. We are thereby able to grasp the tangible influences brought both by the sampling density as well as the propagation direction of plane wave on dispersion error. In our dispersion analysis, we find that the LFE-9 formulation achieves the sixth-order accuracy: the theoretical highest order for discretizing elliptic partial differential equations on a compact nine-point stencil. Additionally, the relative dispersion error of LFE-9 is less than 1%, given that sampling density greater than 2.1 points per wavelength. At this point, the sampling density is nearing that of the Nyquist-Shannon sampling limit, and therefore computational efforts can be significantly reduced.

L'évaluation précise de la surface équivalente radar (SER) large bande de cibles complexes et de grande dimension est réalisée par des méthodes numériques rigoureuses. Parmi celles-ci, la méthode des différences finies dans le domaine temporel (FDTD) est bien adaptée pour effectuer ce calcul de SER sur une large bande de fréquence et obtenir une signature temporelle de la cible. Le schéma de Yee, schéma FDTD historique pour la simulation de propagation d'ondes électromagnétiques en régime transitoire, souffre de deux points faibles cruciaux: la dispersion numérique imposant une finesse de maillage; et l'approximation de la géométrie curviligne par un maillage cartésien avec des marches d'escalier détériorant la qualité des résultats. Les schémas FDTD d'ordre supérieur en espace ont été investigués pour la réduction de l'effet de la dispersion numérique. Dans cette thèse, le schéma Conservative Modified FDTD(2,4) a été développé dont les performances, en précision et en ressources, sont très intéressantes pour le calcul de SER. Liés au problème de l'approximation de la géométrie curviligne, le traitement des bords de plaques métalliques reste une difficulté non résolue pour les schémas FDTD(2,4) à stencil élargi. Les techniques conformes sont des approches développées pour le schéma de Yee, lesquelles ont été étudiées pour les schémas FDTD(2,4) afin de prendre en compte correctement la géométrie curviligne. Nous proposons une nouvelle approche reposant sur le modèle des fils obliques pour la modélisation des éléments surfaciques métalliques. Des applications SER de cibles montrent que celle-ci est prometteuse. / Rigorous numerical methods are used to compute an accurate wideband radar cross section (RCS) evaluation of large complex targets. Among these, finite differences in time domain method is appropriated for the wideband characteristic and also to obtain a transient responses of the target. The Yee scheme, known historically as an FDTD scheme for Maxwell equations, is hindered by two crucial weak points: numerical dispersion which imposes a high mesh resolution; and staircase approximation of curve geometry which deteriorates results quality. High-order space differential operator for FDTD schemes have been investigated to limit numerical dispersion errors. In this thesis, the Conservative Modified FDTD(2,4) scheme has been developed and its performance has shown very accurate results with reasonable workload for RCS computation. Relating to curve geometry modeling problem, metallic edges modeling is still an unsolved problem for FDTD(2,4) schemes with enlarged stencil. Conformal techniques have been developed for the Yee scheme and has been studied for FDTD(2,4) to accurately model curve geometry. We propose a new approach based on oblique thin wire model to model metallic surfaces. RCS computations of several targets have shown that this method is promising.

Le devenir des polluants dans les sols constitue un enjeu environnemental majeur. Dans ce travail, nous apportons une contribution à quelques méthodes numériques pour la simulation de l’écoulement et du transfert de polluants en milieu poreux variablement saturés. La propagation d’un contaminant dans les milieux souterrains dépend en premier lieu des caractéristiques de l’écoulement qui le transporte. Dans la première partie de ce travail, nous présentons la méthode des éléments finis mixtes hybrides pour la résolution de l’équation de Richards. Une procédure de condensation de la masse est proposée pour éviter l’apparition d’oscillations non physiques, notamment lors de la simulation de problèmes d’infiltration dans un milieu initialement sec.Dans la deuxième partie de ce travail, la méthode ELLAM (Eulerian-Lagrangian Localized Adjoint Method) est utilisée pour la modélisation du transport réactif en milieux fortement hétérogènes. En effet, les résultats obtenus pour le transport linéaire, décrit par l’équation d’advection-dispersion, avec les ELLAM sont très encourageants. La méthode ELLAM permet (i) de s’affranchir des contraintes de discrétisations spatiale ettemporelle imposées avec les méthodes eulériennes classiques, (ii) de conserver la masse et (iii) de traiter toutes les conditions aux limites. Par ailleurs, nous proposons une nouvelle formulation des ELLAM (C_ELLAM) permettant d’éviter les oscillations numériques et de limiter la diffusion numérique générées parla formulation standard.Dans la dernière partie, le code de calcul élaboré avec la formulation C_ELLAM est utilisé pour la caractérisation de la macrodispersion dans les milieux hétérogènes. Pour ce faire, il est indispensable de disposer d’outils de simulation précis et efficaces car cette étude est basée sur une méthode Monte Carlo nécessitant la réalisation d’un très grand nombre de simulations sur des grilles de calcul de l’ordre du million de mailles. Les résultats obtenus sont comparés avec une étude antérieure basée sur le Random WalkParticle Method. / The fate of contaminants in soils is a major environmental challenge. In this work, we develop efficient and reliable numerical tools for simulation of water flow and distribution prediction of pollutants in variably saturated porous media. In the first part of this document, the mixed hybrid finite element method is presented for solving Richard’s equation. A mass lumping technique is proposed to avoid unphysical oscillations when sharp infiltration fronts are simulated. In the second part of this work, the Eulerian Lagrangian Localized Adjoint Method (ELLAM) is used for modeling reactive transport in highly heterogeneous domains. Solute transport is described mathematically by the advection-dispersion and results obtained with ELLAM are very encouraging. ELLAM allows (i)overcoming spatial and time discretizations constraints imposed by classical Eulerian method, (ii)conserving mass and (iii) treating general boundary conditions naturally in the formulation. Moreover, we introduce a new ELLAM scheme (C_ELLAM) which avoid unphysical oscillations and reduce the numerical dispersion generated by the standard formulation.In the last part of this document, the C_ELLAM scheme is used to characterize the macrodispersion of a nonreactive solute in heterogeneous domains. This study is based on Monte Carlo simulations andtherefore requires highly efficient simulators. Our results are compared with previous work using Random Walk Particle Method to solve the advection-dispersion equation.

Complex structures and multiple interfaces of modern microelectronic packages complicate the interpretation of acoustic data. This study has four novel contributions. 1) Contributions to the finite element method. 2) Novel approaches to reduce computational cost. 3) New post processing technologies to interpret the simulation data. 4) Formation of theoretical guidance for acoustic image interpretation. The impact of simulation resolution on the numerical dispersion error and the exploration of quadrilateral infinite boundaries make up the first part of this thesis's contributions. The former focuses on establishing the convergence score of varying resolution densities in the time and spatial domain against a very high fidelity numerical solution. The latter evaluates the configuration of quadrilateral infinite boundaries in comparison against traditional circular infinite boundaries and quadrilateral Perfectly Matched Layers. The second part of this study features the modelling of a flip chip with a 140µm solder bump assembly, which is implemented with a 230MHz virtual raster scanning transducer with a spot size of 17µm. The Virtual Transducer was designed to reduce the total numerical elements from hundreds of millions to hundreds of thousands. Thirdly, two techniques are invented to analyze and evaluate simulated acoustic data: 1) The C-Line plot is a 2D max plot of specific gate interfaces that allows quantitative characterization of acoustic phenomena. 2) The Acoustic Propagation Map, contour maps an overall summary of intra sample wave propagation across the time domain in one image. Lastly, combining all the developments. The physical mechanics of edge effects was studied and verified against experimental data. A direct relationship between transducer spot size and edge effect severity was established. At regions with edge effect, the acoustic pulse interfacing with the solder bump edge is scattered mainly along the horizontal axis. The edge effect did not manifest in solder bump models without Under Bump Metallization (UBM). Measurements found acoustic penetration improvements of up to 44% with the removal of (UBM). Other acoustic mechanisms were also discovered and explored. Defect detection mechanism was investigated by modelling crack propagation in the solder bump assembly. Gradual progression of the crack was found have a predictable influence on the edge effect profile. By exploiting this feature, the progress of crack propagation from experimental data can be interpreted by evaluating the C-Scan image.

L'interaction dynamique sol-structure a été largement explorée en supposant le comportement linéaire du sol. Néanmoins, pour des séismes d'intensité modérée à forte, la contrainte de cisaillement maximale peut facilement atteindre la limite élastique du sol. Du point de vue de l'interaction sol-structure, les effets non linéaires peuvent modifier la rigidité du sol à la base de la structure ainsi que la quantité d'énergie dissipée dans le sol. En conséquence, ignorer les caractéristiques non linéaires du sol dans l'interaction dynamique sol-structure (IDSS) peut conduire à des prédictions erronées de la réponse de la structure. Le but de ce travail est d'implémenter dans un code numérique une loi de comportement non linéaire pour le sol afin d'examiner l'effet de la nonlinéarité du sol sur l'interaction dynamique sol-structure. De plus, différents aspects sont pris en compte tels que l'effet de la contrainte de confinement sur le module de cisaillement du sol, les conditions statiques initiales, les conditions d'interface entre le sol et la structure, etc. Durant ce travail, une méthode simple de couche absorbante basée sur une formulation de Rayleigh / Caughey pour l'amortissement, qui est généralement disponible dans les logiciels existants d'éléments finis, a également été développée. Les conditions de stabilité des problèmes de propagation d'onde sont étudiées et on montre que les comportements linéaire et non linéaire sont très différents en ce qui concerne la dispersion numérique. La règle habituelle de 10 points par longueur d'onde, recommandée dans la littérature pour les milieux élastiques, apparaît pas suffisante dans le cas non linéaire.Le modèle implémenté est d'abord vérifié numériquement en comparant les résultats avec ceux d'autres codes numériques connus. Après cela, une étude paramétrique est menée pour différents types de structures et des profils de sol variés afin de caractériser les effets non linéaires. Différentes caractéristiques de l'IDSS sont comparées à celles du cas linéaire: modification de l'amplitude et du contenu fréquentiel des ondes se propageant dans le sol, fréquence fondamentale, dissipation de l'énergie dans le sol et réponse du système sol-structure. A travers ces études paramétriques nous montrons qu'en fonction des propriétés du sol, le contenu fréquentiel de la réponse du sol peut changer significativement à cause des nonlinéarités de comportement. Les pics de la fonction de transfert entre le champ libre et le rocher affleurant se décalent vers les basses fréquences et l'amplification se produit dans cette gamme de fréquences. Une réduction de l'amplification pour les hautes fréquences et même une dé-amplification peuvent se produire pour un fort niveau des mouvements d'entrée. Ces changements influencent la réponse de la structure. Ce travail montre également que la proximité des fréquences fondamentales de la structure et du sol influence fortement l'interaction sol-structure. Enfin, l'effet du poids de la structure et du balancement de la superstructure peut être significatif. Finalement, le bassin de Nice est utilisé comme un exemple de propagation d'onde dans un milieu non linéaire hétérogène et d'interaction dynamique sol-structure. La réponse du bassin dépend fortement de la combinaison de la nonlinéarité du sol, des effets topographiques et du contraste d'impédance entre les couches de sol. Pour les structures et les profils de sol sélectionnés dans ce travail, les simulations numériques réalisées montrent que le décalage de la fréquence fondamentale n'est pas un bon indicateur pour distinguer le comportement linéaire du sol du comportement non linéaire / The interaction of the soil with the structure has been largely explored the assumption of material and geometrical linearity of the soil. Nevertheless, for moderate or strong seismic events, the maximum shear strain can easily reach the elastic limit of the soil behavior. Considering soil-structure interaction, the nonlinear effects may change the soil stiffness at the base of the structure and therefore energy dissipation into the soil. Consequently, ignoring the nonlinear characteristics of the dynamic soil-structure interaction (DSSI) this phenomenon could lead toerroneous predictions of structural response. The goal of this work is to implement a fully nonlinear constitutive model for soils into anumerical code in order to investigate the effect of soil nonlinearity on dynamic soil structureinteraction. Moreover, different issues are taken into account such as the effect of confining stress on the shear modulus of the soil, initial static condition, contact elements in the soil-structure interface, etc. During this work, a simple absorbing layer method based on a Rayleigh / Caughey damping formulation, which is often already available in existing. Finite Element softwares, is also presented. The stability conditions of the wave propagation problems are studied and it is shown that the linear and nonlinear behavior are very different when dealing with numerical dispersion. It is shown that the 10 points per wavelength rule, recommended in the literature for the elastic media is not sufficient for the nonlinear case. The implemented model is first numerically verified by comparing the results with other known numerical codes. Afterward, a parametric study is carried out for different types of structures and various soil profiles to characterize nonlinear effects. Different features of the DSSI are compared to the linear case : modification of the amplitude and frequency content of the waves propagated into the soil, fundamental frequency, energy dissipation in the soil and the response of the soil-structure system. Through these parametric studies we show that depending on the soil properties, frequency content of the soil response could change significantly due to the soil nonlinearity. The peaks of the transfer function between free field and outcropping responsesshift to lower frequencies and amplification happens at this frequency range. Amplificationreduction for the high frequencies and even deamplication may happen for high level inputmotions. These changes influence the structural response.We show that depending on the combination of the fundamental frequency of the structureand the the natural frequency of the soil, the effect of soil-structure interaction could be significant or negligible. However, the effect of structure weight and rocking of the superstructurecould change the results. Finally, the basin of Nice is used as an example of wave propagation ona heterogeneous nonlinear media and dynamic soil-structure interaction. The basin response isstrongly dependent on the combination of soil nonlinearity, topographic effects and impedancecontrast between soil layers. For the selected structures and soil profiles of this work, the performed numerical simulations show that the shift of the fundamental frequency is not a goodindex to discriminate linear from nonlinear soil behavior

The interaction of the soil with the structure has been largely explored the assumption of material and geometrical linearity of the soil. Nevertheless, for moderate or strong seismic events, the maximum shear strain can easily reach the elastic limit of the soil behavior. Considering soil-structure interaction, the nonlinear effects may change the soil stiffness at the base of the structure and therefore energy dissipation into the soil. Consequently, ignoring the nonlinear characteristics of the dynamic soil-structure interaction (DSSI) this phenomenon could lead toerroneous predictions of structural response. The goal of this work is to implement a fully nonlinear constitutive model for soils into anumerical code in order to investigate the effect of soil nonlinearity on dynamic soil structureinteraction. Moreover, different issues are taken into account such as the effect of confining stress on the shear modulus of the soil, initial static condition, contact elements in the soil-structure interface, etc. During this work, a simple absorbing layer method based on a Rayleigh / Caughey damping formulation, which is often already available in existing. Finite Element softwares, is also presented. The stability conditions of the wave propagation problems are studied and it is shown that the linear and nonlinear behavior are very different when dealing with numerical dispersion. It is shown that the 10 points per wavelength rule, recommended in the literature for the elastic media is not sufficient for the nonlinear case. The implemented model is first numerically verified by comparing the results with other known numerical codes. Afterward, a parametric study is carried out for different types of structures and various soil profiles to characterize nonlinear effects. Different features of the DSSI are compared to the linear case : modification of the amplitude and frequency content of the waves propagated into the soil, fundamental frequency, energy dissipation in the soil and the response of the soil-structure system. Through these parametric studies we show that depending on the soil properties, frequency content of the soil response could change significantly due to the soil nonlinearity. The peaks of the transfer function between free field and outcropping responsesshift to lower frequencies and amplification happens at this frequency range. Amplificationreduction for the high frequencies and even deamplication may happen for high level inputmotions. These changes influence the structural response.We show that depending on the combination of the fundamental frequency of the structureand the the natural frequency of the soil, the effect of soil-structure interaction could be significant or negligible. However, the effect of structure weight and rocking of the superstructurecould change the results. Finally, the basin of Nice is used as an example of wave propagation ona heterogeneous nonlinear media and dynamic soil-structure interaction. The basin response isstrongly dependent on the combination of soil nonlinearity, topographic effects and impedancecontrast between soil layers. For the selected structures and soil profiles of this work, the performed numerical simulations show that the shift of the fundamental frequency is not a goodindex to discriminate linear from nonlinear soil behavior

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.