La méthode d’imagerie la plus répandue dans l’industrie pétrolière est la RTM (Reverse Time Migration) qui repose sur la simulation de la propagation des ondes dans le sous-sol. Nous nous sommes concentrés sur un propagateur d'ondes élastiques 3D en milieu anisotrope de type TTI (Tilted Transverse Isotropic). Nous avons directement travaillé dans le code de recherche de Total DIVA (Depth Imaging Velocity Analysis), basé sur une discrétisation par la méthode de Galerkin Discontinue et le schéma Leap-Frog, et développé pour le calcul parallèle intensif – HPC (High Performance Computing). Nous avons ciblé plus particulièrement deux contributions possibles qui, si elles supposent des compétences très différentes, ont la même finalité : réduire les coûts de calculs requis pour la simulation. D'une part, les conditions aux limites classiques de type PML (Perfectly Matched Layers) ne sont pas stables dans des milieux TTI. Nous avons proposé de formuler une CLA (Conditions aux Limites Absorbantes) stable dans des milieux anisotropes. La méthode de construction repose sur les propriétés des courbes de lenteur, ce qui donne à notre approche un caractère original. D'autre part, le parallélisme initial, basé sur une décomposition de domaine et des communications par passage de messages à l'aide de la bibliothèque MPI, conduit à un déséquilibrage de charge qui détériore son efficacité parallèle. Nous avons corrigé cela en remplaçant le paradigme parallélisme par l'utilisation de la programmation à base de tâches sur support d'exécution. Cette thèse a été réalisée dans le cadre de l'action de recherche DIP (Depth Imaging Partnership) qui lie la compagnie pétrolière Total et Inria. / The most common method of Seismic Imaging is the RTM (Reverse Time Migration) which depends on wave propagation simulations in the subsurface. We focused on a 3D elastic wave propagator in anisotropic media, more precisely TTI (Tilted Transverse Isotropic). We directly worked in the Total code DIVA (Depth Imaging Velocity Analysis) which is based on a discretization by the Discontinuous Galerkin method and the Leap-Frog scheme, and developed for intensive parallel computing – HPC (High Performance Computing). We choose to especially target two contributions. Although they required very different skills, they share the same goal: to reduce the computational cost of the simulation. On one hand, classical boundary conditions like PML (Perfectly Matched Layers) are unstable in TTI media. We have proposed a formulation of a stable ABC (Absorbing Boundary Condition) in anisotropic media. The technique is based on slowness curve properties, giving to our approach an original side. On the other hand, the initial parallelism, which is based on a domain decomposition and communications by message passing through the MPI library, leads to load-imbalance and so poor parallel efficiency. We have fixed this issue by replacing the paradigm for parallelism by the use of task-based programming through runtime system. This PhD thesis have been done in the framework of the research action DIP (Depth Imaging Partnership) between the Total oil company and Inria.

Made available in DSpace on 2016-04-28T14:16:34Z (GMT). No. of bitstreams: 1 Douglas Guilherme Schmidt.pdf: 995968 bytes, checksum: 1a8192151fc433a54ec422f421249ac4 (MD5) Previous issue date: 2008-09-15 / Secretaria da Educação do Estado de São Paulo / This historical research analyses the initial period of construction of wave mechanics, as propounded by the Austrian physicist Erwin Schrödinger (1887- 1961), emphasizing his work from December 1925 to February 1926. During these months Schrödinger created the basis of his quantum theory and wrote his two earlier papers on this subject. They were published in the Annalen der Physik. The present dissertation analyses those two papers, and some of their precedents and immediate consequences. Special attention is given to the influence of the theory of matter waves of the French physicist Louis de Broglie (1892-1987) upon the development of Schrödinger s theory, as well as other works and relevant circumstances that contributed to the creation of wave mechanics. The dissertation also discusses the interpretation given by Schrödinger to his own theory, comparing it to the approach of other physicists of that time / O presente trabalho histórico investiga a fase inicial da construção da mecânica ondulatória formulada pelo físico austríaco Erwin Schrödinger (1887- 1961), dando especial atenção aos trabalhos por ele realizados de dezembro de 1925 até fevereiro de 1926. Nesse período, Schrödinger concebeu as bases de sua teoria quântica e redigiu os dois primeiros artigos sobre o assunto, publicados na revista Annalen der Physik. Esta dissertação analisa esses dois artigos, bem como alguns de seus precedentes e repercussões. É analisada em especial a influência da teoria de ondas de matéria do físico francês Louis de Broglie (1892- 1987) no desenvolvimento da teoria de Schrödinger, bem como outros trabalhos e circunstâncias importantes que contribuíram para a elaboração da mecânica ondulatória. Discute-se também a interpretação que o próprio Schrödinger deu à sua teoria, comparando-a com o enfoque adotado por outros físicos da época

