Dans cette thèse, on se propose de participer à la popularisation des méthodes de résolution de problèmes de propagation d'onde basées sur des formulations intégrales en fournissant des indicateurs d'erreur a posteriori utilisable dans le cadre d'algorithmes de raffinement autoadaptatif. Le développement de tels indicateurs est complexe du fait de la non-localité des normes associées aux espaces de Sobolev et des opérateurs entrant en jeu. Des indicateurs de la littérature sont étendus au cas de la propagation d'une onde acoustique. On étend les preuves de convergence quasi-optimale (de la littérature) des algorithmes autoadaptatifs associés dans ce cas. On propose alors une nouvelle approche par rapport à la littérature qui consiste à utiliser une technique de localisation des normes, non pas basée sur des inégalités inverses, mais sur l'utilisation d'un opérateur Λ de localisation bien choisi.On peut alors construire des indicateurs d'erreur a posteriori fiables, efficaces, locaux et asymptotiquement exacts par rapport à la norme de Galerkin de l'erreur. On donne ensuite une méthode pour la construction de tels indicateurs. Les applications numériques sur des géométries 2D et 3D confirment l'exactitude asymptotique ainsi que l'optimalité du guidage de l'algorithme autoadaptatif.On étend ensuite ces indicateurs au cas de la propagation d'une onde électromagnétique. Plus précisément, on s'intéresse au cas de l'EFIE. On propose des généralisations des indicateurs de la littérature. On effectue la preuve de convergence quasi-optimale dans le cas d'un indicateur basé sur une localisation de la norme du résidu. On utilise le principe du Λ pour obtenir le premier indicateur d'erreur fiable, efficace et local pour cette équation. On en propose une seconde forme qui est également, théoriquement asymptotiquement exacte. / The aim of this work is to participate to the popularization of methods for the resolution of wave propagation problems based on integral equations formulations by developping a posteriori error estimates in the context of autoadaptive mesh refinement strategies. The development of such estimates is difficult because of the non-locality of the norms associated to the Sobolev spaces and of the involved integral operators. Estimates from the literature are extended in the case of the propagation of an acoustic wave. The proofs of quasi-optimal convergence of the autoadaptive algorithms are established. We then introduce a new approach with respect to the literature which is based on a new norm-localization technique based on the use of a well-chosen Λ operator and not on inverse inequalities as it was the case previously.We then establish new a posteriori error estimates which are reliable, efficient, local and asymptotically exact with respect to the Galerkin norm of the error. We give a method for the construction of such estimates. Numerical applications on 2D and 3D geometries confirm the asymptotic exactness and the optimality of the autoadaptive algorithm.These estimates are extended in the case of the propagation of an electromagnetic wave. More precisely, we are interested in the EFIE. We suggest generalization of the estimates of the literature. A proof for quasi-optimal convergence is given for an estimate based on a localization of the norm of the residual. The principle of Λ is used to construct the first reliable, efficient, local error estimate for this equation. We give a second forme which is eventually theoretically asymptotically exact.

No description available.

Low frequency acoustic room behaviour can be approximated using numerical methods. Traditionally, music studio control rooms are built with complex geometries, making their eigenmodes difficult to predict mathematically. Hence, a summation-by-parts method with simultaneous-approximation-terms is derived to approximate the time dependent acoustic wave equation in three dimensions. The derived model is limited to rectangular prismatic rooms but planned to be expanded to handle complex geometries in the future. Semi-reflecting boundary conditions are used, corresponding to tabulated reflection and absorption properties of real. walls. Two speakers are modeled as omnidirectional point sources placed on a boundary, to mimic common studio setups. Through tests and examination of eigenvalues of the matrices in the method, conditions for stability and reflection coefficients are derived. Simulations of sound pressure distribution produced by the model correlate well to room mode theory, suggesting the model to be accurate in the application of predicting low frequency acoustic room behaviour. However, the convergence rate of the model turns out to be lower than expected when point sources are introduced. Future steps towards applying the model to real music studio control rooms include modeling the walls as changes in density and wave speed rather than boundaries of the domain. This would potentially allow more complex geometries to be modeled within a larger, rectangular domain.

