A control metric, the weighted sum of spatial gradients (WSSG), has been developed for use in active structural acoustic control (ASAC). Previous development of WSSG [1] showed that it was an effective control metric on simply supported plates, while being simpler to measure than other control metrics, such as volume velocity. The purpose of the current work is to demonstrate that the previous research can be generalized to plates with a wider variety of boundary conditions and on less ideal plates. Two classes of plates have been considered: clamped flat plates, and ribbed plates. On clamped flat plates an analytical model has been developed for use in WSSG that assumes the mode shapes are the product of clamped-clamped beam mode shapes. The boundary condition specific weights for use in WSSG have been derived from this formulation and provide a relatively uniform measurement field, as in the case of the simply supported plate. Using this control metric, control of radiated sound power has been simulated. The results show that WSSG provides comparable control to volume velocity on the clamped plate. Results also show, through random placement of the sensors on the plate, that similar control can be achieved regardless of sensor location. This demonstrates that WSSG is an effective control metric on a variety of boundary conditions. Ribbed plates were considered because of their wide use in aircraft and ships. In this case, a finite-element model of the plate has been used to obtain the displacement field on the plate under a variety of boundary conditions. Due to the discretized model involved, a numerical, as opposed to analytical, formulation for WSSG has been developed. Simulations using this model show that ASAC can be performed effectively on ribbed plates. In particular WSSG was found to perform comparable to or better than volume velocity on all boundary conditions examined. The sensor insensitivity property was found to hold within each section (divided by the ribs) of the plate, a slightly modified form of the flat plate insensitivity property where the plates have been shown to be relatively insensitive to sensor location over the entire surface of the plate. Improved control at natural frequencies can be achieved by applying a second control force. This confirms that ASAC is a viable option for the control of radiated sound power on non-ideal physical systems similar to ribbed plates.

This work addresses the question of the efficient numerical solution of time-domain boundary integral equations with retarded potentials arising in the problems of acoustic and electromagnetic scattering. The convolutional form of the time-domain boundary operators allows to discretize them with the help of Runge-Kutta convolution quadrature. This method combines Laplace-transform and time-stepping approaches and requires the explicit form of the fundamental solution only in the Laplace domain to be known. Recent numerical and analytical studies revealed excellent properties of Runge-Kutta convolution quadrature, e.g. high convergence order, stability, low dissipation and dispersion. As a model problem, we consider the wave scattering in three dimensions. The convolution quadrature discretization of the indirect formulation for the three-dimensional wave equation leads to the lower triangular Toeplitz system of equations. Each entry of this system is a boundary integral operator with a kernel defined by convolution quadrature. In this work we develop an efficient method of almost linear complexity for the solution of this system based on the existing recursive algorithm. The latter requires the construction of many discretizations of the Helmholtz boundary single layer operator for a wide range of complex wavenumbers. This leads to two main problems: the need to construct many dense matrices and to evaluate many singular and near-singular integrals. The first problem is overcome by the use of data-sparse techniques, namely, the high-frequency fast multipole method (HF FMM) and H-matrices. The applicability of both techniques for the discretization of the Helmholtz boundary single-layer operators with complex wavenumbers is analyzed. It is shown that the presence of decay can favorably affect the length of the fast multipole expansions and thus reduce the matrix-vector multiplication times. The performance of H-matrices and the HF FMM is compared for a range of complex wavenumbers, and the strategy to choose between two techniques is suggested. The second problem, namely, the assembly of many singular and nearly-singular integrals, is solved by the use of the Huygens principle. In this work we prove that kernels of the boundary integral operators $w_n^h(d)$ ($h$ is the time step and $t_n=nh$ is the time) exhibit exponential decay outside of the neighborhood of $d=nh$ (this is the consequence of the Huygens principle). The size of the support of these kernels for fixed $h$ increases with $n$ as $n^a,a<1$, where $a$ depends on the order of the Runge-Kutta method and is (typically) smaller for Runge-Kutta methods of higher order. Numerical experiments demonstrate that theoretically predicted values of $a$ are quite close to optimal. In the work it is shown how this property can be used in the recursive algorithm to construct only a few matrices with the near-field, while for the rest of the matrices the far-field only is assembled. The resulting method allows to solve the three-dimensional wave scattering problem with asymptotically almost linear complexity. The efficiency of the approach is confirmed by extensive numerical experiments.

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.

