The research presented herein describes the coupling of acoustic and heat release fluctuations in a perfectly-premixed swirl-stabilized combustor by analysis of simultaneous high-repetition-rate laser diagnostics data. Nine cases are studied, varying the thermal power and the equivalence ratio. Proper orthogonal decomposition (POD) of the velocity data shows that cases with higher amplitude thermoacoustic oscillations have flow fields containing helical vortex cores (HVC); these cases are further analysed to determine the driving mechanisms of the oscillations. Flow and flame statistics are compiled as a function of both the phase in the thermoacoustic cycle and a phase representing the azimuthal position of the HVC relative to the measurement plane. These data are used to spatially map the thermoacoustic energy transfer field, as described by the Rayleigh integral. It is found that periodic deformations of the HVC cause large-scale flame motions, resulting in regions of positive and negative energy transfer.

The research presented herein describes the coupling of acoustic and heat release fluctuations in a perfectly-premixed swirl-stabilized combustor by analysis of simultaneous high-repetition-rate laser diagnostics data. Nine cases are studied, varying the thermal power and the equivalence ratio. Proper orthogonal decomposition (POD) of the velocity data shows that cases with higher amplitude thermoacoustic oscillations have flow fields containing helical vortex cores (HVC); these cases are further analysed to determine the driving mechanisms of the oscillations. Flow and flame statistics are compiled as a function of both the phase in the thermoacoustic cycle and a phase representing the azimuthal position of the HVC relative to the measurement plane. These data are used to spatially map the thermoacoustic energy transfer field, as described by the Rayleigh integral. It is found that periodic deformations of the HVC cause large-scale flame motions, resulting in regions of positive and negative energy transfer.

Knowledge of the indoor propagation of sound has many important applications including acoustic prediction in homes, office buildings, stores, and schools, and the design of concert halls, auditoriums, classrooms, and factories. At low frequencies, interior acoustics are analyzed with the wave equation, but significant computational expense imposes an upper frequency limit. Thus, energy methods are often sought for high frequency analysis. However, conventional energy methods are significantly limited by vast simplifications or computational costs. Therefore, new improvements are still being sought. The basis of this dissertation is a recently developed mathematical model for interior acoustics known as the acoustic diffusion model. The model extends statistical methods in high frequency acoustics to predict the spatial distribution of acoustic energy in the volume over time as a diffusion process. Previously, solutions to the acoustic diffusion model have been limited to one dimensional (1-D) analytical solutions and to the use of the finite element method (FEM). This dissertation focuses on a new, efficient method for solving the acoustic diffusion model based on a boundary element method (BEM) solution using the Laplace transform. First, a Laplace domain solution to the diffusion model is obtained using the BEM. Then, a numerical inverse Laplace transform is used to efficiently compute the time domain response. The diffusion boundary element-Laplace transform solution (BE-LTS) is validated through comparisons with Sabine theory, ray tracing, and a diffusion FEM solution. All methods demonstrate excellent agreement for three increasingly complex acoustic volumes and the computational efficiency of the BE-LTS is exposed. Structural coupling is then incorporated in the diffusion BE-LTS using two methods. First, a simple transmission coefficient separating two acoustic volumes is implemented. Second, a structural power flow model represents the coupling partition separating acoustic volumes. The validation of these methods is successfully performed in an example through comparisons with statistical theory, a diffusion FEM solution, ray tracing, and experimental data. Finally, the diffusion model and the BE-LTS are shown to possess capabilities beyond that of room acoustics. The acoustic transmission through a heat exchanger, acoustic foam, and mufflers is successfully modeled using the diffusion BE-LTS and compared to experimental data. / Ph. D.

