<p>Boundary integral equation techniques are useful in thenumerical simulation of scattering problems for wave equations.Their advantage over methods based on partial di.erentialequations comes from the lack of phase errors in the wavepropagation and from the fact that only the boundary of thescattering object needs to be discretized. Boundary integraltechniques are often applied in frequency domain but recentlyseveral time domain integral equation methods are beingdeveloped.</p><p>We study time domain integral equation methods for thescalar wave equation with a Galerkin discretization of twodi.erent integral formulations for a Dirichlet scatterer. The.rst method uses the Kirchho. formula for the solution of thescalar wave equation. The method is prone to get unstable modesand the method is stabilized using an averaging .lter on thesolution. The second method uses the integral formulations forthe Helmholtz equation in frequency domain, and this method isstable. The Galerkin formulation for a Neumann scattererarising from Helmholtz equation is implemented, but isunstable.</p><p>In the discretizations, integrals are evaluated overtriangles, sectors, segments and circles. Integrals areevaluated analytically and in some cases numerically. Singularintegrands are made .nite, using the Du.y transform.</p><p>The Galerkin discretizations uses constant basis functionsin time and nodal linear elements in space. Numericalcomputations verify that the Dirichlet methods are stable, .rstorder accurate in time and second order accurate in space.Tests are performed with a point source illuminating a plateand a plane wave illuminating a sphere.</p><p>We investigate the On Surface Radiation Condition, which canbe used as a medium to high frequency approximation of theKirchho. formula, for both Dirichlet and Neumann scatterers.Numerical computations are done for a Dirichlet scatterer.</p>

Boundary integral equation techniques are useful in thenumerical simulation of scattering problems for wave equations.Their advantage over methods based on partial di.erentialequations comes from the lack of phase errors in the wavepropagation and from the fact that only the boundary of thescattering object needs to be discretized. Boundary integraltechniques are often applied in frequency domain but recentlyseveral time domain integral equation methods are beingdeveloped. We study time domain integral equation methods for thescalar wave equation with a Galerkin discretization of twodi.erent integral formulations for a Dirichlet scatterer. The.rst method uses the Kirchho. formula for the solution of thescalar wave equation. The method is prone to get unstable modesand the method is stabilized using an averaging .lter on thesolution. The second method uses the integral formulations forthe Helmholtz equation in frequency domain, and this method isstable. The Galerkin formulation for a Neumann scattererarising from Helmholtz equation is implemented, but isunstable. In the discretizations, integrals are evaluated overtriangles, sectors, segments and circles. Integrals areevaluated analytically and in some cases numerically. Singularintegrands are made .nite, using the Du.y transform. The Galerkin discretizations uses constant basis functionsin time and nodal linear elements in space. Numericalcomputations verify that the Dirichlet methods are stable, .rstorder accurate in time and second order accurate in space.Tests are performed with a point source illuminating a plateand a plane wave illuminating a sphere. We investigate the On Surface Radiation Condition, which canbe used as a medium to high frequency approximation of theKirchho. formula, for both Dirichlet and Neumann scatterers.Numerical computations are done for a Dirichlet scatterer. / NR 20140805

This thesis is concerned with the study of transient scattering of acoustic waves by an obstacle in an infinite domain, where the scattered wave is represented in terms of time domain boundary layer potentials. The problem of finding the unknown solution of the scattering problem is thus reduced to the problem of finding the unknown density of the time domain boundary layer operators on the obstacle’s boundary, subject to the boundary data of the known incident wave. Using a Galerkin approach, the unknown density is replaced by a piecewise polynomial approximation, the coefficients of which can be found by solving a linear system. The entries of the system matrix of this linear system involve, for the case of a two dimensional scattering problem, integrals over four dimensional space-time manifolds. An accurate computation of these integrals is crucial for the stability of this method. Using piecewise polynomials of low order, the two temporal integrals can be evaluated analytically, leading to kernel functions for the spatial integrals with complicated domains of piecewise support. These spatial kernel functions are generalised into a class of admissible kernel functions. A quadrature scheme for the approximation of the two dimensional spatial integrals with admissible kernel functions is presented and proven to converge exponentially by using the theory of countably normed spaces. A priori error estimates for the Galerkin approximation scheme are recalled, enhanced and discussed. In particular, the scattered wave’s energy is studied as an alternative error measure. The numerical schemes are presented in such a way that allows the use of non-uniform meshes in space and time, in order to be used with adaptive methods that are based on a posteriori error indicators and which modify the computational domain according to the values of these error indicators. The theoretical analysis of these schemes demands the study of generalised mapping properties of time domain boundary layer potentials and integral operators, analogously to the well known results for elliptic problems. These mapping properties are shown for both two and three space dimensions. Using the generalised mapping properties, three types of a posteriori error estimators are adopted from the literature on elliptic problems and studied within the context of the two dimensional transient problem. Some comments on the three dimensional case are also given. Advantages and disadvantages of each of these a posteriori error estimates are discussed and compared to the a priori error estimates. The thesis concludes with the presentation of two adaptive schemes for the two dimensional scattering problem and some corresponding numerical experiments.

