Un schéma compact de haute précision basé sur le résidu (RBC) a été développé au laboratoire SINUMEF pour la simulation numérique d'écoulements compressibles en maillages structurés. Certaines proriétés intéressantes font de ce schéma un bon choix pour les calculs d'écoulements compressibles. L'objectif de cette thèse est donc de développer un schéma basé sur le résidu (RB) en maillages non-structurés avec une précision d'ordre élevé. A cette fin, deux approches ont été explorées. La première est basé sur la méthode des volumes finies en non-structuré et conduit à un schéma basé sur le résidu appelé FV-RB. Le seconde approche s'appuie sur une nouvelle formulation spatiale dite volumes spectraux (SV) et mène au schéma SV-RB. Le schéma FV-RB a été développé à l'ordre 2 et 3. Avec cette version du schéma, de nombreux cas tests sont calculés: écoulement d'un fluide parfait et visqueux, subsonique, transonique et hypersonique, stationnaire et instationnaire, en 2D et en 3D. Une analyse de la précision et du coût de calcul est effectuée pour le schéma FV-RB. Dans la seconde approche, un schéma SV-RB est développé à l'ordre 2 et 3 pour résoudre le problème d'advection pure et les équations d'Euler. A travers quelques cas tests, une comparaison de la précision et l'efficacité est effectuée entre le schéma RB et un solveur de Riemman classique, et entre deux formulations du schéma RB développés ici.

Thèse diffusée initialement dans le cadre d'un projet pilote des Presses de l'Université de Montréal/Centre d'édition numérique UdeM (1997-2008) avec l'autorisation de l'auteur.

Le travail présenté ici avait pour but le développement d'un schéma de type Galerkin discontinu (GD) d'ordre élevé pour la résolution des équations de la dynamique des gaz écrites dans un formalisme Lagrangien total, sur des maillages bi-dimensionnels totalement déstructurés. A cette fin, une méthode progressive a été utilisée afin d'étudier étape par étape les difficultés numériques inhérentes à la discrétisation Galerkin discontinue ainsi qu'aux équations de la dynamique des gaz Lagrangienne. Par conséquent, nous avons développé dans un premier temps des schémas de type Galerkin discontinu jusqu'à l'ordre trois pour la résolution des lois de conservation scalaires mono-dimensionnelles et bi-dimensionnelles sur des maillages déstructurés. La particularité principale de la discrétisation GD présentée est l'utilisation des bases polynomiales de Taylor. Ces dernières permettent, dans le cadre de maillages bi-dimensionnels déstructurés, une prise en compte globale et unifiée des différentes géométries. Une procédure de limitation hiérarchique, basée aux noeuds et préservant les extrema réguliers a été mise en place, ainsi qu'une forme générale des flux numériques assurant une stabilité globale L2 de la solution. Ensuite, nous avons tâché d'appliquer la discrétisation Galerkin discontinue développée aux systèmes mono-dimensionnels de lois de conservation comme celui de l'acoustique, de Saint-Venant et de la dynamique des gaz Lagrangienne. Nous avons noté au cours de cette étude que l'application directe de la limitation mise en place dans le cadre des lois de conservation scalaires, aux variables physiques des systèmes mono-dimensionnels étudiés provoquait l'apparition d'oscillations parasites. En conséquence, une procédure de limitation basée sur les variables caractéristiques a été développée. Dans le cas de la dynamique des gaz, les flux numériques ont été construits afin que le système satisfasse une inégalité entropique globale. Fort de l'expérience acquise, nous avons appliqué la discrétisation GD mise en place aux équations bi-dimensionnelles de la dynamique des gaz, écrites dans un formalisme Lagrangien total. Dans ce cadre, le domaine de référence est fixe. Cependant, il est nécessaire de suivre l'évolution temporelle de la matrice jacobienne associée à la transformation Lagrange-Euler de l'écoulement, à savoir le tenseur gradient de déformation. Dans le travail présent, la transformation résultant de l'écoulement est discrétisée de manière continue à l'aide d'une base Éléments Finis. Cela permet une approximation du tenseur gradient de déformation vérifiant l'identité essentielle de Piola. La discrétisation des lois de conservation physiques sur le volume spécifique, le moment et l'énergie totale repose sur une méthode Galerkin discontinu. Le schéma est construit de sorte à satisfaire de manière exacte la loi de conservation géométrique (GCL). Dans le cas du schéma d'ordre trois, le champ de vitesse étant quadratique, la géométrie doit pouvoir se courber. Pour ce faire, des courbes de Bézier sont utilisées pour la paramétrisation des bords des cellules du maillage. Nous illustrons la robustesse et la précision des schémas mis en place à l'aide d'un grand nombre de cas tests pertinents, ainsi que par une étude de taux de convergence.