京都大学 / 新制・課程博士 / 博士(農学) / 甲第23527号 / 農博第2474号 / 新制||農||1087(附属図書館) / 学位論文||R3||N5358(農学部図書室) / 京都大学大学院農学研究科地域環境科学専攻 / (主査)教授 藤原 正幸, 教授 中村 公人, 准教授 宇波 耕一 / 学位規則第4条第1項該当 / Doctor of Agricultural Science / Kyoto University / DGAM

A stochastic method is developed, implemented and investigated here for solving Laplace, Poisson's, and standard parabolic wave equations. This method is based on the properties of random walk, diffusion process, Ito formula, Dynkin formula and Monte Carlo simulations. The developed method is a local method i:e: it gives the value of the solution directly at an arbitrary point rather than extracting its value from complete field solution and thus is inherently parallel. Field computation by this method is demonstrated for electrostatic and electrodynamic propagation problems by considering simple examples and numerical results are presented to validate this method. Numerical investigations are carried out to understand efficacy and limitations of this method and to provide qualitative understanding of various parameters involved in this method.

In this work the scattering of an incident plane wave propagating along a plane perpendicular to the xy-plane is studied. The wave is scattered from an infinitely long cylindrical object of arbitrary cross-section. Due to the arbitrary geometry of the obstacle, a finite differences numerical method is employed to approximate the solution of the scattering problems. The wave equation is expressed in terms of generalized curvilinear coordinates. Boundary conforming grids are generated using elliptic grid generators. Then, a explicit marching in time scheme is implemented over these grids. It is found that as time grows the numerical solution converges to a wave with harmonic time dependence. The amplitude of these waves is analyzed and graphed over generalized grids for different geometries. An important physical measure of the energy scattered, the differential scattering cross section, is also obtained. In particular, the method is applied to a circular cylindrical obstacle. For this case, the analytical solution can also be obtained by traditional spectral techniques. The method is validated by comparing this exact solution with the numerical approximation obtained from the application of it.

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.

No description available.

A new wave-kinetic numerical method for the propagation of optical waves in weakly inhomogeneous media is discussed, and it is applied to several canonical problems: the propagation of beam and plane waves through a weak 3-D ( or 2-D ) Gaussian eddy. The numerical results are also compared to those from a Monte-Carlo simulation and the first Born approximation. Within the validity of the Liouville approximation, the Wigner distribution function ( WDF ) is conserved along the conventional ray trajectories, and, thus, by discretizing the input WDF with Gaussian beamlets, we can represent the output WDF as a sum of Gaussians, from which irradiance can be obtained by analytical integration of each Gaussian with respect to wavevector. Although each Gaussian beamlet propagates along a geometrical optics ray trajectory, it can correctly describe diffraction effects, and the propagation of optical waves through caustics or ray crossings. The numerical results agree well with either the Monte-Carlo method or the first Born approximation in regions where one or both of these are expected to be valid. / M.S.

