Spelling suggestions: "subject:"équations""
401 |
Inégalités de type Trudinger-Moser et applications / Trudinger-Moser type inequalities and applicationsZghal, Mohamed Khalil 06 February 2016 (has links)
Cette thèse porte sur quelques inégalités de type Trudinger-Moser et leurs applications à l'étude des injections de Sobolev qu'elles induisent dans les espaces d'Orlicz et à l'analyse d'équations aux dérivées partielles non linéaires à croissance exponentielle.Le travail qu'on présente ici se compose de trois parties. La première partie est consacrée à la description du défaut de compacité de l'injection de Sobolev 4D dans l'espace d'Orlicz dansle cadre radial.L'objectif de la deuxième partie est double. D'abord, on caractérise le défaut de compacité de l'injection de Sobolev 2D dans les différentes classes d'espaces d'Orlicz. Ensuite, on étudiel'équation de Klein-Gordon semi-linéaire avec non linéarité exponentielle, où la norme d'Orlicz joue un rôle crucial. En particulier, on aborde les questions d'existence globale, de complétude asymptotique et d'étude qualitative.Dans la troisième partie, on établit des inégalités optimales de type Adams, en étroite relation avec les inégalités de Hardy, puis on fournit une description du défaut de compacité des injections de Sobolev qu'elles induisent / This thesis focuses on some Trudinger-Moser type inequalities and their applications to the study of Sobolev embeddings they induce into the Orlicz spaces, and the investigation of nonlinear partial differential equations with exponential growth.The work presented here includes three parts. The first part is devoted to the description of the lack of compactness of the 4D Sobolev embedding into the Orlicz space in the radialframework.The aim of the second part is twofold. Firstly, we characterize the lack of compactness of the 2D Sobolev embedding into the different classes of Orlicz spaces. Secondly, we undertakethe study of the nonlinear Klein-Gordon equation with exponential growth, where the Orlicz norm plays a crucial role. In particular, issues of global existence, scattering and qualitativestudy are investigated.In the third part, we establish sharp Adams-type inequalities invoking Hardy inequalities, then we give a description of the lack of compactness of the Sobolev embeddings they induce
|
402 |
A quantum approach to dynamical quarkonia suppression in high energy heavy ion collisions / Une approche quantique de la suppression dynamique des quarkonia dans les collisions d’ions lourds à haute énergieKatz, Roland 14 December 2015 (has links)
La chromodynamique quantique (QCD) prédit l'existence d'un nouvel état de la matière : le plasma de quarks et de gluons (PQG). Celui-ci aurait existé dans les premiers instants suivant le Big Bang et peut en principe être produit sous les conditions extrêmes de température et de densité atteintes lors de collisions d'ions lourds à haute énergie (au LHC par exemple). Un des marqueurs de sa présence est la suppression des quarkonia (états liés de quark/antiquark lourds), caractérisée par une production inférieure de ces états dans les collisions d'ions lourds relativement aux collisions proton-proton où le PQG ne pourrait être créé. Cette suppression a bien été observée expérimentalement, mais l'évolution de ses tendances aux énergies du RHIC et du LHC est un véritable défi qui requiert une meilleure compréhension théorique. La présente thèse a pour but d’étudier l’évolution en temps réel de paires corrélées de quark/antiquark lourds considérées comme des systèmes quantiques ouverts en interaction permanente avec un PQG en refroidissement. Explicitement, l'interaction continue entre le milieu et les degrés de liberté internes de la paire est obtenue par 1) un écrantage de couleur dit « de Debye » dû à la présence de charges de couleur dans leur voisinage et 2) des mécanismes de fluctuation/dissipation qui reflètent les collisions permanentes. Cela mène à une image dynamique et continue de la dissociation des quarkonia, de leur recombinaison et des transitions entre états liés. L'étude est transversale à différents cadres théoriques : semi-classique, quantique et quantique des champs. Les prédictions du modèle sont comparées aux résultats expérimentaux et aux résultats d'autres modèles théoriques. / The theory of quantum chromodynamics (QCD) predicts the existence of a new state of matter: the Quark-GluonPlasma (QGP). The latter may have existed at the first moments of the Universe following the Big Bang and can be, in theory, re-produced under the extreme conditions of temperature and density reached in high energy heavy ion collisions (at the LHC for instance). One of the QGP observables is the suppression of the quarkonia (heavy quark/antiquark bound states), characterised by a smaller production of these states in heavy ion collisions in comparison to proton-protoncollisions, in which no QGP production would be possible. This suppression has indeed been observed experimentally, but the puzzling evolution of its trend from RHIC to LHC energies requires a better theoretical understanding. The present thesis aims at studying the real-time evolution of correlated heavy quark/antiquark pairs described as open quantum systems which permanently interact with a cooling QGP. More explicitly, the continuous interaction between the medium and the pair internal degrees of freedom is obtained through 1) a temperature dependent color screening (“Debye” like) due to color charges in the irvicinity and 2) some fluctuation/dissipation mechanisms reflecting the continuous collisions. It leads to a dynamical and continuous picture of the dissociation, recombination and possible transitions to other bound states. This investigation is at the crossroads of different theoretical frameworks: semi-classic, quantum and quantum fields. The deduced predictions are compared to experimental data and to the results of other theoretical models.
|
403 |
Contribution à la théorie des EDP non linéaires avec applications à la méthode des surfaces de niveau, aux fluides non newtoniens et à l'équation de Boltzmann / A contribution to non-linear PDEs with applications to the level set method, non-Newtonian fluid flows and the Boltzmann equationNtovoris, Eleftherios 12 September 2016 (has links)
Cette thèse comporte trois chapitres indépendants, consacrés à l’étude mathématique de trois problèmes physiques distincts, ayant pour modèles trois équations aux dérivées partielles différentes. Ces équations relèvent plus précisément de la méthode des surfaces de niveau, de la théorie de l’écoulement incompressible des matériaux non newtoniens et de la théorie cinétique des gaz raréfiés. Le premier chapitre de la thèse porte sur la dynamique des frontières en mouvement et contient une justification mathématique de la procédure numérique dite de ré-initialisation, dont les applications sont nombreuses dans le contexte de la célèbre méthode des surfaces de niveau. Nous appliquons ces résultats pour une classe d’équations issues de la méthode des surfaces de niveau de premier ordre. Nous écrivons la procédure de ré-initialisation comme un algorithme de décomposition et nous étudions la convergence de l’algorithme en utilisant des techniques d’homogénéisation dans la variable temporelle. Grâce à cette analyse rigoureuse nous introduisons également une nouvelle méthode pour l’approximation de la fonction de distance dans le contexte de la méthode des surfaces de niveau. Dans le cas où l’on cherche seulement une fonction de l’ensemble de niveau avec un gradient minoré proche du niveau zéro, nous proposons une approximation plus simple. Dans le cas général, où le niveau zéro pourrait présenter des changements de topologie, nous introduisons une nouvelle notion de limites relâchées. Dans le deuxième chapitre de la thèse, nous étudions un problème de frontière libre résultant de l’étude de l’écoulement incompressible d’un matériau non-newtonien, avec limite d’élasticité de type Drucker-Prager, sur un plan incliné et sous l’effet de la pesanteur. Nous obtenons une équation sous-différentielle, que nous formulons comme un problème variationnel avec un terme à croissance linéaire de type gradient, et nous étudions le problème dans un domaine non borné. Nous montrons que les équations sont bien posées et satisfont certaines propriétés de régularité. Nous sommes alors capables de relier les paramètres physiques avec le problème abstrait et de prouver des propriétés quantitatives de la solution. En particulier, nous montrons que la solution a un support compact, la limite de ce que nous appelons la frontière libre. Nous construisons également des solutions explicites d’une équation différentielle ordinaire qui peut estimer la frontière libre. Enfin, le troisième et dernier chapitre de la thèse est dédié aux solutions de l’équation de Boltzmann homogène avec molécules maxwelliennes et énergie infinie. Nous obtenons de nouveaux résultats d’existence de solutions éternelles pour cette équation dans un espace de mesures de probabilité d’énergie infinie (i.e. de moment d’ordre deux infini). Elles permettent de décrire le comportement asymptotique en temps d’autres solutions d’énergie infinie, mais elles apparaissent aussi comme des états asymptotiques intermédiaires dans l’étude des solutions d’énergie finie, mais arbitrairement grande. Les méthodes issues de l’analyse harmonique sont utilisées pour étudier l’équation de Boltzmann, où la variable de vitesse est exprimée en Fourier. Enfin, un changement d’échelle logarithmique en la variable temporelle permet de déterminer le bon comportement asymptotique à l’infini des solutions / This thesis consists of three different and independent chapters, concerning the mathematical study of three distinctive physical problems, which are modelled by three non- linear partial differential equations. These equations concern the level set method, the theory of incompressible flow of non-Newtonian materials and the kinetic theory of rare- fied gases. The first chapter of the thesis concerns the dynamics of moving interfaces and contains a rigorous justification of a numerical procedure called re-initialization, for which there are several applications in the context of the level set method. We apply these results for first order level set equations. We write the re-initialization procedure as a splitting algorithm and study the convergence of the algorithm using homogenization techniques in the time variable. As a result of the rigorous analysis, we are also able to introduce a new method for the approximation of the distance function in the context of the level set method. In the case where one only looks for a level set function with gradient bounded from below near the zero level, we propose a simpler approximation. In the general case where the zero level might present changes of topology we introduce a new notion of relaxed limits. In the second chapter of the thesis, we study a free boundary problem arising in the study of the flow of an incompressible non-Newtonian material with Drucker-Prager plasticity on an inclined plane. We derive a subdifferential equation, which we reformulate as a variational problem containing a term with linear growth in the gradient variable, and we study the problem in an unbounded domain. We show that the equations are well posed and satisfy some regularity properties. We are then able to connect the physical parameters with the abstract problem and prove some quantitative properties of the solution. In particular, we show that the solution has compact support and the support is the free boundary. We also construct explicit solutions of an ordinary differential equation, which we use to estimate the free boundary. The last chapter of the thesis is dedicated to the study of infinite energy solutions of the homogeneous Boltzmann equation with Maxwellian molecules. We obtain new results concerning the existence of eternal solutions in the space of probability measure with infinite energy (i.e. the second order moment is infinite). These solutions describe the asymptotic behaviour of other infinite energy solutions but could also be useful in the study of intermediate asymptotic states of solutions with finite but arbitrarily large energy. We use harmonic analysis tools to study the equation, where the velocity variable is expressed in the Fourier space. Finally, a logarithmic scaling of the time variable allows to determine the correct asymptotic scaling of the solutions
|
404 |
Développement d'une approche numérique et expérimentale pour un conduit avec traitement acoustique : application à la validation de modèles d'impédance en propagation multimodale avec écoulement / Development of a numerical and experimental approach for a duct with an acoustic treatment : application for validation of impedance models with multimodal propagation and mean flowBaccouche, Ryan 01 February 2016 (has links)
La prédiction de la réduction passive du bruit nécessite l’utilisation de modèles d’impédance acoustique fiables. La proposition de tels modèles en présence d’écoulement est un problème qui reste ouvert. Les travaux réalisés dans cette thèse sont consacrés à la mise en place d’une méthode expérimentale et numérique permettant de tester la validité de modèles d’impédance acoustique de traitements SDOF en présence d’écoulement. Cette approche repose sur la comparaison de résultats expérimentaux issus d’un banc aéroacoustique, et de résultats numériques issus d’un modèle FEM-PML axisymétrique basé sur l’équation de Galbrun. Ainsi, dans la première partie de cette thèse, le modèle FEM-PML utilisé est présenté. Dans la seconde partie, des modèles d’impédance acoustique de SDOF tenant compte de l’écoulement sont exposés, ainsi que la méthode de validation de ces modèles et les résultats de ces validations. / Prediction of noise reduction requires reliable acoustic impedance models. The proposal of such models in the presence of a mean flow is a problem which remains open. This thesis is devoted to the development of an experimental and a numerical method to test the validity of SDOF treatments acoustic impedance models in the presence of mean flow. This approach is based on the comparison of experimental results from an aeroacoustical duct, and numerical results from a FEM-PML axisymmetric model based on Galbrun’s equation. Thus, in the first part of this thesis, the FEM-PML model is presented. In the second part, some acoustic impedance models are exposed, and the validation method and the results of these validations are presented.
|
405 |
Efficient domain decomposition methods for electromagnetic modeling of scattering from forest environments / Méthodes de décomposition de domaine pour la modélisation électromagnétique de la diffusion par la forêtFenni, Ines 23 October 2014 (has links)
Cette thèse porte sur la modélisation et l’analyse de la propagation électromagnétique dans un milieu forestier dans les bandes VHF et UHF. L’objectif principal est le développement d’un modèle numérique "full-wave" tridimensionnel de diffusion par la forêt permettant de caractériser l’interaction d’une onde électromagnétique avec un milieu forestier. Un tel modèle s’avère, actuellement, un outil indispensable à l’analyse des mesures radar pour l’étude des paramètres caractéristiques de la forêt tels que la biomasse forestière, la hauteur des arbres et leur densité. La complexité numérique de ce modèle a limité son domaine d’application à de petites parcelles de forêt et aux basses fréquences. Pour pouvoir traiter de larges zones forestières tout en montant en fréquence, et s’approcher ainsi des besoins et exigences des utilisateurs potentiels de notre modèle, nous avons intégré à ce modèle une méthode numérique efficace dédiée à l’analyse de larges problèmes électromagnétiques. La méthode en question,connue sous le nom de Characteristic Basis Function Method (CBFM) était récemment développée dans le laboratoire de Communication et Electromagnétisme de l’université PennState dirigé par le Professeur Mittra. Après une optimisation et une adaptation au problème d’intérêt, la CBFM réalise d’excellentes performances et nous permet une diminution considérable du temps de calcul et des besoins en espace mémoire sans pour autant dégrader la qualité des résultats obtenus ou altérer la fidélité du modèle à la réalité du problème électromagnétique traité. / A 3-D full-wave model, based on the integral representation of the electric field and dedicatedto the analysis of bi-static scattering mechanisms by a forest in the VHF and UHF bands wasefficiently enhanced. In order to overcome the limitation of a previous 3D model to small simulationsscenes and low frequencies, we have developed, during this research work, a new model using basis functionsadapted to the problem of interest, in the context of the Characteristic Basis Function Method(CBFM) and we investigated the suitability of this direct method for computing the electromagneticfields inside and outside three-dimensional dielectric objects representing the tree trunks and branches.The CBFM has shown great performances, when applied to the forest scattering modeling, both interms of CPU time and memory resources needed. Once properly set, the CBFM-E is so efficient thatit is able to treat in few minutes electromagnetic problems totally intractable with the classical MoM.Consequently, we have developed a powerful 3D forest electromagnetic scattering tool which allows ustoday to compute large forest electromagnetic problems in few minutes without worrying about theaccuracy of the solution. On the other hand, we have demonstrated the efficiency and accuracy of theCBFM-E when applied to 3D dielectric objects in the context of the electric volumetric integral equation,and have consolidated thus its leading position in the computational electromagnetics, especiallyagainst the iterative solvers based numerical methods.
|
406 |
Application de l'équation PC-SAFT à la capture du dioxyde de carbone et à la désulfuration des essences / Application of PC-SAFT equation of state in the carbon dioxide capture and the gasoline desulfurizationChen, Yushu 26 September 2013 (has links)
Le remplacement des solvants organiques classiques par une nouvelle génération de solvants moins toxiques, moins inflammables et moins polluants est un défi majeur pour l'industrie chimique. Les liquides ioniques, sels liquides qui satisfont ces critères, sont envisagés comme alternatives. Le but de ce travail est d'évaluer le comportement des liquides ioniques en présence de gaz à effet de serre (CO2, N2O et CH4) ou de composés organiques. Dans un premier temps, une étude théorique présente les performances du modèle thermodynamique PC-SAFT sur la représentation des équilibres liquide-vapeur de systèmes constitués de dioxyde de carbone et de liquide ionique. Ensuite, l'étude de la solubilité du méthane, du dioxyde de carbone et du protoxyde d'azote dans divers liquides ioniques a été effectuée sous basse ou haute pression. Ce travail propose un modèle basé sur le concept de contribution de groupes afin de prédire la constante d'Henry du CO2 dans les liquides ioniques. Enfin, une étude sur les équilibres liquide-vapeur des systèmes binaires rencontrés dans la désulfuration des essences a été effectuée. Les données expérimentales ont permis d'évaluer les performances de l'équation PC-SAFT à représenter les équilibres entre phases de systèmes {composés soufrés / aromatiques + liquide ionique} / The replacement of conventional organic solvents by a new generation of solvents less toxic, less flammable and less polluting is a major challenge for the chemical industry. Ionic liquids have been widely promoted as interesting substitutes for traditional solvents. The purpose of this work is to evaluate the behavior of ionic liquids in the presence of greenhouse gases (CO2, CH4 and N2O) or organic compounds. Firstly, a theoretical study presents the performance of the thermodynamic model PC-SAFT in the representation of vapor-liquid equilibrium of systems containing ionic liquids and carbon dioxide. Then, the solubility study of methane, carbon dioxide and nitrous oxide in various ionic liquids was performed at high or low pressure. The group contribution concept is proposed in this study in order to predict the Henry's law constant of carbon dioxide in ionic liquids. Finally, a study on the vapor-liquid equilibrium of binary systems encountered in gasoline desulfurization was carried out. Experimental data were used to evaluate the performance of PC-SAFT equation of state to represent phase equilibrium of systems {sulfur / aromatic compounds + ionic liquid}
|
407 |
Processus auto-stabilisants dans un paysage multi-puits / Self-stabilizing processes in a multi-well landscapeTugaut, Julian 06 July 2010 (has links)
Les processus auto-stabilisants sont définis comme des solutions d'équations différentielles stochastiques dont le terme de dérive contient à la fois le gradient d'un potentiel ainsi qu'un terme non-linéaire au sens de McKean qui attire le processus vers sa propre loi de distribution. On dispose de nombreux résultats lorsque l'environnement est convexe. L'objet de ce travail est de les étendre autant que possible au cas général notamment lorsque le paysage contient plusieurs puits. Des différences fondamentales sont constatées.Le premier chapitre prouve l'existence d'une solution forte. Le second s'intéresse aux lois de probabilités d'une telle solution. En particulier, l'existence et la non-unicité des mesures stationnaires sont mises en évidence sous des hypothèses faibles. Les chapitres trois et quatre sont affectés au comportement de ces mesures lorsque le coefficient de diffusion tend vers 0.Le chapitre cinq met en relation le processus auto-stabilisant avec des systèmes particulaires via une « propagation du chaos ». Il est ainsi possible de transposer certains résultats du système de particules sur le processus non-markovien et réciproquement. Le chapitre six est dédié au dénombrement exact des mesures stationnaires.Le chapitre sept est employé pour l'étude du comportement en temps long. D'une part, un résultat de convergence dans un cas simple est fourni. D'autre part, un principe de grandes déviations est mis en évidence par l'utilisation des résultats de Freidlin et Wentzell / Self-stabilizing processes are defined as the solutions of stochastic differential equations which drift term contains the gradient of a potential and a term nonlinear in the sense of McKean which attracts the process to its own law distribution. There are many results if the landscape is convex. The purpose of this work is to extend these in the general case especially when the landscape contains contains several wells. Essential differences are found.The first chapter proves the strong existence of a solution. The second one deals with the probability measure of the solution. Particularly, the existence and the non-uniqueness of the stationary measures are highlighted under weak assumptions. Chapter three and four are assigned to the asymptotic analysis in the small noise limit of these measures.Chapter five connects the self-stabilizing process and some particle systems via a « propagation of chaos ». It is thus possible to translate some results from the particle systems to the non-markovian process and reciprocally.Chapter seven is used to study the long time behavior. In one hand, a convergence's result is provided in a simple case. In the other hand, a large deviations principle is highlighted by using the results of Freidlin and Wentzell
|
408 |
Prédiction du comportement de phases et des enthalpies de mélange de gaz naturels atypiques contenant de l'argon, du monoxyde de carbone et de l'hélium / Prediction of phase behaviour and enthalpies of mixing of atypical natural gases containing argon, carbon monoxide and heliumPlée, Vincent 17 December 2014 (has links)
Le développement du modèle prédictif E-PPR78, basé sur une méthode de contribution de groupe, a été entrepris depuis plus de dix ans pour prédire le comportement de systèmes multiconstituants. Ce modèle repose sur l'équation d'état de Peng-Robinson dans sa version de 1978 et les règles de mélanges de Van der Waals. Il utilise un seul paramètre d'interaction binaire, kij, qui dépend de la température. Afin de permettre au modèle E-PPR78 de prédire le comportement du gaz naturel, trois nouveaux groupes sont ajoutés : le monoxyde de carbone, l'hélium et l'argon. Pour cela, il a été nécessaire de former une base de données expérimentales la plus large possible contenant les mesures d'équilibres de phase et d'enthalpies de mélange pour les systèmes binaires constitués par ces trois groupes ainsi que ceux définis dans les études précédentes et présents dans le gaz naturel. Après une description de la classification des diagrammes de phase de Van Konynenburg et Scott, le modèle E-PPR78 est présenté. La troisième partie est consacrée à l'ajout des trois nouveaux groupes au sein du modèle. Les résultats sont obtenus avec une précision satisfaisante. Il apparaît clairement que le modèle E PPR78 est capable de prédire le comportement du gaz naturel dans des conditions de températures et de pressions particulièrement larges / The development of the predictive E-PPR78 model, based on a contribution group method, has been undertaken since ten years to predict accurately the behaviour of multi-component systems. This model lies on the Peng-Robinson equation of state with classical Van der Waals mixing rules. It uses a unique binary interaction parameter, kij, which is temperature dependant. To enable the E-PPR78 model to predict the behavior of natural gases, three new groups are added: carbon monoxide, helium and argon. It was necessary to build an experimental database, as exhaustive as possible, containing phase equilibrium and enthalpies of mixing data for binary systems formed by these groups and those defined in previous studies and present in natural gases. After a description of the classification scheme of Van Konynenburg and Scott, the E-PPR78 model is described. The third part is about the addition of the three new groups to the model. It clearly appears that the E-PPR78 model is able to predict the fluid-phase behavior of natural gases over wide ranges of temperatures and pressures
|
409 |
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
|
410 |
Étude mathématique et numérique des résonances dans une micro-cavité optique / Mathematical and numerical study of resonances in optical micro-cavitiesMoitier, Zoïs 03 October 2019 (has links)
Cette thèse est consacrée à l'étude des fréquences de résonance de cavités optiques bidimensionnelles. Plus particulièrement, on s'intéresse aux résonances à modes de galerie (modes localisés au bord de la cavité avec un grand nombre d'oscillations). La première partie traite du calcul numérique des résonances par la méthode des éléments finis à l'aide de couches parfaitement adaptées, et d'une analyse de sensibilité des paramètres de celles-ci dans les trois situations suivantes : un problème unidimensionnel, une réduction du cas bidimensionnel invariant par rotation et le cas général. La deuxième partie porte sur la construction de développements asymptotiques des résonances à modes de galerie quand le nombre d'oscillations le long du bord tend vers l'infini. On considère d'abord le cas d'un problème invariant par rotation pour lequel le nombre d'oscillations s'interprète comme un paramètre semiclassique grâce à la transformée de Fourier angulaire. Ensuite, pour le cas général, la construction utilise un ansatz phase-amplitude de type BKW qui permet de se ramener à un opérateur de Schrödinger généralisé. Enfin, les résonances calculées numériquement dans la première partie sont comparées aux développements asymptotiques explicités par calcul formel. / This thesis is devoted to the study of resonance frequencies of bidimensional optical cavities. More specifically, we are interested in whispering-gallery modes (modes localized along the cavity boundary with a large number of oscillations). The first part deals with the numerical computation of resonances by the finite element method using perfectly matched layers, and with a sensibility analysis in the three following situations: an unidimensional problem, a reduction of the rotationally invariant bidimensional case, and the general case. The second part focuses on the construction of asymptotic expansions of whispering-gallery modes as the number of oscillations along of boundary goes to infinity. We start by considering the case of a rotationally invariant problem for which the number of oscillations can be interpreted as a semiclassical parameter by means of an angular Fourier transform. Next, for the general case, the construction uses a phase-amplitude ansatz of WKB type which leads to a generalized Schrödinger operator. Finally, the numerically computed resonances obtained in the first part are compared to the asymptotic expansions made explicit by the use of a computer algebra software.
|
Page generated in 0.0851 seconds