La simulation numérique est devenue un outil central dans la modélisation de nombreux systèmes physiques tels que la dynamique des fluides, les plasmas, l'électromagnétisme, etc. L'existence de phénomènes multi-échelles rend l'intégration numérique de ces modèles très difficile du point de vue de la précision et du temps de calcul. En effet, dans les méthodes classiques d'intégration temporelle, le pas de temps est limité par la taille des plus petites mailles au travers d'une relation de type CFL. De plus, la forte disparité entre le pas de temps effectif et la condition CFL favorise les phénomènes de diffusion numérique. Dans la littérature, des nombreux algorithmes à pas de temps locaux (LTS) ont été développés. Pour la plupart des algorithmes LTS, les pas de temps locaux doivent être choisis parmi les fractions du pas de temps global. Nous présentons dans cette thèse une méthode asynchrone pour l'intégration explicite des équations différentielles multi-échelles. Cette méthode repose sur l'utilisation de critères de stabilité locaux, critères déterminés non pas globalement mais à partir de conditions CFL locales. De plus, contrairement aux schémas LTS, l'algorithme asynchrone permet la sélection de pas de temps indépendants pour chaque cellule de maillage. Cette thèse comporte plusieurs volets. Le premier concerne l'étude mathématique des propriétés du schéma asynchrone. Le deuxième a pour objectif d'étudier la montée en ordre, à la fois temporelle et spatiale, des méthodes asynchrones. De nombreux développements dans le cadre des méthodes de haute précision en temps ou en espace, telles que les méthodes de type Galerkin Discontinu, peuvent offrir un cadre naturel pour l'amélioration de la précision des méthodes asynchrones. Toutefois, les estimations garantissant l'ordre de précision de ces méthodes peuvent ne pas être directement compatibles avec l'aspect asynchrone. L'objectif de cette thèse est donc de développer un schéma numérique asynchrone d'ordre élevé mais qui permet également de limiter la quantité de calculs à effectuer. Le troisième volet de cette thèse se focalise sur l'application numérique puisqu'il concerne la mise en oeuvre de la méthode asynchrone dans la simulation des cas-tests représentatifs de problèmes multi-échelles. / Numerical simulation has become a central tool for the modeling of many physical systems (Fluid dynamics, plasmas, electromagnetism, etc). Multi-scale phenomena make the integration of these physical systems difficult in terms of accuracy and computational time. Numerical time-stepping integration techniques used for modeling such problems generally fall into two categories : explicit and implicit schemes. In the explicit schemes, all unknown variables are computed at the current time level from quantities already available. The time step used is then limited by the most restrictive CFL condition over the whole computation domain. In the implicit method the time step is no longer limited by the CFL conditions. However the scheme is generally not suitable for strongly coupled problems. To solve such problems, a number of local time-stepping (LTS) approaches have been developed. These methods are restricted by a local CFL condition rather than the traditional global CFL condition. For most of these LTS algorithms, local time steps are usually selected to be fractions of the global time step so that regular meeting points in time exist, and only little work is available on LTS methods with independent time steps. We present in this thesis an asynchronous method for the explicit integration of multi-scale partial differential equations. This method is restricted by a local CFL condition rather than the traditional global CFL condition. Moreover, contrary to other LTS methods, the asynchronous algorithm permits the selection of independent time steps in each mesh element. Our work consists of several components. The first one concerns the mathematical study of the properties of the asynchronous method. the objective of the second part is to study the improvement of the convergence rate for asynchronous methods. Many approaches in the context of high precision methods in time or in space, such as the Discontinuous Galerkin methods, may offer a natural setting to improve the precision of the asynchronous methods. However, the estimates ensuring the order of the accuracy of the method may not be directly compatible with the asynchronous aspect. Then, the objective is to develop a high order asynchronous numerical scheme which also preserves the computational time reduction. Finally, the third part is focused on the implementation of the asynchronous method and illustrate the advantages of the method on test-cases representative of multiscale problems.