The acoustics of small cavities raises interest of the scientific community since it involves particular damping mechanisms. In fluid dynamics, when a small perturbation is propagating within a Newtonian and heat-conducting fluid bounded by a rigid and isothermal surface, viscous and thermal dissipative mechanisms are generated near the walls. Such effects can have significant impact on the acoustic behaviour of the system.<p>Several types of practical applications can be cited, among which: hearing aids, micro-electro-mechanical systems (transducers, microphones and loud-speakers), absorbing materials made of thin capillary nets or small pores, dissipative silencers, thermo-acoustic heat exchangers, or any kind of device bringing into play small resonant cavities filled with a dissipative fluid (micro-acoustics).<p><p>This study focuses on appropriated reductions of the physical equations, in order to enhance the efficiency of the numerical resolution without adversely affecting the accuracy. Moreover, the proposed strategies lead to numerically stable systems as they involve only one scalar partial order differential equation (or equivalent fluid equation). The emphasis is put on the physical aspect of those reductions, their range of applicability, benefits and drawbacks.<p>Two new reduced models are proposed to estimate the visco-thermal acoustic wave propagation. A first extension deals with waveguide geometries and relax the hypothesis of the fluid at rest. The second original formulation addresses visco-thermal acoustics in 3D arbitrary geometries. This model is based on different considerations coming from existing techniques as well as the estimation of a wall-distance field.<p><p>A second part aims at studying the acoustic behaviour of biphasic materials and more specifically poro-elastic materials. This type of acoustic component is widely used in industry because of their good absorbing properties in the medium- and high-frequency <p>ranges.<p>A preliminary bibliographic research deals with the derivation of the set of partial order differential equations that account for both fluid/structure interactions and the anisotropy of a given poro-elastic material. It has been shown that transversely orientated capillary materials (for instance catalyst substrates) can be simulated using the proposed reduction technique.<p>At last, the modelling of the acoustic transmission between two domains separated by perforated or micro-perforated plates or thin plates of poro-elastic materials is discussed. The analogy between the rigid perforated plate models with an equivalent fluid formulation has been presented. As a result, this model has been extended in order to account for flexural effects of the solid part.<p><p><p>Ce travail porte sur l'étude de certains phénomènes d'amortissements intervenant dans l'acoustique des petites cavités. En méchanique des fluides, lorsqu'une petite perturbation se propage au sein d'un fluide newtonien et caloporteur borné par un mur rigide et isotherme, ces mécanismes dissipatifs particuliers se localisent aux abords des parois et jouent un rôle significatif dans certaines situations.<p>Parmi les exemples d'applications pratiques, il est possible de citer les appareils d'aide auditive, les systèmes microélectromécaniques (transducteurs, microphones et haut-parleurs), les matériaux absorbants constitués de fins réseaux capillaires ou de pores aux dimensions réduites, les systèmes de silencieux, d'échangeurs de chaleur thermo-acoustiques ou tout autre appareil mettant en jeu des cavités résonantes aux dimensions réduites (micro-acoustique).<p><p>L'étude proposée ici se focalise sur des stratégies de réduction appropriées des équations physiques, ceci afin d'améliorer l'efficacité du modèle tout en conservant une précision acceptable. Les techniques présentées aboutissent à des systèmes numériquement stables mettant en jeu une seule équation scalaire (ou équation fluide équivalent). Ainsi, l'accent est porté sur l'aspect physique des réductions, leurs domaines d'application, avantages et inconvénients.<p>Deux modèles originaux sont proposés afin de prédire la propagation acoustique visco-thermique. Une première extension permet d'évaluer la pression acoustique au sein de géométries particulières de type guides d'onde en présence d'un écoulement hydrodynamique. La seconde formulation présentée s'intéresse à l'acoustique dans des domaines 3D arbitraires. Cette méthode se base sur des considérations conjointes de modèles réduits existants ainsi que sur l'estimation d'un champ de distance à la plus proche paroi.<p><p>Dans une seconde partie, nous nous proposons d'étudier le comportement acoustique de matériaux biphasique et plus précisément les matériaux poro-élastiques (très utilisés dans l'industrie en raison de leurs caractéristiques absorbantes dans les domaines des moyennes et hautes fréquences).<p>Une étude bibliographique préliminaire nous a permis d'exprimer l'ensemble des équations aux dérivées <p>partielles modélisant à la fois les interactions fluide/structure et l'anisotropie générale des matériaux. <p>Cette réflexion nous a permis d'aboutir à un modèle de matériau isotrope transverse intéressant, combinant le modèle fluide proposé et la formulation acousto-élastique équivalente. Ainsi la modélisation de structures capillaires orientées (comme les matériaux utilisés dans les catalyseurs automobiles) s'en trouve grandement simplifiée.<p>Enfin, la transmission acoustique intervenant entre deux domaines fluides séparés par une plaque perforée ou micro-perforée ou bien une couche de matériau poreux a été étudiée. L'analogie entre les modèlisations existantes et un modèle générique de fluide équivalent a été mise en évidence. Pour finir, cette formulation a été étendue afin de prendre en compte les effets de flexion de la partie solide. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished