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

This thesis presents how Tikhonov’s regularization can be used to solve an inverse problem of Helmholtz equation inside of a rectangle. The rectangle will be met with both Neumann and Dirichlet boundary conditions. A linear operator containing a Fourier series will be derived from the Helmholtz equation. Using this linear operator, an expression for the inverse operator can be formulated to solve the inverse problem. However, the inverse problem will be found to be ill-posed according to Hadamard’s definition. The regularization used to overcome this ill-posedness (in this thesis) is Tikhonov’s regularization. To compare the efficiency of this inverse operator with Tikhonov’s regularization, another inverse operator will be derived from Helmholtz equation in the partial frequency domain. The inverse operator from the frequency domain will also be regularized with Tikhonov’s regularization. Plots and error measurements will be given to understand how accurate the Tikhonov’s regularization is for both inverse operators. The main focus in this thesis is the inverse operator containing the Fourier series. A series of examples will also be given to strengthen the definitions, theorems and proofs that are made in this work.

Wind power is a rapidly increasing resource of electrical power world-wide. With the increasing number of wind turbines installed one major concern is the noise they generate. Sometimes already built wind turbines have to be put down or down-regulated, when certain noise levels are exceeded, resulting in economical and environmental losses. Therefore, accurate sound propagation calculations would be beneficial already in a planning stage of a wind farm. A model that can account for varying wind speeds and complex terrains could therefore be of great importance when future wind farms are planned. In this report an extended version of the classical wave equation that allows for variations in wind speed and terrain is derived which can be used to solve complex terrain and wind settings. The equation are solved with the use of Fourier transforms and Chebyshev polynomials and a numerical code is developed. The numerical code is evaluated against test cases where analytical and simple solutions exist. Tests with no wind for both totally free propagation and with a ground surface is evaluated in both 2D and 3D settings. For these simple cases the developed code shows good agreement to analytical solutions if the computational domain is sufficiently large. More advanced test cases with wind and terrain is not evaluated in this report and needs further validation. If the sound pressure needs to be calculated for a large area, and if the frequency is high, the developed model has problems regarding computational time and memory. These problems could be solved by further development of the numerical code or by using other solution methods.

A parallel domain decomposition Laplace transform Boundary Element Method, BEM, algorithm for the solution of large-scale transient heat conduction problems will be developed. This is accomplished by building on previous work by the author and including several new additions (most note-worthy is the extension to 3-D) aimed at extending the scope and improving the efficiency of this technique for large-scale problems. A Laplace transform method is utilized to avoid time marching and a Proper Orthogonal Decomposition, POD, interpolation scheme is used to improve the efficiency of the numerical Laplace inversion process. A detailed analysis of the Stehfest Transform (numerical Laplace inversion) is performed to help optimize the procedure for heat transfer problems. A domain decomposition process is described in detail and is used to significantly reduce the size of any single problem for the BEM, which greatly reduces the storage and computational burden of the BEM. The procedure is readily implemented in parallel and renders the BEM applicable to large-scale transient conduction problems on even modest computational platforms. A major benefit of the Laplace space approach described herein, is that it readily allows adaptation and integration of traditional BEM codes, as the resulting governing equations are time independent. This work includes the adaptation of two such traditional BEM codes for steady-state heat conduction, in both two and three dimensions. Verification and validation example problems are presented which show the accuracy and efficiency of the techniques. Additionally, comparisons to commercial Finite Volume Method results are shown to further prove the effectiveness.

[Truncated abstract] The study reported in this thesis extends the scaled boundary finite-element method to firstorder and second-order wave diffraction problems. The scaled boundary finite-element method is a newly developed semi-analytical technique to solve systems of partial differential equations. It works by employing a special local coordinate system, called scaled boundary coordinate system, to define the computational field, and then weakening the partial differential equation in the circumferential direction with the standard finite elements whilst keeping the equation strong in the radial direction, finally analytically solving the resulting system of equations, termed the scaled boundary finite-element equation. This unique feature of the scaled boundary finite-element method enables it to combine many of advantages of the finite-element method and the boundaryelement method with the features of its own. ... In this thesis, both first-order and second-order solutions of wave diffraction problems are presented in the context of scaled boundary finite-element analysis. In the first-order wave diffraction analysis, the boundary-value problems governed by the Laplace equation or by the Helmholtz equation are considered. The solution methods for bounded domains and unbounded domains are described in detail. The solution process is implemented and validated by practical numerical examples. The numerical examples examined include well benchmarked problems such as wave reflection and transmission by a single horizontal structure and by two structures with a small gap, wave radiation induced by oscillating bodies in heave, sway and roll motions, wave diffraction by vertical structures with circular, elliptical, rectangular cross sections and harbour oscillation problems. The numerical results are compared with the available analytical solutions, numerical solutions with other conventional numerical methods and experimental results to demonstrate the accuracy and efficiency of the scaled boundary finite-element method. The computed results show that the scaled boundary finite-element method is able to accurately model the singularity of velocity field near sharp corners and to satisfy the radiation condition with ease. It is worth nothing that the scaled boundary finite-element method is completely free of irregular frequency problem that the Green's function methods often suffer from. For the second-order wave diffraction problem, this thesis develops solution schemes for both monochromatic wave and bichromatic wave cases, based on the analytical expression of first-order solution in the radial direction. It is found that the scaled boundary finiteelement method can produce accurate results of second-order wave loads, due to its high accuracy in calculating the first-order velocity field.

