Spelling suggestions: "subject:"helmholtzequation"" "subject:"helmholtzekvation""
51 |
Métodos de Elementos Finitos e Diferenças Finitas para o Problema de Helmholtz / Finite Elements and Finite Difference Methods for the Helmholtz EquationFernandes, Daniel Thomas 02 March 2009 (has links)
Made available in DSpace on 2015-03-04T18:51:06Z (GMT). No. of bitstreams: 1
tese_danieltf.pdf: 1240547 bytes, checksum: d1fac8fed2c288c3581c57065cf2c0c2 (MD5)
Previous issue date: 2009-03-02 / Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / It is well known that classical finite elements or finite difference methods for Helmholtz problem present pollution effects that can severely deteriorate the quality of the approximate solution. To control pollution effects is especially difficult on non uniform meshes. For uniform meshes of square elements pollution effects can be minimized with the Quasi Stabilized Finite Element Method (QSFEM) proposed by Babus\v ska el al, for example. In the present work we initially present two relatively simple Petrov-Galerkin finite element methods, referred here as RPPG (Reduced Pollution Petrov-Galerkin) and QSPG (Quasi Stabilized Petrov-Galerkin), with reasonable robustness to some type of mesh distortion. The QSPG also shows minimal pollution, identical to QSFEM, for uniform meshes with square elements. Next we formulate the QOFD (Quasi Stabilized Finite Difference) method, a finite difference method for unstructured meshes. The QOFD shows great robustness relative to element distortion, but requires extra work to consider non-essential boundary conditions and source terms. Finally we present a Quasi Optimal Petrov-Galerkin (QOPG) finite element method. To formulate the QOPG we use the same approach introduced for the QOFD, leading to the same accuracy and robustness on distorted meshes, but constructed based on consistent variational formulation. Numerical results are presented illustrating the behavior of all methods developed compared to Galerkin, GLS and QSFEM. / É bem sabido que métodos clássicos de elementos finitos e diferenças finitas para o problema de Helmholtz apresentam efeito de poluição, que pode deteriorar seriamente a qualidade da solução aproximada. Controlar o efeito de poluição é especialmente difícil quando são utilizadas malhas não uniformes. Para malhas uniformes com elementos quadrados são conhecidos métodos (p. e. o QSFEM, proposto por Babuska et al) que minimizam a poluição. Neste trabalho apresentamos inicialmente dois métodos de elementos finitos de Petrov-Galerkin com formulação relativamente simples, o RPPG e o QSPG, ambos com razoável robustez para certos tipos de distorções dos elementos. O QSPG apresenta ainda poluição mínima para elementos quadrados. Em seguida é formulado o QOFD, um método de diferenças finitas aplicável a malhas não estruturadas. O QOFD apresenta grande robustez em relação a distorções, mas requer trabalho extra para tratar problemas não homogêneos ou condições de contorno não essenciais. Finalmente é apresentado um novo método de elementos finitos de Petrov-Galerkin, o QOPG, que é formulado aplicando a mesma técnica usada para obter a estabilização do QOFD, obtendo assim a mesma robustez em relação a distorções da malha, com a vantagem de ser um método variacionalmente consistente. Resultados numéricos são apresentados ilustrando o comportamento de todos os métodos desenvolvidos em comparação com os métodos de Galerkin, GLS e QSFEM.
|
52 |
Résolution numérique de quelques problèmes du type Helmholtz avec conditions au bord d'impédance ou des couches absorbantes (PML) / Numerical resolution of some Helmholtz-type problems with impedance boundary condition or PMLTomezyk, Jérôme 02 July 2019 (has links)
Dans cette thèse, nous étudions la convergence de méthode de type éléments finis pour les équations de Maxwell en régime harmonique avec condition au bord d'impédance et l'équation de Helmholtz avec une couche parfaitement absorbante(PML). On étudie en premier, la formulation régularisée de l'équation de Maxwell en régime harmonique avec condition au bord d'impédance (qui consiste à ajouter le term ∇ div à l'équation originale pour avoir un problème elliptique) et on garde la condition d'impédance comme une condition au bord essentielle. Pour des domaines à bord régulier, le caractère bien posé de cette formulation est bien connu mais cela n'est pas le cas pour des domaines polyédraux convexes. On commence alors le premier chapitre par la preuve du caractère bien posé dans le cas du polyèdre convexe, qui est basé sur le fait que l'espace variationnel est inclus dans H¹. Dans le but d'avoir des estimations explicites en le nombre d'onde k de ce problème, il est obligatoire d'avoir des résultats de stabilité explicites en ce nombre d'onde. C'est aussi proposé, pour quelques situations particulières, dans ce chapitre. Dans le second chapitre on décrit les singularités d'arêtes et de coins pour notre problème. On peut alors déduire la régularité de la solution du problème original, ainsi que de son adjoint. On a tous les ingrédients pour proposer une analyse de convergence explicite en k pour une méthode d'éléments finis avec éléments de Lagrange. Dans le troisième chapitre, on considère une méthode d'éléments finis hp non conforme pour un domaine à bord régulier. Pour obtenir des estimations explicites en k, on introduit un résultat de décomposition, qui sépare la solution du problème original (ou de son adjoint) en une partie régulière mais fortement oscillante et une partie moins régulière mais peu oscillante. Ce résultat permet de montrer des estimations explicites en k. Le dernier chapitre est dédié à l'équation de Helmholtz avec une PML. L'équation de Helmholtz dans l'espace entier est souvent utilisée pour modéliser la diffraction d'onde acoustique (en régime harmonique), avec la condition de radiation à l'infini de Sommerfeld. L'ajout d'une PML est une façon pour passer d'un domaine infini à un domaine fini, elle correspond à l'ajout d'une couche autour du domaine de calcul qui absorbe très vite toutes les ondes sortantes. On propose en premier un résultat de stabilité explicite en k. On propose alors deux schémas numériques, une méthode d'éléments finis hp et une méthode multi- échelle basée sur un sous-espace local de correction. Le résultat de stabilité est utilisé pour mettre en relation de choix des paramètres des méthodes numériques considérées avec k. Nous montrons aussi des estimations d'erreur a priori. A la fin de ces chapitres, des tests numériques sont proposés pour confirmer nos résultats théoriques. / In this thesis, we propose wavenumber explicit convergence analyses of some finite element methods for time-harmonic Maxwell's equations with impedance boundary condition and for the Helmholtz equation with Perfectly Matched Layer (PML). We first study the regularized formulation of time-harmonic Maxwell's equations with impedance boundary conditions (where we add a ∇ div-term to the original equation to have an elliptic problem) and keep the impedance boundary condition as an essential boundary condition. For a smooth domain, the wellposedness of this formulation is well-known. But the well-posedness for convex polyhedral domain has been not yet investigated. Hence, we start the first chapter with the proof of the well-posedness in this case, which is based on the fact that the variational space is embedded in H¹. In order to perform a wavenumber explicit error analysis of our problem, a wavenumber explicit stability estimate is mandatory. We then prove such an estimate for some particular configurations. In the second chapter, we describe the corner and edge singularities for such problem. Then we deduce the regularity of the solution of the original and the adjoint problem, thus we have all ingredients to propose a explicit wavenumber convergence analysis for h-FEM with Lagrange element. In the third chapter, we consider a non conforming hp-finite element approximation for domains with a smooth boundary. To perform a wavenumber explicit error analysis, we split the solution of the original problem (or its adjoint) into a regular but oscillating part and a rough component that behaves nicely for large frequencies. This result allows to prove convergence analysis for our FEM, again explicit in the wavenumber. The last chapter is dedicated to the Helmholtz equation with PML. The Helmholtz equation in full space is often used to model time harmonic acoustic scattering problems, with Sommerfeld radiation condition at infinity. Adding a PML is a way to reduce the infinite domain to a finite one. It corresponds to add an artificial absorbing layer surrounding a computational domain, in which scattered wave will decrease very quickly. We first propose a wavenumber explicit stability result for such problem. Then, we propose two numerical discretizations: an hp-FEM and a multiscale method based on local subspace correction. The stability result is used to relate the choice of the parameters in the numerical methods to the wavenumber. A priori error estimates are shown. At the end of each chapter, we perform numerical tests to confirm our theoritical results.
|
53 |
A Variable Resolution Global Spectral Method With Finer Resolution Over The TropicsJanakiraman, S 08 1900 (has links)
Variable resolution method helps to study the local scale phenomenon of interest within the context of global scale atmosphere/ocean dynamics. Global spectral methods based on spherical harmonics as basis functions are known to resolve a given function defined on the sphere, in an uniform manner. Though known for its mathematical elegance and higher order accuracy, global spectral methods are considered to be restrictive for developing mesh-refinement strategies. The only mesh refinement strategy available until now is due to the pioneering work of F. Schmidt. Schmidt transformation can study only one region with higher resolution.
The study of tropical dynamics is an interesting theme due to the presence of teleconnections between various phenomena, especially Indian Monsoon and the El-Nino. The Inter-Tropical Convergence Zone (ITCZ)is a continental scale phenamenon. It is in the ITCZ, many monsoon systems and tropical cyclones do occur. To study such phenomena under variable resolution method, high resolution is required in the entire tropical belt. Hitherto such a kind of mesh refinement strategies were not available in global spectral models. In this work, a new variable resolution method is developed that can help to study the tropical sub-scale phenomena with high resolution, in global spectral models.
A new conformal coordinate transformation named ’High resolution Tropical Belt Transformation(HTBT)’ is developed to generate high resolution in the entire tropical belt. Mathematical demonstrations are given to show the existence of additional conformal transformations available on the sphere, indicating additional degrees of freedom available to create variable resolution global spectral method.
Variable resolution global spectral method with high resolution over tropics is created through HTBT. The restriction imposed by Schmidt’s framework that the map-ping factor of the transformation need to have a finite-decomposition in the spectral space of the transformed domain is relaxed, by introduction of a new framework. The new framework uses transformed spherical harmonics Bnm as basis for spectral computations. With the use of FFT algorithm and Gaussian quadrature, the efficiency of the traditional spectral method is retained with the variable resolution global spectral method.
The newly defined basis functions Bnm are the eigenvalues of the transformed Laplacian operator . This property of Bnm provide an elegant direct solver for the transformed Helmholtz operator on the sphere. The transformed Helmholtz equations are solved accurately with the variable resolution method.
Advection experiments conducted with variable resolution spectral transport scheme on the HTBT variable grid produces near-dispersion free advection on the tropical belt. Transport across homogeneous resolution regions produce very less dispersion errors. Transport of a feature over the poles result in severe grid representation errors. It is shown that an increase in resolution around the poles greatly reduces this error. Transport of a feature from a point close to poles but not over it, does not produce such representation errors. Fourth-order Runge-Kutta scheme improves the accuracy of the transport scheme. The second order Magazenkov time-scheme proves to be better accurate than the leap-frog scheme with Asselin filter.
The non-divergent barotropic vorticity equation is tested with two exact solutions namely Rochas solution and Rossby-Haurwitz wave solutions. Each of the solution tests certain unique and contrasting characteristic of the system. The numerical behaviour of the solutions show non-linear interactions in them.
The singularity at the poles, arising due to the unbounded nature of the latitudinal derivative of the map factor of HTBT, triggers Gibbs phenomena for certain functions. However the recent advances in spectral methods, especially spectral viscosity method and Boyd-Vandeven filtering strategy provide ways to control the Gibbs oscillation and recover higher accuracy; make the variable resolution global spectral method viable for accurate meteorological computations.
|
54 |
Acoustic Wave Scattering From a Rough Seabed With a Continuously Varying Sediment Layer Overlying an Elastic BasementTsai, Sheng-Hsiung 01 August 2002 (has links)
Acoustic plane wave intearctions with a rough seabed with a continuously varying density and sound speed in a fluid-like sediment layer overlying an elastic basement is considered in this thesis. The acoustic properties in the sediment layer possess an exponential type of variation in density and one of the three classes of sound speed profiles, which are constant, k^2-linear, or inverse-square variations. Analytical solutions for the Helmholtz equation in the sediment layer, combined with a formulation based upon boundary perturbation theory, facilitate numerical implementation for the solution of coherent field.
The coherent reflection coefficients corresponding to the aformentioned density and sound speed profiles for various frequencies, roughness parameters, basement stiffness, are numerically generated and analyzed. Physical interpretations are provided for various results. This simple model characterizes three important features of an realistic sea floor, including seabed roughness, sediment inhomogenieties, and basement shear property,%Two dimensions is considered in the seafloor environment and the random roughness is belong to one dimension space.%
, therefore, provides a canonical model for the study of seabed acoustics.
The variation of the acoustic properties takes such a form that it is not only geologically realistic, but also renders analytical solutions for the Helmholtz equation, thus facilitating the formulation of the problem. The computational algorithm for the spatial spectrum of the scattered field due to random seabed has been developed based upon a boundary perturbation method. %About scattering field, only one time reflection from the sediment is taked account of, because the higher numerical order is, the lower scattering energy exist.%
The results have shown that, while the coherent field mainly depends upon the gross structure of the rough seabed represented by the RMS roughness, the scattered field heavily depends upon the details of the roughness structure specialized by the roughness power spectrum and the spatial correlation length of the rough surface. The dependence of the spatial spectrum on the sediment stratification is also carefully examined.
|
55 |
Modeling and design optimization of electromechanical brake actuator using eddy currentsKarakoc, Kerem 21 September 2012 (has links)
A novel electromechanical brake (EMB) based on the eddy current principle is proposed for application in electrical vehicles. The proposed solution is a feasible replacement for the current conventional hydraulic brake (CHB) systems. Unlike CHBs eddy current brakes (ECBs) use eddy currents and their interaction with an externally applied magnetic field to generate braking torque. Due to their pure electrically controllable and contact free nature, ECBs have multiple advantages over the current CHB systems, such as faster response, reduced weight and number of components, ease of implementing various controllers (e.g., anti-lock braking), and reduced noise levels. However, the torque generated by a typical ECB at low speeds is insufficient to effectively and completely stop a moving vehicle. Therefore, an ECB is commonly used as an assistive brake to the CHB system in heavy vehicles, i.e. trains and trucks In order to overcome this shortcoming, the use of AC magnetic fields is proposed to realize a stand-alone ECB system in which sufficient braking torque can be generated at low speeds. To this end, eddy currents are modeled analytically using the governing Maxwell’s equations with the consideration of time varying field application. The analytical model was validated using finite element analysis. Results show that the braking torque increases with the application of a time varying field.
Various forms of time varying fields have been studied. It was found that the
frequency-modulated applied field in triangular waveform results in the highest braking torque. Next, the design was optimized to maximize the braking torque and an optimum configuration was obtained using multiple pole projection areas (PPAs). Optimization results show that the braking torque significantly increases with the introduction of additional PPAs to the configuration, and the braking torque generation for an optimum four-PPA ECB configuration exceeds the braking requirements for current passenger
vehicles.
For control purposes, a dynamic model for a novel stand-alone ECB system using AC fields for automotive applications has been successfully designed and evaluated. Also, a model-based predictive controller has been developed for the optimum ECB
configuration. Finally an experimental test-bed has been designed for experimentation of both DC and AC field application on ECB. / Graduate
|
56 |
Scalar Waves In Spacetimes With Closed Timelike CurvesBugdayci, Necmi 01 December 2005 (has links) (PDF)
The existence and -if exists- the nature of the solutions of the scalar wave equation in spacetimes with closed timelike curves are investigated. The general properties of the solutions on some class of spacetimes are obtained.
Global monochromatic solutions of the scalar wave equation are obtained in flat wormholes of dimensions 2+1 and 3+1. The solutions are in the form of infinite series involving cylindirical and spherical wave functions and they are elucidated by the multiple scattering method. Explicit solutions for some
limiting cases are illustrated as well. The results of 2+1 dimensions are verified by using numerical methods.
|
57 |
Shape Optimization for Acoustic Wave Propagation ProblemsUdawalpola, Rajitha January 2010 (has links)
Boundary shape optimization is a technique to search for an optimal shape by modifying the boundary of a device with a pre-specified topology. We consider boundary shape optimization of acoustic horns in loudspeakers and brass wind instruments. A horn is an interfacial device, situated between a source, such as a waveguide or a transducer, and surrounding space. Horns are used to control both the transmission properties from the source and the spatial power distribution in the far-field (directivity patterns). Transmission and directivity properties of a horn are sensitive to the shape of the horn flare. By changing the horn flare we design transmission efficient horns. However, it is difficult to achieve both controllability of directivity patterns and high transmission efficiency by using only changes in the horn flare. Therefore we use simultaneous shape and so-called topology optimization to design a horn/acoustic-lens combination to achieve high transmission efficiency and even directivity. We also design transmission efficient interfacial devices without imposing an upper constraint on the mouth diameter. The results demonstrate that there appears to be a natural limit on the optimal mouth diameter. We optimize brasswind instruments with respect to its intonation properties. The instrument is modeled using a hybrid method between a one-dimensional transmission line analogy for the slowly flaring part of the instrument, and a finite element model for the rapidly flaring part. An experimental study is carried out to verify the transmission properties of optimized horn. We produce a prototype of an optimized horn and then measure the input impedance of the horn. The measured values agree reasonably well with the predicted optimal values. The finite element method and the boundary element method are used as discretization methods in the thesis. Gradient-based optimization methods are used for optimization, in which the gradients are supplied by the adjoint methods.
|
58 |
A study on block flexible iterative solvers with applications to Earth imaging problem in geophysics / Étude de méthodes itératives par bloc avec application à l’imagerie sismique en géophysiqueFerreira Lago, Rafael 13 June 2013 (has links)
Les travaux de ce doctorat concernent le développement de méthodes itératives pour la résolution de systèmes linéaires creux de grande taille comportant de nombreux seconds membres. L’application visée est la résolution d’un problème inverse en géophysique visant à reconstruire la vitesse de propagation des ondes dans le sous-sol terrestre. Lorsque de nombreuses sources émettrices sont utilisées, ce problème inverse nécessite la résolution de systèmes linéaires complexes non symétriques non hermitiens comportant des milliers de seconds membres. Dans le cas tridimensionnel ces systèmes linéaires sont reconnus comme difficiles à résoudre plus particulièrement lorsque des fréquences élevées sont considérées. Le principal objectif de cette thèse est donc d’étendre les développements existants concernant les méthodes de Krylov par bloc. Nous étudions plus particulièrement les techniques de déflation dans le cas multiples seconds membres et recyclage de sous-espace dans le cas simple second membre. Des gains substantiels sont obtenus en terme de temps de calcul par rapport aux méthodes existantes sur des applications réalistes dans un environnement parallèle distribué. / This PhD thesis concerns the development of flexible Krylov subspace iterative solvers for the solution of large sparse linear systems of equations with multiple right-hand sides. Our target application is the solution of the acoustic full waveform inversion problem in geophysics associated with the phenomena of wave propagation through an heterogeneous model simulating the subsurface of Earth. When multiple wave sources are being used, this problem gives raise to large sparse complex non-Hermitian and nonsymmetric linear systems with thousands of right-hand sides. Specially in the three-dimensional case and at high frequencies, this problem is known to be difficult. The purpose of this thesis is to develop a flexible block Krylov iterative method which extends and improves techniques already available in the current literature to the multiple right-hand sides scenario. We exploit the relations between each right-hand side to accelerate the convergence of the overall iterative method. We study both block deflation and single right-hand side subspace recycling techniques obtaining substantial gains in terms of computational time when compared to other strategies published in the literature, on realistic applications performed in a parallel environment.
|
59 |
Métodos de Elementos Finitos e Diferenças Finitas para o Problema de Helmholtz / Finite Elements and Finite Difference Methods for the Helmholtz EquationDaniel Thomas Fernandes 02 March 2009 (has links)
É bem sabido que métodos clássicos de elementos finitos e diferenças finitas para o problema de Helmholtz apresentam efeito de poluição, que pode deteriorar seriamente a qualidade da solução aproximada. Controlar o efeito de poluição é especialmente difícil quando são utilizadas malhas não uniformes. Para malhas uniformes com elementos quadrados são conhecidos métodos (p. e. o QSFEM, proposto por Babuska et al) que minimizam a poluição. Neste trabalho apresentamos inicialmente dois métodos de elementos finitos de Petrov-Galerkin com formulação relativamente simples, o RPPG e o QSPG, ambos com razoável robustez para certos tipos de distorções dos elementos. O QSPG apresenta ainda poluição mínima para elementos quadrados. Em seguida é formulado o QOFD, um método de diferenças finitas aplicável a malhas não estruturadas. O QOFD apresenta grande robustez em relação a distorções, mas requer trabalho extra para tratar problemas não homogêneos ou condições de contorno não essenciais. Finalmente é apresentado um novo método de elementos finitos de Petrov-Galerkin, o QOPG, que é formulado aplicando a mesma técnica usada para obter a estabilização do QOFD, obtendo assim a mesma robustez em relação a distorções da malha, com a vantagem de ser um método variacionalmente consistente. Resultados numéricos são apresentados ilustrando o comportamento de todos os métodos desenvolvidos em comparação com os métodos de Galerkin, GLS e QSFEM. / It is well known that classical finite elements or finite difference methods for Helmholtz problem present pollution effects that can severely deteriorate the quality of the approximate solution. To control pollution effects is especially difficult on non uniform meshes. For uniform meshes of square elements pollution effects can be minimized with the Quasi Stabilized Finite Element Method (QSFEM) proposed by Babusv ska el al, for example. In the present work we initially present two relatively simple Petrov-Galerkin finite element methods, referred here as RPPG (Reduced Pollution Petrov-Galerkin) and QSPG (Quasi Stabilized Petrov-Galerkin), with reasonable robustness to some type of mesh distortion. The QSPG also shows minimal pollution, identical to QSFEM, for uniform meshes with square elements. Next we formulate the QOFD (Quasi Stabilized Finite Difference) method, a finite difference method for unstructured meshes. The QOFD shows great robustness relative to element distortion, but requires extra work to consider non-essential boundary conditions and source terms. Finally we present a Quasi Optimal Petrov-Galerkin (QOPG) finite element method. To formulate the QOPG we use the same approach introduced for the QOFD, leading to the same accuracy and robustness on distorted meshes, but constructed based on consistent variational formulation. Numerical results are presented illustrating the behavior of all methods developed compared to Galerkin, GLS and QSFEM.
|
60 |
Development of a reference method based on the fast multipole boundary element method for sound propagation problems in urban environments : formalism, improvements & applications / Développement d’une méthode de référence basée sur la méthode par éléments de frontières multipolaires pour la propagation sonore en environnement urbain : formalisme, optimisations & applicationsVuylsteke, Xavier 10 December 2014 (has links)
Décrit comme l'un des algorithmes les plus prometteurs du 20ème siècle, le formalisme multipolaire appliqué à la méthode des éléments de frontière, permet de nos jours de traiter de larges problèmes encore inconcevables il y a quelques années. La motivation de ce travail de thèse est d'évaluer la capacité, ainsi que les avantages concernant les ressources numériques, de ce formalisme pour apporter une solution de référence aux problèmes de propagation sonore tri-dimensionnels en environnement urbain, dans l'objectif d'améliorer les algorithmes plus rapides déjà existants. Nous présentons la théorie nécessaire à l'obtention de l'équation intégrale de frontière pour la résolution de problèmes non bornés. Nous discutons également de l'équation intégrale de frontière conventionnelle et hyper-singulière pour traiter les artefacts numériques liés aux fréquences fictives, lorsque l'on résout des problèmes extérieurs. Nous présentons par la suite un bref aperçu historique et technique du formalisme multipolaire rapide et des outils mathématiques requis pour représenter la solution élémentaire de l'équation de Helmholtz. Nous décrivons les principales étapes, d'un point de vue numérique, du calcul multipolaire. Un problème de propagation sonore dans un quartier, composé de 5 bâtiments, nous a permis de mettre en évidence des problèmes d'instabilités dans le calcul par récursion des matrices de translations, se traduisant par des discontinuités sur le champs de pression de surface et une non convergence du solveur. Ceci nous a conduits à considérer le travail très récent de Gumerov et Duraiswamy en lien avec un processus récursif stable pour le calcul des coefficients des matrices de rotation. Cette version améliorée a ensuite été testée avec succès sur un cas de multi diffraction jusqu'à une taille dimensionnelle de problème de 207 longueur d'ondes. Nous effectuons finalement une comparaison entre un algorithme d'élément de frontière, Micado3D, un algorithme multipolaire et un algorithme basé sur le tir de rayons, Icare, pour le calcul de niveaux de pression moyennés dans une cour ouverte et fermée. L'algorithme multipolaire permet de valider les résultats obtenus par tir de rayons dans la cour ouverte jusqu'à 300 Hz (i.e. 100 longueur d'ondes), tandis que concernant la cour fermée, zone très sensible par l'absence de contribution directes ou réfléchies, des études complémentaires sur le préconditionnement de la matrice semblent requises afin de s'assurer de la pertinence des résultats obtenus à l'aide de solveurs itératifs / Described as one of the best ten algorithms of the 20th century, the fast multipole formalism applied to the boundary element method allows to handle large problems which were inconceivable only a few years ago. Thus, the motivation of the present work is to assess the ability, as well as the benefits in term of computational resources provided by the application of this formalism to the boundary element method, for solving sound propagation problems and providing reference solutions, in three dimensional dense urban environments, in the aim of assessing or improving fast engineering tools. We first introduce the mathematical background required for the derivation of the boundary integral equation, for solving sound propagation problems in unbounded domains. We discuss the conventional and hyper-singular boundary integral equation to overcome the numerical artifact of fictitious eigen-frequencies, when solving exterior problems. We then make a brief historical and technical overview of the fast multipole principle and introduce the mathematical tools required to expand the elementary solution of the Helmholtz equation and describe the main steps, from a numerical viewpoint, of fast multipole calculations. A sound propagation problem in a city block made of 5 buildings allows us to highlight instabilities in the recursive computation of translation matrices, resulting in discontinuities of the surface pressure and a no convergence of the iterative solver. This observation leads us to consider the very recent work of Gumerov & Duraiswamy, related to a ``stable'' recursive computation of rotation matrices coefficients in the RCR decomposition. This new improved algorithm has been subsequently assessed successfully on a multi scattering problem up to a dimensionless domain size equal to 207 wavelengths. We finally performed comparisons between a BEM algorithm, extit{Micado3D}, the FMBEM algorithm and a ray tracing algorithm, Icare, for the calculation of averaged pressure levels in an opened and closed court yards. The fast multipole algorithm allowed to validate the results computed with Icare in the opened court yard up to 300 Hz corresponding, (i.e. 100 wavelengths), while in the closed court yard, a very sensitive area without direct or reflective fields, further investigations related to the preconditioning seem required to ensure reliable solutions provided by iterative solver based algorithms
|
Page generated in 0.0772 seconds