La physique attoseconde est un domaine en pleine expansion, intrinsèquement lié au processus de génération d’harmoniques d’ordre élevé. Cette émission, sous forme d’un train d’impulsions attosecondes ou d’une impulsion attoseconde isolée, constitue une source de lumière dans le domaine spectral extrême-UV (XUV), ultra-brève, cohérente, parfaitement synchrone du champ générateur. Deux thématiques ont été abordées. La première consiste en la caractérisation complète de l’état de polarisation des harmoniques par Polarimétrie Moléculaire en collaboration avec l’ISMO-Orsay. Cette technique est basée sur la mesure de la distribution angulaire des photoélectrons dans le référentiel moléculaire lors de l’ionisation dissociative de la molécule de NO. Nous l’appliquons à trois configurations produisant un rayonnement harmonique polarisé elliptiquement. Nous obtenons ainsi, pour la première fois, la valeur absolue de l’ellipticité harmonique, son signe, ainsi que le taux de dépolarisation.La seconde thématique est la photoionisation attoseconde résonante : nous avons étudié la photoionisation de l’hélium au voisinage de la résonance d'autoionisation 2s2p à 60.15eV, excitée par une impulsion XUV accordable et sondée par une impulsion laser IR en utilisant la technique RABBIT, qui permet la mesure de l’amplitude et de la phase spectrales de la transition résonante à deux photons. Il est ainsi possible de reconstruire dans le domaine temporel, le paquet d'ondes électronique (POE) à 2 photons. Ces mesures ont été complétées par des simulations effectuées par nos collaborateurs à UAM-Madrid et au LCPMR-Paris, qui montrent que, dans nos conditions expérimentales, ce paquet à deux photons est une image fidèle du paquet résonant à un photon. Ceci représente la première reconstruction de la dynamique temporelle d’une résonance non perturbée par le champ laser, avec une résolution attoseconde. / Attosecond physics is an expending field, intrinsically linked to the High Harmonic Generation process. This emission, which can be either an attosecond pulse train or an isolated attosecond pulse, constitutes a light source in the extreme-UV (XUV) spectral domain, coherent, perfectly synchronous of the generating field. Two thematic have been studied. The first one consists in the complete characterization of the harmonic emission through Molecular Polarimetry, in collaboration with ISMO-Orsay. This technique is based on the measurement of the Molecular Frame PhotoElectron Angular Distribution, during the dissociative ionization of NO molecules. We applied this technique to three configurations producing an elliptically polarized light. For the first time, we obtain the absolute value of the ellipticity, its sign and the depolarization rate. The second topic is the resonant attosecond photoionization: we studied the photoionization of helium, close to the 2s2p autoionization resonance at 60.15 eV, excited by a tunable XUV pulse and probed by an IR pulse, using RABBIT technique, enabling the measurement of the spectral amplitude and phase of the two photons resonant transition. From this, we can reconstruct the two-photons electron wave packets (EWP). These measurements have been completed by simulations done by our collaborator from UAM-Madrid and LCPMR-Paris, showing that, in our experimental conditions, this two photons EWP corresponds to the image of the one-photon EWP. This measurement is the first reconstruction of the temporal dynamic of a resonance non-perturbed by a laser field, with an attosecond resolution.