This thesis highlights one of the most significant concerns that has preoccupied drivetrain engineers in recent times, namely drivetrain clonk. Clonk is an unacceptable audible sound, which is accompanied by a tactile drivetrain response. This may occur under several different driving conditions. Many drivetrain NVH concerns are related to impact loading of subsystems down-line of engine. These concerns are induced by power torque surge through engagement and disengagement processes, which may propagate through various transmission paths as structural waves. The coincidence of these waves with the acoustic modes of sub-system components leads to audible responses, referred to as clonk. The approach usually undertaken and reported in literature is either purely theoretical or constitutes experimental observation of vehicle conditions. A few research workers have reported rig-based investigations, but not under fully dynamic conditions with controlled and reproducible impulsive action. The research reported in this thesis combines experimental and numerical investigation of high frequency behaviour of light truck drivetrain systems, when subjected to sudden impulsive action, due to driver behaviour. The problem is treated as a multi-physics interactive phenomenon under transient conditions. The devised numerical method combines multi-body dynamics, structural modal analysis, impact dynamics in lash zones and acoustic analysis within an overall investigation framework. A representative drivetrain system rig is designed and implemented, and controlled tests simulating driver behaviour undertaken. The combined numerical predictions and experimental noise and vibration monitoring has highlighted the fundamental aspects of drivetrain behaviour. Good agreement is' also found between the detailed numerical approach and the experimental findings. Novel methods of measurement such as Laser Doppler Vibrometery have been employed. Simultaneous measurements of vibration and noise radiation confirm significant elasto-acoustic coupling at high impact energy levels. One of the major finds of the thesis is the complex nature of the clonk signal, being a combination of accelerative and ringing noise, with the latter also comprising of many other lower energy content as observed in the case of transmission rattle and bearing-induced responses. Therefore, the link between rattle and clonk, long suspected, but not hitherto shown has been confirmed in the thesis. Another major find of particular commercial interest is the insignificant contribution of torsional damping devices such as dual mass flywheels upon the accelerative component of the clonk response.

Fully turbulent inflow past symmetrically located side branches mounted in a duct can give rise to pronounced flow oscillations due to coupling between separated shear layers and standing acoustic waves. Experimental investigation of acoustically-coupled flows was conducted using digital particle image velocimetry (DPIV) in conjunction with unsteady pressure measurements. Global instantaneous, phase- and time-averaged flow images, as well as turbulence statistics, were evaluated to provide insight into the flow physics during flow tone generation. Onset of the locked-on resonant states was characterized in terms of the acoustic pressure amplitude, frequency and the quality factor of the resonant pressure peak. Structure of the acoustic noise source is characterized in terms of patterns of generated acoustic power. In contrast to earlier work, the present study represents the first application of vortex sound theory in conjunction with global quantitative flow imaging and numerical simulation of the 2D acoustic field. In addition to the basic side branch configuration, the effects of bluff rectangular splitter plates located along the centerline of the main duct was investigated. The first mode of the shear layer oscillation was inhibited by the presence of plates, which resulted in substantial reduction of the amplitude of acoustic pulsations and the strength of the acoustic source. These results can lead to the development of improved control strategies for coaxial side branch resonators. Motivation for the second part of this study stems from the paper manufacturing industry, where air clamp devices utilize high-speed jets to position paper sheets with respect to other equipment. Thus, vibration of the paper sheet and turbulent flow that emerged from a planar curved nozzle between a flexible wall and a solid surface containing a backward-facing step (BFS) were investigated using high-speed photography and DPIV, respectively. The emphasis was on the characterization of the flow physics in the air clamp device, as well as of the shape of the paper sheet. For the control case, that involved a solid wall with a geometry that represented the time-averaged paper profile, hydrodynamic oscillation frequencies were characterized using unsteady pressure measurements. Experimentally obtained frequencies of the paper sheet vibration were compared to the hydrodynamic frequencies corresponding to the oscillations of the shear layer downstream of the BFS.