La méthode des éléments frontières pour l’équation des ondes (BEM) est utilisée en acoustique eten électromagnétisme pour simuler la propagation d’une onde avec une discrétisation en temps(TD). Elle permet d’obtenir un résultat pour plusieurs fréquences à partir d’une seule résolution.Dans cette thèse, nous nous intéressons à l’implémentation efficace d’un simulateur TD-BEM sousdifférents angles. Nous décrivons le contexte de notre étude et la formulation utilisée qui s’exprimesous la forme d’un système linéaire composé de plusieurs matrices d’interactions/convolutions.Ce système est naturellement calculé en utilisant l’opérateur matrice/vecteur creux (SpMV). Nousavons travaillé sur la limite du SpMV en étudiant la permutation des matrices et le comportementde notre implémentation aidé par la vectorisation sur CPU et avec une approche par bloc surGPU. Nous montrons que cet opérateur n’est pas approprié pour notre problème et nous proposonsde changer l’ordre de calcul afin d’obtenir une matrice avec une structure particulière.Cette nouvelle structure est appelée une matrice tranche et se calcule à l’aide d’un opérateur spécifique.Nous décrivons des implémentations optimisées sur architectures modernes du calculhaute-performance. Le simulateur résultant est parallélisé avec une approche hybride (mémoirespartagées/distribuées) sur des noeuds hétérogènes, et se base sur une nouvelle heuristique pouréquilibrer le travail entre les processeurs. Cette approche matricielle a une complexité quadratiquesi bien que nous avons étudié son accélération par la méthode des multipoles rapides (FMM). Nousavons tout d’abord travaillé sur la parallélisation de l’algorithme de la FMM en utilisant différentsparadigmes et nous montrons comment les moteurs d’exécution sont adaptés pour relâcher le potentielde la FMM. Enfin, nous présentons des résultats préliminaires d’un simulateur TD-BEMaccéléré par FMM . / The time-domain BEM for the wave equation in acoustics and electromagnetism is used to simulatethe propagation of a wave with a discretization in time. It allows to obtain several frequencydomainresults with one solve. In this thesis, we investigate the implementation of an efficientTD-BEM solver using different approaches. We describe the context of our study and the TD-BEMformulation expressed as a sparse linear system composed of multiple interaction/convolutionmatrices. This system is naturally computed using the sparse matrix-vector product (SpMV). Wework on the limits of the SpMV kernel by looking at the matrix reordering and the behavior of ourSpMV kernels using vectorization (SIMD) on CPUs and an advanced blocking-layout on NvidiaGPUs. We show that this operator is not appropriate for our problem, and we then propose toreorder the original computation to get a special matrix structure. This new structure is called aslice matrix and is computed with a custom matrix/vector product operator. We present an optimizedimplementation of this operator on CPUs and Nvidia GPUs for which we describe advancedblocking schemes. The resulting solver is parallelized with a hybrid strategy above heterogeneousnodes and relies on a new heuristic to balance the work among the processing units. Due tothe quadratic complexity of this matrix approach, we study the use of the fast multipole method(FMM) for our time-domain BEM solver. We investigate the parallelization of the general FMMalgorithm using several paradigms in both shared and distributed memory, and we explain howmodern runtime systems are well-suited to express the FMM computation. Finally, we investigatethe implementation and the parametrization of an FMM kernel specific to our TD-BEM, and weprovide preliminary results.