Dans cette thèse, on s'intéresse à la propagation d'ondes en milieu fortement hétérogène modélisée par l'équation d'Helmholtz. Les méthodes numériques permettant de résoudre ce problème souffrent de dispersion numérique, en particulier à haute fréquence. Ce phénomène, appelé "effet de pollution", est largement analysé dans la littérature quand le milieu de propagation est homogène et l'utilisation de "méthodes d'ordre élevé" est souvent proposée pour minimiser ce problème. Dans ce travail, on s'intéresse à un milieu de propagation hétérogène, cas pour lequel on dispose de moins de connaissances. On propose d'adapter des méthodes éléments finis d'ordre élevé pour résoudre l'équation d'Helmholtz en milieu hétérogène, afin de réduire l'effet de pollution. Les méthodes d'ordre élevé étant généralement basées sur des maillages "larges", une stratégie multi-échelle originale est développée afin de prendre en compte des hétérogénéités de petite échelle. La convergence de la méthode est démontrée. En particulier, on montre que la méthode est robuste vis-a-vis de l'effet de pollution. D'autre part, on applique la méthode a plusieurs cas-tests numériques. On s'intéresse d'abord à des problèmes académiques, qui permettent de valider la théorie de convergence développée. On considère ensuite des cas-tests "industriels" appliqués à la Géophysique. Ces derniers nous permettent de conclure que la méthode multi-échelle proposée est plus performante que les éléments finis "classiques" et que des problèmes 3D réalistes peuvent être considérés. / The main objective of this work is the design of an efficient numerical strategy to solve the Helmholtz equation in highly heterogeneous media. We propose a methodology based on coarse meshes and high order polynomials together with a special quadrature scheme to take into account fine scale heterogeneities. The idea behind this choice is that high order polynomials are known to be robust with respect to the pollution effect and therefore, efficient to solve wave problems in homogeneous media. In this work, we are able to extend so-called "asymptotic error-estimate" derived for problems homogeneous media to the case of heterogeneous media. These results are of particular interest because they show that high order polynomials bring more robustness with respect to the pollution effect even if the solution is not regular, because of the fine scale heterogeneities. We propose special quadrature schemes to take int account fine scale heterogeneities. These schemes can also be seen as an approximation of the medium parameters. If we denote by h the finite-element mesh step and by e the approximation level of the medium parameters, we are able to show a convergence theorem which is explicit in terms of h, e and f, where f is the frequency. The main theoretical results are further validated through numerical experiments. 2D and 3D geophysica benchmarks have been considered. First, these experiments confirm that high-order finite-elements are more efficient to approximate the solution if they are coupled with our multiscale strategy. This is in agreement with our results about the pollution effect. Furthermore, we have carried out benchmarks in terms of computational time and memory requirements for 3D problems. We conclude that our multiscale methodology is able to greatly reduce the computational burden compared to the standard finite-element method