This thesis studies interactions between coupled acoustic domain(s) and enclosing rigid or elastic boundary. Boundary element-finite element (BE-FE) sound-structure interaction models are developed by coupling frequency domain BE acoustic and FE structural models using linear inviscid acoustic and elasticity theories. Flexibility in analyses is provided by discontinuous triangular and quadrilateral elements in the BE method (BEM), and a rectangular plate and a triangular shell element in the FE method (FEM). An analytical formulation is developed for an extended fundamental sound-structure interaction problem that involves locally reacting sound absorptive treatment on interior elastic boundary. This new formulation is built upon existing analytical solutions for a configuration known as the cavity-backed-plate problem. Results from developed analytical formulation are compared against those from independent BE-FE analyses. Analytical and BE-FE analysis results for a selection of cavity-plate(s) interaction cases are given. Single- and multi-domain BE analyses of cavity-Helmholtz resonator interaction are provided as an alternative to modal method of acoustoelasticity. A discrete-form of the existing BE acoustic particle velocity formulation is presented and demonstrated on a basic case study. Both the existing and the discretized BE acoustic particle velocity formulations could be utilized in acoustic studies. A selection of case studies involving fundamental configurations are studied both analytically and computationally (by BE or BE-FE methods). These studies could provide a basis for benchmark case development in the field of acoustics.

Optical membranes are currently pursued for their ability to replace the conventional mirrors that are used to correct wave front aberration and space-based telescopes. Among some of the many benefits of using optical membranes, is their ability to considerably reduce the weight of the structure. As a secondary effect, the cost of transportation, which is of great interest in space applications, is reduced as well. Given the low density of these thin-film membranes, the lower end dynamics play a greater significant role than their rigid plate-like counterparts in achieving functional mirrors. Space-based mirrors are subjected to a series of disturbances. Among those encountered are thermal radiation, debris impact, and slewing maneuvers. Thus, dynamic control is essential for the adequate performance of thin-film membrane mirrors. With this in mind, the work described herein aims to improve the performance of optical membranes with an innovative, acoustical control approach to suppress vibration of optical membranes backed by an air cavity. This is achieved by using a centralized acoustic source in the cavity as the method of actuation. The acoustic actuation is of great interest since it does not mass load the membrane in the conventional way, as most methods of actuation would. To achieve this end goal, two structural-acoustic coupled models are developed to describe the dynamics of a pressurized optical membrane system. This is done through an impedance based modeling approach where the subsystems are modeled individually, and then coupled at the interface. The control of the membrane is implemented using a positive position feedback approach. The theory is also extended to positive velocity and positive acceleration feedback. Three experiments are carried out to validate the models previously mentioned. Successful implementation of a control experiment is also accomplished leading to considerable attenuations in the coupled membrane's dynamics. / Ph. D.

A very short and disagreeable audible and tactile response from a vehicle driveline may be excited when the throttle is abruptly applied or released, or when the clutch is rapidly engaged. The condition is most noticeable in low gear and in slow moving traffic, when other background engine and road noise levels are low. This phenomenon is known as clonk and is often associated with the first cycle of shuffle response, which is a low frequency longitudinal vehicle movement excited by throttle demand. It is often reported that clonk may coincide with each cycle of the shuffle response, and multiple clonks may then occur. The problem is aggravated by backlash and wear in the drivetrain, and it conveys a perception of low quality to the customer. Hitherto, reported investigations do not reveal or discuss the mechanism and causal factors of clonk in a quantitative manner, which would relate the engine impulsive torque to the elastic response of the driveline components, and in particular to the noise radiating surfaces. Crucially, neither have the issues of sensitivity, variability and non-linearity been addressed and published. It is also of fundamental importance that clonk is seen as a total system response to impulsive torque, in the presence of distributed lash at the vibro-elastic impact sites. In this thesis, the drivetrain is defined as the torque path from the engine flywheel to the road wheels. The drivetrain is a lightly damped and highly non-linear dynamic system. There are many impact and noise emitting locations in the driveline that contribute to clonk, when the system is subjected to shock torque loading. This thesis examines the clonk energy paths, from the initial impact to many driveline lash locations, and to the various noise radiating surfaces. Both experimental and theoretical methods are applied to this complex system. Structural and acoustic dynamics are considered, as well as the very important frequency couplings between elastic structures and acoustic volumes. Preliminary road tests had indicated that the clonk phenomenon was a, very short transient impact event between lubricated contacts and having a high frequency characteristic. This indicated that a multi-body dynamics simulation of the driveline, in conjunction with a high frequency elasto-acoustic coupling analysis, would be required. In addition, advanced methods of signal analysis would be required to handle the frequency content of the very short clonk time histories. These are the main novelties of this thesis. There were many successful outcomes from the investigation, including quantitative agreement between the numerical and experimental investigations. From the experimental work, it was established that vehicle clonk could be accurately reproduced on a driveline rig and also on a vehicle chassis dynamometer, under controlled test conditions. It then enabled Design of Experiments to be conducted and the principal causal factors to be identified. The experimental input and output data was also used to verify the mathematical simulation. The high frequency FE analysis of the structures and acoustic cavities were used to predict the dynamic modal response to a shock input. The excellent correlation between model and empirical data that was achieved, clearly established the clonk mechanism in mathematical physics terms. Localised impact of meshing gears under impulsive loads were found to be responsible for high frequency structural wave propagation, some of which coupled with the acoustics modes of cavities, when the speed of wave propagation reached supersonic levels. This finding, although previously surmised, has been shown in the thesis and constitutes a major contribution to knowledge.