Our knowledge of elastic wave propagation in general heterogeneous media with complex geological structures comes principally from numerical simulations. In this dissertation, I demonstrate through rigorous theoretical analyses and comprehensive numerical experiments that the Fourier method is a suitable method of choice for large scale 3D seismic modeling and imaging problems, due to its high accuracy and computational efficiency. The most attractive feature of the Fourier method is its ability to produce highly accurate solutions on relatively coarser grids, compared with other numerical methods for solving wave equations. To further advance the Fourier method, I identify two aspects of the method to focus on in this work, i.e., its implementation on modern clusters of computers and efficient high-order time stepping schemes. I propose two new parallel algorithms to improve the efficiency of the Fourier method on distributed memory systems using MPI. The first algorithm employs non-blocking all-to-all communications to optimize the conventional parallel Fourier modeling workflows by overlapping communication with computation. With a carefully designed communication-computation overlapping mechanism, a large amount of communication overhead can be concealed when implementing different kinds of wave equations. The second algorithm combines the advantages of both the Fourier method and the finite difference method by using convolutional high-order finite difference operators to evaluate the spatial derivatives in the decomposed direction. The high-order convolutional finite difference method guarantees a satisfactory accuracy and provides the flexibility of using non-blocking point-to-point communications for efficient interprocessor data exchange and the possibility of overlapping communication and computation. As a result, this hybrid method achieves an optimized balance between numerical accuracy and computational efficiency. To improve the overall accuracy of time domain Fourier simulations, I propose a family of new high-order time stepping schemes, based on a novel algorithm for designing time integration operators, to reduce temporal derivative discretization errors in a cost-effective fashion. I explore the pseudo-analytical method and propose high-order formulations to further improve its accuracy and ability to deal with spatial heterogeneities. I also extend the pseudo-analytical method to solve the variable-density acoustic and elastic wave equations. I thoroughly examine the finite difference method by conducting complete numerical dispersion and stability analyses. I comprehensively compare the finite difference method with the Fourier method and provide a series of detailed benchmarking tests of these two methods under a number of different simulation configurations. The Fourier method outperforms the finite difference method, in terms of both accuracy and efficiency, for both the theoretical studies and the numerical experiments, which provides solid evidence that the Fourier method is a superior scheme for large scale seismic modeling and imaging problems. / text

No description available.

Bragg scattering of optical waves in thick gratings at extreme angles, where the scattered wave propagates parallel (extremely asymmetric scattering - EAS) or nearly parallel (grazing angle scattering - GAS) to the grating boundaries, is associated with many unique and practically important resonant phenomena. It has been demonstrated that one of the main physical mechanisms for these resonant phenomena is the diffractional divergence of the scattered wave inside and outside the grating region. This thesis fills the gaps in the theoretical and experimental understanding of Bragg scattering in gratings at extreme angles by investigating EAS and GAS in structures where diffractional divergence of waves is significantly affected by anisotropy and/or non-uniformities of the dielectric permittivity. Unusually high sensitivity of wave scattering in thick periodic gratings to small step-like variations of mean structural parameters at the grating boundaries is predicted and described for the case when the scattered wave (the +1 diffracted order) propagates almost parallel to the front grating boundary (the geometry of GAS). A unusual pattern of strong multiple resonances for bulk electromagnetic waves is predicted and analysed numerically in thick periodic holographic gratings in a guiding slab with mean permittivity that is greater than that of the surrounding media. It is demonstrated that these resonances are related to resonant generation of a new type of eigenmodes in a thick slab with a periodic grating. These eigenmodes are generically related to the grating -- they do exist not if the grating amplitude is zero. A new type of resonant coupling of bulk radiation into the conventional guided modes of a slab with a thick holographic grating is predicted and explained theoretically. It occurs in the presence of strong frequency detunings of the Bragg condition by means of interaction of the strongly non-eigen +1 diffracted order with the slab-grating boundaries. Therefore, it is only in the presence of step-like variations of the mean permittivity at the grating boundaries that this type of resonant coupling can occur. A new method for the analysis of EAS and GAS in anisotropic gratings is developed. This method is based on the consideration of the diffractional divergence of the scattered wave and the two-wave approximation in anisotropic gratings. Special efforts are focused on the analysis of EAS and GAS of extraordinary waves in uniaxial gratings. In particular, it is demonstrated that increasing curvature of the normal surface in the direction of propagation of the scattered wave results in increase of its diffraction divergence and the resonant amplitude. A theoretical model is developed for comparison of the theoretical predictions with data obtained from experimental observations of EAS in a holographic grating written in a photorefractive medium. The developed model is applied for the interpretation of experimental observations of EAS in BaTiO3 photorefractive crystals. Good agreement with the theoretical predictions is demonstrated.