L'adaptation de maillages est un processus itératif qui consiste à changer localement la taille et l'orientation du maillage en fonction du comportement de la solution physique étudiée. Les méthodes d'adaptation de maillages ont prouvé qu'elles pouvaient être extrêmement efficaces en réduisant significativement la taille des maillages pour une précision donnée et en atteignant rapidement une convergence asymptotique d'ordre 2 pour des problèmes contenant des singularités lorsqu'elles sont couplées à des méthodes numériques d'ordre élevé. Dans les techniques d'adaptation de maillages basées sur les métriques, deux approches ont été proposées: les méthodes multi-échelles basées sur un contrôle de l'erreur d'interpolation en norme Lp et les méthodes ciblées à une fonctionnelle qui contrôle l'erreur d'approximation sur une fonctionnelle d'intérêt via l'utilisation de l'état adjoint. Cependant, avec l'émergence de méthodes numériques d'ordre très élevé telles que la méthode de Galerkin discontinue, il devient nécessaire de prendre en compte l'ordre du schéma numérique dans le processus d'adaptation de maillages. Il est à noter que l'adaptation de maillages devient encore plus cruciale pour de tels schémas car ils ne convergent qu'à l'ordre 1 dans les singularités de l'écoulement. Par conséquent, le raffinement du maillage au niveau des singularités de la solution doit être d'autant plus important que l'ordre de la méthode est élevé. L'objectif de cette thèse sera d'étendre les résultats numériques et théoriques obtenus dans le cas de l'adaptation pour des solutions linéaires par morceaux à l'adaptation pour des solutions d'ordre élevé polynomiales par morceaux. Ces solutions sont représentées sur le maillage par des éléments finis de Lagrange d'ordre k ≥ 2. Cette thèse portera sur la modélisation de l'erreur d'interpolation locale, polynôme homogène de degré k ≥ 3 dans le formalisme du maillage continu. Or, les méthodes d'adaptation de maillages basées sur les métriques nécessitent que le modèle d'erreur soit une forme quadratique, laquelle fait apparaître intrinsèquement un espace métrique. Pour pouvoir exhiber un tel espace, il est nécessaire de décomposer le polynôme homogène et de l'approcher par une forme quadratique à la puissance k/2. Cette modélisation permet ainsi de révéler un champ de métriques indispensable pour communiquer avec le générateur de maillages. En deux et trois dimensions, des méthodes de décomposition de tenseurs telles que la décomposition de Sylvester nous permettront de décomposer la fonction exacte d'erreur puis d'en déduire le modèle d'erreur quadratique. Ce modèle d'erreur local est ensuite utilisé pour contrôler globalement l'erreur en norme Lp et le maillage optimal est obtenu en minimisant cette erreur. Dans cette thèse, on s'attachera à démontrer la convergence à l'ordre k de la méthode d'adaptation de maillages pour des fonctions analytiques et pour des simulations numériques utilisant des solveurs d'ordre k ≥ 3.