La mise en œuvre de la méthode des éléments finis de frontière nécessite l'évaluation d'intégrales comportant un intégrand singulier. Un calcul fiable et précis de ces intégrales peut dans certains cas se révéler à la fois crucial et difficile. La méthode que nous proposons consiste en une réduction récursive de la dimension du domaine d'intégration et aboutit à une représentation de l'intégrale sous la forme d'une combinaison linéaire d'intégrales mono-dimensionnelles dont l'intégrand est régulier et qui peuvent s'évaluer numériquement mais aussi explicitement. L'équation de Helmholtz 3-D sert d'équation modèle mais ces résultats peuvent être utilisés pour les équations de Laplace et de Maxwell 3-D. L'intégrand est décomposé en une partie homogène et une partie régulière ; cette dernière peut être traitée par les méthodes usuelles d'intégration numérique. Pour la discrétisation du domaine, des triangles plans sont utilisés ; par conséquent, nous évaluons des intégrales sur le produit de deux triangles. La technique que nous avons développée nécessite de distinguer entre diverses configurations géométriques ; c'est pourquoi nous traitons séparément le cas de triangles coplanaires, dans des plans sécants ou parallèles. Divers prolongements significatifs de la méthode sont présentés : son extension à l'électromagnétisme, l'évaluation de l'intégrale du noyau de Green complet pour les coefficients d'auto-influence, et le calcul de la partie finie d'intégrales hypersingulières. / The implementation of the boundary element method requires the evaluation of integrals with a singular integrand. A reliable and accurate calculation of these integrals can in some cases be crucial and difficult. The proposed method is a recursive reduction of the dimension of the integration domain and leads to a representation of the integral as a linear combination of one-dimensional integrals whose integrand is regular and that can be evaluated numerically and even explicitly. The 3-D Helmholtz equation is used as a model equation, but these results can be used for the Laplace and the Maxwell equations in 3-D. The integrand is decomposed into a homogeneous part and a regular part, the latter can be treated by conventional numerical integration methods. For the discretization of the domain we use planar triangles, so we evaluate integrals over the product of two triangles. The technique we have developped requires to distinguish between several geometric configurations, that's why we treat separately the case of triangles in the same plane, in secant planes and in parallel planes.

Arbitrary high order numerical methods for time-harmonic acoustic scattering problems originally defined on unbounded domains are constructed. This is done by coupling recently developed high order local absorbing boundary conditions (ABCs) with finite difference methods for the Helmholtz equation. These ABCs are based on exact representations of the outgoing waves by means of farfield expansions. The finite difference methods, which are constructed from a deferred-correction (DC) technique, approximate the Helmholtz equation and the ABCs to any desired order. As a result, high order numerical methods with an overall order of convergence equal to the order of the DC schemes are obtained. A detailed construction of these DC finite difference schemes is presented. Details and results from an extension to heterogeneous media are also included. Additionally, a rigorous proof of the consistency of the DC schemes with the Helmholtz equation and the ABCs in polar coordinates is also given. The results of several numerical experiments corroborate the high order convergence of the proposed method. A novel local high order ABC for elastic waves based on farfield expansions is constructed and preliminary results applying it to elastic scattering problems are presented.

In this thesis, we explore a novel approach to solving the Helmholtz equation,the WaveHolz iteration method. This method aims to overcome some ofthe difficulties with solving the Helmholtz equation by providing a highlyparallelizable iterative method based on solving the time-dependent Waveequation. If this method proves reliable and computationally feasible it wouldhave great value for future application. Therefore, it is of interest to evaluatethe performance and properties of this method. To fully evaluate this method,the method was implemented and conclusions were based on results fromsimulations of the method. The method was able to solve problems in threedimensions and it seems that the method is very well suited for parallelized computations. To replicate real-world scenarios simulations of problems ininfinite and curvilinear domains were conducted. Based on the result presentedhere it is possible to further refine the method, especially regarding the setupof the domain and the implementation of boundary conditions for infinitedomains. / I detta examensarbete presenteras en ny metod för att lösa Helmholtz ekvation, WaveHoltz iterativa metod. Målet med denna metod är att undkomma vissa problem som uppstår med andra metoder för att lösa Helmholtz ekvation genom att tillhandahålla iterativ metod som baseras på lösningar av den tidsberoende vågekvationen samt kan parallelliseras effektivt. Om denna metod visar sig vara stabil och effektiv beräkningsmässigt skulle detta medföra stor potential för framtida tillämpningar. Av denna anledning undersöks metoden och dess egenskaper. För att fullt ut kunna evaluera denna method implementerades den vartefter simuleringar genomfördes och slutsatser drogs. Med metoden var att det var möjligt att lösa problem i tre dimensioner och metoden visade sig vara lämplig för parallella beräkningar. För att återskapa verklighetstrogna scenarion beräknades problem i oändliga och kroklinjiga domäner. Baserat på resultaten som presenteras i denna rapport är det möjligt att förfina metoden, speciellt vid konstruktionen av komplicerade beräkningsnät och randvillkoren för de oändliga problemen.