Esta dissertação apresenta um formalismo baseado no Método dos Elementos Finitos (MEF), adequado a análise de descontinuidades coaxiais com simetria axial, entre duas linhas coaxiais quaisquer, incluindo corrugações nos tubos internos e externos. Pelo método de Galerkin-Budnov deduz-se um operador integral bilinear, aplicado ao campo magnético, expandido em todo o domínio da estrutura, coaxial-descontinuidade-coaxial. As portas de entrada e saída da estrutura são posicionadas distantes da descontinuidade, de forma que nelas só haja o modo Transversal Eletromagnético (TEM). O campo magnético procurado e obtido pelo MEF. Os resultados encontrados; perdas de retorno, comportamento do campo magnético e as equi-fases nas portas de entrada e saída da estrutura, foram calculados e confrontados com as do Método de Casamento de Modos (MCM), com um alto grau de concordância. / This paper presents a formalism based on the Finite Element Method (FEM), suitable for analysis of discontinuities with axial symmetry, between any two coaxial lines, including corrugations. The method of Galerkin-Budnov follows a bi-linear integral operator, applied to the magnetic held, expanded in the whole area of the structure, coaxial-discontinuity-coaxial. The input and output ports of the structure are far from of the discontinuity, so that therein there is only the TEM mode. The magnetic held sought is obtained by MEF. The values for the return loss, and the magnetic helds at both ports were calculated and compared with those yielded by Modes Matching Technique, showing good agreement.

Les méthodes d'éléments finis de type Galerkine discontinu (DG FEM) ont démontré précision et efficacité pour résoudre des problèmes d'ondes dans des milieux complexes. Cependant, elles nécessitent un très grand nombre de degrés de liberté, ce qui augmente leur coût de calcul en comparaison du coût des méthodes d'éléments finis continus.Parmi les différentes approches variationnelles pour résoudre les problèmes aux limites, se distingue une famille particulière, basée sur l'utilisation de fonctions tests qui sont des solutions locales exactes des équations à résoudre. L'idée vient de E.Trefftz en 1926 et a depuis été largement développée et généralisée. Les méthodes variationnelles de type Trefftz-DG appliquées aux problèmes d'ondes se réduisent à des intégrales de surface, ce qui devrait contribuer à réduire les coûts de calcul.Les approches de type Trefftz ont été largement développées pour les problèmes harmoniques, mais leur utilisation pour des simulations en domaine transitoire est encore limitée. Quand elles sont appliquées dans le domaine temporel, les méthodes de Trefftz utilisent des maillages qui recouvrent le domaine espace-temps. C'est une des paraticularités de ces méthodes. En effet, les méthodes DG standards conduisent à la construction d'un système semi-discret d'équations différentielles ordinaires en temps qu'on intègre avec un schéma en temps explicite. Mais les méthodes de Trefftz-DG appliquées aux problèmes d'ondes conduisent à résoudre une matrice globale, contenant la discrétisation en espace et en temps, qui est de grande taille et creuse. Cette particularité gêne considérablement le déploiement de cette technologie pour résoudre des problèmes industriels.Dans ce travail, nous développons un environnement Tre#tz-DG pour résoudre des problèmes d'ondes mécaniques, y compris les équations couplées de l'élasto-acoustique. Nous prouvons que les formulations obtenues sont bien posées et nous considérons la difficulté d'inverser la matrice globale en construisant un inverse approché obtenu à partir de la décomposition de la matrice globale en une matrice diagonale par blocs. Cette idée permet de réduire les coûts de calcul mais sa précision est limitée à de petits domaines de calcul. Etant données les limitations de la méthode, nous nous sommes intéressés au potentiel du "Tent Pitcher", en suivant les travaux récents de Gopalakrishnan et al. Il s'agit de construire un maillage espace-temps composé de macro-éléments qui peuvent être traités indépendamment en faisant une hypothèse de causalité. Nous avons obtenu des résultats préliminaires très encourageants qui illustrent bien l'intérêt du Tent Pitcher, en particulier quand il est couplé à une méthode de Trefftz-DG formulée à partir d'intégrales de surface seulement. Dans ce cas, le maillage espace-temps est composé d'éléments qui sont au plus de dimension 3. Il est aussi important de noter que ce cadre se prête à l'utilisation de pas de temps locaux ce qui est un plus pour gagner en précision avec des coûts de calcul réduits. / Discontinuous Finite Element Methods (DG FEM) have proven flexibility and accuracy for solving wave problems in complex media. However, they require a large number of degrees of freedom, which increases the corresponding computational cost compared with that of continuous finite element methods. Among the different variational approaches to solve boundary value problems, there exists a particular family of methods, based on the use of trial functions in the form of exact local solutions of the governing equations. The idea was first proposed by Trefftz in 1926, and since then it has been further developed and generalized. A Trefftz-DG variational formulation applied to wave problems reduces to surface integrals that should contribute to decreasing the computational costs.Trefftz-type approaches have been widely used for time-harmonic problems, while their implementation for time-dependent simulations is still limited. The feature of Trefftz-DG methods applied to time-dependent problems is in the use of space-time meshes. Indeed, standard DG methods lead to the construction of a semi-discrete system of ordinary differential equations in time which are integrated by using an appropriate scheme. But Trefftz-DG methods applied to wave problems lead to a global matrix including time and space discretizations which is huge and sparse. This significantly hampers the deployment of this technology for solving industrial problems.In this work, we develop a Trefftz-DG framework for solving mechanical wave problems including elasto-acoustic equations. We prove that the corresponding formulations are well-posed and we address the issue of solving the global matrix by constructing an approximate inverse obtained from the decomposition of the global matrix into a block-diagonal one. The inversion is then justified under a CFL-type condition. This idea allows for reducing the computational costs but its accuracy is limited to small computational domains. According to the limitations of the method, we have investigated the potential of Tent Pitcher algorithms following the recent works of Gopalakrishnan et al. It consists in constructing a space-time mesh made of patches that can be solved independently under a causality constraint. We have obtained very promising numerical results illustrating the potential of Tent Pitcher in particular when coupled with a Trefftz-DG method involving only surface terms. In this way, the space-time mesh is composed of elements which are 3D objects at most. It is also worth noting that this framework naturally allows for local time-stepping which is a plus to increase the accuracy while decreasing the computational burden.