La caractérisation d’objets enfouis à partir de mesures d’ondes diffractées est un problème présent dans de nombreuses applications comme l’exploration géophysique, le contrôle non-destructif, l’imagerie médicale, etc. Elle peut être obtenue numériquement par la résolution d’un problème inverse. Néanmoins, c’est un problème non linéaire et mal posé, ce qui rend la tâche difficile. Une reconstruction précise nécessite un choix judicieux de plusieurs paramètres très différents, dépendant des données de la méthode numérique d’optimisation choisie.La contribution principale de cette thèse est une étude de la reconstruction complète d’obstacles élastiques immergés à partir de mesures du champ lointain diffracté. Les paramètres à reconstruire sont la frontière, les coefficients de Lamé, la densité et la position de l’obstacle. On établit tout d’abord des résultats d’existence et d’unicité pour un problème aux limites généralisé englobant le problème direct d’élasto-acoustique. On analyse la sensibilité du champ diffracté par rapport aux différents paramètres du solide, ce qui nous conduit à caractériser les dérivées partielles de Fréchet comme des solutions du problème direct avec des seconds membres modifiés. Les dérivées sont calculées numériquement grâce à la méthode de Galerkine discontinue avec pénalité intérieure et le code est validé par des comparaisons avec des solutions analytiques. Ensuite, deux méthodologies sont introduites pour résoudre le problème inverse. Toutes deux reposent sur une méthode itérative de type Newton généralisée et la première consiste à retrouver les paramètres de nature différente indépendamment, alors que la seconde reconstruit tous les paramètre en même temps. À cause du comportement différent des paramètres, on réalise des tests de sensibilité pour évaluer l’influence de ces paramètres sur les mesures. On conclut que les paramètres matériels ont une influence plus faible sur les mesures que les paramètres de forme et, ainsi, qu’une stratégie efficace pour retrouver des paramètres de nature distincte doit prendre en compte ces différents niveaux de sensibilité. On a effectué de nombreuses expériences à différents niveaux de bruit, avec des données partielles ou complètes pour retrouver certains paramètres, par exemple les coefficients de Lamé et les paramètres de forme, la densité, les paramètres de forme et la localisation. Cet ensemble de tests contribue à la mise en place d’une stratégie pour la reconstruction complète des conditions plus proches de la réalité. Dans la dernière partie de la thèse, on étend ces résultats à des matériaux plus complexes, en particulier élastiques anisotropes. / The characterization of hidden objects from scattered wave measurements arises in many applications such as geophysical exploration, non destructive testing, medical imaging, etc. It can be achieved numerically by solving an Inverse Problem. However, this is a nonlinear and ill-posed problem, thus a difficult task. A successful reconstruction requires careful selection of very different parameters depending on the data and the chosen optimization numerical method.The main contribution of this thesis is an investigation of the full reconstruction of immersed elastic scatterers from far-field pattern measurements. The sought-after parameters are the boundary, the Lamé coefficients, the density and the location of the obstacle. First, existence and uniqueness results of a generalized Boundary Value Problem including the direct elasto-acoustic problem are established. The sensitivity of the scattered field with respect to the different parametersdescribing the solid is analyzed and we end up with the characterization of the corresponding partial Fréchet derivatives as solutions to the direct problem with modified right-hand sides. These Fréchet derivatives are computed numerically thanks to the Interior Penalty Discontinuous Galerkin method and the code is validated thanks to comparison with analytical solutions. Then, two solution methodologies are introduced for solving the inverse problem. Both are based on an iterative regularized Newton-type methodology and the first one consists in retrieving the parameters of different nature independently, while the second one reconstructs all parameters together. Due to the different behavior of the parameters, sensitivity tests are performed to assess the impact of the parameters on the measurements. We conclude that material parameters have a weaker influence on the measurements than shape parameters, and therefore, a successful strategy to retrieve parameters of distinct nature should take into account these different levels of sensitivity. Various experiments at different noise levels and with full or limited aperture data are carried out to retrieve some of the physical properties, e.g. Lamé coefficients with shape parameters, density with shape parameters a, density, shape and location. This set of tests contributes to a final strategy for the full reconstruction and in more realistic conditions. In the final part of the thesis, we extend the results to more complex material parameters, in particular anisotropic elastic.