This thesis considers applications of gradient-based optimization algorithms to the design and control of some mechanics systems. The material distribution approach to topology optimization is applied to design two different acoustic devices, a reactive muffler and an acoustic horn, and optimization is used to control a ball pitching robot. Reactive mufflers are widely used to attenuate the exhaust noise of internal combustion engines by reflecting the acoustic energy back to the source. A material distribution optimization method is developed to design the layout of sound-hard material inside the expansion chamber of a reactive muffler. The objective is to minimize the acoustic energy at the muffler outlet. The presence or absence of material is represented by design variables that are mapped to varying coefficients in the governing equation. An anisotropic design filter is used to control the minimum thickness of materials separately in different directions. Numerical results demonstrate that the approach can produce mufflers with high transmission loss for a broad range of frequencies. For acoustic devices, it is possible to improve their performance, without adding extended volumes of materials, by an appropriate placement of thin structures with suitable material properties. We apply layout optimization of thin sound-hard material in the interior of an acoustic horn to improve its far-field directivity properties. Absence or presence of thin sound-hard material is modeled by a surface transmission impedance, and the optimization determines the distribution of materials along a “ground structure” in the form of a grid inside the horn. Horns provided with the optimized scatterers show a much improved angular coverage, compared to the initial configuration. The surface impedance is handled by a new finite element method developed for Helmholtz equation in the situation where an interface is embedded in the computational domain. A Nitschetype method, different from the standard one, weakly enforces the impedance conditions for transmission through the interface. As opposed to a standard finite-element discretization of the problem, our method seamlessly handles both vanishing and non-vanishing interface conditions. We show the stability of the method for a quite general class of surface impedance functions, provided that possible surface waves are sufficiently resolved by the mesh. The thesis also presents a method for optimal control of a two-link ball pitching robot with the aim of throwing a ball as far as possible. The pitching robot is connected to a motor via a non-linear torsional spring at the shoulder joint. Constraints on the motor torque, power, and angular velocity of the motor shaft are included in the model. The control problem is solved by an interior point method to determine the optimal motor torque profile and release position. Numerical experiments show the effectiveness of the method and the effect of the constraints on the performance.

For the purpose of numerically simulating continuum mechanical structures, different types of material may be represented by the extreme values {<img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />,1}, where 0&lt;<img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /><img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cll" />1, of a varying coefficient <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> in the governing equations. The paramter <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> is not allowed to vanish in order for the equations to be solvable, which means that the exact conditions are approximated. For example, for linear elasticity problems, presence of material is represented by the value <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> = 1, while <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> = <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> provides an approximation of void, meaning that material-free regions are approximated with a weak material. For acoustics applications, the value <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> = 1 corresponds to air and <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> = <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> to an approximation of sound-hard material using a dense fluid. Here we analyze the convergence properties of such material approximations as <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />!0, and we employ this type of approximations to perform design optimization. In Paper I, we carry out boundary shape optimization of an acoustic horn. We suggest a shape parameterization based on a local, discrete curvature combined with a fixed mesh that does not conform to the generated shapes. The values of the coefficient <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" />, which enters in the governing equation, are obtained by projecting the generated shapes onto the underlying computational mesh. The optimized horns are smooth and exhibit good transmission properties. Due to the choice of parameterization, the smoothness of the designs is achieved without imposing severe restrictions on the design variables. In Paper II, we analyze the convergence properties of a linear elasticity problem in which void is approximated by a weak material. We show that the error introduced by the weak material approximation, after a finite element discretization, is bounded by terms that scale as <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> and <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />1/2hs, where h is the mesh size and s depends on the order of the finite element basis functions. In addition, we show that the condition number of the system matrix scales inversely proportional to <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />, and we also construct a left preconditioner that yields a system matrix with a condition number independent of <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />. In Paper III, we observe that the standard sound-hard material approximation with <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Calpha" /> = <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> gives rise to ill-conditioned system matrices at certain wavenumbers due to resonances within the approximated sound-hard material. To cure this defect, we propose a stabilization scheme that makes the condition number of the system matrix independent of the wavenumber. In addition, we demonstrate that the stabilized formulation performs well in the context of design optimization of an acoustic waveguide transmission device. In Paper IV, we analyze the convergence properties of a wave propagation problem in which sound-hard material is approximated by a dense fluid. To avoid the occurrence of internal resonances, we generalize the stabilization scheme presented in Paper III. We show that the error between the solution obtained using the stabilized soundhard material approximation and the solution to the problem with exactly modeled sound-hard material is bounded proportionally to <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" />.