La prévention des risques industriels nécessite de simuler la dispersion turbulente de polluants. Cependant, les outils majoritairement utilisés à ce jour ne permettent pas de traiter les champs proches dans le cas de géométries complexes, et il est nécessaire d'utiliser les outils de CFD (“ Computational Fluid Dynamics ”) plus adaptés, mais plus coûteux. Afin de simuler les écoulements atmosphériques avec dispersion de polluants, les modèles CFD doivent modéliser correctement d'une part, les effets de flottabilité, et d'autre part les effets de la turbulence. Plusieurs approches existent, notamment dans la prise en compte des effets de flottabilité et la modélisation de la turbulence, et nécessitent des méthodes numériques adaptées aux spécificités mathématiques de chacune d'entre elles, ainsi que des schémas numériques précis pour ne pas polluer la modélisation. Une formulation d'ordre élevé en volumes finis, sur maillages non structurés, parallélisée, est proposée pour simuler les écoulements atmosphériques avec dispersion de polluants. L'utilisation de schémas d'ordre élevé doit permettre d'une part de réduire le nombre de cellules et diminuer les temps de simulation pour atteindre une précision donnée, et d'autre part de mieux contrôler la viscosité numérique des schémas en vue de simulations LES (Large Eddy Simulation), pour lesquelles la viscosité numérique des schémas peut masquer les effets de la modélisation. Deux schémas d'ordre élevé ont été étudiés et implémentés dans un solveur 3D Navier Stokes incompressible sur des maillages volumes finis non structurés. Nous avons développé un premier schéma d'ordre élevé, correspondant à un schéma Padé volumes finis, et nous avons étendu le schéma de reconstruction polynomiale de Carpentier (2000) aux écoulements incompressibles. Les propriétés numériques des différents schémas implémentés dans le même code de calcul sont étudiées sur différents cas tests bi-dimensionnels (calcul de flux convectifs et diffusifs sur une solution a-priori, convection d'une tâche gaussienne, décroissance d'un vortex de Taylor et cavité entraînée) et tri-dimensionnel (écoulement autour d'un obstacle cubique). Une attention particulière a été portée à l'étude de la précision et du traitement des conditions limites. L'implémentation proposée du schéma polynomial permet d'approcher, pour un maillage identique, les temps de simulation obtenus avec un schéma décentré classique d'ordre 2, mais avec une précision supérieure. Le schéma compact donne la meilleure précision. En utilisant une méthode de Jacobi sans calcul implicite de la matrice pour calculer le gradient, le temps de simulation devient intéressant uniquement lorsque la précision requise est importante. Une alternative est la résolution du système linéaire par une méthode multigrille algébrique. Cette méthode diminue considérablement le temps de calcul du gradient et le schéma Padé devient performant même pour des maillages grossiers. Enfin, pour réduire les temps de simulation, la parallélisation des schémas d'ordre élevé est réalisée par une décomposition en sous domaines. L'assemblage des flux s'effectue naturellement et différents solveurs proposés par les librairies PETSC et HYPRE (solveur multigrille algébrique et méthode de Krylov préconditionnée) permettent de résoudre les systèmes linéaires issus de notre problème. / The prevention of industrial risks requires simulating turbulent dispersion of pollutants. However, the tools mostly used so far do not allow near fields treated in the case of complex geometries, and it is necessary to utilize the tools of CFD (Computational Fluid Dynamics ") more suitable but more expensive. To simulate atmospheric flows with dispersion of pollutants, the CFD models must correctly model the one hand, the effects of buoyancy, and secondly the effects of turbulence. Several approaches exist, including taking into account the effects of buoyancy and turbulence modeling, and require numerical methods adapted to the specific mathematics of each, and accurate numerical schemes to avoid pollution modeling. A formulation of high order finite volume on unstructured meshes, parallelized, is proposed to simulate the atmospheric flows with dispersion of pollutants. The use of high order schemes allow one hand to reduce the number of cells and decrease the simulation time to achieve a given accuracy, and secondly to better control the viscosity numerical schemes for simulation LES (Large Eddy Simulation), for which the numerical viscosity patterns may mask the effects of modeling. Two high-order schemes have been studied and implemented in a 3D Navier Stokes solver on unstructured mesh finite volume. We developed the first high-order scheme, corresponding to a Padé finite volume scheme, and we have extended the scheme of reconstruction polynomial Carpentier (2000) for incompressible flows. The numerical properties of the various schemes implemented in the same computer code are studied different two-dimensional test cases (calculation of diffusive and convective flow on a solution a priori, a task Gaussian convection, decay of a vortex of Taylor and driven cavity) and tri-dimensional (flow past an obstacle cubic). Particular attention has been paid to the study of the accuracy and treatment of boundary conditions. The implementation of the polynomial allows to obtain quasi identical simulation time compared to a classical upwind scheme of order 2, but with higher accuracy. The compact layout gives the best accuracy. Using a Jacobi method without calculation implied matrix to calculate the gradient, the simulation time becomes interesting only when the required accuracy is important. An alternative is the resolution of linear system by an algebraic multigrid method. This method significantly reduces the computation time of the gradient and the Padé scheme is effective even for coarse meshes. Finally, to reduce simulation time, the parallelization schemes of high order is achieved by a decomposition into subdomains. The assembly flow occurs naturally and different solvers provided by PETSc libraries and HYORE (algebraic multigrid solver and preconditioned Krylov method) used to solve linear systems from our problem. The work was to identify and determine the parameters that lead to lowest time resolution simulation. Various tests of speed-up and scale-up were used to determine the most effective and optimal parameters for solving linear systems in parallel from our problem. The results of this work have been the subject of a communication in an international conference "Parallel CFD 2008" and an article submitted to "International Journal for Numerical Methods in Fluids" (Analysis of high-order finite volume schemes for the incompressible Navier Stokes equations)