Σκοπός της παρούσας διδακτορικής διατριβής είναι η ανάπτυξη αριθμητικής μεθόδου, η οποία επιλύει προβλήματα δισδιάστατης στατικής ελαστικότητας, καθώς και δυναμικής ελαστικότητας στο πεδίο των συχνοτήτων και στο πεδίο του χρόνου. Το κύριο χαρακτηριστικό της είναι ότι η προσέγγιση του άγνωστου πεδίου γίνεται με την τοποθέτηση σημείων και όχι με τη χρήση κάποιου πλέγματος όπως γίνεται στις μέχρι τώρα κλασικές μεθοδολογίες των πεπερασμένων ή συνοριακών στοιχείων. Μέρος της παρούσας διατριβής αποτελεί και η ανάπτυξη προγράμματος ηλεκτρονικού υπολογιστή, ο οποίος υποστηρίζει πλήρως τα όσα αναφέρονται στην παρούσα εργασία. Η παρούσα διατριβή αποτελείται από δύο ενότητες. Στην πρώτη ενότητα, η οποία περιλαμβάνει τα πρώτα τρία κεφάλαια, παρατίθεται το θεωρητικό υπόβαθρο της μεθοδολογίας. Στη δεύτερη ενότητα περιγράφονται διάφορες τεχνικές λεπτομέρειες, όπως ολοκληρώσεις και προσέγγιση πεδίου και δίνονται αρκετά παραδείγματα, τα οποία πιστοποιούν την ακρίβεια και την αξιοπιστία της. / -

京都大学 / 0048 / 新制・課程博士 / 博士(情報学) / 甲第20513号 / 情博第641号 / 新制||情||111(附属図書館) / 京都大学大学院情報学研究科複雑系科学専攻 / (主査)教授 西村 直志, 教授 磯 祐介, 准教授 吉川 仁 / 学位規則第4条第1項該当 / Doctor of Informatics / Kyoto University / DFAM

<p>This investigation is motivated by the problem of optimal design of cooling elements in modern battery systems. We consider a simple model of two-dimensional steady state heat conduction described by elliptic partial differential equations (PDEs) and involving a one dimensional cooling element represented by an open contour. The problem consists in finding an optimal shape of the cooling element which will ensure that the solution in a given region is close (in the least square sense) to some prescribed target distribution. We formulate this problem as PDE-constrained optimization and the locally optimal contour shapes are found using the conjugate gradient algorithm in which the Sobolev shape gradients are obtained using methods of the shape-differential calculus combined with adjoint analysis. The main novelty of this work is an accurate and efficient approach to the evaluation of the shape gradients based on a boundary integral formulation. A number of computational aspects of the proposed approach is discussed and optimization results obtained in several test problems are presented.</p> / Master of Science (MSc)

No description available.