Cette thèse étudie le raffinement local de maillage à la fois en espace et en temps pour l’équation de l’elastodynamique du second ordre pour le calcul haute performance. L’objectif est de mettre en place des méthodes numériques pour traiter des hétérogénéités de petite taille ayant un impact important sur la propagation des ondes. Nous utilisons une approche par éléments finis de Galerkin discontinus avec pénalisation pour leur flexibilité et facilité de parallélisation. La formulation éléments finis que nous proposons a pour particularité d’être élasto-acoustique, pour pouvoir prendre en compte des hétérogénéités acoustiques de petite taille. Par ailleurs, nous proposons un terme de pénalisation optimisé qui est mieux adapté à l’équation de l’élastodynamique, conduisant en particulier à une meilleure condition CFL. Nous avons aussi amélioré une formulation PML du second ordre pour laquelle nous avons proposé une nouvelle discrétisation temporelle qui rend la formulation plus stable. En tirant parti de la p-adaptivité et des maillages non-conformes des méthodes de Galerkin discontinues combiné à une méthode de pas de temps local, nous avons grandement réduit le coût du raffinement local. Ces méthodes ont été implémentées en C++, en utilisant des techniques de template metaprogramming, au sein d’un code parallèle à mémoire distribuée (MPI) et partagée (OpenMP). Enfin, nous montrons le potentiel de notre approche sur des cas tests de validation et sur des cas plus réalistes avec des milieux présentant des hydrofractures. / This thesis studies local mesh refinement both in time and space for the second order elastodynamic equation in a high performance computing context. The objective is to develop numerical methods to treat small heterogeneities that have global impact on wave propagation. We use an internal penalty discontinuous Galerkin finite element approach for its flexibity and parallelization capabilities. The elasto-acoustic finite element formulation we discuss is elasto-acoustic in order to handle local acoustic heterogeneities. We also propose an optimized penalty term more suited to the elastodynamic equation that results in better CFL condition. We improve a second order PML formulation with an original time discretization that results in a more stable formulation. Using the p-adaptivity and nonconforming mesh capabilities of discontinuous Galerkin methods combined with a local time stepping method, we greatly reduce the high computational cost of local refinements. These methods have been implemented in C++, using template metaprogramming, in a distributed memory (MPI) and shared memory (OpenMP) parallel code. Finally, we show the potential of our methods on validation test cases and on more realistic test cases with medium including hydrofractures.