Résumé : La synthèse de champs sonores est un domaine de recherche actif trouvant de nombreuses applications musicales, multimédias ou encore industrielles. Dans ce dernier cas, la re- construction précise du champ sonore est souhaitée, ce qui implique de répondre à un certains nombre de questionnements scientifiques. À l’aide de réseaux de microphones et de haut-parleurs, la captation, la synthèse et la reconstruction précise de champs sonores sont théoriquement possibles. Seulement, pour des applications pratiques, la disposition des haut-parleurs et l’influence acoustique du lieu de restitution sont des facteurs cruciaux à prendre en compte pour s’assurer de la bonne reconstruction du champ sonore. Dans ce contexte, cette thèse de doctorat propose des méthodes et des techniques pour la captation, la transformation et la reconstruction précise de champs sonores en trois dimen- sions en se basant sur la méthode ambisonique d’ordre élevé. Une configuration sphérique pour le réseau de microphones et de haut-parleurs est proposée. Elle suit un maillage de Lebedev à cinquante points qui permet la captation et la reconstruction du champ sonore jusqu’à l’ordre 5 avec le formalisme ambisonique. Les limitations de cette approche, tel le repliement spatial, sont étudiés en détails. De plus, une opération de transformation du champ sonore est présentée. Elle est établie dans le domaine des harmoniques sphériques et permet d’effectuer un filtrage directionnel avant le décodage pour privilégier certaines directions dans le champ sonore, suivant une fonction de directivité choisie. Pour la re- construction, une approche originale, également établie dans le domaine des harmoniques sphériques, permet de prendre en compte l’influence acoustique du lieu de restitution, ainsi que les défauts du système de restitution. Ce traitement permet alors d’adapter la synthèse de champs sonores au lieu de restitution, en conservant le formalisme théorique établi en champ libre. Finalement, une validation expérimentale des méthodes et des tech- niques développées au cours de la thèse est faite. Dans ce contexte, une suite logicielle de synthèse et traitement en temps-réel des champs sonore est développée. / Abstract : Sound field synthesis is an active research domain with various musical, multimedia or industrial applications. In the latter case, the accurate reconstruction of the sound field is targeted, which involves answering several scientific questions. Using arrays of microphones and loudspeakers, the capture, synthesis and accurate reconstruction of sound fields are theoretically possible. However, for practical applications, the arrangement of the loud- speakers and the acoustic influence of the restitution room are critical factors to consider in order to ensure the accurate reconstruction of the sound field. In this context, this thesis proposes methods and techniques for the capture, transforma- tions and accurate reconstruction of sound fields in three dimensions based on the Higher Order Ambisonics (HOA) method. A spherical configuration for the array of microphones and loudspeakers is proposed. It follows a fifty-node Lebedev grid that enables the capture and reconstruction of the sound field up to order 5 with HOA formalism. The limitations of this approach, such as the spatial aliasing, are studied in detail. A transformation op- eration of the sound field is also proposed. The formulation is established in the spherical harmonics domain and enables a directional filtering on the sound field prior to the decod- ing step. For the reconstruction of the sound field, an original approach, also established in the spherical harmonics domain, can take into account the acoustic influence of the restitution room and the defects of the playback system. This treatment then adapts the synthesis of sound fields to the restitution room, maintaining the theoretical formalism established in free field. Finally, an experimental validation of methods and techniques developed in the thesis is made. In this context, a digital signal processing toolkit is de- veloped. It process in real-time the microphones, ambisonics, and loudspeaker signals for the sound field capture, transformations, and decoding.

Cette thèse traite de l'étude numérique des propriétés et des applications des impulsions spatio-temporellement couplées, paquets d'ondes coniques et filaments laser, dans les processus fortement non-linéaires, comme la génération d'harmoniques d'ordre élevé. Nous étudions la redistribution de l'énergie au sein de ces paquets d'ondes en propagation linéaire et non-linéaire. Le flux d'énergie constitue un diagnostic des couplages spatio-temporels que nous avons appliqué à des résultats expérimentaux réels. Nous analysons l'évolution spectrale des filaments dans un gaz et nous obtenons les conditions pour la génération d'impulsions de quelques cycles dans le spectre UV. Nous étudions la génération d'harmoniques d'ordre élevé par des ondes coniques ultra-courtes. En particulier, nous montrons comment leurs propriétés de propagation influencent le champ généré dans la région X-UV. Nous étudions aussi l'interférence des différents chemins quantiques correspondant aux trajectoires électroniques. Enfin, nous obtenons la forme des faisceaux d'Airy stationnaires dans le régime non-linéaire. Pour chaque sujet, nous présentons des résultats expérimentaux qui ont motivé nos travaux ou ont été motivés par nos simulations.