Dans cette thèse, nous nous intéressons à la modélisation mathématique des hétérogénéités de longueurs caractéristiques beaucoup plus petites que la longueur d'ondes. La thèse consiste en deux parties. La partie théorique est dédiée à l'obtention d'un développement asymptotique raccordé: la solution est décrite à l'aide d'un développement de champ proche au voisinage de l'obstacle et par un développement de champ lointain hors de ce voisinage. Le développement de champ lointain met en jeu des solutions singulières de l'équation des ondes tandis que le champ proche lui est régi par un modèle quasi-statique. Ces deux développements sont alors raccordés dans une zone intermédiaire dite de raccord. Nous obtenons alors des estimations d'erreurs permettant de rendre rigoureux ce développement asymptotique formel. La deuxième partie est numérique. Elle décrit à la fois la méthode de Galerkine discontinue, une méthode de raffinement de maillage espace-temps et propose une discrétisation des modèles asymptotiques obtenues précédemment. Elle est illustrée par un certain nombre de tests numériques. / In this thesis, we focus our attention on the modeling of heterogeneities which are smaller than the wavelength. The document is decomposed into two parts : a theoretical one and a numerical one. In the first part, we derive a matched asymptotic expansion composed of a far-field expansion and a near-field expansion. The terms of the far-field expansion are singular solutions of the wave equation whereas the terms of the near-field expansion satisfy quasistatic problems. These expansions are matched in an intermediate region. We justify mathematically this theory by proving error estimates. In the second part, we describe the Discontinuous Galerkin method, a local time stepping method and the implementation of the matched asymptotic method. Numerical simulations illustrate these results.

Dans cette thèse, nous nous intéressons à résoudre certains problèmes inverses pour des équations différentielles aux dérivées fractionnaires. Un problème inverse est généralement mal posé. Un problème mal posé est un problème qui ne répond pas à l’un des trois critères de Hadamard pour être bien posé, c’est-à-dire, soit l’existence, l’unicité ou une dépendance continue aux données n'est plus vraie, à savoir, des petits changements dans les données de mesure entraînent des changements indéfiniment importants dans la solution. La plupart des difficultés à résoudre des problèmes mal posés sont causées par l’instabilité de la solution. D’autre part, les équations différentielles fractionnaires deviennent un outil important dans la modélisation de nombreux problèmes de la vie réelle et il y a eu donc un intérêt croissant pour l’étude des problèmes inverses avec des équations différentielles fractionnaires. Le calcul fractionnaire est une branche des mathématiques qui fait référence à l’extension du concept de dérivation classique à la dérivation d’ordre non entier. Calculer une dérivée fractionnaire à un certain moment exige tous les processus précédents avec des propriétés de mémoire. C’est l’avantage principal du calcul fractionnaire d’expliquer les processus associés aux systèmes physiques complexes qui ont une mémoire à long terme et / ou des interactions spatiales à longue distance. De plus, les équations différentielles fractionnaires peuvent nous aider à réduire les erreurs découlant de paramètres négligés dans la modélisation des phénomènes physiques. / In this thesis, we are interested in solving some inverse problems for fractional differential equations. An inverse problem is usually ill-posed. The concept of an ill-posed problem is not new. While there is no universal formal definition for inverse problems, Hadamard [1923] defined a problem as being ill-posed if it violates the criteria of a well-posed problem, that is, either existence, uniqueness or continuous dependence on data is no longer true, i.e., arbitrarily small changes in the measurement data lead to indefinitely large changes in the solution. Most difficulties in solving ill-posed problems are caused by solution instability. Inverse problems come into various types, for example, inverse initial problems where initial data are unknown and inverse source problems where the source term is unknown. These unknown terms are to be determined using extra boundary data. Fractional differential equations, on the other hand, become an important tool in modeling many real-life problems and hence there has been growing interest in studying inverse problems of time fractional differential equations. The Non-Integer Order Calculus, traditionally known as Fractional Calculus is the branch of mathematics that tries to interpolate the classical derivatives and integrals and generalizes them for any orders, not necessarily integer order. The advantages of fractional derivatives are that they have a greater degree of flexibility in the model and provide an excellent instrument for the description of the reality. This is because of the fact that the realistic modeling of a physical phenomenon does not depend only on the instant time, but also on the history of the previous time, i.e., calculating timefractional derivative at some time requires all the previous processes with memory and hereditary properties.