Cette thèse s'inscrit dans le sujet des opérateurs intégraux de frontière définis sur le disque unitaire en trois dimensions, leurs relations avec les problèmes externes de Laplace et Helmholtz, et leurs applications au préconditionnement des systèmes linéaires obtenus en utilisant la méthode des éléments finis de frontière.On décrit d'abord les méthodes intégrales pour résoudre les problèmes de Laplace et de Helmholtz en dehors des objets à frontière régulière lipschitzienne, et en dehors des surfaces bidimensionnelles ouvertes dans un espace tridimensionnel. La formulation intégrale des problèmes de Laplace et de Helmholtz pour ces cas est décrite formellement. La mise en oeuvre d'une méthode numérique utilisant la méthode des éléments finis de frontière est également décrite. Les avantages et les défis inhérents à la méthode sont abordés : la complexité du calcul numérique (de mémoire et algorithmique) et le mal conditionnement inhérentes à des systèmes linéaires associés.Dans une deuxième partie on expose une technique optimale de préconditionnement (indépendante de la discrétisation) sur la base des opérateurs intégraux de frontière. On montre comment cette technique est facilement réalisable dans le cas de problèmes définis en dehors d'un objet régulier à frontière lipschitzienne, mais qu'elle pose des problèmes quand ils sont définis en dehors d'une surface ouverte dans un espace tridimensionnel. On montre que les opérateurs intégraux de frontière associés à la résolution des problèmes de Dirichlet et Neumann définis en dehors des surfaces ont des inverses bien définis. On montre également que ceux-ci pourraient conduire à des techniques de préconditionnement optimales, mais que ses formes explicites ne sont pas faciles à obtenir. Ensuite, on montre une méthode pour obtenir de tels opérateurs inverses de façon explicite lorsque la surface sur laquelle ils sont définis est un disque unitaire dans un espace tridimensionnel. Ces opérateurs inverses explicites seront, cependant, en forme des séries, et n'auront pas une adaptation immédiate pour leur utilisation dans des méthodes des éléments finis de frontière.Dans une troisième partie on montre comment certaines modifications aux opérateurs inverses mentionnés permettent d'obtenir des expressions variationnelles explicites et fermées, non plus sous la forme des séries, en conservant certaines caractéristiques importantes pour l'effet de préconditionnement cherché. Ces formes explicites sont en effet applicables aux méthodes des éléments finis frontière. On obtient des expressions variationnelles précises et on propose des calculs numériques pour leur utilisation avec des éléments finis frontière. Ces méthodes numériques sont testées en utilisant différentes identités obtenues dans la théorie développée, et sont ensuite utilisées pour produire des matrices préconditionnantes. Leur performance en tant que préconditionneurs de systèmes linéaires associés à des problèmes de Laplace et Helmholtz à l'extérieur du disque est testée. Enfin, on propose extension de cette méthode pour couvrir les cas de surfaces diverses. Ceci est illustré et étudié dans les cas précis des problèmes extérieurs à des surfaces en forme de carré et en forme de L dans un espace tridimensionnel. / This thesis is about some boundary integral operators defined on the unit disk in a three-dimensional spaces, their relation with the exterior Laplace and Helmholtz problems, and their application to the preconditioning of the systems arising when solving these problems using the boundary element method.We begin by describing the so-called integral method for the solution of the exterior Laplace and Helmholtz problems defined on the exterior of objects with Lipschitz-regular boundaries, or on the exterior of open two-dimensional surfaces in a three-dimensional space. We describe the integral formulation for the Laplace and Helmholtz problem in these cases, their numerical implementation using the boundary element method, and we discuss its advantages and challenges: its computational complexity (both algorithmic and memory complexity) and the inherent ill-conditioning of the associated linear systems.In the second part we show an optimal preconditioning technique (independent of the chosen discretization) based on operator preconditioning. We show that this technique is easily applicable in the case of problems defined on the exterior of objects with Lipschitz-regular boundary surfaces, but that its application fails for problems defined on the exterior of open surfaces in three-dimensional spaces. We show that the boundary integral operators associated with the resolution of the Dirichlet and Neumann problems defined on the exterior of open surfaces have inverse operators, and that these operators would provide optimal preconditioners, but that they are not easily obtainable. Then we show a technique to explicitly obtain such inverse operators for the particular case when the open surface is the unit disk in a three-dimensional space. Their explicit inverse operators will be given, however, in the form of series, and will not be immediately applicable in the use of boundary element methods.In the third part we show how some modifications to these inverse operators allow us to obtain variational explicit and closed form expressions, no longer expressed as series, but also conserve nonetheless some characteristics that are relevant for their preconditioning effect. These explicit and closed forms expressions are applicable in boundary element methods. We obtain precise variational expressions for them and propose numerical schemes to compute the integrals needed for their use with boundary elements. The proposed numerical methods are tested using identities available within the developed theory and then used to build preconditioning matrices. Their performance as preconditioners for linear systems is tested for the case of the Laplace and Helmholtz problems for the unit disk. Finally, we propose an extension of this method that allows for the treatment of cases of open surfaces other that the disk. We exemplify and study this extension in its use on a square-shaped and an L-shaped surface screen in a three-dimensional space.

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

Cette thèse aborde deux thématiques différentes. On s’intéresse d’abord au développement et à l’analyse de méthodes pour le sens électrique appliqué à la robotique. On considère en particulier la méthode des réflexions permettant, à l’image de la méthode de Schwarz, de résoudre des problèmes linéaires à partir de sous-problèmes plus simples. Ces deniers sont obtenus par décomposition des frontières du problème de départ. Nous en présentons des preuves de convergence et des applications. Dans le but d’implémenter un simulateur du problème direct d’électrolocation dans un robot autonome, on s’intéresse également à une méthode de bases réduites pour obtenir des algorithmes peu coûteux en temps et en place mémoire. La seconde thématique traite d’un problème inverse dans le domaine de la chimie quantique. Nous cherchons ici à déterminer les caractéristiques d’un système quantique. Celui-ci est éclairé par un champ laser connu et fixé. Dans ce cadre, les données du problème inverse sont les états avant et après éclairage. Un résultat d’existence locale est présenté, ainsi que des méthodes de résolution numériques. / This thesis tackle two different topics.We first design and analyze algorithms related to the electrical sense for applications in robotics. We consider in particular the method of reflections, which allows, like the Schwartz method, to solve linear problems using simpler sub-problems. These ones are obtained by decomposing the boundaries of the original problem. We give proofs of convergence and applications. In order to implement an electrolocation simulator of the direct problem in an autonomous robot, we build a reduced basis method devoted to electrolocation problems. In this way, we obtain algorithms which satisfy the constraints of limited memory and time resources. The second topic is an inverse problem in quantum chemistry. Here, we want to determine some features of a quantum system. To this aim, the system is ligthed by a known and fixed Laser field. In this framework, the data of the inverse problem are the states before and after the Laser lighting. A local existence result is given, together with numerical methods for the solving.