Esta dissertação apresenta um formalismo baseado no Método dos Elementos Finitos (MEF), adequado a análise de descontinuidades coaxiais com simetria axial, entre duas linhas coaxiais quaisquer, incluindo corrugações nos tubos internos e externos. Pelo método de Galerkin-Budnov deduz-se um operador integral bilinear, aplicado ao campo magnético, expandido em todo o domínio da estrutura, coaxial-descontinuidade-coaxial. As portas de entrada e saída da estrutura são posicionadas distantes da descontinuidade, de forma que nelas só haja o modo Transversal Eletromagnético (TEM). O campo magnético procurado e obtido pelo MEF. Os resultados encontrados; perdas de retorno, comportamento do campo magnético e as equi-fases nas portas de entrada e saída da estrutura, foram calculados e confrontados com as do Método de Casamento de Modos (MCM), com um alto grau de concordância. / This paper presents a formalism based on the Finite Element Method (FEM), suitable for analysis of discontinuities with axial symmetry, between any two coaxial lines, including corrugations. The method of Galerkin-Budnov follows a bi-linear integral operator, applied to the magnetic held, expanded in the whole area of the structure, coaxial-discontinuity-coaxial. The input and output ports of the structure are far from of the discontinuity, so that therein there is only the TEM mode. The magnetic held sought is obtained by MEF. The values for the return loss, and the magnetic helds at both ports were calculated and compared with those yielded by Modes Matching Technique, showing good agreement.

Made available in DSpace on 2016-06-02T20:28:25Z (GMT). No. of bitstreams: 1 3166.pdf: 902639 bytes, checksum: ea05b6d6e2b4c76c819c3abd8b7bd595 (MD5) Previous issue date: 2010-07-08 / Financiadora de Estudos e Projetos / Neste trabalho, estudamos o resultado de boa-colocação para a equação da onda cúbica u +uR3 = 0 em R3, devido a H. Bahouri e J.-Y. Chemin, no qual os dados de Cauchy estão no espaço de Sobolev homogêneo H3/4 (R3) H-1/4 (R3). A prova utiliza um método de interpolação não-linear, decomposição de Bony e desigualdade logarítmica de Strichartz, todas formuladas na Teoria de Littlewood-Paley.