High-frequency wave phenomena is observed in many physical settings, most notably in acoustics, electromagnetics, and elasticity. In all of these fields, numerical simulation and modeling of the forward propagation problem is important to the design and analysis of many systems; a few examples which rely on these computations are the development of metamaterial technologies and geophysical prospecting for natural resources. There are two modes of modeling the forward problem: the frequency domain and the time domain. As the title states, this work is concerned with the former regime. The difficulties of solving the high-frequency wave propagation problem accurately lies in the large number of degrees of freedom required. Conventional wisdom in the computational electromagnetics commmunity suggests that about 10 degrees of freedom per wavelength be used in each coordinate direction to resolve each oscillation. If K is the width of the domain in wavelengths, the number of unknowns N grows at least by O(K^2) for surface discretizations and O(K^3) for volume discretizations in 3D. The memory requirements and asymptotic complexity estimates of direct algorithms such as the multifrontal method are too costly for such problems. Thus, iterative solvers must be used. In this dissertation, I will present fast algorithms which, in conjunction with GMRES, allow the solution of the forward problem in O(N) or O(N log N) time. / text

Ce mémoire porte sur l’étude des équations d’évolution sur des variétés à coubure non nulle, plus particulièrement l’équation des ondes sur les espaces symétriques riemanniens de type non compact.Des propriétés de dispersion des solutions du problème de Cauchy homogène sont démontrées. Ces propriétés sont ensuite utilisées pour établir des estimations dites estimations de Strichartz. L’examen de ces estimées permet de déduire que le problème de Cauchy non linéaire avec des non-linéarités de type puissance est globalement bien posé pour des données initiales petites et localement bien posé pour des données arbitraires.Après un chapitre introductif dédié aux définitions, propriétés algébriques et géométriques des espaces symétriques et à quelques aspects élémentaires d’analyse harmonique sphérique sur ces espaces, un article est présenté : Wave equation on Riemannian symmetric spaces. Cet article contient nos résultats principaux. Dans le dernier chapitre nous présentons en détail deux problèmes ouverts qui prolongent nos travaux. Il s’agit respectivement d’établir le lien entre le comportement asymptotique des estimées et les orbites nilpotentes, et l’étude de l’équation des ondes pour les formes différentielles sur les espaces symétriques. / In this memoir we study evolution equations on curved manifolds. In particular we are interested in the wave equation on Riemannian symmetric spaces of the noncompact type.Dispersive properties of solutions of homogeneous Cauchy problem are proved. These properties are then used to establish Strichartz-type estimates. A closer study of these estimates shows that the nonlinear Cauchy problem with power-like nonlinearities is globally well posed for small initial data and locally well posed for arbitrary initial data.The first chapter is devoted to definitions, algebraic and geometric properties of symmetric spaces and to few elementary aspects of spherical analysis on these spaces. Then our main results are represented in an article : Wave equation on Riemannian symmetric spaces. In the last chapter we present in detail two open problems for future work. One issue is to establish a link between the asymptotic behavior of the estimates and nilpotent orbits, while another issue is the study of wave equation for differential forms on symmetric spaces.

Nous étudions l’analyse asymptotique en relativité générale sous deux aspects: le peeling et le scattering (diffusion) conforme. Le peeling est construit pour les champs scalaires linéaire et non-linéaires et pour les champs de Dirac en espace-temps de Kerr (qui est non-stationnaire et à symétrie simplement axiale), généralisant les travaux de L. Mason et J-P. Nicolas (2009, 2012). La méthode des champs de vecteurs (estimations d’énergie géométriques) et la technique de compactification conforme sont développées. Elles nous permettent de formuler les définitions du peeling à tous ordres et d’obtenir les données initiales optimales qui assurent ces comportements. Une théorie de la diffusion conforme pour les équations de champs sans masse de spîn n/2 dans l’espace-temps de Minkowski est construite.En effectuant les compactifications conformes (complète et partielle), l’espace-temps est complété en ajoutant une frontière constituée de deux hypersurfaces isotropes représentant respectivement les points limites passés et futurs des géodésiques de type lumière. Le comportement asymptotique des champs s’obtient en résolvant le problème de Cauchy pour l’équation rééchelonnée et en considérant les traces des solutions sur ces bords. L’inversibilité des opérateurs de trace, qui associent le comportement asymptotique passé ou futur aux données initiales, s’obtient en résolvant le problème de Goursat sur le bord conforme. L’opérateur de diffusion conforme est alors obtenu par composition de l’opérateur de trace futur avec l’inverse de l’opérateur de trace passé. / This work explores two aspects of asymptotic analysis in general relativity: peeling and conformal scattering.On the one hand, the peeling is constructed for linear and nonlinear scalar fields as well as Dirac fields on Kerr spacetime, which is non-stationary and merely axially symmetric. This generalizes the work of L. Mason and J-P. Nicolas (2009, 2012). The vector field method (geometric energy estimates) and the conformal technique are developed. They allow us to formulate the definition of the peeling at all orders and to obtain the optimal space of initial data which guarantees these behaviours. On the other hand, a conformal scattering theory for the spin-n/2 zero rest-mass equations on Minkowski spacetime is constructed. Using the conformal compactifications (full and partial), the spacetime is completed with two null hypersurfaces representing respectively the past and future end points of null geodesics. The asymptotic behaviour of fields is then obtained by solving the Cauchy problem for the rescaled equation and considering the traces of the solutions on these hypersurfaces. The invertibility of the trace operators, that to the initial data associate the future or past asymptotic behaviours, is obtained by solving the Goursat problem on the conformal boundary. The conformal scattering operator is then obtained by composing the future trace operator with the inverse of the past trace operator.

Dans cette thèse, nous présentons des résultats d’homogénéisation périodique d’un problème spectral et de l’équation d’ondes de Bloch. Il permet de modéliser les ondes à basse et haute fréquences. La partie modèle à basse fréquence est bien connu et n’est pas donc abordée dans ce travail. A contrario ; la partie à haute fréquence du modèle, qui représente des oscillations aux échelles microscopiques et macroscopiques, est un problème laissé ouvert. En particulier, les conditions aux limites de l’équation macroscopique à hautes fréquences établies dans [36] n’étaient pas connues avant le début de la thèse. Ce dernier travail apporte trois contributions principales. Les deux premières contributions, portent sur le comportement asymptotique du problème d’homogénéisation périodique du problème spectral et de l’équation des ondes en une dimension. La troisième contribution consiste en une extension du modèle du problème spectral posé dans une bande bi dimensionnelle et bornée. Le résultat d’homogénéisation comprend des effets de couche limite qui se produisent dans les conditions aux limites de l’équation macroscopique à haute fréquence. / In this dissertation, we present the periodic homogenization of a spectral problem and the waveequation with periodic rapidly varying coefficients in a bounded domain. The asymptotic behavioris addressed based on a method of Bloch wave homogenization. It allows modeling both the lowand high frequency waves. The low frequency part is well-known and it is not a new point here.In the opposite, the high frequency part of the model, which represents oscillations occurringat the microscopic and macroscopic scales, was not well understood. Especially, the boundaryconditions of the high-frequency macroscopic equation established in [36] were not known prior to thecommencement of thesis. The latter brings three main contributions. The first two contributions, areabout the asymptotic behavior of the periodic homogenization of the spectral problem and waveequation in one-dimension. The third contribution consists in an extension of the model for thespectral problem to a thin two-dimensional bounded strip Ω = (0; _) _ (0; ") _ R2. The homogenizationresult includes boundary layer effects occurring in the boundary conditions of the high-frequencymacroscopic equation.

Cette thèse porte sur l’étude des propriétés qualitatives des solutions des équations des ondes et de la chaleur semi-linéaires. Les résultats qui y sont décrits sont les suivants. Les deux premiers concernent l’existence et la description de dynamiques explosives de concentration en temps fini de l’état stationnaire à symétrie radiale dans le régime dit énergie surcritique ; en outre, pour l’équation des ondes la stabilité de ces phénomènes est étudiée dans le cas radial, et pour l’équation de la chaleur le cas plus général d’un domaine borné avec conditions de Dirichlet au bord est considéré. Le troisième porte sur la classification des dynamiques possibles près de l’état stationnaire radial pour l’équation de la chaleur dans le régime dit énergie critique, trois scénarios ayant lieu : la stabilisation, l’instabilité par explosion auto-similaire à profil explosif constant en espace, et l’instabilité par dissipation vers la solution nulle. Enfin, le quatrième a pour objet l’existence et la stabilité de profils explosifs auto-similaires non constants en espace pour l’équation de la chaleur dans le cas énergie surcritique / This thesis is devoted to the study of qualitative properties for solutions to the semilinear heat and wave equations. The results that are described are the following. The first two concern the existence and description of blow-up dynamics in which the radially symmetric stationary state is concentrated in finite time in the so-called energy supercritical regime; in addition, for the wave equation the stability of these phenomena is studied in the radial case, and for the heat equation the more general case of a bounded domain with Dirichlet condition at the boundary is considered. The third one deals with the classification of the possible dynamics near the radial stationary state for the heat equation in the so-called energy critical regime, where three scenarii occur: stabilization, instability by blow-up with the constant in space blow-up profile, and instability by dissipation to the null solution. Eventually, in the forth result we investigate the existence and the stability of self-similar blow-up profiles that are not constant in space, for the heat equation in the energy supercritical case